1. 引言
随着一带一路建设的推进,高效而精确地测量海底地形就显得尤为重要。这就需要将有限的等深线数据,在科学的插值方法下,得到尽可能准确的连续水深数据 [1] 。水深测量数据具有数据大、离散化、网格化、范围广、时空域等特点。传统的水深插值方法如反距离加权插值、全局多项式插值法等,运用这些方法可以较好的内插出水深值,得到较好的海底地形模型。但由于时代的限制,反距离加权插值算法虽简便易行,但计算量大,速度慢;全局多项式插值的适用范围相对有限,易出现病态方程组,对精度影响较大,计算复杂 [2] [3] [4] 。而阵列代数可以有效地将网格状大规模的水深数据块分解为小数据块进行高效率的处理和分析,大大提高了效率;并与最小二乘拟合思想相结合,借助正规方程求出最小二乘解,在一定程度上提高了精度,在水深测量内插计算中具有较好的优势 [5] [6] 。
2. 基于阵列代数的双二次多项式
传统双二次多项式可表示为:
(1)
将(1)通过阵列代数的思想变换为矩阵形式:
(2)
线状测深模式是海洋测深的基本作业模式 [2] 。船载测深仪在海面沿计划测线边进行航行测量,获得多条测线上的水深数据,在测区按一定间距加测垂直交于计划测线的检查测线,以便获得交叉点附近的差值数据来检核测深成果。这样,测线和检查线便形成网状,如图1所示。
假设取九条主测线,垂直于主测线选择九条检查测线,如图所示。可以得到
个测点
及对应的水深值
,
。则得到:
(3)
对于形如
形式的矩阵方程。设观测阵列方程为式(4),即求其最小二乘解最为关键。
(4)
3. 基于正规方程的最小二乘解法
本文用
表示所有
阶实矩阵全体。
表示矩阵
与
的Kronecker乘积,即若
,
,则
。
表示矩阵
的按行拉直,即若
,其中
为
的行向量,而
。
设
,考虑线性矩阵方程:
(5)
运用Kronecker乘积,方程(5)等价于
(6)
从而得到方程(6)的正规方程:
(7)
由(7)进一步化简得到:
(8)
将(8)还原可得:
(9)
由于
与
均为范德蒙矩阵,所以可知
与
皆为可逆矩阵,于是由方程(9)得到其最小二乘解为:
(10)
由图2,可以看出,由于
为
的矩阵,
也为
的矩阵,因此,计算速度优于普通的接线性方程组。
4. 实验对比分析
为验证插值效率和所构造的双二次插值最小二乘内插模型的性能,取双曲抛物面函数:
(11)
写成矩阵形式为:
![](//html.hanspub.org/file/10-2620467x49_hanspub.png)
Figure 2. Process of decomposing large matrices into small matrix blocks by array algebra
图2. 阵列代数将大矩阵分解为小矩阵块的过程
(12)
取
并等分为8份,得到结点
;取
并等分为8份,得到结点
。可得到样本点
。采用基于阵列代数的双二次最小二乘拟合模型进行测试,用时0.010秒,拟合结果如下:
(13)
取
并等分为49份,得到结点
;取
并等分为49份,得到结点
。可得到检查点
,求其拟合精度的公式为:
(14)
其中,
为标准差,
为真值,
为拟合值,
为真值
与拟合值
的差值,
为试样总数
,
为各次拟合值。
由式(14)可得精度结果:
(15)
如图3所示,图中曲面为基于阵列代数的双二次插值最小二乘拟合模型及原函数图像的叠加结果图,可以误差很小,肉眼不可见。
用Matlab编程,使用反距离加权插值法(IDW),同样的采样点,当
时得到内插标准差较小,标准差为:
(16)
同时,将内插出的点与原函数点绘制成图,进行比较。可以看出,当权重为4时,内插出的图像有很多小的凹陷,不平整,如图3(b)所示。用Matlab编程得到的反距离加权插值法的插值效果一般,该方法总共需要花费时长为0.209 s,远远超出基于阵列代数的双二次插值最小二乘内插模型的时间。
采用传统不完全双二次函数进行拟合。不完全双二次函数如下式(17)所示。
(17)
(a) (b)
Figure 3. Image comparison between different fitting method (a) Array Algebra; (b) IDW)
图3. (a) 本文方法图像与原函数图像对比图(b) IDW(
)与原函数图像对比图
![](Images/Table_Tmp.jpg)
Table 1. Comparison of time and precision of three methods
表1. 三种方法用时与精度比较
用Matlab编程用时为0.072秒,得到的拟合结果为:
(18)
即:
(19)
标准差为:
(20)
比较三种方法的用时和精度,可以得到表1。
由表1可以看出,相比较传统双二次多项式最小二乘拟合,基于阵列代数的双二次多项式最小二乘拟合算法能够有效地分解大型稀疏矩阵,不需要复杂的奇异值分解等分解形式,即可求出拟合值,同时具有较好的精度和效率。相比于常用水深内插方法反距离加权插值法,在用时和精度方面,本论文方法都有较大优势。
5. 结束语
双二次插值最小二乘内插模型属于整体内插方法,容易受区域内极值的影响,适用于较平缓地区;同时,阵列代数的应用要求所给的数据必须是格网数据。基于阵列代数的双二次插值最小二乘内插模型,将水深测量的大数据块分割成若干小数据块,同时将双二次多项式与最小二乘拟合相结合得出基于正规方程的最小二乘解,一定程度上提高了计算精度,同时有效得提高了计算机的运算效率,减小了计算机存储量,具有非常重要的实用价值。
基金项目
国家自然科学基金项目(41576105, 41631072)。