CN101718560B - 基于单轴四位置转停方案的捷联系统的误差抑制方法 - Google Patents

Abstract

本发明提供的是一种基于单轴四位置转停方案的捷联系统误差抑制方法。(1)通过GPS确定载体的初始位置参数；(2)采集光纤陀螺仪输出和加速度计输出的数据并对数据进行处理；(3)惯性测量单元绕着载体方位轴固定的四个位置正反转停；(4)将IMU旋转后光纤陀螺仪和石英加速度计生成的数据转换到导航坐标系下，得到惯性器件常值偏差的调制形式；(5)利用光纤陀螺的输出值ω_is ^s对捷联矩阵T_s ⁿ进行更新；(6)计算IMU旋转调制后载体的速度和位置；本发明水平方向上的惯性器件常值偏差进行调制，提高导航定位精度。

Description

基于单轴四位置转停方案的捷联系统的误差抑制方法

技术领域

本发明涉及的是一种捷联系统误差抑制方法，尤其涉及的是一种基于单轴四位置转停方案的捷联系统误差抑制方法。

背景技术

在捷联式惯导系统中，所有的惯性元件都直接安装在载体上，惯性元件的输出就是载体相对惯性空间的加速度和角速度，有计算机将载体坐标系下测得的数据变换到导航坐标系再进行导航解算，这相当于在计算机内虚拟了一个平台作为导航计算的参考。由于取消了平台结构，因而系统的体积大大缩小，成本和维护费也大大降低。

惯性敏感元件的误差是惯性导航系统误差的主要决定因素。从工艺上提高惯性敏感元件的精度，存在技术难度大，周期长的问题，而且也不能超越现阶段实现的可能。惯性元件误差补偿方法有两种：一是采用高精度的漂移误差补偿软件，另一种是采用系统的反转或者是旋转补偿法。

自从光学陀螺出现以后，在旋转自动补偿技术和捷联惯性导航技术的基础上出现了一种新型的惯导系统，即旋转式惯性导航系统。它具有与传统平台系统类似的框架和转轴，然而却没有稳定的平台。简而言之，此类惯导系统相当于在捷联惯导系统的外贸加上转动机构和测角装置，导航解算也采用了捷联惯导算法，这样导航计算出来的依然是载体的位置和速度信息，而导航直接计算出的姿态信息只是惯性测量单元(IMU)的姿态，这时候需要加上IMU相对载体的转动角度，就得到了载体的姿态信息。

旋转式惯性导航系统中转动机构的目的是为了在一个转动周期内让惯性元件的常值漂移对导航所产生的误差能够自动抵消，相应的技术也被称为系统误差平均技术，被广泛的应用到当前的船用光学陀螺惯性导航系统中。

CNKI库中与本发明相关的公开报道有：1、《光学陀螺旋转惯导系统原理探讨》(袁保伦、饶谷音；国防科技大学学报；2006年第6期)，该文章主要以静电陀螺的旋转调制方式为例讨论了单轴、双轴旋转调制方法，但是没有涉及到本文提到的单轴旋转调制方法。2、《旋转IMU在光纤捷联航姿系统中的应用》(王其、徐晓苏；中国惯性技术学报；2007年第3期)，该文章主要介绍了单轴、双轴旋转方式，并在理论上给予证明。3、《捷联惯导动基座对准新方法及导航误差抑制技术研究》(练军想，国防科学技术大学博士论文，出版日2009年7月)，该论文主要介绍了单轴旋转调制技术的误差补偿原理，但是文章中提到的是单轴连续旋转，没有提及本发明提出的单轴四位置转停方案。

发明内容

本发明的目的在于提供一种能够提高导航定位精度的基于单轴四位置转停方案的捷联系统误差抑制方法。

本发明的目的是这样实现的：包括以下步骤：

(1)利用全球定位系统GPS确定载体的初始位置参数，将它们装订至导航计算机中；

(2)光纤陀螺捷联惯性导航系统进行预热后采集光纤陀螺仪和石英加速度计输出的数据；

(3)IMU采用8个转停次序为一个旋转周期的转位方案；

(4)将惯性测量单元旋转后陀螺仪生成的数据转换到导航坐标系下，得到惯性器件常值偏差的调制形式；

(5)将IMU坐标系下陀螺仪的输出值

带入惯导系统中采用四元数法对捷联矩阵

进行更新：

${\mathrm{\ω}}_{\mathrm{ns}}^s={\mathrm{\ω}}_{\mathrm{is}}^s-{\left(T_s^n\right)}^T({\mathrm{\ω}}_{\mathrm{ie}}^n+{\mathrm{\ω}}_{\mathrm{en}}^n)$

其中：为地球自转角速度在导航系下的分量；

为导航坐标系相对地球坐标系的运动角速度在导航系下的分量；

为IMU相对导航坐标系的运动角速度在载体坐标系上的分量；

(6)利用石英加速度计的输出值和步骤(5)计算的姿态矩阵计算出经过IMU旋转调制后载体的速度和位置。

本发明还可以包括：

1、所述的IMU采用8个转停次序为一个旋转周期的转位方案为：

次序1，IMU从A点出发顺时针转动180°到达位置C，停止时间T_t；次序2，IMU从C点出发顺时针转动90°到达位置D，停止时间T_t；次序3，IMU从D点出发逆时针转动180°到达位置B，停止时间T_t；次序4，IMU从B点出发逆时针转动90°到达位置A，停止时间T_t；次序5，IMU从A点出发逆时针转动180°到达位置C，停止时间T_t；次序6，IMU从C点出发逆时针转动90°到达位置B，停止时间T_t；次序7，IMU从B点出发顺时针转动180°到达位置D，停止时间T_t；次序8，IMU从D点出发顺时针转动90°到达位置A，停止时间T_t；IMU按照此转动顺序循环进行；水平东向轴上的IMU停顿点p3、p8与p4、p7对称于转轴中心；北向轴上的停顿点p1、p5与p2、p6对称于转轴中心；四位置转停方案仍然是转动角度为180°或90°间隔进行。

2、所述的将惯性测量单元旋转后陀螺仪生成的数据转换到导航坐标系下，得到惯性器件常值偏差的调制形式包括：

IMU水平方向上的陀螺常值漂移分别为ε_x和ε_y，载体静止条件下，IMU停顿的A、B、C、D四个位置相对于导航坐标系对称，在一个循环周期内四个固定位置上器件偏差产生的水平误差角在导航系下满足：

$\left\{\begin{array}{c}2{\left({\∫}_0^{T_t}{\mathrm{\ϵ}}_E^n\mathrm{dt}\right)}_A+2{\left({\∫}_0^{T_t}{\mathrm{\ϵ}}_E^n\mathrm{dt}\right)}_B+2{\left({\∫}_0^{T_t}{\mathrm{\ϵ}}_E^n\mathrm{dt}\right)}_C+2{\left({\∫}_0^{T_t}{\mathrm{\ϵ}}_E^n\mathrm{dt}\right)}_D=0\\ 2{\left({\∫}_0^{T_t}{\mathrm{\ϵ}}_N^n\mathrm{dt}\right)}_A+2{\left({\∫}_0^{T_t}{\mathrm{\ϵ}}_N^n\mathrm{dt}\right)}_B+2{\left({\∫}_0^{T_t}{\mathrm{\ϵ}}_N^n\mathrm{dt}\right)}_C+2{\left({\∫}_0^{T_t}{\mathrm{\ϵ}}_N^n\mathrm{dt}\right)}_D=0\end{array}\right.$

IMU正反转动的8个环节内，导航系水平方向上的误差角度分别表示为：

将一个完整转动周期内的IMU停顿位置和转动过程的东向和北向的水平误差角分别相加，得到导航系下累计的水平误差角度，

$\left\{\begin{array}{c}\∫{\mathrm{\ϵ}}_E^n\mathrm{dt}=0\\ \∫{\mathrm{\ϵ}}_N^n\mathrm{dt}=0\end{array}\right.$

导航坐标系上水平方向上陀螺仪常值偏差得到调制，该常值偏差对系统的作用为零，水平方向上加速度计零偏也得到类似的调制。

3、所述的将IMU坐标系下陀螺仪的输出值

带入惯导系统中采用四元数法对捷联矩阵

进行更新包括：

更新四元数和姿态矩阵：

设载体坐标系相对导航坐标系的转动四元数为：

Q＝q₀+q₁i_s+q₂j_s+q₃k_s

其中：i_s、j_s、k_s分别表示IMU坐标系ox_s轴、oy_s轴、oz_s轴上的单位方向向量；

四元数的即时修正通过解四元数微分方程

来实现：

$\left[\begin{array}{c}{\stackrel{\·}{q}}_0\\ {\stackrel{\·}{q}}_1\\ {\stackrel{\·}{q}}_2\\ {\stackrel{\·}{q}}_3\end{array}\right]=\frac{1}{2}\left[\begin{array}{cccc}0& {-\mathrm{\ω}}_{\mathrm{ns}}^{\mathrm{sx}}& -{\mathrm{\ω}}_{\mathrm{ns}}^{\mathrm{sy}}& -{\mathrm{\ω}}_{\mathrm{ns}}^{\mathrm{sz}}\\ {\mathrm{\ω}}_{\mathrm{ns}}^{\mathrm{sx}}& 0& {\mathrm{\ω}}_{\mathrm{ns}}^{\mathrm{sz}}& -{\mathrm{\ω}}_{\mathrm{ns}}^{\mathrm{sy}}\\ {\mathrm{\ω}}_{\mathrm{ns}}^{\mathrm{sy}}& {-\mathrm{\ω}}_{\mathrm{ns}}^{\mathrm{sz}}& 0& {\mathrm{\ω}}_{\mathrm{ns}}^{\mathrm{sx}}\\ {\mathrm{\ω}}_{\mathrm{ns}}^{\mathrm{sz}}& {\mathrm{\ω}}_{\mathrm{ns}}^{\mathrm{sy}}& -{\mathrm{\ω}}_{\mathrm{ns}}^{\mathrm{sx}}& 0\end{array}\right]\left[\begin{array}{c}{q}_0\\ q_1\\ q_2\\ q_3\end{array}\right]$

其中：

分别表示载体相对导航系的运动角速度在IMU坐标系ox_s轴、oy_s轴、oz_s轴上的分量；

姿态矩阵

的更新过程为：

$T_s^n=\left[\begin{array}{ccc}{q}_0^2+q_1^2-q_2^2-q_3^2& 2(q_1{q}_2-q_0{q}_3)& 2(q_1{q}_3+q_0{q}_2)\\ 2(q_1{q}_2+q_0{q}_3)& q_0^2-q_1^2+q_2^2-q_3^2& 2(q_2{q}_3-q_0{q}_1)\\ 2(q_1{q}_3-q_0{q}_2)& 2(q_2{q}_3+q_0{q}_1)& q_0^2-q_1^2-q_2^2+q_3^2\end{array}\right]$

4、所述的计算出经过IMU旋转调制后载体的速度和位置地方法为：

1)计算导航系下加速度：

$\left[\begin{array}{c}{f}^{\mathrm{nx}}\\ f^{\mathrm{ny}}\\ f^{\mathrm{nz}}\end{array}\right]=T_s^n\left[\begin{array}{c}{f}_{\mathrm{is}}^{\mathrm{sx}}\\ f_{\mathrm{is}}^{\mathrm{sy}}\\ f_{\mathrm{is}}^{\mathrm{sz}}\end{array}\right]$

2)计算载体的水平速度和位置：

根据t₁时刻的载体东向水平速度V_x(t₁)和北向水平速度V_y(t₁)，求取t₁时刻载体水平速度的变化率为：

$\left\{\begin{array}{c}{\stackrel{\·}{V}}_x\left(t_1\right)=f^{\mathrm{nx}}+({2\mathrm{\ω}}_{\mathrm{ie}}^{\mathrm{nz}}+{\mathrm{\ω}}_{\mathrm{en}}^{\mathrm{nz}})V_y\left(t_1\right)\\ {\stackrel{\·}{V}}_y\left(t_1\right)=f^{\mathrm{ny}}-({2\mathrm{\ω}}_{\mathrm{ie}}^{\mathrm{nz}}+{\mathrm{\ω}}_{\mathrm{en}}^{\mathrm{nz}})V_x\left(t_1\right)\end{array}\right.$

在t₂时刻水平速度和载体位置分别为：

$\left\{\begin{array}{c}{V}_x\left(t_2\right)=V_x\left(t_1\right)+{\stackrel{\·}{V}}_x\left(t_1\right)(t_2-t_1)\\ V_y\left(t_2\right)=V_y\left(t_1\right)+{\stackrel{\·}{V}}_y\left(t_1\right)(t_2-t_1)\end{array}\right.$

3)计算载体速度误差和位置误差：

$\left\{\begin{array}{c}{\mathrm{\ΔV}}_x=V_x\left(t_2\right)-V_{x0}\\ {\mathrm{\ΔV}}_y=V_y\left(t_2\right)-V_{y0}\end{array}\right.$

其中：V_x0、V_y0分别表示初始时刻载体的东向和北向速度；ΔV_x、ΔV_y分别表示载体东向、北向速度的变化量；

λ₀分别表示初始时刻载体所处位置的经度和纬度；Δλ分别表示载体的纬度、经度的变化量；R_xp、R_yp分别表示地球子午圈、卯酉圈的曲率半径；t₁、t₂为惯导系统的解算过程中两个相邻的时间点。

本发明与现有技术相比的优点在于：本发明打破了传统捷联惯导系统中IMU与载体固连导致系统导航精度受到惯性器件偏差影响的约束，提出一种将IMU绕载体方位轴固定的四个位置正反转停的惯性器件常值偏差调制方案，该方法可以将水平方向上惯性器件常值偏差进行调制，有效地提高导航定位精度。

对本发明有益的效果说明如下：

在VC++仿真条件下，对该方法进行仿真实验：

载体处于静止状态，IMU四位置转停方案的误差模型参数：

四个位置的停顿时间：T_t＝5min；

转动180°和90°时消耗的时间：T_z＝12s；

转动180°和90°的过程中，每一个转位中的加减速时间各为4s；

载体初始位置：北纬45.7796°，东经126.6705°；

初始姿态误差角：三个初始姿态误差角均为零；

赤道半径：R_e＝6378393.0m；

椭球度：e＝3.367e-3；

由万有引力可得的地球表面重力加速度：g₀＝9.78049；

地球自转角速度(弧度/秒)：7.2921158e-5；

陀螺仪常值漂移：0.01度/小时；

加速度计零偏：10^-4g₀；

常数：π＝3.1415926；

利用发明所述方法得到载体速度误差曲线和位置误差曲线分别如图3、图4所示。结果表明IMU四位置转停条件下，采用本发明方法可以获得较高的定位精度。

附图说明

图1为本发明的一种基于单轴四位置转停方案的捷联系统误差抑制方法流程图；

图2为本发明的IMU单轴四位置转停方案图；

图3为本发明的基于IMU单轴四位置转停方案的陀螺仪常值漂移在导航系下的调制曲线；

图4为本发明的基于IMU单轴四位置转停方案的载体速度误差与IMU静止状态时载体速度误差的对比实验曲线；

图5为本发明的基于IMU单轴四位置转停方案的载体位置误差与IMU静止状态时载体定位误差的对比实验曲线。

具体实施方式

下面结合附图举例对本发明做更详细地描述：

(1)利用全球定位系统GPS确定载体的初始位置参数，将它们装订至导航计算机中；

(2)光纤陀螺捷联惯性导航系统进行预热后采集光纤陀螺仪和石英加速度计输出的数据；

(3)IMU采用8个转停次序为一个旋转周期的转位方案；

次序1，IMU从A点出发顺时针转动180°到达位置C，停止时间T_t；次序2，IMU从C点出发顺时针转动90°到达位置D，停止时间T_t；次序3，IMU从D点出发逆时针转动180°到达位置B，停止时间T_t；次序4，IMU从B点出发逆时针转动90°到达位置A，停止时间T_t；次序5，IMU从A点出发逆时针转动180°到达位置C，停止时间T_t；次序6，IMU从C点出发逆时针转动90°到达位置B，停止时间T_t；次序7，IMU从B点出发顺时针转动180°到达位置D，停止时间T_t；次序8，IMU从D点出发顺时针转动90°到达位置A，停止时间T_t；IMU按照此转动顺序循环进行。为了有效地对水平方向上的惯性器件偏差在对称位置上进行正负平均，定义水平东向轴上的IMU停顿点p3、p8与p4、p7对称于转轴中心；北向轴上的停顿点p1、p5与p2、p6对称于转轴中心。改进的四位置转停方案仍然是转动角度为180°或90°间隔进行。

(4)将惯性测量单元旋转后陀螺仪生成的数据转换到导航坐标系下，得到惯性器件常值偏差的调制形式；

假定IMU水平方向上的陀螺常值漂移分别为ε_x和ε_y。载体静止条件下，由于IMU停顿的A、B、C、D四个位置相对于导航坐标系对称，因此在一个循环周期内四个固定位置上器件偏差产生的水平误差角在导航系下必然满足：

$\left\{\begin{array}{c}2{\left({\∫}_0^{T_t}{\mathrm{\ϵ}}_E^n\mathrm{dt}\right)}_A+2{\left({\∫}_0^{T_t}{\mathrm{\ϵ}}_E^n\mathrm{dt}\right)}_B+2{\left({\∫}_0^{T_t}{\mathrm{\ϵ}}_E^n\mathrm{dt}\right)}_C+2{\left({\∫}_0^{T_t}{\mathrm{\ϵ}}_E^n\mathrm{dt}\right)}_D=0\\ 2{\left({\∫}_0^{T_t}{\mathrm{\ϵ}}_N^n\mathrm{dt}\right)}_A+2{\left({\∫}_0^{T_t}{\mathrm{\ϵ}}_N^n\mathrm{dt}\right)}_B+2{\left({\∫}_0^{T_t}{\mathrm{\ϵ}}_N^n\mathrm{dt}\right)}_C+2{\left({\∫}_0^{T_t}{\mathrm{\ϵ}}_N^n\mathrm{dt}\right)}_D=0\end{array}\right.---\left(1\right)$

IMU正反转动的8个环节内，导航系水平方向上的误差角度可以分别表示为：

将一个完整转动周期内的IMU停顿位置和转动过程的东向和北向的水平误差角分别相加，得到导航系下累计的水平误差角度。

$\left\{\begin{array}{c}\∫{\mathrm{\ϵ}}_E^n\mathrm{dt}=0\\ \∫{\mathrm{\ϵ}}_N^n\mathrm{dt}=0\end{array}\right.---\left(4\right)$

至此，导航坐标系上水平方向上陀螺仪常值偏差得到调制，该常值偏差对系统的作用为零，水平方向上加速度计零偏也得到类似的调制。

(5)将IMU坐标系下陀螺仪的输出值

带入惯导系统中采用四元数法对捷联矩阵

进行更新：

${\mathrm{\ω}}_{\mathrm{ns}}^s={\mathrm{\ω}}_{\mathrm{is}}^s-{\left(T_s^n\right)}^T({\mathrm{\ω}}_{\mathrm{ie}}^n+{\mathrm{\ω}}_{\mathrm{en}}^n)---\left(5\right)$

其中：

为地球自转角速度在导航系下的分量；

为导航坐标系相对地球坐标系的运动角速度在导航系下的分量；

为IMU相对导航坐标系的运动角速度在载体坐标系上的分量。

更新四元数和姿态矩阵：

设载体坐标系相对导航坐标系的转动四元数为：

Q＝q₀+q₁i_s+q₂j_s+q₃k_s    (6)

其中：i_s、j_s、k_s分别表示IMU坐标系ox_s轴、oy_s轴、oz_s轴上的单位方向向量。

四元数的即时修正可以通过解四元数微分方程

来实现：

$\left[\begin{array}{c}{\stackrel{\·}{q}}_0\\ {\stackrel{\·}{q}}_1\\ {\stackrel{\·}{q}}_2\\ {\stackrel{\·}{q}}_3\end{array}\right]=\frac{1}{2}\left[\begin{array}{cccc}0& {-\mathrm{\ω}}_{\mathrm{ns}}^{\mathrm{sx}}& -{\mathrm{\ω}}_{\mathrm{ns}}^{\mathrm{sy}}& -{\mathrm{\ω}}_{\mathrm{ns}}^{\mathrm{sz}}\\ {\mathrm{\ω}}_{\mathrm{ns}}^{\mathrm{sx}}& 0& {\mathrm{\ω}}_{\mathrm{ns}}^{\mathrm{sz}}& -{\mathrm{\ω}}_{\mathrm{ns}}^{\mathrm{sy}}\\ {\mathrm{\ω}}_{\mathrm{ns}}^{\mathrm{sy}}& {-\mathrm{\ω}}_{\mathrm{ns}}^{\mathrm{sz}}& 0& {\mathrm{\ω}}_{\mathrm{ns}}^{\mathrm{sx}}\\ {\mathrm{\ω}}_{\mathrm{ns}}^{\mathrm{sz}}& {\mathrm{\ω}}_{\mathrm{ns}}^{\mathrm{sy}}& -{\mathrm{\ω}}_{\mathrm{ns}}^{\mathrm{sx}}& 0\end{array}\right]\left[\begin{array}{c}{q}_0\\ q_1\\ q_2\\ q_3\end{array}\right]---\left(7\right)$

其中：

分别表示载体相对导航系的运动角速度在IMU坐标系ox_s轴、oy_s轴、oz_s轴上的分量。

姿态矩阵

的更新过程为：

$T_s^n=\left[\begin{array}{ccc}{q}_0^2+q_1^2-q_2^2-q_3^2& 2(q_1{q}_2-q_0{q}_3)& 2(q_1{q}_3+q_0{q}_2)\\ 2(q_1{q}_2+q_0{q}_3)& q_0^2-q_1^2+q_2^2-q_3^2& 2(q_2{q}_3-q_0{q}_1)\\ 2(q_1{q}_3-q_0{q}_2)& 2(q_2{q}_3+q_0{q}_1)& q_0^2-q_1^2-q_2^2+q_3^2\end{array}\right]---\left(8\right)$

(6)利用石英加速度计的输出值

和步骤(5)计算的姿态矩阵计算出经过IMU旋转调制后载体的速度和位置。

1)计算导航系下加速度：

$\left[\begin{array}{c}{f}^{\mathrm{nx}}\\ f^{\mathrm{ny}}\\ f^{\mathrm{nz}}\end{array}\right]=T_s^n\left[\begin{array}{c}{f}_{\mathrm{is}}^{\mathrm{sx}}\\ f_{\mathrm{is}}^{\mathrm{sy}}\\ f_{\mathrm{is}}^{\mathrm{sz}}\end{array}\right]---\left(9\right)$

2)计算载体的水平速度和位置：

根据t₁时刻的载体东向水平速度V_x(t₁)和北向水平速度V_y(t₁)，求取t₁时刻载体水平速度的变化率为：

$\left\{\begin{array}{c}{\stackrel{\·}{V}}_x\left(t_1\right)=f^{\mathrm{nx}}+({2\mathrm{\ω}}_{\mathrm{ie}}^{\mathrm{nz}}+{\mathrm{\ω}}_{\mathrm{en}}^{\mathrm{nz}})V_y\left(t_1\right)\\ {\stackrel{\·}{V}}_y\left(t_1\right)=f^{\mathrm{ny}}-({2\mathrm{\ω}}_{\mathrm{ie}}^{\mathrm{nz}}+{\mathrm{\ω}}_{\mathrm{en}}^{\mathrm{nz}})V_x\left(t_1\right)\end{array}\right.---\left(10\right)$

在t₂时刻水平速度和载体位置分别为：

$\left\{\begin{array}{c}{V}_x\left(t_2\right)=V_x\left(t_1\right)+{\stackrel{\·}{V}}_x\left(t_1\right)(t_2-t_1)\\ V_y\left(t_2\right)=V_y\left(t_1\right)+{\stackrel{\·}{V}}_y\left(t_1\right)(t_2-t_1)\end{array}\right.---\left(11\right)$

3)计算载体速度误差和位置误差：

$\left\{\begin{array}{c}{\mathrm{\ΔV}}_x=V_x\left(t_2\right)-V_{x0}\\ {\mathrm{\ΔV}}_y=V_y\left(t_2\right)-V_{y0}\end{array}\right.---\left(13\right)$

其中：V_x0、V_y0分别表示初始时刻载体的东向和北向速度；ΔV_x、ΔV_y分别表示载体东向、北向速度的变化量；

λ₀分别表示初始时刻载体所处位置的经度和纬度；Δλ分别表示载体的纬度、经度的变化量；R_xp、R_yp分别表示地球子午圈、卯酉圈的曲率半径；t₁、t₂为惯导系统的解算过程中两个相邻的时间点。

Claims (3)

1.一种基于单轴四位置转停方案的捷联系统的误差抑制方法，其特征在于包括以下步骤：

(1)利用全球定位系统GPS确定载体的初始位置参数，将它们装订至导航计算机中；

(2)捷联系统进行预热后采集光纤陀螺仪和石英加速度计输出的数据；

(3)IMU采用8个转停次序为一个旋转周期的转位方案，方案为：

次序1，IMU从A点出发顺时针转动180°到达位置C，停止时间T_t；次序2，IMU从C点出发顺时针转动90°到达位置D，停止时间T_t；次序3，IMU从D点出发逆时针转动180°到达位置B，停止时间T_t；次序4，IMU从B点出发逆时针转动90°到达位置A，停止时间T_t；次序5，IMU从A点出发逆时针转动180°到达位置C，停止时间T_t；次序6，IMU从C点出发逆时针转动90°到达位置B，停止时间T_t；次序7，IMU从B点出发顺时针转动180°到达位置D，停止时间T_t；次序8，IMU从D点出发顺时针转动90°到达位置A，停止时间T_t；IMU按照此转动顺序循环进行；水平东向轴上的IMU停顿点p3、p8与p4、p7对称于转轴中心；北向轴上的停顿点p1、p5与p2、p6对称于转轴中心；四位置转停方案仍然是转动角度为180°或90°间隔进行；

(4)将IMU旋转后光纤陀螺仪生成的数据转换到导航坐标系下，得到惯性器件常值偏差的调制形式，包括：

IMU水平方向上的陀螺常值漂移分别为ε_x和ε_y，载体静止条件下，IMU停顿的A、B、C、D四个位置相对于导航坐标系对称，在一个循环周期内四个固定位置上器件偏差产生的水平误差角在导航坐标系下满足： 

IMU正反转动的8个环节内，导航坐标系水平方向上的误差角度分别表示为：

将一个完整转动周期内的IMU停顿位置和转动过程的东向和北向的水平误差角分别相加，得到导航坐标系下累计的水平误差角度，

导航坐标系上水平方向上光纤陀螺仪常值偏差得到调制，该常值偏差对系统的作用为零；

(5)将IMU坐标系下光纤陀螺仪的输出值 

带入捷联系统中采用四元数法对捷联矩阵 

进行更新：

其中： 

为地球自转角速度在导航坐标系下的分量； 

为导航坐标系相对地球坐标系的运动角速度在导航坐标系下的分量； 

为IMU相对导航坐标系的运动角速度在载体坐标系上的分量；

(6)利用石英加速度计的输出值 

和步骤(5)计算的姿态矩阵 

计算出经过IMU旋转调制后载体的速度和位置。

2.根据权利要求1所述的基于单轴四位置转停方案的捷联系统的误差抑制方法，其特征在于所述的将IMU坐标系下光纤陀螺仪的输出值 带入捷联系统中采用四元数法对捷联矩阵 

进行更新包括：

更新四元数和姿态矩阵：

设载体坐标系相对导航坐标系的转动四元数为：

Q＝q₀+q₁i_s+q₂j_s+q₃k_s

其中：i_s、j_s、k_s分别表示IMU坐标系ox_s轴、oy_s轴、oz_s轴上的单位方向向量；

四元数的即时修正通过解四元数微分方程 

来实现： 

其中： 

分别表示载体相对导航坐标系的运动角速度在IMU坐标系ox_s轴、oy_s轴、oz_s轴上的分量；

姿态矩阵 的更新过程为：

。

3.根据权利要求2所述的一种基于单轴四位置转停方案的捷联系统的误差抑制方法，其特征在于所述的计算出经过IMU旋转调制后载体的速度和位置地方法为：

1)计算导航坐标系下加速度：

2)计算载体的水平速度和位置：

根据t₁时刻的载体东向水平速度V_x(t₁)和北向水平速度V_y(t₁)，求取t₁时刻载体水平速度的变化率为：

在t₂时刻水平速度和载体位置分别为：

3)计算载体速度误差和位置误差： 

其中：V_x0、V_y0分别表示初始时刻载体的东向和北向速度；ΔV_x、ΔV_y分别表示载体东向、北向速度的变化量； 

λ₀分别表示初始时刻载体所处位置的经度和纬度； 

Δλ分别表示载体的纬度、经度的变化量；R_xp、R_yp分别表示地球子午圈、卯酉圈的曲率半径；t₁、t₂为捷联系统的解算过程中两个相邻的时间点。 

Images