CN105300387A - 一种火星大气进入段非线性非高斯秩滤波方法 - Google Patents
一种火星大气进入段非线性非高斯秩滤波方法 Download PDFInfo
- Publication number
- CN105300387A CN105300387A CN201510738010.XA CN201510738010A CN105300387A CN 105300387 A CN105300387 A CN 105300387A CN 201510738010 A CN201510738010 A CN 201510738010A CN 105300387 A CN105300387 A CN 105300387A
- Authority
- CN
- China
- Prior art keywords
- gamma
- martian atmosphere
- detector
- cos
- martian
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/24—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
- Navigation (AREA)
Abstract
一种火星大气进入段非线性非高斯秩滤波方法,该方法有四大步骤:步骤一、建立基于火星大气进入段探测器的工程实际方程:步骤二、给定初始值:P0及系统噪声wk和vk的统计特性即均值,方差;步骤三、非线性非高斯秩方法:步骤四、令k=k+1,返回步骤三继续循环,直到k等于火星大气进入段时间截止所对应的时刻T即超音速降落伞打开所对应的时间为止,至此完成火星大气进入段非线性非高斯秩滤波方法。本发明通过应用非线性非高斯秩滤波方法,探测器在火星大气进入段,可以减少非高斯噪声对系统状态估计的影响,提高状态估计的精度,很好的满足未来火星探测任务对着陆精度的要求。
Description
技术领域
本发明涉及一种火星大气进入段非线性非高斯秩滤波方法。属于航天导航技术领域。
背景技术
火星大气进入段是整个火星进入、下降、着陆阶段历时时间最长,气动环境最复杂、最恶劣,导航误差最大,对着陆精度影响最大的阶段。未来火星探测需要完成采样返回与载人登陆等着陆误差不超过1km或小于100m的任务,这就对火星大气进入段的自主导航精度提出了更高的要求。而影响火星大气进入段自主导航精度的因素主要包括以下三个方面:
一、火星大气进入段精确的动力学模型;
二、火星大气进入段精确的传感器量测模型;
三、火星大气进入自主导航算法。
要实现精确的自主导航,传统的Kalman滤波方法对动力学模型和量测模型要求较为苛刻。然而,在火星大气进入段,气动环境复杂且恶劣,动力学模型具有强非线性,且其噪声的统计特性很难了解,所以要建立精确的动力学模型较为困难。因而,在工程实际中,常将动力学模型的噪声统计特性假设为服从高斯分布的白噪声,因此,这种假设将给实际的火星大气进入段自主导航带来潜在的导航误差。传统的Kalman滤波方法只适用于高斯噪声情况的线性滤波问题;扩展Kalman滤波方法主要适用于弱非线性系统,且其存在线性化导致的截断误差;而无迹Kalman滤波方法虽可用于强非线性系统,但其主要用于高斯噪声系统的非线性滤波。
基于以上的实际情况,建立了一种火星大气进入段非线性非高斯秩滤波方法,该方法适用于非高斯噪声系统的非线性滤波,在火星大气进入段,可以很好地减少非高斯噪声系统对自主导航精度的影响,并提高了火星大气进入段探测器的状态估计精度以满足未来火星精确着陆任务的要求。
发明内容
1、目的:本发明的目的是提供一种火星大气进入段非线性非高斯秩滤波方法,减少非高斯噪声对导航精度的影响,以减小探测器状态估计误差,提高其状态估计的精度。
2、技术方案:本发明的目的是通过以下技术方案来实现的。
一种火星大气进入段非线性非高斯秩滤波方法,该方法具体步骤如下:
步骤一、建立基于火星大气进入段探测器的工程实际方程:离散时间下的动力学系统和量测系统
xk+1=f(xk,uk)+wk(1)
zk=h(xk)+vk(2)
式中,xk表示系统状态量;zk是量测系统测量值;uk是状态方程的确定性控制项;非线性函数f(·)和h(·)分别为非线性状态转移函数,表示满足探测器动力学原理的函数和非线性量测函数,表示满足量测敏感器测量原理的函数;wk和vk分别是动力学系统噪声向量和量测噪声向量,其方差阵分别为Qk和Rk,且满足
式中,δkj是克罗内克函数,在数学中,克罗内克函数δkj是一个二元函数,克罗内克函数的自变量(输入值)一般是两个整数,如果两者相等,则其输出值为1,否则为0。
步骤二、给定初始值:P0及系统噪声wk和vk的统计特性(均值,方差)其中,为初始状态的估计值,P0为初始状态的误差方差阵。
步骤三、非线性非高斯秩方法:
(1)、秩采样点集
在4个采样点的情况下,采样策略和对称分布情况的秩采样点集{χk-1,i}为
式中,χk-1,i为xk-1的第i个采样点,共有4n个样本点;n为状态向量xk-1的维数;为Pk-1平方根的第i列向量;为标准正态偏量,用中位秩计算pj=(j+2.7)/5.4(也可用平均秩同样计算),p1=0.6852,p2=0.8704,
(2)、更新
(a)、时间更新
状态一步预测
xk/k-1,i=f(χk-1,i)i=1,2,…,4n(6)式中,xk/k-1,i为第i个采样点χk-1,i从tk-1时刻到tk时刻的一步预测。
一步预测误差的方差阵
式中,协方差权重系数
(b)、量测更新
重新秩采样
状态估计
式中
zk/k-1,i=h(χk/k-1,i)i=1,2,…,4n(11)
而Κk为tk时刻的状态滤波增益,zk为tk时刻的量测向量,为tk-1时刻到tk时刻的量测一步预测,为tk时刻的状态估计,zk/k-1,i为第i个重新采样点χk/k-1,i从tk-1时刻到tk时刻的量测一步预测。
估计误差的方差阵
式中,Pzz为系统量测输出变量的方差阵。
滤波增益矩阵
式中
而Pxz为系统状态一步预测与系统量测输出变量的协方差阵。
步骤四、令k=k+1,返回步骤三继续循环。直到k等于火星大气进入段时间截止所对应的时刻T(即超音速降落伞打开所对应的时间)为止。至此完成火星大气进入段非线性非高斯秩滤波方法。
其中,在步骤一中所述的“建立基于火星大气进入段探测器的工程实际方程”,其步骤如下:
探测器沿飞行轨迹进入火星大气,其对应的简化动力学系统为如下方程。
其中,r表示探测器到火星中心的距离,v是探测器的速度,θ是经度,λ是纬度,γ是飞行路径角,ψ是航向角,σ是滚转角(是控制量);火星引力加速度为火星引力常数μ=4.28221×1013m3/s2;气动升力加速度l和阻力加速度d分别为:其中CL和CD为升力系数和阻力系数;火星大气密度ρ近似满足指数表达形式,其表达形式为而ρ0为参考密度,hs为火星大气标高,大小为7500m,rs为距离火星表面40km的火星参考径向半径,大小为3437.2km。
目前火星进入段的量测敏感器主要有惯性量测单元中的加速度计和甚高频无线电测距,量测模型分别为
其中,am为火星固联坐标系的加速度,表示导航坐标系到火星固联坐标系的转换矩阵,表示速度坐标系到导航坐标系的转换矩阵,av为惯性量测单元中的加速度计在速度坐标系下输出的加速度,表达式分别为
av=[-d-lsinσlcosσ]T(21)
Ri(i=1,2,3)为探测器与在轨卫星或火星表面信标通过无线电通信测得二者之间的距离,其表达式为
其中,rl是探测器的位置矢量,表示第i个火星表面导航信标的位置矢量。
1)火星大气进入段对应的离散动力学系统可以改写为如下形式:
xk+1=f(xk,uk)+wk(23)式中,xk=[rvγθλψ]T为火星大气进入段动力学系统中的状态变量包含式(17)中
左边的各分量,控制量uk对应于滚转角σ。
2)对应的离散量测方程为
zk=h(xk)+vk(24)
式中
h(xk)=[am,R1,R2,R3]T(25)
其中,步骤二中根据实际情况估计初值。在火星大气进入段之前,探测器状态估计值和估计均方误差由大气层外探测器飞行段末端得到。需要特别强调的是,探测器初始状态估计值要是无偏的。
3、优点和功效:
本发明统筹考虑了火星实际大气进入过程中,非线性非高斯随机系统的探测器状态估计问题。通过应用非线性非高斯秩滤波方法,探测器在火星大气进入段,可以减少非高斯噪声对系统状态估计的影响,提高状态估计的精度。通过计算机仿真验证,该方法对系统的状态可以进行精确的估计,其中探测器的位置估计误差精度可以达到1.11m,速度估计误差精度可以达到0.01m/s。因此该方法可以很好的满足未来火星探测任务对着陆精度的要求。
附图说明
图1为火星大气进入段导航示意图。
图2a为火星大气进入段无迹Kalman滤波方法与非线性非高斯秩滤波方法高度估计误差示意图。
图2b为火星大气进入段无迹Kalman滤波方法与非线性非高斯秩滤波方法速度估计误差示意图。
图2c为火星大气进入段无迹Kalman滤波方法与非线性非高斯秩滤波方法飞行路径角估计误差示意图。
图2d为火星大气进入段无迹Kalman滤波方法与非线性非高斯秩滤波方法经度估计误差示意图。
图2e为火星大气进入段无迹Kalman滤波方法与非线性非高斯秩滤波方法纬度估计误差示意图。
图2f为火星大气进入段无迹Kalman滤波方法与非线性非高斯秩滤波方法航向角估计误差示意图。
图3为火星大气进入段非线性非高斯秩滤波方法具体实施的流程图。
图中的代号、符号说明如下:
UKF为无迹Kalman滤波方法。
RF为非线性非高斯秩滤波方法。
xk为状态向量。
zk是量测系统测量值。
uk是状态方程的确定性控制项。
f(·)和h(·)分别为非线性状态转移函数,表示满足探测器动力学原理的函数和非线性量测函数,表示满足量测敏感器测量原理的函数zk为量测向量。
wk和vk分别是动力学系统噪声向量和量测噪声向量。
为初始状态的估计值。
P0为初始状态的误差方差阵。
具体实施方式
见图1-图3,本发明涉及一种火星大气进入段非线性非高斯秩滤波方法,具体实施步骤如下:
探测器沿飞行轨迹进入火星大气,其对应的简化动力学系统为如下方程。
其中,r表示探测器到火星中心的距离,v是飞行器的速度,θ是经度,λ是纬度,γ是飞行路径角,ψ是航向角,σ是滚转角(是控制量);火星引力加速度为火星引力常数μ=4.28221×1013m3/s2;气动升力加速度l和阻力加速度d分别为:其中CL和CD为升力系数和阻力系数;火星大气密度ρ近似满足指数表达形式,其表达形式为而ρ0为参考密度,hs为火星大气标高,大小为7500m,rs为距离火星表面40km的火星参考径向半径,大小为3437.2km。
目前火星进入段的量测敏感器主要有惯性量测单元中的加速度计和甚高频无线电测距(考虑如图1所示的三颗火星表面信标),量测模型分别为
其中,am为火星固联坐标系的加速度,表示导航坐标系到火星固联坐标系的转换矩阵,表示速度坐标系到导航坐标系的转换矩阵,av为惯性量测单元中的加速度计在速度坐标系下输出的加速度,表达式分别为
av=[-d-lsinσlcosσ]T(30)
Ri(i=1,2,3)为探测器与在轨卫星或火星表面信标的无线电通信测得二者之间的距离
其中,rl是探测器的位置矢量,表示第i个火星表面导航信标的位置矢量。
步骤一:建立工程实际方程:火星大气进入段对应的离散动力学系统可以改写为如下形式:
xk+1=f(xk,uk)+wk(32)
式中,xk=[rvγθλψ]T为火星大气进入段动力学系统中的状态变量包含式(26)中左边的各分量,控制量uk对应于滚转角σ。
对应的离散量测方程为
zk=h(xk)+vk(33)
式中
h(xk)=[am,R1,R2,R3]T(34)
步骤二、给定初始值:
初始状态真实值和估计值,如表一所示
表一火星大气进入段探测器状态初始估计值及真实值
其中真实值为提前规划的火星大气进入点;估计值实际上与真实值存在一定误差。初始状态估计均方误差 而动力学系统噪声服从均值为0,方差阵为 的极值分布;量测噪声服从均值为0,方差阵为 的高斯分布。
三个火星表面导航信标的位置如表二所示。
表二火星表面导航信标的位置
步骤三、非线性非高斯秩滤波:
(1)、秩采样点集
按照公式(4)进行秩采样,得到下一步所需的秩采样点集。
(2)、更新
(a)、时间更新
按照式(5)-式(8)对系统状态进行时间更新,得到tk-1时刻到tk时刻状态一步预测及一步预测误差的方差阵。
(b)、量测更新
按照式(9)-式(16)对系统进行量测更新,得到tk时刻状态估计及估计误差的方差阵。
步骤四、令k=k+1,返回步骤三继续循环。直到k等于火星大气进入段时间截止所对应的时刻T(即超音速降落伞打开所对应的时间)为止。至此完成火星大气进入段非线性非高斯秩滤波方法。
其中截止时间主要取决于探测器与火星表面的高度及速度,能否满足超音速降落伞打开。
从图2a-f可以看出,在相同的初始条件下,非线性非高斯秩滤波方法得到的状态估计精度比无迹Kalman滤波方法的精度高,能够减少系统非高斯噪声对滤波的影响,提高状态估计的精度。由于火星大气进入段探测器与火星表面的高度及速度对于超音速降落伞打开的重要性,所以有必要对火星大气进入段探测器与火星表面的高度及速度的估计精度进行进一步的研究。从图2的子图a和子图b可以看出,非线性非高斯秩滤波方法给出的火星大气进入段探测器高度及速度的估计误差精度相对于由无迹Kalman滤波方法得到的估计误差精度提高了很多。可以进一步看出,非线性非高斯秩滤波方法相对于无迹Kalman滤波方法的优越性。
以上所述仅为本发明较佳的实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化和替换都应涵盖在本发明的保护范围之内,另外本发明提供的方法可以集成到火星大气进入探测器位置速度估计软件中。
Claims (3)
1.一种火星大气进入段非线性非高斯秩滤波方法,其特征在于:该方法具体步骤如下:
步骤一、建立基于火星大气进入段探测器的工程实际方程:离散时间下的动力学系统和量测系统
xk+1=f(xk,uk)+wk(1)
zk=h(xk)+vk(2)
式中,xk表示系统状态量;zk是量测系统测量值;uk是状态方程的确定性控制项;非线性函数f(·)和h(·)分别为非线性状态转移函数,表示满足探测器动力学原理的函数和非线性量测函数,表示满足量测敏感器测量原理的函数;wk和vk分别是动力学系统噪声向量和量测噪声向量,其方差阵分别为Qk和Rk,且满足
式中,δkj是克罗内克函数,在数学中,克罗内克函数δkj是一个二元函数,克罗内克函数的自变量即输入值是两个整数,如果两者相等,则其输出值为1,否则为0;
步骤二、给定初始值:P0及系统噪声wk和vk的统计特性即均值,方差其中,为初始状态的估计值,P0为初始状态的误差方差阵;
步骤三、非线性非高斯秩方法:
(1)、秩采样点集
在4个采样点的情况下,采样策略和对称分布情况的秩采样点集{χk-1,i}为
式中,χk-1,i为xk-1的第i个采样点,共有4n个样本点;n为状态向量xk-1的维数;为Pk-1平方根的第i列向量;为标准正态偏量,用中位秩计算pj=(j+2.7)/5.4,也能用平均秩同样计算,p1=0.6852,up1=0.4822,p2=0.8704,up2=1.1281;
(2)、更新
(a)、时间更新
状态一步预测
xk/k-1,i=f(χk-1,i)i=1,2,…,4n(6)
式中,xk/k-1,i为第i个采样点χk-1,i从tk-1时刻到tk时刻的一步预测;
一步预测误差的方差阵
式中,协方差权重系数
(b)、量测更新
重新秩采样
状态估计
式中
zk/k-1,i=h(χk/k-1,i)i=1,2,…,4n(11)
而Kk为tk时刻的状态滤波增益,zk为tk时刻的量测向量,为tk-1时刻到tk时刻的量测一步预测,为tk时刻的状态估计,zk/k-1,i为第i个重新采样点χk/k-1,i从tk-1时刻到tk时刻的量测一步预测;
估计误差的方差阵
式中,Pzz为系统量测输出变量的方差阵,
滤波增益矩阵
式中
而Pxz为系统状态一步预测与系统量测输出变量的协方差阵;
步骤四、令k=k+1,返回步骤三继续循环,直到k等于火星大气进入段时间截止所对应的时刻T即超音速降落伞打开所对应的时间为止,至此完成火星大气进入段非线性非高斯秩滤波方法。
2.根据权利要求1所述的一种火星大气进入段非线性非高斯秩滤波方法,其特征在于:在步骤一中所述的“建立基于火星大气进入段探测器的工程实际方程”,其步骤如下:
探测器沿飞行轨迹进入火星大气,其对应的简化动力学系统为如下方程:
其中,r表示探测器到火星中心的距离,v是探测器的速度,θ是经度,λ是纬度,γ是飞行路径角,ψ是航向角,σ是滚转角即控制量;火星引力加速度为火星引力常数μ=4.28221×1013m3/s2;气动升力加速度l和阻力加速度d分别为:其中CL和CD为升力系数和阻力系数;火星大气密度ρ近似满足指数表达形式,其表达形式为而ρ0为参考密度,hs为火星大气标高,大小为7500m,rs为距离火星表面40km的火星参考径向半径,大小为3437.2km;
目前火星进入段的量测敏感器有惯性量测单元中的加速度计和甚高频无线电测距,量测模型分别为
其中,am为火星固联坐标系的加速度,表示导航坐标系到火星固联坐标系的转换矩阵,表示速度坐标系到导航坐标系的转换矩阵,av为惯性量测单元中的加速度计在速度坐标系下输出的加速度,表达式分别为
av=[-d-lsinσlcosσ]T(21)
Ri(i=1,2,3)为探测器与在轨卫星或火星表面信标通过无线电通信测得二者之间的距离,其表达式为
其中,rl是探测器的位置矢量,rbi表示第i个火星表面导航信标的位置矢量;
1)火星大气进入段对应的离散动力学系统改写为如下形式:
xk+1=f(xk,uk)+wk(23)
式中,xk=[rvγθλψ]T为火星大气进入段动力学系统中的状态变量包含式(17)中左边的各分量,控制量uk对应于滚转角σ;
2)对应的离散量测方程为
zk=h(xk)+vk(24)
式中
h(xk)=[am,R1,R2,R3]T(25)。
3.根据权利要求1所述的一种火星大气进入段非线性非高斯秩滤波方法,其特征在于:步骤二中根据实际情况估计初值,在火星大气进入段之前,探测器状态估计值和估计均方误差由大气层外探测器飞行段末端得到,强调的是,探测器初始状态估计值要是无偏的。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510738010.XA CN105300387B (zh) | 2015-11-03 | 2015-11-03 | 一种火星大气进入段非线性非高斯秩滤波方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510738010.XA CN105300387B (zh) | 2015-11-03 | 2015-11-03 | 一种火星大气进入段非线性非高斯秩滤波方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105300387A true CN105300387A (zh) | 2016-02-03 |
CN105300387B CN105300387B (zh) | 2018-04-10 |
Family
ID=55197924
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510738010.XA Expired - Fee Related CN105300387B (zh) | 2015-11-03 | 2015-11-03 | 一种火星大气进入段非线性非高斯秩滤波方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105300387B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106441309A (zh) * | 2016-11-14 | 2017-02-22 | 郑州轻工业学院 | 基于协方差交叉融合的火星进入段分布式自主导航方法 |
CN106772524A (zh) * | 2016-11-25 | 2017-05-31 | 安徽科技学院 | 一种基于秩滤波的农业机器人组合导航信息融合方法 |
CN107547066A (zh) * | 2017-09-15 | 2018-01-05 | 北京航空航天大学 | 一种自校准秩滤波方法 |
CN109765141A (zh) * | 2018-12-18 | 2019-05-17 | 电子科技大学 | 一种基于swarm-c卫星提取大气密度的方法 |
CN111158382A (zh) * | 2020-01-19 | 2020-05-15 | 郑州轻工业大学 | 基于无线超宽带网络的无人车定位模型构建方法及系统 |
CN111189454A (zh) * | 2020-01-09 | 2020-05-22 | 郑州轻工业大学 | 基于秩卡尔曼滤波的无人车slam导航方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103323009A (zh) * | 2013-07-10 | 2013-09-25 | 北京航空航天大学 | 火星大气进入段的非线性三步滤波方法 |
CN103335653A (zh) * | 2013-06-06 | 2013-10-02 | 北京航空航天大学 | 火星大气进入段的自适应增量粒子滤波方法 |
CN103344244A (zh) * | 2013-07-09 | 2013-10-09 | 北京航空航天大学 | 火星大气进入段消除测量数据中系统误差的两步滤波方法 |
CN103411627A (zh) * | 2013-08-07 | 2013-11-27 | 北京航空航天大学 | 火星动力下降段非线性三步滤波方法 |
CN103411614A (zh) * | 2013-08-07 | 2013-11-27 | 北京航空航天大学 | 火星动力下降段多源信息组合导航的迭代skf方法 |
-
2015
- 2015-11-03 CN CN201510738010.XA patent/CN105300387B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103335653A (zh) * | 2013-06-06 | 2013-10-02 | 北京航空航天大学 | 火星大气进入段的自适应增量粒子滤波方法 |
CN103344244A (zh) * | 2013-07-09 | 2013-10-09 | 北京航空航天大学 | 火星大气进入段消除测量数据中系统误差的两步滤波方法 |
CN103323009A (zh) * | 2013-07-10 | 2013-09-25 | 北京航空航天大学 | 火星大气进入段的非线性三步滤波方法 |
CN103411627A (zh) * | 2013-08-07 | 2013-11-27 | 北京航空航天大学 | 火星动力下降段非线性三步滤波方法 |
CN103411614A (zh) * | 2013-08-07 | 2013-11-27 | 北京航空航天大学 | 火星动力下降段多源信息组合导航的迭代skf方法 |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106441309A (zh) * | 2016-11-14 | 2017-02-22 | 郑州轻工业学院 | 基于协方差交叉融合的火星进入段分布式自主导航方法 |
CN106772524A (zh) * | 2016-11-25 | 2017-05-31 | 安徽科技学院 | 一种基于秩滤波的农业机器人组合导航信息融合方法 |
CN106772524B (zh) * | 2016-11-25 | 2019-10-18 | 安徽科技学院 | 一种基于秩滤波的农业机器人组合导航信息融合方法 |
CN107547066A (zh) * | 2017-09-15 | 2018-01-05 | 北京航空航天大学 | 一种自校准秩滤波方法 |
CN109765141A (zh) * | 2018-12-18 | 2019-05-17 | 电子科技大学 | 一种基于swarm-c卫星提取大气密度的方法 |
CN109765141B (zh) * | 2018-12-18 | 2021-09-14 | 电子科技大学 | 一种基于swarm-c卫星提取大气密度的方法 |
CN111189454A (zh) * | 2020-01-09 | 2020-05-22 | 郑州轻工业大学 | 基于秩卡尔曼滤波的无人车slam导航方法 |
CN111158382A (zh) * | 2020-01-19 | 2020-05-15 | 郑州轻工业大学 | 基于无线超宽带网络的无人车定位模型构建方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN105300387B (zh) | 2018-04-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105300387A (zh) | 一种火星大气进入段非线性非高斯秩滤波方法 | |
Mulder et al. | Non-linear aircraft flight path reconstruction review and new advances | |
Rhudy et al. | Onboard wind velocity estimation comparison for unmanned aircraft systems | |
Rhudy et al. | UAV attitude, heading, and wind estimation using GPS/INS and an air data system | |
CN106885570A (zh) | 一种基于鲁棒sckf滤波的紧组合导航方法 | |
Rhudy et al. | Aircraft model-independent airspeed estimation without pitot tube measurements | |
CN105136145A (zh) | 一种基于卡尔曼滤波的四旋翼无人机姿态数据融合的方法 | |
CN102937449A (zh) | 惯性导航系统中跨音速段气压高度计和gps信息两步融合方法 | |
CN105865446A (zh) | 基于大气辅助的惯性高度通道阻尼卡尔曼滤波方法 | |
Dorobantu et al. | An airborne experimental test platform: From theory to flight | |
CN103708045A (zh) | 一种探月飞船跳跃式再入的在线参数辨识方法 | |
CN106441291A (zh) | 一种基于强跟踪sdre滤波的组合导航系统及导航方法 | |
RU2539140C1 (ru) | Интегрированная бесплатформенная система навигации средней точности для беспилотного летательного аппарата | |
Yan et al. | High-precision simulator for strapdown inertial navigation systems based on real dynamics from GNSS and IMU integration | |
CN103344245A (zh) | 火星进入段imu和甚高频无线电组合导航的ud-skf方法 | |
EP4033205A1 (en) | Systems and methods for model based vehicle navigation | |
CN104765373A (zh) | 一种星上相对运动状态获取方法 | |
CN106597428A (zh) | 一种海面目标航向航速估算方法 | |
CN103344244B (zh) | 火星大气进入段消除测量数据中系统误差的两步滤波方法 | |
CN103616024B (zh) | 一种行星探测进入段自主导航系统可观测度确定方法 | |
US20050159857A1 (en) | System and method for determining aircraft tapeline altitude | |
CN104297525A (zh) | 基于火箭橇试验的惯性测量系统加速度计标定方法 | |
Huang et al. | Observability analysis of flight state estimation for UAVs and experimental validation | |
CN103323009B (zh) | 火星大气进入段的非线性三步滤波方法 | |
CN103344246B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180410 Termination date: 20191103 |
|
CF01 | Termination of patent right due to non-payment of annual fee |