CN111238532B - 一种适用于晃动基座环境的惯性测量单元标定方法 - Google Patents
一种适用于晃动基座环境的惯性测量单元标定方法 Download PDFInfo
- Publication number
- CN111238532B CN111238532B CN201911340104.6A CN201911340104A CN111238532B CN 111238532 B CN111238532 B CN 111238532B CN 201911340104 A CN201911340104 A CN 201911340104A CN 111238532 B CN111238532 B CN 111238532B
- Authority
- CN
- China
- Prior art keywords
- measurement unit
- order
- axis
- parameter
- accelerometer
- 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
Images
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
本发明公开了一种适用于晃动基座环境的惯性测量单元标定方法,涉及航空航天捷联惯性导航技术中的惯性测量组合测试技术领域,该方法包括利用一阶中间参数和一阶误差参数间的关系建立联立方程,基于所建立的联立方程,求解得到一阶误差参数;利用一阶误差参数计算得到第i个位置的二阶中间参数,根据二阶中间参数和二阶误差参数间的关系建立方程,求解得到二阶误差参数,所述二阶误差参数为陀螺零偏;判断一阶误差参数和二阶误差参数的残差是否小于设定阈值,若是,则结束。本发明采用迭代方法逐步消除在初始对准过程中晃动基座环境引入的姿态误差,从而更好的适用于晃动基座环境下的标定。
Description
技术领域
本发明涉及航空航天捷联惯性导航技术中的惯性测量组合测试技术领域,具体涉及一种适用于晃动基座环境的惯性测量单元标定方法。
背景技术
惯性导航***能够提供完全自主的导航信息,具有反应时间短、可靠性高、体积小、重量轻等优点,广泛应用于飞机、舰船等军用和民用导航领域,具有重要的国防意义和巨大的经济效益。但惯性导航***在使用的过程中,其误差会随着时间不断积累,当前,两条措施可以提高惯导***的导航精度:一条是设计和制造出更高精度的惯性元件,如采用新型材料、新工艺或新型惯性器件;另一条就是采用***技术,如运用先进的控制策略与算法。但当前光学陀螺的研制已经达到较高的水平,惯性器件精度的提升空间已然有限,为了进一步提升光学陀螺惯导***的精度,国内外研究人员开始转向更合理、有效的***补偿技术的研究。
标定技术是惯性导航领域的核心技术之一,其是一种误差辨识技术,即首先建立惯性测量单元误差模型,然后通过一系列试验求解出误差模型中的误差项,进而通过软件算法对误差进行补偿。惯性测量单元的标定结果好坏直接影响到惯性导航***的精度。
当前广泛使用的惯性测量单元标定方法需要在静基座环境下应用,即在整个标定过程中,惯性测量单元在除翻转外,其它时间均处于静止状态,但在实际工况下,难以避免因发动机振动、人员走动或自然环境(强风、浪涌等)引起的载体晃动,因而极大限制了标定技术的应用。
发明内容
针对现有技术中存在的缺陷,本发明的目的在于提供一种适用于晃动基座环境的惯性测量单元标定方法,采用迭代方法逐步消除在初始对准过程中晃动基座环境引入的姿态误差,从而更好的适用于晃动基座环境下的标定。
为达到以上目的,本发明采取的技术方案是,包括:
沿惯性测量单元的轴向方向依次转动设定角度,且相邻两次转动过程间,惯性测量单元静止设定时间,并采集整个转动过程中惯性测量单元输出的数据;
基于采集的数据,在惯性测量单元的相邻两次转动过程中,以及第i个位置上的静止过程中进行连续导航,获取转动到达第i个位置时的瞬间速度和瞬间天向转角,以及转动完成后在第i个位置静止过程中的实时速度和实时天向转角;
基于获取的瞬间速度、瞬间天向转角、实时速度和实时天向转角,进行拟合,得到第i个位置的一阶中间参数;
利用一阶中间参数和一阶误差参数间的关系建立联立方程,基于所建立的联立方程,求解得到一阶误差参数,所述一阶误差参数包括加速度计零偏、加速度计标度因数、加速度计失准角、加速度计二次项系数、陀螺标度因数和陀螺失准角;
利用一阶误差参数计算得到第i个位置的二阶中间参数,根据二阶中间参数和二阶误差参数间的关系建立方程,求解得到二阶误差参数,所述二阶误差参数为陀螺零偏;
判断一阶误差参数和二阶误差参数的残差是否小于设定阈值,若是,则结束,若否,则使用得到的一阶误差参数和二阶误差参数的残差对一阶误差参数和二阶误差参数进行补偿,再次计算得到一阶误差参数和二阶误差参数,以此类推,直至计算得到的一阶误差参数和二阶误差参数的残差小于设定阈值。
在上述技术方案的基础上,所述沿惯性测量单元的轴向方向依次转动设定角度,且相邻两次转动过程间,惯性测量单元静止设定时间,具体步骤包括:
将惯性测量单元安装于双轴转动设备上,从惯性测量单元的初始位置开始,沿惯性测量单元的轴向方向依次转动设定角度,且每转动一次设定角度后,到达一个位置,然后惯性测量单元在当前位置静止设定时间,再转动设定角度,达到下一个位置,所述i为不小于0的整数,当i取值为0时,表示惯性测量单元处于初始位置。
在上述技术方案的基础上,所述惯性测量单元的惯性测量单元坐标系具体为:x轴与x加速度计输入轴方向相同,y轴位于x加速度计和y加速度计输入轴构成的平面内,并接近y加速度计输入轴方向,z轴方向由右手定则确定。
在上述技术方案的基础上,所述基于采集的数据,在惯性测量单元的相邻两次转动过程中,以及第i个位置上的静止过程中进行连续导航,获取转动到达第i个位置时的瞬间速度和瞬间天向转角,以及转动完成后在第i个位置静止过程中的实时速度和实时天向转角,具体步骤包括:
当i等于0时,基于采集的数据,使用双矢量定姿法进行解析粗对准,然后进行精对准,得到初始位置上导航起始时刻的瞬间速度和瞬间天向转角,然后基于精对准结果和初始位置上惯性测量单元输出的数据,进行导航解算,得到初始位置上导航过程中的实时速度和实时天向转角;
当i大于0时,在第i个位置上使用上一次转动过程中的惯性测量单元输出的数据进行粗对准,然后进行精对准,然后在第i-1到第i个位置的转动过程中,以及在第i个位置上的静止过程中进行连续导航,获取转动到达第i个位置时的瞬间速度和瞬间天向转角,以及转动完成后在第i个位置静止过程中的实时速度和实时天向转角。
在上述技术方案的基础上,所述基于获取的瞬间速度、瞬间天向转角、实时速度和实时天向转角,进行拟合,得到第i个位置的一阶中间参数,进行拟合的公式为:
其中,g表示重力加速度,T表示实时时间,表示到达第i个位置时在x轴方向上的瞬间速度,表示到达第i个位置时在y轴方向上的瞬间速度,表示到达第i个位置时在z轴方向上的瞬间速度,表示转动到达第i个位置时在x轴方向上的瞬间速度,表示转动到达第i个位置时在y轴方向上的瞬间速度,表示转动到达第i个位置时在z轴方向上的瞬间速度,θn(i)在第i个位置静止过程中的实时天向转角,转动到达第i个位置时的瞬间天向转角,ωv表示降低系数,ωv包括ωvx、ωvy和ωvz,ωvx表示ωv在x轴方向上的分量,ωvy表示ωv在y轴方向上的分量,ωvz表示ωv在z轴方向上的分量,Δg表示一阶中间参数,Δg包括Δgx、Δgy和Δgz,Δgx表示Δg在x轴方向上投影的标量,Δgy表示Δg在y轴方向上投影的标量,Δgz表示Δg在z轴方向上投影的标量。
在上述技术方案的基础上,
在所述惯性测量单元坐标系中,所述加速度计的误差模型为:
所述误差模型的向量表示形式为:
其中,fb表示载体坐标系下加速度计测得的比力,表示载体坐标系下的加速度计零偏,Ka表示加速度计标度因数误差和加速度计失准角,Ka2表示加速度计二次项系数,δfb表示载体坐标系下加速度计测得的比力误差,fx、fy和fz分别表示载体坐标系下加速度计测得的x轴、y轴和z轴方向上的比力,δfx、δfy和δfz分别表示载体坐标系下加速度计测得的x轴、y轴和z轴方向上的比力误差,Bax、Bay和Baz表示加速度计零偏,Kaxx、Kayy和Kazz表示加速度计标度因数,Kax2、Kay2、Kaz2表示加速度计二次项系数。
在上述技术方案的基础上,
所述惯性测量单元的陀螺的误差模型为:
所述陀螺的误差模型的向量表示形式为:
其中,εb表示载体坐标系下陀螺测得的角速度误差,表示载体坐标系下陀螺零偏,Kg表示陀螺标度因数误差和陀螺失准角,ωb表示载体坐标系下陀螺测得的角速度,εx、εy和εz分别表示载体坐标系下陀螺测得的x轴、y轴和z轴方向上的角速度误差,Bgx、Bgy和Bgz分别表示载体坐标系下x轴、y轴和z轴方向上的陀螺零偏,Kgxx、Kgyy和Kgzz表示陀螺标度因数,Kgxy、Kgxz、Kgyz、Kgyx、Kgzx和Kgzy表示陀螺失准角,ωx、ωy和ωz分别表示载体坐标系下陀螺测得的x轴、y轴和z轴方向上的角速度。
在上述技术方案的基础上,所述利用一阶中间参数和一阶误差参数间的关系建立联立方程,所述联立方程为Δg=AK,其中,A表示系数矩阵,它反映了一阶中间参数Δg与一阶误差参数K的系数关系,K表示一阶误差参数。
在上述技术方案的基础上,所述利用一阶误差参数计算得到第i个位置的二阶中间参数,具体过程包括:
在上述技术方案的基础上,所述根据二阶中间参数和二阶误差参数间的关系建立方程,求解得到二阶误差参数,其中,所述方程为:
与现有技术相比,本发明的优点在于:通过计算一阶中间参数和一阶误差参数,进而计算得到二阶中间参数和二阶误差参数,对晃动基座环境下引入的额外误差进行了建模,采用迭代方法逐步消除在初始对准过程中晃动基座环境引入的姿态误差,从而更好的适用于晃动基座环境下的标定。
附图说明
图1为本发明实施例中一种适用于晃动基座环境的惯性测量单元标定方法的流程图。
具体实施方式
以下结合附图及实施例对本发明作进一步详细说明。
参见图1所示,本发明实施例提供的一种适用于晃动基座环境的惯性测量单元标定方法,包括以下步骤:
S1:沿惯性测量单元的轴向方向依次转动设定角度,且相邻两次转动过程间,惯性测量单元静止设定时间,并采集整个转动过程中惯性测量单元输出的数据。
本发明实施例中,沿惯性测量单元的轴向方向依次转动设定角度,且相邻两次转动过程间,惯性测量单元静止设定时间,具体步骤包括:
将惯性测量单元安装于双轴转动设备上,从惯性测量单元的初始位置开始,沿惯性测量单元的轴向方向依次转动设定角度,且每转动一次设定角度后,到达一个位置,然后惯性测量单元在当前位置静止设定时间,再转动设定角度,达到下一个位置,所述i为不小于0的整数,当i取值为0时,表示惯性测量单元处于初始位置。
惯性测量单元的初始位置保证x轴朝上或朝下,优选的,导航坐标系选取北天东坐标系,惯性测量单元通电预热30分钟开始采集输出的数据,采样周期为0.01s,同时惯性测量单元开始旋转,惯性测量单元的旋转顺序为:沿z轴顺时针旋转90度、沿z轴顺时针旋转90度、沿z轴顺时针旋转90度、沿z轴逆时针旋转90度、沿z轴逆时针旋转90度、沿z轴逆时针旋转90度、沿x轴顺时针旋转90度、沿y轴顺时针旋转90度、沿y轴顺时针旋转90度、沿y轴顺时针旋转90度、沿y轴逆时针旋转90度、沿y轴逆时针旋转90度、沿y轴逆时针旋转90度、沿x轴顺时针旋转90度、沿x轴顺时针旋转90度、沿x轴逆时针旋转90度、沿x轴逆时针旋转90度、沿x轴逆时针旋转90度。惯性测量单元每旋转一次到达位置后,静止3~5分钟,再开始下一次的旋转,旋转全部结束后,停止采集惯性测量单元输出的数据,并关闭惯性测量单元。
本发明实施例中,惯性测量单元的惯性测量单元坐标系具体为:x轴与x加速度计输入轴方向相同,y轴位于x加速度计和y加速度计输入轴构成的平面内,并接近y加速度计输入轴方向,z轴方向由右手定则确定。
S2:基于采集的数据,在惯性测量单元的相邻两次转动过程中,以及第i个位置上的静止过程中进行连续导航,获取转动到达第i个位置时的瞬间速度和瞬间天向转角,以及转动完成后在第i个位置静止过程中的实时速度和实时天向转角。
本发明实施例中,基于采集的数据,在惯性测量单元的相邻两次转动过程中,以及第i个位置上的静止过程中进行连续导航,获取转动到达第i个位置时的瞬间速度和瞬间天向转角,以及转动完成后在第i个位置静止过程中的实时速度和实时天向转角,具体步骤包括:
当i等于0时,基于采集的数据,使用双矢量定姿法进行解析粗对准,然后进行精对准,得到初始位置上导航起始时刻的瞬间速度和瞬间天向转角,然后基于精对准结果和初始位置上惯性测量单元输出的数据,进行导航解算,得到初始位置上导航过程中的实时速度和实时天向转角;
当i大于0时,在第i个位置上使用上一次转动过程中的惯性测量单元输出的数据进行粗对准,然后进行精对准,然后在第i-1到第i个位置的转动过程中,以及在第i个位置上的静止过程中进行连续导航,获取转动到达第i个位置时的瞬间速度和瞬间天向转角,以及转动完成后在第i个位置静止过程中的实时速度和实时天向转角。
本发明实施例中,使用的精对准算法为:
(1)参数修正:
(2)使用修正后的姿态与地速,在每个标定位置的静止段进行导航与姿态更新;
(3)更新迭代中间参数说明
(4)迭代中间参数解算方法
使用最小二乘拟合方法计算第j个标定位置的中间参数,其中F1为二次多项式拟合方法的广义逆运算。脚标C1表示矩阵的第1列。
其中脚标1、3分别表示零次项、一次项系数、二次项系数,那么迭代中间参数更新如下:
S3:基于获取的瞬间速度、瞬间天向转角、实时速度和实时天向转角,进行拟合,得到第i个位置的一阶中间参数;
本发明实施例中,基于获取的瞬间速度、瞬间天向转角、实时速度和实时天向转角,进行拟合,得到第i个位置的一阶中间参数,进行拟合的公式为:
其中,g表示重力加速度,T表示实时时间,表示到达第i个位置时在x轴方向上的瞬间速度,表示到达第i个位置时在y轴方向上的瞬间速度,表示到达第i个位置时在z轴方向上的瞬间速度,表示转动到达第i个位置时在x轴方向上的瞬间速度,表示转动到达第i个位置时在y轴方向上的瞬间速度,表示转动到达第i个位置时在z轴方向上的瞬间速度,θn(i)在第i个位置静止过程中的实时天向转角,转动到达第i个位置时的瞬间天向转角,ωv表示降低系数,ωv包括ωvx、ωvy和ωvz,ωvx表示ωv在x轴方向上的分量,ωvy表示ωv在y轴方向上的分量,ωvz表示ωv在z轴方向上的分量,Δg表示一阶中间参数,Δg包括Δgx、Δgy和Δgz,Δgx表示Δg在x轴方向上投影的标量,Δgy表示Δg在y轴方向上投影的标量,Δgz表示Δg在z轴方向上投影的标量。
在所述惯性测量单元坐标系中,所述加速度计的误差模型为:
所述误差模型的向量表示形式为:
其中,fb表示载体坐标系下加速度计测得的比力,表示载体坐标系下的加速度计零偏,Ka表示加速度计标度因数误差和加速度计失准角,Ka2表示加速度计二次项系数,δfb表示载体坐标系下加速度计测得的比力误差,fx、fy和fz分别表示载体坐标系下加速度计测得的x轴、y轴和z轴方向上的比力,δfx、δfy和δfz分别表示载体坐标系下加速度计测得的x轴、y轴和z轴方向上的比力误差,Bax、Bay和Baz表示加速度计零偏,Kaxx、Kayy和Kazz表示加速度计标度因数,Kax2、Kay2、Kaz2表示加速度计二次项系数。
惯性测量单元的陀螺的误差模型为:
所述陀螺的误差模型的向量表示形式为:
其中,εb表示载体坐标系下陀螺测得的角速度误差,表示载体坐标系下陀螺零偏,Kg表示陀螺标度因数误差和陀螺失准角,ωb表示载体坐标系下陀螺测得的角速度,εx、εy和εz分别表示载体坐标系下陀螺测得的x轴、y轴和z轴方向上的角速度误差,Bgx、Bgy和Bgz分别表示载体坐标系下x轴、y轴和z轴方向上的陀螺零偏,Kgxx、Kgyy和Kgzz表示陀螺标度因数,Kgxy、Kgxz、Kgyz、Kgyx、Kgzx和Kgzy表示陀螺失准角,ωx、ωy和ωz分别表示载体坐标系下陀螺测得的x轴、y轴和z轴方向上的角速度。
S4:利用一阶中间参数和一阶误差参数间的关系建立联立方程,基于所建立的联立方程,求解得到一阶误差参数,一阶误差参数包括加速度计零偏、加速度计标度因数、加速度计失准角、加速度计二次项系数、陀螺标度因数和陀螺失准角。
本发明实施例中,利用一阶中间参数和一阶误差参数间的关系建立联立方程,所述联立方程为Δg=AK,其中,A表示系数矩阵,它反映了一阶中间参数Δg与一阶误差参数K的系数关系,K表示一阶误差参数。
S5:利用一阶误差参数计算得到第i个位置的二阶中间参数,根据二阶中间参数和二阶误差参数间的关系建立方程,求解得到二阶误差参数,所述二阶误差参数为陀螺零偏。
本发明实施例中,利用一阶误差参数计算得到第i个位置的二阶中间参数,具体过程包括:
本发明实施例中,根据二阶中间参数和二阶误差参数间的关系建立方程,求解得到二阶误差参数,其中,所述方程为:
S6:判断一阶误差参数和二阶误差参数的残差是否小于设定阈值,若是,则结束,若否,则使用得到的一阶误差参数和二阶误差参数的残差对一阶误差参数和二阶误差参数进行补偿,再次计算得到一阶误差参数和二阶误差参数,以此类推,直至计算得到的一阶误差参数和二阶误差参数的残差小于设定阈值。
晃动基座条件主要影响初始对准过程,由于方位误差已放在二阶误差参数中进行估计,为提高初始对准精度,在各个静止位置上需要估计初始水平姿态误差,并使用迭代方法进一步消除随机误差的影响。
本发明实施例的适用于晃动基座环境的惯性测量单元标定方法,通过计算一阶中间参数和一阶误差参数,进而计算得到二阶中间参数和二阶误差参数,对晃动基座环境下引入的额外误差进行了建模,采用迭代方法逐步消除在初始对准过程中晃动基座环境引入的姿态误差,从而更好的适用于晃动基座环境下的标定。
本发明不局限于上述实施方式,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也视为本发明的保护范围之内。本说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
Claims (9)
1.一种适用于晃动基座环境的惯性测量单元标定方法,其特征在于,包括:
沿惯性测量单元的轴向方向依次转动设定角度,且相邻两次转动过程间,惯性测量单元静止设定时间,并采集整个转动过程中惯性测量单元输出的数据;
基于采集的数据,在惯性测量单元的相邻两次转动过程中,以及第i个位置上的静止过程中进行连续导航,获取转动到达第i个位置时的瞬间速度和瞬间天向转角,以及转动完成后在第i个位置静止过程中的实时速度和实时天向转角;
基于获取的瞬间速度、瞬间天向转角、实时速度和实时天向转角,进行拟合,得到第i个位置的一阶中间参数;
利用一阶中间参数和一阶误差参数间的关系建立联立方程,基于所建立的联立方程,求解得到一阶误差参数,所述一阶误差参数包括加速度计零偏、加速度计标度因数、加速度计失准角、加速度计二次项系数、陀螺标度因数和陀螺失准角;
利用一阶误差参数计算得到第i个位置的二阶中间参数,根据二阶中间参数和二阶误差参数间的关系建立方程,求解得到二阶误差参数,所述二阶误差参数为陀螺零偏;
判断一阶误差参数和二阶误差参数的残差是否小于设定阈值,若是,则结束,若否,则使用得到的一阶误差参数和二阶误差参数的残差对一阶误差参数和二阶误差参数进行补偿,再次计算得到一阶误差参数和二阶误差参数,以此类推,直至计算得到的一阶误差参数和二阶误差参数的残差小于设定阈值;
其中,所述沿惯性测量单元的轴向方向依次转动设定角度,且相邻两次转动过程间,惯性测量单元静止设定时间,具体步骤包括:
将惯性测量单元安装于双轴转动设备上,从惯性测量单元的初始位置开始,沿惯性测量单元的轴向方向依次转动设定角度,且每转动一次设定角度后,到达一个位置,然后惯性测量单元在当前位置静止设定时间,再转动设定角度,达到下一个位置,所述i为不小于0的整数,当i取值为0时,表示惯性测量单元处于初始位置。
2.如权利要求1所述的一种适用于晃动基座环境的惯性测量单元标定方法,其特征在于,所述惯性测量单元的惯性测量单元坐标系具体为:x轴与x加速度计输入轴方向相同,y轴位于x加速度计和y加速度计输入轴构成的平面内,并接近y加速度计输入轴方向,z轴方向由右手定则确定。
3.如权利要求2所述的一种适用于晃动基座环境的惯性测量单元标定方法,其特征在于,所述基于采集的数据,在惯性测量单元的相邻两次转动过程中,以及第i个位置上的静止过程中进行连续导航,获取转动到达第i个位置时的瞬间速度和瞬间天向转角,以及转动完成后在第i个位置静止过程中的实时速度和实时天向转角,具体步骤包括:
当i等于0时,基于采集的数据,使用双矢量定姿法进行解析粗对准,然后进行精对准,得到初始位置上导航起始时刻的瞬间速度和瞬间天向转角,然后基于精对准结果和初始位置上惯性测量单元输出的数据,进行导航解算,得到初始位置上导航过程中的实时速度和实时天向转角;
当i大于0时,在第i个位置上使用上一次转动过程中的惯性测量单元输出的数据进行粗对准,然后进行精对准,然后在第i-1到第i个位置的转动过程中,以及在第i个位置上的静止过程中进行连续导航,获取转动到达第i个位置时的瞬间速度和瞬间天向转角,以及转动完成后在第i个位置静止过程中的实时速度和实时天向转角。
4.如权利要求3所述的一种适用于晃动基座环境的惯性测量单元标定方法,其特征在于,所述基于获取的瞬间速度、瞬间天向转角、实时速度和实时天向转角,进行拟合,得到第i个位置的一阶中间参数,进行拟合的公式为:
其中,g表示重力加速度,T表示实时时间,表示到达第i个位置时在x轴方向上的瞬间速度,表示到达第i个位置时在y轴方向上的瞬间速度,表示到达第i个位置时在z轴方向上的瞬间速度,表示转动到达第i个位置时在x轴方向上的瞬间速度,表示转动到达第i个位置时在y轴方向上的瞬间速度,表示转动到达第i个位置时在z轴方向上的瞬间速度,θn(i)在第i个位置静止过程中的实时天向转角,转动到达第i个位置时的瞬间天向转角,ωv表示降低系数,ωv包括ωvx、ωvy和ωvz,ωvx表示ωv在x轴方向上的分量,ωvy表示ωv在y轴方向上的分量,ωvz表示ωv在z轴方向上的分量,Δg表示一阶中间参数,Δg包括Δgx、Δgy和Δgz,Δgx表示Δg在x轴方向上投影的标量,Δgy表示Δg在y轴方向上投影的标量,Δgz表示Δg在z轴方向上投影的标量。
5.如权利要求4所述的一种适用于晃动基座环境的惯性测量单元标定方法,其特征在于,
在所述惯性测量单元坐标系中,所述加速度计的误差模型为:
所述误差模型的向量表示形式为:
6.如权利要求5所述的一种适用于晃动基座环境的惯性测量单元标定方法,其特征在于,
所述惯性测量单元的陀螺的误差模型为:
所述陀螺的误差模型的向量表示形式为:
7.如权利要求6所述的一种适用于晃动基座环境的惯性测量单元标定方法,其特征在于,所述利用一阶中间参数和一阶误差参数间的关系建立联立方程,所述联立方程为Δg=AK,其中,A表示系数矩阵,它反映了一阶中间参数Δg与一阶误差参数K的系数关系,K表示一阶误差参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911340104.6A CN111238532B (zh) | 2019-12-23 | 2019-12-23 | 一种适用于晃动基座环境的惯性测量单元标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911340104.6A CN111238532B (zh) | 2019-12-23 | 2019-12-23 | 一种适用于晃动基座环境的惯性测量单元标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111238532A CN111238532A (zh) | 2020-06-05 |
CN111238532B true CN111238532B (zh) | 2022-02-01 |
Family
ID=70863892
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911340104.6A Active CN111238532B (zh) | 2019-12-23 | 2019-12-23 | 一种适用于晃动基座环境的惯性测量单元标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111238532B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112033438B (zh) * | 2020-08-18 | 2022-09-02 | 湖北航天技术研究院总体设计所 | 一种基于速度拟合的晃动基座自对准方法 |
CN113566849B (zh) * | 2021-07-29 | 2024-03-05 | 深圳元戎启行科技有限公司 | 惯性测量单元的安装角度标定方法、装置和计算机设备 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101131311A (zh) * | 2007-10-15 | 2008-02-27 | 北京航空航天大学 | 一种智能化机载导弹动基座对准及标定方法 |
CN103604442A (zh) * | 2013-11-14 | 2014-02-26 | 哈尔滨工程大学 | 应用于捷联惯导***在线标定的可观测性分析方法 |
CN103983274A (zh) * | 2014-04-11 | 2014-08-13 | 湖北航天技术研究院总体设计所 | 一种适用于低精度无方位基准双轴转位设备的惯性测量单元标定方法 |
CN103994775A (zh) * | 2014-04-11 | 2014-08-20 | 湖北航天技术研究院总体设计所 | 一种适用于低精度有方位基准双轴转位设备的惯性测量单元标定方法 |
CN104121928A (zh) * | 2014-05-29 | 2014-10-29 | 湖北航天技术研究院总体设计所 | 一种适用于低精度有方位基准单轴转位设备的惯性测量单元标定方法 |
CN108458725A (zh) * | 2017-11-17 | 2018-08-28 | 北京计算机技术及应用研究所 | 捷联惯导***晃动基座上的***级标定方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105203129B (zh) * | 2015-10-13 | 2019-05-07 | 上海华测导航技术股份有限公司 | 一种惯导装置初始对准方法 |
-
2019
- 2019-12-23 CN CN201911340104.6A patent/CN111238532B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101131311A (zh) * | 2007-10-15 | 2008-02-27 | 北京航空航天大学 | 一种智能化机载导弹动基座对准及标定方法 |
CN103604442A (zh) * | 2013-11-14 | 2014-02-26 | 哈尔滨工程大学 | 应用于捷联惯导***在线标定的可观测性分析方法 |
CN103983274A (zh) * | 2014-04-11 | 2014-08-13 | 湖北航天技术研究院总体设计所 | 一种适用于低精度无方位基准双轴转位设备的惯性测量单元标定方法 |
CN103994775A (zh) * | 2014-04-11 | 2014-08-20 | 湖北航天技术研究院总体设计所 | 一种适用于低精度有方位基准双轴转位设备的惯性测量单元标定方法 |
CN104121928A (zh) * | 2014-05-29 | 2014-10-29 | 湖北航天技术研究院总体设计所 | 一种适用于低精度有方位基准单轴转位设备的惯性测量单元标定方法 |
CN108458725A (zh) * | 2017-11-17 | 2018-08-28 | 北京计算机技术及应用研究所 | 捷联惯导***晃动基座上的***级标定方法 |
Non-Patent Citations (1)
Title |
---|
激光陀螺捷联惯导的动基座标定;杨玉孔等;《光学与光电技术》;20111031(第05期);第79-82页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111238532A (zh) | 2020-06-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108168574B (zh) | 一种基于速度观测的8位置捷联惯导***级标定方法 | |
CN108827299B (zh) | 一种基于改进四元数二阶互补滤波的飞行器姿态解算方法 | |
CN107655493B (zh) | 一种光纤陀螺sins六位置***级标定方法 | |
CN105180968B (zh) | 一种imu/磁强计安装失准角在线滤波标定方法 | |
CN113029199B (zh) | 一种激光陀螺惯导***的***级温度误差补偿方法 | |
CN110160554B (zh) | 一种基于寻优法的单轴旋转捷联惯导***标定方法 | |
CN106969783B (zh) | 一种基于光纤陀螺惯性导航的单轴旋转快速标定技术 | |
CN110398257A (zh) | Gps辅助的sins***快速动基座初始对准方法 | |
CN105806367B (zh) | 无陀螺惯性***误差标定方法 | |
CN103575299A (zh) | 利用外观测信息的双轴旋转惯导***对准及误差修正方法 | |
CN103245359A (zh) | 一种惯性导航***中惯性传感器固定误差实时标定方法 | |
CN108458725A (zh) | 捷联惯导***晃动基座上的***级标定方法 | |
CN103759729B (zh) | 采用捷联惯导的月球软着陆地面试验用初始姿态获取方法 | |
CN113503894B (zh) | 基于陀螺基准坐标系的惯导***误差标定方法 | |
CN111238532B (zh) | 一种适用于晃动基座环境的惯性测量单元标定方法 | |
CN106940193A (zh) | 一种基于Kalman滤波的船舶自适应摇摆标定方法 | |
CN111121773A (zh) | 一种mems惯性测量组合 | |
CN109708663B (zh) | 基于空天飞机sins辅助的星敏感器在线标定方法 | |
CN104121927A (zh) | 一种适用于低精度无方位基准单轴转位设备的惯性测量单元标定方法 | |
CN104121928A (zh) | 一种适用于低精度有方位基准单轴转位设备的惯性测量单元标定方法 | |
CN111912427B (zh) | 一种多普勒雷达辅助捷联惯导运动基座对准方法及*** | |
CN114877915A (zh) | 一种激光陀螺惯性测量组件g敏感性误差标定装置及方法 | |
CN109029499B (zh) | 一种基于重力视运动模型的加速度计零偏迭代寻优估计方法 | |
CN111044082A (zh) | 一种基于星敏感器辅助的陀螺误差参数在轨快速标定方法 | |
CN103983274A (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 |