面向深空探测捕获段的深度组合导航方法技术

本发明专利技术提供一种面向深空探测捕获段的深度组合导航方法，包括预备阶段和滤波阶段；所述预备阶段建立深空探测器的轨道动力学模型、测向模型、测距模型、测速模型；所述滤波阶段利用扩展卡尔曼滤波器滤波，在导航滤波器中的状态转移模型为轨道动力学模型，导航滤波器中的测量模型选择包括在脉冲观测周期内，未获得测距信息时选择测向模型或测速模型，并对X射线敏感器接收的脉冲信号进行补偿；当脉冲信号累积完毕获得测距信息时，选择测距模型；导航滤波器利用测距模型，根据获取导航所需的位置和速度矢量。本发明专利技术抑制了脉冲到达时间多普勒偏差，滤波器收敛，定位精度高，并且对传感器要求很低。因此，本发明专利技术对航天器自主导航具有重要的实际意义。

【技术实现步骤摘要】

本专利技术属于航天器自主导航领域，特别涉及一种面向深空探测捕获段的测向/测距/测速深度组合导航方法。
技术介绍
导航信息是制导的前提，对于深空探测的成败至关重要。受超远距离及其带来的长时延影响，地面站无法提供实时高精度的导航信息，特别是在捕获段。自1990年以来，深空探测任务共失败7次，其中4次与捕获段有关。航天器自主导航系统通过测量天体信息并加以解算，可为深空探测器提供实时高精度的自主导航。因此，对于捕获段而言，天文自主导航是极其重要的。捕获段是一种高动态环境，其轨道动力学模型是强非线性，时变的。在此环境下实现高精度自主导航也是极其困难的。目前，在深空探测领域，主要有以下几种自主导航测量方法：(1)X射线脉冲星测距导航。X射线脉冲星测距导航通过观测脉冲星辐射信号，并按脉冲周期累积，可获得高信噪比累积轮廓，将其与标准轮廓对比，即可获得脉冲TOA(到达时间，time-of-arrival)，通过对TOA解算即可获得高精度测距信息。但是，在捕获段这一高动态环境下，由于航天器的变加速飞行，脉冲星辐射信号受多普勒效应影响极大，且难以补偿。脉冲轮廓发生较大畸变，脉冲到达时间存在较大偏差。这将严重影响导航性能。(2)测向导航。测向导航是传统的天文导航方式，通过测量近天体获得航天器相对于该天体的方位信息。但是，该方法无法提供高精度的航天器与近天体之间距离信息。特别是针对金星这种无天然卫星的行星，航天器仅能获得相对于一个天体的方位信息。航天器与近天体之间距离信息完全无从获得。这导航滤波器极易发散。(3)测速导航。测速导航通过测量太阳的光谱频移来获得航天器相对于恒星的速度信息。其测速精度较高。但是，测速导航方法无法直接提供位置信息，位置信息是通过积分速度信息获得。因此速度误差长期累积，位置信息必存在较大积分误差。一般来说，测速导航无法单独工作，常作为其他导航方式的辅助手段。以上三种方法各有优劣。已有学者将他们进行组合，如：测角/测距组合导航(马杰，刘劲，田金文.一种脉冲星/CNS组合导航方法，国家专利技术专利，ZL 2009100632674)，测距/测速组合导航(Liu J,Kang Z W,White P,Ma J,Tian J W.Doppler/XNAV-integrated navigation system using small-area X-ray sensor,IET Radar,Sonar and Navigation.5(9):1010-1017)。但这都
不是深度组合，没有考虑脉冲星测距导航多普勒效应的影响，更加无法适应捕获段这一高动态环境。综上所述，在深空探测捕获段这一高动态环境下，X射线脉冲星测距导航系统受多普勒效应影响严重，其测量信息存在较大偏差；测速导航无法提供高精度的定位信息；测向导航切向精度尚可，但是在径向上精度极低。即单独一种导航方式无法胜任深空探测捕获段高精度自主导航这一任务。
技术实现思路
本专利技术提出了一种面向深空探测捕获段的测向/测距/测速深度组合导航方法，旨在深空探测捕获段为航天器提供高精度的定位，定速自主导航信息。“深度”二字体现在利用测向和测速信息辅助测距信息的获取，具体为利用测向和测速信息提供高精度的三维速度估计信息，并用其补偿航天器高速飞行引起的脉冲星光子到达时间误差，从而达到抑制多普勒偏差的效果。组合导航则是利用测向，测距，测速信息更新导航滤波器状态，实现多源信息融合。本专利技术提供一种面向深空探测捕获段的深度组合导航方法，包括预备阶段和滤波阶段，所述预备阶段，包括建立导航滤波所需的各种模型，包括以下步骤，步骤A1，建立深空探测器的轨道动力学模型，实现如下，设深空探测器的状态矢量X为，其中，r＝[x,y,z]T和v＝[vx,vy,vz]T分别为深空探测器的位置和速度矢量，x,y,z分别为深空探测器的位置在三轴上的分量，vx,vy,vz分别为深空探测器的速度在三轴上的分量；则深空探测器的轨道动力学模型为，其中，分别为x,y,z，vx,vy,vz的导数，式(2)表示为，其中，是状态矢量X的导数，为时刻t的f(X,t)为深空探测器的状态转移模型，[x1,y1,z1]和[x2,y2,z2]分别是金星和地球相对于太阳系质心的相对位置矢量，μs,μv,μe分别是太阳、金星和地球的引力常数，rps,rpv,rpe分别是深空探测器到太阳质心，金星质心以及地球质心之间的距离；rsv,rse分别是金星质心、地球质心分别到太阳质心之间的距离；ω(t)为时刻t深空探测器的导航系统噪声；步骤A2，建立测向模型；步骤A3，建立测距模型；步骤A4，建立测速模型；所述滤波阶段利用扩展卡尔曼滤波器滤波，包括在导航滤波器中的状态转移模型为轨道动力学模型，导航滤波器中的测量模型选择包括以下步骤，所述滤波阶段利用扩展卡尔曼滤波器滤波，包括在导航滤波器中的状态转移模型为轨道动力学模型，导航滤波器中的测量模型选择包括以下步骤，步骤B1，在当前脉冲观测周期内，未获得当前脉冲观测周期的测距信息时，选择测向模型或测速模型，并利用基于多普勒补偿的历元叠加方法对X射线敏感器接收的脉冲信号进行补偿，补偿实现如下，步骤B11，X射线敏感器记录单个X射线光子的到达时间；步骤B12，对X射线光子到达时间进行多普勒补偿，过程如下，(a)估计航天器当前速度当滤波器有反馈时，该数值采用反馈值；否则，通过积分式三获得；(b)利用按式四补偿X射线光子到达时间；第i个子脉冲多普勒补偿量表示如下，其中，第i个脉冲周期为Pi，脉冲星观测周期内的脉冲数为N，ti为第i个子脉冲的到达时间，n为脉冲星方位矢量，T表示转置，c为光速第k个脉冲周期为Pk，是第k个脉冲周期中的航天器速度矢量；步骤B13，将光子按照预测脉冲周期进行叠加，获得脉冲TOA，得到测距信息；步骤B2，当脉冲信号累积完毕，获得测距信息时，选择测距模型；步骤B3，导航滤波器利用测距模型，根据接收到的脉冲TOA、金星方位和多普勒速度进行处理，得到状态矢量，获取导航所需的位置和速度矢量；当前脉冲观测周期结束后，返回步骤B1，继续下一脉冲观测周期的导航。而且，步骤A2中，建立测向模型如下，其中，Z是测向值，rv是金星位置矢量，υ是测向噪声；而且，步骤A3中，建立测距模型如下，其中，t和tb分别是脉冲到达航天器和太阳系质心的时间；n是脉冲星方位矢量；D0是脉冲星到太阳系质心的距离，b是太阳系质心相对于太阳质心的位置矢量，c是光速，σ是TOA测量噪声，|·|表示矢量的模。而且，步骤A4中，建立测速模型如下，其中，V是测速值，υ是测速噪声。本专利技术抑制了脉冲到达时间多普勒偏差，滤波器收敛，定位精度高，并且对传感器要求很低。因此，本专利技术对航天器自主导航具有重要的实际意义。本专利技术与现有技术相比的优点在于：(1)航天器自主导航系统对传感器要求较低。一般而言，高精度导航需要高精度的传感器。而本专利技术所需的传感器(天文光学导航相机，X射线敏感器，光谱仪)都是已有的，无需重新研制或采集，节约了成本和时间。此外，本方法即使采用低分辨率的天文光学导航相机，小面积的X射线敏感器，低精度的光谱仪也能获得高精度的导航定位信息，能满足深空探测的需求。以上各种探测器仅需一个即可，降低了载重。(2)本文档来自技高网...

【技术保护点】
一种面向深空探测捕获段的深度组合导航方法，其特征在于，包括预备阶段和滤波阶段，所述预备阶段，包括建立导航滤波所需的各种模型，包括以下步骤，步骤A1，建立深空探测器的轨道动力学模型，实现如下，设深空探测器的状态矢量X为，其中，r＝[x,y,z]T和v＝[vx,vy,vz]T分别为深空探测器的位置和速度矢量，x,y,z分别为深空探测器的位置在三轴上的分量，vx,vy,vz分别为深空探测器的速度在三轴上的分量；则深空探测器的轨道动力学模型为，其中，分别为x,y,z，vx,vy,vz的导数，式(2)表示为，其中，是状态矢量X的导数，为时刻t的f(X,t)为深空探测器的状态转移模型，[x1,y1,z1]和[x2,y2,z2]分别是金星和地球相对于太阳系质心的相对位置矢量，μs,μv,μe分别是太阳、金星和地球的引力常数，rps,rpv,rpe分别是深空探测器到太阳质心，金星质心以及地球质心之间的距离；rsv,rse分别是金星质心、地球质心分别到太阳质心之间的距离；ω(t)为时刻t深空探测器的导航系统噪声；步骤A2，建立测向模型；步骤A3，建立测距模型；步骤A4，建立测速模型；所述滤波阶段利用扩展卡尔曼滤波器滤波，包括在导航滤波器中的状态转移模型为轨道动力学模型，导航滤波器中的测量模型选择包括以下步骤，步骤B1，在当前脉冲观测周期内，未获得当前脉冲观测周期的测距信息时，选择测向模型或测速模型，并利用基于多普勒补偿的历元叠加方法对X射线敏感器接收的脉冲信号进行补偿，补偿实现如下，步骤B11，X射线敏感器记录单个X射线光子的到达时间；步骤B12，对X射线光子到达时间进行多普勒补偿，过程如下，(a)估计航天器当前速度当滤波器有反馈时，该数值采用反馈值；否则，通过积分式三获得；(b)利用按式四补偿X射线光子到达时间；第i个子脉冲多普勒补偿量表示如下，其中，第i个脉冲周期为Pi，脉冲星观测周期内的脉冲数为N，ti为第i个子脉冲的到达时间，n为脉冲星方位矢量，T表示转置，c为光速第k个脉冲周期为Pk，是第k个脉冲周期中的航天器速度矢量；步骤B13，将光子按照预测脉冲周期进行叠加，获得脉冲TOA，得到测距信息；步骤B2，当脉冲信号累积完毕，获得测距信息时，选择测距模型；步骤B3，导航滤波器利用测距模型，根据接收到的脉冲TOA、金星方位和多普勒速度进行处理，得到状态矢量，获取导航所需的位置和速度矢量；当前脉冲观测周期结束后，返回步骤B1，继续下一脉冲观测周期的导航。...

【技术特征摘要】
1.一种面向深空探测捕获段的深度组合导航方法，其特征在于，包括预备阶段和滤波阶段，所述预备阶段，包括建立导航滤波所需的各种模型，包括以下步骤，步骤A1，建立深空探测器的轨道动力学模型，实现如下，设深空探测器的状态矢量X为，其中，r＝[x,y,z]T和v＝[vx,vy,vz]T分别为深空探测器的位置和速度矢量，x,y,z分别为深空探测器的位置在三轴上的分量，vx,vy,vz分别为深空探测器的速度在三轴上的分量；则深空探测器的轨道动力学模型为，其中，分别为x,y,z，vx,vy,vz的导数，式(2)表示为，其中，是状态矢量X的导数，为时刻t的f(X,t)为深空探测器的状态转移模型，[x1,y1,z1]和[x2,y2,z2]分别是金星和地球相对于太阳系质心的相对位置矢量，μs,μv,μe分别是太阳、金星和地球的引力常数，rps,rpv,rpe分别是深空探测器到太阳质心，金星质心以及地球质心之间的距离；rsv,rse分别是金星质心、地球质心分别到太阳质心之间的距离；ω(t)为时刻t深空探测器的导航系统噪声；步骤A2，建立测向模型；步骤A3，建立测距模型；步骤A4，建立测速模型；所述滤波阶段利用扩展卡尔曼滤波器滤波，包括在导航滤波器中的状态转移模型为轨道动力学模型，导航滤波器中的测量模型选择包括以下步骤，步骤B1，在当前脉冲观测周期内，未获得当前脉冲观测周期的测距信息时，选择测向模型或测速模型，并利用基于多普勒补偿的历元叠加方法对X射线敏感器接收的脉冲信号进行补偿，补偿实现如下，步骤B11，X射线敏感器记录单个X射线光子...

【专利技术属性】
技术研发人员：刘劲，吴谨，李娟，邓慧萍，王文武，李富年，
申请(专利权)人：武汉科技大学，
类型：发明
国别省市：湖北;42