阅读说明：本技术 一种基于协方差分析的弹类惯导导航误差分析方法及装置 (A kind of bullet class inertial navigation navigation error analysis method and device based on the analysis of covariance ) 是由 夏家和 李伟 张金亮 吉翠萍 于 2019-08-16 设计创作，主要内容包括：本发明提供一种基于协方差分析的弹类惯导导航误差分析方法及装置,包括：获取关于惯导建模误差状态的高阶惯导误差传播模型；对高阶惯导误差传播模型进行离散化处理,得到上一时刻和当前时刻的惯导建模误差状态的转移关系；根据状态转移关系,弹道递推计算导航参数误差统计值。本发明解决仿真时间较长和仿真样板需求大的问题。(The present invention provides a kind of bullet class inertial navigation navigation error analysis method and device based on the analysis of covariance, comprising: obtains the high-order ins error propagation model about inertial navigation modeling error state；Sliding-model control is carried out to high-order ins error propagation model, obtains the transfer relationship of the inertial navigation modeling error state at last moment and current time；According to state transfer relationship, trajectory recurrence calculation navigational parameter error statistics value.The present invention solves the problems, such as that simulation time is big compared with long and emulation template demand.)

一种基于协方差分析的弹类惯导导航误差分析方法及装置

技术领域

本发明属于惯性导航技术领域，涉及一种基于协方差分析的弹类惯导导航误差分析方法及装置。

背景技术

在惯导方案设计阶段，为实现最终系统精度需对各种误差进行预估及分配计算。导航误差的传播分析主要可通过两种方法来完成，一个是蒙特卡洛(Monte Carlo)仿真方法，这种方法可以应用于任意系统类型(适用于非线性系统)，但是需要较长的时间来完成计算，特别是惯导系统误差源很多，为获得较准确统计特性，需很大量仿真样本。

发明内容

本发明解决的技术问题：提供一种基于协方差分析的弹类惯导导航误差仿真分析方法，解决仿真时间较长和仿真样板需求大的问题。

本发明的技术方案：

本发明提供一种基于协方差分析的弹类惯导导航误差分析方法，包括：

获取关于惯导建模误差状态的高阶惯导误差传播模型；

对高阶惯导误差传播模型进行离散化处理，得到上一时刻和当前时刻的惯导建模误差状态的转移关系；

根据状态转移关系，弹道递推计算导航参数误差统计值。

进一步的，所述方法还包括：

建立高阶惯导误差传播模型；高阶惯导误差传播模型是关于惯导建模误差状态的误差传播微分模型；

给定高阶惯导误差传播模型中各误差源的初始统计特性值。

进一步的，对高阶惯导误差传播模型进行离散化处理，得到上一时刻和当前时刻的惯导建模误差状态的转移关系包括：

对高阶惯导误差传播模型进行离散化处理，得到离散后惯导建模误差状态的转移矩阵的递推公式，递推公式为：

这里，Φ_k|k-1为离散后惯导建模误差状态的状态转移矩阵，F为连续系统转移阵，Δt为离散化时间步长；Q_k-1为离散化后的系统噪声阵，Q为连续系统噪声阵。

进一步的，给定高阶惯导误差传播模型中各误差源的初始统计特性值之前，所述方法还包括：

建立初始的误差协方差阵，误差协方差阵包括数学平台偏角的协方差阵、速度误差的协方差阵、位置误差的协方差阵。

进一步的，其特征在于，误差协方差阵还包括陀螺漂移的协方差阵、加计零位的协方差阵、陀螺刻度系数安装偏角的协方差阵、加计刻度系数安装偏角相应的协方差阵。

进一步的，根据状态转移关系，弹道递推计算导航参数误差统计值包括：

根据状态转移关系，弹道递推计算误差协方差阵；

从误差协方差阵提取导航参数误差统计值。

进一步的，根据状态转移关系，弹道递推计算误差协方差阵包括：

通过弹道递推公式计算误差协方差阵，所述弹道递推公式为：

这里，P_K、P_K-1分别为k、k-1时刻的协方差阵。

本发明提供一种基于协方差分析的弹类惯导导航误差分析装置，包括：

获取模块，用于获取关于惯导建模误差状态的高阶惯导误差传播模型；

离散化处理模块，用于对高阶惯导误差传播模型进行离散化处理，得到上一时刻和当前时刻的惯导建模误差状态的转移关系；

计算模块，用于根据状态转移关系，弹道递推计算导航参数误差统计值。

进一步的，所述装置还包括：

建立模块，用于建立高阶惯导误差传播模型；高阶惯导误差传播模型是关于惯导建模误差状态的误差传播微分模型；

给定模块，用于给定高阶惯导误差传播模型中各误差源的初始统计特性值。

本发明提供一种计算机可读的存储介质，存储有计算机能够运行的程序，计算机程序被处理器执行以实现上述方法。

本发明的有益效果：本发明的弹类惯导导航误差仿真分析方法基于高阶捷联惯导误差模型和预设弹道剖面，通过协方差分析方法获得各导航误差源对惯性导航误差的贡献。方案根据各项惯性器件的统计特性及初始对准精度，基于高阶捷联惯导误差传播模型，通过协方差分析方法获得特性弹道剖面下的纯惯性误差统计特性。本发明可应用于各种惯性导航系统设计阶段的误差预估，给出了一种弹类惯导导航误差仿真分析方法。方法提供了一种系统统计性能的仿真手段，可避免通过大量的单次仿真进行统计而获得系统性能的预估。通过协方差仿真可以评估在特定任务场景下各项误差的贡献，据此可分析哪些误差在系统中是主要误差源，并且在后续系统实现中通过技术手段予以降低，以获得满足需求的导航系统性能。

附图说明

图1为本发明的一种基于协方差分析的弹类惯导导航误差仿真分析方法的流程图。

具体实施方式

根据附图和具体实施方式对本发明作进一步详细说明。

本发明方法是协方差分析方法，基于这种方法如果限定误差的分布模型可以快速得到误差传播的分布规律。对于线性系统，协方差分析方法不但省时间，而且是精确的，同时结果功能直观。

协方差分析是上个世纪八十年代由美国科学分析公司提出的一种系统分析方法，可以直接确定具有随机输入的线性系统统计特性，并用于大型导弹系统的精度分析中，对其非线性采用近似分析法。这个方法的主要优点是大大节省计算时间。本专利设计一种基于协方差分析的弹类惯导导航误差仿真分析方法，该方法只需要运行一次即可获得与蒙特卡罗仿真相似的结果从而大大节省仿真时间。

误差协方差仿真分析是评估降阶卡尔曼滤波算法在真实环境下工作性能的有效方法[9]。协方差分析与仿真方法在卡尔曼滤波器初始设计阶段是非常有效的。协方差仿真提供了一种系统统计性能的仿真手段，通过协方差仿真可以避免通过大量的单次仿真进行统计而获得系统性能的预估。通过协方差仿真可以评估在特定任务场景下各项误差的贡献，据此可分析哪些误差在系统中是主要误差源，并且在后续系统实现中通过技术手段予以降低，以获得满足需求的导航系统性能。

本发明提供一种基于协方差分析的弹类惯导导航误差仿真分析方法，如图1所示，包括如下步骤：

步骤S1、建立高阶惯导误差传播模型；

步骤S2、给定各误差源的初始统计特性值；

步骤S3、对高阶连续模型进行离散化处理；

步骤S4、根据弹道递推计算误差协方差阵，提取导航参数误差统计值。

进一步的，在步骤1中，采用的高阶误差模型为：

这里，

x惯导建模误差状态；F为惯导误差系统矩阵；w为驱动噪声。

为数学平台偏角；

为地速东北天分量；

δp＝[δL δλ δh]^T，δL,δλ,δh分别为纬度、经度、高度误差；

为机体系陀螺漂移；

为机体系加计零位；

δKg＝[δk_gxx δk_gxy δk_gxz δk_gyx δk_gyy δk_gyz δk_gzx δk_gzy δk_gzz]^T，陀螺刻度系数误差及安装偏角；

δKa＝[δk_axx δk_axy δk_axz δk_ayy δk_ayz δk_azz]^T，加计刻度系数误差及安装偏角。

进一步的，在步骤2中，给定各误差源的初始统计特性值，具体计算公式如下：

P_φ＝diag(var(φ_E) var(φ_N) var(φ_U))

P_δv＝diag(var(δv_E) var(δv_N) var(δv_U))

P_δp＝diag(var(δL) var(δλ) var(δh))

这里，diag表示对角矩阵，var表示方差；P(0)为计算起始时刻的误差协方差阵；P_φ、P_δv、P_δp、P_ε、P_δKa、P_δKg分别为数学平台偏角、速度误差、位置误差、陀螺漂移、加计零位、陀螺刻度系数安装偏角、加计刻度系数安装偏角相应的协方差阵；

进一步的，步骤3中，进行离散化处理，具体计算公式如下：

这里，Φ_k|k-1为离散后的状态转移矩阵，F为连续系统转移阵，Δt为离散化时间步长。Q_k-1为离散化后的系统噪声阵，Q为连续系统噪声阵，即为w的统计特性。

进一步的，步骤4中根据弹道递推计算误差协方差，具体计算公式如下：

这里，P_K、P_K-1分别为k、k-1时刻的协方差阵。从P_K可以提取当前时刻导航参数的协方差值。

本发明采取的方案为基于高阶捷联惯导误差模型和预设弹道剖面，通过协方差分析方法获得各导航误差源对惯性导航误差的贡献。方案根据各项惯性器件的统计特性及初始对准精度，基于高阶捷联惯导误差传播模型，通过协方差分析方法获得特性弹道剖面下的纯惯性误差统计特性。本发明可应用于各种惯性导航系统设计阶段的误差预估，给出了一种弹类惯导导航误差仿真分析方法。方法提供了一种系统统计性能的仿真手段，可避免通过大量的单次仿真进行统计而获得系统性能的预估。