CN103389095A - 一种用于捷联惯性/多普勒组合导航***的自适应滤波方法 - Google Patents
一种用于捷联惯性/多普勒组合导航***的自适应滤波方法 Download PDFInfo
- Publication number
- CN103389095A CN103389095A CN2013103136494A CN201310313649A CN103389095A CN 103389095 A CN103389095 A CN 103389095A CN 2013103136494 A CN2013103136494 A CN 2013103136494A CN 201310313649 A CN201310313649 A CN 201310313649A CN 103389095 A CN103389095 A CN 103389095A
- Authority
- CN
- China
- Prior art keywords
- constantly
- expression
- sins
- phi
- navigation
- 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.)
- Pending
Links
Images
Landscapes
- Navigation (AREA)
Abstract
本发明公开了一种用于捷联惯性/多普勒组合导航***的自适应滤波方法,其目的是提高滤波器在动态条件下的响应速度和滤波精度,并提高捷联惯性/多普勒组合导航***的定位精度。该方法通过引入关于新息协方差的限定窗口平滑器,利用该新息协方差平滑值直接修正滤波器中的增益矩阵,并引入调节因子对一步预测均方误差进行修正,进而达到改善滤波器动态响应速度以及提高滤波精度的目的。本发明所设计的自适应滤波方法用于捷联惯性/多普勒组合导航***中,有效地提高了组合***在动态条件下的导航定位精度。
Description
技术领域
本发明属于组合导航***技术领域,尤其涉及的是一种捷联惯性/多普勒组合导航***的自适应滤波方法。
背景技术
常规卡尔曼滤波可以估计组合导航***的位置、速度误差以及平台失准角等状态变量,因此在组合导航领域中获得了广泛的应用。但是常规卡尔曼滤波在工作时需要满足模型精确和噪声统计特性准确已知的要求,对于捷联惯性与多普勒组合导航***而言,***中不可避免的含有不确定性噪声,例如陀螺随机漂移、随机风速以及海流等干扰因素,此时常规卡尔曼滤波将失去最优估计特性;而且当常规卡尔曼滤波在特定条件下达到稳态时,滤波参数也会收敛于稳态值,若此时***的外部运动发生剧烈变化,滤波参数不能立刻做出调整,使得滤波精度下降,甚至引起滤波发散。
为解决常规卡尔曼滤波在噪声统计信息不准确以及动态条件下导致滤波精度下降等问题,近年来学者们提出了许多自适应滤波算法,主要集中在新息自适应滤波和渐消记忆衰减自适应滤波两方面。文献“基于新息协方差的自适应渐消卡尔曼滤波器,***工程与电子技术,2012,第13卷第12期”,提出了一种利用新息协方差计算渐消因子的方法,该方法利用新息序列的协方差自适应的改变渐消因子调整新息的权重,减小了陈旧量测值对估计的影响,而且通过渐消因子自适应地调整误差协方差,补偿不完整信息的影响;但是该方法需要同时求取两个调节因子,不仅增加了计算量,还由于调节因子本身会带有一定的误差,可能会引起估计精度下降,甚至导致滤波器的发散。
发明内容
为了克服现有技术中的缺陷,解决上述技术问题,本发明提供一种用于捷联惯性/多普勒组合导航***的自适应滤波方法,该方法可以明显改善由于噪声统计信息不准导致常规卡尔曼滤波发散的情况,提高滤波器在动态条件下的响应速度和滤波精度,可有效改进捷联惯性/多普勒组合导航***的定位精度,其技术方案如下:
一种用于捷联惯性/多普勒组合导航***的自适应滤波方法,包括以下步骤:
步骤一:利用捷联惯性导航***中惯性测量单元测得沿载体系相对惯性空间的角速率和加速度分量信息,并与***的初始位置、速度以及姿态信息一起进行捷联惯导解算,得到载体实时的位置、速度和姿态信息;
步骤二:根据捷联惯导***的误差方程选取状态变量,建立***的状态方程,并选取经纬度误差、速度误差以及失准角作为状态变量;利用多普勒导航***与捷联惯导***所提供速度的差值作为量测变量,建立组合导航***的量测方程;
步骤三:经过连续***离散化处理后,对组合导航***进行滤波估计;首先更新状态一步预测值及其均方误差Pk,k-1,并利用k时刻量测信息Zk和滤波器的一步预测值计算k时刻的新息序列εk;再确定平滑窗口的宽度N,并借助量测新息序列εk建立关于新息协方差的N步限定窗口平滑器
所述的一步预测值:
所述的一步预测均方误差:
所述的k时刻新息序列:
所述的N步新息协方差平滑器:
其中,Pk-1表示k-1时刻的滤波均方误差,Qk-1表示k-1时刻***噪声协方差,表示k-1时刻状态转移矩阵的转置,表示k-1时刻***噪声系数矩阵的转置,εk表示k时刻的新息序列,和分别是k和k-1时刻的新息协方差N步平滑值,N=10,且当k<N时,
步骤六:将步骤三至步骤五中得到的滤波参数带入到自适应滤波器中,并不断重复以上过程完成对导航参数误差的估计;利用输出校正实时修正捷联惯导***输出的导航信息,从而得到更高精度的导航参数信息,即补偿后的位置、速度和姿态信息,直至捷联惯性/多普勒组合导航过程结束。
所述的方法,所述步骤二中,
所述的状态方程:Xk=ΦkXk-1+ΓkWk
所述的量测方程:Zk=HkXk+Vk
状态一步转移矩阵
Φk=
***噪声系数阵
Γk=
量测系数阵
其中,wie表示地球自转角速度,L表示载体所在纬度,g表示地球重力加速度,R表示地球半径,T表示滤波器数据处理周期,vE、vN分别表示载体相对地球的东北向速度,fE、fN表示加速度计测得的比力信息,表示载体坐标系到导航坐标系之间的转换矩阵。
所述的方法,所述步骤四的具体方法为:
(2)修正k时刻的滤波增益矩阵Kk,其计算公式写作:
其中,Pk,k-1表示k时刻一步预测均方误差,Hk表示k时刻量测系数矩阵。
所述的方法,所述步骤五的具体方法为:
(1)得到k时刻的调节因子ρk,其表达式写作:
其中,max{□}表示取内部所有元素的最大值,tr[□]表示对内部矩阵求迹;
(2)修正一步预测均方误差Pk,k-1,其表达式写作:
其中,Φk-1表示k-1时刻的***状态转移矩阵,Pk-1表示k-1时刻的滤波估计均方误差,Γk-1表示k-1时刻的***噪声系数阵,Qk-1表示k-1时刻的***噪声协方差。
所述的方法,所述步骤六的具体方法为:
(1)自适应滤波器对导航参数误差的估计值:
k时刻状态估计值:
k时刻估计均方误差:Pk=(I-KkHk)Pk,k-1
(2)得到导航参数误差的估计值后,利用输出校正对捷联惯导***输出的导航参数进行补偿:
①位置修正: λ=λSINS-δλ
L=LSINS-δL
其中,λ、L分别表示修正后的经纬度信息,λSINS、LSINS分别表示捷联惯导***输出的经纬度信息,δλ、δL分别表示滤波器估计的经纬度误差值;
②速度修正:vE=vSINS-E-δvE
vN=vSINS-N-δvN
其中,vE、vN分别表示修正后的东北向速度信息,vSINS-E、vSINS-N分别表示捷联惯导***输出的东北向速度信息,δvE、δvN分别表示滤波器估计的东北向速度误差值;
③姿态修正:
其中,φx、φy和φz分别表示计算导航系与真实导航系坐标轴之间的失准角;
且修正后的姿态角矩阵可以表示成矩阵的形式:
c)计算补偿后的姿态角,再根据
本发明的有益效果:
1、本发明通过对新近的新息协方差序列进行窗口平滑,使得新息变化更加平稳,并直接对滤波器的增益矩阵进行修正,提高了滤波器的鲁棒性。
2、本发明通过引入调节因子对滤波过程进行修正,使得滤波器能够根据量测信息的变化及时修正滤波估计值,提高了滤波器在动态条件下的响应速度以及滤波精度。
3、本发明的有益效果可以通过Matlab仿真实验进行验证,在捷联惯性/多普勒组合导航***中仿真参数设置如下:
①自适应滤波器的初始参数设置如下:
P(0)=diag{(100m/Re)2,(100m/Re)2,(0.5m/s)2,(0.5m/s)2,(0.1°)2,(0.1°)2,(0.1°)2}
Q(0)=diag{(100ug)2,(100ug)2,(0.01°/h)2,(0.01°/h)2,(0.01°/h)2}
R(0)=10R=10×diag{(0.1m/s)2,(0.1m/s)2}
②初始位置姿态:经度为126.67°,纬度为45.77°,初始速度(6m/s,8m/s,0),
姿态角(0,0,30°);仿真时间3600s,计算步长设为0.1s。
由仿真结果图3-图5可以得出,本发明提出的自适应滤波方法通过对新息协方差序列进行限定窗口平滑处理,并直接调节增益矩阵,而且通过引入调节因子来修正一步预测均方误差,不仅改善了滤波器的鲁棒性,而且进一步提高了滤波精度;本发明提出的自适应滤波方法用于捷联惯性/多普勒组合导航***中,提高了组合导航***的定位精度,具有广阔的工程应用前景。
本发明提出一种自适应卡尔曼滤波的设计方法,通过引入关于新息协方差的限定窗口平滑器,对滤波器的增益矩阵进行实时修正,通过引入调节因子对一步预测均方误差进行修正,使得滤波器能够根据新近的新息序列自适应的调节增益矩阵,提高滤波器在动态条件下的响应速度和滤波精度,进而提高捷联惯性/多普勒组合导航***的定位精度,具有实际的工程应用意义。
附图说明
图1是本发明自适应滤波方法在组合导航***中的工作流程图;
图2是本发明自适应滤波方法用于捷联惯性/多普勒组合***的方框图;
图3是本发明自适应滤波方法在组合导航***中的速度误差仿真图;
图4是本发明自适应滤波方法在组合导航***中的位置误差仿真图;
图5是本发明自适应滤波方法在组合导航***中的姿态误差角仿真图;
具体实施方式
以下结合具体实施例,对本发明进行详细说明。
一种用于捷联惯性/多普勒组合导航***的自适应滤波方法,如图1-图2所示,通过下列步骤实现自适应滤波过程:
步骤一:利用捷联惯性导航***中惯性测量单元测得沿载体系相对惯性空间的角速率和加速度分量信息,并与***的初始位置、速度以及姿态信息一起进行捷联惯导解算,得到载体实时的位置、速度和姿态信息;
步骤二:根据捷联惯导***的误差方程选取状态变量,建立***的状态方程,并选取经纬度误差、速度误差以及失准角作为状态变量;利用多普勒导航***与捷联惯导***所提供速度的差值作为量测变量,建立组合导航***的量测方程;
状态方程:Xk=ΦkXk-1+ΓkWk
量测方程:Zk=HkXk+Vk
状态一步转移矩阵
Φk=
***噪声系数阵
Γk=
量测系数阵
其中,wie表示地球自转角速度,L表示载体所在纬度,g表示地球重力加速度,R表示地球半径,T表示滤波器数据处理周期,vE、vN分别表示载体相对地球的东北向速度,fE、fN表示加速度计测得的比力信息,表示载体坐标系到导航坐标系之间的转换矩阵。
步骤三:经过连续***离散化处理后,对组合导航***进行滤波估计。首先更新状态一步预测值及其均方误差Pk,k-1,并利用k时刻量测信息Zk和滤波器的一步预测值计算k时刻的新息序列εk;再确定平滑窗口的宽度N,并借助量测新息序列εk建立关于新息协方差的N步限定窗口平滑器
一步预测值:
一步预测均方误差:
k时刻新息序列:
N步新息协方差平滑器:
其中,Pk-1表示k-1时刻的滤波均方误差,Qk-1表示k-1时刻***噪声协方差,表示k-1时刻状态转移矩阵的转置,表示k-1时刻***噪声系数矩阵的转置,εk表示k时刻的新息序列,和分别是k和k-1时刻的新息协方差N步平滑值,N=10,且当k<N时,
(1)建立k时刻新息协方差的限定窗口平滑器其递推形式:
(2)修正k时刻的滤波增益矩阵Kk,其计算公式写作:
其中,Pk,k-1表示k时刻一步预测均方误差,Hk表示k时刻量测系数矩阵。
(1)得到k时刻的调节因子ρk,其表达式写作:
其中,max{□}表示取内部所有元素的最大值,tr[□]表示对内部矩阵求迹;
(2)修正一步预测均方误差Pk,k-1,其表达式写作:
其中,Φk-1表示k-1时刻的***状态转移矩阵,Pk-1表示k-1时刻的滤波估计均方误差,Γk-1表示k-1时刻的***噪声系数阵,Qk-1表示k-1时刻的***噪声协方差。
步骤六:将步骤三至步骤五中得到的滤波参数带入到自适应滤波器中,并不断重复以上过程完成对导航参数误差的估计;利用输出校正实时修正捷联惯导***输出的导航信息,从而得到更高精度的导航参数信息,即补偿后的位置、速度和姿态信息,直至捷联惯性/多普勒组合导航过程结束。
(1)自适应滤波器对导航参数误差的估计值:
k时刻状态估计值:
k时刻估计均方误差: Pk=(I-KkHk)Pk,k-1
(2)得到导航参数误差的估计值后,利用输出校正对捷联惯导***输出的导航参数进行补偿:
①位置修正:λ=λSINS-δλ
L=LSINS-δL
其中,λ、L分别表示修正后的经纬度信息,λSINS、LSINS分别表示捷联惯导***输出的经纬度信息,δλ、δL分别表示滤波器估计的经纬度误差值。
②速度修正: vE=vSINS-E-δvE
vN=vSINS-N-δvN
其中,vE、vN分别表示修正后的东北向速度信息,vSINS-E、vSINS-N分别表示捷联惯导***输出的东北向速度信息,δvE、δvN分别表示滤波器估计的东北向速度误差值。
③姿态修正:
其中,φx、φy和φz分别表示计算导航系与真实导航系坐标轴之间的失准角。
b)对捷联惯导***解算的姿态角矩阵进行修正
c)计算补偿后的姿态角,再根据
至此,得到的θ,γ和即为通过自适应卡尔曼滤波修正后的姿态角信息。
应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。
Claims (5)
1.一种用于捷联惯性/多普勒组合导航***的自适应滤波方法,其特征在于,包括以下步骤:
步骤一:利用捷联惯性导航***中惯性测量单元测得沿载体系相对惯性空间的角速率和加速度分量信息,并与***的初始位置、速度以及姿态信息一起进行捷联惯导解算,得到载体实时的位置、速度和姿态信息;
步骤二:根据捷联惯导***的误差方程选取状态变量,建立***的状态方程,并选取经纬度误差、速度误差以及失准角作为状态变量;利用多普勒导航***与捷联惯导***所提供速度的差值作为量测变量,建立组合导航***的量测方程;
步骤三:经过连续***离散化处理后,对组合导航***进行滤波估计;首先更新状态一步预测值及其均方误差Pk,k-1,并利用k时刻量测信息Zk和滤波器的一步预测值计算k时刻的新息序列εk;再确定平滑窗口的宽度N,并借助量测新息序列εk建立关于新息协方差的N步限定窗口平滑器
所述的一步预测值:
所述的一步预测均方误差:
所述的N步新息协方差平滑器:
其中,Pk-1表示k-1时刻的滤波均方误差,Qk-1表示k-1时刻***噪声协方差,表示k-1时刻状态转移矩阵的转置,表示k-1时刻***噪声系数矩阵的转置,εk表示k时刻的新息序列,和分别是k和k-1时刻的新息协方差N步平滑值,N=10,且当k<N时,
步骤四:借助步骤三中得到的关于新息协方差的限定窗口平滑器对k时刻滤波器的增益矩阵Kk进行修正,使得增益矩阵能够根据新息序列的变化直接做出调整;
步骤六:将步骤三至步骤五中得到的滤波参数带入到自适应滤波器中,并不断重复以上过程完成对导航参数误差的估计;利用输出校正实时修正捷联惯导***输出的导航信息,从而得到更高精度的导航参数信息,即补偿后的位置、速度和姿态信息,直至捷联惯性/多普勒组合导航过程结束。
4.根据权利要求3所述的方法,其特征在于,所述步骤五的具体方法为:
(1)得到k时刻的调节因子ρk,其表达式写作:
其中,max{□}表示取内部所有元素的最大值,tr[□]表示对内部矩阵求迹;
(2)修正一步预测均方误差Pk,k-1,其表达式写作:
其中,Φk-1表示k-1时刻的***状态转移矩阵,Pk-1表示k-1时刻的滤波估计均方误差,Γk-1表示k-1时刻的***噪声系数阵,Qk-1表示k-1时刻的***噪声协方差。
5.根据权利要求4所述的方法,其特征在于,所述步骤六的具体方法为:
(1)自适应滤波器对导航参数误差的估计值:
k时刻状态估计值:
k时刻估计均方误差: Pk=(I-KkHk)Pk,k-1
(2)得到导航参数误差的估计值后,利用输出校正对捷联惯导***输出的导航参数进行补偿:
①位置修正: λ=λSINS-δλ
L=LSINS-δL
其中,λ、L分别表示修正后的经纬度信息,λSINS、LSINS分别表示捷联惯导***输出的经纬度信息,δλ、δL分别表示滤波器估计的经纬度误差值;
②速度修正: vE=vSINS-E-δvE
vN=vSINS-N-δvN
其中,vE、vN分别表示修正后的东北向速度信息,vSINS-E、vSINS-N分别表示捷联惯导***输出的东北向速度信息,δvE、δvN分别表示滤波器估计的东北向速度误差值;
③姿态修正:
其中,φx、φy和φz分别表示计算导航系与真实导航系坐标轴之间的失准角;
c)计算补偿后的姿态角,再根据
由于航向角θ,横摇角γ,纵摇角的定义域分别为:θ∈[0,2π],
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2013103136494A CN103389095A (zh) | 2013-07-24 | 2013-07-24 | 一种用于捷联惯性/多普勒组合导航***的自适应滤波方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2013103136494A CN103389095A (zh) | 2013-07-24 | 2013-07-24 | 一种用于捷联惯性/多普勒组合导航***的自适应滤波方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103389095A true CN103389095A (zh) | 2013-11-13 |
Family
ID=49533440
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2013103136494A Pending CN103389095A (zh) | 2013-07-24 | 2013-07-24 | 一种用于捷联惯性/多普勒组合导航***的自适应滤波方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103389095A (zh) |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103697887A (zh) * | 2013-12-05 | 2014-04-02 | 东南大学 | 基于捷联惯性制导和多普勒计程仪的优化导航方法 |
CN104020671A (zh) * | 2014-05-30 | 2014-09-03 | 哈尔滨工程大学 | 一种量测干扰下用于飞行器姿态估计的鲁棒递推滤波方法 |
CN104061930A (zh) * | 2013-12-05 | 2014-09-24 | 东南大学 | 基于捷联惯性制导和多普勒计程仪的导航方法 |
CN104215244A (zh) * | 2014-08-22 | 2014-12-17 | 南京航空航天大学 | 基于发射惯性坐标系的空天飞行器组合导航鲁棒滤波方法 |
CN107764257A (zh) * | 2017-09-14 | 2018-03-06 | 中国电子科技集团公司第五十四研究所 | 一种惯性器件数值模拟方法 |
CN108088498A (zh) * | 2017-11-29 | 2018-05-29 | 安徽省通信息科技有限公司 | 一种航向角与转动阻力***联合估计方法 |
CN108663068A (zh) * | 2018-03-20 | 2018-10-16 | 东南大学 | 一种应用在初始对准中的svm自适应卡尔曼滤波方法 |
CN108896036A (zh) * | 2018-05-09 | 2018-11-27 | 中国人民解放军国防科技大学 | 一种基于新息估计的自适应联邦滤波方法 |
CN108957495A (zh) * | 2018-05-03 | 2018-12-07 | 广州中海达卫星导航技术股份有限公司 | Gnss与mimu组合导航方法 |
CN109142779A (zh) * | 2018-08-09 | 2019-01-04 | 东莞市诺丽电子科技有限公司 | 一种车辆adas及dsm产品车速采集***及采集方法 |
CN110006427A (zh) * | 2019-05-20 | 2019-07-12 | 中国矿业大学 | 一种低动态高振动环境下的bds/ins紧组合导航方法 |
CN110146075A (zh) * | 2019-06-06 | 2019-08-20 | 哈尔滨工业大学(威海) | 一种增益补偿自适应滤波的sins/dvl组合定位方法 |
CN110763228A (zh) * | 2019-10-08 | 2020-02-07 | 哈尔滨工程大学 | 基于海底油气管节点位置的组合导航***误差修正方法 |
CN110763872A (zh) * | 2019-11-21 | 2020-02-07 | 中国船舶重工集团公司第七0七研究所 | 一种多普勒测速仪多参数在线标定方法 |
CN110823213A (zh) * | 2018-08-14 | 2020-02-21 | 北京自动化控制设备研究所 | 一种提高sins/dr组合导航***相对航向角精度的方法 |
CN111024074A (zh) * | 2019-12-12 | 2020-04-17 | 北京遥测技术研究所 | 一种基于递推最小二乘参数辨识的惯导速度误差确定方法 |
CN112636719A (zh) * | 2020-12-17 | 2021-04-09 | 郑州轻工业大学 | 数据丢失和信道噪声干扰下的ilc***输入信号滤波方法 |
CN112710304A (zh) * | 2020-12-17 | 2021-04-27 | 西北工业大学 | 一种基于自适应滤波的水下自主航行器导航方法 |
CN113175926A (zh) * | 2021-04-21 | 2021-07-27 | 哈尔滨工程大学 | 一种基于运动状态监测的自适应水平姿态测量方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH02116711A (ja) * | 1988-10-27 | 1990-05-01 | Tech Res & Dev Inst Of Japan Def Agency | ハイブリッド航法の誤差回避機能をもつ慣性航法装置 |
CN101393025A (zh) * | 2008-11-06 | 2009-03-25 | 哈尔滨工程大学 | Auv组合导航***无迹切换方法 |
CN101464152A (zh) * | 2009-01-09 | 2009-06-24 | 哈尔滨工程大学 | 一种sins/gps组合导航***自适应滤波方法 |
-
2013
- 2013-07-24 CN CN2013103136494A patent/CN103389095A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH02116711A (ja) * | 1988-10-27 | 1990-05-01 | Tech Res & Dev Inst Of Japan Def Agency | ハイブリッド航法の誤差回避機能をもつ慣性航法装置 |
CN101393025A (zh) * | 2008-11-06 | 2009-03-25 | 哈尔滨工程大学 | Auv组合导航***无迹切换方法 |
CN101464152A (zh) * | 2009-01-09 | 2009-06-24 | 哈尔滨工程大学 | 一种sins/gps组合导航***自适应滤波方法 |
Non-Patent Citations (4)
Title |
---|
WANG XIN-LONG等: "An innovative scheme for SINS/GPS ultra-tight integration system with low-grade IMU", 《AEROSPACE SCIENCE AND TECHNOLOGY》, 31 December 2012 (2012-12-31) * |
唐康华等: "基于自适应滤波的水下SINS/相控阵DVL组合导航算法设计", 《中国惯性技术学报》, vol. 21, no. 1, 1 March 2013 (2013-03-01) * |
李旦等: "组合导航自适应卡尔曼滤波改进算法研究", 《测控技术》, vol. 30, no. 3, 31 December 2011 (2011-12-31) * |
范科: "自适应滤波在组合导航和初始对准中的应用研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》, no. 2, 15 December 2011 (2011-12-15) * |
Cited By (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103697887A (zh) * | 2013-12-05 | 2014-04-02 | 东南大学 | 基于捷联惯性制导和多普勒计程仪的优化导航方法 |
CN104061930A (zh) * | 2013-12-05 | 2014-09-24 | 东南大学 | 基于捷联惯性制导和多普勒计程仪的导航方法 |
CN103697887B (zh) * | 2013-12-05 | 2017-03-01 | 东南大学 | 一种基于捷联惯导***和多普勒计程仪的优化导航方法 |
CN104061930B (zh) * | 2013-12-05 | 2017-06-16 | 东南大学 | 一种基于捷联惯性制导和多普勒计程仪的导航方法 |
CN104020671A (zh) * | 2014-05-30 | 2014-09-03 | 哈尔滨工程大学 | 一种量测干扰下用于飞行器姿态估计的鲁棒递推滤波方法 |
CN104020671B (zh) * | 2014-05-30 | 2017-01-11 | 哈尔滨工程大学 | 一种量测干扰下用于飞行器姿态估计的鲁棒递推滤波方法 |
CN104215244A (zh) * | 2014-08-22 | 2014-12-17 | 南京航空航天大学 | 基于发射惯性坐标系的空天飞行器组合导航鲁棒滤波方法 |
CN107764257B (zh) * | 2017-09-14 | 2019-10-22 | 中国电子科技集团公司第五十四研究所 | 一种惯性器件数值模拟方法 |
CN107764257A (zh) * | 2017-09-14 | 2018-03-06 | 中国电子科技集团公司第五十四研究所 | 一种惯性器件数值模拟方法 |
CN108088498A (zh) * | 2017-11-29 | 2018-05-29 | 安徽省通信息科技有限公司 | 一种航向角与转动阻力***联合估计方法 |
CN108663068A (zh) * | 2018-03-20 | 2018-10-16 | 东南大学 | 一种应用在初始对准中的svm自适应卡尔曼滤波方法 |
CN108957495A (zh) * | 2018-05-03 | 2018-12-07 | 广州中海达卫星导航技术股份有限公司 | Gnss与mimu组合导航方法 |
CN108957495B (zh) * | 2018-05-03 | 2020-12-29 | 广州中海达卫星导航技术股份有限公司 | Gnss与mimu组合导航方法、装置及计算机设备 |
CN108896036A (zh) * | 2018-05-09 | 2018-11-27 | 中国人民解放军国防科技大学 | 一种基于新息估计的自适应联邦滤波方法 |
CN108896036B (zh) * | 2018-05-09 | 2021-01-22 | 中国人民解放军国防科技大学 | 一种基于新息估计的自适应联邦滤波方法 |
CN109142779A (zh) * | 2018-08-09 | 2019-01-04 | 东莞市诺丽电子科技有限公司 | 一种车辆adas及dsm产品车速采集***及采集方法 |
CN109142779B (zh) * | 2018-08-09 | 2021-06-15 | 东莞市诺丽电子科技有限公司 | 一种车辆adas及dsm产品车速采集***及采集方法 |
CN110823213B (zh) * | 2018-08-14 | 2022-07-08 | 北京自动化控制设备研究所 | 一种提高sins/dr组合导航***相对航向角精度的方法 |
CN110823213A (zh) * | 2018-08-14 | 2020-02-21 | 北京自动化控制设备研究所 | 一种提高sins/dr组合导航***相对航向角精度的方法 |
CN110006427B (zh) * | 2019-05-20 | 2020-10-27 | 中国矿业大学 | 一种低动态高振动环境下的bds/ins紧组合导航方法 |
CN110006427A (zh) * | 2019-05-20 | 2019-07-12 | 中国矿业大学 | 一种低动态高振动环境下的bds/ins紧组合导航方法 |
CN110146075B (zh) * | 2019-06-06 | 2022-06-21 | 哈尔滨工业大学(威海) | 一种增益补偿自适应滤波的sins/dvl组合定位方法 |
CN110146075A (zh) * | 2019-06-06 | 2019-08-20 | 哈尔滨工业大学(威海) | 一种增益补偿自适应滤波的sins/dvl组合定位方法 |
CN110763228A (zh) * | 2019-10-08 | 2020-02-07 | 哈尔滨工程大学 | 基于海底油气管节点位置的组合导航***误差修正方法 |
CN110763872A (zh) * | 2019-11-21 | 2020-02-07 | 中国船舶重工集团公司第七0七研究所 | 一种多普勒测速仪多参数在线标定方法 |
CN111024074A (zh) * | 2019-12-12 | 2020-04-17 | 北京遥测技术研究所 | 一种基于递推最小二乘参数辨识的惯导速度误差确定方法 |
CN112710304A (zh) * | 2020-12-17 | 2021-04-27 | 西北工业大学 | 一种基于自适应滤波的水下自主航行器导航方法 |
CN112636719A (zh) * | 2020-12-17 | 2021-04-09 | 郑州轻工业大学 | 数据丢失和信道噪声干扰下的ilc***输入信号滤波方法 |
CN112710304B (zh) * | 2020-12-17 | 2022-12-13 | 西北工业大学 | 一种基于自适应滤波的水下自主航行器导航方法 |
CN112636719B (zh) * | 2020-12-17 | 2023-10-13 | 郑州轻工业大学 | 数据丢失和信道噪声干扰下的ilc***输入信号滤波方法 |
CN113175926A (zh) * | 2021-04-21 | 2021-07-27 | 哈尔滨工程大学 | 一种基于运动状态监测的自适应水平姿态测量方法 |
CN113175926B (zh) * | 2021-04-21 | 2022-06-21 | 哈尔滨工程大学 | 一种基于运动状态监测的自适应水平姿态测量方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103389095A (zh) | 一种用于捷联惯性/多普勒组合导航***的自适应滤波方法 | |
CN111024064B (zh) | 一种改进Sage-Husa自适应滤波的SINS/DVL组合导航方法 | |
CN110579740B (zh) | 一种基于自适应联邦卡尔曼滤波的无人船组合导航方法 | |
US20220404152A1 (en) | Motion constraint-aided underwater integrated navigation method employing improved sage-husa adaptive filtering | |
WO2020062791A1 (zh) | 一种深海潜航器的sins/dvl水下抗晃动对准方法 | |
CN103389506A (zh) | 一种用于捷联惯性/北斗卫星组合导航***的自适应滤波方法 | |
CN102608596B (zh) | 一种用于机载惯性/多普勒雷达组合导航***的信息融合方法 | |
CN112097763B (zh) | 一种基于mems imu/磁力计/dvl组合的水下运载体组合导航方法 | |
CN102829777B (zh) | 自主式水下机器人组合导航***及方法 | |
CN102519450B (zh) | 一种用于水下滑翔器的组合导航装置及方法 | |
CN103955218B (zh) | 一种基于非线性控制理论的无人艇轨迹跟踪控制装置及方法 | |
CN103759742B (zh) | 基于模糊自适应控制技术的捷联惯导非线性对准方法 | |
CN109000640B (zh) | 基于离散灰色神经网络模型的车辆gnss/ins组合导航方法 | |
CN100575879C (zh) | 自适应闭环h∞滤波器对北斗/捷联导航***的修正方法 | |
CN103063212B (zh) | 一种基于非线性映射自适应混合Kalman/H∞滤波器的组合导航方法 | |
CN105891863B (zh) | 一种基于高度约束的扩展卡尔曼滤波定位方法 | |
CN101900573B (zh) | 一种实现陆用惯性导航***运动对准的方法 | |
CN102508278A (zh) | 一种基于观测噪声方差阵估计的自适应滤波方法 | |
CN103017755A (zh) | 一种水下导航姿态测量方法 | |
CN103278837A (zh) | 基于自适应滤波的sins/gnss多级容错组合导航方法 | |
CN101949703A (zh) | 一种捷联惯性/卫星组合导航滤波方法 | |
CN104567930A (zh) | 一种能够估计和补偿机翼挠曲变形的传递对准方法 | |
CN104215259A (zh) | 一种基于地磁模量梯度和粒子滤波的惯导误差校正方法 | |
CN105424036A (zh) | 一种低成本水下潜器地形辅助惯性组合导航定位方法 | |
CN103822633A (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 | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20131113 |