CN102878995B - 一种静止轨道卫星自主导航方法 - Google Patents
一种静止轨道卫星自主导航方法 Download PDFInfo
- Publication number
- CN102878995B CN102878995B CN201210413994.0A CN201210413994A CN102878995B CN 102878995 B CN102878995 B CN 102878995B CN 201210413994 A CN201210413994 A CN 201210413994A CN 102878995 B CN102878995 B CN 102878995B
- Authority
- CN
- China
- Prior art keywords
- satellite
- orbit
- cos
- sin
- earth
- 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
Landscapes
- Navigation (AREA)
Abstract
本发明涉及一种静止轨道卫星自主导航方法,属于卫星自主导航研究领域。基于星敏感器和地球敏感器的kalman滤波算法,实时得到卫星相对于定点的位置偏差值;对位置偏差值用最小二乘方法进行数据处理,获得一天内的轨道平面内的平均轨道根数,再用平均滤波方法获得一天内的轨道平面外的平均轨道根数;以平均轨道根数作为星上轨道解析外推算法的输入值,外推计算一天内卫星轨道位置,提供连续导航定位信息,实现卫星自主导航功能。本发明的方法已经在中星2A上成功应用,经过在轨标定后自主轨道确定精度优于10km,该方法可以推广应用于所有要求具备自主功能的地球静止轨道卫星。
Description
技术领域
本发明涉及一种静止轨道卫星自主导航方法,属于卫星自主导航研究领域,可以应用于静止轨道卫星自主导航任务。
背景技术
传统卫星在寿命周期内完全依赖地面测控系统的支持,与之相比,中星2A卫星是我国第一颗具有自主生存能力的地球静止轨道通信卫星,其所采用的自主导航方法为该技术首次工程应用,所用的方法属于天文导航方法,这一大类方法从算法角度已经较为成熟,并也产生了多种新型敏感器,如空间六分仪等,如中星2A上用到的基于地球敏感器和星敏感器的导航滤波算法是比较适合目前发展趋势的方法,首先是方法较为简单(已经获得我国专利保护),其次是方法所用到的部件属于星上常用姿态确定部件,能够长期在轨不间断工作,即导航方案不对星上资源产生新的需求;在原有导航滤波方法的基础上,中星2A针对静止轨道卫星的运行要求,尤其是轨道控制要求,创造性将星上轨道解析外推算法和导航滤波后的处理结果有机结合起来,既使得卫星自主获取了实时位置,同时也获得了适合于轨道控制的轨道信息。
发明内容
本发明的目的是为了提出一种静止轨道卫星自主导航方法,该方法能够给出静止轨道卫星实时的轨道信息。
本发明的目的是通过以下技术方案实现的。
本发明的一种静止轨道卫星自主导航方法,该方法的步骤为:
1)基于星敏感器和地球敏感器的kalman滤波算法,实时得到卫星相对于定点的位置偏差值;
2)对步骤1)得到的位置偏差值用最小二乘方法进行数据处理,获得一天内的轨道平面内的平均轨道根数:即半长轴、偏心率、近地点幅角和真近点角;对步骤1)得到的位置偏差值用平均滤波方法获得一天内的轨道平面外的平均轨道根数:即轨道倾角和升交点赤经;
3)以步骤2)得到的平均轨道根数作为星上轨道解析外推算法的输入值,外推计算一天内卫星轨道位置,提供连续导航定位信息,实现卫星自主导航功能。
有益效果
本发明的方法已经在中星2A上成功应用,经过在轨标定后自主轨道确定精度优于10km,该方法可以推广应用于所有要求具备自主功能的地球静止轨道卫星。
附图说明
图1为卫星自主导航测量原理示意图;
图2为轨道平面参数确定示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明。
实施例
1)卫星自主导航所需的星上测量部件为星敏感器和地球敏感器,基于星敏感器和地球敏感器的kalman滤波算法
以此标称轨道定点位置建立Hill方程,即卫星相对标称轨道定点位置的相对运动方程,作为kalman滤波算法的状态方程,方程形式如下:
其中:
α为经度偏差,单位为rad;
β为纬度偏差,单位为rad;
γ为径向偏差,单位为rad;
r为静止轨道卫星理论半径,单位为km;
ω0为标称静止轨道角速度,单位为rad/s;
ax、ay、az是外界加速度在标称轨道定点位置为原点的坐标系的分量;
标称轨道定点位置定义为:不考虑卫星漂移且轨道倾角为零的情况下的理论定点位置,与地球自转完全同步;
卫星自主导航测量原理如图1所示,包括地球、卫星和标称轨道定点位置;其中,表示星敏感器给出的第一个方向矢量,表示星敏感器给出的第二个方向矢量,代表地球敏感器给出的地心方向矢量。
方向矢量在卫星本体系下的坐标位置由星敏的安装矩阵决定;根据星敏感器的测量原理,其测量输出为方向矢量在惯性坐标系中的方位坐标,分别记为:和具体表达为
红外地球敏感器的测量输出为滚动角φh和俯仰角θh。根据测量得到的滚动角和俯仰角,经过两次旋转变换即可得到在本体坐标系下卫星轨道坐标系原点到地心的位置矢量,记为:
从标称轨道定点位置的指向地心的单位矢量记为
在卫星本体坐标系得到两个夹角η1、η2:
在惯性坐标系下计算得到另外两个夹角η1′、η2′:
可得到测量方程g:
g=(cosη1-cosη′1 cosη2-cosη′2)T
为此构建一个自主导航的模型,其对象方程为:
其中:
b=[0 0 0 1 1 1]T
U为外部输入。
将测量方程在参考轨道处线性化:
其中:
h11=-sin(αT+α)cos(δT-β)S1ix+cos(αT+α)cos(δT-β)S1iy
h12=cos(αT+α)sin(δT-β)S1ix+sin(αT+α)sin(δT-β)S1iy-cos(δT-β)S1iz
h21=-sin(αT+α)cos(δT-β)S2ix+cos(αT+α)cos(δT-β)S2iy
h22=cos(αT+α)sin(δT-β)S2ix+sin(αT+α)sin(δT-β)S2iy-cos(δT-β)S2iz
αT为标称轨道定点位置的赤经,δT为标称轨道定点位置的赤纬
根据上述状态方程和测量方程,采用扩展Kalman滤波方法,可获得当前卫星轨道相对于标称轨道定点位置的偏移量,即:
2)对步骤1)得到的位置偏差值用最小二乘方法进行数据处理,获得一天内的轨道平面内的平均轨道根数:即半长轴、偏心率、近地点幅角和真近点角;对步骤1)得到的位置偏差值用平均滤波方法获得一天内的轨道平面外的平均轨道根数:即轨道倾角和升交点赤经;
根据步骤1)中kalman滤波后获得的真实卫星相对于标称轨道定点位置的相对位置偏差值,为了获得便于轨道控制的控制变量以及提高导航精度,取滤波结果中的东西偏差α和南北偏差β,基于静止轨道卫星的受摄漂移原理,通过最小二乘法和平均滤波方法获得卫星的平均轨道六根数,包括:半长轴、偏心率、轨道倾角、升交点赤经、近地点幅角和真近点角,其中平均的意义为:利用一天的α和β数据拟合获得。
根据α和β,可以转换成卫星的赤经和赤纬,即赤经为αT+α:,赤纬为δT-β,可以获得当前k时刻的地心单位矢量其表达式为:
由若干组这样的单位矢量即可确定卫星的轨道平面参数,即轨道倾角i和升交点赤经Ω,采用的方法为平均滤波方法。
在一个轨道周期T内,首先每隔设定时间dt,可获得一个地心单位矢量,分别记为:(k=1...N),其中N=T/dt,如图2所示,则通过叉乘可得一个倾角矢量:
最终形成一个总的倾角矢量:
如此在一个轨道周期内通过一系列倾角矢量得到一个总的倾角矢量,根据倾角矢量在惯性坐标系中的XOY平面内投影定义:[sin i cos(Ω-90°),sin i sin(Ω-90°)],可解算出一组轨道倾角i和升交点赤经Ω。由于该倾角矢量是由一个轨道周期内的数据平均后获得,所以可以认为求出的i和Ω为上一个轨道周期内轨道倾角和升交点赤经的平均值,即一天数据可获得一组和
采用递推最小二乘法对平面内轨道参数进行估计,即卫星东西方向运行规律可简化为:
α=a0+a1torb+b1cos(ωetorb)+b2sin(ωetorb)
其中,torb一个轨道周期的时间变量,ωe为静止轨道角速度,待估参数为a0、a1、b1、b2。
通过递推最小二乘法可实时估计参数a0、a1、b1、b2,由此在一个轨道周期结束时,可解算出卫星轨道平面内的平均轨道根数:
半长轴可由a1直接解出:其中rs为静止轨道标称半径,a1对应的单位为rad/天。
根据数学模型2esin(ωetorb+θ),其中θ=atan2(b1,b2)为相位偏移量,由b1、b2可解算出偏心率即:
并且,令 与当前的卫星地心单位矢量进行内积,再求反余弦获得卫星幅角u,另外,真近点角则:
3)以步骤2)得到的六个平均轨道根数作为星上轨道解析外推算法的输入值,外推计算一天内卫星轨道位置,提供连续导航定位信息,实现卫星自主导航功能。
由于地球同步轨道卫星偏心率e和轨道倾角i量级都比较小,故在进行解析外推计算时需采用适用于0≤e<1和0≤i<180°的第二类六个无奇点根数:
σ=[a,h,k,ξ,η,λ]T
其表达是为:a为半长轴,h=sin i cosΩ,k=sin i sinΩ, 其中
将步骤2)得到的六个平均轨道根数代入上式,令f≈M,可获得第二类无奇点根数的平均值: 其中的即为作为星上轨道解析外推算法的初始拟平均根数,且根据步骤2)的结果每天更新一次,进而获得自主导航期间的实时轨道根数。
采用拟平均根数法构造的分析解可以获得瞬时轨道根数σ(t):
其中:
σ(t)为瞬时轨道根数
为拟平均轨道根数
σ1为一阶长期项系数,包含长期项和长周期项两部分
为一阶短周期项
为初始拟平均根数
t0为初始拟平均根数所对应的轨道历元
为轨道平均角速度
地球非球形长期项系数σ1:
a1=0
λ1=M1+ω1+Ω1
其中Ω1,ω1和M1由下式表达:
上述各式中的a,e,i均为且有:
地球非球形短周期项
上述短周期项各式右端出现的根数均为而u和a/r由下列相应公式进行迭代计算获得:
其中u=f+ω,E是偏近点角。
地球非球形长周期项系数:
上述各式中的a,h,k等量均为各式中的谐系数J2,2以及等量的意义如下:
其中,ne为标称轨道角速度,C2,2,S2,2是非归一化的谐系数,d为J2000.0起算的儒略日,SG是格林尼治恒星时,λ2,2=-14.545°。
日月引力摄动部分:
长期项系数σ1:
a1=0
各式右端出现的a,e,i和n等量均为轨道根数和日、月的有关量β定义如下:
(其中M为地球质量,S为日或月质量,r′为日或月至地心的距离)
日、月位置量可采用t0时刻的值。式中出现的几个量A,B,Ai,Bi,AΩ,,BΩ定义如下,其中的i′,Ω′,ω′,u′分别为日或月在惯性坐标系下的轨道根数:
短周期项
式中变量含义同上。
Claims (4)
1.一种静止轨道卫星自主导航方法,其特征在于该方法的步骤为:
1)基于星敏感器和地球敏感器的kalman滤波算法,实时得到卫星相对于定点的位置偏差值;
2)对步骤1)得到的位置偏差值用最小二乘方法进行数据处理,获得一天内的轨道平面内的平均轨道根数:即半长轴、偏心率、近地点幅角和真近点角;对步骤1)得到的位置偏差值用平均滤波方法获得一天内的轨道平面外的平均轨道根数:即轨道倾角和升交点赤经;
3)以步骤2)得到的平均轨道根数作为星上轨道解析外推算法的输入值,外推计算一天内卫星轨道位置,提供连续导航定位信息,实现卫星自主导航功能。
2.根据权利要求1所述的一种静止轨道卫星自主导航方法,其特征在于:步骤1)的具体步骤为:
以此标称轨道定点位置建立Hill方程,即卫星相对标称轨道定点位置的相对运动方程,作为kalman滤波算法的状态方程,方程形式如下:
其中:
α为经度偏差,单位为rad;
β为纬度偏差,单位为rad;
γ为径向偏差,单位为rad;
r为静止轨道卫星理论半径,单位为km;
ω0为标称静止轨道角速度,单位为rad/s;
ax、ay、az是外界加速度在标称轨道定点位置为原点的坐标系的分量;
标称轨道定点位置定义为:不考虑卫星漂移且轨道倾角为零的情况下的理论定点位置,与地球自转完全同步;
卫星自主导航测量要用到地球、卫星和标称轨道定点位置;其中,表示星敏感器给出的第一个方向矢量,表示星敏感器给出的第二个方向矢量,代表地球敏感器给出的地心方向矢量;
方向矢量在卫星本体系下的坐标位置由星敏的安装矩阵决定;根据星敏感器的测量原理,其测量输出为方向矢量在惯性坐标系中的方位坐标,分别记为:和具体表达为
红外地球敏感器的测量输出为滚动角φh和俯仰角θh;根据测量得到的滚动角和俯仰角,经过两次旋转变换即可得到在本体坐标系下卫星轨道坐标系原点到地心的位置矢量,记为:
从标称轨道定点位置的指向地心的单位矢量记为
在卫星本体坐标系得到两个夹角η1、η2:
在惯性坐标系下计算得到另外两个夹角η1′、η2′:
可得到测量方程g:
为此构建一个自主导航的模型,其对象方程为:
其中:
b=[0 0 0 1 1 1]T
U为外部输入;
将测量方程在参考轨道处线性化:
其中:
h11=-sin(αT+α)cos(δT-β)S1ix+cos(αT+α)cos(δT-β)S1iy
h12=cos(αT+α)sin(δT-β)S1ix+sin(αT+α)sin(δT-β)S1iy-cos(δT-β)S1iz
h21=-sin(αT+α)cos(δT-β)S2ix+cos(αT+α)cos(δT-β)S2iy
h22=cos(αT+α)sin(δT-β)S2ix+sin(αT+α)sin(δT-β)S2iy-cos(δT-β)S2iz
αT为标称轨道定点位置的赤经,δT为标称轨道定点位置的赤纬
根据上述状态方程和测量方程,采用扩展Kalman滤波方法,可获得当前卫星轨道相对于标称轨道定点位置的偏移量,即:
3.根据权利要求1所述的一种静止轨道卫星自主导航方法,其特征在于:步骤2)的具体步骤为:对步骤1)得到的位置偏差值用最小二乘方法进行数据处理,获得一天内的轨道平面内的平均轨道根数:即半长轴、偏心率、近地点幅角和真近点角;对步骤1)得到的位置偏差值用平均滤波方法获得一天内的轨道平面外的平均轨道根数:即轨道倾角和升交点赤经;
根据步骤1)中kalman滤波后获得的真实卫星相对于标称轨道定点位置的相对位置偏差值,为了获得便于轨道控制的控制变量以及提高导航精度,取滤 波结果中的东西偏差α和南北偏差β,基于静止轨道卫星的受摄漂移原理,通过最小二乘法和平均滤波方法获得卫星的平均轨道六根数,包括:半长轴、偏心率、轨道倾角、升交点赤经、近地点幅角和真近点角,其中平均的意义为:利用一天的α和β数据拟合获得;
根据α和β,可以转换成卫星的赤经和赤纬,即赤经为αT+α,赤纬为δT-β,可以获得当前k时刻的地心单位矢量其表达式为:
由若干组这样的单位矢量即可确定卫星的轨道平面参数,即轨道倾角i和升交点赤经Ω,采用的方法为平均滤波方法;
在一个轨道周期T内,首先每隔设定时间dt,可获得一个地心单位矢量,分别记为:(k=1...N),其中N=T/dt,则通过叉乘可得一个倾角矢量:
最终形成一个总的倾角矢量:
如此在一个轨道周期内通过一系列倾角矢量得到一个总的倾角矢量,根据倾角矢量在惯性坐标系中的XOY平面内投影定义:[sinicos(Ω-90°),sinisin(Ω-90°)],可解算出一组轨道倾角i和升交点赤经Ω;由于该倾角矢量是由一个轨道周期内的数据平均后获得,所以可以认为求出的i和Ω为上一个轨道周期内轨道倾角和升交点赤经的平均值,即一天数据可获得一组和
采用递推最小二乘法对平面内轨道参数进行估计,即卫星东西方向运行规律可简化为:
α=a0+α1torb+b1cos(ωetorb)+b2sin(ωetorb)
其中,torb一个轨道周期的时间变量,ωe为静止轨道角速度,待估参数为a0、a1、b1、b2;
通过递推最小二乘法可实时估计参数a0、a1、b1、b2,由此在一个轨道周期结束时,可解算出卫星轨道平面内的平均轨道根数:
半长轴可由a1直接解出:其中rs为静止轨道标称半径,a1对应的单位为rad/天;
根据数学模型2esin(ωetorb+θ),其中θ=a tan2(b1,b2)为相位偏移量,由b1、b2可解算出偏心率即:
并且,令与当前的卫星地心单位矢量进行内积,再求反余弦获得卫星幅角u,另外,真近点角则:
4.根据权利要求1所述的一种静止轨道卫星自主导航方法,其特征在于:步骤3)的具体步骤为:以步骤2)得到的六个平均轨道根数作为星上轨道解析外推算法的输入值,外推计算一天内卫星轨道位置,提供连续导航定位信息,实现卫星自主导航功能;
由于地球同步轨道卫星偏心率e和轨道倾角i量级都比较小,故在进行解析外推计算时需采用适用于0≤e<1和0≤i<180°的第二类六个无奇点根数:
σ=[a,h,k,ξ,η,λ]T
其表达式为:a为半长轴,h=sinicosΩ,k=sinisinΩ, 其中
将步骤2)得到的六个平均轨道根数代入上式,令f≈M,可获得第二类无奇点根数的平均值:其中的即为作为星上轨道解析外推算法的初始拟平均根数,且根据步骤2)的结果每天更新一次,进而获得自主导航期间的实时轨道根数;
采用拟平均根数法构造的分析解可以获得瞬时轨道根数σ(t):
其中:
σ(t)为瞬时轨道根数
为拟平均轨道根数
σ1为一阶长期项系数,包含长期项和长周期项两部分
为一阶短周期项
为初始拟平均根数
t0为初始拟平均根数所对应的轨道历元
为轨道平均角速度
地球非球形长期项系数σ1:
a1=0
λ1=M1+ω1+Ω1
其中Ω1,ω1和M1由下式表达:
上述各式中的a,e,i均为且有:
地球非球形短周期项
上述短周期项各式右端出现的根数均为而u和a/r由下列相应公式进行迭代计算获得:
其中u=f+ω,E是偏近点角;
地球非球形长周期项系数:
上述各式中的a,h,k等量均为各式中的谐系数J2,2以及等量的意义如下:
其中,ne为标称轨道角速度,d为J2000.0起算的儒略日,SG是格林尼治恒星时,λ2,2=-14.545°;
日月引力摄动部分:
长期项系数σ1:
a1=0
各式右端出现的a,e,i和n等量均为轨道根数和日、月的有关量β定义如下:
(其中M为地球质量,S为日或月质量,r′为日或月至地心的距离)
日、月位置量可采用t0时刻的值;式中出现的几个量A,B,Ai,Bi,AΩ,BΩ定义如下,其中的i′,Ω′,ω′,u′分别为日或月在惯性坐标系下的轨道根数:
短周期项
式中变量含义同上。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210413994.0A CN102878995B (zh) | 2012-10-24 | 2012-10-24 | 一种静止轨道卫星自主导航方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210413994.0A CN102878995B (zh) | 2012-10-24 | 2012-10-24 | 一种静止轨道卫星自主导航方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102878995A CN102878995A (zh) | 2013-01-16 |
CN102878995B true CN102878995B (zh) | 2014-12-17 |
Family
ID=47480415
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210413994.0A Active CN102878995B (zh) | 2012-10-24 | 2012-10-24 | 一种静止轨道卫星自主导航方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102878995B (zh) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103279127B (zh) * | 2013-05-22 | 2016-06-22 | 上海新跃仪表厂 | 一种仅用角度信息的geo轨道卫星自主控制方法 |
CN105737834B (zh) * | 2014-12-09 | 2018-06-26 | 上海新跃仪表厂 | 一种基于轨道平根数的相对导航鲁棒滤波方法 |
CN104574422B (zh) * | 2015-01-30 | 2017-04-19 | 北京控制工程研究所 | 一种低信噪比红外地球图像信息处理方法 |
CN106197425B (zh) * | 2016-06-30 | 2019-04-19 | 中国电子科技集团公司第五十四研究所 | 基于卫星姿态角的地面目标点位置的计算方法 |
CN106628257B (zh) * | 2016-09-28 | 2018-01-19 | 西北工业大学 | 地球摄动引力场中近地航天器相对运动轨道的保持方法 |
CN106697333B (zh) * | 2017-01-12 | 2019-09-06 | 北京理工大学 | 一种航天器轨道控制策略的鲁棒性分析方法 |
CN107031868B (zh) * | 2017-03-23 | 2019-06-18 | 北京空间飞行器总体设计部 | 一种低轨遥感卫星自主轨道控制方法 |
CN107656293A (zh) * | 2017-08-30 | 2018-02-02 | 清华大学 | 基于无奇点根数的14参数广播星历卫星位置估计方法 |
CN107797130B (zh) * | 2017-10-16 | 2021-01-05 | 中国西安卫星测控中心 | 低轨航天器多点多参数轨道上行数据计算方法 |
CN108279426B (zh) * | 2018-01-24 | 2020-05-26 | 北京电子工程总体研究所 | 一种测控站至卫星星下点航路捷径的解析计算方法 |
CN109000666B (zh) * | 2018-06-05 | 2021-11-23 | 北京电子工程总体研究所 | 一种基于中心天体矢量观测的自主定轨方法及其系统 |
CN109579829B (zh) * | 2018-11-29 | 2022-06-14 | 天津津航技术物理研究所 | 一种小视场星敏感器短波导航星识别方法 |
CN109506645B (zh) * | 2018-12-13 | 2020-09-18 | 上海航天控制技术研究所 | 一种星敏感器安装矩阵地面精测方法 |
CN110068846B (zh) * | 2019-04-30 | 2022-01-07 | 上海微小卫星工程中心 | 一种基于星载gnss接收机在星上自主确定轨道平根数的方法 |
CN113740887A (zh) * | 2019-04-30 | 2021-12-03 | 上海微小卫星工程中心 | 一种卫星注入轨道外推及卫星理论轨道确定方法 |
CN111428365B (zh) * | 2020-03-24 | 2021-12-28 | 中国人民解放军32035部队 | 利用天文测量数据判别geo目标的方法 |
CN113984069B (zh) * | 2021-04-30 | 2023-06-06 | 北京临近空间飞行器系统工程研究所 | 基于人造卫星的星光定位导航方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101226062A (zh) * | 2007-12-26 | 2008-07-23 | 北京控制工程研究所 | 一种星上实时计算环月轨道的方法 |
CN101893712A (zh) * | 2010-07-09 | 2010-11-24 | 中国科学院测量与地球物理研究所 | 用于地球静止卫星精密定轨的选权拟合方法 |
CN102305630A (zh) * | 2011-05-17 | 2012-01-04 | 哈尔滨工业大学 | 基于扩展卡尔曼滤波的sar卫星自主定轨方法 |
CN102322862A (zh) * | 2011-06-29 | 2012-01-18 | 航天东方红卫星有限公司 | 一种编队飞行卫星绝对和相对轨道确定方法 |
CN102519472A (zh) * | 2011-12-08 | 2012-06-27 | 北京控制工程研究所 | 一种利用偏航机动的自主导航敏感器系统误差校正方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0434310A (ja) * | 1990-05-31 | 1992-02-05 | Mitsubishi Electric Corp | 宇宙機の航法装置 |
US7668629B2 (en) * | 2006-12-05 | 2010-02-23 | The Boeing Company | Ultra-tightly coupled global navigation satellite system space borne receiver system |
EP2603769A4 (en) * | 2010-08-12 | 2015-05-13 | Us Gov Sec Navy | IMPROVED SYSTEM AND METHOD OF CARRYING UP COVARION MEASUREMENT AND ANALYSIS |
-
2012
- 2012-10-24 CN CN201210413994.0A patent/CN102878995B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101226062A (zh) * | 2007-12-26 | 2008-07-23 | 北京控制工程研究所 | 一种星上实时计算环月轨道的方法 |
CN101893712A (zh) * | 2010-07-09 | 2010-11-24 | 中国科学院测量与地球物理研究所 | 用于地球静止卫星精密定轨的选权拟合方法 |
CN102305630A (zh) * | 2011-05-17 | 2012-01-04 | 哈尔滨工业大学 | 基于扩展卡尔曼滤波的sar卫星自主定轨方法 |
CN102322862A (zh) * | 2011-06-29 | 2012-01-18 | 航天东方红卫星有限公司 | 一种编队飞行卫星绝对和相对轨道确定方法 |
CN102519472A (zh) * | 2011-12-08 | 2012-06-27 | 北京控制工程研究所 | 一种利用偏航机动的自主导航敏感器系统误差校正方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102878995A (zh) | 2013-01-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102878995B (zh) | 一种静止轨道卫星自主导航方法 | |
CN101788296B (zh) | 一种sins/cns深组合导航系统及其实现方法 | |
CN100501331C (zh) | 基于x射线脉冲星的导航卫星自主导航系统与方法 | |
Konopliv et al. | A global solution for the Mars static and seasonal gravity, Mars orientation, Phobos and Deimos masses, and Mars ephemeris | |
CN103076015B (zh) | 一种基于全面最优校正的sins/cns组合导航系统及其导航方法 | |
CA2973741C (en) | Inertial sensing augmentation for navigation of spacecraft | |
CN104298647B (zh) | 基于低轨道地球卫星的地影时刻预报的星上确定方法 | |
CN102880184A (zh) | 一种静止轨道卫星自主轨道控制方法 | |
CN103017760B (zh) | 一种大椭圆轨道火星探测器自主对火定向方法 | |
CN106595674A (zh) | 基于星敏感器和星间链路的heo卫星编队飞行自主导航方法 | |
CN108562295B (zh) | 一种基于同步卫星二体模型的三站时差定轨方法 | |
Zurek et al. | Application of MAVEN accelerometer and attitude control data to Mars atmospheric characterization | |
CN103900611A (zh) | 一种惯导天文高精度复合两位置对准及误差标定方法 | |
CN103868514A (zh) | 一种在轨飞行器自主导航系统 | |
CN106643741A (zh) | 一种卫星相对小行星视觉自主导航方法 | |
CN104457705A (zh) | 基于天基自主光学观测的深空目标天体初定轨方法 | |
CN104006813A (zh) | 一种高轨卫星的脉冲星/星光角距组合导航方法 | |
CN104048664A (zh) | 一种导航卫星星座自主定轨的方法 | |
CN106595673A (zh) | 面对地球静止轨道目标操作的空间多机器人自主导航方法 | |
CN103968834A (zh) | 一种近地停泊轨道上深空探测器的自主天文导航方法 | |
CN112014869A (zh) | 基于天文导航的星间链路自主导航方法及系统 | |
CN113589832A (zh) | 对地表固定区域目标稳定观测覆盖的星座快速设计方法 | |
CN114684389B (zh) | 考虑再入约束的月地转移窗口及精确转移轨道确定方法 | |
CN115392540A (zh) | 一种用于月球轨道交会制导的快速预报方法 | |
CN112013834B (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |