matlab n维插值,简单调研多维插值方法

原标题：简单调研多维插值方法

■2020-11-07 11:36:50

以前的时候用过二维的插值, 见 二维三次卷积插值算法及Fortran代码 [1], 也用过matlab自带的插值方法, 见 Matlab插值曲面及其曲率的计算 [2]. 最近因为需要自己实现多维插值方法, 就进行了简单的调研, 简单地罗列一下所得.

中文世界

先简单搜索一下中文的资料, 看看中文世界能找到什么信息, 发现大部分都是地学相关的. 虽然参考意义不大, 但至少熟悉一下相关的名词.

■"自然地理要素空间插值的几个问题", 地理研究, 23(4), 2004■"多维线性插值方法的构建及应用分析", 南水北调与水利科技, 11(4), 2013■"复杂地形影响下四川省逐月气温空间插值方法研究", 生态与农村环境学报, 35(6), 2019

常用的空间插值方法可以分为几何插值法、数学插值法以及地学插值方法。

几何插值方法的基本原理是距离越近相似性越高, 常用方法包括泰森多边形法和反距离权重插值法, 适用于下垫面特征较为一致的地区, 而在地形起伏较大的区域其插值误差普遍较大。

数学插值方法是假定插值要素在研究区域呈一定的曲面分布, 通过拟合曲面方程对未知站点进行插值, 常见方法有趋势面法和样条插值法。

地学插值法主要指克里格插值法, 其原理是利用区域化变量数据和半变异函数的结构特点, 对未知点的区域化变量进行最优无偏估计。

数学插值方法和地学插值方法均是根据数学原理通过已知点的信息拟合区域曲面方程或者半变异函数对未知点变量进行无偏估计, 因对区域下垫面差异性较大的区域具有一定的适用性, 且通常已知站点密度越大插值结果精度越高。然而由于区域下垫面的差异性的分布特征, 其空间变量场的分布是复杂多样的, 如何确定合理的曲面方程和半变异函数来准确描述空间变量场的分布特征是十分困难的。实际上, 区域变量场的分布在空间维度的上通常满足一定的递变规律, 因此在空间任一维度上, 均可通过线性插值对未知点