- 10 -
高 新 技 术
０ 引言
潜艇在水下执行情报搜集、军事打击以及潜伏等任务时，必须要保证它的隐蔽性和导航定位的准确性，这样才能满足潜艇对水下安全航行及精确打击的需求。惯性导航系统具有自主、高精度、高隐蔽性以及不易受到外界信息干扰等优点，因此潜艇主要依靠惯性导航系统在水下航行。但是惯性导航系统存在系统误差，并且该误差随时间的累积而增大，这就会降低潜艇在水下运行时的精度。针对惯性导航系统存在的不足之处，许多学者提出潜艇在水下航行时，需要修正潜艇惯性导航系统的累积误差；将重力信息与惯性导航系统联合起来的无源辅助导航系统可以对惯性导航的累积误差进行校正，它是目前国际研究的热点。重力匹配导航系统不需要水下潜器浮出或接近水面，测量时也不会向外辐射能量，因此有望实现水下潜器连续长航时精确、自主定位。 １ 重力匹配导航系统原理
重力匹配导航系统指的是潜艇通过自身搭载的重力仪来获
得潜艇目前位置重力场的特征信息，将获得的实测重力值与海洋重力基准图进行匹配，按照合适的匹配算法最终确定潜艇所在的位置。海洋重力基准图是重力匹配导航的重要组成部分，因此，构建全球高精度和高空间分辨率的海洋重力基准图是非常必要的。船载重力测量、航空重力测量以及卫星重力/测高反演是获取全球海洋重力场的主要技术途径，全球高精度和高空间分辨率的海洋重力基准图能有效提高水下重力匹配导航的精度；然而目前全球海洋重力基准图的质量较低（精度为3 mGal~8 mGal，空间分辨率为1′×1′），因此需要对其进行一定程度的插值加密，将其重构为高精度和高空间分辨率的海洋重力基准图。 ２ 海洋重力异常图插值方法
目前常见的重力场插值方法有距离加权反比插值法、径向基函数插值法、样条函数插值法、 克里金插值法以及改进的Shepard 插值法。
２．１ 距离加权反比插值法
距离加权反比插值法是最简单的插值方法，该方法根据未知点到已知点距离的加权平均值来计算未知点的值。距离加权反比插值法的具体计算如公式（1）所示[1]。
（1）
式中：
G 0为预测点的重力值；G i 为第i 个已知点的重力值；h i 为第i 个已知点到预测点的距离；
μ为指数参数。２．２ 径向基函数插值法
径向基函数（Radial Basis Function， RBF）插值法能够应
用于各个领域，例如全球地形、水文、地球物理科学和图像重
建，其原理是将已知点作为位置中心，求解得出插值函数的系数，从而获得预测点的值，径向基函数插值法的插值函数的表
达式如公式（2）所示[2]
彭真。
（2）
式中：
为插值函数；
n 为已知点的总数；为与第i 个点相关的权系数；
刘招华
为基函数。
２．３ 样条函数插值法
样条函数插值法是由工程制图中的样条绘图工具发展而
来的，一般指的是三次样条函数插值法（Cubic spline function interpolation）。三次样条函数插值法是由几段三次多项式拼合而成的曲线，如公式（3）所示[3]。
S i （x ）=a i +b i x +c i x 2+d i x 3 （3）式中：
S i （x ）为三次样条函数；a i 、b i 、c i 和d i 为待求系数。２．４ 克里金插值法
克里金插值法（Kriging）主要采用地质学中的统计理论和
插值方法，以变异函数为插值基础和基本工具，选择合适的变异函数，从而对空间结构和参数进行最优线性无偏估计。克里金插值法分为以下5步：1） 重力场数据的预处理。2） x 和y 2个方向的实验方差图。3） 拟合最优变异函数计算变异函数
（半方差函数）的系数。4） 克里金交叉验证。5） 未知点的克里金估计值。变异函数描述了各个点的空间关系，其表达式如公
式（4）所示[4]
。
（4）
式中：
h 为点x i 到点x i +h 之间的距离；n 为已知点的总数；G （x i ）为点x i 的重力值；γ（h ）为变异函数。
２．５ 改进的二次曲面Ｓｈｅｐａｒｄ插值法
改进的二次曲面Shepard （Modified Quadratic Shepard， MQS）插值法可用于对大规模散乱数据的插值拟合，与距离加权反比插值法相似。改进的二次曲面插值法对领域内外进行主要
电网谐波治理装置重力插值重构前沿与进展
王 彤
（中国空间技术研究院，北京 100094）
摘 要：水下导航一直是导航领域重要的研究内容，重力匹配导航凭借具有自主、高隐蔽性等优点，在导航领域拥有良好的应用前景。该文详细分析并比较了目前常用的海洋重力插值计算方法，例如距离加权反比插值法、径向基函数插值法、样条函数插值法、克里金插值法和改进的Shepard 插值法。通过仿真模拟，比较了5种方法的插值精度，结果表明克里金插值法的精度最高。关键词：水下导航；重力匹配导航；插值方法中图分类号：P 631.1 文献标志码：A
- 11 -
2011年高考数学
高 新 技 术
区分，对进入领域的点进行二次曲面加权运算，领域外的权系数为0，不进行运算。二次曲面函数如公式（5）所示[5]。
Q i （x ）=G i +a i 1（x -x i ）2+a i 2（x -x i ）（y -y i ）+a i 3（y -y i ）2+a i 4
（x -x i ）+a i 5（y -y i ） （5）
式中：
x i 、y i 为第i 个点的经纬度坐标；a i 1、a i 2、a i 3、a i 4和a i 5为二次曲面函数的待求系数；
G i 为第i 个点的重力值；Q i （x ）为二次曲面函数。
可以通过最小二乘逼近法对已知点的重力值进行拟合求出。该文对上述5种方法进行深入分析，根据插值、精度计算复杂度、计算时间以及适用范围等评价指标对常见的插值方法进行统计，各个方法的优缺点见表1。
表1 插值法的优缺点对比分析
插值方法
优点缺点
距离加权反比插值法计算简单、运行时间短精度较差径向基函数插值法精度高、计算时间短适合数据集变
化平缓的情况
数据集变化剧烈时插值效果不佳
样条函数插值法精度高、计算时间较短
重力场高频区域插值时误差较大克里金插值法精度高
计算时间长，边缘数据的插值
效果不好改进的Ｓｈｅｐａｒｄ插值法
计算简单、计算时间短
精度略高于距离加权反比插值法
除了上述5种插值方法，目前很多学者将其他领域的插值方法也应用到重力场插值的过程中，提出了一些改进的插值方法，例如广义抽样扩展插值法、双三次孔斯曲面插值法、不等权最小二乘预测插值和窗口移除恢复技术—克里金插值法等，这些方法的插值精度较高，能有效提高重力场的空间分辨率。
３ 团队研究进展
研究团队从事航天、导航、测绘以及海洋等交叉学科的前
沿探索与军事应用，对卫星重力学、卫星测高学以及水下导航学开展研究，构建了全球高精度和高空间分辨率的海洋重力基准图，进而为实现自主、隐蔽和高精度的全球水下导航提供了理论和方法支持。该文将上述5种传统插值法进行实验验证和对比分析，选取我国南海区域（东经115 °E~116 °E、北纬20 °N~21 °N）的格网化海洋重力异常数据作为实验数据（空间分辨率为1′×1′），如图1所示。
图1（a）和图1（b）分别表示研究区域空间分辨率为1′×1′的二维和三维海洋重力异常图。该区域内海洋异常变化范围较大，因此不同插值方法的误差变化也较大。首先，将空间分辨率为1′×1′的实验数据稀疏为4′×4′；其次，分别用距离加权反比插值法、径向基函数插值法、样条函数插值法、克里金插值法和改进的Shepard 插值法将空间分辨率为4′×4′的海洋重力异常插值为空间分辨率为1′×1′的海洋重力异常；最后，求得与原始1′×1′的海洋重力异常的误差，如图2所示。图2为用这5种方法对我国南海区域插值后的误差散点图，可以看出克里金插值法误差分布较为集中，误差最小。插值重构时，为了与原始图片对比，分析插值后海洋重力异常误差的大小，在判断插值方法精度时，一般选用
平均误差（Mean error， ME）和均方根误差（Root mean square error， RMSE）作为评价指标，计算公式如公式（6）所示[6]。
（6）
RMSE ME 式中：
为x i 预测点的重力值；
为x i 实测的重力值。
平均误差和均方根误差越大，表明精度越低，插值结果越差。因此，平均误差和均方根误差反映了插值精度的大小。
图2的统计结果见表2。根据表2可知，与其他4种插值方法相比，距离加权反比插值法的误差较大，其平均误差为1.07×10-1 mGal，均方根误差为2.96 mGal，对应图2（a）误差的变化范围较大，误差点分布离散。因此，距离加权反比插值法的插值精度较低，插值效果不好。克里金插值法精度最高，平均误差为3.16×10-2 mGal，均方根误差为1.113 mGal，其误差较小，插值效果较佳，适用于对海洋重力场的插值重构。
表2 5种插值方法对我国南海区域插值误差统计
插值方法
平均误差／ｍＧａｌ均方根误差／ｍＧａｌ
距离加权反比插值法１．０７×１０－１２．９６径向基函数插值法３．１７×１０－２２．５８样条
函数插值法３．５２×１０－２１．１６克里金插值法３．１６×１０－２１．１３改进的Ｓｈｅｐａｒｄ插值法
西溪海小海７．７８×１０－２
１．５３
立方氮化硼砂轮
４ 结语
重力匹配导航是1种自主无源的导航方式，对潜艇的水下航行具有极其重要的意义。海洋重力基准图对重力匹配导航精度具有一定的影响。该文对目前常见的重力场插值方法的拟合过程进行研究，通过分析对现有的距离加权反比插值法、径向基函数插值法、样条函数插值法、克里金插值法和改进的Shepard 插值法进行对比，体现出克里金插值法和样条函数数插值法对海洋重力异常插值重构的可行性，也为未来构建全球海洋重力场基准图提供了可行性建议。
参考文献
[1]王博，周明龙.水下重力辅助导航适配区选取的研究进展[J].导航定位学报，2020，8（3）：36-43.
[2]王先泽，李忠科，张晓娟，等.特征保持的基于紧支径向基函数的点云简化[J].计算机工程与设计，2013，34（1）：201-206.[3]谭文，方淼，段峰，等.基于三次样条插值函数的3D 激光
（a） 二维（2D）
（b） 三维（3D）
图1 分辨率为1′×1′的海洋重力异常图
重力异常/m G a l
经度
/°E
经度/°E
纬度
/°N
纬度/°N
- 12 -
高 新 技 术
图像去噪算法研究[J].中南大学学报：自然科学版，2020，51（9）：2496-2503.
[4]Oliver M，Webster R.A tutorial guide to geostatistics ：Computi ng and modelling variograms and kriging[J].Catena， 2014，113（2）：56-69.
[5]吴太旗，黄谟涛，陆秀平，等.重力场匹配导航的重力图生成技术[J].中国惯性技术学报，2007，15（4）：438-441.
[6]Vasat R，Pavlu L，Boruvka L，et al.Mapping the topsoil pH and hu-mus quality of forest soils in the North Bohemian Jizersk hory Mts.Region with ordinary，universal，and regression kriging ： Cross-Validation comparison[J].Soil & Water Research， 2013，8（3）：97-104.
图 2 不同方法插值后误差散点图
（a） 距离加权反比插值法的误差点序列
（b） 径向基函数插值法的误差点序列
（c） 样条插值法的误差点序列（d） 克里金插值法的误差点序列
（e） 改进的Shepard 插值法的误差点序列
重力异常误差/m G a l
重力异常误差/m G a l
重力异常误差/m G a l
重力异常误差/m G a l
重力异常误差/m G a l
1 0001 000 1 0001 0001 000
2 0002 000 2 0002 0002 000
3 0003 000 3 0003 0003 000
4 000
4 000
4 000
4 000
4 000
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议