CN114216480B - 一种惯性导航系统精对准方法 - Google Patents
一种惯性导航系统精对准方法 Download PDFInfo
- Publication number
- CN114216480B CN114216480B CN202111410043.3A CN202111410043A CN114216480B CN 114216480 B CN114216480 B CN 114216480B CN 202111410043 A CN202111410043 A CN 202111410043A CN 114216480 B CN114216480 B CN 114216480B
- Authority
- CN
- China
- Prior art keywords
- matrix
- carrier
- ins
- moment
- initial
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
- G01C25/005—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass initial alignment, calibration or starting-up of inertial devices
Abstract
本发明公开了一种惯性导航系统精对准方法,包括初始化INS/GPS系统;采集INS和GPS的输出数据;对惯性组件输出进行导航解算,以位置误差为观测量,INS导航误差为状态量,进行卡尔曼滤波时间更新;根据时间更新结果计算新息二阶矩,与判决条件进行比较,当满足判决条件时,使用自适应因子计算k时刻量测噪声阵,否则,基于Allan方差估计k时刻量测噪声阵;依据k时刻量测噪声阵进行卡尔曼滤波量测更新,对INS导航误差进行估计,补偿INS失准角,提高初始对准精度。本发明解决了受工况影响测量噪声发生变化时姿态角估计不准确问题,消除了由变化测量噪声引起的初始对准误差,提高INS定位精度。
Description
技术领域
本发明属于惯性技术领域,涉及一种惯性导航系统精对准方法,特别是一种基于新息检测进行自适应方法调整的惯性导航系统精对准方法。
背景技术
初始对准是为惯性导航系统(Inertial Navigation System,INS)提供初始姿态角,对准精度和对准时间是影响INS系统后续工作性能的两项重要指标。初始对准分粗对准和精对准两个阶段:粗对准阶段粗略估计初始姿态角;精对准在粗对准基础上提高对准精度,因此初始对准最终精度主要与精对准精度有关。基于GPS辅助的INS精对准方法以GPS位置信息为基准,获取INS位置误差,根据INS误差模型,通过卡尔曼滤波对INS姿态角校正。然而,受工况影响,INS测量噪声统计特性随时间发生变化,导致卡尔曼滤波的量测噪声阵不适用,影响卡尔曼滤波性能,降低INS失准角的估计精度,进而影响INS初始对准精度。
授权号为CN107479076B的中国发明专利在2017年8月8日公开的《一种动基座下联合滤波初始对准方法》中,提出了一种根据失准角判决门限实现无迹卡尔曼滤波和卡尔曼滤波算法切换的INS初始对准方法,缩短了对准时间。公开号为CN112857398A的中国发明专利在2021年05月28日公开的《一种系泊状态下舰船的快速初始对准方法》中,对初始对准误差模型进行了改进,利用四元数无味估计器进行滤波估计,提升对准精度。《武汉大学学报·信息科学版》期刊2020年45卷7期,由查峰等人撰写的《一种改进的捷联惯导参数辨识初始对准方法》一文中,采用泰勒展开方式对地球自转角的三角函数进行高阶近似,减少因地球自转角的正弦函数一阶近似对参数辨识模型的影响,提高初始对准精度。《中国惯性技术学报》期刊2020年28卷6期,由徐博等人撰写的《基于粒子群优化的舰船捷联惯导初始对准方法》一文中,提出了一种基于惯性系的粒子群算法优化方法,降低舰船基座摇摆对对准性能影响,提高航向角误差收敛速度,缩短对准时间。《IEEE TRANSACTIONS ONINSTRUMENTATION AND MEASUREMENT》期刊2020年69卷10期,由Yusen Lin等人撰写的《AnImproved MCMC-Based Particle Filter for GPS-Aided SINS In-Motion InitialAlignment》一文中,提出了一种改进的基于马尔科夫链蒙特卡罗的粒子滤波算法,通过避免粒子滤波中的粒子退化和粒子贫化现象,提高滤波精度,减少计算量,缩短对准时间。《IEEE TRANSACTIONS ON INSTRUMENTATION AND MEASUREMENT》期刊2020年70卷,由Xiaoren Zhou等人撰写的《A Robust Quaternion Kalman Filter Method for MIMU/GPSIn-Motion Alignment》一文中,提出了一种新的鲁棒四元数卡尔曼滤波算法,消除GPS异常值对对准精度的影响。以上技术都通过改进滤波模型或改变滤波方法的角度,提高初始对准精度、缩短对准时间,增加了算法的复杂度。
发明内容
针对上述现有技术,本发明要解决的技术问题是提供一种基于新息检测进行自适应方法调整的惯性导航系统精对准方法,解决测量噪声发生变化情况下的姿态角估计不准确问题,消除由变化测量噪声引起的初始对准误差。
为解决上述技术问题,本发明的一种惯性导航系统精对准方法,包括以下步骤:
步骤一:初始化INS/GPS组合系统;
步骤二:采集惯性组件输出的角速度ωb、加速度fb和GPS输出的位置pgps;
步骤三:对步骤二中采集的ωb和fb进行导航解算,得到载体的位置速度/>姿态/>
步骤四:构造卡尔曼滤波器,滤波状态量为:X=[δpT δvT δφT ΔT εT]T,其中δp为INS位置误差,δp=[δpx δpy δpz]T,δv为速度误差,δv=[δvx δvy δvz]T,δφ为姿态失准角,δφ=[δφx δφy δφz]T,Δ为加速度计零偏,Δ=[Δx Δy Δz]T,ε为陀螺零偏,ε=[εx εyεz]T,以GPS输出的位置pgps和步骤三中解算位置之差为观测量Z(k),进行卡尔曼滤波时间更新,计算状态一步预测阵X(k/k-1)和一步预测均方误差阵P(k/k-1);
步骤五:根据步骤四的Zk和X(k/k-1),计算k时刻新息二阶矩C(k),即
C(k)=E((Z(k)-HX(k/k-1))(Z(k)-HX(k/k-1))T)
步骤六:自适应方法选取:比较C(k)与HP(k/k-1)HT+R(k-1),同时比较k与M,当C(k)=HP(k/k-1)HT+R(k-1),且k≥M时,则转至步骤七,否则转至步骤十,其中H表示量测矩阵,R(k-1)表示k-1时刻量测噪声阵;
步骤七:计算(k-M+1)~k时刻M个点新息二阶矩的均值即
其中,C(j)表示j时刻新息二阶矩;
步骤八:根据步骤五计算的C(k)和步骤七计算的计算自适应因子α,即
其中,tr表示对矩阵求迹;
步骤九:根据步骤八中的α,计算k时刻R(k),即:
R(k)=(1+α)R(k-1)
计算完成后,转至步骤十二。
步骤十:计算k时刻指数渐消记忆加权因子β(k),即:
其中,b表示遗忘因子;
步骤十一:根据步骤四的Z(k)和步骤十的β(k),基于Allan方差估计k时刻R(k),即:
计算完成后,转至步骤十二;
步骤十二:进行卡尔曼滤波测量更新,计算k时刻滤波增益K(k)、估计均方误差P(k)和状态估计值X(k);
步骤十三:根据步骤十二中得到的状态估计值X(k),对步骤三中得到的载体姿态进行补偿,得到载体姿态角φ;存储并输出载体姿态角。
进一步的,步骤一中初始化INS/GPS组合系统包括:
初始化INS导航解算初值:载体的位置信息λ0和h0分别为初始纬度、初始经度和初始高度,三维速度信息vx0、vy0、vz0,以及三个粗对准姿态角信息φx0、φy0、φz0,初始转换矩阵/>初始化四元数q0;
初始转换矩阵满足:
其中,b表示载体坐标系,n表示导航坐标系,表示b系到n系的转换矩阵;
初始化四元数q0满足:
其中,q0=[q00 q10 q20 q30]T,cbnij(i,j=1,2,3)表示中第i行、第j列矩阵元素;
初始化常值参量,包括新息二阶矩均值窗长M和遗忘因子b;
初始化指数渐消记忆加权因子初值β0;
初始化卡尔曼滤波器初值:状态变量初值X0=[δpT δvT δφT ΔT εT]T,均方误差阵P0,系统噪声方差阵Q,量测噪声方差阵R0,量测阵Hk,Hk=[Ι3×3 Ο3×3 Ο3×3 Ο3×3 Ο3×3],其中,Ι表示单位矩阵;Ο表示零矩阵。
进一步的,步骤三中对步骤二中采集的ωb和fb进行导航解算,得到载体的位置速度/>姿态/>包括:
四元数姿态矩阵更新运算:
其中,表示INS解算纬度;ωie表示地球自转角速度;/>分别表示INS解算载体速度在导航系oxn轴、oyn轴轴上的分量;Re表示地球半径;T为采样时间;/>当k=1时,q(1)为步骤1中的载体初始四元数;
根据q(k)中元素qi(k),i=0,1,2,3,更新捷联矩阵
其中,上式中的qi,i=0,1,2,3为式(9)中qi(k)(i=0,1,2,3);
更新姿态信息:
其中,i,j=1,2,3表示/>中第i行、第j列矩阵元素;则姿态/> 和/>分别表示横滚角、俯仰角、航向角解算值;
速度更新运算:
其中,gn=[0 0 -g]T,g表示当地重力加速度,速度和/>分别表示解算载体速度在导航系oxn轴、oyn轴、ozn轴上的分量;当k=1时,/>为步骤一中初始化系统时获得的载体初始速度;
位置更新计算:
位置表示解算载体纬度,/>表示解算载体经度,/>表示解算载体高度;当t=1时,/>为步骤一中载体初始位置。
进一步的,步骤四中计算状态一步预测阵X(k/k-1)和一步预测均方误差阵P(k/k-1)具体为:
X(k/k-1)=F(k/k-1)X(k-1)
其中,X(k/k-1)表示状态一步预测阵;F(k/k-1)表示k-1时刻到k时刻的状态转移矩阵;X(k-1)表示k-1时刻状态估计值;
P(k/k-1)=F(k/k-1)P(k-1)FT(k/k-1)+Q
其中,P(k/k-1)表示一步预测均方误差阵;P(k-1)表示k-1时刻状态估计均方误差;Q表示系统噪声。
进一步的,步骤十二中进行卡尔曼滤波测量更新,计算k时刻滤波增益K(k)、估计均方误差P(k)和状态估计值X(k)具体为:
滤波增益K(k)满足:
K(k)=P(k,k-1)HT(HP(k,k-1)HT+R(k))
估计均方误差P(k)满足:
P(k)=(I15×15-K(k)H)P(k/k-1)
状态估计值X(k)满足:
X(k)=X(k/k-1)+K(k)(Z(k)-HX(k/k-1))
其中,X(k)表示k时刻状态估计值。
进一步的,步骤十三中根据步骤十二中得到的状态估计值X(k),对步骤三中得到的四元数进行补偿,得到载体姿态角φ具体为:
对步骤三中得到的四元数进行补偿:
根据上式中元素,计算更新捷联矩阵更新姿态信息:
其中,cij(i,j=1,2,3)表示中第i行、第j列矩阵元素;,φi(i=x,y,z)分别表示横滚角、俯仰角、航向角解算值。
本发明的有益效果:本发明针对初始精对准过程中由测量噪声统计特性变化引起的对准误差增大问题,提出了一种基于新息检测进行自适应方法调整的INS精对准方法,以位置误差为观测量,INS导航误差为状态量,进行卡尔曼滤波时间更新;根据时间更新结果计算新息二阶矩,与判决条件进行比较,当满足判决条件时,使用自适应因子计算k时刻量测噪声阵,否则,基于Allan方差估计k时刻量测噪声阵;依据k时刻量测噪声阵进行卡尔曼滤波量测更新,对INS导航误差进行估计,补偿INS失准角,提高初始对准精度。
本发明在进行卡尔曼滤波时,根据新息矩阵进行自适应R阵方法调整,对卡尔曼滤波测量噪声阵进行估计,经过进一步解算得到姿态失准角估计值,对姿态角进行校正,消除了由变化测量噪声引起的姿态角估计不准确问题,提高INS精对准精度。其优点在于:(1)解决了受工况影响,测量噪声发生变化情况下的姿态角估计不准确问题,消除了由变化测量噪声引起的初始对准误差;(2)不需要任何外界辅助信息就可提高定位精度;(3)计算量小,简单易操作。
附图说明
图1为本发明的方法流程图;
图2为利用本发明进行仿真试验结果;
图3为利用本发明进行实测试验结果。
具体实施方式
下面结合说明书附图和具体实施例对本发明做进一步说明。
结合图1,本发明包括以下步骤:
步骤一:上电,初始化INS/GPS组合系统。导航初始时刻,初始化系统:
(1)初始化INS导航解算初值:船的位置信息λ0(单位均为度),h0(单位为m)分别为初始纬度、初始经度和初始高度,三维速度信息vx0、vy0、vz0(单位均为m/s),以及三个粗对准姿态角信息φx0、φy0、φz0(单位均为度),初始转换矩阵/>初始化四元数q0;
(2)初始化常值参量:新息二阶矩均值窗长M;遗忘因子b;
(3)初始化指数渐消记忆加权因子初值:β0;
(4)初始化卡尔曼滤波器初值:状态变量初值X0=]δpT δvT δφT ΔT εT]T,均方误差阵P0,系统噪声方差阵Q,量测噪声方差阵R0,量测阵Hk,Hk=[Ι3×3 Ο3×3 Ο3×3 Ο3×3Ο3×3],其余初值根据实际情况设定。其中,Ι表示单位矩阵;Ο表示零矩阵。
初始转换矩阵计算如下:
其中,b表示载体坐标系,n表示导航坐标系,表示b系到n系的转换矩阵。
初始化四元数q0计算如下:
其中,q0=[q00 q10 q20 q30]T,cbnij(i,j=1,2,3)表示中第i行、第j列矩阵元素;
将以上初始化信息装订至导航计算机中。
导航过程中,利用该初始信息进行更新,得到任意时刻载体的位置、速度和姿态信息。
步骤二:实时采集惯性组件加速度计输出的三轴加速度陀螺输出的三轴角速度/>和GPS输出的位置信息/>其中,fi b(i=x,y,z)表示加速度fb在载体系oxb轴、oyb轴、ozb轴上的分量(单位均为m/s2);分别表示角速度ωb在载体系oxb轴、oyb轴、ozb轴上的分量(单位均为rad/s);表示GPS输出纬度,λgps表示GPS输出经度,hgps表示GPS输出高度;角标T表示矩阵转置。
步骤三:对步骤二中获取的惯性组件输出数据进行导航解算:
四元数姿态矩阵更新运算:
其中,表示INS解算纬度;ωie表示地球自转角速度;/>分别表示INS解算载体速度在导航系oxn轴、oyn轴轴上的分量;Re表示地球半径;T为采样时间;/>当k=1时,q(1)为步骤1中的载体初始四元数。
根据q(k)中元素qi(k)(i=0,1,2,3),更新捷联矩阵
其中,上式中的qi(i=0,1,2,3)为式(3)中qi(k)(i=0,1,2,3),更新姿态信息:
其中,表示/>中第i行、第j列矩阵元素;/> 分别表示横滚角、俯仰角、航向角解算值。
速度更新运算:
其中,分别表示解算载体速度在导航系oxn轴、oyn轴、ozn轴上的分量;gn=[0 0 -g]T,g表示当地重力加速度。当k=1时,/>为步骤一中初始化系统时获得的载体初始速度。
位置更新计算:
其中,表示解算载体纬度,/>表示解算载体经度,/>表示解算载体高度;当t=1时,/>为步骤一中载体初始位置。
至此,根据式(5)、(6)、(7)得到INS解算的载体姿态角、速度和位置。
步骤四:构造卡尔曼滤波器,滤波状态量:X=[δpT δvT δφT ΔT εT]T,以步骤二中GPS位置信息pgps和步骤三中解算位置信息之差为观测量,即:/>进行卡尔曼滤波时间更新。其中,δp表示INS位置误差,δp=[δpx δpy δpz]T;δv表示INS速度误差,δv=[δvx δvy δvz]T;δφ表示INS姿态失准角,δφ=[δφx δφy δφz]T;Δ表示加速度计零偏,Δ=[Δx Δy Δz]T;ε表示陀螺零偏,ε=[εx εy εz]T,具体更新过程如下:
状态一步预测阵运算:
X(k/k-1)=F(k/k-1)X(k-1) (8)
其中,X(k/k-1)表示状态一步预测阵;F(k/k-1)表示k-1时刻到k时刻的状态转移矩阵;X(k-1)表示k-1时刻状态估计值。
一步预测均方误差阵运算:
P(k/k-1)=F(k/k-1)P(k-1)FT(k/k-1)+Q (9)
其中,P(k/k-1)表示一步预测均方误差阵;P(k-1)表示k-1时刻状态估计均方误差;Q表示系统噪声。
步骤五:根据步骤四的Z(k)和X(k/k-1),计算k时刻新息二阶矩C(k),即
C(k)=E((Z(k)-HX(k/k-1))(Z(k)-HX(k/k-1))T) (10)
步骤六:自适应方法选取。比较C(k)与HP(k/k-1)HT+R(k-1),同时比较k与M,当C(k)=HP(k/k-1)HT+R(k-1),且k≥M时,则转至步骤七,否则转至步骤十。其中,H表示量测矩阵,R(k-1)表示k-1时刻量测噪声阵。
步骤七:计算(k-M+1)~k时刻,M个点新息二阶矩均值即
其中,C(j)表示j时刻新息二阶矩。
步骤八:根据步骤五的C(k)和步骤七的计算自适应因子α,即
其中,tr表示对矩阵求迹。
步骤九:根据步骤八中的自适应因子,计算k时刻R(k),即:
R(k)=(1+α)R(k-1) (13)
计算完成后,转至步骤十二。
步骤十:计算k时刻指数渐消记忆加权因子β(k),即:
其中,b表示遗忘因子。
步骤十一:根据步骤四的Z(k)和步骤十的β(k),基于Allan方差估计k时刻R(k),即:
计算完成后,转至步骤十二。
步骤十二:根据步骤四的Z(k)、X(k/k-1)、P(k,k-1)和步骤九或步骤十一的R(k),进行卡尔曼滤波测量更新。具体计算过程如下:
滤波增益计算:
K(k)=P(k,k-1)HT(HP(k,k-1)HT+R(k)) (16)
估计均方误差计算:
P(k)=(I15×15-K(k)H)P(k/k-1) (17)
状态估计值计算:
X(k)=X(k/k-1)+K(k)(Z(k)-HX(k/k-1)) (18)
其中,X(k)表示k时刻状态估计值。
步骤十三:根据步骤十二中得到的状态估计值X(k),对步骤三中得到的四元数进行补偿,得到载体姿态角φ。具体为:
根据上式中元素,结合公式计算更新捷联矩阵更新姿态信息:
其中,cij(i,j=1,2,3)表示中第i行、第j列矩阵元素;,φi(i=x,y,z)分别表示横滚角、俯仰角、航向角解算值;存储并输出载体姿态角。
下面结合具体参数对本发明的有益效果如下方式得以验证:
船舶航行路线:直线运动
(1)系统初始化参数如下:
舰船运动参数:舰船速度:vx0=0m/s,vy0=0m/s,vz0=0m/s;
舰船姿态:φx0=0°,φy0=0°,φz0=30°;
舰船姿态角误差:0.5°,0.5°,5°;
当地经纬度:λ0=126.6705°;
陀螺参数:陀螺常值漂移:0.01°/h
陀螺白噪声误差:0.001°/h
加速度计参数:加速度计零偏:10-5g
加速度计白噪声误差:10-6g
(2)常值参量:
采样时间:T=0.005s
新息二阶矩均值窗长:M=5
遗忘因子:b=0.95
(3)指数渐消记忆加权因子初值设置:β0=1
(4)卡尔曼滤波器参数设置:
状态变量初值:X0=[O1×3 O1×3 O1×3 O1×3 O1×3]T
均方误差阵:
Pa0=diag[2.12×10-8 2.12×10-8 3.58×10-4 10-2 10-2 10-4 2.46×10-14 2.46×10-14 1 2.12×10-14 2.12×10-14 2.12×10-14 9.57×10-7 9.57×10-7 9.57×10-7]
其中,diag表示矩阵对角线上元素。
系统噪声方差阵:
Q=diag[8.46×10-14 8.46×10-14 8.46×10-14 9.57×10-9 9.57×10-9 9.57×10-9 0 0 0 0 0 0 0 0 0]
量测噪声方差阵:Ra=diag[2.46×10-12 2.46×10-12 2.46×10-12]
利用本发明方法,得到应用此方法的失准角曲线,图2表示失准角曲线。因实测试验没有外姿态基准,因此采用在初始对准完成后进行组合导航,用组合导航的定位误差反映对准精度;图3表示有、无自适应调整R阵进行精对准再进行组合导航的定位误差曲线。结果表明本发明较好抑制了由测量噪声变化引起的初始对准精度下降,可以满足实际需求。
Claims (6)
1.一种惯性导航系统精对准方法,其特征在于,包括以下步骤:
步骤一:初始化INS/GPS组合系统;
步骤二:采集惯性组件输出的角速度ωb、加速度fb和GPS输出的位置pgps;
步骤三:对步骤二中采集的ωb和fb进行导航解算,得到载体的位置速度/>姿态/>
步骤四:构造卡尔曼滤波器,滤波状态量为:X=[δpT δvT δφT ΔT εT]T,其中δp为INS位置误差,δp=[δpx δpy δpz]T,δv为速度误差,δv=[δvx δvy δvz]T,δφ为姿态失准角,δφ=[δφx δφy δφz]T,Δ为加速度计零偏,Δ=[Δx Δy Δz]T,ε为陀螺零偏,ε=[εx εy εz]T,以GPS输出的位置pgps和步骤三中解算位置之差为观测量Z(k),进行卡尔曼滤波时间更新,计算状态一步预测阵X(k/k-1)和一步预测均方误差阵P(k/k-1);
步骤五:根据步骤四的Zk和X(k/k-1),计算k时刻新息二阶矩C(k),即
C(k)=E((Z(k)-HX(k/k-1))(Z(k)-HX(k/k-1))T)
步骤六:自适应方法选取:比较C(k)与HP(k/k-1)HT+R(k-1),同时比较k与M,当C(k)=HP(k/k-1)HT+R(k-1),且k≥M时,则转至步骤七,否则转至步骤十,其中H表示量测矩阵,R(k-1)表示k-1时刻量测噪声阵;
步骤七:计算(k-M+1)~k时刻M个点新息二阶矩的均值即
其中,C(j)表示j时刻新息二阶矩;
步骤八:根据步骤五计算的C(k)和步骤七计算的计算自适应因子α,即
其中,tr表示对矩阵求迹;
步骤九:根据步骤八中的α,计算k时刻R(k),即:
R(k)=(1+α)R(k-1)
计算完成后,转至步骤十二;
步骤十:计算k时刻指数渐消记忆加权因子β(k),即:
其中,b表示遗忘因子;
步骤十一:根据步骤四的Z(k)和步骤十的β(k),基于Allan方差估计k时刻R(k),即:
计算完成后,转至步骤十二;
步骤十二:进行卡尔曼滤波测量更新,计算k时刻滤波增益K(k)、估计均方误差P(k)和状态估计值X(k);
步骤十三:根据步骤十二中得到的状态估计值X(k),对步骤三中得到的载体姿态进行补偿,得到载体姿态角φ;存储并输出载体姿态角。
2.根据权利要求1所述的一种惯性导航系统精对准方法,其特征在于:步骤一所述初始化INS/GPS组合系统包括:
初始化INS导航解算初值:载体的位置信息 λ0和h0分别为初始纬度、初始经度和初始高度,三维速度信息vx0、vy0、vz0,以及三个粗对准姿态角信息φx0、φy0、φz0,初始转换矩阵/>初始化四元数q0;
初始转换矩阵满足:
其中,b表示载体坐标系,n表示导航坐标系,表示b系到n系的转换矩阵;
初始化四元数q0满足:
其中,q0=[q00 q10 q20 q30]T,cbnij(i,j=1,2,3)表示中第i行、第j列矩阵元素;
初始化常值参量,包括新息二阶矩均值窗长M和遗忘因子b;
初始化指数渐消记忆加权因子初值β0;
初始化卡尔曼滤波器初值:状态变量初值X0=[δpT δvT δφT ΔT εT]T,均方误差阵P0,系统噪声方差阵Q,量测噪声方差阵R0,量测阵Hk,Hk=[Ι3×3 Ο3×3 Ο3×3 Ο3×3 Ο3×3],其中,Ι表示单位矩阵;Ο表示零矩阵。
3.根据权利要求1所述的一种惯性导航系统精对准方法,其特征在于:步骤三所述对步骤二中采集的ωb和fb进行导航解算,得到载体的位置速度/>姿态/>包括:
四元数姿态矩阵更新运算:
其中, 表示INS解算纬度;ωie表示地球自转角速度; 分别表示INS解算载体速度在导航系oxn轴、oyn轴轴上的分量;Re表示地球半径;T为采样时间;/>当k=1时,q(1)为步骤1中的载体初始四元数;
根据q(k)中元素qi(k),i=0,1,2,3,更新捷联矩阵
其中,上式中的qi,i=0,1,2,3为式(9)中qi(k)(i=0,1,2,3);
更新姿态信息:
其中,表示/>中第i行、第j列矩阵元素;则姿态/> 和/>分别表示横滚角、俯仰角、航向角解算值;
速度更新运算:
其中,gn=[0 0 -g]T,g表示当地重力加速度,速度 和/>分别表示解算载体速度在导航系oxn轴、oyn轴、ozn轴上的分量;当k=1时,/>为步骤一中初始化系统时获得的载体初始速度;
位置更新计算:
位置 表示解算载体纬度,/>表示解算载体经度,/>表示解算载体高度;当t=1时,/>为步骤一中载体初始位置。
4.根据权利要求1所述的一种惯性导航系统精对准方法,其特征在于:步骤四所述计算状态一步预测阵X(k/k-1)和一步预测均方误差阵P(k/k-1)具体为:
X(k/k-1)=F(k/k-1)X(k-1)
其中,X(k/k-1)表示状态一步预测阵;F(k/k-1)表示k-1时刻到k时刻的状态转移矩阵;X(k-1)表示k-1时刻状态估计值;
P(k/k-1)=F(k/k-1)P(k-1)FT(k/k-1)+Q
其中,P(k/k-1)表示一步预测均方误差阵;P(k-1)表示k-1时刻状态估计均方误差;Q表示系统噪声。
5.根据权利要求1所述的一种惯性导航系统精对准方法,其特征在于:步骤十二所述进行卡尔曼滤波测量更新,计算k时刻滤波增益K(k)、估计均方误差P(k)和状态估计值X(k)具体为:
滤波增益K(k)满足:
K(k)=P(k,k-1)HT(HP(k,k-1)HT+R(k))
估计均方误差P(k)满足:
P(k)=(I15×15-K(k)H)P(k/k-1)
状态估计值X(k)满足:
X(k)=X(k/k-1)+K(k)(Z(k)-HX(k/k-1))
其中,X(k)表示k时刻状态估计值。
6.根据权利要求1所述的一种惯性导航系统精对准方法,其特征在于:步骤十三所述根据步骤十二中得到的状态估计值X(k),对步骤三中得到的载体姿态进行补偿,得到载体姿态角φ具体为:
对步骤三中得到的四元数进行补偿:
根据上式中元素,计算更新捷联矩阵更新姿态信息:
其中,cij(i,j=1,2,3)表示中第i行、第j列矩阵元素;φi(i=x,y,z)分别表示横滚角、俯仰角、航向角解算值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111410043.3A CN114216480B (zh) | 2021-11-25 | 2021-11-25 | 一种惯性导航系统精对准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111410043.3A CN114216480B (zh) | 2021-11-25 | 2021-11-25 | 一种惯性导航系统精对准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114216480A CN114216480A (zh) | 2022-03-22 |
CN114216480B true CN114216480B (zh) | 2023-07-21 |
Family
ID=80698241
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111410043.3A Active CN114216480B (zh) | 2021-11-25 | 2021-11-25 | 一种惯性导航系统精对准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114216480B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108663068A (zh) * | 2018-03-20 | 2018-10-16 | 东南大学 | 一种应用在初始对准中的svm自适应卡尔曼滤波方法 |
CN111024070A (zh) * | 2019-12-23 | 2020-04-17 | 哈尔滨工程大学 | 一种基于航向自观测的惯性足绑式行人定位方法 |
CN111795708A (zh) * | 2020-06-16 | 2020-10-20 | 湖南跨线桥航天科技有限公司 | 晃动基座条件下陆用惯性导航系统的自适应初始对准方法 |
WO2021110055A1 (zh) * | 2019-12-02 | 2021-06-10 | 东南大学 | 一种用于最优估计精对准的严格逆向导航方法 |
-
2021
- 2021-11-25 CN CN202111410043.3A patent/CN114216480B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108663068A (zh) * | 2018-03-20 | 2018-10-16 | 东南大学 | 一种应用在初始对准中的svm自适应卡尔曼滤波方法 |
WO2021110055A1 (zh) * | 2019-12-02 | 2021-06-10 | 东南大学 | 一种用于最优估计精对准的严格逆向导航方法 |
CN111024070A (zh) * | 2019-12-23 | 2020-04-17 | 哈尔滨工程大学 | 一种基于航向自观测的惯性足绑式行人定位方法 |
CN111795708A (zh) * | 2020-06-16 | 2020-10-20 | 湖南跨线桥航天科技有限公司 | 晃动基座条件下陆用惯性导航系统的自适应初始对准方法 |
Non-Patent Citations (1)
Title |
---|
自适应Kalman滤波在SINS初始对准中的应用;苏宛新;黄春梅;刘培伟;马明龙;;中国惯性技术学报(第01期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114216480A (zh) | 2022-03-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108827310B (zh) | 一种船用星敏感器辅助陀螺仪在线标定方法 | |
JP5237723B2 (ja) | 動的に較正されるセンサデータと、ナビゲーションシステム内の繰り返し拡張カルマンフィルタとを使用する、ジャイロコンパスの整合用のシステム及び方法 | |
CN111024064B (zh) | 一种改进Sage-Husa自适应滤波的SINS/DVL组合导航方法 | |
CN107655493B (zh) | 一种光纤陀螺sins六位置系统级标定方法 | |
CN109916395B (zh) | 一种姿态自主冗余组合导航算法 | |
CN110398257A (zh) | Gps辅助的sins系统快速动基座初始对准方法 | |
CN113063429B (zh) | 一种自适应车载组合导航定位方法 | |
CN110631574B (zh) | 一种惯性/里程计/rtk多信息融合方法 | |
CN111024070A (zh) | 一种基于航向自观测的惯性足绑式行人定位方法 | |
CN103822633A (zh) | 一种基于二阶量测更新的低成本姿态估计方法 | |
CN111238535B (zh) | 一种基于因子图的imu误差在线标定方法 | |
CN110954102B (zh) | 用于机器人定位的磁力计辅助惯性导航系统及方法 | |
CN112880669B (zh) | 一种航天器星光折射和单轴旋转调制惯性组合导航方法 | |
CN112945271A (zh) | 磁力计信息辅助的mems陀螺仪标定方法及标定系统 | |
CN116007620A (zh) | 一种组合导航滤波方法、系统、电子设备及存储介质 | |
CN110285830B (zh) | 基于mems传感器的sins/gps速度匹配对准方法 | |
CN113503892A (zh) | 一种基于里程计和回溯导航的惯导系统动基座初始对准方法 | |
CN111649747A (zh) | 一种基于imu的自适应ekf姿态测量改进方法 | |
CN109084755B (zh) | 一种基于重力视速度与参数辨识的加速度计零偏估计方法 | |
CN113008229B (zh) | 一种基于低成本车载传感器的分布式自主组合导航方法 | |
CN114216480B (zh) | 一种惯性导航系统精对准方法 | |
CN116222551A (zh) | 一种融合多种数据的水下导航方法及装置 | |
CN113551669B (zh) | 基于短基线的组合导航定位方法及装置 | |
CN114184190A (zh) | 一种惯性/里程计组合导航系统及方法 | |
CN114993346A (zh) | 适用于跨海空介质的捷联式惯性导航系统的初始对准方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |