CN110440830A - 动基座下车载捷联惯导***自对准方法 - Google Patents
动基座下车载捷联惯导***自对准方法 Download PDFInfo
- Publication number
- CN110440830A CN110440830A CN201910766894.8A CN201910766894A CN110440830A CN 110440830 A CN110440830 A CN 110440830A CN 201910766894 A CN201910766894 A CN 201910766894A CN 110440830 A CN110440830 A CN 110440830A
- Authority
- CN
- China
- Prior art keywords
- error
- matrix
- equation
- coordinate system
- carrier
- 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
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Navigation (AREA)
Abstract
本发明公开了一种动基座下车载捷联惯导***自对准方法,包括以下步骤:步骤1,获得捷联惯导***所在地理位置信息,其中地理位置信息包括经度λ、纬度重力加速度g;步骤2,利用经度、纬度、重力加速度信息及惯性器件输出信息,采用凝固坐标系,利用多矢量定姿方法进行粗对准,求得初始姿态矩阵;步骤3,根据捷联惯性***误差方程及惯性器件误差方程构建kalman滤波状态方程和量测方程;步骤4,设计混合校准kalman滤波方法估计***误差量和惯性器件误差量并修正,实现精对准过程,完成初始对准。本发明具有较强的抗干扰能力,保证了动基座条件下的自对准精度。
Description
技术领域
本发明特别涉及一种动基座下车载捷联惯导***自对准方法,涉及动基座下基于凝固坐标系粗对准、精对准kalman滤波状态空间模型的建立及基于混合校准方法的kalman滤波方法。本发明适合应用于车体晃动、人员上下、自然干扰、车辆发动机启动等车载环境下及拥有相同工作状态情况的其他载体的自对准过程。
背景技术
对于捷联惯导***(strapdown inertial navigation system)而言,初始对准是其开展导航任务的重要前提,其目的是建立起精确的初始坐标基准,即载体系相对于导航系的姿态矩阵。对准精度和对准速度是初始对准的两项重要技术指标。初始对准的精度对惯导***的导航精度有着重要影响,初始对准的速度在很大程度上决定了武器***的快速反应能力,因此要求对准精度高、对准时间短,既精又快。
初始对准按不同的划分标准可以把惯性导航***初始对准方式作如下分类:按对外部信息的依赖程度,可分为主动式对准和非主动式对准:按是否需要更高精度主惯导提供匹配参数,可分为传递对准和自对准;按对准的阶段,一般可分为粗对准和精对准;按基座的运动状态,可分为静基座对准和动基座对准。
常用的自对准方案是在忽略扰动的条件下利用加速度计和陀螺仪测量两个在空间不共线的矢量实现粗对准,求取粗略的姿态矩阵,然后利用基于小失准角线性误差模型的卡尔曼滤波算法完成初始姿态失准角的精确估计。但当载体处于晃动干扰条件下,惯性测量组件输出中包含干扰加速度信息和干扰角速度信息,鉴于陀螺输出信号的信噪比较低,而干扰信号频带较宽,此时难以获得准确的地球自转角速度,进而导致粗对准精度不高甚至得不到有效的对准结果,对精对准结果影响较大。
在精对准过程中,常用kalman滤波进行***状态变量估计,解算精确姿态转换矩阵,实现精对准过程。然而在实际***中,由于测量过程中动态干扰条件的影响,导致***模型的线性近似度变弱使得滤波效果变差,甚至可能直接导致滤波器发散。针对上述滤波环境,有必要提出新的滤波方法,能够修正***内部的状态误差,能保证长时间滤波的稳定性,实现高精度的精对准过程。
发明内容
本发明的目的在于,针对动基座条件下***噪声变大的情况,提出一种动基座下车载捷联惯导***自对准方法,首先采用凝固法粗对准,在时间可以精确已知的基础上,计算重力加速度矢量在惯性空间中旋转的角度,提取重力加速度矢量在惯性空间方向改变中所包含的真北信息,实现粗对准;然后在完成粗对准基础上,利用惯导***误差方程及惯性器件误差方程构建kalman滤波状态方程和量测方程;最后再利用混合校准kalman滤波方法估计惯导***速度误差、位置误差、姿态误差及惯性器件误差等并对其进行修正,实现高精度自对准。本发明具有较强的抗干扰能力,保证了动基座条件下的自对准精度。
为解决上述技术问题,本发明所采用的技术方案是:
一种动基座下车载捷联惯导***自对准方法,其特点是包括以下步骤:
步骤1,获得捷联惯导***所在地理位置信息,其中地理位置信息包括经度λ、纬度重力加速度g;
步骤2,利用经度、纬度、重力加速度信息及惯性器件输出信息,采用凝固坐标系,利用多矢量定姿方法进行粗对准,求得初始姿态矩阵;
步骤3,根据捷联惯性***误差方程及惯性器件误差方程构建kalman滤波状态方程和量测方程;
步骤4,设计混合校准kalman滤波方法估计***误差量和惯性器件误差量并修正,实现精对准过程,完成初始对准。
作为一种优选方式,所述步骤2中,初始姿态矩阵分解为以下公式(1)中的形式:
其中,
为地理坐标系e系与导航坐标系n系之间的变换矩阵,表示载体与地球的相对位置;
为地心惯性坐标系i系与地理坐标系e系之间的变换矩阵,反映了载体的自转,由时间间隔Δt=t-t0确定;
记录了载体坐标系b系与载体惯性凝固坐标系ib0系之间的姿态变化信息;
为载体惯性凝固坐标系ib0系与地心惯性坐标系i系之间的变换矩阵。
作为一种优选方式,
其中,为初始单位阵,为陀螺仪输出,表示载体坐标系b系相对于地心惯性坐标系i系的载体运动角速度在载体坐标系b系内的投影,为向量构造的叉乘反对称矩阵;
计算过程包括:
选择m个时刻的速度作为参考矢量,构造公式(5)中的目标函数:
其中,
公式(7)中,fb为加速度计输出;
整理变形公式(5)得:
L(C)=L'-tr(CBT) (8)
其中,且与B无关; 是使L取最小值的正交矩阵C;
对B进行奇异值分解,得:
B=USVT (9)
式(9)中,U和V为正交矩阵;S=diag(s1,s2,s3),s1≥s2≥s3;
求得:
作为一种优选方式,所述步骤3中,捷联惯性***误差方程及惯性器件误差方程如下:
其中,δve、δvn、δvu分别为惯导***东向、北向、天向速度误差;φe、φn、φu分别为俯仰、横滚、方位角误差;分别为X、Y、Z三轴向加速度计的零偏;τa为加速度计零偏相关时间;ε为陀螺的零漂;τg为陀螺零漂相关时间;ωie为地球自转角速度;RM为地球纬度圈半径,RN为地球经度圈半径;wa和wg分别为加速度计和陀螺仪零均值的白噪声;
根据捷联惯导***误差方程,建立***状态方程为:
其中,A(t)为状态转移矩阵,由式(11)得到;略去惯性器件的刻度系数误差和安装误差,又因为惯性器件偏置量的重复性误差对***精度的影响最大,因此将陀螺漂移和加速度零偏的随机常数部分列入状态,则滤波状态矢量及***噪声分别为:
假设载***置不变,且已经精确已知,选取位置误差作为观测量,则量测方程:
其中,分别为捷联惯导***解算的经纬度值;V为量测白噪声,H为量测矩阵,H=[I2×2 02×12]。
将式(12)和式(15)进行离散化,可得到***离散化的kalman滤波状态方程和量测方程。
作为一种优选方式,所述步骤4包括步骤41和步骤42,其中:
步骤41,设计混合校准kalman滤波方法:
输出校准采用公式(16)中的kalman滤波方程:
其中,表示一步预测估计值;表示Xk的状态估计值;Kk为k时刻的增益矩阵;Pk,k-1表示一步预测估计方差阵;Pk-1表示状态估计方差阵;Rk表示量测噪声方差阵;Qk-1表示***噪声方差阵;具有反馈控制的***状态方程和量测方程为:
对应的滤波方程为:
将控制矢量Uk-1选取为:
反馈校准的滤波方程为:
步骤42,精对准过程误差姿态误差修正:
姿态捷联矩阵修正如下:
其中,
kalman滤波估计出失准角φe、φn、φu后,通过混合校准方法对姿态矩阵进行校准,即可精确确定出载体姿态,完成自对准过程。
本发明具有以下有益效果:
1、引入惯性凝固思想:在时间可以精确已知的基础上,地球自转角速度可以精确计算。通过选取多个中间时刻计算重力加速度矢量在惯性空间中旋转的角度,提取重力加速度矢量在惯性空间方向改变中所包含的真北信息,实现粗对准。该方法只利用了地球自转角速度和重力加速度矢量在惯性空间的运动,避免了载体处于动基座条件下惯性器件测量受到的干扰信息的影响,有效提高了粗对准精度。
2、本发明提供混合校准kalman滤波方法,即在滤波初始阶段先进行输出校准,滤波稳定后进行反馈校准。进行反馈校准时将估计误差量反馈到***中,可以防止滤波结果发散,使滤波收敛速度更快,自对准精度更高。
3、反馈校准在滤波稳定后进行,更新后的状态协方差矩阵Pk反应了滤波估计的稳定状态,以前八个对角元素作为判断滤波稳定的参考量,通过与设定值比较,若小于设定值则认为此时滤波处于稳定状态可对***误差及惯性器件误差进行修正,修正后的误差参数直接进行下一周期计算,重复此过程可将误差逐步削减直至消失。
附图说明
图1为动基座下捷联惯导***自对准过程流程图。
图2为混合校准kalman滤波方法计算流程图。
具体实施方式
下面结合图1和图2对本发明的具体实施步骤进行详细的描述:
在下文实施步骤的的详细描述中坐标系定义如下:
地球坐标系e系:原点选取地球中心,X轴位于赤道平面内,从地心指向载体所在点经线,Z轴沿地球自转轴方向,随地球自转而转动,X轴、Y轴和Z轴构成右手坐标系,随地球自转而转动。
地心惯性坐标系i系:惯性坐标系是指在空间静止或匀速直线运动的参考坐标系。地心惯性坐标系是原点取在地球地心Oe时的右手直角惯性坐标系。zi轴沿地轴指向北极,yi轴,zi轴在地球赤道平面内,并指向空间的两颗恒星。
导航坐标系n系:是在导航时根据导航***工作的需要而选取的作为导航基准的坐标系。本发明选取东北天坐标系作为导航坐标系,即Ox指向水平东向(e),Oy轴指向水平北向(n),Oz沿垂线方向指向天(u)。
载体坐标系b系:坐标原点O位于载体的重心处,纵轴Oyb沿其首尾线方向并指向前,横轴Oxb指向其右,Ozb垂直于Oxbyb,沿载体竖轴指向上。
惯性凝固坐标系ib0系:在t0时刻将载体坐标系b经惯性凝固后得到,其中t0为粗对准的起始时刻。
动基座下车载捷联惯导***自对准方法包括以下步骤:
步骤1,获得捷联惯导***所在地理位置信息,其中地理位置信息包括经度λ、纬度重力加速度g。
步骤2,利用经度、纬度、重力加速度信息及惯性器件输出信息,采用凝固坐标系,利用多矢量定姿方法进行粗对准,求得初始姿态矩阵。
初始姿态矩阵分解为以下公式(1)中的形式:
其中,
为地理坐标系e系与导航坐标系n系之间的变换矩阵,表示载体与地球的相对位置;
为地心惯性坐标系i系与地理坐标系e系之间的变换矩阵,反映了载体的自转,由时间间隔Δt=t-t0确定;
记录了载体坐标系b系与载体惯性凝固坐标系ib0系之间的姿态变化信息;根据坐标系的相关定义,在载体纬度精确及地球自转角速度精确已知时,和可以精确求得;由定义可知,在t0时刻,可利用陀螺仪输出的b系相对于ib0系的角运动信息通过姿态跟踪算法实时求解,计算公式如下:
其中,为初始单位阵,为陀螺仪输出,表示载体坐标系b系相对于地心惯性坐标系i系的载体运动角速度在载体坐标系b系内的投影,为向量构造的叉乘反对称矩阵。
为载体惯性凝固坐标系ib0系与地心惯性坐标系i系之间的变换矩阵。由定义可知该矩阵不随时间变化且与载体的运动状态无关,为一常值矩阵。选择其中m个时刻的速度作为参考矢量,构造公式(5)中的目标函数:
其中,
公式(7)中,fb为加速度计输出;
整理变形公式(5)得:
L(C)=L'-tr(CBT) (8)
其中,且与B无关; 是使L取最小值的正交矩阵C(行列式为1)。
对B进行奇异值分解,得:
B=USVT (9)
式(9)中,U和V为正交矩阵;S=diag(s1,s2,s3),s1≥s2≥s3;
求得:
由式(6)计算t时刻的vi(t),由式(7)分别计算t时刻的再由式(10)即可计算通过式(2)~式(4)分别计算 代入式(1)即可得到粗略的初始姿态矩阵,完成粗对准。
步骤3,根据捷联惯性***误差方程及惯性器件误差方程构建kalman滤波状态方程和量测方程。其中,捷联惯性***误差方程及惯性器件误差方程如下:
其中,δve、δvn、δvu分别为惯导***东向、北向、天向速度误差;φe、φn、φu分别为俯仰、横滚、方位角误差;分别为X、Y、Z三轴向加速度计的零偏;τa为加速度计零偏相关时间;ε为陀螺的零漂;τg为陀螺零漂相关时间;ωie为地球自转角速度;RM为地球纬度圈半径,RN为地球经度圈半径;wa和wg分别为加速度计和陀螺仪零均值的白噪声。
根据捷联惯导***误差方程,建立***状态方程为:
其中,A(t)为状态转移矩阵,由式(11)得到;略去惯性器件的刻度系数误差和安装误差,又因为惯性器件偏置量的重复性误差对***精度的影响最大,因此将陀螺漂移和加速度零偏的随机常数部分列入状态,则滤波状态矢量及***噪声分别为:
假设载***置不变,且已经精确已知,选取位置误差作为观测量,则量测方程:
其中,分别为捷联惯导***解算的经纬度值;V为量测白噪声,H为量测矩阵,H=[I2×2 02×12]。
将式(12)和式(15)进行离散化,可得到***离散化的kalman滤波状态方程和量测方程。
步骤4,设计混合校准kalman滤波方法估计***误差量和惯性器件误差量并修正,实现精对准过程,完成初始对准。具体实现方法如下:
(1)设计混合校准kalman滤波方法
目前,应用kalman滤波方法对***进行校准常用的两种方法为:输出校准和反馈校准。进行输出校准时,产品进行滤波估计但是估计值不对***状态量进行修正;当应用反馈校准时,它将导航***输出导航参数的误差值的滤波估计反馈到内部计算中,直接参与每次的导航解算,在理想情况下,这些误差量会在不断的迭代过程中被渐进消除。但kalman滤波器从开始到稳定需要一定的时间,且在滤波刚开始的时候,估计误差很大,并且有研究显示滤波估计的精度会受到载体运动状态的影响,因而若是在此时采用反馈校准将估计误差较大的导航误差补偿到***内部计算中会更容易导致导航***精度的降低,因此一般情况下反馈校正不能在滤波不稳定时实时进行。基于对二者适用性的分析,本发明在滤波开始阶段使用输出校准方式进行滤波估计,滤波稳定后进行反馈校准即用滤波估计值对***速度误差、姿态误差、位置误差及惯性器件误差进行修正,修正后的***误差直接参与下一周期的计算,重复此过程直至自对准结束。
上述滤波稳定状态通过取状态估计方差阵Pk中对应姿态、速度、位置对角线元素作为度量滤波估值程度的参考量,这八个对角元素反映了相应估值的稳定程度,将其与设定值(根据实际***及实验条件进行设定)相比较,若小于设定值则认为滤波稳定,此时将滤波估计值对状态变量相应***误差及惯性器件误差进行修正,修正后误差参数重新参与下一周期的计算,否则只进行输出校准。
输出校准采用公式(16)中的一般kalman滤波方程即可实现:
其中,表示一步预测估计值;表示Xk的状态估计值;Kk为k时刻的增益矩阵;Pk,k-1表示一步预测估计方差阵;Pk-1表示状态估计方差阵;Rk表示量测噪声方差阵;Qk-1表示***噪声方差阵。
具有反馈控制的***状态方程和量测方程为:
与基本滤波方程相比,反馈校准仅在状态方程中增加控制项Bk-1Uk-1,其中Bk-1是控制矢量Uk-1的系数矩阵。其对应的滤波方程为:
根据反馈校准的原理可知,校准所要达到的效果就要将被估算出的误差量消减直至没有估计误差,因此可将控制矢量Uk-1选取为:
反馈校准的滤波方程为:
(2),精对准过程误差姿态误差修正
粗对准确定的姿态矩阵是精对准过程姿态更新的初值。由于粗对准确定的姿态矩阵误差较大,实际建立的导航坐标系n'与要求建立的理想导航坐标系n系间存在偏差角,该偏差角即为对准失准角。设该失准角为φe、φn、φu,粗对准后该偏差角都为小角,所以φ=[φe φn φu]T可看做是n系至n'系的等效旋转矢量。则姿态捷联矩阵修正如下:
其中,
由式(23)可知,kalman滤波估计出失准角φe、φn、φu后,通过混合校准方法对姿态矩阵进行校准,即可精确确定出载体姿态,完成自对准过程。
上面结合附图对本发明的实施例进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,而不是局限性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还可做出很多形式,这些均属于本发明的保护范围之内。
Claims (5)
1.一种动基座下车载捷联惯导***自对准方法,其特征在于,包括以下步骤:
步骤1,获得捷联惯导***所在地理位置信息,其中地理位置信息包括经度λ、纬度重力加速度g;
步骤2,利用经度、纬度、重力加速度信息及惯性器件输出信息,采用凝固坐标系,利用多矢量定姿方法进行粗对准,求得初始姿态矩阵;
步骤3,根据捷联惯性***误差方程及惯性器件误差方程构建kalman滤波状态方程和量测方程;
步骤4,设计混合校准kalman滤波方法估计***误差量和惯性器件误差量并修正,实现精对准过程,完成初始对准。
2.如权利要求1所述的动基座下车载捷联惯导***自对准方法,其特征在于,所述步骤2中,初始姿态矩阵分解为以下公式(1)中的形式:
其中,
为地理坐标系e系与导航坐标系n系之间的变换矩阵,表示载体与地球的相对位置;
为地心惯性坐标系i系与地理坐标系e系之间的变换矩阵,反映了载体的自转,由时间间隔Δt=t-t0确定;
记录了载体坐标系b系与载体惯性凝固坐标系ib0系之间的姿态变化信息;
为载体惯性凝固坐标系ib0系与地心惯性坐标系i系之间的变换矩阵。
3.如权利要求2所述的动基座下车载捷联惯导***自对准方法,其特征在于,
其中,为初始单位阵,为陀螺仪输出,表示载体坐标系b系相对于地心惯性坐标系i系的载体运动角速度在载体坐标系b系内的投影,为向量构造的叉乘反对称矩阵;
计算过程包括:
选择m个时刻的速度作为参考矢量,构造公式(5)中的目标函数:
其中,
公式(7)中,fb为加速度计输出;
整理变形公式(5)得:
L(C)=L'-tr(CBT) (8)
其中,且与B无关; 是使L取最小值的正交矩阵C;
对B进行奇异值分解,得:
B=USVT (9)
式(9)中,U和V为正交矩阵;S=diag(s1,s2,s3),s1≥s2≥s3;
求得:
4.如权利要求1或2所述的动基座下车载捷联惯导***自对准方法,其特征在于,所述步骤3中,捷联惯性***误差方程及惯性器件误差方程如下:
其中,δve、δvn、δvu分别为惯导***东向、北向、天向速度误差;φe、φn、φu分别为俯仰、横滚、方位角误差;分别为X、Y、Z三轴向加速度计的零偏;τa为加速度计零偏相关时间;ε为陀螺的零漂;τg为陀螺零漂相关时间;ωie为地球自转角速度;RM为地球纬度圈半径,RN为地球经度圈半径;wa和wg分别为加速度计和陀螺仪零均值的白噪声;
根据捷联惯导***误差方程,建立***状态方程为:
其中,A(t)为状态转移矩阵,由式(11)得到;略去惯性器件的刻度系数误差和安装误差,又因为惯性器件偏置量的重复性误差对***精度的影响最大,因此将陀螺漂移和加速度零偏的随机常数部分列入状态,则滤波状态矢量及***噪声分别为:
假设载***置不变,且已经精确已知,选取位置误差作为观测量,则量测方程:
其中,分别为捷联惯导***解算的经纬度值;V为量测白噪声,H为量测矩阵,H=[I2×2 02×12]。
将式(12)和式(15)进行离散化,可得到***离散化的kalman滤波状态方程和量测方程。
5.如权利要求1或2所述的动基座下车载捷联惯导***自对准方法,其特征在于,所述步骤4包括步骤41和步骤42,其中:
步骤41,设计混合校准kalman滤波方法:
输出校准采用公式(16)中的kalman滤波方程:
其中,表示一步预测估计值;表示Xk的状态估计值;Kk为k时刻的增益矩阵;Pk,k-1表示一步预测估计方差阵;Pk-1表示状态估计方差阵;Rk表示量测噪声方差阵;Qk-1表示***噪声方差阵;
具有反馈控制的***状态方程和量测方程为:
对应的滤波方程为:
将控制矢量Uk-1选取为:
反馈校准的滤波方程为:
步骤42,精对准过程误差姿态误差修正:
姿态捷联矩阵修正如下:
其中,
kalman滤波估计出失准角φe、φn、φu后,通过混合校准方法对姿态矩阵进行校准,即可精确确定出载体姿态,完成自对准过程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910766894.8A CN110440830B (zh) | 2019-08-20 | 2019-08-20 | 动基座下车载捷联惯导***自对准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910766894.8A CN110440830B (zh) | 2019-08-20 | 2019-08-20 | 动基座下车载捷联惯导***自对准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110440830A true CN110440830A (zh) | 2019-11-12 |
CN110440830B CN110440830B (zh) | 2023-03-28 |
Family
ID=68436614
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910766894.8A Active CN110440830B (zh) | 2019-08-20 | 2019-08-20 | 动基座下车载捷联惯导***自对准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110440830B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110906927A (zh) * | 2019-12-06 | 2020-03-24 | 中国空空导弹研究院 | 一种凝固坐标系下重力加速度简化算法 |
CN111795696A (zh) * | 2020-06-28 | 2020-10-20 | 中铁第一勘察设计院集团有限公司 | 一种基于零速修正的多惯导冗余***初始结构优选方法 |
CN112762961A (zh) * | 2020-12-28 | 2021-05-07 | 厦门华源嘉航科技有限公司 | 一种车载惯性里程计组合导航在线标定方法 |
CN112882118A (zh) * | 2020-12-30 | 2021-06-01 | 中国人民解放军海军工程大学 | 地固坐标系下动基座重力矢量估计方法、***及存储介质 |
CN113483756A (zh) * | 2021-07-13 | 2021-10-08 | 北京信息科技大学 | 数据的处理方法和***、存储介质及电子设备 |
CN113834499A (zh) * | 2021-08-26 | 2021-12-24 | 北京航天发射技术研究所 | 一种车载惯组与里程计行进间对准方法及*** |
CN116539029A (zh) * | 2023-04-03 | 2023-08-04 | 中山大学 | 一种水下动基座的基座定位方法、装置、存储介质及设备 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090089001A1 (en) * | 2007-08-14 | 2009-04-02 | American Gnc Corporation | Self-calibrated azimuth and attitude accuracy enhancing method and system (SAAAEMS) |
CN103245360A (zh) * | 2013-04-24 | 2013-08-14 | 北京工业大学 | 晃动基座下的舰载机旋转式捷联惯导***自对准方法 |
CN103557876A (zh) * | 2013-11-15 | 2014-02-05 | 山东理工大学 | 一种用于天线跟踪稳定平台的捷联惯导初始对准方法 |
CN103727940A (zh) * | 2014-01-15 | 2014-04-16 | 东南大学 | 基于重力加速度矢量匹配的非线性初始对准方法 |
FR3013829A1 (fr) * | 2013-11-22 | 2015-05-29 | Sagem Defense Securite | Procede d'alignement d'une centrale inertielle |
US20150285835A1 (en) * | 2013-12-23 | 2015-10-08 | InvenSense, Incorporated | Systems and methods for sensor calibration |
CN106123921A (zh) * | 2016-07-10 | 2016-11-16 | 北京工业大学 | 动态干扰条件下捷联惯导***的纬度未知自对准方法 |
CN106871928A (zh) * | 2017-01-18 | 2017-06-20 | 北京工业大学 | 基于李群滤波的捷联惯性导航初始对准方法 |
US20180274940A1 (en) * | 2015-10-13 | 2018-09-27 | Shanghai Huace Navigation Technology Ltd | A method for initial alignment of an inertial navigation apparatus |
CN108731702A (zh) * | 2018-07-03 | 2018-11-02 | 哈尔滨工业大学 | 一种基于Huber方法的大失准角传递对准方法 |
-
2019
- 2019-08-20 CN CN201910766894.8A patent/CN110440830B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090089001A1 (en) * | 2007-08-14 | 2009-04-02 | American Gnc Corporation | Self-calibrated azimuth and attitude accuracy enhancing method and system (SAAAEMS) |
CN103245360A (zh) * | 2013-04-24 | 2013-08-14 | 北京工业大学 | 晃动基座下的舰载机旋转式捷联惯导***自对准方法 |
CN103557876A (zh) * | 2013-11-15 | 2014-02-05 | 山东理工大学 | 一种用于天线跟踪稳定平台的捷联惯导初始对准方法 |
FR3013829A1 (fr) * | 2013-11-22 | 2015-05-29 | Sagem Defense Securite | Procede d'alignement d'une centrale inertielle |
US20150285835A1 (en) * | 2013-12-23 | 2015-10-08 | InvenSense, Incorporated | Systems and methods for sensor calibration |
CN103727940A (zh) * | 2014-01-15 | 2014-04-16 | 东南大学 | 基于重力加速度矢量匹配的非线性初始对准方法 |
US20180274940A1 (en) * | 2015-10-13 | 2018-09-27 | Shanghai Huace Navigation Technology Ltd | A method for initial alignment of an inertial navigation apparatus |
CN106123921A (zh) * | 2016-07-10 | 2016-11-16 | 北京工业大学 | 动态干扰条件下捷联惯导***的纬度未知自对准方法 |
CN106871928A (zh) * | 2017-01-18 | 2017-06-20 | 北京工业大学 | 基于李群滤波的捷联惯性导航初始对准方法 |
CN108731702A (zh) * | 2018-07-03 | 2018-11-02 | 哈尔滨工业大学 | 一种基于Huber方法的大失准角传递对准方法 |
Non-Patent Citations (3)
Title |
---|
GAO WEI 等: "The application in SINS"s initial alignment based on combination of LQ optimal control and Kalman filter", 《2008 IEEE INTERNATIONAL CONFERENCE ON MECHATRONICS AND AUTOMATION》, 6 March 2009 (2009-03-06), pages 1 - 5 * |
赵小明等: "车载单轴旋转激光捷联惯导抗晃动初始对准和零速修正方法", 《中国惯性技术学报》, no. 03, 15 June 2013 (2013-06-15), pages 302 - 307 * |
黄湘远: "5阶CKF在捷联惯导非线性对准中的应用研究", 《***工程与电子技术》, vol. 37, no. 3, 31 March 2015 (2015-03-31), pages 633 - 638 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110906927A (zh) * | 2019-12-06 | 2020-03-24 | 中国空空导弹研究院 | 一种凝固坐标系下重力加速度简化算法 |
CN110906927B (zh) * | 2019-12-06 | 2023-04-14 | 中国空空导弹研究院 | 一种凝固坐标系下重力加速度简化算法 |
CN111795696A (zh) * | 2020-06-28 | 2020-10-20 | 中铁第一勘察设计院集团有限公司 | 一种基于零速修正的多惯导冗余***初始结构优选方法 |
CN112762961A (zh) * | 2020-12-28 | 2021-05-07 | 厦门华源嘉航科技有限公司 | 一种车载惯性里程计组合导航在线标定方法 |
CN112762961B (zh) * | 2020-12-28 | 2022-06-03 | 厦门华源嘉航科技有限公司 | 一种车载惯性里程计组合导航在线标定方法 |
CN112882118A (zh) * | 2020-12-30 | 2021-06-01 | 中国人民解放军海军工程大学 | 地固坐标系下动基座重力矢量估计方法、***及存储介质 |
CN113483756A (zh) * | 2021-07-13 | 2021-10-08 | 北京信息科技大学 | 数据的处理方法和***、存储介质及电子设备 |
CN113834499A (zh) * | 2021-08-26 | 2021-12-24 | 北京航天发射技术研究所 | 一种车载惯组与里程计行进间对准方法及*** |
CN116539029A (zh) * | 2023-04-03 | 2023-08-04 | 中山大学 | 一种水下动基座的基座定位方法、装置、存储介质及设备 |
CN116539029B (zh) * | 2023-04-03 | 2024-02-23 | 中山大学 | 一种水下动基座的基座定位方法、装置、存储介质及设备 |
Also Published As
Publication number | Publication date |
---|---|
CN110440830B (zh) | 2023-03-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Bian et al. | Inertial navigation | |
CN109556632B (zh) | 一种基于卡尔曼滤波的ins/gnss/偏振/地磁组合导航对准方法 | |
CN110440830A (zh) | 动基座下车载捷联惯导***自对准方法 | |
US9541392B2 (en) | Surveying system and method | |
Fang et al. | Predictive iterated Kalman filter for INS/GPS integration and its application to SAR motion compensation | |
Fakharian et al. | Adaptive Kalman filtering based navigation: An IMU/GPS integration approach | |
De Agostino et al. | Performances comparison of different MEMS-based IMUs | |
CN112505737B (zh) | 一种gnss/ins组合导航方法 | |
Stančić et al. | The integration of strap-down INS and GPS based on adaptive error damping | |
JP2007536499A (ja) | 測定及び記憶された重力傾度を用いて慣性航法測定値の精度を改善する方法及びシステム | |
Oh | Multisensor fusion for autonomous UAV navigation based on the Unscented Kalman Filter with Sequential Measurement Updates | |
Xue et al. | In-motion alignment algorithm for vehicle carried SINS based on odometer aiding | |
CN109708663B (zh) | 基于空天飞机sins辅助的星敏感器在线标定方法 | |
CN110849360A (zh) | 面向多机协同编队飞行的分布式相对导航方法 | |
CN114877915B (zh) | 一种激光陀螺惯性测量组件g敏感性误差标定装置及方法 | |
CN111220151B (zh) | 载体系下考虑温度模型的惯性和里程计组合导航方法 | |
Zorina et al. | Enhancement of INS/GNSS integration capabilities for aviation-related applications | |
Omerbashich | Integrated INS/GPS navigation from a popular perspective | |
Mu et al. | Improved decentralized GNSS/SINS/odometer fusion system for land vehicle navigation applications | |
CN112683265B (zh) | 一种基于快速iss集员滤波的mimu/gps组合导航方法 | |
Emel’yantsev et al. | Use of maneuvering to improve the accuracy of ship autonomous SINS | |
Akimov et al. | Compensation for errors in determining the angle in the flight-navigation complex of the aircraft in case of failure of the satellite navigation system | |
CN111473786A (zh) | 一种基于局部反馈的两层分布式多传感器组合导航滤波方法 | |
Ebrahim et al. | Initial alignment of strap-down inertial navigation system on stationary base for high-speed flying vehicle | |
Gothard et al. | Lessons learned on a low-cost global navigation system for the Surrogate Semiautonomous Vehicle (SSV) |
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 |