CN115265590B - 一种双轴旋转惯导动态误差抑制方法 - Google Patents

一种双轴旋转惯导动态误差抑制方法 Download PDF

Info

Publication number
CN115265590B
CN115265590B CN202210835839.1A CN202210835839A CN115265590B CN 115265590 B CN115265590 B CN 115265590B CN 202210835839 A CN202210835839 A CN 202210835839A CN 115265590 B CN115265590 B CN 115265590B
Authority
CN
China
Prior art keywords
axis
accelerometer
error
inertial navigation
oriented
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
Application number
CN202210835839.1A
Other languages
English (en)
Other versions
CN115265590A (zh
Inventor
涂勇强
蔡庆中
杨功流
李晶
尹洪亮
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN202210835839.1A priority Critical patent/CN115265590B/zh
Publication of CN115265590A publication Critical patent/CN115265590A/zh
Application granted granted Critical
Publication of CN115265590B publication Critical patent/CN115265590B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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

本发明公开了一种双轴旋转惯导动态误差抑制方法，步骤为：S1、构建坐标系，确定陀螺和加速度计的安装误差矩阵；S2、确定双轴旋转惯导在导航坐标系下的误差方程；S3、构建三十六维卡尔曼滤波器；S4、双轴转位机构根据设计的转位路径进行转动，将该过程中的角速率数据和比力数据作为标定数据；S5、利用卡尔曼滤波算法和卡尔曼滤波器处理标定数据，以得到卡尔曼滤波器中的状态量估计值；S6、利用标定参数得到标定后输出的角速率和比力、以及经过加速度计内杆臂补偿的加速度，进而解算得到导航结果；该方法全面考虑零偏、标度因素、安装误差、加速度计二阶非线性系数和加速度计内杆臂导致的双轴旋转惯导动态误差，应用简便、精度高，实用性佳。

Description

一种双轴旋转惯导动态误差抑制方法
技术领域
本发明涉及双轴旋转惯导误差标定与抑制技术领域，特别涉及一种双轴旋转惯导动态误差抑制方法。
背景技术
双轴旋转惯导由惯性测量单元(Initial Measurement Unit，简称为IMU)和转位机构组成，工作原理为：将IMU安装在转位机构上，使IMU在导航过程中可以相对固定的坐标系轴旋转，由于IMU通过特定的旋转，将惯性器件的常值误差调制成为均值为零的周期变化量，可以在现有的惯性器件精度水平上，大大的提高双轴旋转惯导长航时的导航精度。因此，双轴旋转惯导被广泛地应用在需要高精度惯性导航***的应用场合中，例如：舰艇、战车、飞机等运载体中。
在双轴旋转惯导的实际使用中，合理的双轴转位方案虽然可以完全调制IMU中惯性器件的常值零偏，但是无法调制由运载体运动激励惯性器件标度因数误差、安装误差、加速度计二阶非线性系数和加速度计内杆臂引起的动态误差。因此，为提高双轴旋转惯导在动态环境下的精度，需要对双轴旋转惯导的惯性器件标度因数误差、安装误差、加速度计二阶非线性系数和加速度计内杆臂进行精确标定，并在标定后进行补偿以抑制双轴旋转惯导的动态误差。
为了达到这个目的，已公开发明专利CN103575296B提供了一种双轴旋转惯导***自标定方法，该方法利用遗传算法优化标定试验路径并利用最小二乘法处理标定数据标定出了双轴旋转惯导IMU中惯性器件的零偏、标度因素和安装误差；已公开发明专利CN104165638B提出了一种双轴旋转惯导***多位置自主标定方法，该方法利用卡尔曼滤波和最小二乘法标定出了陀螺和加速度计的零偏、标度因素和安装误差；已公开发明专利CN108981751A公开了一种双轴旋转惯导***的八位置在线标定方法，该方法标定出了IMU的零偏、标度因素和安装误差。上述现有方法能较好地标定出IMU的零偏、标度因素和安装误差；然而，上述方案均未考虑加速度计二阶非线性系数和加速度计内杆臂引起的双轴旋转惯导动态误差，导致了双轴旋转惯导依然存在动态环境下出现较大误差的问题。
因此，为了克服上述现有方法中存在着的未考虑加速度计二阶非线性系数和加速度计内杆臂引起的双轴旋转惯导动态误差导致双轴旋转惯导在动态环境下误差较大的问题，需要提供一种简便、高精度的，全面考虑IMU的零偏、标度因素、安装误差、加速度计二阶非线性系数和加速度计内杆臂的双轴旋转惯导动态误差抑制方法。
发明内容
本发明的目的是提供一种解决目前现有的双轴旋转惯导误差标定与抑制方法中存在着的未考虑加速度计二阶非线性系数和加速度计内杆臂引起的双轴旋转惯导动态误差导致双轴旋转惯导在动态环境下误差较大的问题的双轴旋转惯导动态误差抑制方法。
为此，本发明技术方案如下：
一种双轴旋转惯导动态误差抑制方法，步骤如下：
S1、构建陀螺敏感轴系、加速度计敏感轴系、IMU坐标系、地心惯性坐标系、地球坐标系和导航坐标系，并确定陀螺的安装误差矩阵和加速度计的安装误差矩阵；
S2、基于IMU的零偏、标度因素、安装误差、加速度计二阶非线性系数和加速度计内杆臂引发的误差，确定双轴旋转惯导在导航坐标系下的误差方程；
S3、基于双轴旋转惯导的误差方程，构建三十六维卡尔曼滤波器，包括三十六维卡尔曼滤波器的状态量、三十六维卡尔曼滤波器滤波器的状态方程、三十六维卡尔曼滤波器滤波器的观测方程和卡尔曼滤波估计结果的反馈补偿形式；其中，三十六维卡尔曼滤波器的状态量包含双轴旋转惯导的姿态误差角速度误差δV，位置误差δP，陀螺的标定参数误差Xg，加速度计的标定参数误差Xa，加速度计二阶非线性系数标定误差XKa2和内杆臂误差Xr；
S4、标定试验：开机预热完成后，双轴转位机构根据设计的转位路径进行转动，保存该过程中由双轴旋转惯导的陀螺输出的角速率数据和加速度计输出的比力数据，作为标定数据；其中，转位路径应设计应满足：双轴旋转惯导的惯性测量单元(IMU)的每个轴均经过单方向的两次180°旋转，以及IMU的每个轴均历经指天、指地两个位置；
S5、利用卡尔曼滤波算法和步骤S3构建的卡尔曼滤波器对步骤S4的标定数据进行处理，使卡尔曼滤波器中的状态量从初值开始随步骤S4的多次转位而得到不断迭代和反馈补偿，最终获得卡尔曼滤波器中的状态量估计值；其中，三十六维状态量的初值基于双轴旋转的初对准结果和仪器参数进行赋值；在最优的卡尔曼滤波器中的状态量估计值中，关于惯性器件误差的状态量分量为标定参数；
S6、将标定参数代入陀螺的标定模型和加速度计的标定模型中，得到标定后输出的角速率和比力，并基于标定后的角速率和比力得到经过加速度计内杆臂补偿的加速度，进而解算得到双轴旋转惯导经动态误差实时补偿后的导航结果。
进一步地，步骤S1的具体实施步骤为：
S101、构建陀螺敏感轴系，即g系，该坐标系以三个陀螺的测量中心点为坐标原点，Xg轴与X向陀螺的敏感轴方向一致，Yg轴与Y向陀螺的敏感轴方向一致，Zg轴与Z向陀螺的敏感轴方向一致；
S102、构建加速度计敏感轴系，即a系，该坐标系以g系的坐标原点为坐标原点，Xa轴与X向加速度计的敏感轴方向一致，Ya轴与Y向加速度计的敏感轴方向一致，Za轴与Z向加速度计的敏感轴方向一致；
S103、构建IMU坐标系，即m系，该坐标系以g系的坐标原点为坐标原点，Xm轴与陀螺敏感轴Xg轴一致，Ym轴在陀螺敏感轴Xg轴与Yg轴构成的平面中，且与Yg轴相差小角度γyz，Zm轴为陀螺敏感轴Zg轴绕Xm轴旋转小角度γzx，再绕Ym轴旋转小角度γzy；
S104、构建地心惯性坐标系，即i系，该坐标系以地球中心为坐标原点，Xi和Yi轴在地球赤道平面内，且Xi轴指向春分点，Zi轴为地球自转轴，并指向北极；
S105、构建地球坐标系，即e系，该坐标系以地球中心为坐标原点，Xe和Ye轴在地球赤道平面内，且Xe指向本初子午线，Ze轴为地球自转轴，并指向北极；
S106、构建导航坐标系，即n系，该坐标系为东北天地理坐标系，其坐标原点为IMU坐标系的坐标原点，Xn轴指向东，Yn轴指向北，Zn轴指向天；
S107、基于陀螺敏感轴系和IMU坐标系，陀螺的安装误差矩阵的表达式为：
式中，为陀螺的安装误差矩阵，γyz、γzy和γzx为陀螺敏感轴系相对于IMU坐标系的三个安装误差角，具体地，γyz表示Y向陀螺敏感轴Yg轴绕m系的Zm轴的偏转角度，γzy表示Z向陀螺敏感轴Zg轴绕m系的Ym轴的偏转角度，γzx表示Z向陀螺敏感轴Zg轴绕m系的Xm轴的偏转角度；
S108、加速度计的安装误差矩阵由六个加速度计安装误差角构成，其表达式为：
式中，为加速度计的安装误差矩阵，αxz、αxy、αyz、αyx、αzy和αzx为加速度计敏感轴系相对于IMU坐标系的六个安装误差角，具体地，αxz表示X向加速度计敏感轴Xa轴绕m系的Zm轴的偏转角度，αxy表示X向加速度计敏感轴Xa轴绕m系的Ym轴的偏转角度，αyz表示Y向加速度计敏感轴Ya轴绕m系的Zm轴的偏转角度，αyx表示Y向加速度计敏感轴Ya轴绕m系的Xm轴的偏转角度，αzy表示Z向加速度计敏感轴Za轴绕m系的Ym轴的偏转角度，αzx表示Z向加速度计敏感轴Za轴绕m系的Xm轴的偏转角度。
进一步地，步骤S2的具体实施步骤为：
双轴旋转惯导在导航坐标系(n系)下的误差方程，其表达式为：
式中，为/>相对于时间的导数，/>为双轴旋转惯导的姿态误差角，和/>分别为双轴旋转惯导的东向误差角、北向误差角和天向误差角；/>为n系相对于i系的转动角速率，其由地球自转与载体运动产生；/>为导航解算中/>的估计误差；/>为m系到n系的坐标转换矩阵；/>为陀螺的测量误差；/>为δVn相对于时间的导数，δVn为双轴旋转惯导的速度误差，δVn＝[δVE δVN δVU]T，δVE、δVN和δVU分别为东向速度误差、北向速度误差和天向速度误差；fn为n系下的比力；/>为地球自转角速率；为载体绕地球运动产生的角速率；；Vn为双轴旋转惯导的速度，Vn＝[VE VN VU]T，VE、VN和VU分别为双轴旋转惯导的东向速度、北向速度和天向速度；/>为/>的误差；/>为/>的误差；δfm为加速度计的测量误差；/>为加速度计内杆臂误差；δg为重力矢量误差；/>为δL相对于时间的导数，δL为纬度误差；/>为δλ相对于时间的导数，δλ为经度误差；/>为δh相对于时间的导数，δh为高度误差；L、λ和h分别为当地地理纬度、经度和高度；RM和RN分别为当地地球子午圈和卯酉圈半径；
其中，陀螺的测量误差的表达式为：/>式中，/>为陀螺在m系下的测量值；δKG为陀螺的标度因数和安装误差的误差矩阵，
和/>分别为X向陀螺标度因数/>Y向陀螺标度因数/>和Z向陀螺标度因数/>的误差；δγyz、δγzy和δγzx分别为陀螺安装误差角γyz、γzy和γzx的误差；εm为陀螺的零偏误差，其表达式为：εm＝[εx εy εz]T，式中，εx、εy和εz分别为X向陀螺、Y向陀螺和Z向陀螺的零偏误差；
加速度计的测量误差δfm的表达式为：式中，/>为加速度计在m系下的测量值；δKA为加速度计的标度因数和安装误差的误差矩阵，和/>分别为X向加速度计标度因数/>Y向加速度计标度因数/>和Z向加速度计标度因数/>的误差；δαxz、δαxy、δαyz、δαyx、δαzy和δαzx分别为加速度计安装误差角αxz、αxy、αyz、αyx、αzy和αzx的误差；δKA2为加速度计的二阶非线性系数的标定误差矩阵，/>δKax2、δKay2和δKaz2分别为X向加速度计的二阶非线性系数Kax2的标定误差、Y向加速度计的二阶非线性系数Kay2的标定误差和Z向加速度计的二阶非线性系数Kaz2的标定误差；/>为加速度计的零偏误差，和/>分别为X向加速度计、Y向加速度计和Z向加速度计的零偏误差；
加速度计内杆臂误差的表达式为：/>式中，/> 和/>分别为内杆臂误差/>在m系中Xm轴、Ym轴和Zm轴上的分量，其中，和/>分别为陀螺在m系下的测量值/>在Xm轴、Ym轴和Zm轴上的分量；δrx、δry和δrz分别为X向加速度计内杆臂rx、Y向加速度计内杆臂ry和Z向加速度计内杆臂rz的误差。
进一步地，步骤S3的具体实施步骤如下：
S301、构建三十六维卡尔曼滤波器的状态量X36：
式中，式中，为双轴旋转惯导的姿态误差角，/>δV为速度误差，δV＝[δVE δVN δVU]T；δP为位置误差，δP＝[δL δλ δh]T；Xg为陀螺的标定参数误差，Xa为加速度计的标定参数误差，XKa2为加速度计二阶非线性系数标定误差，XKa2＝[δKax2 δKay2 δKaz2]T；Xr为内杆臂误差，Xr＝[δrx δry δrz]T；
S302、构建三十六维卡尔曼滤波器滤波器的状态方程：
式中，在F33中，
F36中，
G36为测量噪声输入矩阵为，u为测量噪声，/>其中，ug为陀螺的测量噪声，ua为加速度计的测量噪声；
S303、构建三十六维卡尔曼滤波器滤波器的观测方程：Z＝H36X36+v，
式中，Z为观测量，以速度为观测量且在标定时双轴旋转惯导保持静止状态，则Z＝[VE VN VU]T；H36为观测矩阵，H36＝[03×3 I3 03×30]，I3为三行三列的单位矩阵，即03×30为三行三十列的零矩阵，03×3为三行三列的零矩阵；v是观测噪声；
S304、确定卡尔曼滤波估计结果的反馈补偿形式：
式中，为/>的反对称矩阵，/>为双轴旋转惯导的/>的解算结果；/>为双轴旋转惯导的Vn的解算结果；/>为双轴旋转惯导的L的解算结果；/>为双轴旋转惯导的λ的解算结果；/>为双轴旋转惯导的h的解算结果；，/>为卡尔曼滤波估计结果反馈补偿前的KG；/>为卡尔曼滤波估计结果反馈补偿前的εm；，/>为卡尔曼滤波估计结果反馈补偿前的KA；/>为卡尔曼滤波估计结果反馈补偿前的/>为卡尔曼滤波估计结果反馈补偿前的KA2；/>为卡尔曼滤波估计结果反馈补偿前的rx；/>为卡尔曼滤波估计结果反馈补偿前的ry；/>为卡尔曼滤波估计结果反馈补偿前的rz。
进一步地，在步骤S4中，双轴转位机构的转位路径具体设计为：
基于次序0对应初始位置：双轴旋转惯导以其Xm轴朝东向，Ym轴朝北向，Zm轴朝天向的姿态，静止10分钟；
次序1：以中框轴为旋转轴，顺时针旋转90°；该状态下，双轴旋转惯导的Xm轴朝东向，Ym轴朝天向，Zm轴朝南向；转位完成后，静止180s；
次序2：以中框轴为旋转轴，顺时针旋转180°；该状态下，双轴旋转惯导的Xm轴朝东向，Ym轴朝地向，Zm轴朝北向；转位完成后，静止180s；
次序3：以中框轴为旋转轴，顺时针旋转180°；该状态下，双轴旋转惯导的Xm轴朝东向，Ym轴朝天向，Zm轴朝南向；转位完成后，静止180s；
次序4：以内框轴为旋转轴，顺时针旋转90°；该状态下，双轴旋转惯导的Xm轴朝天向，Ym轴朝西向，Zm轴朝南向；转位完成后，静止180s；
次序5：以内框轴为旋转轴，顺时针旋转180°；该状态下，双轴旋转惯导的Xm轴朝地向，Ym轴朝东向，Zm轴朝南向；转位完成后，静止180s；
次序6：以内框轴为旋转轴，顺时针旋转180°；该状态下，双轴旋转惯导的Xm轴朝天向，Ym轴朝西向，Zm轴朝南向；转位完成后，静止180s；
次序7：以中框轴为旋转轴，顺时针旋转90°；该状态下，双轴旋转惯导的Xm轴朝南向，Ym轴朝西向，Zm轴朝地向；转位完成后，静止180s；
次序8：以中框轴为旋转轴，顺时针旋转180°；该状态下，双轴旋转惯导的Xm轴朝北向，Ym轴朝西向，Zm轴朝天向；转位完成后，静止180s；
次序9：以中框轴为旋转轴，顺时针旋转180°；该状态下，双轴旋转惯导的Xm轴朝南向，Ym轴朝西向，Zm轴朝地向；转位完成后，静止180s；
次序10：以中框轴为旋转轴，顺时针旋转90°；该状态下，双轴旋转惯导的Xm轴朝地向，Ym轴朝西向，Zm轴朝北向；转位完成后，静止180s；
次序11：以中框轴为旋转轴，顺时针旋转90°；该状态下，双轴旋转惯导的Xm轴朝北向，Ym轴朝西向，Zm轴朝天向；转位完成后，静止180s；
次序12：以中框轴为旋转轴，顺时针旋转90°；该状态下，双轴旋转惯导的Xm轴朝天向，Ym轴朝西向，Zm轴朝南向；转位完成后，静止180s；
次序13：以内框轴为旋转轴，顺时针旋转90°；该状态下，双轴旋转惯导的Xm轴朝西向，Ym轴朝地向，Zm轴朝南向；转位完成后，静止180s；
次序14：以内框轴为旋转轴，顺时针旋转90°；该状态下，双轴旋转惯导的Xm轴朝地向，Ym轴朝东向，Zm轴朝南向；转位完成后，静止180s；
次序15：以内框轴为旋转轴，顺时针旋转90°；该状态下，双轴旋转惯导的Xm轴朝东向，Ym轴朝天向，Zm轴朝南向；转位完成后，静止180s；
次序16：以中框轴为旋转轴，顺时针旋转90°；该状态下，双轴旋转惯导的Xm轴朝东向，Ym轴朝南向，Zm轴朝地向；转位完成后，静止180s；
次序17：以中框轴为旋转轴，顺时针旋转90°；该状态下，双轴旋转惯导的Xm轴朝东向，Ym轴朝地向，Zm轴朝北向；转位完成后，静止180s；
次序18：以中框轴为旋转轴，顺时针旋转90°；该状态下，双轴旋转惯导的Xm轴朝东向，Ym轴朝北向，Zm轴朝天向；转位完成后，静止180s。
进一步地，在步骤S4中，一次标定实验中按照设计的转位路径进行1～3次双轴转位机构的转位。
进一步地，步骤S5的具体实施步骤为：
S501、对三十六维状态量进行初值赋值：和δV的初值分别由步骤S4中双轴旋转的初对准结果中的初始姿态角和初始速度进行赋值；δP的初值由步骤S4中双轴旋转惯导完成初对准后输出的初始位置减去实际实验位置得到的差值进行赋值；Xg、Xa、XKa2和Xr的初值根据双轴旋转惯导中陀螺和加速度计的仪器实际情况进行赋值；
S502、基于步骤S501确定的三十六维状态量初值作为卡尔曼滤波的状态量的初值，将步骤S301构建的三十六维状态量分别带入步骤S302构建的三十六维卡尔曼滤波器状态方程和步骤S303构建的三十六维卡尔曼滤波器观测方程中，采用卡尔曼滤波算法和步骤S3构建的卡尔曼滤波器处理标定数据，通过对卡尔曼滤波器中的状态量从初值开始不断迭代和反馈补偿，获得卡尔曼滤波器中的状态量估计值；其中，状态量估计值中关于惯性器件误差的状态量分量即准确的标定参数，包括：X向陀螺标度因数Y向陀螺标度因数/>和Z向陀螺标度因数/>陀螺安装误差角γyz、γzy和γzx；X向陀螺零偏/>Y向陀螺零偏和Z向陀螺零偏/>X向加速度计标度因数/>Y向加速度计标度因数/>和Z向加速度计标度因数/>的误差；加速度计安装误差角αxz、αxy、αyz、αyx、αzy和αzx；X向加速度计的二阶非线性系数Kax2、Y向加速度计的二阶非线性系数Kay2和Z向加速度计的二阶非线性系数Kaz2；X向加速度计零偏/>Y向加速度计零偏/>和Z向加速度计零偏/>X向加速度计内杆臂rx、Y向加度计内杆臂ry和Z向加速度计内杆臂rz；
进一步地，步骤S6的具体实施步骤为：
S601、采用现有的六十四次序旋转调制方案作为双轴旋转惯导导航过程中的转位机构旋转方案，并在每转动至一个位置后停止10s时间，获得旋转调制后的惯性器件输出结果，包括陀螺输出的角速率和加速度计输出的比力/>
S602、将由步骤S5获得的标定参数分别代入至陀螺标定模型和加速度计标定模型中，以获得陀螺标定后的输出和加速度计标定后的输出；其中，
(1)陀螺的标定模型为：
式中，
(2)加速度计的标定模型为：
式中，
S603、利用经过步骤S602得到的陀螺标定后的输出和加速度计标定后的输出，利用X向加速度计内杆臂rx、Y向加速度计内杆臂ry和Z向加速度计内杆臂rz对加速度数据进行实时补偿，其具体公式为：
式中，axb、ayb和azb为加速度计内杆臂补偿后的加速度实时结果在m系中Xm轴、Ym轴和Zm轴的分量；axc、ayc和azc为加速度计内杆臂补偿前的加速度实时结果在m系中Xm轴、Ym轴和Zm轴的分量；ωxc、ωyc和ωzc为角速率在m系中Xm轴、Ym轴和Zm轴的分量；
S604、将步骤S603中的axb、ayb、azb和ωxc、ωyc、ωzc带入惯导解算方程，获得双轴旋转惯导的导航结果，即采用本申请实现的对双轴旋转惯导动态误差进行实时补偿后的导航结果。
与现有技术相比，该双轴旋转惯导动态误差抑制方法的有益效果在于：
(1)本申请的双轴旋转惯导动态误差抑制方法简便且精度高，全面考虑IMU的零偏、标度因素、安装误差、加速度计二阶非线性系数和加速度计内杆臂导致的双轴旋转惯导动态误差；
(2)采用本申请的方法进行标定补偿后，由动态引起的速度误差明显减少，其效果可达90％，且能提高动态环境下双轴旋转惯导的定位精度约8米，证明本申请方法的正确性和准确性，即能很好抑制双轴旋转惯导动态误差，实用性佳。
附图说明
图1为本发明的双轴旋转惯导动态误差抑制方法的流程示意图；
图2为本发明在步骤S1中构建的陀螺敏感轴系和IMU坐标系间坐标系转换的示意图；
图3(a)为本发明实施例中双轴旋转惯导在进行动态误差补偿前后的东向速度误差对比的示意图；
图3(b)为本发明实施例中双轴旋转惯导在进行动态误差补偿前后的北向速度误差对比的示意图；
图4(a)为本发明实施例中双轴旋转惯导在进行动态误差补偿前后的东向位置误差对比的示意图；
图4(b)为本发明实施例中双轴旋转惯导在进行动态误差补偿前后的北向位置误差对比的示意图。
具体实施方式
下面结合附图及具体实施例对本发明做进一步的说明，但下述实施例绝非对本发明有任何限制。
如图1所示，该双轴旋转惯导动态误差抑制方法的具体实施步骤如下：
S1、构建陀螺敏感轴系(g系)、加速度计敏感轴系(a系)、IMU坐标系(m系)、地心惯性坐标系(i系)、地球坐标系(e系)和导航坐标系(n系)，并确定陀螺的安装误差矩阵和加速度计的安装误差矩阵；
具体地，在该步骤S1中，坐标系构建如下：
(1)陀螺敏感轴系，又简称为g系，其表示为：o-XgYgZg；该坐标系以三个陀螺的测量中心点为坐标原点(o点)，其Xg轴与X向陀螺的敏感轴方向一致，其Yg轴与Y向陀螺的敏感轴方向一致，其Zg轴与Z向陀螺的敏感轴方向一致；由于双轴旋转惯导的IMU中，三个陀螺无法保证严格正交安装，因此，g系不是正交坐标系；
(2)加速度计敏感轴系，又简称为a系，其表示为：o-XaYaZa；该坐标系的坐标原点与g系的坐标原点重合，即以g系的坐标原点(o点)为坐标原点，Xa轴与X向加速度计的敏感轴方向一致，Ya轴与Y向加速度计的敏感轴方向一致，Za轴与Z向加速度计的敏感轴方向一致；同样地，由于双轴旋转惯导的IMU中，三个加速度计无法保证严格正交安装，因此，a系也不是正交坐标系；
(3)IMU坐标系，又简称为m系，其表示为：o-XmYmZm；由于惯性导航解算需要将各惯性器件的输出投影到同一个正交坐标系中，因此，该坐标系为用于将陀螺和加速度计的输出投影为正交向量的正交坐标系；如图2所示，该坐标系的坐标原点与g系的坐标原点重合，即以g系的坐标原点(o点)为坐标原点，Xm轴与陀螺敏感轴Xg轴一致，Ym轴在陀螺敏感轴Xg轴与Yg轴构成的平面中，且与Yg轴相差小角度γyz，Zm轴为陀螺敏感轴Zg轴绕Xm轴旋转小角度γzx，再绕Ym轴旋转小角度γzy；
(4)地心惯性坐标系，又简称为i系，其表示为oe-XiYiZi，该坐标系作为惯性器件的输出的参考基准，其坐标原点oe为地球中心，Xi和Yi轴在地球赤道平面内，且Xi轴指向春分点(即赤道面与黄道面的交线再与天球相交的交点之一)，Zi轴为地球自转轴，并指向北极；
(5)地球坐标系，又简称为e系，其表示为oe-XeYeZe；该坐标系的坐标原点oe为地球中心，Xe和Ye轴在地球赤道平面内，且Xe指向本初子午线，Ze轴为地球自转轴，并指向北极；该地球坐标系与地球固连，e系相对于i系的角速率为地球自转角速率；
(6)导航坐标系，又简称为n系，其表示为o-XnYnZn；该坐标系作为惯导在求解导航参数时采用的坐标系，具体为东北天地理坐标系，其坐标原点o为IMU坐标系的坐标原点，Xn轴指向东，Yn轴指向北，Zn轴指向天，n系相对于e系的方位关系就是载体的地理位置；
基于上述构建的陀螺敏感轴系(g系)和IMU坐标系(m系)，可以确定陀螺的安装误差矩阵的表达式为：
式中，为陀螺的安装误差矩阵，γyz、γzy和γzx为陀螺敏感轴系相对于IMU坐标系的三个安装误差角，具体地，γyz表示Y向陀螺敏感轴Yg轴绕m系的Zm轴的偏转角度，γzy表示Z向陀螺敏感轴Zg轴绕m系的Ym轴的偏转角度，γzx表示Z向陀螺敏感轴Zg轴绕m系的Xm轴的偏转角度；
另外，由于陀螺组件与加速度计组件之间存在小角度安装误差，因此，为了将加速度计的测量信息转换到IMU坐标系下，加速度计的安装误差矩阵由六个加速度计安装误差角构成，其表达式为：
式中，为加速度计的安装误差矩阵，αxz、αxy、αyz、αyx、αzy和αzx为加速度计敏感轴系相对于IMU坐标系的六个安装误差角，具体地，αxz表示X向加速度计敏感轴Xa轴绕m系的Zm轴的偏转角度，αxy表示X向加速度计敏感轴Xa轴绕m系的Ym轴的偏转角度，αyz表示Y向加速度计敏感轴Ya轴绕m系的Zm轴的偏转角度，αyx表示Y向加速度计敏感轴Ya轴绕m系的Xm轴的偏转角度，αzy表示Z向加速度计敏感轴Za轴绕m系的Ym轴的偏转角度，αzx表示Z向加速度计敏感轴Za轴绕m系的Xm轴的偏转角度；
S2、基于IMU的零偏、标度因素、安装误差、加速度计二阶非线性系数和加速度计内杆臂引发的误差，确定双轴旋转惯导在导航坐标系(n系)下的误差方程，其表达式为：
/>
式中，为/>相对于时间的导数，/>为双轴旋转惯导的姿态误差角，和/>分别为双轴旋转惯导的东向误差角、北向误差角和天向误差角；/>为n系相对于i系的转动角速率，其由地球自转与载体运动产生；/>为导航解算中/>的估计误差；/>为m系到n系的坐标转换矩阵；/>为陀螺的测量误差；/>为δVn相对于时间的导数，δVn为双轴旋转惯导的速度误差，δVn＝[δVE δVN δVU]T，δVE、δVN和δVU分别为东向速度误差、北向速度误差和天向速度误差；fn为n系下的比力；/>为地球自转角速率；为载体绕地球运动产生的角速率；；Vn为双轴旋转惯导的速度，Vn＝[VE VN VU]T，VE、VN和VU分别为双轴旋转惯导的东向速度、北向速度和天向速度；/>为/>的误差；/>为/>的误差；δfm为加速度计的测量误差；/>为加速度计内杆臂误差；δg为重力矢量误差；/>为δL相对于时间的导数，δL为纬度误差；/>为δλ相对于时间的导数，δλ为经度误差；/>为δh相对于时间的导数，δh为高度误差；L、λ和h分别为当地地理纬度、经度和高度；RM和RN分别为当地地球子午圈和卯酉圈半径；
其中，陀螺的测量误差的表达式为：/>
式中，为陀螺在m系下的测量值；δKG为陀螺的标度因数和安装误差的误差矩阵，/>和/>分别为X向陀螺标度因数Y向陀螺标度因数/>和Z向陀螺标度因数/>的误差；δγyz、δγzy和δγzx分别为陀螺安装误差角γyz、γzy和γzx的误差；
εm为陀螺的零偏误差，其表达式为：εm＝[εx εy εz]T，
式中，εx、εy和εz分别为X向陀螺、Y向陀螺和Z向陀螺的零偏误差；
加速度计的测量误差δfm的表达式为：
式中，为加速度计在m系下的测量值；δKA为加速度计的标度因数和安装误差的误差矩阵，/>和/>分别为X向加速度计标度因数/>Y向加速度计标度因数/>和Z向加速度计标度因数/>的误差；δαxz、δαxy、δαyz、δαyx、δαzy和δαzx分别为加速度计安装误差角αxz、αxy、αyz、αyx、αzy和αzx的误差；δKA2为加速度计的二阶非线性系数的标定误差矩阵，/>δKax2、δKay2和δKaz2分别为X向加速度计的二阶非线性系数Kax2的标定误差、Y向加速度计的二阶非线性系数Kay2的标定误差和Z向加速度计的二阶非线性系数Kaz2的标定误差；/>为加速度计的零偏误差，/>和/>分别为X向加速度计、Y向加速度计和Z向加速度计的零偏误差；
加速度计内杆臂误差的表达式为：/>
式中，和/>分别为内杆臂误差/>在m系中Xm轴、Ym轴和Zm轴上的分量，/> 其中，/>和/>分别为陀螺在m系下的测量值/>在Xm轴、Ym轴和Zm轴上的分量；δrx、δry和δrz分别为X向加速度计内杆臂rx、Y向加速度计内杆臂ry和Z向加速度计内杆臂rz的误差；
S3、基于步骤S2确定的双轴旋转惯导的误差方程，构建三十六维卡尔曼滤波器，其具体实施步骤如下：
S301、构建三十六维卡尔曼滤波器的状态量X36为：
式中，为双轴旋转惯导的姿态误差角，/>δV为速度误差，δV＝[δVE δVN δVU]T；δP为位置误差，δP＝[δL δλ δh]T；Xg为陀螺的标定参数误差，Xa为加速度计的标定参数误差，XKa2为加速度计二阶非线性系数标定误差，XKa2＝[δKax2 δKay2 δKaz2]T；Xr为内杆臂误差，Xr＝[δrx δry δrz]T；
S302、构建三十六维卡尔曼滤波器滤波器的状态方程：
式中，其中，在F33中，每个非零矩阵元素为：/>
/>
F36中，
G36为测量噪声输入矩阵为，u为测量噪声，/>其中，ug为陀螺的测量噪声，ua为加速度计的测量噪声；
S303、构建三十六维卡尔曼滤波器滤波器的观测方程为：
Z＝H36X36+v，
式中，Z为观测量，以速度为观测量且在标定时双轴旋转惯导保持静止状态，则Z＝[VE VN VU]T；H36为观测矩阵，H36＝[03×3 I3 03×30]，I3为三行三列的单位矩阵，即03×30为三行三十列的零矩阵，03×3为三行三列的零矩阵；v是观测噪声；
S304、确定卡尔曼滤波估计结果的反馈补偿形式：
式中，为/>的反对称矩阵，即/>为双轴旋转惯导的/>的解算结果；/>为双轴旋转惯导的Vn的解算结果；/>为双轴旋转惯导的L的解算结果；/>为双轴旋转惯导的λ的解算结果；/>为双轴旋转惯导的h的解算结果；，/>为卡尔曼滤波估计结果反馈补偿前的KG；/>为卡尔曼滤波估计结果反馈补偿前的εm；，/>为卡尔曼滤波估计结果反馈补偿前的KA；/>为卡尔曼滤波估计结果反馈补偿前的/>为卡尔曼滤波估计结果反馈补偿前的KA2；/>为卡尔曼滤波估计结果反馈补偿前的rx；/>为卡尔曼滤波估计结果反馈补偿前的ry；/>为卡尔曼滤波估计结果反馈补偿前的rz；
S4、设计标定试验中双轴转位机构的转位路径并进行标定试验；
该步骤S4的具体步骤为：
S401、设计标定试验中双轴转位机构的转位路径；
双轴旋转惯导***的双轴转位机构包括有中框轴与内框轴，二者相互垂直并在同一平面内；通过控制中框轴和内框轴的旋转来控制双轴转位机构的转位路径，使得IMU坐标系不断变换姿态；设定双轴转位机构的转速为5°/s；转位路径如下表1所示，具体为：
次序0：双轴转位机构在初始位置上，此时，其中框轴朝东向，内框轴朝天向，对应地，双轴旋转惯导的Xm轴朝东向，Ym轴朝北向，Zm轴朝天向；双轴转位机构初始静止10分钟，使双轴旋转惯导在静止过程中进行初对准，并获得初始的包括姿态角和速度在内的导航解算结果，之后即进入导航解算状态；
次序1：以中框轴为旋转轴，顺时针旋转90°，旋转后的转位机构姿态记为次序1，对应双轴旋转惯导的具体姿态为：Xm轴朝东向，Ym轴朝天向，Zm轴朝南向；并在旋转后的次序1的位置静止180s；
次序2：以中框轴为旋转轴，顺时针旋转180°，旋转后的转位机构姿态记为次序2，对应双轴旋转惯导的具体姿态为：Xm轴朝东向，Ym轴朝地向，Zm轴朝北向；并在旋转后的次序2的位置静止180s；
次序3：以中框轴为旋转轴，顺时针旋转180°，旋转后的转位机构姿态记为次序3，对应双轴旋转惯导的具体姿态为：Xm轴朝东向，Ym轴朝天向，Zm轴朝南向；并在旋转后的次序3的位置静止180s；
次序4：以内框轴为旋转轴，顺时针旋转90°，旋转后的转位机构姿态记为次序4，对应双轴旋转惯导的具体姿态为：Xm轴朝天向，Ym轴朝西向，Zm轴朝南向；并在旋转后的次序4的位置静止180s；
次序5：以内框轴为旋转轴，顺时针旋转180°，旋转后的转位机构姿态记为次序5，对应双轴旋转惯导的具体姿态为：Xm轴朝地向，Ym轴朝东向，Zm轴朝南向；并在旋转后的次序5的位置静止180s；
次序6：以内框轴为旋转轴，顺时针旋转180°，旋转后的转位机构姿态记为次序6，对应双轴旋转惯导的具体姿态为：Xm轴朝天向，Ym轴朝西向，Zm轴朝南向；并在旋转后的次序6的位置静止180s；
次序7：以中框轴为旋转轴，顺时针旋转90°，旋转后的转位机构姿态记为次序7，对应双轴旋转惯导的具体姿态为：Xm轴朝南向，Ym轴朝西向，Zm轴朝地向；并在旋转后的次序7的位置静止180s；
次序8：以中框轴为旋转轴，顺时针旋转180°，旋转后的转位机构姿态记为次序8，对应双轴旋转惯导的具体姿态为：Xm轴朝北向，Ym轴朝西向，Zm轴朝天向；并在旋转后的次序8的位置静止180s；
次序9：以中框轴为旋转轴，顺时针旋转180°，旋转后的转位机构姿态记为次序9，对应双轴旋转惯导的具体姿态为：Xm轴朝南向，Ym轴朝西向，Zm轴朝地向；并在旋转后的次序9的位置静止180s；
次序10：以中框轴为旋转轴，顺时针旋转90°，旋转后的转位机构姿态记为次序10，对应双轴旋转惯导的具体姿态为：Xm轴朝地向，Ym轴朝西向，Zm轴朝北向；并在旋转后的次序10的位置静止180s；
次序11：以中框轴为旋转轴，顺时针旋转90°，旋转后的转位机构姿态记为次序11，对应双轴旋转惯导的具体姿态为：Xm轴朝北向，Ym轴朝西向，Zm轴朝天向；并在旋转后的次序11的位置静止180s；
次序12：以中框轴为旋转轴，顺时针旋转90°，旋转后的转位机构姿态记为次序12，对应双轴旋转惯导的具体姿态为：Xm轴朝天向，Ym轴朝西向，Zm轴朝南向；并在旋转后的次序12的位置静止180s；
次序13：以内框轴为旋转轴，顺时针旋转90°，旋转后的转位机构姿态记为次序13，对应双轴旋转惯导的具体姿态为：Xm轴朝西向，Ym轴朝地向，Zm轴朝南向；并在旋转后的次序13的位置静止180s；
次序14：以内框轴为旋转轴，顺时针旋转90°，旋转后的转位机构姿态记为次序14，对应双轴旋转惯导的具体姿态为：Xm轴朝地向，Ym轴朝东向，Zm轴朝南向；并在旋转后的次序14的位置静止180s；
次序15：以内框轴为旋转轴，顺时针旋转90°，旋转后的转位机构姿态记为次序15，对应双轴旋转惯导的具体姿态为：Xm轴朝东向，Ym轴朝天向，Zm轴朝南向；并在旋转后的次序15的位置静止180s；
次序16：以中框轴为旋转轴，顺时针旋转90°，旋转后的转位机构姿态记为次序16，对应双轴旋转惯导的具体姿态为：Xm轴朝东向，Ym轴朝南向，Zm轴朝地向；并在旋转后的次序16的位置静止180s；
次序17：以中框轴为旋转轴，顺时针旋转90°，旋转后的转位机构姿态记为次序17，对应双轴旋转惯导的具体姿态为：Xm轴朝东向，Ym轴朝地向，Zm轴朝北向；并在旋转后的次序17的位置静止180s；
次序18：以中框轴为旋转轴，顺时针旋转90°，旋转后的转位机构姿态记为次序18，对应双轴旋转惯导的具体姿态为：Xm轴朝东向，Ym轴朝北向，Zm轴朝天向；并在旋转后的次序18的位置静止180s；
表1：
整个转位过程中共涉及18次转序，且每次转动至指定位置后停止180s，即整个转位路径可以在1小时内完成；
在上述转位路径中，前九次转序(次序1～次序9)的目的在于激励陀螺标度因数误差和安装误差，因此共包含了IMU的每个轴单方向的两次180°旋转；后九次转序(次序10～次序18)的目的在于激励加速度计标度因数和安装误差，因此共包含了IMU的每个轴的指天、指地两个位置；与此同时，上述标定路径还可以满足加速度计二阶非线性系数的激励和解耦；由于实际标定中，陀螺的随机噪声比较大，卡尔曼滤波器中陀螺零偏误差的估计通常需要较长时间，因此在一次标定中进行两次的标定路径转位，既可以保证各标定参数误差的估计曲线完全收敛，最终获得卡尔曼滤波器中的状态量估计值；；
S402、开机预热至少四个小时，以减小温度变化对标定的影响；
S403、完成标定实验，即按照步骤S401设计的标定转位路径控制双轴转位机构转动，并保存由陀螺输出的角速率数据和加速度计输出的比力数据，作为标定数据；
S5、利用步骤S3构建的卡尔曼滤波器处理由步骤S4得到标定数据，以获得标定参数，其具体实施过程如下：
S501、给定三十六维卡尔曼滤波器中三十六维状态量的初值，即给定的初值；其中，
和δV＝[δVE δVN δVU]T的初值由双轴旋转的初对准结果中的初始姿态角和初始速度直接赋值，具体为：将通过步骤S4在初始状态下转位机构静止10分钟进行双轴旋转惯导初对准后所得到姿态的初值：/>和速度的初值：V0＝[VE0 VN0 VU0]T作为/>和δV的初值，即初始状态下，/>δV＝V0；
δP的初值由步骤S4的双轴旋转惯导完成初对准后输出的初始位置减去实际实验室位置的差值进行赋值，具体为：1)利用GPS获得实验室位置为Pt＝[Lt λt ht]T，其中，Lt、λt和ht分别为实验室的纬度、经度和高度；2)步骤S4的双轴旋转惯导完成初对准后输出的初始位置为P0＝[L0 λ0 h0]T；3)δP的初值设置为：[L0-Lt λ0-λt h0-ht]T；
Xg、Xa、XKa2和Xr的初值根据具体双轴旋转惯导所用的陀螺和加速度计的仪器实际情况后进行赋值(即通过仪器类型、精度水平以及机械结构设计特点进行人为考量后确定初值)；
S502、基于步骤S501确定的三十六维状态量初值作为卡尔曼滤波的状态量的初值，将步骤S301构建的三十六维状态量带入步骤S302构建的三十六维卡尔曼滤波器状态方程和步骤S303构建的三十六维卡尔曼滤波器观测方程，利用卡尔曼滤波算法和步骤S3构建的卡尔曼滤波器，处理标定数据，使卡尔曼滤波器中的状态量从初值开始随步骤S4的多次转位而得到不断迭代和反馈补偿，反馈补偿采用S304确定的估计结果的反馈补偿补偿方式，最终获得卡尔曼滤波器中的状态量估计值；
其中，在状态量估计值中，关于惯性器件误差的状态量分量即为经过步骤S5获得准确的标定参数，包括：X向陀螺标度因数Y向陀螺标度因数/>和Z向陀螺标度因数陀螺安装误差角γyz、γzy和γzx；X向陀螺零偏/>Y向陀螺零偏/>和Z向陀螺零偏X向加速度计标度因数/>Y向加速度计标度因数/>和Z向加速度计标度因数/>的误差；加速度计安装误差角αxz、αxy、αyz、αyx、αzy和αzx；X向加速度计的二阶非线性系数Kax2、Y向加速度计的二阶非线性系数Kay2和Z向加速度计的二阶非线性系数Kaz2；X向加速度计零偏Y向加速度计零偏/>和Z向加速度计零偏/>X向加速度计内杆臂rx、Y向加速度计内杆臂ry和Z向加速度计内杆臂rz；
S6、将标定参数代入标定模型中，以补偿双轴旋转惯导的输出；其具体实施步骤为：
S601、采用常用的转位机构旋转方案(例如：六十四次序旋转调制方案)作为双轴旋转惯导导航过程中的转位机构旋转方案，并在每转动至一个位置后停止10s时间，获得旋转调制后的惯性器件输出结果，包括陀螺输出的角速率和加速度计输出的比力/>
其中，在该步骤S603中，六十四次序旋转调制方案是一种针对双轴旋转惯导使用过程中的公知的转位方案，即安装在运载体上进行导航时候的转位方案，因此在本实施例中不作详细说明；
S602、将由步骤S5获得的标定参数分别代入至陀螺标定模型和加速度计标定模型中，以获得陀螺标定后的输出和加速度计标定后的输出；其中，
(1)陀螺的标定模型为：/>
式中，
(2)加速度计的标定模型为：
式中，
S603、利用经过步骤S602得到的陀螺标定后的输出和加速度计标定后的输出，利用X向加速度计内杆臂rx、Y向加速度计内杆臂ry和Z向加速度计内杆臂rz对加速度数据进行实时补偿，其具体公式为：
式中，axb、ayb和azb为加速度计内杆臂补偿后的加速度实时结果在m系中Xm轴、Ym轴和Zm轴的分量；axc、ayc和azc为加速度计内杆臂补偿前的加速度实时结果fm在m系中Xm轴、Ym轴和Zm轴的分量；ωxc、ωyc和ωzc为角速率在m系中Xm轴、Ym轴和Zm轴的分量；
S604、将步骤S603中的axb、ayb、azb和ωxc、ωyc和ωzc带入惯导解算方程，获得双轴旋转惯导的导航结果；
经过步骤S604得到的导航结果，即为采用本申请实现的对双轴旋转惯导动态误差进行实时补偿后的导航结果。
为验证本发明提出的双轴旋转惯导动态误差抑制方法的正确性和准确性，选用一套双轴旋转惯导进行了标定实验，选用的双轴旋转惯导内的IMU由三个零偏稳定性为0.01°/h的激光陀螺仪和三个零偏稳定性为10μg的加速度计组成；选用的双轴旋转惯导的转位机构的姿态控制精度为5″(1σ)。
首先，采用本发明提供的方法进行误差参数的标定，X36中Xg、Xa、XKa2和Xr的初值根据所选用的双轴旋转惯导所使用的陀螺、加速度计类型和精度水平以及机械结构设计确定为：X向陀螺标度因数Y向陀螺标度因数/>和Z向陀螺标度因数/>的初值均为-0.84″/pulse；陀螺安装误差角γyz、γzy和γzx的初值均为10′；X向陀螺零偏/>Y向陀螺零偏/>和Z向陀螺零偏/>的初值均为0.01°/h；X向加速度计标度因数/>Y向加速度计标度因数/>和Z向加速度计标度因数/>的初值均为450μm/s/pulse；加速度计安装误差角αxz、αxy、αyz、αyx、αzy和αzx的初值均为10′；X向加速度计的二阶非线性系数Kax2、Y向加速度计的二阶非线性系数Kay2和Z向加速度计的二阶非线性系数Kaz2的初值均为20μg/g2；X向加速度计零偏/>Y向加速度计零偏/>和Z向加速度计零偏/>的初值均为1000μg；X向加速度计内杆臂rx、Y向加速度计内杆臂ry和Z向加速度计内杆臂rz的初值均为3cm；
根据设置的初值和本发明提供的方法获得的标定参数如下：
表2：
在获得标定参数后，将选用的双轴惯导安装在试验车上进行车载动态实验，实验包括10分钟静态初始对准和10分钟导航解算，首先在试验车进行情况下进行10分钟静态初始对准，在对准完成后进入导航状态，试验车开始以15m/s的速度直线开动模拟动态环境，并在双轴旋转惯导进入导航状态约4分钟后，试验车180°转弯，在转弯后继续以15m/s的速度直线开动；在导航状态约8分钟后，试验车再次180°转弯，两次旋转方向相同。
分别采用未考虑加速度计二次非线性参数、加速度计内杆臂误差的标定方法和本发明提出的方法的两组参数进行对同一组实验数据进行导航解算，将导航解算结果与安装于试验车上的GPS输出的速度和位置进行对比。
如图3(a)所示为双轴旋转惯导在进行动态误差补偿前、后的东向速度误差对比的示意图；如图3(b)所示为双轴旋转惯导在进行动态误差补偿前、后的北向速度误差对比的示意图。从上述两幅图中对比可以看出，在未补偿加速度计二次非线性参数、加速度计内杆臂误差的情况下，双轴旋转惯导180°转向将导致在东向和北向产生约0.025m/s的速度误差，而采用本发明提出的方法进行标定补偿后，由动态引起的速度误差明显减少，效果可达90％；
如图4(a)所示为双轴旋转惯导在进行动态误差补偿前、后的东向位置误差对比的示意图；如图4(b)所示为双轴旋转惯导在进行动态误差补偿前、后的北向位置误差对比的示意图。从上述两幅图中对比可以看出，对比未补偿加速度计二次非线性参数、加速度计内杆臂误差和本发明得到的导航结果的位置误差可见，在短航时的导航环境中，通过本发明提供的方法可以提高动态环境下双轴旋转惯导的定位精度约8米。
综上所述，上述实验证明了本发明提供的双轴旋转惯导动态误差抑制方法的正确性和准确性，具有能够能很好地抑制双轴旋转惯导动态误差的效果，有很好的实用性。
本发明未详细公开的部分属于本领域的公知技术。尽管上面对本发明说明性的具体实施方式进行了描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化时显而易见的，一切利用本发明构思的发明创造均为保护之列。

Claims (8)

1.一种双轴旋转惯导动态误差抑制方法，其特征在于，步骤如下：
S1、构建陀螺敏感轴系、加速度计敏感轴系、IMU坐标系、地心惯性坐标系、地球坐标系和导航坐标系，并确定陀螺的安装误差矩阵和加速度计的安装误差矩阵；
S2、基于IMU的零偏、标度因素、安装误差、加速度计二阶非线性系数和加速度计内杆臂引发的误差，确定双轴旋转惯导在导航坐标系下的误差方程；
S3、基于双轴旋转惯导的误差方程，构建三十六维卡尔曼滤波器，包括三十六维卡尔曼滤波器的状态量、三十六维卡尔曼滤波器滤波器的状态方程、三十六维卡尔曼滤波器滤波器的观测方程和卡尔曼滤波估计结果的反馈补偿形式；其中，三十六维卡尔曼滤波器的状态量包含双轴旋转惯导的姿态误差角速度误差δV，位置误差δP，陀螺的标定参数误差Xg，加速度计的标定参数误差Xa，加速度计二阶非线性系数标定误差XKa2和内杆臂误差Xr；
S4、标定试验：开机预热完成后，双轴转位机构根据设计的转位路径进行转动，保存该过程中由双轴旋转惯导的陀螺输出的角速率数据和加速度计输出的比力数据，作为标定数据；其中，转位路径应设计应满足：双轴旋转惯导惯性测量单元的每个轴均经过单方向的两次180°旋转，以及惯性测量单元的每个轴均历经指天、指地两个位置；
S5、利用卡尔曼滤波算法和步骤S3构建的卡尔曼滤波器对步骤S4的标定数据进行处理，使卡尔曼滤波器中的状态量从初值开始随步骤S4的多次转位而得到不断迭代和反馈补偿，最终获得卡尔曼滤波器中的状态量估计值；其中，三十六维状态量的初值基于双轴旋转的初对准结果和仪器参数进行赋值；在最优的卡尔曼滤波器中的状态量估计值中，关于惯性器件误差的状态量分量为标定参数；
S6、将标定参数代入陀螺的标定模型和加速度计的标定模型中，得到标定后输出的角速率和比力，并基于标定后的角速率和比力得到经过加速度计内杆臂补偿的加速度，进而解算得到双轴旋转惯导经动态误差实时补偿后的导航结果。
2.根据权利要求1所述的双轴旋转惯导动态误差抑制方法，其特征在于，步骤S1的具体实施步骤为：
S101、构建陀螺敏感轴系，即g系，该坐标系以三个陀螺的测量中心点为坐标原点，Xg轴与X向陀螺的敏感轴方向一致，Yg轴与Y向陀螺的敏感轴方向一致，Zg轴与Z向陀螺的敏感轴方向一致；
S102、构建加速度计敏感轴系，即a系，该坐标系以g系的坐标原点为坐标原点，Xa轴与X向加速度计的敏感轴方向一致，Ya轴与Y向加速度计的敏感轴方向一致，Za轴与Z向加速度计的敏感轴方向一致；
S103、构建IMU坐标系，即m系，该坐标系以g系的坐标原点为坐标原点，Xm轴与陀螺敏感轴Xg轴一致，Ym轴在陀螺敏感轴Xg轴与Yg轴构成的平面中，且与Yg轴相差小角度γyz，Zm轴为陀螺敏感轴Zg轴绕Xm轴旋转小角度γzx，再绕Ym轴旋转小角度γzy；
S104、构建地心惯性坐标系，即i系，该坐标系以地球中心为坐标原点，Xi和Yi轴在地球赤道平面内，且Xi轴指向春分点，Zi轴为地球自转轴，并指向北极；
S105、构建地球坐标系，即e系，该坐标系以地球中心为坐标原点，Xe和Ye轴在地球赤道平面内，且Xe指向本初子午线，Ze轴为地球自转轴，并指向北极；
S106、构建导航坐标系，即n系，该坐标系为东北天地理坐标系，其坐标原点为IMU坐标系的坐标原点，Xn轴指向东，Yn轴指向北，Zn轴指向天；
S107、基于陀螺敏感轴系和IMU坐标系，陀螺的安装误差矩阵的表达式为：
式中，为陀螺的安装误差矩阵，γyz、γzy和γzx为陀螺敏感轴系相对于IMU坐标系的三个安装误差角，具体地，γyz表示Y向陀螺敏感轴Yg轴绕m系的Zm轴的偏转角度，γzy表示Z向陀螺敏感轴Zg轴绕m系的Ym轴的偏转角度，γzx表示Z向陀螺敏感轴Zg轴绕m系的Xm轴的偏转角度；
S108、加速度计的安装误差矩阵由六个加速度计安装误差角构成，其表达式为：
式中，为加速度计的安装误差矩阵，αxz、αxy、αyz、αyx、αzy和αzx为加速度计敏感轴系相对于IMU坐标系的六个安装误差角，具体地，αxz表示X向加速度计敏感轴Xa轴绕m系的Zm轴的偏转角度，αxy表示X向加速度计敏感轴Xa轴绕m系的Ym轴的偏转角度，αyz表示Y向加速度计敏感轴Ya轴绕m系的Zm轴的偏转角度，αyx表示Y向加速度计敏感轴Ya轴绕m系的Xm轴的偏转角度，αzy表示Z向加速度计敏感轴Za轴绕m系的Ym轴的偏转角度，αzx表示Z向加速度计敏感轴Za轴绕m系的Xm轴的偏转角度。
3.根据权利要求2所述的双轴旋转惯导动态误差抑制方法，其特征在于，步骤S2的具体实施步骤为：
双轴旋转惯导在导航坐标系下的误差方程，其表达式为：
式中，为/>相对于时间的导数，/>为双轴旋转惯导的姿态误差角，/> 和/>分别为双轴旋转惯导的东向误差角、北向误差角和天向误差角；/>为n系相对于i系的转动角速率，其由地球自转与载体运动产生；/>为导航解算中/>的估计误差；为m系到n系的坐标转换矩阵；/>为陀螺的测量误差；/>为δVn相对于时间的导数，δVn为双轴旋转惯导的速度误差，δVn＝[δVE δVN δVU]T，δVE、δVN和δVU分别为东向速度误差、北向速度误差和天向速度误差；fn为n系下的比力；/>为地球自转角速率；/>为载体绕地球运动产生的角速率；；Vn为双轴旋转惯导的速度，Vn＝[VE VN VU]T，VE、VN和VU分别为双轴旋转惯导的东向速度、北向速度和天向速度；/>为/>的误差；/>为/>的误差；δfm为加速度计的测量误差；/>为加速度计内杆臂误差；δg为重力矢量误差；/>为δL相对于时间的导数，δL为纬度误差；/>为δλ相对于时间的导数，δλ为经度误差；/>为δh相对于时间的导数，δh为高度误差；L、λ和h分别为当地地理纬度、经度和高度；RM和RN分别为当地地球子午圈和卯酉圈半径；
其中，陀螺的测量误差的表达式为：
式中，为陀螺在m系下的测量值；δKG为陀螺的标度因数和安装误差的误差矩阵， 和/>分别为X向陀螺标度因数/>Y向陀螺标度因数/>和Z向陀螺标度因数/>的误差；δγyz、δγzy和δγzx分别为陀螺安装误差角γyz、γzy和γzx的误差；
εm为陀螺的零偏误差，其表达式为：
εm＝[εx εy εz]T，
式中，εx、εy和εz分别为X向陀螺、Y向陀螺和Z向陀螺的零偏误差；
加速度计的测量误差δfm的表达式为：
式中，为加速度计在m系下的测量值；δKA为加速度计的标度因数和安装误差的误差矩阵，/> 和/>分别为X向加速度计标度因数/>Y向加速度计标度因数/>和Z向加速度计标度因数/>的误差；δαxz、δαxy、δαyz、δαyx、δαzy和δαzx分别为加速度计安装误差角αxz、αxy、αyz、αyx、αzy和αzx的误差；δKA2为加速度计的二阶非线性系数的标定误差矩阵，/>δKax2、δKay2和δKaz2分别为X向加速度计的二阶非线性系数Kax2的标定误差、Y向加速度计的二阶非线性系数Kay2的标定误差和Z向加速度计的二阶非线性系数Kaz2的标定误差；/>为加速度计的零偏误差， 和/>分别为X向加速度计、Y向加速度计和Z向加速度计的零偏误差；
加速度计内杆臂误差的表达式为：
式中，和/>分别为内杆臂误差/>在m系中Xm轴、Ym轴和Zm轴上的分量， 其中，/>和/>分别为陀螺在m系下的测量值/>在Xm轴、Ym轴和Zm轴上的分量；δrx、δry和δrz分别为X向加速度计内杆臂rx、Y向加速度计内杆臂ry和Z向加速度计内杆臂rz的误差。
4.根据权利要求3所述的双轴旋转惯导动态误差抑制方法，其特征在于，步骤S3的具体实施步骤如下：
S301、构建三十六维卡尔曼滤波器的状态量X36：
式中，式中，为双轴旋转惯导的姿态误差角，/>δV为速度误差，δV＝[δVE δVN δVU]T；δP为位置误差，δP＝[δL δλ δh]T；Xg为陀螺的标定参数误差，Xa为加速度计的标定参数误差，XKa2为加速度计二阶非线性系数标定误差，XKa2＝[δKax2 δKay2 δKaz2]T；Xr为内杆臂误差，Xr＝[δrx δry δrz]T；
S302、构建三十六维卡尔曼滤波器滤波器的状态方程：
式中，在F33中，每个非零矩阵元素为：
F36中，
G36为测量噪声输入矩阵为，u为测量噪声，/>其中，ug为陀螺的测量噪声，ua为加速度计的测量噪声；
S303、构建三十六维卡尔曼滤波器滤波器的观测方程：
Z＝H36X36+v，
式中，Z为观测量，以速度为观测量且在标定时双轴旋转惯导保持静止状态，则Z＝[VEVN VU]T；H36为观测矩阵，H36＝[03×3 I3 03×30]，I3为三行三列的单位矩阵，即03×30为三行三十列的零矩阵，03×3为三行三列的零矩阵；v是观测噪声；
S304、确定卡尔曼滤波估计结果的反馈补偿形式：
式中，为/>的反对称矩阵，/>为双轴旋转惯导的/>的解算结果；/>为双轴旋转惯导的Vn的解算结果；/>为双轴旋转惯导的L的解算结果；/>为双轴旋转惯导的λ的解算结果；/>为双轴旋转惯导的h的解算结果；，/>为卡尔曼滤波估计结果反馈补偿前的KG；/>为卡尔曼滤波估计结果反馈补偿前的εm；，/>为卡尔曼滤波估计结果反馈补偿前的KA；/>为卡尔曼滤波估计结果反馈补偿前的/> 为卡尔曼滤波估计结果反馈补偿前的KA2；/>为卡尔曼滤波估计结果反馈补偿前的rx；/>为卡尔曼滤波估计结果反馈补偿前的ry；/>为卡尔曼滤波估计结果反馈补偿前的rz。
5.根据权利要求1所述的双轴旋转惯导动态误差抑制方法，其特征在于，在步骤S4中，双轴转位机构的转位路径具体设计为：
基于次序0对应初始位置：双轴旋转惯导以其Xm轴朝东向，Ym轴朝北向，Zm轴朝天向的姿态，静止10分钟；
次序1：以中框轴为旋转轴，顺时针旋转90°；该状态下，双轴旋转惯导的Xm轴朝东向，Ym轴朝天向，Zm轴朝南向；转位完成后，静止180s；
次序2：以中框轴为旋转轴，顺时针旋转180°；该状态下，双轴旋转惯导的Xm轴朝东向，Ym轴朝地向，Zm轴朝北向；转位完成后，静止180s；
次序3：以中框轴为旋转轴，顺时针旋转180°；该状态下，双轴旋转惯导的Xm轴朝东向，Ym轴朝天向，Zm轴朝南向；转位完成后，静止180s；
次序4：以内框轴为旋转轴，顺时针旋转90°；该状态下，双轴旋转惯导的Xm轴朝天向，Ym轴朝西向，Zm轴朝南向；转位完成后，静止180s；
次序5：以内框轴为旋转轴，顺时针旋转180°；该状态下，双轴旋转惯导的Xm轴朝地向，Ym轴朝东向，Zm轴朝南向；转位完成后，静止180s；
次序6：以内框轴为旋转轴，顺时针旋转180°；该状态下，双轴旋转惯导的Xm轴朝天向，Ym轴朝西向，Zm轴朝南向；转位完成后，静止180s；
次序7：以中框轴为旋转轴，顺时针旋转90°；该状态下，双轴旋转惯导的Xm轴朝南向，Ym轴朝西向，Zm轴朝地向；转位完成后，静止180s；
次序8：以中框轴为旋转轴，顺时针旋转180°；该状态下，双轴旋转惯导的Xm轴朝北向，Ym轴朝西向，Zm轴朝天向；转位完成后，静止180s；
次序9：以中框轴为旋转轴，顺时针旋转180°；该状态下，双轴旋转惯导的Xm轴朝南向，Ym轴朝西向，Zm轴朝地向；转位完成后，静止180s；
次序10：以中框轴为旋转轴，顺时针旋转90°；该状态下，双轴旋转惯导的Xm轴朝地向，Ym轴朝西向，Zm轴朝北向；转位完成后，静止180s；
次序11：以中框轴为旋转轴，顺时针旋转90°；该状态下，双轴旋转惯导的Xm轴朝北向，Ym轴朝西向，Zm轴朝天向；转位完成后，静止180s；
次序12：以中框轴为旋转轴，顺时针旋转90°；该状态下，双轴旋转惯导的Xm轴朝天向，Ym轴朝西向，Zm轴朝南向；转位完成后，静止180s；
次序13：以内框轴为旋转轴，顺时针旋转90°；该状态下，双轴旋转惯导的Xm轴朝西向，Ym轴朝地向，Zm轴朝南向；转位完成后，静止180s；
次序14：以内框轴为旋转轴，顺时针旋转90°；该状态下，双轴旋转惯导的Xm轴朝地向，Ym轴朝东向，Zm轴朝南向；转位完成后，静止180s；
次序15：以内框轴为旋转轴，顺时针旋转90°；该状态下，双轴旋转惯导的Xm轴朝东向，Ym轴朝天向，Zm轴朝南向；转位完成后，静止180s；
次序16：以中框轴为旋转轴，顺时针旋转90°；该状态下，双轴旋转惯导的Xm轴朝东向，Ym轴朝南向，Zm轴朝地向；转位完成后，静止180s；
次序17：以中框轴为旋转轴，顺时针旋转90°；该状态下，双轴旋转惯导的Xm轴朝东向，Ym轴朝地向，Zm轴朝北向；转位完成后，静止180s；
次序18：以中框轴为旋转轴，顺时针旋转90°；该状态下，双轴旋转惯导的Xm轴朝东向，Ym轴朝北向，Zm轴朝天向；转位完成后，静止180s。
6.根据权利要求5所述的双轴旋转惯导动态误差抑制方法，其特征在于，在步骤S4中，一次标定实验中按照设计的转位路径进行两次双轴转位机构的转位。
7.根据权利要求3所述的双轴旋转惯导动态误差抑制方法，其特征在于，步骤S5的具体实施步骤为：
S501、对三十六维状态量进行初值赋值：和δV的初值分别由步骤S4中双轴旋转的初对准结果中的初始姿态角和初始速度进行赋值；δP的初值由步骤S4中双轴旋转惯导完成初对准后输出的初始位置减去实际实验位置得到的差值进行赋值；Xg、Xa、XKa2和Xr的初值根据双轴旋转惯导中陀螺和加速度计的仪器实际情况进行赋值；
S502、基于步骤S501确定的三十六维状态量初值作为卡尔曼滤波的状态量的初值，将步骤S301构建的三十六维状态量分别带入步骤S302构建的三十六维卡尔曼滤波器状态方程和步骤S303构建的三十六维卡尔曼滤波器观测方程中，采用卡尔曼滤波算法和步骤S3构建的卡尔曼滤波器处理标定数据，通过对卡尔曼滤波器中的状态量从初值开始不断迭代和反馈补偿，获得卡尔曼滤波器中的状态量估计值；其中，状态量估计值中关于惯性器件误差的状态量分量即准确的标定参数，包括：X向陀螺标度因数Y向陀螺标度因数/>和Z向陀螺标度因数/>陀螺安装误差角γyz、γzy和γzx；X向陀螺零偏/>Y向陀螺零偏/>和Z向陀螺零偏/>X向加速度计标度因数/>Y向加速度计标度因数/>和Z向加速度计标度因数/>的误差；加速度计安装误差角αxz、αxy、αyz、αyx、αzy和αzx；X向加速度计的二阶非线性系数Kax2、Y向加速度计的二阶非线性系数Kay2和Z向加速度计的二阶非线性系数Kaz2；X向加速度计零偏/>Y向加速度计零偏/>和Z向加速度计零偏/>X向加速度计内杆臂rx、Y向加度计内杆臂ry和Z向加速度计内杆臂rz。
8.根据权利要求7所述的双轴旋转惯导动态误差抑制方法，其特征在于，步骤S6的具体实施步骤为：
S601、采用现有的六十四次序旋转调制方案作为双轴旋转惯导导航过程中的转位机构旋转方案，并在每转动至一个位置后停止10s时间，获得旋转调制后的惯性器件输出结果，包括陀螺输出的角速率和加速度计输出的比力/>
S602、将由步骤S5获得的标定参数分别代入至陀螺标定模型和加速度计标定模型中，以获得陀螺标定后的输出和加速度计标定后的输出；其中，
(1)陀螺的标定模型为：
式中，
(2)加速度计的标定模型为：
式中，
S603、利用经过步骤S602得到的陀螺标定后的输出和加速度计标定后的输出，利用X向加速度计内杆臂rx、Y向加速度计内杆臂ry和Z向加速度计内杆臂rz对加速度数据进行实时补偿，其具体公式为：
式中，axb、ayb和azb为加速度计内杆臂补偿后的加速度实时结果在m系中Xm轴、Ym轴和Zm轴的分量；axc、ayc和azc为加速度计内杆臂补偿前的加速度实时结果在m系中Xm轴、Ym轴和Zm轴的分量；ωxc、ωyc和ωzc为角速率在m系中Xm轴、Ym轴和Zm轴的分量；
S604、将步骤S603中的axb、ayb、azb和ωxc、ωyc、ωzc带入惯导解算方程，经解算得到的双轴旋转惯导的导航结果，即为对双轴旋转惯导动态误差进行实时补偿后的导航结果。
CN202210835839.1A 2022-07-15 2022-07-15 一种双轴旋转惯导动态误差抑制方法 Active CN115265590B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210835839.1A CN115265590B (zh) 2022-07-15 2022-07-15 一种双轴旋转惯导动态误差抑制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210835839.1A CN115265590B (zh) 2022-07-15 2022-07-15 一种双轴旋转惯导动态误差抑制方法

Publications (2)

Publication Number Publication Date
CN115265590A CN115265590A (zh) 2022-11-01
CN115265590B true CN115265590B (zh) 2024-04-30

Family

ID=83765352

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210835839.1A Active CN115265590B (zh) 2022-07-15 2022-07-15 一种双轴旋转惯导动态误差抑制方法

Country Status (1)

Country Link
CN (1) CN115265590B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116481564B (zh) * 2023-03-11 2024-02-23 中国人民解放军国防科技大学 基于Psi角误差修正模型的极地双惯导协同标定方法
CN116242397B (zh) * 2023-03-11 2024-01-30 中国人民解放军国防科技大学 一种速度误差修正模型下的双惯导协同标定方法
CN117782001B (zh) * 2024-02-28 2024-05-07 爱瑞克(大连)安全技术集团有限公司 一种papi助航灯动态角度监测预警方法及***

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109029500A (zh) * 2018-07-24 2018-12-18 中国航空工业集团公司西安飞行自动控制研究所 一种双轴旋转调制***全参数自标定方法
CN113029199A (zh) * 2021-03-15 2021-06-25 中国人民解放军国防科技大学 一种激光陀螺惯导***的***级温度误差补偿方法
WO2021227011A1 (zh) * 2020-05-11 2021-11-18 中国科学院地质与地球物理研究所 一种陀螺随钻测量***及方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110487301B (zh) * 2019-09-18 2021-07-06 哈尔滨工程大学 一种雷达辅助机载捷联惯性导航***初始对准方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109029500A (zh) * 2018-07-24 2018-12-18 中国航空工业集团公司西安飞行自动控制研究所 一种双轴旋转调制***全参数自标定方法
WO2021227011A1 (zh) * 2020-05-11 2021-11-18 中国科学院地质与地球物理研究所 一种陀螺随钻测量***及方法
CN113029199A (zh) * 2021-03-15 2021-06-25 中国人民解放军国防科技大学 一种激光陀螺惯导***的***级温度误差补偿方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于高阶卡尔曼滤波的激光捷联惯导***级标定方法;刘冰;任继山;白焕旭;王盛;陈鸿跃;;导弹与航天运载技术;20170810(第04期);全文 *
激光陀螺捷联惯导***多位置***级标定方法;石文峰;王省书;郑佳兴;战德军;王以忠;;红外与激光工程;20161125(第11期);全文 *

Also Published As

Publication number Publication date
CN115265590A (zh) 2022-11-01

Similar Documents

Publication Publication Date Title
CN115265590B (zh) 一种双轴旋转惯导动态误差抑制方法
CN108318052B (zh) 一种基于双轴连续旋转的混合式平台惯导***标定方法
CN113029199B (zh) 一种激光陀螺惯导***的***级温度误差补偿方法
CN111678538B (zh) 一种基于速度匹配的动态水平仪误差补偿方法
CN106969783B (zh) 一种基于光纤陀螺惯性导航的单轴旋转快速标定技术
CN110160554B (zh) 一种基于寻优法的单轴旋转捷联惯导***标定方法
CN103630146B (zh) 一种离散解析与Kalman滤波结合的激光陀螺IMU标定方法
CN111561948B (zh) 四轴冗余捷联惯导的***级标定方法
CN112595350B (zh) 一种惯导***自动标定方法及终端
CN104181574A (zh) 一种捷联惯导***/全球导航卫星***组合导航滤波***及方法
CN108132060B (zh) 一种捷联惯导***无基准的***级标定方法
CN113503894B (zh) 基于陀螺基准坐标系的惯导***误差标定方法
GB2378765A (en) Error compensation in an inertial navigation system
CN112710328B (zh) 一种四轴冗余惯导***的误差标定方法
CN115143993A (zh) 基于三轴转台的激光陀螺惯导***g敏感性误差标定方法
CN107246883A (zh) 一种高精度星敏感器安装矩阵在轨实时校准方法
CN110749338A (zh) 一种惯性测量单元偏轴-旋转复合转位误差标定方法
CN114877915B (zh) 一种激光陀螺惯性测量组件g敏感性误差标定装置及方法
CN112697143B (zh) 高精度载体动态姿态测量方法及***
CN110940357B (zh) 一种用于旋转惯导单轴自对准的内杆臂标定方法
CN115683155A (zh) 一种旋转捷联惯导***误差***级标定方法
CN111089606B (zh) 一种三自激光惯组关键参数快速自标定方法
Liang et al. A novel calibration method between two marine rotational inertial navigation systems based on state constraint Kalman filter
CN111964671B (zh) 一种基于双轴旋转调制的惯性天文组合导航***及方法
CN110220534B (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