一种大方位失准角线性对准方法

摘要

本发明涉及一种大方位失准角线性对准方法，步骤如下：(1)根据粗对准系统的状态方程和观测方程，利用GPS观测信息，通过线性卡尔曼滤波，实现粗对准过程，直至航向误差角满足小角度阈值条件；(2)保留粗对准收敛时的系统协方差阵，并将之作为精对准过程的初始条件；(3)延续粗对准的系统状态变量及其位置、速度误差方程和观测方程，使用本发明给出姿态误差方程，并将保存的协方差阵作为初始条件，进行精对准，至系统收敛至预期水平。本发明使系统误差的协方差阵能从粗对准模型直接传递给精对准模型，实现平稳的模型切换过程，改善精对准收敛过程。

权利要求

1.一种大方位失准角线性对准方法，其特征在于步骤如下：(1)选择经度误差，纬度误差，高程误差，东向、北向和天向速度误差，俯仰角误差，滚动角误差，以及航向角误差的正弦项和余弦项作为状态变量，根据状态方程和观测方程，利用GPS观测信息，通过线性卡尔曼滤波，进行粗对准；其中状态方程由位置误差方程、线性化的速度误差方程和姿态误差方程组成；

(2)当粗对准的航向角误差收敛到满足阈值条件时，保存协方差阵；

(3)延续步骤(1)的状态变量及其位置误差方程、速度误差方程和观测方程，更新姿态误差方程，并将保存的协方差阵作为初始条件，进行精对准，至姿态角收敛至预期水平；

所述步骤(2)中的阈值条件是指当航向角误差 收敛至使 成立的阈值

2.根据权利要求1所述的一种大方位失准角线性对准方法，其特征在于：所述步骤(3)中，姿态误差方程为：其中θ、γ和 分别为俯仰、滚动和航向角误差， 为计算坐标系相对惯性坐标系的角速率在计算坐标系中的投影，εx、εy和εz为陀螺仪三轴漂移。

说明书

一种大方位失准角线性对准方法

技术领域

[0001] 本发明涉及一种大方位失准角线性对准方法，属于惯性导航技术领域。

背景技术

[0002] 初始对准是惯性导航的关键技术，也是INS/GNSS组合导航关键技术之一。在基于MEMS-INS/GNSS组合导航系统中，由于MEMS器件特别是陀螺仪的限制，导致无法通过自对准实现方位失准角的初始化，从而造成了大方位失准角问题，其中一种解决方法是直接进行动基座的对准。

[0003] 建立准确的INS误差传播方程和采用适当的滤波技术是进行初始对准的主要问题。在大方位失准角的情况下的动基座对准模型本质上式非线性的，而非线性滤波不适宜工程应用，故多采用线性化方法。

[0004] 现有的技术方案将对准分为粗对准和精对准两个过程，在粗对准精度达到某一阈值条件时实现切换。在用欧拉角法描述失准角的条件下，粗对准将游动方位坐标系中的姿态误差状态替换为正弦和余弦项，从而实现方程线性化。同时，针对游动方位的三角函数可观度较弱的问题，还提出一种显著增加了可观度的改进状态方程。在精对准过程中，由于经过粗对准的失准角已足够小，可以通过sinα≈α，cosα≈1实现线性化。

[0005] 现有对线性滤波方案中，在粗对准和精对准中使用不同的状态变量和状态方程，这造成了由粗到精过渡时的模型切换问题。由于系统状态变量不一致，粗对准得到的协方差阵不能直接用于精对准，相应的，精对准过程也需要重新对设定协方差阵。由于协方差阵的损失，模型切换往往使滤波过渡无法平稳实现，影响精对准过程的收敛速度。

发明内容

[0006] 本发明技术解决问题：克服现有技术在粗对准和精对准模型切换中的过渡不平稳问题，提供了一种新型大方位失准角线性对准方法，使系统误差的协方差阵能从粗对准模型直接传递给精对准模型，实现平稳的模型切换过程，改善精对准收敛过程。

[0007] 本发明的技术要点：

[0008] 1.粗对准过程延续现有方案；

[0009] 2.精对准过程中，延续粗对准系统的状态变量，即系统状态变量定义为：

[0010]

[0011] 其中，L、λ和h分别为纬度、经度和高度，δVE、δVN和δVU分别为东向、北向和天向的速度误差，θ、γ和 分别为俯仰、滚动和航向角误差。

[0012] 粗对准系统的姿态误差方程为：

[0013]

[0014] 其中 为计算坐标系相对惯性坐标系的角速率在计算坐标系中的投影，εx、εy和εz为陀螺仪漂移项，Fs和Fc为非线性项，定义为：

[0015]

[0016] 粗对准中，Fs和Fc近似为0，精对准则不能这样简化。精对准时，式(2)中θ、γ、εz均是小量， 也是随着方位角的减小而趋于零的小量，当方位角减小到8-10度时，的值减小到0.01以内。因此 与θ、γ、εz的乘积为高阶小量，消去高阶小量近似得到：

[0017]

[0018] 此时，得到精对准的系统姿态误差方程：

[0019]

[0020] 精对准中，除姿态误差方程外，其他误差方程保持不变。对比粗对准的姿态误差方程，本发明提出的方案只需对传递函数部分进行调整，并增加陀螺漂移εz即可。

[0021] 本发明技术解决方案：一种大方位失准角线性对准方法，步骤如下：

[0022] (1)粗对准采用现有的技术方案，选择经度误差，纬度误差，高程误差，东向、北向和天向速度误差，俯仰角误差，滚动角误差，以及航向角误差的正弦项和余弦项作为系统状态变量，根据粗对准系统的状态方程和观测方程，利用GPS观测信息，通过线性卡尔曼滤波，进行粗对准。其中，状态方程由位置误差方程、线性化的速度和姿态误差方程组成；

[0023] (2)当粗对准的航向角误差收敛到满足阈值条件 时，保存滤波系统的协方差阵，其中 为设定的阈值，使 成立。

[0024] (3)保持粗对准的系统状态变量及其位置、速度误差方程和观测方程不变，使用本发明给出姿态误差方程，并将保存的协方差阵作为初始条件，进行精对准，至系统收敛至预期水平。本发明给出的系统的姿态误差方程为：

[0025]

[0026] 本发明与现有技术相比的优点在于：现有的粗对准至精对准的切换方案，粗对准和精对准使用不同的系统变量，相应的需要不同的系统模型。而本发明给出的切换方案的系统状态量保持一致，就状态方程而言，只需要将姿态误差方程的传递函数作微调，并添加陀螺漂移分量εz，其他方程保持变即可。

[0027] 系统状态量保持不变的优势在于，衡量对准水平的系统协方差阵可以直接从粗对准过渡到精对准，使精对准过程有准确的初始条件。故而本发明给出的切换方案能够保持对准模型的一致性，实现平稳过渡。

附图说明

[0028] 图1为本发明方法实现流程图。

具体实施方式

[0029] 如图1所示，本发明具体实现如下：

[0030] (1)粗对准过程

[0031] 定义系统状态变量为：

[0032]

[0033] 其中，L、λ和h分别为纬度、经度和高度，δVE、δVN和δVU分别为东向、北向和天向的速度误差，θ、γ和 分别为俯仰、滚动和航向角误差。另外，‘x’表示任意变量x的计算值，‘δx’表示任意变量x的误差。

[0034] 速度误差方程：

[0035]

[0036] 其中， 和 分别为地球自转角速率和计算坐标系相对地球坐标系的角速率在计算坐标系中的投影，fc为真实比力在计算坐标系上的投影， 分别为三个方向上加速度计的零偏。式中F为非线性项，将模型线性化认为此项在车载应用中是小量，在粗对准时可以省略。

[0037]

[0038] 姿态误差方程：

[0039]

[0040] 其中， 分别计算坐标系的角速率在计算坐标系中的投影，εx和εy为陀螺仪漂移项。

[0041] 位置误差方程：

[0042]

[0043] 其中RM,RN分别是子午圈半径和卯酉圈半径。

[0044] 姿态误差方程、速度误差方程和位置误差方程共同构成状态方程。

[0045] 动基座对准中采用GPS输出的速度和位置作为观测量，观测方程如下：

[0046]

[0047] 式中

[0048]

[0049] pGPS和vGPS分别为GPS提供的位置和速度；δpGPS和δvGPS为GPS的位置和速度误差，构成观测噪声n；pIMU和vIMU分别为INS提供的位置和速度；0m×n和Ik分别表示大小为m×n的零矩阵和k×k的单位矩阵。

[0050] 最后，通过状态方程和观测方程构成线性卡尔曼滤波器，实现粗对准过程。

[0051] 当航向误差角收敛至满足小角度条件时，结束粗对准过程。

[0052] (2)协方差阵传递

[0053] 保留粗对准收敛时的系统协方差阵，将其作为精对准过程的协方差阵初始条件。

[0054] (3)精对准过程

[0055] 当粗对准的航向误差收敛到满足阈值条件 时，切换到精对准过程。精对准过程的系统状态变量、状态方程和观测方程延续粗对准过程，只将姿态误差方程替换为如下方程即可。

[0056]