首页 / 专利库 / 天文学 / 深空网 / 一种动压辅助的火星大气进入段组合导航方法

一种动压辅助的火星大气进入段组合导航方法

阅读:434发布:2020-07-31

专利汇可以提供一种动压辅助的火星大气进入段组合导航方法专利检索,专利查询,专利分析的服务。并且本 发明 公开的一种动压辅助的火星大气进入段组合导航方法,涉及火星大气进入段组合导航方法,属于深空探测技术领域。本发明在航位递推导航方法的 基础 上,建立等离子模型实时估计探测器驻点的 等离子体 密度 。当密度低于通信临界 电子 密度时,通过测量探测器与无线电信标之间的相对距离信息对探测器的状态估计偏差进行修正;当密度大于通信临界电子密度时,把惯性测量单元输出的三轴 加速 度和火星进入大气数据系统测量得到的动压值视为外部测量,并结合非线性滤波方法对探测器的状态进行估计,提高大气进入段导航方法可靠性及状态估计 精度 。本发明可以有效抑制无线电测量长时间中断时采用单一航位递推导航方法探测器的状态估计偏差存在发散的难题。,下面是一种动压辅助的火星大气进入段组合导航方法专利的具体信息内容。

1.一种动压辅助的火星大气进入段组合导航方法,其特征在于:在航位递推导航方法的基础上,通过建立的等离子模型实时估计探测器驻点的等离子体密度,根据在线实时估计得到的等离子体密度自主切换导航方法;当等离子体密度低于通信临界电子密度时,无线电通信正常,采用无线电和惯性测量输出的组合导航方法,通过测量探测器与无线电信标之间的相对距离信息对探测器的状态估计偏差进行修正;当等离子密度大于无线电通信临界电子密度时,无线电通信“黑障”,采用动压辅助的导航方法,把惯性测量单元输出的三轴加速度和火星进入大气数据系统测量得到的动压值视为外部测量,并结合非线性滤波方法对探测器的状态进行估计,进而提高星大气进入段导航方法的可靠性及状态估计精度
2.如权利要求1所述的一种动压辅助的火星大气进入段组合导航方法,其特征在于:
所述的非线性滤波方法采用扩展卡尔曼滤波(Extend Kalman Filter,EKF),无迹卡尔曼滤波(Unscented Kalman Filter,UKF)算法
3.如权利要求1所述的一种动压辅助的火星大气进入段组合导航方法,其特征在于:
具体实现方法包括如下步骤,
步骤1:建立等离子密度模型,实时估计探测器驻点的等离子体密度;
在火星大气进入高超声速段,探测器与火星大气剧烈摩擦产生的气动热会导致探测器周围的空气电离产生等离子体,当等离子体的密度超过探测器的无线电通信频率对应的临界电子密度时,会导致无线电通信中断,产生通信“黑障”;因此,实时估计出探测器周围的等离子体密度是动压辅助导航方法自主切换的前提条件;探测器周围的等离子体密度主要由探测器相对火星表面的速度以及火星大气密度决定,通过式(1)获得探测器驻点等离子密度的估计值:
a b
ne=cρ V (1)
其中,a,b,c为模型系数,V为探测器相对于火星表面的速度,单位为km/s,ρ为火星大气密度,假设火星大气密度ρ呈指数形式分布,由式(2)计算得到:
3
其中,ρ0=2e-4kg/m 表示火星大气参考密度,r0=3,437,200m表示参考高度,hs=
7500m火星大气标高;
步骤2:根据步骤1所述的等离子密度模型能够实时估计出探测器周围的等离子体密度;当等离子密度低于探测器无线电通信频率对应的临界电子密度时,无线电通信正常,系统采用无线电和惯性测量输出的组合导航方法,当等离子体密度大于无线电通信临界电子密度时,无线电通信“黑障”,探测器采用动压辅助的导航方法;所述的动压测量辅助的导航方法指把火星进入大气数据系统测量得到的驻点动压值和惯性测量单元输出的三轴加速度值视为外部测量,通过非线性滤波算法对探测器的状态进行估计;根据无线电导航方法或动压辅助组合导航方法分别建立对应的无线电测量模型和动压测量模型;
步骤3:根据步骤1得到的等离子密度模型,以及步骤2建立的无线电测量模型和动压测量模型,通过非线性滤波算法对探测器的状态进行估计,进而提高火星大气进入段导航状态估计精度。
4.如权利要求3所述的一种动压辅助的火星大气进入段组合导航方法,其特征在于:
所述的步骤2具体实现方法为,
当等离子密度低于通信临界电子密度时,无线电通信正常,通过测量探测器与无线电信标之间的相对距离对探测器的状态偏差进行修正;
所述的无线电信标的数量及位置根据导航精度需要而定,无线电信标越多,测量信息越丰富,导航精度越高;假设探测器在进入过程中有三颗信标在可见范围内:一颗轨道器和两颗火星表面信标;无线电信标的位置在探测器进入火星大气前由地面深空网准确测定;探测器与第i颗信标之间的距离由式(3)计算得到:
其中,r,ri分别表示探测器和第i(i=1,2,3)颗信标在火星惯性系下的位置矢量,εRi表示无线电测量噪声,服从高斯分布;
当等离子体密度大于无线电通信临界电子密度时,即通信“黑障”时,探测器通过惯性测量输出的三轴加速度值及火星进入大气数据系统获得飞行过程中探测器驻点动压值,通过非线性估计算法估计探测器的状态;所述的动压测量辅助的导航方法的测量模型建立如式(4):
其中,εDA表示测量噪声,服从高斯分布;a表示探测器三轴加速度值,由式(5)确定:
q∞表示探测器驻点动压值,由式(6)确定:
其中, 表示速度坐标系单位矢量在火星惯性系下的表示;σ表示倾侧
,为系统控制变量,假设倾侧角σ为常值;D、L表示阻加速和升力加速度,由式(7)确定:
(7)
其中,CD,CL表示探测器的阻力系数和升力系数,S表示探测器参考面积,m表示探测器质量,参数CD,CL,S,m都视为常值;
速度坐标系三轴定义由式(8)和式(9)确定:
T
其中, 为探测器相对火星表面的速度,Ω=[0, 0, 7.095e-5]rad/s
表示火星自转角速度在火星惯性系下的表示;
则整个导航方法的量测模型如下:
等离子密度低于无线电通信临界电子密度时,即通信正常时:
等离子体密度大于无线电通信临界电子密度时,即通信“黑障”时:
火星大气进入段导航系统根据系统在线实时估计得到的等离子体密度自主切换导航方法,确保进入过程中火星大气进入段导航系统能够获得的测量信息最大化。
5.如权利要求3或4所述的一种动压辅助的火星大气进入段组合导航方法,其特征在于:所述的非线性滤波方法采用扩展卡尔曼滤波(Extend Kalman Filter,EKF),无迹卡尔曼滤波(Unscented Kalman Filter,UKF)算法。

说明书全文

一种动压辅助的火星大气进入段组合导航方法

技术领域

[0001] 本发明涉及一种火星大气进入段组合导航方法,尤其涉及一种动压辅助的火星大气进入段组合导航方法,属于深空探测技术领域。

背景技术

[0002] 未来火星采样返回及载人火星探测任务要求探测器具有定点着陆的能,定点着陆要求在大气进入段对探测器进行精确的导航、制导与控制,确保探测器能精确达到预开伞点。在火星大气进入段精确确定探测器的状态信息是定点着陆探测的前提条件。
[0003] 以往成功着陆的火星探测任务中,火星大气进入段都采用基于惯性测量单元(Inertial Measurement Unit,IMU)输出的航位递推(Dead Reckoning,DR)导航方法。由于航位递推导航方法不能对初始状态偏差进行修正,以及难以消除敏感器系统偏差对状态估计精度的影响,导致状态估计偏差随时间发散,难以满足未来定点着陆的导航精度需求。而目前基于无线电和惯性测量单元输出的组合导航方法中,由于无线电信号对外界环境敏感,尤其在高超声速段,探测器与火星大气剧烈摩擦产生的气动热会导致周围空气电离。等离子体密度超过无线电通信频率对应的临界电子密度时会导致无线电通信中断,产生通信“黑障”。长时间的通信“黑障”会导致进入段状态估计偏差随时间发散。
[0004] 通过在线实时估计探测器周围的等离子体密度,把火星大气进入段导航分为通信正常和通信“黑障”两部分。等离子体密度低于无线电通信临界电子密度时,无线电通信正常,通过测量探测器与无线电信标之间的相对距离,可对探测器的状态进行估计;当等离子密度高于无线电通信临界电子密度时,无线电通信“黑障”,采用动压辅助的(Dynamic Pressure Measurement Aid,DA)导航方法,把驻点动压测量值及惯性测量单元三轴加速度信息视为外部测量,对探测器的状态进行估计。

发明内容

[0005] 针对目前无线电和惯性测量单元输出的组合导航方法可能存在无线电通信“黑障”,及“黑障”段采用单一航位递推导航方法探测器的状态估计偏差存在发散难题。本发明公开的一种动压辅助的火星大气进入段组合导航方法要解决的技术问题是,解决现有组合导航方法中存在的上述技术难题,进而提高火星大气进入段导航方法的可靠性及状态估计精度。
[0006] 本发明的目的是通过下述技术方案实现的:
[0007] 本发明公开的一种动压辅助的火星大气进入段组合导航方法,在航位递推导航方法的基础上,通过建立的等离子模型实时估计探测器驻点的等离子体密度,根据在线实时估计得到的等离子体密度自主切换导航方法。当等离子体密度低于通信临界电子密度时,即无线电通信正常时,采用无线电和惯性测量单元输出的组合导航方法,通过测量探测器与无线电信标之间的相对距离信息对探测器的状态估计偏差进行修正;当等离子密度大于无线电通信临界电子密度时,即通信“黑障”时,采用动压辅助的导航方法,把惯性测量单元输出的三轴加速度和火星进入大气数据系统测量得到的动压值视为外部测量,并结合非线性滤波方法对探测器的状态进行估计,进而提高星大气进入段导航方法的可靠性及状态估计精度。
[0008] 本发明公开的一种动压辅助的火星大气进入段组合导航方法,可以有效抑制无线电测量长时间中断时采用单一航位递推导航方法探测器的状态估计偏差存在发散的难题,保证了导航方法的可靠性及状态估计精度。
[0009] 所述的非线性滤波方法采用扩展卡尔曼滤波(Extend Kalman Filter,EKF),无迹卡尔曼滤波(Unscented Kalman Filter,UKF)算法
[0010] 本发明公开的一种动压辅助的火星大气进入段组合导航方法,包括如下步骤:
[0011] 步骤1:建立等离子密度模型,实时估计探测器驻点的等离子体密度。
[0012] 在火星大气进入高超声速段,探测器与火星大气剧烈摩擦产生的气动热会导致探测器周围的空气电离产生等离子体,当等离子体的密度超过探测器的无线电通信频率对应的临界电子密度时,会导致无线电通信中断,产生通信“黑障”。因此,实时估计出探测器周围的等离子体密度是动压辅助导航方法自主切换的前提条件。探测器周围的等离子体密度主要由探测器相对火星表面的速度以及火星大气密度决定,通过式(1)获得探测器驻点等离子密度的估计值:
[0013] ne=cρaVb (1)
[0014] 其中,a,b,c为模型系数,V为探测器相对于火星表面的速度,单位为km/s,ρ为火星大气密度,假设火星大气密度ρ呈指数形式分布,由式(2)计算得到:
[0015]
[0016] 其中,ρ0=2e-4kg/m3表示火星大气参考密度,r0=3,437,200m表示参考高度,hs=7500m火星大气标高。
[0017] 步骤2:根据步骤1所述的等离子密度模型能够实时估计出探测器周围的等离子体密度。当等离子密度低于探测器无线电通信频率对应的临界电子密度时,无线电通信正常,系统采用无线电和惯性测量输出的组合导航方法,当等离子体密度大于无线电通信临界电子密度时,无线电通信“黑障”,探测器采用动压辅助的导航方法。所述的动压测量辅助的导航方法指把火星进入大气数据系统测量得到的驻点动压值和惯性测量单元输出的三轴加速度值视为外部测量,通过非线性滤波算法对探测器的状态进行估计。根据无线电导航方法或动压辅助组合导航方法分别建立对应的无线电测量模型和动压测量模型。
[0018] 当等离子密度低于通信临界电子密度时,无线电通信正常,通过测量探测器与无线电信标之间的相对距离对探测器的状态偏差进行修正。
[0019] 所述的无线电信标的数量及位置根据导航精度需要而定,无线电信标越多,测量信息越丰富,导航精度越高。假设探测器在进入过程中有三颗信标在可见范围内:一颗轨道器和两颗火星表面信标。无线电信标的位置在探测器进入火星大气前由地面深空网准确测定。探测器与第i颗信标之间的距离由式(3)计算得到:
[0020]
[0021] 其中,r,ri分别表示探测器和第i(i=1,2,3)颗信标在火星惯性系下的位置矢量,εRi表示无线电测量噪声,服从高斯分布。
[0022] 当等离子体密度大于无线电通信临界电子密度时,无线电通信“黑障”,探测器通过惯性测量输出的三轴加速度值及火星进入大气数据系统获得飞行过程中探测器驻点动压值,通过非线性估计算法估计探测器的状态。所述的动压测量辅助的导航方法的测量模型建立如式(4):
[0023]
[0024] 其中,εDA表示测量噪声,服从高斯分布;a表示探测器三轴加速度值,由式(5)确定:
[0025]
[0026] q∞表示探测器驻点动压值,由式(6)确定:
[0027]
[0028] 其中, 表示速度坐标系单位矢量在火星惯性系下的表示;σ表示倾侧,为系统控制变量,假设倾侧角σ为常值;D、L表示阻力加速和升力加速度,由式(7)确定:
[0029]
[0030] 其中,CD,CL表示探测器的阻力系数和升力系数,S表示探测器参考面积,m表示探测器质量,参数CD,CL,S,m都视为常值。
[0031] 速度坐标系三轴定义由式(8)和式(9)确定:
[0032]
[0033]
[0034] 其中, 为探测器相对火星表面的速度,Ω=[0,0,7.095e-5]Trad/s表示火星自转角速度在火星惯性系下的表示。
[0035] 则整个导航方法的量测模型如下:
[0036] 等离子密度低于无线电通信临界电子密度时,即通信正常时:
[0037]
[0038] 等离子体密度大于无线电通信临界电子密度时,即通信“黑障”时:
[0039]
[0040] 火星大气进入段导航系统根据系统在线实时估计得到的等离子体密度自主切换导航方法,确保进入过程中火星大气进入段导航系统能够获得的测量信息最大化。
[0041] 步骤3:根据步骤1得到的等离子密度模型,以及步骤2建立的无线电测量模型和动压测量模型,通过非线性滤波算法对探测器的状态进行估计,进而提高火星大气进入段导航系统的可靠性及状态估计精度。
[0042] 有益效果:
[0043] 1、本发明公开的一种动压辅助的火星大气进入段组合导航方法,可以实时估计探测器驻点等离子密度,通过对比驻点等离子体密度与无线电通信临界电子密度的大小自主切换无线电导航和动压辅助导航方法。
[0044] 2、本发明公开的一种动压辅助的火星大气进入段组合导航方法,能够增加当等离子体密度大于无线电通信临界电子密度时火星大气进入段测量信息,抑制单一航位递推导航方法探测器的状态估计偏差存在发散,提高火星大气进入段导航方法的可靠性及状态估计精度。
[0045] 3、本发明公开的一种动压辅助的火星大气进入段组合导航方法,采用非线性滤波算法,可以提高火星大气进入段导航数据解算速度,满足未来火星大气进入段导航实时性需求。附图说明
[0046] 图1为动压测量辅助导航示意图;
[0047] 图2为火星大气进入段导航解算流程图
[0048] 图3 I-VI分别表示动压测量辅助的组合导航与传统无线电和惯性测量单元输出的组合导航位置、速度、航迹角、方位角、经度及纬度估计偏差对比图。图中符号a所指的实线代表本发明组合导航方法获得的状态估计偏差;符号b所指的点线代表传统无线电和惯性测量单元输出的组合导航的状态偏差;符号c所指的虚线表示本发明状态估计3sigma曲线,符号d所指的点画线表示传统无线电和惯性测量单元输出的组合导航状态估计3sigma曲线。

具体实施方式

[0049] 为了更好的说明本发明的目的和优点,下面结合附图和实例对本发明的内容做进一步说明。
[0050] 本实例针对火星大气进入段动压测量辅助的组合导航方法,当等离子体密度低于无线电通信临界电子密度式时,无线电通信正常,探测器采用无线电和惯性测量单元输出的组合导航方法;当等离子密度大于无线电通信临界电子密度时,无线电通信“黑障”,采用动压辅助的组合导航方法。并采用EKF滤波器对探测器的状态进行解算,增加了无线电测量长时间中断时导航系统的可靠性及状态估计精度。本实例的具体实施方法如下:
[0051] 步骤1:建立等离子密度模型,实时估计探测器驻点的等离子体密度。
[0052] 在火星大气进入高超声速段,探测器与火星大气剧烈摩擦产生的气动热会导致探测器周围的空气电离产生等离子体,当等离子体的密度超过探测器的无线电通信频率对应的临界电子密度时,会导致无线电通信中断,产生通信“黑障”。因此,实时估计出探测器周围的等离子体密度是动压辅助导航方法自主切换的前提条件。探测器周围的等离子体密度主要由探测器相对火星表面的速度以及火星大气密度决定,通过式(12)获得探测器驻点等离子密度的估计值:
[0053] ne=cρaVb (12)
[0054] 其中,a,b,c为模型系数,V为探测器相对于火星表面的速度,单位为km/s,ρ为火星大气密度,假设火星大气密度ρ呈指数形式分布,由式(13)计算得到:
[0055]
[0056] 其中,ρ0=2e-4kg/m3表示火星大气参考密度,r0=3,437,200m表示参考高度,hs=7500m火星大气标高。
[0057] 步骤2:根据步骤1所述的等离子密度模型能够实时估计出探测器周围的等离子体密度。当等离子密度低于探测器无线电通信频率对应的临界电子密度时,即通信正常时,系统采用无线电和惯性测量输出的组合导航方法,当等离子体密度大于无线电通信临界电子密度时,即无线电通信“黑障”时,探测器采用动压辅助的导航方法。所述的动压测量辅助的导航方法指把火星进入大气数据系统测量得到的驻点动压值和惯性测量单元输出的三轴加速度值视为外部测量,通过非线性滤波算法对探测器的状态进行估计。根据无线电导航方法或动压辅助组合导航方法分别建立对应的无线电测量模型和动压测量模型。
[0058] 当等离子密度低于通信临界电子密度时,无线电通信正常,通过测量探测器与无线电信标之间的相对距离对探测器的状态偏差进行修正。
[0059] 所述的无线电信标的数量及位置根据导航精度需要而定,无线电信标越多,测量信息越丰富,导航精度越高。假设探测器在进入过程中有三颗信标在可见范围内:一颗轨道器和两颗火星表面信标。无线电信标的位置在探测器进入火星大气前由地面深空网准确测定。探测器与第i颗信标之间的距离由式(14)计算得到:
[0060]
[0061] 其中,r,ri分别表示探测器和第i(i=1,2,3)颗信标在火星惯性系下的位置矢量,εRi表示无线电测量噪声,服从高斯分布。
[0062] 当等离子体密度大于无线电通信临界电子密度时,即通信“黑障”时,探测器通过惯性测量输出的三轴加速度值及火星进入大气数据系统获得飞行过程中探测器驻点动压值,通过非线性估计算法估计探测器的状态。所述的动压测量辅助的导航方法的测量模型建立如式(15):
[0063]
[0064] 其中,εDA表示测量噪声,服从高斯分布;a表示探测器三轴加速度值,由式(16)确定:
[0065]
[0066] q∞表示探测器驻点动压值,由式(17)确定:
[0067]
[0068] 其中, 表示速度坐标系单位矢量在火星惯性系下的表示;σ表示倾侧角,为系统控制变量,假设倾侧角σ为常值;D、L表示阻力加速和升力加速度,可由式(18)确定:
[0069]
[0070] 其中,CD,CL表示探测器的阻力系数和升力系数,S表示探测器参考面积,m表示探测器质量,参数CD,CL,S,m都视为常值。
[0071] 速度坐标系三轴定义由式(19)和式(20)确定:
[0072]
[0073]
[0074] 其中, 为探测器相对火星表面的速度,Ω=[0,0,7.095e-5]Trad/s表示火星自转角速度在火星惯性系下的表示。
[0075] 则整个导航方法的量测模型如下:
[0076] 等离子密度低于无线电通信临界电子密度时,即通信正常时:
[0077]
[0078] 等离子体密度大于无线电通信临界电子密度时,即通信“黑障”时:
[0079]
[0080] 火星大气进入段导航系统根据系统在线实时估计得到的等离子体密度自主切换导航方法,确保进入过程中火星大气进入段导航系统能够获得的测量信息最大化。
[0081] 步骤3:根据步骤1得到的等离子密度模型,以及步骤2建立的无线电测量模型和动压测量模型,通过非线性滤波算法对探测器的状态进行估计,进而提高火星大气进入段导航状态估计精度。
[0082] 仿真中探测器初始状态及偏差如表1所示,其它物理参数如表2所示。
[0083] 表1、探测器初始状态及偏差
[0084]
[0085]
[0086] 表2、探测器物理参数
[0087]
[0088] 加速度敏感器性能指标如表3所示:
[0089] 表3、LN-200惯性单元性能指标
[0090]
[0091] 表4、“黑障”结束点两种导航方法状态估计偏差对比
[0092]
[0093] 表4给出了本发明与传统无线电和惯性测量输出的组合导航方法开伞点精度对比,当“黑障”时间为40s时,本发明相比传统无线电和惯性测量输出的组合导航方法在“黑障”结束点能提高一个数量级的估计精度。从图3可以看出,采用本发明的动压测量辅助的组合导航方法在无线电长时间中断时状态估计偏差收敛性更好,“黑障”段状态估计精度更高,可靠性更好。
[0094] 本发明保护范围不仅局限于实施例,实施例用于解释本发明,凡与本发明在相同原理和构思条件下的变更或修改均在本发明公开的保护范围之内。
高效检索全球专利

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

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

申请试用

分析报告

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

申请试用

QQ群二维码
意见反馈