实验名称 插值法2.doc

探索实验5 插值法 实验目的 了解插值问题及其适用的场合，理解并掌握常用的插值算法的构造和计算，了解差商概念、Runge现象及样条插值方法，学习用计算机求近似函数的一些科学计算方法和简单的编程技术。 二、概念与结论 1. 插值问题与插值函数： 由实验或测量的方法得到所求函数 y=f(x) 在互异点x0 , x1, ... , xn 处的值 y0 , y1 , … , yn ,构造一个简单函数 ((x) 作为函数 y=f(x) 的近似表达式 y= f(x) ( ((x) 使 ((x0)=y0 , ((x1)=y1 , (, ((xn)=yn , (1) 这类问题称为插值问题。 f(x) 称为被插值函数，((x) 称为插值函数， x0 , x1, ... , xn 称为插值节点。(1)式称为插值条件。 常用的插值函数是多项式函数。且当n=1时是称为线性插值，n=2时称为Simpson插值或抛物线插值。 2．插值定理： 假设x0 ,x1,…,xn 是n+1个互异节点,函数f(x)在这组节点的值f(xk)(k=0,1,…,n)是给定的，那么存在唯一的n 次次多项式pn (x)满足 pn (xk)=f(xk), k=0,1,…,n 3.插值的截断误差 设(n(x)是过点x0 ，x1 ，x2 ，…xn的 n 次插值多项式， f(n+1)(x)在（a，b）上存在，其中[a，b]是包含点x0 ，x1 ，x2 ，…,xn的任一区间，则对任意给定的x([a，b]，总存在一点((（a，b）（依赖于x）使 其中 (n+1 (x)=(x –x0) (x - x1)…(x-xn ) ，f(n+1)(() 是f(x)的n+1阶微商在 ( 的值。 4. 差商： 　给定一个函数表 x | x0 x1 ….... xn -------- --------------------------------------------------------- y | y0 ,y1 ……. yn 其中当i (j时 ,xi(xj 记 f[xi]=f(xi) ,定义f(x)关于xi,xj的一 阶差商 一般的, f(x)关于xi,xi+1,…,xi+k的k 阶差商定义为: 5. 分段线性插值： 在区间[a,b]上给定一组节点: a=x0x1x2(((xn=b 且在节点上的函数值f(xi )=yi, i=0,1,2,(((,n, 如果函数((x) 满足： ( (x)在[a , b]上连续 ( (xi )=yi, i=0,1,2,(((,n ( (x)在每个子区间[xi , xi+1](i=0,1,2,(((,n-1)上是次数为1的多项式 则称( (x)是定义在[a ,b]上的分段线性插值函数。 6. m次样条插值函数： 设f(x)是区间[a,b]上的一个连续可微函数,在区间[a,b]上给定一组节点: a=x0x1x2(((xn=b 且y = f(x)在点 x0，x1，x2, ((( xn的值为y0，y1，y2, (((yn，如果函数s(x)满足条件 (1) s(x)在每个子区间[xi , xi+1](i=0,1,2,(((,n-1)上是次数不超过m的多项式; s(x)在区间[a , b]上有m-1阶连续导数; s(xi)=f(xi) =yi , i=0,1,2,(((,n 则称s(x) 为函数f(x)的在[a ,b]上的m次样条插值函数。x0，x1，x2, (((称为样条结点,其中x1,(((,xn-1称为内结点, x0 , xn 称为边界结点。当m=3时,则称s(x)为函数f(x)的三次样条插值函数。 三、程序中Mathematica语句解释： Product [f(i) , {i ，imin，imax，h }] 计算乘积 f(imin)(f(imin +h) (f(imin +2h) (……(f(imin +nh)} imax – h ( imin + nh ( imax , h0 2．Simplify[expr] 对expr进行化简 3. Append[list, elem] 在表list的后面加一个元素 elem 4. Prepend[list, elem] 在表list的前面加一个元素 elem 四、方法与程序 插值法是函数逼近的重要方法之一,