CN103063217B - 一种基于星历修正的深空探测器天文/无线电组合导航方法 - Google Patents
一种基于星历修正的深空探测器天文/无线电组合导航方法 Download PDFInfo
- Publication number
- CN103063217B CN103063217B CN201310006575.XA CN201310006575A CN103063217B CN 103063217 B CN103063217 B CN 103063217B CN 201310006575 A CN201310006575 A CN 201310006575A CN 103063217 B CN103063217 B CN 103063217B
- Authority
- CN
- China
- Prior art keywords
- err
- navigation
- celestial body
- subsystem
- error
- 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
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明涉及一种基于目标天体星历修正的深空探测器捕获段天文/无线电组合导航方法,首先建立火星探测器状态模型、天文导航和无线电导航子***量测模型,然后获取天文和无线电导航子***的量测量,滤波估计得到探测器在日心和以目标天体为中心惯性坐标系中的位置和速度;在此基础上,建立目标天体星历误差的状态模型和量测模型,并由天文和无线电两个导航子***的估计状态向量获得目标天体星历误差量测量,利用卡尔曼滤波方法估计目标天体星历误差,并反馈至导航***模型中,再进行信息融合。本发明属于航天导航技术领域,可在线估计天体星历误差,修正导航***的模型误差,适用于探测器捕获段。
Description
技术领域
本发明涉及在深空探测器处于捕获段时,基于目标天体图像和无线电信号时间延迟和频移信息的组合导航方法,是一种非常适用于深空探测器捕获段的组合导航方法。
背景技术
深空探测技术作为一个国家综合国力和科学技术发展水平的重要特征与标志,已引起世界各国的极大关注。新一轮深空探测的争夺战已经拉开了序幕。21世纪初,各航天大国纷纷将目光聚焦至距离地球38万公里以外的深空宇宙。美国、欧空局、俄罗斯、日本以及印度在内的世界主要航天集团都提出了未来的深空探测计划,要对各大行星及其卫星进行载人或基于机器人的无人探测。随着我国运载火箭和其他深空探测技术的发展和经济实力的提高,我国已具备探测火星甚至更远太阳系行星的能力。
深空探测器飞行过程主要包括地球逃逸、日心转移、目标捕获、环绕、着陆、巡视探测等过程。其中捕获段是指从深空探测器进入目标影响球开始到点火制动的全过程,处于该阶段的深空探测器飞行速度快,飞行弧段短,控制精度要求高,且机会唯一。深空探测器减速制动的入轨点距离目标行星表面非常近(近地点),被捕获是整个深空探测任务的一个关键时间节点,捕获阶段相对目标天体的相对导航精度和相对于太阳或者地球的绝对导航精度对后续探测任务有直接影响。然而深空探测器在捕获段航行速度快、空间电离环境未知、大气环境复杂,这些因素都对深空探测器的入轨精度有很大的影响,也制约着深空探测器捕获后绕飞、着陆等阶段的导航精度,当入轨精度不能达到指标要求时,甚至无法完成科学探测任务,导致整个任务的失败。
目标天体星历数据是影响深空探测器捕获段导航性能的主要因素之一。目标天体星历数据是描述目标天***置、速度等特征的一类天体数据库,由长期天文观测拟合而成的。若一段时间没有天文观测信息进行修正,目标天体星历数据误差会随时间的递推而增加。太阳系内行星(除地球外)目前的星历误差约为200m~100km。
现有的深空探测导航方式主要包括地面无线电导航和天文导航两种导航方式,其中无线电导航的主要观测量为探测器相对于地面站的距离与距离变化率等,这种导航方式可以测量相对于太阳的导航信息;而天文导航是测量目标天体及其背景恒星图像信息,从而从图像中获得相对于目标天体(如火星)的导航信息。由于目标天体星历误差大、天文敏感器测量精度有限等因素的影响,这两种方法仍存在以下不足:无线电导航方法仅能获取相对太阳的导航信息,因此由该方法获得的相对目标天体导航信息受到目标天体星历误差的影响大,相对目标天体导航精度低;天文导航方法由于受到天文导航敏感器精度有限这一因素的影响,其获得的相对目标天体导航信息精度有限,且受到目标天体星历误差影响,由该方法获得的相对太阳导航精度低。因此如何有效利用无线电导航和天文导航所提供的两种测量信息、减小星历误差对导航性能的影响,是捕获段高精度导航的一个关键问题。
传统的天文/无线电组合导航方法有两种:一种方法为将无线电导航的结果转换至目标天体为中心的惯性坐标系中后,将天文导航子***和无线电子***的导航结果进行信息融合。这种方法计算简单,易于实现,但由于在信息融合过程中没有考虑星历误差,因此无法克服目标天体星历误差对状态模型精度的影响,限制了导航精度的提高,导航精度较低。另一种方法为将无线电导航信息与天文导航信息组合进行几何解析,求解目标天体星历及星历误差,此后将无线电导航导航结果转换至目标天体为中心的惯性坐标系中,获得较高精度的导航信息。但是这种方法由于没有实时修正状态模型的模型误差,仍然无法消除目标天体星历误差对探测器状态模型精度的影响,导航精度仍较低。
发明内容
本发明要解决的技术问题是:克服星历误差对无线电导航信息相对目标天体精度的影响,弥补现有方法难以消除目标天体星历误差对探测器状态模型的影响这一不足,合理有效利用无线电导航和天文导航所提供的导航信息,为深空探测器捕获段提供一种高精度的组合导航方法。
本发明解决其技术问题所采用的技术方案为:建立目标天体为中心惯性坐标系中基于太阳和八大行星引力的深空探测器状态模型,通过光学导航敏感器获得目标天体及其卫星与恒星之间角度信息量测量,建立天文导航***的量测模型,并使用Unscented卡尔曼滤波方法获得深空探测器相对目标天体的位置和速度参数;建立日心惯性坐标系中基于太阳和八大行星引力的深空探测器状态模型,通过无线电接收机获得探测器相对太阳的位置、速度量测量,建立无线电导航***的量测模型,并使用Unscented卡尔曼滤波方法获得深空探测器相对日心的位置和速度参数;根据目标天体星历误差特性,建立目标天体星历误差状态模型和量测模型,利用无线电子***和天文导航子***的导航信息,获得目标天体星历误差量测量,并使用卡尔曼滤波方法获得目标天体星历误差的估计状态向量与估计均方误差阵,并将其反馈回无线电导航***及天文导航***的状态模型中,修正状态模型误差;利用星历修正***获得的星历误差,将无线电导航***的导航信息转换至以目标天体为中心的惯性坐标系中,与天文导航***的导航信息进行融合,最终为处于捕获段的深空探测器提供位置、速度等导航参数。
具体包括以下步骤:
1.建立基于太阳和八大行星引力轨道动力学的深空探测器导航***状态模型
A.在目标天体为中心的惯性坐标系中建立深空探测器基于太阳和八大行星引力轨道动力学的第一状态模型,即天文导航子***的状态模型;
式中,X′(t)=[x′,y′,z′,vx′,v′y,vz′]T为状态向量,f1(X′(t),t)为***非线性连续状态转移函数,w′(t)=[wx′,w′y,wz′]T为状态模型误差。
B.在日心惯性坐标系中建立深空探测器基于太阳和八大行星引力轨道动力学的第二状态模型,即无线电导航子***的状态模型;
式中,X(t)=[x,y,z,vx,vy,vz]T为状态向量,f2(X(t),t)为***非线性连续状态转移函数,w(t)=[wx,wy,wz]T为状态模型误差。
2.分别建立天文导航子***和无线电导航子***量测模型
(1)天文导航子***量测模型
目标天体和两个卫星与三颗背景恒星的角度信息θ1i、θ2i和θ3i(i=1,2,3)表达式为:
式中,为在惯性坐标系中目标天体到探测器的单位矢量,为在惯性坐标系中目标天体图像中第i个恒星的单位矢量;为在惯性坐标系中目标天体的第一个卫星到探测器的单位矢量,为在惯性坐标系中目标天体第一个卫星图像中i个恒星的单位矢量;为在惯性坐标系中第二个目标天体的卫星到探测器的单位矢量,为在惯性坐标系中目标天体第二个卫星图像中i个恒星的单位矢量。
设天文导航子***量测量Z1=[θ11,θ12,θ13,θ21,θ22,θ23,θ31,θ32,θ33]T,天文导航子***量测噪声 分别为测量θ11,θ12,θ13,θ21,θ22,θ23,θ31,θ32,θ33的观测误差,由于各变量都是与时间t有关的变量,则可根据式(3)建立天文导航子***量测模型的表达式为:
Z1(t)=h1[X(t),t]+v1(t) (4)
式中,h1[X′(t),t]为天文导航子***非线性连续量测函数。
(2)无线电导航子***量测模型
无线电测控站可根据无线电信号的时延与频移测量探测器到测控站的距离与距离变化率,一般主测控站附近还配有多个副站,根据多个测控站测量的距离与距离变化率可以获得探测器相对于地心的位置与速度,无线电导航所获得的距离与距离变化率表达式为:
式中,Ri为第i个测控站的测距量测信息(探测器与测控站之间的距离),为第i个测控站的距离变化率量测信息(探测器与测控站之间的距离变化率),(xf,yf,zf)为地面测控站在地心惯性坐标系中的位置坐标。
设无线电导航子***量测量无线电导航子***量测噪声 分别为三个测控站所测量得到的探测器与测控站之间的距离和距离变化率,分别为测量的观测误差,由于各变量都是与时间t有关的变量,则可根据式(5)和式(6)建立无线电导航子***量测模型的表达式为:
Z2(t)=h2[X(t),t]+v2(t) (7)
式中,h2[X(t),t]为无线电导航子***非线性连续量测函数。
3.对步骤1和步骤2中的状态模型和量测模型进行离散化
X′(k)=F1(X′(k-1),k-1)+W′(k-1) (8)
X(k)=F2(X(k-1),k-1)+W(k-1) (9)
Z′(k)=H1(X′(k),k)+V1(k) (10)
Z(k)=H2(X(k),k)+V2(k) (11)
式中,k=1,2,…,F1(X′(k-1),k-1)和F2(X(k-1),k-1)分别为f1(X′(t),t)和f2(X(t),t)离散后从第k-1时刻到第k时刻的非线性状态转移函数,H1(X′(k),k)和H2(X(k),k)分别为h1(X′(t),t)和h2(X(t),t)离散后第k时刻的非线性量测函数,W′(k),W(k),V1(k),V2(k)分别为w′(t),w(t),v1(t)和v2(t)离散后第k时刻的等效噪声,且W′(k)和V1(k)、W(k)和V2(k)互不相关。
4.天文导航和无线电导航量测量的获取及处理
(1)天文导航子***量测量的获取及处理
由天文导航敏感器获得目标天体的图像,利用图像处理技术,确定天体质心的像元像线坐标;经过从像元像线坐标系到二维成像平面坐标系、从二维成像平面坐标系到敏感器测量坐标系的三次转换,确定天体及其背景恒星在敏感器坐标系中的单位矢量;最后计算天体单位矢量与背景恒星单位矢量间的星光角距。
(2)无线电导航子***量测量的获取及处理
探测器与测控站之间的距离变化率是通过测量探测器接收的载波信号频率f′rec与地面站发射的载波信号频率f0相比的频率偏移来获得的;探测器与测控站之间的距离是通过测量无线电信号由地面测控站发射至探测器再返回的时间延迟而获得的
5.对天文导航子***进行Unscented卡尔曼滤波
根据第一状态模型、天文导航子***量测模型、天文导航敏感器获得的量测量,进行天文导航子***Unscented卡尔曼滤波,获得在目标天体惯性坐标系中表示深空探测器位置和速度的估计状态向量和估计均方误差阵Pk′。
6.对无线电导航子***进行Unscented卡尔曼滤波
根据第二状态模型、无线电导航子***量测模型、由无线电接收装置获得的量测量,进行无线电导航子***Unscented卡尔曼滤波,获得在日心惯性坐标系中表示深空探测器位置和速度的估计状态向量和估计均方误差阵Pk。
7.判定是否需要进行目标天体星历校正
当天文导航子***Unscented卡尔曼滤波的估计均方误差阵Pk′大于已有目标天体星历误差均方误差阵时,即Pk′>Peph,则不进行星历校正,直接进行第9步信息融合;当Pk′<Peph时,则进行星历误差校正,执行第8步;
8.对目标天体星历误差进行建模、估计并反馈校正
A.建立目标天体星历误差状态模型
将捕获段内其误差特性考虑为常值误差,在日心惯性坐标系中建立目标天体星历误差状态模型为:
式中, 为日心惯性坐标系中目标天体星历的三轴位置误差的微分,可离散化后简写为:
Xerr(k)=Ferr(Xerr(k-1),k-1)+Werr(k-1) (13)
式中,状态转移函数Ferr(Xerr(k-1),k-1)=Φerr,k,k-1Xerr,k-1,Φerr,k,k-1为第k-1时刻到第k时刻的状态转移矩阵,Xerr(k)为第k时刻目标天体星历误差状态向量,且Xerr(k)=Xerr,k,Werr(k-1)为第k-1时刻目标天体星历误差状态模型误差。
B.建立目标天体星历误差量测模型
目标天体星历误差的量测模型可以表示为:
Zerr=H3(Xerr(k),k)+V3 (14)
式中,H3(Xerr(k),k)为k时刻的量测函数,V3为目标天体星历误差量测噪声。
C.获取目标天体星历误差量测量
目标天体星历误差量测量Zerr可以表示为:
式中,为无线电导航子***获得的相对于太阳的位置和速度,为天文导航子***获得的相对于目标天体的位置和速度,为目标天体相对于太阳的位置和速度,从天体星历数据库中获得。
D.对目标天体星历误差进行卡尔曼滤波估计
根据步骤A和步骤B建立的目标天体星历误差状态模型和量测模型,以及步骤C所获取的目标天体星历误差量测量,利用卡尔曼滤波方法,对目标天体星历误差进行估计,获得目标天体星历误差估计状态向量与估计均方误差阵,具体如下:
目标天体星历误差状态向量的一步预测
式中,为k-1时刻火星星历误差一步预测状态向量。
目标天体星历误差状态向量的估计均方误差阵一步预测
Perr,k/k-1=Φerr,k,k-1Perr,k-1Φerr,k,k-1 T+Qerr,k (17)
式中,Perr,k-1为k-1时刻目标天体星历误差状态向量的估计均方误差阵,Qerr,k为k时刻目标天体星历误差状态模型误差协方差阵。
卡尔曼滤波增益
Kerr,k=Perr,k/k-1Herr,k T(Herr,kPerr,k/k-1Herr,k T+Rerr,k)-1 (18)
式中,Herr,k为k时刻目标天体星历误差量测矩阵,Herr,kXerr,k=H3(Xerr,k),Rerr,k为k时刻目标天体星历误差量测模型误差协方差阵。
目标天体星历误差估计状态向量
式中,Zerr,k为k时刻目标天体星历误差量测量。
目标天体星历误差估计均方误差阵
Perr,k=(I-Kerr,kHerr,k)Perr,k/k-1 (20)
式中,I为单位阵。
E.对目标天体星历误差进行反馈校正
将步骤D中获得的目标天体星历误差和目标天体星历误差估计均方误差阵Perr,k反馈回深空探测器的第一状态模型和第二状态模型中,并重新确定第一状态模型和第二状态模型的模型误差协方差阵,最后将校正后的模型误差协方差阵输入至天文导航子***Unscented卡尔曼滤波和无线电导航子***Unscented卡尔曼滤波中,修正下一时刻的导航结果。
9.信息融合
利用星历修正***获得的星历误差,将无线电导航***的导航信息转换至以目标天体为中心的惯性坐标系中,与天文导航***的导航信息进行融合。当探测器不在无线电测控范围内时,无线电导航子***没有输入的无线电导航量测量,此时对天文导航子***进行Unscented卡尔曼滤波,无线电导航子***只利用第二状态模型进行时间更新;当探测器处于无线电测控范围内时,无线电导航子***有输入的无线电导航量测量,对两个导航子***同时进行Unscented卡尔曼滤波;
最终输出第k时刻在目标天体为中心的惯性坐标系中表示探测器位置和速度的估计状态向量和估计均方误差阵,并根据修正后的目标行星星历,将该结果转换至日心惯性坐标系中,输出在日心惯性坐标系中表示探测器位置和速度的估计状态向量和估计均方误差阵,将这些导航信息分别返回天文导航子***和无线电导航子***中,用于k+1时刻的位置、速度导航信息的估计,k=1,2,...;
本发明的原理是:首先选择基于太阳和八大行星引力的轨道动力学模型作为***状态模型,根据天文导航和无线电导航测量的不同特点,分别建立在目标天体为中心的惯性坐标系和日心惯性坐标系中的两个状态模型;之后根据天文导航敏感器和无线电导航接收装置的工作原理,获取并处理的天文导航敏感器和无线电导航接收装置所测量的量测量;然后,建立天文导航子***和无线电导航子***的量测模型;此后,由于探测器的状态模型和量测模型都呈现非线性特性,且***噪声为非高斯噪声,因此采用Unscented卡尔曼滤波方法,对两个子***的探测器导航信息进行估计;其次,由于天文导航可以测量高精度的相对目标天体导航信息,而无线电导航可以测量高精度的相对于太阳导航信息;结合目标天体、探测器、太阳的几何关系,可以确定目标天体的星历,与原有的目标天体星历相比较,即可得到目标天体星历误差量测量;为了获得更为准确的目标天体星历误差,结合捕获段持续时间短、目标天体星历误差在捕获段变化小的特点,建立捕获段目标天体星历误差的状态模型和量测模型,并根据目标天体星历误差状态模型和量测模型都为线性模型的特点,采用卡尔曼滤波方法,估计目标天体星历误差,并将所估计的目标天体星历误差反馈至第一状态模型和第二状态模型中,修正目标天体星历数据,提高状态模型的模型精度;最后通过信息融合方法,有效利用天文导航子***和无线电子***的测量信息,提高对相对于目标天体和相对于日心的探测器导航信息的估计精度。
本发明与现有技术相比的优点在于:利用星历校正过程所估计的目标天体星历误差,一方面实现了对目标天体星历误差的高精度估计,同时获得了相对目标天体和相对日心的探测器高精度导航信息,另一方面修正了导航***的状态模型,减小了目标天体星历误差对状态模型精度的影响,进一步提高了深空探测器的导航精度。
附图说明
图1为本发明基于星历修正的天文无线电组合导航流程图。
图2为本发明天文导航子***所用天文导航敏感器的成像原理图。
具体实施方式
如图1所示,前述技术方案中所涉及的目标天体可以为火星、金星、木星、土星等太阳系内的天体,以下以火星作为实施例,说明本发明的具体实施过程:
1.建立深空探测器基于太阳和八大行星引力轨道动力学的状态模型
首先初始化位置、速度,设X=[x y z vx vy vz]T为在日心惯性坐标系中的状态向量,x,y,z,vx,vy,vz分别为探测器在日心惯性坐标系中三轴的位置和速度,X′=[x′ y′ z′ vx′ vy′ vz′]T为在火心惯性坐标系中的状态向量,x′,y′,z′,vx′,v′y,vz′分别为探测器在火心惯性坐标系中三轴的位置和速度,上述变量都是与t有关的函数,根据探测器的轨道设计,选取探测器的位置和速度初值为X(0)和X′(0);其次初始化火星星历误差初值为Xerr(0)=[xerr yerr zerr vxerr vyerr vzerr]T,xerr,yerr,zerr,vxerr,vyerr,vzerr分别为日心坐标系中火星星历三轴的位置误差和速度误差。
A.在火心惯性坐标系中建立深空探测器基于太阳和八大行星引力轨道动力学的第一状态模型,即天文导航子***的状态模型;
考虑太阳和火星、地球等太阳和八大行星对探测器的引力作用,选取火心惯性坐标系,可得深空探测器在火心惯性坐标系中第一状态模型,即天文导航子***的状态模型为:
式中,x′,y′,z′为探测器在火心惯性坐标系中三轴位置,vx′,v′y,vz′为探测器在火心惯性坐标系中三轴速度,为探测器在火心惯性坐标系中三轴位置的微分,为探测器在火心惯性坐标系中三轴速度的微分,μs、μm和μic分别为太阳、火星和第ic颗行星的引力常数;r′ps为日心到探测器的距离;r′pm为火星到探测器的距离;r′ms为日心到火心的距离;为第ic颗行星到探测器的距离;r′mi为第ic颗行星质心到火心的距离;(xs′,ys′,zs′),分别为太阳、第ic颗行星在火心惯性坐标系中的三轴位置坐标,可根据时间由行星星历表获得,wx′,wy′,wz′分别为第一状态模型中探测器三轴的状态模型误差;ic表示太阳和八大行星中从内至外的第ic颗行星,ic=1,2,3...,N(ic≠4),N=8,由于ic=4表示第4颗行星(火星),因此不包含在求和表达式中。
式(1)中的各变量都是与时间t有关的变量,可简写为
式中,X′(t)=[x′,y′,z′,vx′,v′y,vz′]T为第一状态模型的状态向量,为X′(t)的微分,f1(X′(t),t)为第一状态模型的***非线性连续状态转移函数,w′(t)=[wx′,w′y,wz′]T为第一状态模型的状态模型误差。
B.在日心惯性坐标系中建立深空探测器基于太阳和八大行星引力轨道动力学的第二状态模型,即无线电导航子***的状态模型;
考虑太阳和火星、地球等太阳和八大行星对探测器的引力作用,选取日心惯性坐标系,可得深空探测器在日心惯性坐标系中展开为分量形式的第二状态模型,即无线电导航子***的状态模型为:
式中,x,y,z为探测器在日心惯性坐标系中三轴位置,vx,vy,vz为探测器在日心惯性坐标系中三轴速度,为探测器在日心惯性坐标系中三轴位置的微分,为探测器在日心惯性坐标系中三轴速度的微分,μs为太阳引力常数,为第ic个行星的引力常数;rps为日心到探测器的距离;为第ic个行星到探测器的距离;为第ic个行星质心到日心的距离;分别为第ic个行星在日心惯性坐标系中的坐标,可根据时间由行星星历表获得,wx,wy,wz分别为第二状态模型中探测器三轴的状态模型误差;
式(3)中的各变量都是与时间t有关的变量,可简写为:
式中,X(t)=[x,y,z,vx,vy,vz]T为第二状态模型的状态向量,为X(t)的微分,f2(X(t),t)为第二状态模型***非线性连续状态转移函数,w(t)=[wx,wy,wz]T为第二状态模型的状态模型误差。
2.分别建立天文导航子***和无线电导航子***量测模型;
(1)天文导航子***量测模型
火星与第i颗背景恒星之间角度信息θmi的大小在不同坐标系中是相同的,因此其表达式为:
式中,为在火星敏感器测量坐标系中从火星至探测器的单位矢量,为在惯性坐标系中火星到探测器的单位矢量,可表示为:
式中,(x′,y′,z′)为探测器在火心惯性坐标系中三轴位置坐标,为火星图像中第i颗背景恒星在火星敏感器测量坐标系中的单位矢量,为在惯性坐标系中第i颗背景恒星星光的单位矢量,i=1,2,3,可由恒星星历数据库直接得到,
同理可得火卫一和火卫二与其第i颗背景恒星之间角度信息θpi和θdi的表达式为:
设天文导航子***量测量Z1=[θm1,θm2,θm3,θp1,θp2,θp3,θd1,θd2,θd3]T,天文导航子***量测噪声 分别为测量θm1,θm2,θm3,θp1,θp2,θp3,θd1,θd2,θd3的观测误差,由于各变量都是与时间t有关的变量,则可根据式(6)~式(8)建立天文导航子***量测模型的表达式为
Z1(t)=h1[X′(t),t]+v1(t) (9)
式中,h1[X′(t),t]为天文导航子***非线性连续量测函数。
(2)无线电导航子***量测模型
无线电测控站可根据无线电信号的时延与多普勒频移测量探测器到测控站的距离与距离变化率,一般主测控站附近还配有多个副站,根据多个测控站测量的距离与距离变化率可以获得探测器相对于地心的位置与速度,无线电测控导航所获得的距离与距离变化率表达式为:
式中,Ri为第i个测控站的测距量测信息(探测器与测控站之间的距离),为第i个测控站的距离变化率量测信息(探测器与测控站之间的距离变化率),i=1,2,3,(xf,yf,zf)为地面测控站在地心惯性坐标系中的位置坐标。
设无线电导航子***量测量无线电导航子***量测噪声 分别为三个测控站所测量得到的探测器与测控站之间的距离和距离变化率,分别为测量的观测误差,由于各变量都是与时间t有关的变量,则可根据式(10)和式(11)建立无线电导航子***量测模型的表达式为:
Z2(t)=h2[X(t),t]+v2(t) (12)
式中,h2[X(t),t]为无线电导航子***非线性连续量测函数。
3.对步骤1和步骤2中的状态模型和量测模型进行离散化
X′(k)=F1(X′(k-1),k-1)+W′(k-1) (13)
X(k)=F2(X(k-1),k-1)+W(k-1) (14)
Z′(k)=H1(X′(k),k)+V1(k) (15)
Z(k)=H2(X(k),k)+V2(k) (16)
式中,k=1,2,…,F1(X′(k-1),k-1)和F2(X(k-1),k-1)分别为f1(X′(t),t)和f2(X(t),t)离散后从第k-1时刻到第k时刻的非线性状态转移函数,H1(X′(k),k)和H2(X(k),k)分别为h1(X′(t),t)和h2(X(t),t)离散后第k时刻的非线性量测函数,W′(k),W(k),V1(k),V2(k)分别为w′(t),w(t),v1(t)和v2(t)离散后第k时刻的等效噪声,且W′(k)和V1(k)、W(k)和V2(k)互不相关。
4.天文导航和无线电导航量测量的获取及处理
(1)天文导航子***量测量的获取及处理
图2给出了天文导航子***所用火星敏感器的成像原理图,火卫一敏感器、火卫二敏感器成像过程与之类似。火星敏感器主要由光学透镜和二维成像面阵组成,在敏感器测量坐标系O′XcYcZc中沿火星到探测器的矢量方向火星反射太阳光线射向火星敏感器,此时,火星在火星敏感器测量坐标系中的坐标为(xc,yc,zc);火星敏感器的光学透镜以焦距f将火星的光线折射后成像在二维成像面阵上,二维成像面阵将照在每个成像单元上的图像亮度信号储存。根据敏感器的成像原理,天文导航子***量测量的处理过程如下所示:
A.天体图像的处理
由于火星在二维成像面阵上的图像并不是一个点,而是一个圆,通过质心识别等图像处理技术确定火星图像在二维成像平面坐标系OX2dY2d的质心(x2d,y2d),这一中心可以用像元像线坐标系OplXplYpl中的像元像线(p,l)表示。
B.质心坐标从像元像线坐标系转换至二维成像平面坐标系的转换
根据像元像线坐标系与二维成像平面坐标系之间的关系,将火星质心坐标从像元像线坐标系转换至二维成像平面坐标系:
式中,(x2d,y2d)为火星在二维成像平面OX2dY2d中的坐标,p和l分别为火星在火星敏感器二维成像平面的像元和像线, 为由毫米转为像素的敏感器转换矩阵,p0和l0分别为火星敏感器中心在像元像线坐标系OXplYpl中的像元和像线。
C.质心坐标从二维成像平面坐标系转换至敏感器测量坐标系的转换
根据透镜成像原理,将火星质心坐标从二维成像平面坐标系转换至敏感器测量坐标系:
式中,f为火星敏感器的焦距,为在火星敏感器测量坐标系中从火星至探测器的单位矢量。
同理可获得火星图像中第i颗背景恒星在火星敏感器测量坐标系中的单位矢量i=1,2,3。
D.计算星光角距
根据在火星敏感器测量坐标系中火星至探测器的单位矢量与火星图像中第i颗背景恒星的单位矢量计算两个矢量之间的星光角距θmi
同理可获得火卫一与其背景恒星、火卫二与其背景恒星之间的星光角距θpi,θdi。
(2)无线电导航子***量测量的获取及处理
选取探测器与测控站之间的距离变化率和探测器与测控站之间的距离R作为无线电导航子***的量测量。
探测器与测控站之间的距离变化率是通过测量探测器接收的载波信号频率f′rec与地面站发射的载波信号频率f0相比的频率偏移来获得的,具体可表示为:
式中,为探测器与地面站间的相对距离变化率,c为光在真空中的传播速度,f0为地面站发射的无线电信号的固有频率,f′rec为探测器上的接收机接收到的无线电信号的频率,δfatm为大气层对信号的时延,δf0为由信号源本振频率的漂移引起的误差,由于目前地面站多采用USO(Ultra Stable Oscillators),该误差的量级很小。
探测器与测控站之间的距离R是通过测量无线电信号由地面测控站发射至探测器再返回的时间延迟而获得的,具体可表示为:
R=cτ/2+cδtR-cδtT+δρatm (21)
式中,c为光在真空中的传播速度,τ为无线电信号由地面测控站发射至探测器再返回的时间延迟,δtR为接收装置时钟同步误差,δtT为发射装置时钟同步误差,δρatm为由原子钟差引起的距离测量误差。
5.对天文导航子***进行Unscented卡尔曼滤波
根据第一状态模型式(13)、天文导航子***量测模型式(15)、由天文导航敏感器获得的量测量式(19),进行天文导航子***Unscented卡尔曼滤波。具体步骤如下:
A.初始化
式中,为第0时刻(初始时刻)在火心惯性坐标系中探测器的三轴位置和速度估计值,x0′为第0时刻(初始时刻)在火心惯性坐标系中探测器的三轴位置和速度真实值,P0′为状态向量的初始均方误差阵。
B.计算采样点
在天文导航子***第k-1时刻状态向量附近选取一系列样本点,这些样本点的均值和均方误差阵分别为和设状态向量为6×1维,那么13个天文导航子***的样本点χ′0,k,...,χ′12,k及其权重W0′…W′12分别如下:
式中,当P′k-1=A′TA′时,取A′的第j行,当P′k-1=A′A′T时,取A′的第j列,得第k-1时刻采样点χ′k-1的统一表达式为:
C.时间更新
天文导航子***状态向量的一步预测χ′k|k-1为:
χ′k|k-1=F1(χ′k-1,k-1) (25)
天文导航子***所有采样点状态向量的一步预测加权后结果为:
式中,Wj′为第j个采样点的权值;
天文导航子***状态向量的估计均方误差阵一步预测Pk′-为:
式中,Qk′为k时刻天文导航子***的状态模型误差协方差阵;
天文导航子***采样点对应的量测估计向量Z′k|k-1
Z′k|k-1=H1(χ′k|k-1,k) (28)
天文导航子***所有采样点量测估计加权向量
D.量测更新
天文导航子***量测均方误差阵为:
式中,Rk′为k时刻天文导航子***的量测噪声协方差阵;
天文导航子***状态向量量测量均方误差阵
天文导航子***滤波增益Kk′为:
天文导航子***的估计状态向量和估计均方误差阵Pk′为
6.对无线电导航子***进行Unscented卡尔曼滤波
根据第二状态模型式(14)、无线电导航子***量测模型式(16)、由无线电接收装置获得的量测量式(20)和式(21),进行无线电导航子***Unscented卡尔曼滤波。具体步骤如下:
A.初始化
式中,为第0时刻(初始时刻)无线电导航子***在日心惯性坐标系中探测器的三轴位置和速度估计值,x0为第0时刻(初始时刻)无线电导航子***在日心惯性坐标系中探测器的三轴位置和速度真实值,P0为无线电导航子***状态向量的初始均方误差阵。
B.计算采样点
在无线电导航子***第k-1时刻状态向量附近选取一系列样本点,这些样本点的均值和均方误差阵分别为和Pk-1。设状态向量为6×1维,那么13个无线电导航子***的样本点χ0,k,...,χ12,k及其权重W0…W12分别如下:
式中,当Pk-1=ATA时,取A的第j行,当Pk-1=AAT时,取A的第j列,得第k-1时刻采样点χk-1的统一表达式为:
C.时间更新
无线电导航子***状态向量的一步预测χk|k-1为:
χk|k-1=F2(χk-1,k-1) (38)
无线电导航子***所有采样点状态向量的一步预测加权后结果为:
式中,Wj为第j个采样点的权值;
无线电导航子***状态向量的估计均方误差阵一步预测为:
式中,Qk为k时刻无线电导航子***的状态模型误差协方差阵;
无线电导航子***采样点对应的量测估计向量Zk|k-1
Zk|k-1=H2(χk|k-1,k) (41)
无线电导航子***所有采样点量测估计加权向量
D.量测更新
无线电导航子***量测均方误差阵为:
式中,Rk为无线电导航子***量测噪声协方差阵;
无线电导航子***状态向量量测量均方误差阵为:
无线电导航子***滤波增益Kk为:
无线电导航子***估计状态向量和估计均方误差阵Pk为:
7.判定是否需要进行火星星历校正
当天文导航子***Unscented卡尔曼滤波的估计均方误差阵Pk′大于已有火星星历误差均方误差阵时,即Pk′>Peph,则不进行星历校正,直接进行第9步信息融合;当Pk′<Peph时,则进行星历误差校正,执行第8步;
8.对火星星历误差进行建模、估计并反馈校正
无线电测控导航主要的观测量为探测器相对于地面站的距离与距离变化率等,这种导航方式可以直接测量相对于太阳的导航信息,天文导航可以直接测量相对于目标天体(如火星)的导航信息。由于目标天体的星历存在误差(200m~100km),而无线电导航方法仅能获取相对太阳的高精度导航信息,因此由该方法获得的相对目标天体导航信息精度低;虽然天文导航方法可以获得较为准确的相对目标天体导航信息,因此由该方法无法获得相对太阳高精度的导航信息。因此需要对目标天体星历误差进行估计,并反馈校正由目标天体星历误差引起的导航误差。
天文导航子***获得的相对于火星的位置和速度为天文导航子***获得的相对于太阳的位置和速度为为火星相对于太阳的位置和速度,可从天体星历数据库获得);无线电导航子***获得的相对于太阳的位置和速度为无线电导航子***获得的相对于火星的位置和速度为由此可以看出,单一导航***会受到火星星历误差的影响,无法同时满足探测器相对太阳和火星高精度导航的需求。因此可利用天文导航子***Unscented卡尔曼滤波的结果和无线电导航子***Unscented卡尔曼滤波的结果对火星星历误差进行估计,具体步骤如下:
A.建立火星星历误差状态模型
考虑到捕获段持续时间短(约40天)这一特点,火星星历误差变化较小,将捕获段内火星的星历误差特性考虑为常值误差,在日心惯性坐标系中建立火星星历误差状态模型为:
式中,为日心惯性坐标系中火星星历的三轴位置误差的微分,离散化后简写为:
Xerr(k)=Ferr(Xerr(k-1),k-1)+Werr(k-1) (49)
式中,状态转移函数Ferr(Xerr(k-1),k-1)=Φerr,k,k-1Xerr,k-1,Φerr,k,k-1为第k-1时刻到第k时刻的状态转移矩阵,Xerr(k)为第k时刻火星星历误差状态向量,且Xerr(k)=Xerr,k,Werr(k-1)为第k-1时刻火星误差状态模型误差。
B.建立火星星历误差量测模型
因此火星星历误差的量测模型可以表示为:
Zerr=H3(Xerr(k),k)+V3 (50)
式中,H3(Xerr(k),k)为k时刻的量测函数,V3为火星星历误差量测噪声。
C.获取火星星历误差量测量
火星星历误差量测量Zerr可以表示为:
式中,为无线电导航子***获得的相对于太阳的位置和速度,为天文导航子***获得的相对于火星的位置和速度,为火星相对于太阳的位置和速度,从天体星历数据库中获得。
D.对火星星历误差进行卡尔曼滤波估计
根据步骤A和步骤B建立的火星星历误差状态模型式(49)和量测模型式(50),以及步骤C所获取的火星星历误差量测量式(51),利用卡尔曼滤波方法,对火星星历误差进行估计,获得火星星历误差估计状态向量与估计均方误差阵,具体如下:
火星星历误差状态向量的一步预测
式中,为k-1时刻火星星历误差一步预测状态向量。
火星星历误差状态向量的估计均方误差阵一步预测
Perr,k/k-1=Φerr,k,k-1Perr,k-1Φerr,k,k-1 T+Qerr,k (53)
式中,Perr,k-1为k-1时刻火星星历误差状态向量的估计均方误差阵,Qerr,k为k时刻火星星历误差状态模型误差均方误差阵。
卡尔曼滤波增益
Kerr,k=Perr,k/k-1Herr,k T(Herr,kPerr,k/k-1Herr,k T+Rerr,k)-1 (54)
式中,Herr,k为k时刻火星星历误差量测矩阵,Herr,kXerr,k=H3(Xerr,k),Rerr,k为k时刻火星星历误差量测模型误差协方差阵。
火星星历误差估计状态向量
式中,Zerr,k为k时刻火星星历误差量测量。
火星星历误差估计均方误差阵
Perr,k=(I-Kerr,kHerr,k)Perr,k/k-1 (56)
式中,I为单位阵。
E.对火星星历误差进行反馈校正
将步骤D中获得的火星星历误差和火星星历估计均方误差阵Perr,k反馈回深空探测器的第一状态模型和第二状态模型中,并重新确定第一状态模型和第二状态模型的模型误差协方差阵Qk′和Qk,最后将校正后的模型误差协方差阵Qk′和Qk输入至天文导航子***Unscented卡尔曼滤波和无线电导航子***Unscented卡尔曼滤波中,修正下一时刻的导航结果。
9.信息融合
当探测器不在无线电测控范围内时,无线电导航子***没有输入的无线电导航量测量,此时对天文导航子***进行Unscented卡尔曼滤波,包括时间更新和量测更新(第5步的步骤C和步骤D)等步骤,无线电导航子***只进行时间更新(第6步的步骤C);当探测器处于无线电测控范围内时,无线电导航子***有输入的无线电导航量测量,对两个子***同时进行Unscented卡尔曼滤波,都进行时间更新和量测更新。
在滤波过程中得到的两个局部估计状态向量 两个估计均方误差阵 按下式进行融合,得到全局的估计状态向量和全局估计均方误差阵分别为:
将全局估计结果反馈给两个导航子***,作为k时刻两个导航子***的估计结果:
式中,为信息分配因子。
信息分配因子选择的基本原则是在满足信息守恒公式的前提下与局部导航滤波精度成正比,为了使导航***具有更强的自适应能力和容错能力,使用基于估计均方误差阵范数的动态分配信息因子的算法。
令
式中,||·||F为Frobenius范数,即对于任意矩阵A,
最终将式(59)和式(60)获得的第k时刻在火心惯性坐标系中和在日心惯性坐标系中的估计状态向量和估计均方误差阵P1(k),P2(k)输出,估计状态向量分别表示在火心惯性坐标系中和日心惯性坐标系中探测器的位置、速度信息,输出的估计均方误差阵P1(k),P2(k)表示了滤波估计的性能,并将这些导航信息分别返回天文导航子***和无线电导航子***中,用于k+1时刻的位置、速度导航信息,k=1,2,...。
本发明说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
Claims (1)
1.一种基于星历修正的深空探测器天文与无线电组合导航方法,其特征在于:首先建立深空探测器的状态模型和量测模型,利用无线电导航子***和天文导航子***分别获取相对太阳和相对目标天体的量测量,通过Unscented卡尔曼滤波估计得到探测器在日心和目标天体为中心惯性坐标系中的位置和速度;在此基础上,建立目标天体星历误差的状态模型和量测模型,并由天文和无线电两个子***的导航信息获得目标天体星历误差量测量,利用卡尔曼滤波方法对目标天体星历误差进行估计,并将目标天体星历误差反馈至导航***模型中,对***模型进行修正,此后将两个子***的结果进行信息融合,获得校正星历误差后的相对于目标天体和相对于日心的探测器位置和速度;具体包括以下步骤:
①建立深空探测器基于太阳和八大行星引力轨道动力学状态模型;
A.在目标天体为中心的惯性坐标系中,建立深空探测器基于太阳和八大行星引力轨道动力学的第一状态模型,即天文导航子***的状态模型;
B.在日心惯性坐标系中,建立深空探测器基于太阳和八大行星引力轨道动力学的第二状态模型,即无线电导航子***的状态模型;
②分别建立天文导航子***和无线电导航子***量测模型;
③对步骤①和步骤②中的状态模型和量测模型进行离散化;
④天文导航和无线电导航量测量的获取及处理;
⑤对天文导航子***进行Unscented卡尔曼滤波;
根据第一状态模型、天文导航子***量测模型、天文导航敏感器获得的量测量,进行天文导航子***Unscented卡尔曼滤波,获得在目标天体惯性坐标系中深空探测器位置和速度的估计状态向量和估计均方误差阵P′k;
⑥对无线电导航子***进行Unscented卡尔曼滤波;
根据第二状态模型、无线电导航子***量测模型、由无线电接收装置获得的量测量,进行无线电导航子***Unscented卡尔曼滤波,获得在日心惯性坐标系中深空探测器位置和速度的估计状态向量与估计均方误差阵Pk;
⑦判定是否需要进行目标天体星历校正;
当天文导航子***Unscented卡尔曼滤波的估计均方误差阵P′k大于已有目标天体星历误差均方差阵Peph时,即P′k>Peph,则不进行星历校正,直接进行第⑨步信息融合;当P′k<Peph时,进行目标天体星历误差校正,执行步骤⑧;
⑧对目标天体星历误差进行建模、估计并反馈校正;
A.建立目标天体星历误差状态模型
在日心惯性坐标系中建立目标天体星历误差状态模型为:
式中, 为日心惯性坐标系中目标天体星历的三轴位置误差的微分,离散化后为:
Xerr(k)=Ferr(Xerr(k-1),k-1)+Werr(k-1)
式中,状态转移函数Ferr(Xerr(k-1),k-1)=Φerr,k,k-1Xerr,k-1,其中Φerr,k,k-1为第k-1时刻到第k时刻的状态转移矩阵,Xerr(k)为第k时刻目标天体星历误差状态向量,且Xerr(k)=Xerr,k,Werr(k-1)为第k-1时刻目标天体星历误差状态模型误差,k=1,2,...;
B.建立目标天体星历误差量测模型
建立目标天体星历误差的量测模型为:
Zerr=H3(Xerr(k),k)+V3
式中,H3(Xerr(k),k)为k时刻的量测函数,V3为目标天体星历误差量测噪声;
C.获取目标天体星历误差量测量
目标天体星历误差量测量Zerr表示为:
式中,为无线电导航子***获得的相对于太阳的位置和速度,为天文导航子***获得的相对于目标天体的位置和速度,为目标天体相对于太阳的位置和速度,从天体星历数据库中获得;
D.对目标天体星历误差进行Unscented卡尔曼滤波估计
根据目标天体星历误差的状态模型、量测模型和获取的目标天体星历误差量测量利用Unscented卡尔曼滤波方法,对目标天体星历误差进行估计,获得目标天体星历误差估计状态向量与估计均方误差阵;
E.对目标天体星历误差进行反馈校正
将步骤D中获得的目标天体星历误差估计状态向量及估计均方误差阵反馈回深空探测器的第一状态模型和第二状态模型中,并重新确定第一状态模型和第二状态模型的模型误差协方差阵,最后将校正目标天体星历后的状态模型和模型误差协方差阵输入至天文导航子***Unscented卡尔曼滤波和无线电导航子***Unscented卡尔曼滤波中,修正下一时刻的导航结果;
⑨信息融合;
利用星历修正***获得的星历误差,将无线电导航***的导航信息转换至以目标天体为中心的惯性坐标系中,与天文导航***的导航信息进行融合,当探测器不在无线电测控范围内时,无线电导航子***没有输入的无线电导航量测量,此时对天文导航子***进行Unscented卡尔曼滤波,无线电导航子***只利用第二状态模型进行时间更新;当探测器处于无线电测控范围内时,无线电导航子***有输入的无线电导航量测量,对两个导航子***同时进行Unscented卡尔曼滤波;
最终输出第k时刻在目标天体为中心的惯性坐标系中表示探测器位置和速度的估计状态向量和估计均方误差阵,并根据修正后的目标行星星历,将该结果转换至日心惯性坐标系中,输出在日心惯性坐标系中表示探测器位置和速度的估计状态向量和估计均方误差阵,将这些导航信息分别返回天文导航子***和无线电导航子***中,用于k+1时刻的位置、速度导航信息的估计,k=1,2,...。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310006575.XA CN103063217B (zh) | 2013-01-08 | 2013-01-08 | 一种基于星历修正的深空探测器天文/无线电组合导航方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310006575.XA CN103063217B (zh) | 2013-01-08 | 2013-01-08 | 一种基于星历修正的深空探测器天文/无线电组合导航方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103063217A CN103063217A (zh) | 2013-04-24 |
CN103063217B true CN103063217B (zh) | 2015-04-29 |
Family
ID=48105946
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310006575.XA Active CN103063217B (zh) | 2013-01-08 | 2013-01-08 | 一种基于星历修正的深空探测器天文/无线电组合导航方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103063217B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103344245B (zh) * | 2013-07-09 | 2015-11-18 | 北京航空航天大学 | 火星进入段imu和甚高频无线电组合导航的ud-skf方法 |
US10262272B2 (en) * | 2014-12-07 | 2019-04-16 | Microsoft Technology Licensing, Llc | Active machine learning |
CN104501804B (zh) * | 2014-12-17 | 2017-06-13 | 深圳航天东方红海特卫星有限公司 | 一种基于gps测量数据的卫星在轨轨道预报方法 |
CN104764449B (zh) * | 2015-04-23 | 2017-07-11 | 北京航空航天大学 | 一种基于星历修正的捕获段深空探测器自主天文导航方法 |
CN105301958B (zh) * | 2015-11-03 | 2018-01-02 | 北京理工大学 | 一种基于气动力辅助的平衡点周期轨道捕获方法 |
CN105509750B (zh) * | 2015-11-27 | 2018-02-06 | 上海卫星工程研究所 | 一种天文测速与地面无线电组合的火星捕获段导航方法 |
CN107024212B (zh) * | 2017-06-22 | 2019-11-22 | 北京航空航天大学 | 一种深空探测器x射线脉冲星/时间差分天文多普勒组合导航方法 |
CN108413969B (zh) * | 2018-01-31 | 2021-02-09 | 上海航天控制技术研究所 | 一种采用卫星影像辅助和无线通信网络的定位方法 |
CN110376618B (zh) * | 2019-08-30 | 2020-08-28 | 北京航天宏图信息技术股份有限公司 | 基于北斗三号卫星星基增强的定位方法、装置及终端 |
CN110672105B (zh) * | 2019-11-22 | 2021-04-20 | 北京理工大学 | 一种小天体接近段双探测器高精度协同光学导航方法 |
CN110940333B (zh) * | 2019-12-12 | 2021-07-13 | 中南大学 | 一种基于在线估计的深空探测器测角及时延组合导航方法 |
CN111024121B (zh) * | 2019-12-13 | 2023-03-31 | 中国科学院光电技术研究所 | 一种光电设备自主精度鉴定的***和方法 |
CN111323020B (zh) * | 2020-02-25 | 2021-11-16 | 上海航天控制技术研究所 | 一种基于火星边缘及中心多矢量观测的自主定轨方法 |
CN111637896B (zh) * | 2020-05-12 | 2021-12-07 | 北京控制工程研究所 | 一种基于星历约束辅助的自主天文导航方法 |
CN117055074B (zh) * | 2023-10-13 | 2024-01-23 | 中国电子科技集团公司第十五研究所 | 相对精度综合量化评估方法、服务器及存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101270993A (zh) * | 2007-12-12 | 2008-09-24 | 北京航空航天大学 | 一种远程高精度自主组合导航定位方法 |
CN101285687A (zh) * | 2008-05-30 | 2008-10-15 | 中国科学院上海技术物理研究所 | 地空天一体化自主导航***设计方法 |
CN102168981A (zh) * | 2011-01-13 | 2011-08-31 | 北京航空航天大学 | 一种深空探测器火星捕获段自主天文导航方法 |
CN102168980A (zh) * | 2011-01-13 | 2011-08-31 | 北京航空航天大学 | 一种基于小行星交会的深空探测器自主天文导航方法 |
CN102607564A (zh) * | 2012-03-09 | 2012-07-25 | 北京航空航天大学 | 一种基于星光/地磁组合信息的小卫星自主导航***及其导航方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8825399B2 (en) * | 2008-07-24 | 2014-09-02 | Raytheon Company | System and method of passive and autonomous navigation of space vehicles using an extended Kalman filter |
-
2013
- 2013-01-08 CN CN201310006575.XA patent/CN103063217B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101270993A (zh) * | 2007-12-12 | 2008-09-24 | 北京航空航天大学 | 一种远程高精度自主组合导航定位方法 |
CN101285687A (zh) * | 2008-05-30 | 2008-10-15 | 中国科学院上海技术物理研究所 | 地空天一体化自主导航***设计方法 |
CN102168981A (zh) * | 2011-01-13 | 2011-08-31 | 北京航空航天大学 | 一种深空探测器火星捕获段自主天文导航方法 |
CN102168980A (zh) * | 2011-01-13 | 2011-08-31 | 北京航空航天大学 | 一种基于小行星交会的深空探测器自主天文导航方法 |
CN102607564A (zh) * | 2012-03-09 | 2012-07-25 | 北京航空航天大学 | 一种基于星光/地磁组合信息的小卫星自主导航***及其导航方法 |
Non-Patent Citations (1)
Title |
---|
基于UKF的航天器多普勒/天文组合导航方法研究;宋伟等;《载人航天》;20071231;第20-28页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103063217A (zh) | 2013-04-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103063217B (zh) | 一种基于星历修正的深空探测器天文/无线电组合导航方法 | |
CN105203101A (zh) | 一种基于目标天体星历修正的深空探测器捕获段天文导航方法 | |
Williams et al. | OSIRIS-REx flight dynamics and navigation design | |
CN102175241B (zh) | 一种火星探测器巡航段自主天文导航方法 | |
CN102168981B (zh) | 一种深空探测器火星捕获段自主天文导航方法 | |
CN102168980B (zh) | 一种基于小行星交会的深空探测器自主天文导航方法 | |
CN103674032B (zh) | 融合脉冲星辐射矢量和计时观测的卫星自主导航***及方法 | |
CN104764449B (zh) | 一种基于星历修正的捕获段深空探测器自主天文导航方法 | |
CN104298647B (zh) | 基于低轨道地球卫星的地影时刻预报的星上确定方法 | |
Cui et al. | X-ray pulsars/Doppler integrated navigation for Mars final approach | |
US6622970B2 (en) | Method and apparatus for autonomous solar navigation | |
CN103148856B (zh) | 一种基于自适应尺度变化的借力飞行探测器自主天文导航方法 | |
CN103968834B (zh) | 一种近地停泊轨道上深空探测器的自主天文导航方法 | |
Franzese et al. | Deep-space optical navigation for M-ARGO mission | |
Kai et al. | Autonomous navigation for a group of satellites with star sensors and inter-satellite links | |
Liu et al. | X-ray pulsar/Doppler difference integrated navigation for deep space exploration with unstable solar spectrum | |
CN106643741A (zh) | 一种卫星相对小行星视觉自主导航方法 | |
CN107144283A (zh) | 一种用于深空探测器的高可观度光学脉冲星混合导航方法 | |
CN104457705A (zh) | 基于天基自主光学观测的深空目标天体初定轨方法 | |
CN105547303A (zh) | 一种平动点星座的自主导航方法 | |
Wang et al. | Absolute navigation for Mars final approach using relative measurements of X-ray pulsars and Mars orbiter | |
CN103017772A (zh) | 一种基于可观性分析的光学和脉冲星融合自主导航方法 | |
Feng et al. | Autonomous orbit determination and its error analysis for deep space using X-ray pulsar | |
Bu et al. | A novel interplanetary optical navigation algorithm based on Earth–Moon group photos by Chang’e-5T1 probe | |
Guinn et al. | The deep-space positioning system concept: Automating complex navigation operations beyond the earth |
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 |