CN106005482B - 一种适用于导航倾斜轨道卫星的零偏持续天数确定方法 - Google Patents
一种适用于导航倾斜轨道卫星的零偏持续天数确定方法 Download PDFInfo
- Publication number
- CN106005482B CN106005482B CN201610344877.1A CN201610344877A CN106005482B CN 106005482 B CN106005482 B CN 106005482B CN 201610344877 A CN201610344877 A CN 201610344877A CN 106005482 B CN106005482 B CN 106005482B
- Authority
- CN
- China
- Prior art keywords
- satellite
- plane
- ecliptic
- zero
- orbit
- 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.)
- Active
Links
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Radar, Positioning & Navigation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Navigation (AREA)
Abstract
一种适用于导航倾斜轨道卫星的零偏持续天数确定方法,根据所述导航倾斜轨道卫星的轨道升交点赤经和倾角计算得到黄道面与轨道面的夹角,根据黄道面与轨道面的夹角、动偏转零偏条件中的太阳矢量与轨道面夹角得到卫星零偏期间太阳在黄道上的弧段;根据卫星零偏期间太阳在黄道上的弧段确定卫星零偏持续天数。本发明通过几何分析,根据卫星的两个轨道根数在轨道设计期间和设计完成后均可计算卫星零偏置的持续天数。该方法不仅计算简单,需要的轨道信息少,同时也可以反向应用于零偏置持续天数约束下的轨道设计。
Description
技术领域
本发明涉及一种卫星零偏持续天数的确定方法,适用于偏航姿态控制采用动态偏置和零偏置相结合策略的导航倾斜轨道卫星。
背景技术
导航IGSO\MEO卫星轨道为倾角为55°的倾斜轨道,为满足整星供电要求,需要进行偏航姿态连续测量与控制,即偏航动态偏置。但当太阳矢量与轨道面的夹角较小时,偏航角速度较大,需要反作用轮提供的偏航控制力矩也较大。由于轮子控制力矩有限,在偏航快速变化期间控制精度将降低,进而影响帆板跟踪太阳的精度要求,为此在太阳矢量与轨道面的夹角较小时采用偏航零偏置。整星偏航姿态控制采用动态偏置和零偏置相结合的策略。
对于导航定位系统,倾斜轨道卫星的姿态控制方法不同,精密定轨的太阳光压动力学模型也将不同,直接影响卫星的定轨精度。目前,国内定轨软件的动力学模型适用于动态偏置姿态,在零偏置期间的定轨精度会降低。因此,分析导航倾斜轨道卫星零偏置的持续天数可以用于分析姿态控制对定轨影响的持续时间。
目前,已有的导航倾斜轨道卫星零偏置的持续天数计算方法是根据卫星的六个轨道根数(表征了轨道在空间的位置和卫星在轨道上的位置)通过数值计算,得到卫星零偏置的持续天数,即在卫星轨道设计完成后,卫星在轨道上的位置完全确定后方可计算得到,因此,原有方法计算复杂,需要的轨道信息多,且计算中需要进行轨道外推。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提供一种适用于导航倾斜轨道卫星的零偏持续天数确定方法,该方法不仅计算简单,需要的轨道信息少,同时不需要进行轨道外推,可以反向应用于零偏置持续天数约束下的轨道设计。
本发明的技术解决方案是:一种适用于导航倾斜轨道卫星的零偏持续天数确定方法,步骤如下:
(1)根据导航倾斜轨道卫星的轨道升交点赤经和倾角计算黄道面与轨道面的夹角;
(2)根据黄道面与轨道面的夹角以及动偏转零偏条件中的太阳矢量与轨道面夹角得到卫星零偏期间太阳在黄道上的弧段;
(3)根据卫星零偏期间太阳在黄道上的弧段确定卫星零偏持续天数。
所述步骤(1)中黄道面与轨道面的夹角β'依据下式确定:
其中,Ω为卫星轨道升交点赤经,i为卫星轨道倾角,ε是黄道面与赤道面的夹角,β'为黄道面与轨道面的夹角,所述卫星轨道升交点赤经和卫星轨道倾角为预设值。
所述步骤(2)中卫星零偏期间太阳在黄道上的弧段us满足如下条件:
us=asin(sinus),us>=0
其中,ρ为动偏转零偏条件中的太阳矢量与轨道面夹角,为预设值,β'为黄道面与轨道面的夹角。
所述步骤(3)中确定卫星零偏持续天数T的公式如下:
其中,us为卫星零偏期间太阳在黄道上的弧段,ns为太阳运动平均角速度。
本发明与现有技术相比的优点在于:
(1)现有数值方法需要获取轨道全部信息即六个轨道根数,才可以进行轨道外推,计算零偏置的天数;本发明通过球面几何的分析方法,仅需要卫星的两个轨道根数(倾角和升交点赤经)即可确定零偏持续天数,计算简单,不需要轨道外推。
(2)现有数值方法由于需要获得全部轨道信息,因此通常为轨道设计已经完成后才可以使用。本发明由于需要的轨道信息少,所以可在卫星轨道设计过程中和设计完成后均可计算卫星零偏置的持续天数。
(3)现有数值方法通过轨道外推后,计算得到不同时刻太阳与轨道面的夹角,进而计算太阳与轨道面的夹角小于某个值的时间长度,该时间长度即为零偏置的天数。因此现有数值方法在已知零偏置持续天数时,无法反向计算轨道应该满足的条件,不可以用于轨道设计。本发明可以在已知零偏置持续天数时,按照步骤(3)得到卫星零偏期间太阳在黄道上的弧段,据此依照步骤(2)中的公式(2)可以得到黄道面与轨道面的夹角,再按照步骤(1)中的公式(1)计算得到轨道倾角和升交点赤经需满足的条件。因此,本方法可以反向应用于零偏持续天数约束下的轨道参数设计。
附图说明
图1为本发明方法流程图;
图2为本发明的零偏持续天数分析示意图;其中图2中(a)为0°≤Ω<180°的零偏持续天数分析示意图,图2中(b)为180°≤Ω<360°的零偏持续天数分析示意图。
具体实施方式
本发明提出一种适用于导航倾斜轨道卫星的零偏持续天数确定方法,采用几何分析法,根据卫星的两个轨道根数(倾角和升交点赤经),在轨道设计期间和设计完成后均可确定卫星零偏置的持续天数。该方法不仅计算简单,需要的轨道信息少,同时也可以反向应用于零偏置持续天数约束下的轨道设计。
本发明步骤如下:
(1)如图2所示,根据所述导航倾斜轨道卫星的轨道升交点赤经和倾角计算得到黄道面与轨道面的夹角β'。所述轨道升交点赤经和倾角为预设值。具体计算方法为:
其中,Ω为卫星轨道升交点赤经;i为卫星轨道倾角;ε是黄道面与赤道面的夹角;β'为黄道面与轨道面的夹角。其中图2中(a)为0°≤Ω<180°的示意图,图2中(b)为180°≤Ω<360°的示意图。
(2)根据黄道面与轨道面的夹角、动偏转零偏条件中的太阳矢量与轨道面夹角得到卫星零偏期间太阳在黄道上的弧段。所述动偏转零偏条件中的太阳矢量与轨道面夹角为预设值。具体为:
其中,ρ为动偏转零偏条件中的太阳矢量与轨道面夹角;β'为黄道面与轨道面的夹角;us为卫星零偏期间太阳在黄道上的弧段。
(3)根据卫星零偏期间太阳在黄道上的弧段得到卫星零偏持续天数。具体为:
其中,us为卫星零偏期间太阳在黄道上的弧段;ns为太阳运动平均角速度;T为卫星零偏持续天数。
本方法通过几何分析,仅需要卫星的两个轨道根数(表征轨道在空间的位置),在卫星轨道设计期间和设计完成后均可计算卫星零偏置的持续天数,并且计算中不需要进行轨道外推。因此,该方法不仅计算简单,需要的轨道信息少,同时也可以反向应用于零偏置持续天数约束下的轨道设计。该方法得到的零偏持续天数,一方面可以用于地面应用系统分析零偏置姿态控制对定轨精度影响的持续时间,以及用于卫星可用性的辅助分析;另一方面,当地面应用系统反向要求零偏置姿态对定轨精度影响的持续时间小于某个值时,也可以采用本方法中步骤(3)(2)(1)进行轨道参数的设计。
实施例:
针对某型号导航卫星,预设卫星轨道升交点赤经Ω=90°,卫星轨道倾角i=55°,动偏转零偏条件中的太阳矢量与轨道面夹角ρ=5°。
黄道面与赤道面的夹角以及太阳运动平均角速度为常值,即黄道面与赤道面的夹角ε=23.5°,太阳运动平均角速度ns=0.986°/day。
根据步骤(1)计算黄道面与轨道面的夹角:
cosβ'=cos(23.5°)cos(55°)+sin(23.5°)sin(55°)cos(90°)=0.526;
β'=acos(0.526)=58.264°
根据步骤(2)计算卫星零偏期间太阳在黄道上的弧段:
sinus=sin(5°)/sin(58.264°)=0.1025;
us=asin(0.1025)=5.882°
根据步骤(3)确定卫星零偏持续天数:
T=2*5.882°/(0.986°/day)=11.931day
根据得到的零偏置持续天数对单颗卫星的可用性影响进行分析,以1年为考察周期,卫星将经历2个阶段的零偏置,1个阶段的零偏置时间为11.931day,2个阶段的零偏置时间共计23.862day,考虑最恶劣情况,零偏置期间地面定轨精度均不满足要求,卫星均不可用,则1年中由零偏置引起的单颗卫星不可用时间约为23.862day,卫星可用性约为23.862/365=6.538%。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (4)
1.一种适用于导航倾斜轨道卫星的零偏持续天数确定方法,其特征在于步骤如下:
(1)根据导航倾斜轨道卫星的轨道升交点赤经和倾角计算黄道面与轨道面的夹角;
(2)根据黄道面与轨道面的夹角以及动偏转零偏条件中的太阳矢量与轨道面夹角得到卫星零偏期间太阳在黄道上的弧段;
(3)根据卫星零偏期间太阳在黄道上的弧段确定卫星零偏持续天数。
2.根据权利要求1所述的一种适用于导航倾斜轨道卫星的零偏持续天数确定方法,其特征在于:所述步骤(1)中黄道面与轨道面的夹角β'依据下式确定:
β'=arccos(cosβ'),0°=<β'<=90°
其中,Ω为卫星轨道升交点赤经,i为卫星轨道倾角,ε是黄道面与赤道面的夹角,β'为黄道面与轨道面的夹角,所述卫星轨道升交点赤经和卫星轨道倾角为预设值。
3.根据权利要求1所述的一种适用于导航倾斜轨道卫星的零偏持续天数确定方法,其特征在于:所述步骤(2)中卫星零偏期间太阳在黄道上的弧段us满足如下条件:
us=arcsin(sinus),us>=0
<mrow>
<mi>sin</mi>
<mi> </mi>
<msub>
<mi>u</mi>
<mi>s</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&rho;</mi>
</mrow>
<mrow>
<msup>
<mi>sin&beta;</mi>
<mo>&prime;</mo>
</msup>
</mrow>
</mfrac>
</mrow>
其中,ρ为动偏转零偏条件中的太阳矢量与轨道面夹角,为预设值,β'为黄道面与轨道面的夹角。
4.根据权利要求1所述的一种适用于导航倾斜轨道卫星的零偏持续天数确定方法,其特征在于:所述步骤(3)中确定卫星零偏持续天数T的公式如下:
<mrow>
<mi>T</mi>
<mo>=</mo>
<mn>2</mn>
<mfrac>
<msub>
<mi>u</mi>
<mi>s</mi>
</msub>
<msub>
<mi>n</mi>
<mi>s</mi>
</msub>
</mfrac>
</mrow>
其中,us为卫星零偏期间太阳在黄道上的弧段,ns为太阳运动平均角速度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610344877.1A CN106005482B (zh) | 2016-05-23 | 2016-05-23 | 一种适用于导航倾斜轨道卫星的零偏持续天数确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610344877.1A CN106005482B (zh) | 2016-05-23 | 2016-05-23 | 一种适用于导航倾斜轨道卫星的零偏持续天数确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106005482A CN106005482A (zh) | 2016-10-12 |
CN106005482B true CN106005482B (zh) | 2018-05-01 |
Family
ID=57095824
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610344877.1A Active CN106005482B (zh) | 2016-05-23 | 2016-05-23 | 一种适用于导航倾斜轨道卫星的零偏持续天数确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106005482B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108657467B (zh) * | 2018-05-16 | 2019-05-21 | 清华大学 | 一种利用虚拟太阳矢量的航天器偏航机动控制方法及系统 |
CN109552670B (zh) * | 2018-12-03 | 2021-11-02 | 中科星图(西安)测控技术有限公司 | 一种小推力控制在地球静止同步卫星轨道倾角保持中的应用 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5794891A (en) * | 1995-01-06 | 1998-08-18 | Matra Marconi Space France | Method of controlling the attitude control for satellites on an orbit inclined relative to the equator |
JP3947121B2 (ja) * | 2003-03-20 | 2007-07-18 | 富士通株式会社 | 衛星投入方法,投入ルート提案装置及び投入ルート提案プログラム |
CN101786505A (zh) * | 2010-02-26 | 2010-07-28 | 哈尔滨工业大学 | 基于气动力姿态控制的低轨飞行器及其姿态控制方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101949274B1 (ko) * | 2012-08-24 | 2019-02-19 | 한국전자통신연구원 | 정지 궤도 위성의 궤도 결정 장치 및 방법 |
-
2016
- 2016-05-23 CN CN201610344877.1A patent/CN106005482B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5794891A (en) * | 1995-01-06 | 1998-08-18 | Matra Marconi Space France | Method of controlling the attitude control for satellites on an orbit inclined relative to the equator |
JP3947121B2 (ja) * | 2003-03-20 | 2007-07-18 | 富士通株式会社 | 衛星投入方法,投入ルート提案装置及び投入ルート提案プログラム |
CN101786505A (zh) * | 2010-02-26 | 2010-07-28 | 哈尔滨工业大学 | 基于气动力姿态控制的低轨飞行器及其姿态控制方法 |
Non-Patent Citations (1)
Title |
---|
《姿态模式切换期间QZSS卫星轨道及其钟差产品特性分析》;周佩元等;《测绘学报》;20160331;第45卷(第3期);第275-281页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106005482A (zh) | 2016-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103983254B (zh) | 一种新型敏捷卫星机动中成像方法 | |
CN103917850B (zh) | 一种惯性导航系统的运动对准方法 | |
CN101788296B (zh) | 一种sins/cns深组合导航系统及其实现方法 | |
CN103235328B (zh) | 一种gnss与mems组合导航的方法 | |
WO2017113567A1 (zh) | 火星探测器自主导航方法 | |
CN105160125B (zh) | 一种星敏感器四元数的仿真分析方法 | |
CN100476360C (zh) | 一种基于星敏感器标定的深综合组合导航方法 | |
CN105184002B (zh) | 一种数传天线指向角度的仿真分析方法 | |
CN104567930A (zh) | 一种能够估计和补偿机翼挠曲变形的传递对准方法 | |
CN104792340A (zh) | 一种星敏感器安装误差矩阵与导航系统星地联合标定与校正的方法 | |
CN104848860A (zh) | 一种敏捷卫星成像过程姿态机动规划方法 | |
CN102426025B (zh) | 遥感卫星姿态机动时偏流修正角的仿真分析方法 | |
CN107085634B (zh) | 快速计算太阳光与太阳同步卫星星敏感器最小夹角的方法 | |
CN102944241B (zh) | 基于多胞型线性微分包含的航天器相对姿态确定方法 | |
CN104406583B (zh) | 双星敏感器联合确定载体姿态方法 | |
CN108303079B (zh) | 一种水下usbl反向应用的数据平滑方法 | |
CN109489661B (zh) | 一种卫星初始入轨时陀螺组合常值漂移估计方法 | |
CN101214860A (zh) | 轨控过程中自主选取定姿方式的方法 | |
CN102495950A (zh) | 一种适用于太阳同步轨道的倾角偏置量获取方法 | |
CN109708663B (zh) | 基于空天飞机sins辅助的星敏感器在线标定方法 | |
CN103968834A (zh) | 一种近地停泊轨道上深空探测器的自主天文导航方法 | |
CN106005482B (zh) | 一种适用于导航倾斜轨道卫星的零偏持续天数确定方法 | |
CN103968844A (zh) | 基于低轨平台跟踪测量的大椭圆机动航天器自主导航方法 | |
CN112649006A (zh) | 一种太阳同步圆轨道的轨道规划方法 | |
CN116105730A (zh) | 基于合作目标卫星甚短弧观测的仅测角光学组合导航方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |