几种常用的插值方法.pptxVIP

几种常用的插值方法 数学系 信息与计算科学 1 班 李平 指导老师：唐振先 摘要:插值在诸如机械加工等工程技术和数据处理等科学研究中有许多直接的应用,在很多 领域都要用插值的办法找出表格和中间值,插值还是数值积分微分方程数值解等数值计算的 基础。本文归纳了几种常用的插值方法,并简单分析了其各自的优缺点。 关键词:任意阶多项式插值,分段多项式插值。 引言:所谓插值,通俗地说就是在若干以知的函数值之间插入一些未知函数值,而 插值函数的类型最简单的选取是代数多项式。用多项??建立插值函数的方法主要 用两种：一种是任意阶的插值多项式,它主要有三种基本的插值公式：单项式,拉 格朗日和牛顿插值；另一种是分段多项式插值,它有Hermite 和 spine 插值和分段 线性插值。 一．任意阶多项式插值： 用单项式基本插值公式进行多项式插值： 多项式 插值是 求通 过几个 已知 数据点 的那个 n-1 阶多项 式, 即 Pn-1(X)=A1+A2X+…AnXn-1,它是一个单项式基本函数 X0,X1…Xn-1 的集合来定义多 项式,由已知 n 个点（X,Y）构成的集合,可以使多项式通过没数据点,并为 n 个未 知系数 Ai 写出n 个方程,这n 个方程组成的方程组的系数矩阵为 Vandermonde 矩 阵。 虽然这个过程直观易懂, 但它都不是建立插值多项式最好的办法, 因为 Vandermonde 方程组有可能是病态的,这样会导致单项式系数不确定。另外,单项 式中的各项可能在大小上有很大的差异,这就导致了多项式计算中的舍入误差。 拉格朗日基本插值公式进行插值：;n.容易看出 n 次多项式 Li（x）满足 Li（x）=1，（i=j）；Li（ x ）=0，（i≠j），其;结点时,它不必重新计算,只要在后面再计算一项均插即可,减少了计算量,不用计 算全部系数,节约了大量人力,物力,财力。 增加插值多项式的阶数并不一定能增加插值的精度,据定义,插值式,F(x)可以与 结点(xi，yi)，i=1，…，n 处的实际函数匹配,但却不能保证支点之间求F(x),还能 很好的逼近产生(xi，yi)数据的实际函数 F(x)。例如,如果 F(x)为一个已知的解析 函数,而且定义 F(x)的节点集合中数据点的数目可以增加（多项式 F(x)的阶数也 增加）,但是,由于 F(x)的起伏增加,那么插值式就可能在节点见振带,基于当实际函 数F(x)平滑时,这种多项式摆动也可能发生,这种振荡不是由多项式摆动引起的,而 是由多项式的项相加来求插值多项式时发生舍入误差造成的。有时多项式摆动可 通过谨慎选择基础函数的取样来成为,但如果数据是由不容易重复实验取得的,就 不能这么做了,这会司会用下面介绍分段插值法。 二、分段插值多项式 1、分段线性插值： 分段线性插值最简单的插值方案,只要将每个相邻的节点用直线接起来,如此形 成的一条新的折线就是分段线性插值函数,记作In( x j)=yi 而且In(x)在每个区间[ x j x j+1]上是线性函数（j=0,1…n-1）;4;;在[a,b]上取 n+1 个插值结点 a=x0x1…xn=b 已知函数 f(x)在这 n+1 个点的函数 值为yk=f(xk)则在[a,b]上函数y=f(x)的m 次样条插值函数S(x)满足: (1)S(x)在(a,b)上直到m-1 阶导数连续; (2)S(xk)=yk,(k=0 1…n); (3)在区间[xk,xk+1](k=0 1 …n-1)上,S(x)是m 次多项式。 2.三次样条函数 在[a,b]上函数y=f(x)的三次样条插值函数S(x)满足: (1)在(a,b)上 0、1、2 阶导数连续,即: s(xk-0)=s(xk+0)，s″(xk-0)=s″(xk+0) (k=0 1…n-1) (2)S(xk)=yk (k=0，1，…n); (3)在区间[xk xk+1](k=0，1…n-1)上S(x)是三次多项式。 3.三次样条函数的计算 由二阶导数连续,设 s″(xk )=mk,(k=0,1, …,n),mk 是未知待定的数。因 S(x)是分段 三次多项式,则在每个区间[xk xk+1]内,S″(x)是分段一次多项式,记 hk=xk+1-xk 则:;而s(xk+0)=s(xk-0)联立上述两式,得到关于m 的方程: