实验名称插值法技术总结.docVIP

探索实验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 四、方法与程序 插值法是函数逼近的重要方法之一,