基于恒星时角的姿态控制方法和系统
专利汇可以提供基于恒星时角的姿态控制方法和系统专利检索，专利查询，专利分析的服务。并且本 发明 公开一种基于 恒星 时 角 的 姿态 控制方法，包含：S1、获取卫星的当前星上时间对应的儒略世纪；S2、计算恒星时角；S3、计算得到的轨道 坐标系 与惯性坐标系之间的姿态转换四元数；S4、计算得到卫星本体系相对轨道系的姿态四元数；S5、计算 俯仰 ‑滚动‑ 偏航 转序下卫星本体系下的欧拉角，得到基于星敏感器信息的姿态确定角；S6、执行PD控制 算法 ，计算得到控制 力 矩；重复S1至S6，实现卫星的自主连续导航和姿态控制。本发明仅通过恒星时角的计算，进行 航天器 位置 自主导航和姿态控制，解决异常状态下因无法获取地面或GNSS测定轨信息的问题，从而实现卫星的自主管理，提高卫星的可靠性和生存能力。，下面是基于恒星时角的姿态控制方法和系统专利的具体信息内容。
1.一种基于恒星时角的姿态控制方法，其特征在于，该姿态控制方法包含：
S1、获取卫星的当前星上时间对应的儒略世纪；
所述S1中，获取卫星的当前星上时间对应的儒略世纪T，如式(1)：
式(1)中：t为相对2000年1月1日12时UTC时间的星上当前时刻，单位为秒；
S2、利用当前星上时间对应的儒略世纪和卫星在地球同步轨道的定点位置信息，计算恒星时角；
所述S2中，利用当前星上时间对应的儒略世纪和卫星在地球同步轨道的定点位置信息，根据式(2)计算恒星时角SG包含：
SG＝((18h.697374558+879000h.051336907T)×15+P)×π/180(2)
式(2)中：T为星上时间对应的儒略世纪；P为卫星在地球同步轨道定点位置，单位为度；
S3、利用恒星时角，计算得到的轨道坐标系与惯性坐标系之间的姿态转换四元数；
所述S3中，利用恒星时角根据式(3)计算得到的轨道坐标系与惯性坐标系之间的姿态转换四元数qoi包含：
式(3)中：SG为恒星时角，T为星上时间对应的儒略世纪；
S4、利用轨道坐标系与惯性坐标系之间的姿态转换四元数，结合星敏感器测量得到的通过安装误差修正后的本体系相对惯性系姿态信息，计算得到卫星本体系相对轨道系的姿态四元数；
所述S4中，计算得到卫星本体系相对轨道系的姿态四元数包含：
将星敏感器输出姿态四元数qoi经安装偏差修正后的值qbi变换为轨道计算时刻卫星本体坐标系相对轨道坐标系的姿态四元数，用qbo表示；
姿态四元数qbo计算过程如式(4)：
式(4)中：qbi为星敏感器输出姿态四元数qoi经安装偏差修正后的值；qio表示轨道坐标系变换到惯性坐标系的姿态四元数；
S5、利用卫星本体系相对轨道系的姿态四元数，计算俯仰-滚动-偏航转序下卫星本体系下的欧拉角，得到基于星敏感器信息的姿态确定角，作为系统控制输入；
所述S5中，计算俯仰-滚动-偏航转序下卫星本体系下的欧拉角，得到基于星敏感器信息的姿态确定角包含：
将卫星本体系相对轨道系的姿态四元数qbo转换成姿态确定角，姿态确定角用
表示，转换方法如式(5)：
θST＝atan2(2(qbo,1qbo,3+qbo,2qbo,4),-(qbo,1)2-(qbo,2)2+(qbo,3)2+(qbo,4)2)ψST＝atan2(2(qbo,1qbo,2+qbo,3qbo,4),-(qbo,1)2+(qbo,2)2-(qbo,3)2+(qbo,4)2)(5)式(5)中， 为基于星敏感器姿态确定角的滚动角，θST为姿态确定角的俯仰角，ψST为姿态确定角的偏航角；qbo＝[qbo,1 qbo,2 qbo,3 qbo,4]；qbo,4为四元素标量部分，qbo,1,qbo,2,qbo,3为四元素矢量部分；
函数atan2(Y,X)计算方法如式(6)：
S6、利用由前后拍星敏感器确定的姿态确定角微分计算得到的姿态确定角速度，结合由基于星敏感器信息的姿态确定角确定的姿态确定角估值，执行PD控制算法，计算得到控制力矩；
所述S6包含：
利用基于星敏感器的姿态确定角计算姿态确定角估值与姿态确定角速度，姿态确定角估值用 表示，姿态确定角速度用 表示，如式(7)：
式(7)中， 为基于星敏感器姿态确定角的滚动角，θST为姿态确定角的俯仰角，ψST为姿态确定角的偏航角；Ts为控制周期；
执行PD控制算法，控制力矩用Tci，i＝x,y,z表示，如式(8)：
式(8)中，KPi,KDi，i＝x,y,z为控制器参数；qPi,qDi，i＝x,y,z为整定系数，缺省为1；
重复S1至S6，实现卫星的自主连续导航和姿态控制。
2.一种基于恒星时角的姿态控制系统，其特征在于，该姿态控制系统包含：
儒略世纪获取模块，其获取卫星的当前星上时间对应的儒略世纪；获取卫星的当前星上时间对应的儒略世纪T，如式(1)：
式(1)中：t为相对2000年1月1日12时UTC时间的星上当前时刻，单位为秒；
恒星时角计算模块，其利用当前星上时间对应的儒略世纪和卫星在地球同步轨道的定点位置信息，计算恒星时角；利用当前星上时间对应的儒略世纪和卫星在地球同步轨道的定点位置信息，根据式(2)计算恒星时角SG包含：
SG＝((18h.697374558+879000h.051336907T)×15+P)×π/180(2)
式(2)中：T为星上时间对应的儒略世纪；P为卫星在地球同步轨道定点位置，单位为度；
姿态转换四元数计算模块，其利用恒星时角，计算得到的轨道坐标系与惯性坐标系之间的姿态转换四元数；利用恒星时角根据式(3)计算得到的轨道坐标系与惯性坐标系之间的姿态转换四元数qoi包含：
式(3)中：SG为恒星时角，T为星上时间对应的儒略世纪；
姿态四元数计算模块，其利用轨道坐标系与惯性坐标系之间的姿态转换四元数，结合星敏感器测量得到的通过安装误差修正后的本体系相对惯性系姿态信息，计算得到卫星本体系相对轨道系的姿态四元数；计算得到卫星本体系相对轨道系的姿态四元数包含：
将星敏感器输出姿态四元数qoi经安装偏差修正后的值qbi变换为轨道计算时刻卫星本体坐标系相对轨道坐标系的姿态四元数，用qbo表示；姿态四元数qbo计算过程如式(4)：
式(4)中：qbi为星敏感器输出姿态四元数qoi经安装偏差修正后的值；qio表示轨道坐标系变换到惯性坐标系的姿态四元数；
姿态确定角获取模块，其利用卫星本体系相对轨道系的姿态四元数，计算俯仰-滚动-偏航转序下卫星本体系下的欧拉角，得到基于星敏感器信息的姿态确定角，作为系统控制输入；计算俯仰-滚动-偏航转序下卫星本体系下的欧拉角，得到基于星敏感器信息的姿态确定角包含：
将卫星本体系相对轨道系的姿态四元数qbo转换成姿态确定角，姿态确定角用
表示，转换方法如式(5)：
θST＝atan2(2(qbo,1qbo,3+qbo,2qbo,4),-(qbo,1)2-(qbo,2)2+(qbo,3)2+(qbo,4)2)ψST＝atan2(2(qbo,1qbo,2+qbo,3qbo,4),-(qbo,1)2+(qbo,2)2-(qbo,3)2+(qbo,4)2)(5)式(5)中， 为基于星敏感器姿态确定角的滚动角，θST为姿态确定角的俯仰角，ψST为姿态确定角的偏航角；qbo＝[qbo,1 qbo,2 qbo,3 qbo,4]；qbo,4为四元素标量部分，qbo,1,qbo,2,qbo,3为四元素矢量部分；
函数atan2(Y,X)计算方法如式(6)：
控制力矩获取模块，其利用由前后拍星敏感器确定的姿态确定角微分计算得到的姿态确定角速度，结合由基于星敏感器信息的姿态确定角确定的姿态确定角估值，执行PD控制算法，计算得到控制力矩；利用基于星敏感器的姿态确定角计算姿态确定角估值与姿态确定角速度，姿态确定角估值用 表示，姿态确定角速度用 表示，如式(7)：
式(7)中， 为基于星敏感器姿态确定角的滚动角，θST为姿态确定角的俯仰角，ψST为姿态确定角的偏航角；Ts为控制周期；
执行PD控制算法，控制力矩用Tci，i＝x,y,z表示，如式(8)：
式(8)中，KPi,KDi，i＝x,y,z为控制器参数；qPi,qDi，i＝x,y,z为整定系数，缺省为1。
基于恒星时角的姿态控制方法和系统
技术领域
背景技术
发明内容
具体实施方式
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 一种基于恒星地面观测系统的焦距测量方法 | 2020-05-15 | 107 |
| 一种恒星敏感器的模拟器 | 2020-05-12 | 702 |
| 恒星遥感影像的点扩散函数提取方法 | 2020-05-12 | 103 |
| 一种手持式恒星和星座指示装置 | 2020-05-12 | 776 |
| 一种恒星陀螺及其实现方法 | 2020-05-13 | 589 |
| 一种恒星陀螺敏感器 | 2020-05-11 | 701 |
| 一种恒星轴调差差速器 | 2020-05-11 | 722 |
| 基于恒星时角的姿态控制方法和系统 | 2020-05-12 | 837 |
| 一种恒星时估算方法 | 2020-05-11 | 938 |
| 恒星日与太阳日示教板 | 2020-05-14 | 139 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
