[发明专利]一种基于PC-104的高精度惯导系统扰动重力补偿方法

说明书

本发明涉及一种基于PC‑104的高精度惯导系统扰动重力补偿方法，通过使用惯导系统将位置信息输入至PC‑104，PC‑104将输入位置信息中的大地经纬度转换为地心经纬度后，使用地心经纬度进行Legendre函数及其导数的计算，把Legendre函数计算结果和完全规格化位系数导入全球重力模型EGM2008计算垂线偏差东西向和南北向分量，通过垂线偏差东西向和南北向分量计算惯导系统重力补偿。本方法能够精确计算垂线偏差东西向和南北向分量，并通过计算得到的垂线偏差东西向和南北向分量对惯导系统进行重力补偿，满足了惯导系统的垂线偏差低频补偿需求，提高了高精度惯导系统的精确性，为高精度惯导系统的发展奠定了基础。

技术领域

本发明属于重力计算领域，尤其是一种基于PC-104的高精度惯导系统扰动重力补偿方法。

背景技术

CHAMP、GRACE和GOCE等卫星重力计划的成功实施使得全球重力场模型向着超高阶次发展，精度也得到较大提高提升。目前常用的全球重力场模型最高阶次已达到2190阶，综合重力观测数据和地形数据获得的全球网格重力数据的最高分辨率已达到7.2″×7.2″，这使我们对地球重力场的认识不断深化，通过重力场模型算得的重力扰动也越来越精确，给我们提供了补偿重力扰动的理论基础。

其中全球重力场的精度、稳定性和效率主要取决于缔合勒让德函数的精度、稳定性和效率，现阶段缔合勒让德函数的计算方法主要有两类：直接计算法和递推计算法。对于超高阶重力场模型，直接计算法在计算精度和效率方面均不能满足要求，所以计算方式主要采用四种方法：标准向前行递推法、标准向前列递推法、跨阶次递推法和Belikov递推法。研究证明标准向前列递推法和标准向前行递推法适用范围在1900阶以内，并且适用效果随着阶数的增加而减小，超过1900阶则不适用。跨阶次递推法和Belikov递推法在3000阶以内计算均适用。

但是在惯性器件精度不高的系统中，与原件误差相比，扰动重力对系统的测数影响可以忽略，随着惯性器件精度的不断提高与高精度导航系统需求的发展，扰动重力成为高精度惯导系统的一项主要误差。重力扰动分为垂线偏差和重力异常连部分，其中垂线偏差的影响远大于单纯的重力异常的影响。由于垂线偏差在系统中的影响可以等同为加速度计零偏，那么如果知道测量点的垂线偏差，补偿将会方便计算。

发明内容

本发明的目的在于克服现有技术的不足，提出一种基于PC-104的高精度惯导系统扰动重力补偿方法，精确计算惯导系统的重力补偿。

本发明解决其技术问题是采取以下技术方案实现的：

一种基于PC-104的高精度惯导系统扰动重力补偿方法，包括以下步骤：

步骤1、用户下载完全规格化位系数的文本文档并与惯导系统位置信息输入至PC-104；

步骤2、PC-104将惯导系统位置信息中的大地经纬度转换为地心经纬度；

步骤3、PC-104使用地心经纬度进行Legendre函数及其导数的计算；

步骤4、PC-104将计算得到的Legendre函数及其导数与完全规格化位系数文本文档中完全规格化位系数导入全球重力模型EGM2008进行垂线偏差计算，得到垂线偏差东西向和南北向分量；

步骤5、PC-104将计算得到垂线偏差东西向和南北向分量输入至惯导系统，惯导系统根据垂线偏差东西向和南北向分量对扰动重力进行补偿。

而且，所述步骤2中的具体方法为：

Lat_Earthcore＝arctan(Lat_Geo)/((1-f)^2)

其中，Lat_Geo为大地经纬度，Lat_Earthcore为地心经纬度，f为地球扁率。