CN109460049B - 基于惯性指向模式的地球同步轨道卫星远地点变轨方法 - Google Patents
基于惯性指向模式的地球同步轨道卫星远地点变轨方法 Download PDFInfo
- Publication number
- CN109460049B CN109460049B CN201811354941.XA CN201811354941A CN109460049B CN 109460049 B CN109460049 B CN 109460049B CN 201811354941 A CN201811354941 A CN 201811354941A CN 109460049 B CN109460049 B CN 109460049B
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- satellite
- inertial
- orbit
- transfer matrix
- 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
Images
Landscapes
- Navigation (AREA)
Abstract
基于惯性指向模式的地球同步轨道卫星远地点变轨方法,首先建立J2000地心惯性系、J2000准惯性坐标系、卫星本体坐标系、东南地坐标系、轨道坐标系,然后测量计算得到各个坐标系之间的转换矩阵,并计算地球同步轨道卫星在J2000准惯性系下的三轴目标姿态,最后计算当前三轴目标姿态、目标三轴目标姿态的差,并进行地球同步轨道卫星远点点火期间的姿态控制。本发明利用星敏感器进行远地点点火姿态建立,不要求偏航太阳可见,降低了对发射窗口和入轨精度的要求,提高了点火时机选择的灵活性,具有很好的使用价值。
Description
技术领域
本发明涉及卫星控制领域,特别是基于惯性指向模式(IPM,Inertial PointingMode)的地球同步轨道卫星远地点变轨方法。
背景技术
当前,东三、东四平台等地球同步轨道卫星,火箭把卫星送入转移轨道,需要在星箭分离后进行多次轨道控制。轨道控制任务是卫星飞控工作的重要内容。
在远地点点火前,一般在地球指向模式下利用地球敏感器和数字太阳敏感器标定陀螺,这要求地球、卫星和太阳夹角满足一定要求,从而限制了卫星发射窗口。利用星敏建立远地点姿态对地球、卫星和太阳夹角没有要求,减少了卫星发射窗口约束条件,利用星敏进行远地点姿态建立和点火控制,星上完全可以自主进行,大大减少了当前飞控的工作量,因此需要提出一种利用星敏来完成地球同步轨道卫星远地点变轨的方法。
发明内容
本发明解决的技术问题是:克服现有技术的不足,提供了基于惯性指向模式的地球同步轨道卫星远地点变轨方法,解决了现有技术利用地球敏感器和数字太阳敏感器标定陀螺,要求地球、卫星和太阳夹角满足一定要求的缺陷,通过利用星敏建立远地点姿态对地球、卫星和太阳夹角没有要求,减少了卫星发射窗口约束条件。
本发明的技术解决方案是:基于惯性指向模式的地球同步轨道卫星远地点变轨方法,包括如下步骤:
(1)以地心为坐标原点,历元J2000.0时刻的平赤道面为XY坐标面,右手定则确定Z轴建立J2000地心惯性系OXiYiZi;以卫星质心为坐标系原点,Y轴垂直于J2000赤道平面指向正南,Z平行于J2000地心惯性坐标系的X轴,X轴由右手系确定建立J2000准惯性坐标系;以卫星质心为原点,三轴与卫星机械坐标系的三轴平行,建立卫星本体坐标系OXbYbZb;以卫星质心为原点,Z轴从卫星所处的轨道位置指向地心,X轴沿当地水平面正东方向,Y轴沿当地子午线指向当地正南,Z轴与轨道坐标系Z轴重合建立东南地坐标系OXEYEZE;以卫星质心为原点,Z轴从卫星所处的轨道位置指向地心,Y轴为轨道平面的负法线方向,X轴与Y轴、Z轴构成右手直角坐标系建立轨道坐标系OXoYoZo;
(2)测量计算得到J2000地心惯性系到轨道坐标系的转移矩阵Coi、轨道坐标系到东南地坐标系的转移矩阵Cdo、东南地坐标系到卫星本体坐标系的转移矩阵Cbd、轨道坐标系到卫星本体坐标系的转移矩阵Cbo;设J2000地心惯性坐标系到卫星本体系的转移矩阵为Cbi,J2000准惯性坐标系到卫星本体系的转移矩阵为CbI;
(3)根据J2000地心惯性系到轨道坐标系的转移矩阵Coi、轨道坐标系到东南地坐标系的转移矩阵Cdo、东南地坐标系到卫星本体坐标系的转移矩阵Cbd计算得到J2000地心惯性坐标系到卫星本体系的转移矩阵为Cbi或者根据J2000地心惯性系到轨道坐标系的转移矩阵Coi、轨道坐标系到卫星本体坐标系的转移矩阵Cbo计算得到J2000地心惯性坐标系到卫星本体系的转移矩阵为Cbi;根据轨道坐标系到卫星本体坐标系的转移矩阵Cbo、J2000地心惯性系到轨道坐标系的转移矩阵Coi计算得到J2000准惯性坐标系到卫星本体系的转移矩阵CbI;
计算得到
θbI(t)=θt0+θrate(ta-t0)
ψbI(t)=ψt0+ψrate(ta-t0)
所述的J2000地心惯性系到J2000准惯性系的转换矩阵为
所述的根据轨道坐标系到卫星本体坐标系的转移矩阵Cbo、J2000地心惯性系到轨道坐标系的转移矩阵Coi计算得到J2000准惯性坐标系到卫星本体系的转移矩阵CbI的方法为:
所述的根据J2000地心惯性系到轨道坐标系的转移矩阵Coi、轨道坐标系到东南地坐标系的转移矩阵Cdo、东南地坐标系到卫星本体坐标系的转移矩阵Cbd计算得到J2000地心惯性坐标系到卫星本体系的转移矩阵为Cbi的方法为:
Cbi=CbdCdoCoi。
所述的根据J2000地心惯性系到轨道坐标系的转移矩阵Coi、轨道坐标系到卫星本体坐标系的转移矩阵Cbo计算得到J2000地心惯性坐标系到卫星本体系的转移矩阵为Cbi的方法为
Cbi=CboCoi。
(3)计算姿态变化率为
一种计算机可读存储介质,所述的计算机可读存储介质存储有计算机程序,所述的计算机程序被处理器执行时实现如权利要求1-权利要求6任一所述方法的步骤。
一种基于惯性指向模式的地球同步轨道卫星远地点变轨终端设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,所述的处理器执行所述的计算机程序时实现如权利要求1-权利要求6任一所述方法的步骤。
本发明与现有技术相比的优点在于:
(1)本发明利用星敏感器进行远地点点火姿态建立,不要求偏航太阳可见,降低了对发射窗口和入轨精度的要求,提高了点火时机选择的灵活性;
(2)本发明利用星敏感器,自主实时标定陀螺漂移,相比较传统的利用地敏和太敏数据通过地面计算的方法,既提高了效率,又提高了标定精度;
(3)本发明利用星敏输出的惯性系姿态,实时拟合点火弧段,和传统通过注入等价陀螺漂移参数分段拟合点火弧段相比,提高了弧段拟合的灵活性和拟合精度;
(4)本发明利用星敏星上自主进行远地点姿态建立和点火控制,精简了飞控流程,减少了当前飞控的工作量。
附图说明
图1为本发明坐标系说明;
图2为本发明的流程图;
图3为本发明点火过程中姿态拟合曲线。
具体实施方式
本发明克服现有技术的不足,提供了基于惯性指向模式的地球同步轨道卫星远地点变轨方法,解决了现有技术利用地球敏感器和数字太阳敏感器标定陀螺,要求地球、卫星和太阳夹角满足一定要求的缺陷,通过利用星敏建立远地点姿态对地球、卫星和太阳夹角没有要求,减少了卫星发射窗口约束条件。下面结合附图对本发明方法进行详细说明。如图2所示为本发明的流程图,本发明方法包括如下步骤。
1、坐标系定义
(1)J2000地心惯性系OXiYiZi:坐标原点在地心,XY坐标面是历元J2000.0时刻的平赤道面,X轴方向是该历元的平春分点。
(2)J2000准惯性坐标系OXIYIZI:坐标系原点在卫星质心随卫星运动,Y轴垂直于J2000赤道平面指向正南,Z平行于J2000地心惯性坐标系的X轴,X轴由右手系确定。J2000准惯性坐标系与J2000地心惯性坐标系只有相对平动,没有转动。
按照定义,J2000地心惯性系到J2000准惯性系的转换矩阵为
J2000地心惯性系到J2000准惯性系的变换关系如附图1所示。
(3)卫星本体坐标系OXbYbZb:原点在卫星质心,三轴与卫星机械坐标系的三轴平行。
(4)东南地坐标系OXEYEZE:原点在卫星质心,ZE轴从卫星所处的轨道位置指向地心,XE轴沿当地水平面正东方向,YE轴沿当地子午线指向当地正南,构成右手坐标系。东南地坐标系ZE轴与轨道坐标系Zo轴重合。
(5)轨道坐标系OXoYoZo:原点在卫星质心,Zo轴从卫星所处的轨道位置指向地心,Y轴为轨道平面的负法线方向,X轴在轨道平面内,与Y轴、Z轴构成右手直角坐标系。
2、惯性指向模式建立远地点姿态控制系数计算
(a)已知东南地坐标系下的卫星点火姿态,计算IPM远地点姿态
东三、东四平台的卫星习惯在东南地坐标系描述远地点姿态,为了在使用上兼容,本文给出如何从东南系的点火姿态计算IPM所需的姿态偏置量。
定义:
Coi为J2000地心惯性系到轨道坐标系的转移矩阵,根据测轨结果和轨道动力学可得;
Cdo为轨道坐标系到东南地坐标系的转移矩阵,根据测轨结果和轨道动力学可得;
Cbd为东南地坐标系到卫星本体坐标系的转移矩阵,根据点火策略可得;
Cbo为轨道坐标系到卫星本体坐标系的转移矩阵,根据点火策略可得;
Cbi为J2000地心惯性坐标系到卫星本体系的转移矩阵,由Coi、Cdo和Cbd计算获得;
已知Cbd或Cbo、Cdo和Coi,可以求得Cbi:
根据Cbi=CbICIi和CIi,可求得CbI
可以计算出三个姿态目标值:
卫星在J2000准惯性系下,θbI,ψbI目标值和卫星在轨道上的位置相关,也就是是时间t的函数,改定义为θbI(t),ψbI(t)。点火弧段较长时,考虑点火弧段的影响,点火过程中,星上采用一阶线性拟合该弧段,变化规律为:
θbI(t)=θt0+θrate(ta-t0)
ψbI(t)=ψt0+ψrate(ta-t0)
其中,θt0,ψt0为点火开始时刻J2000准惯性系下三轴姿态,θrate,ψrate为点火期间J2000准惯性系三轴姿态变化率,t0为点火开始时刻,点火结束时刻为t1。7个参数为需要地面计算后遥控注入星上。
t0、t1时刻的轨道和姿态偏置已知,系数计算步骤如下:
(3)计算姿态变化率
(b)已知轨道坐标系下的卫星点火姿态,计算IPM远地点姿态
和已知东南地坐标系下的卫星点火姿态相比,就计算Cbi使用的公式不同,其他计算公式和方法完全相同。
t0、t1时刻的轨道和姿态偏置已知,系数计算步骤如下:
(3)计算姿态变化率
3、远地点点火期间姿态控制
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (8)
1.基于惯性指向模式的地球同步轨道卫星远地点变轨方法,其特征在于包括如下步骤:
(1)以地心为坐标原点,历元J2000.0时刻的平赤道面为XY坐标面,右手定则确定Z轴建立J2000地心惯性系OXiYiZi;以卫星质心为坐标系原点,Y轴垂直于J2000赤道平面指向正南,Z平行于J2000地心惯性坐标系的X轴,X轴由右手系确定建立J2000准惯性坐标系;以卫星质心为原点,三轴与卫星机械坐标系的三轴平行,建立卫星本体坐标系OXbYbZb;以卫星质心为原点,Z轴从卫星所处的轨道位置指向地心,X轴沿当地水平面正东方向,Y轴沿当地子午线指向当地正南,Z轴与轨道坐标系Z轴重合建立东南地坐标系OXEYEZE;以卫星质心为原点,Z轴从卫星所处的轨道位置指向地心,Y轴为轨道平面的负法线方向,X轴与Y轴、Z轴构成右手直角坐标系建立轨道坐标系OXoYoZo;
(2)测量计算得到J2000地心惯性系到轨道坐标系的转移矩阵Coi、轨道坐标系到东南地坐标系的转移矩阵Cdo、东南地坐标系到卫星本体坐标系的转移矩阵Cbd、轨道坐标系到卫星本体坐标系的转移矩阵Cbo;设J2000地心惯性坐标系到卫星本体坐标系的转移矩阵为Cbi,J2000准惯性坐标系到卫星本体坐标系的转移矩阵为CbI;
(3)根据J2000地心惯性系到轨道坐标系的转移矩阵Coi、轨道坐标系到东南地坐标系的转移矩阵Cdo、东南地坐标系到卫星本体坐标系的转移矩阵Cbd计算得到J2000地心惯性坐标系到卫星本体坐标系的转移矩阵为Cbi或者根据J2000地心惯性系到轨道坐标系的转移矩阵Coi、轨道坐标系到卫星本体坐标系的转移矩阵Cbo计算得到J2000地心惯性坐标系到卫星本体坐标系的转移矩阵为Cbi;根据轨道坐标系到卫星本体坐标系的转移矩阵Cbo、J2000地心惯性系到轨道坐标系的转移矩阵Coi计算得到J2000准惯性坐标系到卫星本体坐标系的转移矩阵CbI;
计算得到
θbI(t)=θt0+θrate(ta-t0)
ψbI(t)=ψt0+ψrate(ta-t0)
4.根据权利要求1所述的基于惯性指向模式的地球同步轨道卫星远地点变轨方法,其特征在于:所述的根据J2000地心惯性系到轨道坐标系的转移矩阵Coi、轨道坐标系到东南地坐标系的转移矩阵Cdo、东南地坐标系到卫星本体坐标系的转移矩阵Cbd计算得到J2000地心惯性坐标系到卫星本体坐标系的转移矩阵为Cbi的方法为:
Cbi=CbdCdoCoi。
5.根据权利要求1所述的基于惯性指向模式的地球同步轨道卫星远地点变轨方法,其特征在于:所述的根据J2000地心惯性系到轨道坐标系的转移矩阵Coi、轨道坐标系到卫星本体坐标系的转移矩阵Cbo计算得到J2000地心惯性坐标系到卫星本体坐标系的转移矩阵为Cbi的方法为
Cbi=CboCoi。
7.一种计算机可读存储介质,所述的计算机可读存储介质存储有计算机程序,其特征在于,所述的计算机程序被处理器执行时实现如权利要求1-权利要求6任一所述方法的步骤。
8.一种基于惯性指向模式的地球同步轨道卫星远地点变轨终端设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其特征在于:所述的处理器执行所述的计算机程序时实现如权利要求1-权利要求6任一所述方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811354941.XA CN109460049B (zh) | 2018-11-14 | 2018-11-14 | 基于惯性指向模式的地球同步轨道卫星远地点变轨方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811354941.XA CN109460049B (zh) | 2018-11-14 | 2018-11-14 | 基于惯性指向模式的地球同步轨道卫星远地点变轨方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109460049A CN109460049A (zh) | 2019-03-12 |
CN109460049B true CN109460049B (zh) | 2021-11-16 |
Family
ID=65610450
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811354941.XA Active CN109460049B (zh) | 2018-11-14 | 2018-11-14 | 基于惯性指向模式的地球同步轨道卫星远地点变轨方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109460049B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111427000A (zh) * | 2020-03-19 | 2020-07-17 | 上海卫星工程研究所 | 适用于地面测站天线对卫星指向的目标视向量确定方法 |
CN112357121B (zh) * | 2020-10-30 | 2022-04-29 | 哈尔滨工业大学 | 一种基于预置姿态的卫星入轨快速成像方法 |
CN113465570B (zh) * | 2021-06-17 | 2023-05-12 | 北京控制工程研究所 | 一种基于高精度imu的气浮台初始对准方法及系统 |
CN113761664B (zh) * | 2021-09-15 | 2024-01-30 | 长光卫星技术股份有限公司 | 一种遥感卫星成像期间飞轮转速优化方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4943014A (en) * | 1979-10-30 | 1990-07-24 | Trw Inc. | Soft ride method for changing the altitude or position of a spacecraft in orbit |
US5716029A (en) * | 1994-03-25 | 1998-02-10 | Hughes Electronics | Constant sun angle transfer orbit sequence and method using electric propulsion |
EP0640524B1 (en) * | 1993-08-24 | 1999-10-13 | Orbital Sciences Corporation | Method for injecting payloads into orbit |
CN1239458A (zh) * | 1997-09-12 | 1999-12-22 | 国家航空工业公司 | 把人造卫星安置到对地静止轨道上的方法 |
CN101186236A (zh) * | 2007-12-26 | 2008-05-28 | 北京控制工程研究所 | 一种减少航天器重力损失的变轨方法 |
CN102424116A (zh) * | 2011-12-08 | 2012-04-25 | 中国空间技术研究院 | 一种静止轨道卫星变轨策略优化方法 |
CN102880184A (zh) * | 2012-10-24 | 2013-01-16 | 北京控制工程研究所 | 一种静止轨道卫星自主轨道控制方法 |
CN105334731A (zh) * | 2015-10-30 | 2016-02-17 | 北京控制工程研究所 | 一种远地点点火高精度解析轨道自主预报方法 |
CN107826269A (zh) * | 2017-09-18 | 2018-03-23 | 北京控制工程研究所 | 一种适用于地球静止轨道卫星平台的近地点变轨方法 |
CN107985629A (zh) * | 2017-10-17 | 2018-05-04 | 上海宇航系统工程研究所 | 火星探测方法及火星探测器 |
CN108791957A (zh) * | 2018-06-14 | 2018-11-13 | 上海卫星工程研究所 | 一种卫星变轨期间横向质心高精度实时辨识方法 |
-
2018
- 2018-11-14 CN CN201811354941.XA patent/CN109460049B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4943014A (en) * | 1979-10-30 | 1990-07-24 | Trw Inc. | Soft ride method for changing the altitude or position of a spacecraft in orbit |
EP0640524B1 (en) * | 1993-08-24 | 1999-10-13 | Orbital Sciences Corporation | Method for injecting payloads into orbit |
US5716029A (en) * | 1994-03-25 | 1998-02-10 | Hughes Electronics | Constant sun angle transfer orbit sequence and method using electric propulsion |
CN1239458A (zh) * | 1997-09-12 | 1999-12-22 | 国家航空工业公司 | 把人造卫星安置到对地静止轨道上的方法 |
CN101186236A (zh) * | 2007-12-26 | 2008-05-28 | 北京控制工程研究所 | 一种减少航天器重力损失的变轨方法 |
CN102424116A (zh) * | 2011-12-08 | 2012-04-25 | 中国空间技术研究院 | 一种静止轨道卫星变轨策略优化方法 |
CN102880184A (zh) * | 2012-10-24 | 2013-01-16 | 北京控制工程研究所 | 一种静止轨道卫星自主轨道控制方法 |
CN105334731A (zh) * | 2015-10-30 | 2016-02-17 | 北京控制工程研究所 | 一种远地点点火高精度解析轨道自主预报方法 |
CN107826269A (zh) * | 2017-09-18 | 2018-03-23 | 北京控制工程研究所 | 一种适用于地球静止轨道卫星平台的近地点变轨方法 |
CN107985629A (zh) * | 2017-10-17 | 2018-05-04 | 上海宇航系统工程研究所 | 火星探测方法及火星探测器 |
CN108791957A (zh) * | 2018-06-14 | 2018-11-13 | 上海卫星工程研究所 | 一种卫星变轨期间横向质心高精度实时辨识方法 |
Non-Patent Citations (2)
Title |
---|
基于冲量变轨原理的地球同步卫星有限推力变轨策略;彭坤 等;《空间控制技术与应用》;20081231;第34卷(第6期);第48-51,57页 * |
高分四号卫星控制方案设计特点及在轨实现;黎康 等;《空间控制技术与应用》;20161231;第42卷(第6期);第1-8页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109460049A (zh) | 2019-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109460049B (zh) | 基于惯性指向模式的地球同步轨道卫星远地点变轨方法 | |
US9073648B2 (en) | Star tracker rate estimation with kalman filter enhancement | |
WO2022007437A1 (zh) | 传感器安装偏差角的标定方法、组合定位系统和车辆 | |
US9643740B2 (en) | Gyroless three-axis sun acquisition using sun sensor and unscented kalman filter | |
EP3346234A1 (en) | Autonomous navigation method for mars probe | |
US11142351B2 (en) | Earth satellite attitude data fusion system and method thereof | |
CN110435930B (zh) | 低轨光学卫星匀降速推扫姿态规划方法 | |
CN108663052B (zh) | 一种星上自主空间非合作目标相对导航相机指向控制方法 | |
Wolf et al. | Toward improved landing precision on Mars | |
CN110096726B (zh) | 基于月球借力的geo卫星应急转移轨道快速优化设计方法 | |
CN103063216A (zh) | 一种基于星像坐标建模的惯性与天文组合导航方法 | |
CN116105730A (zh) | 基于合作目标卫星甚短弧观测的仅测角光学组合导航方法 | |
CN114993305A (zh) | 一种基于发射坐标系的制导炮弹组合导航方法 | |
CN109407696B (zh) | 一种无人机航向角动态校定方法 | |
CN114879709A (zh) | 一种面向运动目标跟踪观测的卫星姿态控制方法及装置 | |
CN112985391B (zh) | 一种基于惯性和双目视觉的多无人机协同导航方法和装置 | |
CN108981691B (zh) | 一种天空偏振光组合导航在线滤波与平滑方法 | |
CN112945242B (zh) | 一种卫星在轨自主规划任务最佳时间与姿态的方法 | |
CN107883925B (zh) | 一种导航星座星间观测目标卫星图像模拟方法 | |
CN112525204A (zh) | 一种航天器惯性和太阳多普勒速度组合导航方法 | |
CN106250684B (zh) | 基于地固系数据的卫星过境时间快速计算方法 | |
CN115392540A (zh) | 一种用于月球轨道交会制导的快速预报方法 | |
CN112329202B (zh) | 一种火星车对环绕器天线指向算法的优化实现方法 | |
CN113447025A (zh) | 基于克雷洛夫角的惯性导航高精度姿态角解算方法和系统 | |
KR101510435B1 (ko) | 단일차분을 이용한 3축 각속도센서의 오정렬 및 축척계수 오차 보정 방법 |
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 |