插值法与数值微分教学幻灯片讲义.ppt

第二章 插值法与数值微分;求近似函数的方法:由实验或测量的方法得到所求函数 y=f(x) 在互异点x0 , x1, ... , xn 处的值 y0 , y1 , … , yn , 构造一个简单函数 ?(x) 作为函数 y=f(x) 的近似表达式 y= f(x) ? ?(x) 使 ?(x0)=y0 , ?(x1)=y1 , ?, ?(xn)=yn , (a) 这类问题称为插值问题。 f(x) 称为被插值函数，?(x) 称为插值函数， x0 , x1, ... , xn 称为插值节点。 (a)式称为插值条件。常用的插值函数是多项式。 插值的任务就是由已知的观测点,为物理量(未知量)建立一个简单的、连续的解析模型，以便能根据该模型推测该物理量在非观测点处的特性。; ; ; 2.1 线性插值 先从最简单的线性插值(n=1)开始。这时插值问题(2)就是求一次多项式 P1(x)=a0+a1x 使它满足条件 P1(x0)=y0 , P1(x1)=y1 , 令P1(x)=l0(x)y0+l1(x)y1 ,由于 l0(x0)=1, l0(x1)=0, l1(x0)=0, l1(x1)=1. 这样l0(x)含有因子x-x1, 令 l0(x)=λ(x-x1), 再利用 l0(x0)=1确定其中的系数，结果得到 x-x1 l0(x)=------------ , x0-x1 类似的可得到 x-x0 l1(x)=------------ , x1-x0 这样 x-x1 x-x0 P1(x)=---------y0 + --------y1 , x0-x1 x1-x0 l0(x), l1(x)称为以x0 , x1 为节点的插值基函数。;Newton插值;Newton插值;因此, 是f(x) 在 x1 , x0处的一阶 均差， 利用均差的对称性， （＊） 式可以表示为： 这种形式的插值叫做牛顿（Newton）插值． ;;; 构造l0(x):由于l0(x)有x1和x2二个零点，因此有因子(x-x1)(x-x2)，又因有l0(x）是一个次数不高于二次的多项式，所以，还可能相差一个常数因子，于是把l0(x）写成：;x0; 利用条件 ，可求得： 那么x0点的二次插值基函数为： 同理构造 为： ;插值函数： 牛顿二次插值多项式： 假设过 的二次插值多项式具有下面的形式： 确定系数A,B,C： 利用 我们有： ;再利用 有： B = 最后 确定C： （为一阶均差的均差） 其中 ;f [xi, xj , xk ]为f (x)在点xi, xj, xk处的二阶差商 一般的称 为f (x)在点x0, x1,…, xn处的n阶差商。 于是得到二次牛顿插值多项式： ; 三、例题： 例1: 已给sin0.32=0.314567,sin0.34=0.333487, sin0.36=0.352274, 用线性插值及抛物插值计算 sin0.3367 的值并估计截断误差。 解： 由题意取x0=0.32, y0=0.314567 , x1=0.34 , y1=0.333487 , x2=0.36 , y2=0.352274 。 用线性插值及抛物插值计算，取 x0=0.32 及 x1=0.34 , 又由公式得 ;其截断误差为 其中 ，因 f(x)=sinx，f//(x)= -sinx， 可取 ，于是 ?R1(0.3367)?=?sin 0.3367 –P1(0.3367)? ?1/2(0.3335)(0.0167)(0.0033)?0.92?10–5， 若取x1=0.34，x2=0.