插值方法比较Word版
数值分析论文 ――几种插值方法的比较课程论文8(学院+专业+学号)
数值分析论文——几种插值方法的比较1．插值法概述插值法是函数逼近的重要方法之一，有着广泛的应用 。
在生产和实验中，函数或者其表达式不便于计算复杂或者无表达式而只有函数在给定点的函()x f 数值(或其导数值) ，此时我们希望建立一个简单的而便于计算的函数，使()x ϕ其近似的代替，有很多种插值法,其中以拉格朗日(Lagrange)插值和牛顿()x f (Newton)插值为代表的多项式插值最有特点,常用的插值还有Hermite 插值,分段插值和样条插值.这里主要介绍拉格朗日(Lagrange)插值和牛顿(Newton)插值和埃尔米特插值（Hermite 插值）。
2．插值方法的比较2．1拉格朗日插值2.1.1基本原理构造次多项式，这是n ()()()()()x l y x l y x l y x l y x P n n k nk k n +⋅⋅⋅++==∑=11000不超过次的多项式，其中基函数：n()x l k =)...()()...()(()...()()...()(()1110)1110n k k k k k k k n k k x x x x x x x x x x x x x x x x x x x x ----------+-+-显然满足 =()x l k ()i k x l ⎩⎨⎧≠=)(0)(1k i k i 此时,误差()()x f x P n ≈()()()=-=x P x f x R n n (x))!1()(1)1(+++n n n f ωξ其中∈且依赖于，.ξ()b a ,x ()()()()n n x x x x x x x -⋅⋅⋅--=+101ω很显然，当，插值节点只有两个,时1=n k x 1+k x ()()()x l y x l y x P k k k k i 11+++=其中基函数 = ， =()x l k 11++--k k k x x x x ()x l k 1+kk kx x x x --+12.1.2优缺点可对插值函数选择多种不同的函数类型，由于代数多项式具有简单和一些良好的特性，故常选用代数多项式作为插值函数。
空间插值方法对比整理版
• 由于建立在统计学的基础上，因此不仅可 以产生预测曲面，而且可以产生误差和不 确定性曲面，用来评估预测结果的好坏
• 多种 kriging 方法
a
9
3、精确插值和近似插值
• 精确插值：产生通过所有观测点的曲面。
• 在精确插值中，插值点落在观测点上，内插值等 于估计值。
• 近似插值：插值产生的曲面不通过所有观测 点。
a
11
插值方法选择的原则
① 精确性：
② 参数的敏感性：许多的插值方法都涉及到一个或多个参数， 如距离反比法中距离的阶数等。有些方法对参数的选择相当 敏感，而有些方法对变量值敏感。后者对不同的数据集会有 截然不同的插值结果。希望找到对参数的波动相对稳定，其 值不过多地依赖变量值的插值方法。
③ 耗时：一般情况下，计算时间不是很重要，除非特别费时。
空间插值 Spatial Interpolation
• 空间插值的概念 • 空间插值的类型 • 空间插值的方法
a
1
空间插值概念
空间插值——空间插值常用于将离散点的测量数据转换为连 续的数据曲面，以便与其它空间现象的分布模式进行比较， 它包括了空间内插和外推两种算法。空间内插算法：通过已 知点的数据推求同一区域未知点数据。空间外推算法：通过 已知区域的数据，推求其它区域数据。
• 典型例子是：全局趋势面分析 、Fourier Series （周期序列）
a
4
局部内插法
➢ 局部内插法只使用邻近的数据点来估计未知点的值，步骤如 下： • 定义一个邻域或搜索范围； • 搜索落在此邻域范围的数据点； • 选择能表达这有限个点空间变化的数学函数； • 为未知的数据点赋值。
➢ 局部内插方法： • 样条函数插值法 • 距离倒数插值 • Kriging插值（空间自由协方差最佳内插）
[转载]插值算法（一）：各种插值方法比较
![[转载]插值算法（一）：各种插值方法比较](https://img.taocdn.com/s3/m/807a27741fd9ad51f01dc281e53a580216fc500f.png)
[转载]插值算法（⼀）：各种插值⽅法⽐较原⽂地址：插值算法（⼀）：各种插值⽅法⽐较作者：稻草⼈确定性随机性确定性随机性趋势⾯（⾮精确）回归（⾮精确）泰森（精确）克⾥⾦（精确）密度估算（⾮精确）反距离权重（精确）薄板样条（精确）整体拟合利⽤现有的所有已知点来估算未知点的值。
局部插值使⽤已知点的样本来估算位置点的值。
确定性插值⽅法不提供预测值的误差检验。
随机性插值⽅法则⽤估计变异提供预测误差的评价。
对于某个数据已知的点，精确插值法在该点位置的估算值与该点已知值相同。
也就是，精确插值所⽣成的⾯通过所有控制点，⽽⾮精确插值或叫做近似插值，估算的点值与该点已知值不同。
1、反距离加权法（Inverse Distance Weighted）反距离加权法是⼀种常⽤⽽简单的空间插值⽅法，IDW是基于“地理第⼀定律”的基本假设:即两个物体相似性随他们见的距离增⼤⽽减少。
它以插值点与样本点间的距离为权重进⾏加权平均，离插值点越近的样本赋予的权重越⼤，此种⽅法简单易⾏，直观并且效率⾼，在已知点分布均匀的情况下插值效果好，插值结果在⽤于插值数据的最⼤值和最⼩值之间，但缺点是易受极值的影响。
2、样条插值法（Spline）样条插值是使⽤⼀种数学函数，对⼀些限定的点值，通过控制估计⽅差，利⽤⼀些特征节点，⽤多项式拟合的⽅法来产⽣平滑的插值曲线。
这种⽅法适⽤于逐渐变化的曲⾯，如温度、⾼程、地下⽔位⾼度或污染浓度等。
该⽅法优点是易操作，计算量不⼤，缺点是难以对误差进⾏估计，采样点稀少时效果不好。
样条插值法⼜分为张⼒样条插值法（Spline with Tension）规则样条插值法（Regularized Spline）薄板样条插值法 (Thin-Plate Splin)3、克⾥⾦法（Kriging）克⾥⾦⽅法最早是由法国地理学家Matheron和南⾮矿⼭⼯程师Krige提出的，⽤于矿⼭勘探。
这种⽅法认为在空间连续变化的属性是⾮常不规则的，⽤简单的平滑函数进⾏模拟将出现误差，⽤随机表⾯函数给予描述会⽐较恰当。
(完整版)几种插值法比较与应用
多种插值法比较与应用（一）Lagrange 插值 1． Lagrange 插值基函数 n+1个n 次多项式∏≠=--=nkj j j kjk x xx x x l 0)( n k ,,1,0 =称为Lagrange 插值基函数 2． Lagrange 插值多项式设给定n+1个互异点))(,(k k x f x ，n k ,,1,0 =，j i x x ≠，j i ≠，满足插值条件)()(k k n x f x L =，n k ,,1,0 =的n 次多项式∏∏∏=≠==--==nk nkj j jk j k k nk k n x x x x x f x l x f x L 000))(()()()(为Lagrange 插值多项式，称∏=+-+=-=nj j x n n x x n f x L x f x E 0)1()()!1()()()()(ξ 为插值余项，其中),()(b a x x ∈=ξξ （二）Newton 插值 1．差商的定义 )(x f 关于i x 的零阶差商)(][i i x f x f = )(x f 关于i x ，j x 的一阶差商ij i j j i x x x f x f x x f --=][][],[依次类推，)(x f 关于i x ，1+i x ，……，k i x +的k 阶差商ik i k i i k i i k i i i x x x x f x x f x x x f --=+-+++++],,[],,[],,,[1112． Newton 插值多项式设给定的n+1个互异点))(,(k k x f x ，n k ,,1,0 =，j i x x ≠，j i ≠， 称满足条件)()(k k n x f x N =，n k ,,1,0 =的n 次多项式)()](,,,[)](,[][)(10100100---++-+=n n n x x x x x x x f x x x x f x f x N为Newton 插值多项式，称],[,)(],,,[)()()(010b a x x x x x x f x N x f x E nj j n n ∈-=-=∏=为插值余项。
插值方法比较范文
插值方法比较范文插值方法是数值计算中常用的一种数值逼近技术，用于通过已知数据点之间的关系来估计未知数据点的值。
在插值过程中，根据不同的插值方法，可以得到不同的近似函数，从而得到不同的结果。
常见的插值方法包括拉格朗日插值、牛顿插值、埃尔米特插值和样条插值等。
下面将对这些插值方法进行比较，包括优缺点。
首先是拉格朗日插值法，它是通过使用已知数据点的函数值来构建一个多项式，再利用这个多项式来估算未知数据点的函数值。
拉格朗日插值法的优点是简单易懂、计算简便，而且在已知数据点分布较为均匀的情况下效果较好。
然而，拉格朗日插值法的缺点是对于较多数据点的情况，构建的多项式会非常复杂，容易导致插值结果的振荡。
此外，拉格朗日插值法对于增加或减少一个数据点都需要重新计算，不够灵活。
其次是牛顿插值法，它也是通过已知数据点的函数值来构建一个多项式，但是与拉格朗日插值法不同，牛顿插值法利用差商的概念来简化多项式的计算。
牛顿插值法的优点是可以递推计算差商，避免了重复计算，因此对于增加或减少一个数据点时比较方便。
此外，牛顿插值法的插值多项式在已知数据点分布较为稀疏的情况下效果较好。
缺点是对于较多数据点的情况，插值多项式同样会变得复杂，容易导致插值结果的振荡。
再者是埃尔米特插值法，它是拉格朗日插值法的一种改进方法。
埃尔米特插值法不仅利用已知数据点的函数值，还利用已知数据点的导数值来构建插值函数，从而提高了插值的精度。
埃尔米特插值法的优点是可以通过已知数据点的导数值来更好地拟合函数的特点，从而得到更准确的插值结果。
缺点是在计算过程中需要求解一系列线性方程组，计算量较大。
最后是样条插值法，它是常用的插值方法之一、样条插值法通过将插值区间划分为若干小区间，在每个小区间上构建一个低次多项式，通过满足一定的光滑性条件来保证插值函数的平滑性。
样条插值法的优点是插值函数的平滑性较好，能够解决拉格朗日插值法和牛顿插值法的振荡问题。
缺点是在计算过程中需要求解大规模的线性方程组，计算量较大。
（完整word版）几种插值法的应用和比较
（完整word版）⼏种插值法的应⽤和⽐较插值法的应⽤与⽐较信科1302 万贤浩 132710381格朗⽇插值法在数值分析中，拉格朗⽇插值法是以法国⼗⼋世纪数学家约瑟夫·路易斯·拉格朗⽇命名的⼀种多项式插值⽅法.许多实际问题中都⽤函数来表⽰某种内在联系或规律，⽽不少函数都只能通过实验和观测来了解.如对实践中的某个物理量进⾏观测，在若⼲个不同的地⽅得到相应的观测值，拉格朗⽇插值法可以找到⼀个多项式，其恰好在各个观测的点取到观测到的值.这样的多项式称为拉格朗⽇（插值）多项式.数学上来说，拉格朗⽇插值法可以给出⼀个恰好穿过⼆维平⾯上若⼲个已知点的多项式函数.拉格朗⽇插值法最早被英国数学家爱德华·华林于1779年发现，不久后由莱昂哈德·欧拉再次发现.1795年，拉格朗⽇在其著作《师范学校数学基础教程》中发表了这个插值⽅法，从此他的名字就和这个⽅法联系在⼀起.1.1拉格朗⽇插值多项式图1已知平⾯上四个点：(?9, 5), (?4, 2), (?1, ?2), (7, 9)，拉格朗⽇多项式：)(x L （⿊⾊）穿过所有点.⽽每个基本多项式：)(00x l y ,)(11x l y , )(22x l y 以及)(x l y ??各穿过对应的⼀点，并在其它的三个点的x 值上取零.对于给定的若1+n 个点),(00y x ,),(11y x ,………),(n n y x ，对应于它们的次数不超过n 的拉格朗⽇多项式L 只有⼀个.如果计⼊次数更⾼的多项式，则有⽆穷个，因为所有与L 相差))((10x x x x --λ……)(n x x -的多项式都满⾜条件.对某个多项式函数，已知有给定的1+k 个取值点：),(00y x ,……,),(k k y x ,其中i x 对应着⾃变量的位置，⽽i y 对应着函数在这个位置的取值.假设任意两个不同的i x 都互不相同，那么应⽤拉格朗⽇插值公式所得到的拉格朗⽇插值多项式为：)()(0x l y x L j kj j ∑==，其中每个)(x l j 为拉格朗⽇基本多项式（或称插值基函数），其表达式为：)()()()()()()()()(111100,0k j k j j j j j j j kj i i ij i j x x x x x x x x x x x x x x x x x x x x x l --------=--=++--≠=∏ΛΛ，拉格朗⽇基本多项式()x l i 的特点是在j x 上取值为1，在其它的点i x ，j i ≠ 上取值为0. 例：设有某个多项式函数f ，已知它在三个点上的取值为：10)4(=f ， ? 25.5)5(=f ， ?1)6(=f ，要求)18(f 的值.⾸先写出每个拉格朗⽇基本多项式：())64)(54()6)(5(0----=x x x l ；())65)(45()6)(4(1----=x x x l ；())56)(46()5)(4(2----=x x x l ；然后应⽤拉格朗⽇插值法，就可以得到p 的表达式（p 为函数f 的插值函数）：)()6()()5()()4()(210x l f x l f x l f x p ++=)56)(46()5)(4(1)65)(45()6)(4(25.5)64)(54()6)(5(10----?+----?+----?=x x x x x x)13628(412+-=x x ，此时数值18就可以求出所需之值：11)18()18(-==p f .1.2插值多项式的存在性与唯⼀性存在性对于给定的1+k 个点：),(),,(00k k y x y x K 拉格朗⽇插值法的思路是找到⼀个在⼀点j x 取值为1，⽽在其他点取值都是0的多项式)(x l j .这样，多项式)(x l y j j 在点j x 取值为j y ，⽽在其他点取值都是0.⽽多项式()∑==kj jj x ly x L 0)(就可以满⾜∑==++++==ki j j j i y y x l y x L 0000)()(ΛΛ，在其它点取值为0的多项式容易找到，例如：)())(()(110k j j x x x x x x x x ----+-ΛΛ，它在点j x 取值为：)()()(10k j j j i x x x x x x ---+ΛΛ.由于已经假定i x 两两互不相同，因此上⾯的取值不等于0.于是，将多项式除以这个取值，就得到⼀个满⾜“在j x 取值为1，⽽在其他点取值都是0的多项式”：)()()()()()()()(111100k j k j j j j j j j i j j x x x x x x x x x x x x x x x x x x xx l --------=--=++--∏ΛΛ，这就是拉格朗⽇基本多项式. 唯⼀性次数不超过k 的拉格朗⽇多项式⾄多只有⼀个，因为对任意两个次数不超过k 的拉格朗⽇多项式：1p 和2p ，它们的差21p p -在所有1+k 个点上取值都是0,因此必然是多项式)())((10k x x x x x x ---Λ的倍数.因此，如果这个差21p p -不等于0，次数就⼀定不⼩于1+k .但是21p p -是两个次数不超过k 的多项式之差，它的次数也不超过k ，所以021=-p p 也就是说21p p =.这样就证明了唯⼀性.1.3性质拉格朗⽇插值法中⽤到的拉格朗⽇基本多项式n l l l ,,,10Λ（由某⼀组n x x x <<<Λ10 确定）可以看做是由次数不超过n 的多项式所组成的线性空间：[]X n K 的⼀组基底.⾸先，如果存在⼀组系数：n λλλ,,,10Λ使得，01100=+++=n n l l l P λλλΛ，那么，⼀⽅⾯多项式p 是满⾜n n x P x P x P λλλ===)(,,)(,)(1100Λ的拉格朗⽇插值多项式，另⼀⽅⾯p 是零多项式，所以取值永远是0.所以010====n λλλΛ，这证明了n l l l ,,,10Λ是线性⽆关的.同时它⼀共包含1+n 个多项式，恰好等于[]X n K 的维数.所以n l l l ,,,10Λ构成了[]X n K 的⼀组基底.拉格朗⽇基本多项式作为基底的好处是所有的多项式都是齐次的（都是n 次多项式）.1.4优点与缺点拉格朗⽇插值法的公式结构整齐紧凑，在理论分析中⼗分⽅便，然⽽在计算中，当插值点增加或减少⼀个时，所对应的基本多项式就需要全部重新计算，于是整个公式都会变化，⾮常繁琐.这时可以⽤重⼼拉格朗⽇插值法或⽜顿插值法来代替.此外，当插值点⽐较多的时候，拉格朗⽇插值多项式的次数可能会很⾼，因此具有数值不稳定的特点，也就是说尽管在已知的⼏个点取到给定的数值，但在附近却会和“实际上”的值之间有很⼤的偏差.这类现象也被称为龙格现象，解决的办法是分段⽤较低次数的插值多项式.2 重⼼拉格朗⽇插值法重⼼拉格朗⽇插值法是拉格朗⽇插值法的⼀种改进.在拉格朗⽇插值法中，运⽤多项式)())(()(10k x x x x x x x l ---=Λ，图（2）拉格朗⽇插值法的数值稳定性：如图(2)，⽤于模拟⼀个⼗分平稳的函数时，插值多项式的取值可能会突然出现⼀个⼤的偏差（图中的14⾄15中间）可以将拉格朗⽇基本多项式重新写为：∏≠=--=kji i i j jj x x x x x l x l ,0)(1)()(，定义重⼼权∏≠=-=k ji i i j j x x ,0)(1ω，上⾯的表达式可以简化为：jjj x x x l x l -=ω)()(，于是拉格朗⽇插值多项式变为：j kj jjy xx x l x L ∑=-=0)()(ω，（1）即所谓的重⼼拉格朗⽇插值公式（第⼀型）或改进拉格朗⽇插值公式.它的优点是当插值点的个数增加⼀个时，将每个j ω都除以)(1+-k j x x ，就可以得到新的重⼼权1+k ω，计算复杂度为)(n O ，⽐重新计算每个基本多项式所需要的复杂度)(2n O 降了⼀个量级.将以上的拉格朗⽇插值多项式⽤来对函数1)(≡x g 插值，可以得到：∑=-=?kj jjx x x l x g x 0)()(,ω，因为1)(≡x g 是⼀个多项式. 因此，将)(x L 除以)(x g 后可得到：∑∑==--=k j jjk j jjx x x x x L 00)(ωω，（2）这个公式被称为重⼼拉格朗⽇插值公式（第⼆型）或真正的重⼼拉格朗⽇插值公式.它继承了（1）式容易计算的特点，并且在代⼊x 值计算)(x L 的时候不必计算多项式)(x l 它的另⼀个优点是，结合切⽐雪夫节点进⾏插值的话，可以很好地模拟给定的函数，使得插值点个数趋于⽆穷时，最⼤偏差趋于零.同时，重⼼拉格朗⽇插值结合切⽐雪夫节点进⾏插值可以达到极佳的数值稳定性.第⼀型拉格朗⽇插值是向后稳定的，⽽第⼆型拉格朗⽇插值是向前稳定的，并且勒贝格常数很⼩.3.分段线性插值对于分段线性插值，我们看⼀下下⾯的情况.3.1问题的重诉已知211)(xx g +=,66≤≤-x ⽤分段线性插值法求插值，绘出插值结果图形，并观察插值误差.1.在[-6，6]中平均选取5个点作插值；2.在[-6，6]中平均选取11个点作插值；3.在[-6，6]中平均选取21个点作插值；4.在[-6，6]中平均选取41个点作插值.3.2问题的分析在数值计算中，已知数据通常是离散的，如果要得到这些离散点以外的其他点的函数值，就需要根据这些已知数据进⾏插值.⽽本题只提供了取样点和原函数)(x g .分析问题求解⽅法如下：（1）利⽤已知函数式211)(xx g +=计算取样点X 对应的函数值Y ；将Y X ,作为两个等长的已知向量，分别描述采样点和样本值.因此被插值函数是⼀个单变量函数，可利⽤⼀维插值处理该数据插值问题.⼀维插值采⽤的⽅法通常有拉格朗⽇多项式插值（本题采⽤3次多项式插值），3次样条插值法和分段线性插值.（2）分别利⽤以上插值⽅法求插值.以0.5个单位为步长划分区间[-6，6]，并将每⼀点作为插值函数的取样点.再根据插值函数计算所选取样点的函数值.最后再利⽤所得函数值画出相应的函数图象，并与原函数)(x g 的图象进⾏对⽐.3.3问题的假设为了解决上述分析所提到的问题，本题可以作出如下假设：（1）假设原函数)(x g 仅作为求解取样点对应的样点值的函数关系式.⽽其他各点的函数值都是未知量，叙⽤插值函数计算.（2）为了得到理想的对⽐函数图象，假设)(x g 为已知的标准函数.可以选取0.5个单位为步长划分区间[-6，6]，分别计算插值函数和标准函数)(x g 在该区间的取样点的函数值.画出函数图象进⾏对⽐.3.4分段线性插值原理给定区间[]b a ,, 将其分割成b x x x a n =<<<=Λ10，已知函数)(x f y =在这些插值结点的函数值为),1,0)((n k x f y k k Λ==；求⼀个分段函数)(x I k ,使其满⾜：(1) k k h y x I =)(，),1,0(n k Λ=；(2) 在每个区间[]1,+k k x x 上, )(x I h 是个⼀次函数.易知,)(x I h 是个折线函数, 在每个区间[]1,+k k x x 上,),1,0(n k Λ=1111)(++++--+--=k kk kk k k k k h y x x x x y x x x x x I ,于是, )(x I h 在[]b a ,上是连续的，但其⼀阶导数是不连续的. 于是即可得到如下分段线性插值函数：)()(0x l y x I ni i i n ∑==,其中=≤≤--=≤≤--=+++---.,0;,;0,111111其他时舍去时，且当时舍去时，且当n i x x x x x x x i x x x xx x x l i i i i i i i i ii i3.5问题的求解在MATLAB 中实现分段线性插值，最近点插值，3次多项式插值，3次样条插值的命令为interp 1,其调⽤格式为: Y 1=interp 1(X ,Y ,X 1，’method ’)函数根据X ,Y 的值，计算函数在X 1处的值.X ,Y 是两个等长的已知向量，分别描述采样点和样本值，X 1是⼀个向量或标量，描述欲插值点，Y 1是⼀个与X 1等长的插值结果.method 是插值⽅法，包括：linear ：分段线性插值.它是把与插值点靠近的两个数据点⽤直线连接，然后在直线让选取对应插值点的数.nearest ：近点插值法.根据已知两点间的插值点与这两点间的位置远近插值.当插值点距离前点远时,取前点的值,否则取后点的值.cubic ：3次多项式插值.根据已知数据求出⼀个3次多项式，然后根据多项式进⾏插值. spline ：3次样条插值.在每个分段（⼦区间）内构造⼀个3次多项式，使其插值函数除满⾜插值条件外，还要求个节点处具有光滑条件.再根据已知数据求出样条函数后,按照样条函数插值.运⽤Matlab ⼯具软件编写代码，并分别画出图形如下： (⼀)在[-6，6]中平均选取5个点作插值：-10-5051000.20.40.60.81分段线性插值g(x)y1-10-50510-0.500.513次样条插值g(x)y2-10-5051000.20.40.60.81最近点插值g(x)y3-10-5051000.20.40.60.813次多项式插值g(x)y4（⼆）在[-6，6]中平均选取11个点作插值：-10-5051000.20.40.60.81-10-5051000.20.40.60.81-10-5051000.20.40.60.81最近点插值-10-5051000.20.40.60.813次多项式插值g(x )y1g(x )y2g(x )y3g(x )y4（三）在[-6，6]中平均选取21个点作插值：-10-5051000.20.40.60.81分段线性插值-10-551000.20.40.60.813次样条插值-10-551000.20.40.60.81最近点插值-10-551000.20.40.60.813次多项式插值g(x )y1g(x )y2g(x )y3g(x )y4（四）在[-6，6]中平均选取41个点作插值-10-5051000.20.40.60.81g(x )y1-10-5051000.20.40.60.81g(x )y2-10-5051000.20.40.60.81最近点插值g(x )y3-10-5051000.20.40.60.813次多项式插值g(x )y43.6 分段插值⽅法的优劣性分析从以上对⽐函数图象可以看出，分段线性插值其总体光滑程度不够.在数学上，光滑程度的定量描述是函数(曲线) 的k 阶导数存在且连续，则称该曲线具有k 阶光滑性.⼀般情况下，阶数越⾼光滑程度越好.分段线性插值具有零阶光滑性，也就是不光滑.3次样条插值就是较低次数的多项式⽽达到较⾼阶光滑性的⽅法.总体上分段线性插值具有以下特点：优点: 1.分段线性插值在计算上具有简洁⽅便的特点.2.分段线性插值与3次多项式插值函数在每个⼩区间上相对于原函数都有很强的收敛性，（舍⼊误差影响不⼤），数值稳定性好且容易在计算机上编程实现等优点缺点: 分段线性插值在节点处具有不光滑性的缺点(不能保证节点处插值函数的导数连续)，从⽽不能满⾜某些⼯程技术上的要求.⽽3次样条插值却具有在节点处光滑的特点.。
插值方法比较
1. 克里金法（Kriging）克里金法是通过一组具有z 值的分散点生成估计表面的高级地统计过程。
与其他插值方法不同，选择用于生成输出表面的最佳估算方法之前应对由z 值表示的现象的空间行为进行全面研究。
克里金插值与IDW插值的区别在于权重的选择，IDW仅仅将距离的倒数作为权重,而克里金考虑到了空间相关性的问题。
它首先将每两个点进行配对，这样就能产生一个自变量为两点之间距离的函数。
对于这种方法，原始的输入点可能会发生变化.在数据点多时，结果更加可靠。
该方法通常用在土壤科学和地质中。
2. 反距离权重法(Inverse Distance Weighted，IDW）反距离权重法（反距离权重法）工具所使用的插值方法可通过对各个待处理像元邻域中的样本数据点取平均值来估计像元值。
点到要估计的像元的中心越近，则其在平均过程中的影响或权重越大.此方法假定所映射的变量因受到与其采样位置间的距离的影响而减小。
例如,为分析零售网点而对购电消费者的表面进行插值处理时,在较远位置购电影响较小，这是因为人们更倾向于在家附近购物.反距离权重法主要依赖于反距离的幂值。
幂参数可基于距输出点的距离来控制已知点对内插值的影响。
幂参数是一个正实数，默认值为2。
通过定义更高的幂值,可进一步强调最近点.因此，邻近数据将受到最大影响，表面会变得更加详细（更不平滑）。
随着幂数的增大，内插值将逐渐接近最近采样点的值。
指定较小的幂值将对距离较远的周围点产生更大影响，从而导致更加平滑的表面。
由于反距离权重公式与任何实际物理过程都不关联，因此无法确定特定幂值是否过大。
作为常规准则，认为值为30 的幂是超大幂，因此不建议使用.此外还需牢记一点，如果距离或幂值较大，则可能生成错误结果。
3。
含障碍的样条函数（Spline with Barriers)含障碍的样条函数工具使用的方法类似于样条函数法工具中使用的技术，其主要差异是此工具兼顾在输入障碍和输入点数据中编码的不连续性.含障碍的样条函数工具应用了最小曲率方法，其实现方式为通过单向多格网技术，以初始的粗糙格网(在本例中是已按输入数据的平均间距进行初始化的格网)为起点在一系列精细格网间移动，直至目标行和目标列的间距足以使表面曲率接近最小值为止。
常见插值方法及其介绍
常见插值方法及其介绍常见的插值方法有最邻近插值、双线性插值、双三次插值和基于样条的插值方法。
下面将对这些方法进行介绍。
1.最邻近插值：最邻近插值是最简单也是最直观的插值方法之一、该方法的原理是将待插值点附近最近的一个已知像素的灰度值赋给待插值点。
这种插值方法的优点是计算简单且实时性好，但缺点是结果较为粗糙，会出现明显的锯齿状边缘。
2.双线性插值：双线性插值是一种基于线性插值的方法，它考虑了待插值点附近四个已知像素的灰度值来生成新的像素值。
具体而言，对于一个待插值点，首先在水平方向上计算它上下两个已知像素的插值，然后在竖直方向上计算其左右两个已知像素的插值，最后再在这两次插值的基础上进行一次线性插值。
这种插值方法的优点是计算相对简单，效果较好，但仍然会存在锯齿状边缘。
3.双三次插值：双三次插值是一种更为复杂的插值方法，它通过分析待插值点周围的16个已知像素的灰度值来生成新的像素值。
具体而言，双三次插值首先根据已知像素的位置与待插值点的距离计算出一个权重系数矩阵，然后将这个系数矩阵与对应的已知像素灰度值相乘并相加。
这种插值方法的优点是结果较为平滑，点缺失问题较少，但计算量较大。
4.基于样条的插值方法：基于样条的插值方法主要包括线性样条插值、三次样条插值和B样条插值。
这些方法是基于插值函数的一种改进，通过选取合适的插值函数形式来拟合已知像素点，从而实现待插值点的灰度值推测。
这些方法计算量较大，但插值效果相对较好，具有高度灵活性。
总结：常见的插值方法包括最邻近插值、双线性插值、双三次插值和基于样条的插值方法。
最邻近插值计算简单且实时性好，但结果较为粗糙；双线性插值效果较好，但仍然存在锯齿状边缘；双三次插值平滑度较高，但计算量较大；基于样条的插值方法具有高度灵活性，但计算量较大。
选择适合的插值方法需根据具体需求考虑。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
1. 克里金法（Kriging）
克里金法是通过一组具有z 值的分散点生成估计表面的高级地统计过程。
与其他插值方法不同，选择用于生成输出表面的最佳估算方法之前应对由z 值表示的现象的空间行为进行全面研究。
克里金插值与IDW插值的区别在于权重的选择，IDW仅仅将距离的倒数作为权重，而克里金考虑到了空间相关性的问题。
它首先将每两个点进行配对，这样就能产生一个自变量为两点之间距离的函数。
对于这种方法，原始的输入点可能会发生变化。
在数据点多时，结果更加可靠。
该方法通常用在土壤科学和地质中。
2. 反距离权重法（Inverse Distance Weighted，IDW）
反距离权重法（反距离权重法）工具所使用的插值方法可通过对各个待处理像元邻域中的样本数据点取平均值来估计像元值。
点到要估计的像元的中心越近，则其在平均过程中的影响或权重越大。
此方法假定所映射的变量因受到与其采样位置间的距离的影响而减小。
例如，为分析零售网点而对购电消费者的表面进行插值处理时，在较远位置购电影响较小，这是因为人们更倾向于在家附近购物。
反距离权重法主要依赖于反距离的幂值。
幂参数可基于距输出点的距离来控制已知点对内插值的影响。
幂参数是一个正实数，默认值为2。
通过定义更高的幂值，可进一步强调最近点。
因此，邻近数据将受到最大影响，表面会变得更加详细（更不平滑）。
随着幂数的增大，内插值将逐渐接近最近采样点的值。
指定较小的幂值将对距离较远的周围点产生更大影响，从而导致更加平滑的表面。
由于反距离权重公式与任何实际物理过程都不关联，因此无法确定特定幂值是否过大。
作为常规准则，认为值为30 的幂是超大幂，因此不建议使用。
此外还需牢记一点，如果距离或幂值较大，则可能生成错误结果。
3. 含障碍的样条函数（Spline with Barriers）
含障碍的样条函数工具使用的方法类似于样条函数法工具中使用的技术，其主要差异是此工具兼顾在输入障碍和输入点数据中编码的不连续性。
含障碍的样条函数工具应用了最小曲率方法，其实现方式为通过单向多格网技术，以初始的粗糙格网（在本例中是已按输入数据的平均间距进行初始化的格网）为起点在一系列精细格网间移动，直至目标行和目标列的间距足以使表面曲率接近最小值为止。
4. 地形转栅格（Topo to Raster）
地形转栅格和依据文件实现地形转栅格工具所使用插值技术是旨在用于创建可更准确地表示自然水系表面的表面，而且通过这种技术创建的表面可更好的保留输入等值线数据中的山脊线和河流网络。
5. 样条函数（Spline）
样条函数法工具所使用的插值方法使用可最小化整体表面曲率的数学函数来估计值，以生成恰好经过输入点的平滑表面。
从概念上讲，采样点被拉伸到它们数量上的高度；样条函数折弯一个橡皮页，该橡皮页在最小化表面总曲率的同时穿过这些输入点。
在穿过采样点时，它将一个数学函数与指定数量的最近输入点进行拟合。
此方法最适合生成平缓变化的表面，例如高程、地下水位高度或污染程度。
IDW 插值主要受幂指数和各采样点属性值变化情况的影响，幂指数越高，其局部影响的程度越高，在IDW搜索半径内，若各个采样点属性值变化较小时，内插结果受幂指数的影响较小；Spline 插值主要受插值类型（Regularized 或Tension）和weight 值的影响，一般Regularize 插值结果比Tension插值结果光滑，在Regularized Spline 插值中，weight 值越高生成的表面越光滑，Tension Spline 插值则相反；总体来看，IDW和SPLINE 插值受采样点范围、采样点密度、采样点属性取值变化以及各自的参数影响，当采样点足够密时，使用IDW 插值可以取得良好效果，SPLINE插值则适合那些空间连续变化且光滑的表面的生成。
6. 自然邻域法（Natural Neighbor）
自然邻域法插值可找到距查询点最近的输入样本子集，并基于区域大小按比例对这些样本应用权重来进行插值（Sibson，1981）。
该插值也称为Sibson 或“区域占用(area-stealing)”插值。
该插值方法的基本属性是它具有局部性，仅使用查询点周围的样本子集，且保证插值高度在所使用的样本范围之内。
该插值方法不会推断趋势且不会生成输入样本尚未表示的山峰、凹地、山脊或山谷。
该表面将通过输入样本且在除输入样本位置之外的其他所有位置均是平滑的。
原理是构建voronoi多边形，也就是泰森多边形。
首先将所有的空间点构建成voronoi多边形，然后将待求点也构建一个voronoi多边形，这样就与圆多边形有很多相交的地方，根据每一块的面积按比例设置权重，这样就能够求得待求点的值了。
7. 趋势（trend）
趋势面法是一种可将由数学函数（多项式）定义的平滑表面与输入样本点进行拟合的全局多项式插值法。
趋势表面会逐渐变化，并捕捉数据中的粗尺度模式。
使用趋势插值法可获得表示感兴趣区域表面渐进趋势的平滑表面。
此种插值法适用于以下几种情况
（1）感兴趣区域的表面在各位置间出现渐变时，可将该表面与采样点拟合，例如，工业区的污染情况。
（2）检查或排除长期趋势或全局趋势的影响。
此类情况下，采用的方法通常为趋势面分析。