CN113485095A - 北斗三号卫星处于地影期时预报姿态的方法 - Google Patents
北斗三号卫星处于地影期时预报姿态的方法 Download PDFInfo
- Publication number
- CN113485095A CN113485095A CN202110920677.7A CN202110920677A CN113485095A CN 113485095 A CN113485095 A CN 113485095A CN 202110920677 A CN202110920677 A CN 202110920677A CN 113485095 A CN113485095 A CN 113485095A
- Authority
- CN
- China
- Prior art keywords
- satellite
- beidou
- attitude
- yaw angle
- angle
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 88
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 39
- 238000012545 processing Methods 0.000 claims abstract description 12
- 230000001133 acceleration Effects 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims description 27
- 238000005259 measurement Methods 0.000 claims description 16
- 230000005540 biological transmission Effects 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 7
- 230000009466 transformation Effects 0.000 claims description 6
- 238000005096 rolling process Methods 0.000 claims description 4
- 238000000844 transformation Methods 0.000 claims description 4
- 230000003111 delayed effect Effects 0.000 claims description 2
- 230000007246 mechanism Effects 0.000 claims description 2
- 238000004088 simulation Methods 0.000 claims description 2
- 238000013277 forecasting method Methods 0.000 claims 1
- 230000014509 gene expression Effects 0.000 claims 1
- 238000011084 recovery Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 8
- 230000000712 assembly Effects 0.000 description 2
- 238000000429 assembly Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000013524 data verification Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- DMBHHRLKUKUOEG-UHFFFAOYSA-N diphenylamine Chemical compound C=1C=CC=CC=1NC1=CC=CC=C1 DMBHHRLKUKUOEG-UHFFFAOYSA-N 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B11/00—Automatic controllers
- G05B11/01—Automatic controllers electric
- G05B11/36—Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential
- G05B11/42—Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential for obtaining a characteristic which is both proportional and time-dependent, e.g. P. I., P. I. D.
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明提供了一种北斗三号卫星处于地影期时偏航角正负阶跃时姿态预报的方法，包括：通过PID算法，结合卫星参数，完成地影期偏航正负阶跃时每一状态的迭代，具体算法包括：计算当前时刻指令；指令一阶惯性处理，计算角加速度；计算下一时刻角速度，偏航角；以及循环迭代计算每一时刻偏航角步骤。
Description
技术领域
本发明涉及北斗卫星技术领域，特别涉及一种北斗三号卫星处于地影期偏航角正负阶跃时预报姿态的方法。
背景技术
在现有技术中，为了满足导航载荷与测控天线对地指向的需求，同时保证帆板的对日指向以获取能源，北斗三号导航系统MEO卫星常规运行过程中姿态采用动态偏航模式，在太阳高度角较小时采用连续动偏控制模式。导航卫星偏航姿态模式的变化会引起卫星太阳能帆板和卫星本体受照面的变化，从而影响卫星在轨的受力状态。随着北斗导航卫星的发展，对定轨要求的提高，对光压模型和偏航姿态模式研究分析进而提高北斗的定轨精度成为重要课题。
现有技术具有如下缺陷，当轨道太阳角β处于[-3°,+3°]时，卫星处于深地影区。在一年内，北斗MEO卫星轨道太阳角经历两次符号变换，即两段地影期。偏航角也经历两次符号变换，常规预报算法直接计算会产生偏航角符号变化。实际姿控伺服系统的为减少正负变换的负荷压力，采用连续动态偏航，以使当前状态偏航角与上一状态偏航角尽量保持平稳变化。该阶段分为两部分，其一为姿态角缓慢变化到±5°，其二为±5°以内执行的阶跃变化。
在姿态角正负阶跃的过程内，计算偏航角最大从+5°跳变为-5°，而在轨实际偏航角以连续曲线变换，所以预报误差最大可达到10°。如图5所示，偏航角正负阶跃时，常规预报算法计算偏航角与在轨实际偏航角趋势差异明显。
发明内容
本发明的目的在于提供一种北斗三号卫星处于地影期时预报姿态的方法，以解决现有的北斗三号卫星偏航角正负阶跃时，常规预报算法计算偏航角与在轨实际偏航角趋势差异明显的问题。
为解决上述技术问题，本发明提供一种北斗三号卫星处于地影期时预报姿态的方法，包括：
轨道太阳角指的是卫星到太阳的矢量与卫星轨道面的夹角；卫星轨道坐标系+Z轴指向地心，+Y轴垂直于轨道面，与轨道面法线方向相反，+X轴由右手定则确定；由于卫星天线指地需求，卫星质心坐标系+Z轴与轨道坐标系+Z轴重合，卫星滚动角和俯仰角都为0。
偏航角为绕+Z轴旋转，卫星质心坐标系+X与轨道坐标系+X轴的夹角，在地影期偏航角正负阶跃时计算偏航角需如下步骤：
通过PID算法，结合卫星参数，完成地影期偏航正负阶跃时每一状态的迭代，具体算法包括：
计算当前时刻指令；
指令一阶惯性处理；
计算下一时刻角速度，偏航角；以及
循环迭代计算至计算当前时刻指令步骤；
地影期偏航角正负阶跃时，根据PID算法计算偏航角与在轨实际偏航角趋势。
可选的，在所述的北斗三号卫星处于地影期时预报姿态的方法中，
计算当前时刻指令包括：
eψi＝ψi-ψAim (9)
eωi＝ωi-ωaim (10)
其中Kp、Ki、Kd为PID算法中比例系数、积分系数和微分系数，ψi为当前时刻偏航角，ψaim为目标偏航角，eψi为当前时刻偏航角与目标偏航角差值，ωi当前时刻偏航角速度，ωaim为目标偏航角速度。
可选的，在所述的北斗三号卫星处于地影期时预报姿态的方法中，指令一阶惯性处理包括：
经过一阶惯性处理，Insi变换为Insouti，则角加速度αi为
αi＝Insouti/J (11)
其中J为转动惯量。
可选的，在所述的北斗三号卫星处于地影期时预报姿态的方法中，计算下一时刻角速度，偏航角包括：
下一状态角加速度为
ωi+1＝ωi+αi*T (12)
下一状态偏航角
ψi+1＝ψi+ωi*T+0.5*αi*T2 (13)。
可选的，在所述的北斗三号卫星处于地影期时预报姿态的方法中，循环迭代计算至计算当前时刻指令步骤包括：
可选的，在所述的北斗三号卫星处于地影期时预报姿态的方法中，利用北斗三号MEO卫星处于地影期时预报姿态的方法，以将预报姿态用于星间链路指向计算；
利用北斗三号卫星处于地影期时预报姿态的方法，以解决现有算法在地影期偏航正负阶跃时预报误差极大。
可选的，在所述的北斗三号卫星处于地影期时预报姿态的方法中，星间链路指向计算包括：
根据卫星在地影期时轨道信息和偏航设置，计算北斗三号卫星的当前姿态信息；以及
将北斗三号卫星的当前姿态信息用于星间链路指向，以替代或更正由于星内姿态信息传输失败引起的建链指向错误。
可选的，在所述的北斗三号卫星处于地影期时预报姿态的方法中，还包括：
在地面站无法获取北斗三号遥测时，包括由于境外遥测中断、星间链路误码、星地链路误码情况导致的姿态信息获取失败，用所述方法根据卫星位置信息计算相应时刻的姿态信息以用于地面事后或事前计算，包括定轨、光压估计、指向仿真、卫星建模；
在无权获取北斗三号遥测的单位，包括工厂、软件中，用所述方法进行姿态计算，以完成其所需。
可选的，在所述北斗三号卫星处于地影期偏航角正负阶跃时的预报姿态的方法中，还包括：
当轨道太阳角β处于[-3°,+3°]时，卫星处于深地影区。在一年内，北斗三号MEO卫星轨道太阳角经历两次符号变换，即两段地影期。此时卫星设置太阳高度角为
则卫星偏航角为
当β正负切换时，偏航角将产生巨大的符号变换，这将给卫星的姿控伺服系统造成巨大的负荷压力，是姿控伺服机构承受不了的。所以在卫星姿控系统设置中，除了轨道太阳角度外，连续动态偏航还需考虑上一时刻偏航角与这一时刻偏航角尽量平稳变化，将偏航角正负切换延时到±5°执行。
此时，真正的偏航角需在式(9)的基础上增加惯性处理，具体计算方式为
因此，卫星的偏航角将产生从+5°到-5°(或从-5°到+5°)的阶跃变化。此即为本发明所述方法的进入条件。
在本发明提供的北斗三号卫星处于地影期时预报姿态的方法中，通过PID算法，结合卫星参数，完成地影期偏航正负阶跃时每一状态的迭代，计算当前时刻指令，指令一阶惯性处理，计算下一时刻角速度和偏航角，循环迭代计算至计算当前时刻指令步骤，地影期偏航角正负阶跃时，根据PID算法计算偏航角与在轨实际偏航角趋势，解决了现有算法在地影期偏航正负阶跃时预报误差极大的情况。使用本发明可以使北斗三号MEO卫星预报算法在正负阶跃时预报误差减小为1.7°以下。
附图说明
图1是现有的姿控信息在星内传播路径示意图；
图2是卫星偏航角示意图；
图3是本发明一实施例中的星间链路恢复方法使用预报姿态计算离轴角误差示意图；
图4是本发明一实施例中的星间链路恢复方法使用预报姿态计算旋转角误差示意图；
图5是现有的常规预报算法计算偏航角与在轨实际偏航角趋势示意图；
图6是本发明一实施例中的北斗卫星处于地影期时在偏航角正负阶跃时的预报算法示意图；
图7是本发明一实施例中的地影期偏航角正负阶跃时根据PID算法计算偏航角与在轨实际偏航角趋势示意图。
具体实施方式
下面结合具体实施方式参考附图进一步阐述本发明。
应当指出，各附图中的各组件可能为了图解说明而被夸大地示出，而不一定是比例正确的。在各附图中，给相同或功能相同的组件配备了相同的附图标记。
在本发明中，除非特别指出，“布置在…上”、“布置在…上方”以及“布置在…之上”并未排除二者之间存在中间物的情况。此外，“布置在…上或上方”仅仅表示两个部件之间的相对位置关系，而在一定情况下、如在颠倒产品方向后，也可以转换为“布置在…下或下方”，反之亦然。
在本发明中，各实施例仅仅旨在说明本发明的方案，而不应被理解为限制性的。
在本发明中，除非特别指出，量词“一个”、“一”并未排除多个元素的场景。
在此还应当指出，在本发明的实施例中，为清楚、简单起见，可能示出了仅仅一部分部件或组件，但是本领域的普通技术人员能够理解，在本发明的教导下，可根据具体场景需要添加所需的部件或组件。另外，除非另行说明，本发明的不同实施例中的特征可以相互组合。例如，可以用第二实施例中的某特征替换第一实施例中相对应或功能相同或相似的特征，所得到的实施例同样落入本申请的公开范围或记载范围。
在此还应当指出，在本发明的范围内，“相同”、“相等”、“等于”等措辞并不意味着二者数值绝对相等，而是允许一定的合理误差，也就是说，所述措辞也涵盖了“基本上相同”、“基本上相等”、“基本上等于”。以此类推，在本发明中，表方向的术语“垂直于”、“平行于”等等同样涵盖了“基本上垂直于”、“基本上平行于”的含义。
另外，本发明的各方法的步骤的编号并未限定所述方法步骤的执行顺序。除非特别指出，各方法步骤可以以不同顺序执行。
以下结合附图和具体实施例对本发明提出的北斗三号卫星处于地影期时预报姿态的方法作进一步详细说明。根据下面说明和权利要求书，本发明的优点和特征将更清楚。需说明的是，附图均采用非常简化的形式且均使用非精准的比例，仅用以方便、明晰地辅助说明本发明实施例的目的。
本发明的目的在于提供一种北斗三号卫星处于地影期时预报姿态的方法，以解决现有的北斗三号卫星偏航角正负阶跃时，常规预报算法计算偏航角与在轨实际偏航角趋势差异明显的问题。
为实现上述目的，本发明提供了一种北斗三号卫星处于地影期时预报姿态的方法，包括：通过PID算法，结合卫星参数，完成地影期偏航正负阶跃时每一状态的迭代，具体算法包括：计算当前时刻指令；指令一阶惯性处理；计算下一时刻角速度，偏航角；以及循环迭代计算至计算当前时刻指令步骤；地影期偏航角正负阶跃时，根据PID算法计算偏航角与在轨实际偏航角趋势。
本发明的实施例提供一种北斗三号卫星处于地影期时预报姿态的方法，如图6所示，包括：预报滚动和俯仰角都为0；通过PID算法，结合卫星参数，完成地影期偏航正负阶跃时每一状态的迭代，具体算法包括：计算当前时刻指令；指令一阶惯性处理；计算下一时刻角速度，偏航角；以及循环迭代计算至计算当前时刻指令步骤；地影期偏航角正负阶跃时，根据PID算法计算偏航角与在轨实际偏航角趋势。
在本发明的一个实施例中，在所述的北斗三号卫星处于地影期时预报姿态的方法中，计算当前时刻指令包括：
eψi＝ψi-ψAim (16)
eωi＝ωi-ωaim (17)
其中Kp、Ki、Kd为PID算法中比例系数、积分系数和微分系数，ψi为当前时刻偏航角，ψaim为目标偏航角，eψi为当前时刻偏航角与目标偏航角差值，ωi当前时刻偏航角速度，ωaim为目标偏航角速度。
指令一阶惯性处理包括：
经过一阶惯性处理，Insi变换为Insouti，则角加速度αi为
αi＝Insouti/J (18)
其中J为转动惯量。
计算下一时刻角速度，偏航角包括：
下一状态角速度为
ωi+1＝ωi+αi*T (19)
下一状态偏航角
ψi+1＝ψi+ωi*T+0.5*αi*T2 (20)
循环迭代计算至计算当前时刻指令步骤包括：
重复(1)-(7)步骤，循环计算整个阶跃过程的偏航角。
在本发明提供的北斗三号卫星处于地影期时预报姿态的方法中，通过PID算法，结合卫星参数，完成地影期偏航正负阶跃时每一状态的迭代，计算当前时刻指令，指令一阶惯性处理，计算下一时刻角速度和偏航角，循环迭代计算至计算当前时刻指令步骤，地影期偏航角正负阶跃时，根据PID算法计算偏航角与在轨实际偏航角趋势，解决了现有算法在地影期偏航正负阶跃时预报误差极大的情况。使用本发明可以使北斗三号MEO卫星预报算法在正负阶跃时预报误差减小为1.7°以下。
本发明的北斗三号卫星处于地影期时预报姿态的方法还可以解决现有算法在地影期偏航正负阶跃时预报误差极大的情况。使用本发明可以使北斗三号MEO卫星预报算法在正负阶跃时预报误差减小为1.7°以下。
如图7所示，地影期偏航角正负阶跃时，根据PID算法计算偏航角与在轨实际偏航角趋势。本发明创新性的提出了北斗卫星处于地影期时在偏航角正负阶跃时的预报算法。
在本发明的一个实施例中，在所述的北斗三号卫星处于地影期时预报姿态的方法中，利用北斗三号MEO卫星处于地影期时预报姿态的方法，以将预报姿态用于星间链路指向计算；利用北斗三号卫星处于地影期时预报姿态的方法，以解决现有算法在地影期偏航正负阶跃时预报误差极大。
北斗三号导航系统中，星间链路起到了连接组网的作用。受地面分布测站地理限制，处于境外的卫星需要由境内星中转上行控制信息和下行遥测信息。除了通信功能外，星地/星间的测距下行到地面后，采用联合定轨的方法能有效提高定轨精度，显著降低设备时延对轨道和钟差解算精度的影响。同时星间链路测距值也被用于北斗三号卫星的自主导航。在星间链路建链时，需要上位机能够按照时隙约定计算出正确的星间链路指向角。姿控信息在星内传播时，如总线发生错误或者总线控制单机出现错误，将会影响星间链路上位机收到正确的姿态信息，从而导致建链失败，影响星间测距通信。尤其卫星处于境外时，若发生该类故障，将导致测控站无法监测卫星状态，卫星也无法获取境内更新电文，严重降低系统可靠性。
图1为姿控信息在星内传播路径示意图，如图1所示，在红线标识的星内单机或总线故障时，姿态信息传输失败，将导致建链故障。因此本发明的目的还在于提供一种北斗三号卫星处于地影期时预报姿态的方法，以解决现有的北斗三号卫星由于星内姿态信息传输故障引起的星间链路指向失败问题。本发明通过北斗三号卫星处于地影期时预报姿态的方法实现一种星间链路指向恢复方法，尤其涉及北斗三号导航卫星星内姿态信息传输故障时星间链路恢复指向方法。
为实现上述目的，本发明还提供了一种北斗三号导航卫星星内姿态信息传输故障时星间链路恢复指向方法，包括：确保满足北斗三号卫星的导航载荷与测控天线对地指向的需求，并根据需求确定北斗三号卫星的轨道信息和约束项，根据轨道信息和约束项与太阳的关系，计算北斗三号卫星的当前姿态信息；以及将北斗三号卫星的当前姿态信息用于星间链路指向，以替代或更正由于星内姿态信息传输失败引起的建链指向错误。
本发明的实施例提供一种星间链路恢复方法，包括：确保满足北斗三号卫星的导航载荷与测控天线对地指向的需求，并根据需求确定北斗三号卫星的轨道信息和约束项，根据轨道信息和约束项与太阳的关系，计算北斗三号卫星的当前姿态信息；以及将北斗三号卫星的当前姿态信息用于星间链路指向，以替代或更正由于星内姿态信息传输失败引起的建链指向错误。
在本发明的一个实施例中，在所述的星间链路恢复方法中，还包括：在卫星正常情况下，北斗三号卫星的当前姿态信息对星内传输的姿态信息进行核验。
在本发明的一个实施例中，在所述的星间链路恢复方法中，确保满足北斗三号卫星的导航载荷与测控天线对地指向的需求包括：北斗三号导航系统MEO卫星常规运行过程中姿态采用动态偏航模式，在太阳高度角小于10°时采用连续动偏控制模式，以满足北斗三号卫星的导航载荷与测控天线对地指向的需求，并同时保证帆板的对日指向以获取能源。
在本发明的一个实施例中，在所述的星间链路恢复方法中，计算北斗三号卫星的当前姿态信息包括：根据卫星质心坐标系和轨道坐标系，计算北斗三号卫星的姿态角；所述北斗三号卫星的姿态角包括偏航角、滚动角和俯仰角。轨道太阳角与偏航角示意图如图2所示。
在本发明的一个实施例中，在所述的星间链路恢复方法中，采用连续动偏控制模式包括：所述卫星质心坐标系中，将测控天线部署在卫星质心坐标系的+Z轴方向；其中卫星质心坐标系的+Z轴与轨道坐标系的+Z轴重合，指向地球中心，以满足北斗三号卫星的测控天线对地指向的需求。
在本发明的一个实施例中，在所述的星间链路恢复方法中，采用连续动偏控制模式还包括：北斗三号卫星运动时，使太阳矢量一直保持在卫星质心坐标系的X-O-Z面内，以保证帆板的对日指向以获取能源，满足能源需求。
在本发明的一个实施例中，在所述的星间链路恢复方法中，还包括：设置初始状态，所述初始状态包括：卫星质心坐标系的+X轴与轨道坐标系的+X轴重合，北斗三号卫星绕卫星质心坐标系的+Z轴旋转的角度为北斗三号卫星的姿态角中的偏航角ψ；北斗三号卫星的姿态角中的滚动角和俯仰角都为0度。
在本发明的一个实施例中，在所述的星间链路恢复方法中，还包括：轨道太阳高度角β为太阳矢量与轨道坐标系X-O-Z平面的夹角。
在本发明提供的星间链路恢复方法中，提出了将预报姿态用于星间链路指向计算的方案。根据在轨数据核验，在一年中预报姿态与在轨实际姿态的误差在99.9987％的时间内误差都在0.2°以下，100％时间内误差在1.7°以下，如图3～4所示。使用预报姿态计算星间链路指向角，离轴角误差为0，旋转角误差在1.7°以下。使用本申请的方案，在星间链路系统获取姿态角失败的情况下，能够将星间链路的捕获成功概率从0提高到100％。
本发明还涉及一种北斗三号卫星处于地影期时预报姿态的方法，以实现将预报姿态用于星间链路指向计算的方案。尤其涉及北斗三号MEO卫星处于地影期时预报姿态的方法。
在本发明的一个实施例中，首先为了满足导航载荷与测控天线对地指向的需求，同时保证帆板的对日指向以获取能源，北斗三号导航系统MEO卫星常规运行过程中姿态采用动态偏航模式，在太阳高度角较小时采用连续动偏控制模式。其轨道信息已知，约束项已知，则根据轨道与太阳的关系，可以计算出当时的姿态信息。
其次将计算出的姿态信息用于星间链路指向，可解决由于星内姿态信息传输失败引起的建链指向错误。
另外在卫星正常情况下，计算出的姿态信息也可对星内传输的姿态信息进行核验。
计算姿态角涉及两个坐标系。一是卫星质心坐标系，二是轨道坐标系。由于卫星质心坐标系中，天线部署在+Z轴方向，故而为了指地需求，卫星质心坐标系的+Z轴与轨道坐标系+Z轴重合，指向地球中心。卫星运动时，为了能源需求，使太阳矢量一直保持在质心坐标系X-O-Z面内。设置初始状态卫星本体系+X与轨道坐标系+X轴重合，卫星绕Z轴旋转的角度即为偏航角ψ，滚动角和俯仰角都为0。轨道太阳高度角β为太阳矢量与轨道坐标系X-O-Z平面的夹角。
本发明提出了将预报姿态用于星间链路指向计算的方案。根据在轨数据核验，在一年中预报姿态与在轨实际姿态的误差在99.9987％的时间内误差都在0.2°以下，100％时间内误差在1.7°以下。使用预报姿态计算星间链路指向角，离轴角误差为0，旋转角误差在1.7°以下。使用该方案，在星间链路系统获取姿态角失败的情况下，能够将星间链路的捕获成功概率从0提高到100％。本发明创造性的提出了北斗三号导航卫星星内姿态信息传输故障时星间链路恢复指向的方法。
综上，上述实施例对星间链路恢复方法的不同构型进行了详细说明，当然，本发明包括但不局限于上述实施中所列举的构型，任何在上述实施例提供的构型基础上进行变换的内容，均属于本发明所保护的范围。本领域技术人员可以根据上述实施例的内容举一反三。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
上述描述仅是对本发明较佳实施例的描述，并非对本发明范围的任何限定，本发明领域的普通技术人员根据上述揭示内容做的任何变更、修饰，均属于权利要求书的保护范围。
Claims (10)
1.一种北斗三号卫星处于地影期偏航角正负阶跃时姿态预报的方法，其特征在于，包括：通过PID算法，结合卫星参数，完成地影期偏航正负阶跃时每一状态的迭代，具体算法包括：
计算当前时刻指令；
指令一阶惯性处理，计算角加速度；
计算下一时刻角速度，偏航角；以及
循环迭代计算每一时刻偏航角步骤。
2.如权利要求1所述的北斗三号卫星处于地影期偏航角正负阶跃时预报姿态的方法，其特征在于，
轨道太阳角指的是卫星到太阳的矢量与卫星轨道面的夹角；卫星轨道坐标系+Z轴指向地心，+Y轴垂直于轨道面，与轨道面法线方向相反，+X轴由右手定则确定；由于卫星天线指地需求，卫星质心坐标系+Z轴与轨道坐标系+Z轴重合，卫星滚动角和俯仰角都为0。
3.如权利要求1所述的北斗三号卫星处于地影期偏航角正负阶跃时预报姿态的方法，其特征在于，
计算偏航角需如下算法，
计算当前时刻指令包括：
eψi＝ψi-ψAim (2)
eωi＝ωi-ωaim (3)
其中Kp、Ki、Kd为PID算法中比例系数、积分系数和微分系数，ψi为当前时刻偏航角，ψaim为目标偏航角，eψi为当前时刻偏航角与目标偏航角差值，ωi当前时刻偏航角速度，ωaim为目标偏航角速度。
4.如权利要求1所述的北斗三号卫星处于地影期偏航角正负阶跃时预报姿态的方法，其特征在于，指令一阶惯性处理，计算角加速度包括：
经过一阶惯性处理，Insi变换为Insouti，则角加速度αi为
αi＝Insouti/J (4)
其中J为转动惯量。
5.如权利要求1所述的北斗三号卫星处于地影期偏航角正负阶跃时预报姿态的方法，其特征在于，计算下一时刻角速度，偏航角包括：
下一状态角速度为
ωi+1＝ωi+αi*T (5)
下一状态偏航角
ψi+1＝ψi+ωi*T+0.5*αi*T2 (6)。
6.如权利要求1所述的北斗三号卫星处于地影期时预报姿态的方法，其特征在于，循环迭代计算每一时刻偏航角步骤包括：
将下一时刻偏航角和偏航角速度迭代至当前时刻，并根据所述步骤计算其下一时刻偏航角和偏航角速度。
之后按式(1)至(6)进行下一步计算，循环(1)-(7)直至偏航角正负阶跃完成。
7.如权利要求1所述的北斗三号卫星处于地影期偏航角正负阶跃时预报姿态的方法，其特征在于，利用北斗三号MEO卫星处于地影期偏航角正负阶跃时预报姿态的方法，以将预报姿态用于星间链路指向计算；
利用北斗三号卫星处于地影期偏航角正负阶跃时预报姿态的方法，以解决现有预报姿态算法在地影期偏航正负阶跃时预报误差极大。
8.如权利要求7所述的北斗三号卫星处于地影期偏航角正负阶跃时预报姿态的方法，其特征在于，星间链路指向计算包括：
根据卫星在地影期时轨道信息和偏航设置，计算北斗三号卫星的当前姿态信息；以及
将北斗三号卫星的当前姿态信息用于星间链路指向，以替代或更正由于星内姿态信息传输失败引起的建链指向错误。
9.如权利要求1所述的北斗三号卫星处于地影期偏航角正负阶跃预报姿态的方法，其特征在于，还包括：
在地面站无法获取北斗三号遥测时，包括由于境外遥测中断、星间链路误码、星地链路误码情况导致的姿态信息获取失败，用所述方法根据卫星位置信息计算相应时刻的姿态信息以用于地面事后或事前计算，包括定轨、光压估计、指向仿真、卫星建模；
在无权获取北斗三号遥测的单位，包括工厂、软件中，用所述方法进行姿态计算，以完成其所需。
10.如权利要求1所述的北斗三号卫星处于地影期偏航角正负阶跃预报姿态的方法，其特征在于，还包括：
当轨道太阳角β处于[-3°,+3°]时，卫星处于深地影区。在一年内，北斗三号MEO卫星轨道太阳角经历两次符号变换，即两段地影期。此时卫星设置太阳高度角为
则卫星偏航角为
当β正负切换时，偏航角将产生巨大的符号变换，这将给卫星的姿控伺服系统造成巨大的负荷压力，是姿控伺服机构承受不了的。所以在卫星姿控系统设置中，除了轨道太阳角度外，连续动态偏航还需考虑上一时刻偏航角与这一时刻偏航角尽量平稳变化，将偏航角正负切换延时到±5°执行。
此时，真正的偏航角需在式(9)的基础上增加惯性处理，具体计算方式为
因此，卫星的偏航角将产生从+5°到-5°(或从-5°到+5°)的阶跃变化。此即为本发明所述方法的进入条件。同时式(8)-(10)也是地影期偏航角未发生阶跃时期的预报方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110920677.7A CN113485095B (zh) | 2021-08-11 | 2021-08-11 | 北斗三号卫星处于地影期时预报姿态的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110920677.7A CN113485095B (zh) | 2021-08-11 | 2021-08-11 | 北斗三号卫星处于地影期时预报姿态的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113485095A true CN113485095A (zh) | 2021-10-08 |
| CN113485095B CN113485095B (zh) | 2022-09-13 |
Family
ID=77946255
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110920677.7A Active CN113485095B (zh) | 2021-08-11 | 2021-08-11 | 北斗三号卫星处于地影期时预报姿态的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113485095B (zh) |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE3032571A1 (de) * | 1980-08-29 | 1982-03-11 | Messerschmitt-Bölkow-Blohm GmbH, 8000 München | Verfahren zum betrieb eines satelliten ohne batterie |
| EP0174715A2 (en) * | 1984-09-13 | 1986-03-19 | Mitsubishi Denki Kabushiki Kaisha | Attitude angle calculation apparatus for a geostationary satellite |
| US5556058A (en) * | 1994-05-16 | 1996-09-17 | Hughes Electronics | Spacecraft attitude determination using sun sensor, earth sensor, and space-to-ground link |
| DE69606275D1 (de) * | 1995-08-11 | 2000-02-24 | Daimler Chrysler Ag | Verfahren zur Lageregelung eines dreiachsenstabilisierten, erdorientierten trägheitsmomentstabilisierten Raumfahrzeuges |
| CN102591349A (zh) * | 2012-03-12 | 2012-07-18 | 北京控制工程研究所 | 高轨道卫星大初始角速率情况的无陀螺太阳捕获控制方法 |
| CN104181930A (zh) * | 2014-09-02 | 2014-12-03 | 上海新跃仪表厂 | 一种倾斜轨道卫星偏航机动的自主控制方法 |
| CN104267735A (zh) * | 2014-09-02 | 2015-01-07 | 上海新跃仪表厂 | 一种倾斜轨道卫星偏航机动时机判断方法 |
| CN108710379A (zh) * | 2018-06-14 | 2018-10-26 | 上海卫星工程研究所 | 静止卫星成像偏航导引角计算方法 |
| CN108803642A (zh) * | 2018-06-14 | 2018-11-13 | 上海卫星工程研究所 | 光学成像卫星相机的太阳保护姿控关联设计方法 |
| CN112014869A (zh) * | 2020-08-12 | 2020-12-01 | 中国科学院微小卫星创新研究院 | 基于天文导航的星间链路自主导航方法及系统 |
| CN112213747A (zh) * | 2020-09-29 | 2021-01-12 | 中国科学院微小卫星创新研究院 | 粗轨航天器与北斗三号卫星星间链路信号双向捕获方法 |
| CN112632756A (zh) * | 2020-12-07 | 2021-04-09 | 上海卫星工程研究所 | 基于太阳敏感器的卫星地影自主预报方法及系统 |
| CN112710298A (zh) * | 2020-12-02 | 2021-04-27 | 惠州学院 | 基于动力学模型辅助的旋转弹用地磁卫星组合导航方法 |
-
2021
- 2021-08-11 CN CN202110920677.7A patent/CN113485095B/zh active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE3032571A1 (de) * | 1980-08-29 | 1982-03-11 | Messerschmitt-Bölkow-Blohm GmbH, 8000 München | Verfahren zum betrieb eines satelliten ohne batterie |
| EP0174715A2 (en) * | 1984-09-13 | 1986-03-19 | Mitsubishi Denki Kabushiki Kaisha | Attitude angle calculation apparatus for a geostationary satellite |
| US5556058A (en) * | 1994-05-16 | 1996-09-17 | Hughes Electronics | Spacecraft attitude determination using sun sensor, earth sensor, and space-to-ground link |
| DE69606275D1 (de) * | 1995-08-11 | 2000-02-24 | Daimler Chrysler Ag | Verfahren zur Lageregelung eines dreiachsenstabilisierten, erdorientierten trägheitsmomentstabilisierten Raumfahrzeuges |
| CN102591349A (zh) * | 2012-03-12 | 2012-07-18 | 北京控制工程研究所 | 高轨道卫星大初始角速率情况的无陀螺太阳捕获控制方法 |
| CN104267735A (zh) * | 2014-09-02 | 2015-01-07 | 上海新跃仪表厂 | 一种倾斜轨道卫星偏航机动时机判断方法 |
| CN104181930A (zh) * | 2014-09-02 | 2014-12-03 | 上海新跃仪表厂 | 一种倾斜轨道卫星偏航机动的自主控制方法 |
| CN108710379A (zh) * | 2018-06-14 | 2018-10-26 | 上海卫星工程研究所 | 静止卫星成像偏航导引角计算方法 |
| CN108803642A (zh) * | 2018-06-14 | 2018-11-13 | 上海卫星工程研究所 | 光学成像卫星相机的太阳保护姿控关联设计方法 |
| CN112014869A (zh) * | 2020-08-12 | 2020-12-01 | 中国科学院微小卫星创新研究院 | 基于天文导航的星间链路自主导航方法及系统 |
| CN112213747A (zh) * | 2020-09-29 | 2021-01-12 | 中国科学院微小卫星创新研究院 | 粗轨航天器与北斗三号卫星星间链路信号双向捕获方法 |
| CN112710298A (zh) * | 2020-12-02 | 2021-04-27 | 惠州学院 | 基于动力学模型辅助的旋转弹用地磁卫星组合导航方法 |
| CN112632756A (zh) * | 2020-12-07 | 2021-04-09 | 上海卫星工程研究所 | 基于太阳敏感器的卫星地影自主预报方法及系统 |
Non-Patent Citations (1)
| Title |
|---|
| 陈婷婷等: "基于星间链路的星上时间自主完好性监测方法", 《应用科学学报》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113485095B (zh) | 2022-09-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ahn et al. | Trajectory‐keeping in satellite formation flying via robust periodic learning control | |
| CN104332707A (zh) | 一种用于低轨星载天线跟踪地面站的方法 | |
| CN103279127B (zh) | 一种仅用角度信息的geo轨道卫星自主控制方法 | |
| Roth et al. | Flight results from the CanX-4 and CanX-5 formation flying mission | |
| Klesh et al. | Cyber-physical challenges for space systems | |
| CN112607065B (zh) | 一种基于电推进系统的高精度相位控制方法 | |
| CN103413006A (zh) | 一种空间惯性定向姿态卫星数传天线波束的设计方法 | |
| CN112014869B (zh) | 基于天文导航的星间链路自主导航方法及系统 | |
| Hill et al. | A lunar L2 navigation, communication, and gravity mission | |
| Miotto | Designing and validating proximity operations rendezvous and approach trajectories for the cygnus mission | |
| CN110077627B (zh) | 一种空间激光干涉引力波探测器轨道修正方法及系统 | |
| Scharf et al. | Flight-like ground demonstrations of precision maneuvers for spacecraft formations—Part II | |
| CN112468211B (zh) | 跟踪天线在轨指向误差修正方法及天基测控通信系统 | |
| CN113485095B (zh) | 北斗三号卫星处于地影期时预报姿态的方法 | |
| EP4032816A1 (en) | Manipulation method suitable for non-redundant sgcmg group | |
| CN103514362A (zh) | 基于模型误差补偿的两行根数生成方法 | |
| US20210403182A1 (en) | Fail-Safe Vehicle Rendezvous in Case of Total Control Failure | |
| CN112762925A (zh) | 一种基于地磁计和陀螺仪的低轨卫星定姿方法 | |
| CN102394688B (zh) | 一种天线固联式对地通信方法 | |
| CN109918706B (zh) | 一种基于广义动力学的卫星-天线耦合系统路径规划算法 | |
| CN113624238B (zh) | 星间链路指向恢复方法 | |
| Sun et al. | Modeling and Calibrating the Ground-Surface Beam Pointing of GEO Satellite | |
| CN114413911A (zh) | 掩星探测自主任务姿态导引方法及系统 | |
| CN112013834A (zh) | 基于天文导航的星间链路自主恢复方法及系统 | |
| Jin et al. | Smart Communication Satellite (SCS) Project Overview |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |