计算方法-4插值方法.pptVIP

1 分段线性插值 定义：设f(x)在给定的n+1个互异节点上的函数值为f(xi)=yi,(i=0,1,2,…,n),S1(x)为区间[a,b]上的函数，如果S1(x)满足下列条件： S1(x)是在区间内地每个小区间[xi, xi+1], (i=0,1,2,…,n)的线性函数； S1(xi)=yi, (i=0,1,2,…,n) S1(xi)是插值区间[a,b]上的连续函数 则称S1(x)是插值区间[a,b]上的分段线性插值函数 分段线性插值 几何意义：在每个小区间[xi, xi+1], (i=0,1,2,…,n) 上，作连接点(xi,yi)与(xi+1,yi+1)的直线。 函数S1(x)在区间[xi-1,xi) (i=0,1,2,…,n)上的表达式为： 函数S1(x)在区间[xi-1,xi) (i=0,1,2,…,n)上的表达式为 分段线性插值 如果用基函数表示， S1(x)在插值区间[a,b]上的表达式为： 分段线性插值余项 定理：设f(x)在给定的互异节点xi上的函数值为f(xi)=yi,(i=0,1,2,…,n), f(x)属于C1[a,b], f’’(x)在[a,b]上存在， S1(x)为插值区间[a,b]上的、由数据(xi,yi)(i=0,1,2,…,n)构成的分段线性插值函数。则, 分段线性插值余项 例：设f(x)=1/(1+x2),在[-5,5]上取n=10，按等间距节点求线性插值函数S1(x)与f(x)的值，并估计误差。 解： h=(5-(-5))/10=1 在区间[xi-1,xi] (i=0,1,2,…,n)上的表达式为 分段线性插值余项 分段线性插值余项 4.5 样条插值函数 高次插值：计算复杂,可能出现不收敛现象,Runge现 象; 分段线性插值: 计算简单,一致收敛,但光滑性不足; 样条函数插值: 分段多项式,光滑连接 4.5.1 三次样条插值函数 定义：设f(x)在区间[a,b]上给定n+1个节点a=x0x1…xn=b.若函数S(x)满足: S(x)在每一个小区间 [xi, xi+1], (i=0,1,2,…,n-1)上是三次多项式； S(x)在[a,b]上有连续二阶导数. 则称S(x)为三次样条函数 三次样条插值函数 设三次样条函数表达式为: 其中系数A,B,C,D为待定,需要满足的条件为: 三次样条插值函数 定义：设f(x)在区间[a,b]上给定n+1个节点a=x0x1…xn=b及函数y=f(x)在这些节点上的值yi=f(xi)(i=0,1,…,n).若3次样条函数S(x)满足一下插值条件: 则称S(x)为三次样条插值函数 三次样条插值函数 求解条件：需求4n个未知数,目前有4n-2个条件,通常还需加上边界条件,通常有下面两种: 当n大时,计算量非常大 4.5.2三次样条插值函数求解 求解条件：需求4n个未知数,目前有4n-2个条件,通常还需加上边界条件,通常有下面三种: 已知端点一阶导数 已知端点二阶导数 特别地,若 称为自然边界条件 周期性边界条件 三次样条插值函数求解 例：已知函数f(x)在三个点处的值为f(-1)=1,f(0)=0, f(1)=1,在区间[-1,1]上,求f(x)在自然边界条件下的三次样条插值函数 三次样条插值函数求解 三次样条插值函数求解 三次样条插值函数求解 三次样条插值函数求解 三次样条插值函数求解 三次样条插值函数求解 三次样条插值函数求解 三次样条插值函数求解 三次样条插值函数求解 三次样条插值函数求解 * * * * * * * * * * * * * 为了确定未知数，还需要补充两个边界条件。下面分两种边界条件讨论： (1)若要求满足第一边界条件： 方程4.2可以转化为n-1阶方程组 (2)若要求S(x)满足第二边界条件： 分别在[x0,x1]和[xn-1,xn]中建立关于m0,m1,…,mn的方程。 在式4.1中，令i=0，x=x0 在式4.1中，令i=n-1，x=xn 将上面两式联立4.2，得到n+1阶方程组： 其中di与4.3式相同，d0和dn的值为： 若要求S(x)满足第二边界条件： 分别在[x0,x1]和[xn-1,xn]中建立关于m0,m1,…,mn的方程。 在式4.1中，令i=0，x=x0 在式4.1中，