首页 / 专利库 / 光学 / 多普勒光谱 / 多普勒差分测速模型及与X射线脉冲星组合导航的方法

多普勒差分测速模型及与X射线脉冲星组合导航的方法

阅读:391发布:2020-05-17

专利汇可以提供多普勒差分测速模型及与X射线脉冲星组合导航的方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了一种多普勒差分测速模型及与 X射线 脉冲星组合导航的方法,其多普勒差分测速模型为Ys(t)=hs(X,t)+ωv,其中,ωv为多普勒测量噪声,Ys和hs分别表示测量值和测量方程,Ys=vD2-vD1,组合导航的方法包括以下步骤:步骤1:建立深空探测器的轨道动 力 学模型;步骤2:建立X射线脉冲星导航测距模型;步骤3:建立多普勒差分测速模型;步骤4:利用 扩展卡 尔曼 滤波器 滤波。本发明与 现有技术 相比的优点在于:(1)与现有脉冲星/多普勒组合导航方法相比,本发明在不稳定的太阳 光谱 下仍能正常稳定工作;(2)与脉冲星导航方法和多普勒导航方法相比,本发明充分利用了多种导航信息,具有完全可观测性,能提供更高 精度 的 定位 信息。,下面是多普勒差分测速模型及与X射线脉冲星组合导航的方法专利的具体信息内容。

1.一种多普勒差分测速模型,其特征在于:
首先假设深空探测器在t1时刻和t时刻的位置分别为r1和r,火星在t1时刻的位置和速度分别为rM和vM,深空探测器上搭载两个光谱仪,第一光谱仪在t1时刻测量太阳直射光的光谱频移,第二光谱仪在t时刻测量火星反射太阳光的光谱频移,且t1和t的关系如下所示:
(式壹)
其中,c为光速;
基于上述假设,多普勒差分测速模型的建立包括以下步骤:
步骤A1:在t1时刻之前,根据t和式壹确定t1的值;
步骤A2:在t1时刻,利用第一光谱仪获取太阳径向上的多普勒速度vD1;在t时刻,利用第二光谱仪获取火星径向上的多普勒速度vD2;
步骤A3:建立多普勒差分测速模型:
s s
Y(t)=h(X,t)+ωv (式贰)
s s s
其中,ωv为多普勒测量噪声,Y 和h 分别表示测量值和测量方程,Y =vD2-vD1,
2.一种利用权利要求1所述的多普勒差分测速模型与X射线脉冲星组合导航的方法,其特征在于,包括以下步骤:
步骤B1:建立深空探测器的轨道动学模型;
步骤B2:建立X射线脉冲星导航测距模型;
步骤B3:建立多普勒差分测速模型;
步骤B4:利用扩展卡尔曼滤波器滤波。
3.根据权利要求2所述的多普勒差分测速模型与X射线脉冲星组合导航的方法,其特征在于:步骤B1中所述的建立深空探测器的轨道动力学模型,其具体实现过程为:
因为深空探测器的状态矢量X为:
(式叁)
其中,r=[x,y,z]T和v=[vx,vy,vz]T分别为深空探测器的位置和速度矢量,x,y,z分别为深空探测器的位置在三轴上的分量,vx,vy,vz分别为深空探测器的速度在三轴上的分量;
则深空探测器的轨道动力学模型为:
(式肆)
其中, 分别为x,y,z,vx,vy,vz的导数,
式肆可表示为:
(式伍)
其中, 是X的导数, 为时刻t的 ,f(X,t)为深空探测器的状态
转移模型,[x1,y1,z1]和[x2,y2,z2]分别是火星和地球相对于太阳系质心的相对位置矢量,μs,μm,μe分别是太阳,火星和地球的引力常数,rps,rpm,rpe分别是深空探测器到太阳质心,火星质心以及地球质心之间的距离,其计算公式 为 :
分别是火星
质心、地球质心分别到太阳质心之间的距离;深空探测器的导航系统噪声ω=[0,0,0,T
ΔFx,ΔFy,ΔFz],其中,ΔFx,ΔFy和ΔFz是摄动力,ω(t)为时刻t深空探测器的导航系统噪声。
4.根据权利要求3所述的多普勒差分测速模型与X射线脉冲星组合导航的方法,其特征在于:步骤B2中所述的建立X射线脉冲星导航测距模型,其具体实现包括以下子步骤:
步骤2.1:确定用于导航的X射线脉冲星数量及其方位参数;
步骤2.2:建立X射线脉冲到达时间转换模型为:
(式
陆)
(式柒)
其中, 为深空探测器与太阳系质心之间的距离在第j颗脉冲星方向上的投j j j
影;n 是第j颗脉冲星的方向矢量,j=1,2,…I,I为采用的X射线脉冲星数量,α 和δ分别为第j颗脉冲星的赤经和赤纬, 和 分别为脉冲到达深空探测器时间与到达太阳系质心时间,c为光速, 为第j颗脉冲星到太阳系质心的距离,b为太阳系质心相对于太阳的位置矢量,|b|为位置矢量b的长度;μSun为太阳引力常数,r是深空探测器相对于太阳系质心的位置矢量,|r|为位置矢量r的长度;
X
步骤2.3:确定X射线脉冲星导航观测模型h(X(t),t);
X
首先假设X射线脉冲星导航观测量Y 为:
(式捌)
其对应的量测噪声为V,则X射线脉冲星导航观测模型可表示为:
YX=hX(X(t),t)+V(t) (式玖)
其中,V(t)为时刻t的量测噪声,观测模型hX(X(t),t)如下:
(式拾)
其中,第j颗脉冲星的相应项hj(X(t),t)的表达式如下:

(式十一)。
5.根据权利要求4所述的多普勒差分测速模型与X射线脉冲星组合导航的方法,其特征在于:步骤B4中所述的利用扩展卡尔曼滤波器滤波,其具体实现过程为:首先建立脉冲星中滤波器的测量模型,在脉冲星观测期间,其测量模型h(X,t)和测量值Y为:
s
h(X,t)=h(X,t) (式十二)
s
Y=Y (式十三)
一旦获得了脉冲到达时间,测量模型h(X,t)和测量值Y则表示为:
(式十四)
(式十五)
其中,hs(X(t),t)即hs(X,t),hI(X(t),t)即hI(X,t)。

说明书全文

多普勒差分测速模型及与X射线脉冲星组合导航的方法

技术领域

[0001] 本发明属于深空探测器自主导航技术领域,特别涉及一种多普勒差分测速模型及与X射线脉冲星组合导航的方法。

背景技术

[0002] 对于深空探测器而言,天文自主导航技术是至关重要的。天文自主导航系统通过测量天文数据获得导航信息。天文导航测量方法包括测、测距、测速三种方法。在深空探测转移轨道,深空探测器与近天体之间的距离很远,这将导致测角导航方法定位精度很低,无法满足深空探测任务定位高精度的要求。测距导航方法(即X射线脉冲星导航方法)的性能受X射线敏感器面积的限制,精度有限。测速导航方法则利用太阳光谱频移反演深空探测器速度,测速精度很高,可达到1cm/s量级。但该系统不具备完全可观性,无法长时间单独工作,只能作为辅助手段。
[0003] 2011年《IET Radar,Sonar&Navigation(雷达,声纳与导航)》上发表的学术论文《Doppler/XNAV-integrated navigation system using small-area X-ray sensor(采用小面积X射线敏感器的多普勒/脉冲星组合导航系统)》中,将脉冲星导航与多普勒导航相结合,提出了一种脉冲星与多普勒组合导航方法。该方法能在稳定的太阳光源下提供高精度定位和定速信息。但在实际中,受太阳耀斑、日珥、黑子等干扰因素的影响,太阳光谱并不稳定。这将会导致在某些时段导航测速信息存在较大偏差,进而影响整个组合导航系统的性能。

发明内容

[0004] 为了解决上述的技术问题,本发明提出了一种多普勒差分测速模型及将多普勒差分测速模型与X射线脉冲星导航相结合的导航方法,旨在为深空探测器提供高精度的位置和速度信息。
[0005] 本发明的多普勒差分测速模型所采用的技术方案是:一种多普勒差分测速模型,其特征在于:
[0006] 首先假设深空探测器在t1时刻和t时刻的位置分别为r1和r,火星在t1时刻的位置和速度分别为rM和vM,深空探测器上搭载两个光谱仪,第一光谱仪在t1时刻测量太阳直射光的光谱频移,第二光谱仪在t时刻测量火星反射太阳光的光谱频移,且t1和t的关系如下所示:
[0007] (式壹)
[0008] 其中,c为光速;
[0009] 基于上述假设,多普勒差分测速模型的建立包括以下步骤:
[0010] 步骤A1:在t1时刻之前,根据t和式壹确定t1的值;
[0011] 步骤A2:在t1时刻,利用第一光谱仪获取太阳径向上的多普勒速度vD1;在t时刻,利用第二光谱仪获取火星径向上的多普勒速度vD2;
[0012] 步骤A3:建立多普勒差分测速模型:
[0013] Ys(t)=hs(X,t)+ωv (式贰)
[0014] 其中,ωv为多普勒测量噪声,Ys和hs分别表示测量值和测量方程,Ys=vD2-vD1,[0015] 本发明的多普勒差分测速模型与X射线脉冲星组合导航的方法所采用的技术方案是:一种多普勒差分测速模型与脉冲星组合导航的方法,其特征在于,包括以下步骤:
[0016] 步骤B1:建立深空探测器的轨道动学模型;
[0017] 步骤B2:建立X射线脉冲星导航测距模型;
[0018] 步骤B3:建立多普勒差分测速模型;
[0019] 步骤B4:利用扩展卡尔曼滤波器滤波。
[0020] 作为优选,步骤B1中所述的建立深空探测器的轨道动力学模型,其具体实现包括以下子步骤:
[0021] 步骤1.1:确定深空探测器的状态矢量X为:
[0022] (式叁)
[0023] 其中,r=[x,y,z]T和v=[vx,vy,vz]T分别为深空探测器的位置和速度矢量,x,y,z分别为深空探测器的位置在三轴上的分量,vx,vy,vz分别为深空探测器的速度在三轴上的分量;
[0024] 步骤1.2:确定深空探测器的轨道动力学模型为:
[0025] (式肆)
[0026] 其中, 分别为x,y,z,vx,vy,vz的导数,
[0027] 式肆可表示为:
[0028] (式伍)
[0029] 其 中, 是X的 导 数, 为 时 刻t的 ,f(X,t) 为 深 空 探 测 器的状 态 转 移 模型,[x1,y1,z1]和[x2,y2,z2]分 别 是火 星 和 地 球相 对 于 太 阳系质心的相对位置矢量,μs,μm,μe分别是太阳,火星和地球的引力常数,rps,rpm,rpe分别是深空探测器到太阳质心,火星质心以及地球质心之间的距离,其 计 算 公 式 为:分别是火星
质心、地球质心分别到太阳质心之间的距离;深空探测器的导航系统噪声ω=[0,0,0,ΔFx,ΔFy,ΔFz]T,其中,ΔFx,ΔFy和ΔFz是摄动力,ω(t)为时刻t深空探测器的系统噪声。
[0030] 作为优选,步骤B2中所述的建立X射线脉冲星导航测距模型,其具体实现包括以下子步骤:
[0031] 步骤2.1:确定用于导航的X射线脉冲星数量及其方位参数;
[0032] 步骤2.2:建立X射线脉冲到达时间转换模型为:
[0033](式陆)
[0034] (式柒)
[0035] 其中, 为深空探测器与太阳系质心之间的距离在第j颗脉冲星方向上的j j投影;n 是第j颗脉冲星的方向矢量,j=1,2,…I,I为采用的X射线脉冲星数量,α 和j
δ 分别为第j颗脉冲星的赤经和赤纬, 和 分别为脉冲到达深空探测器时间与到达太阳系质心时间,c为光速, 为第j颗脉冲星到太阳系质心的距离,b为太阳系质心相对于太阳的位置矢量,|b|为位置矢量b的长度;μSun为太阳引力常数,r是深空探测器相对于太阳系质心的位置矢量,|r|为位置矢量r的长度;
X
[0036] 步骤2.3:确定X射线脉冲星导航观测模型h(X(t),t);X
[0037] 首先假设X射线脉冲星导航观测量Y 为:
[0038] (式捌)
[0039] 其对应的量测噪声为V,则X射线脉冲星导航观测模型可表示为:
[0040] YX=hX(X(t),t)+V(t) (式玖)
[0041] 其中,V(t)为时刻t的量测噪声,观测模型hX(X(t),t)如下:
[0042] (式捌)
[0043] 其中,第j颗脉冲星的相应项hj(X(t),t)的表达式如下:
[0044](式十一)。
[0045] 作为优选,步骤B4中所述的利用扩展卡尔曼滤波器滤波,其具体实现过程为:首先建立脉冲星中滤波器的测量模型,在脉冲星观测期间,其测量模型h(X,t)和测量值Y为:
[0046] h(X,t)=hs(X,t) (式十二)
[0047] Y=Ys (式十三)
[0048] 一旦获得了脉冲到达时间,测量值h(X,t)和测量模型Y则表示为:
[0049] (式十四)
[0050] (式十五)s s I I
[0051] 其中,h(X(t),t)即h(X,t),h(X(t),t)即h(X,t)。
[0052] 本发明与现有技术相比的优点在于:
[0053] (1)与现有脉冲星和多普勒组合导航方法相比,本发明在不稳定的太阳光谱下仍能正常稳定工作;
[0054] (2)与脉冲星导航方法和多普勒导航方法相比,本发明充分利用了多种导航信息,具有完全可观测性,能提供更高精度的定位信息。附图说明
[0055] 图1:本发明实施例的多普勒差分测速原理示意图。

具体实施方式

[0056] 为了便于本领域普通技术人员理解和实施本发明,下面结合附图及实施例对本发明作进一步的详细描述,应当理解,此处所描述的实施示例仅用于说明和解释本发明,并不用于限定本发明。
[0057] 本实施例首先给出火星探测器轨道,如表1所示。
[0058] 表1始轨道参数
[0059]
[0060] 请见图1,实施例提供了一种多普勒差分测速模型,首先假设深空探测器在t1时刻和t时刻的位置分别为r1和r,火星在t1时刻的位置和速度分别为rM和vM,深空探测器上搭载两个光谱仪,第一光谱仪在t1时刻测量太阳直射光的光谱频移,第二光谱仪在t时刻测量火星反射太阳光的光谱频移,且t1和t的关系如下所示:
[0061] (式壹)
[0062] 其中,c为光速;
[0063] 基于上述假设,多普勒差分测速模型的建立包括以下步骤:
[0064] 步骤A1:在t1时刻之前,根据t和式壹确定t1的值;
[0065] 步骤A2:在t1时刻,利用第一光谱仪获取太阳径向上的多普勒速度vD1;在t时刻,利用第二光谱仪获取火星径向上的多普勒速度vD2;
[0066] 步骤A3:建立多普勒差分测速模型:
[0067] Ys(t)=hs(X,t)+ωv (式贰)
[0068] 其中,ωv为多普勒测量噪声,Ys和hs分别表示测量值和测量方程,Ys=vD2-vD1,[0069] 本实施例基于以上多普勒差分测速模型,提供一种多普勒差分测速模型与X射线脉冲星组合导航的方法,包括以下步骤:
[0070] 步骤B1:建立深空探测器的轨道动力学模型;其具体实现包括以下子步骤:
[0071] 步骤1.1:确定深空探测器的状态矢量X为:
[0072] (式叁)
[0073] 其中,r=[x,y,z]T和v=[vx,vy,vz]T分别为深空探测器的位置和速度矢量,x,y,z分别为深空探测器的位置在三轴上的分量,vx,vy,vz分别为深空探测器的速度在三轴上的分量;
[0074] 步骤1.2:确定深空探测器的轨道动力学模型为:
[0075] (式肆)
[0076] 其中, 分别为x,y,z,vx,vy,vz的导数,
[0077] 式肆可表示为:
[0078] (式伍)
[0079] 其中, 是X的 导数, 为 时刻t的 ,f(X,t) 为深空探测器 的状态转移模型,[x1,y1,z1]和[x2,y2,z2]分别是火星和地球相对于太阳系质心的相对位置矢量,μs,μm,μe分别是太阳,火星和地球的引力常数,rps,rpm,rpe分别是深空探测器到太阳质心,火星质心以及地球质心之间的距离,其计算公 式 为 :分别是火星
质心、地球质心分别到太阳质心之间的距离;深空探测器的导航系统噪声ω=[0,0,0,ΔFx,ΔFx,ΔFz]T,其中,ΔFx,ΔFy和ΔFz是摄动力,ω(t)为时刻t深空探测器的系统噪声。
[0080] 步骤B2:建立X射线脉冲星导航测距模型;其具体实现包括以下子步骤:
[0081] 步骤2.1:确定用于导航的X射线脉冲星数量及其方位参数;
[0082] 本实施例的脉冲星采用了三颗脉冲星,脉冲星及其方位参数如表2所示。
[0083] 表2脉冲星方位参数
[0084]
[0085] 步骤2.2:建立X射线脉冲到达时间转换模型为:
[0086](式陆)
[0087] (式柒)
[0088] 其中, 为深空探测器与太阳系质心之间的距离在第j颗脉冲星方向上的投影;nj是第j颗脉冲星的方向矢量,j=1,2,…I,I为采用的X射线脉冲星数量,本实施例中I=3,αj和δj分别为第j颗脉冲星的赤经和赤纬, 和 分别为脉冲到达深空探测器时间与到达太阳系质心时间,c为光速, 为第j颗脉冲星到太阳系质心的距离,b为太阳系质心相对于太阳的位置矢量,|b|为位置矢量b的长度;μSun为太阳引力常数,r是深空探测器相对于太阳系质心的位置矢量,|r|为位置矢量r的长度;
[0089] 步骤2.3:确定X射线脉冲星导航观测模型hX(X(t),t);
[0090] 首先假设X射线脉冲星导航观测量YX为:
[0091] (式捌)
[0092] 其对应的量测噪声为V,则X射线脉冲星导航观测模型可表示为:
[0093] YX=hX(X(t),t)+V(t) (式玖)
[0094] 其中,V(t)为时刻t的量测噪声,观测模型hX(X(t),t)如下:
[0095] (式拾)
[0096] 其中,第j颗脉冲星的相应项hj(X(t),t)的表达式如下:
[0097](式十一)。
[0098] 步骤B3:建立多普勒差分测速模型;
[0099] 步骤B4:利用扩展卡尔曼滤波器滤波。其具体实现过程为:首先建立脉冲星中滤波器的测量模型,在脉冲星观测期间,其测量模型h(X,t)和测量值Y为:
[0100] h(X,t)=hs(X,t) (式十二)
[0101] Y=Ys (式十三)
[0102] 一旦获得了脉冲到达时间,测量模型h(X,t)和测量值Y则表示为:
[0103] (式十四)
[0104] (式十五)
[0105] 其中,hs(X(t),t)即hs(X,t),hI(X(t),t)即hI(X,t)。
[0106] 本实施例的滤波器参数如表3所示:
[0107] 表3滤波器参数
[0108]
[0109] 其中,δX(0)为深空探测器的状态初始误差,P(0)为初始状态误差矩阵,Q为状态噪声协方差, 即q1的平方, 即q2的平方。
[0110] 应当理解的是,上述针对较佳实施例的描述较为详细,并不能因此而认为是对本发明专利保护范围的限制,本领域的普通技术人员在本发明的启示下,在不脱离本发明权利要求所保护的范围情况下,还可以做出替换或变形,均落入本发明的保护范围之内,本发明的请求保护范围应以所附权利要求为准。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈