数值计算方法课件.pptxVIP

第1章插值方法;1.1拉格朗日插值公式;图1－1插值多项式;1.n=1的情况

已知函数y=f(x)在点x0,x1上的值为y0,y1，要求多项式y=p1(x)，使p1(x0)=y0,p1(x1)=y1。其几何意义，就是通过两点A(x0,y0),B(x1,y1)的一条直线，如图1－2所示。;图1－2一次插值多项式;由直线两点式可知，通过A，B的直线方程为

它也可变形为p1(x)=l0(x)y0+l1(x)y1

显然有：l0(x0)=l1(x1)=1，l0(x1)=l1(x0)=0，p1(x0)=y0，p1(x1)=y1;我们称l0(x)为点x0的一次插值基函数，l1(x)为点x1的一次插值基函数。它们在对应的插值点上取值为1，而在另外的插值点上取值为0。插值函数p1(x)是这两个插值基函数的线性组合，其组合系数就是对应点上的函数值。这种形式的插值称作为拉格朗日（Lagrange）插值。;2.n=2的情况

线性插值只利用两对值(x0,y0)及(x1,y1)求得y=f(x)的近似值，误差较大。

p2(x0)=y0,p2(x1)=y1,p2(x2)=y2

p2(x)是x的二次函数，称为二次插值多项式。通过三点的插值问题称为二次插值或抛物插值。; 3.一般情况

我们看到，两个插值点可求出一次插值多项式p1(x)，而三个插值点可求出二次插值多项式p2(x)。当插值点增加到n+1个时，我们可以利用Lagrange插值方法写出n次插值多项

式pn(x)，如下所示：

;1.2牛顿插值公式;Newton插值算法如下：

inputx,(xi,yi),i=0,1,…,n。

y=y0,t=1。

forj=1,…,ndo

t=t*(x-xj-1)

fori=0,…,n-jdo



end

y=y+y0*t

end

output(x,y),(xi,yi),i=0,1,…,n。;Newton插值算法中的j循环由三部分组成：计算(x-xj)的累积，存入t单元；内套一个i循环用来依次计算差商表中的各阶差商，存入yi单元；y单元用于存放Newton公式中各项累加之和。;［例3］已知f(-1)=2,f(1)=1,f(2)=1，求f(x)的Newton插值多项式。

解：设x0=-1,x1=1,x2=2,则;1.3埃特金插值公式;x

; 从Aitken插值公式向算法转化要考虑的问题是：

(1)插值公式右端n-1次多项式应如何处理；

(2)插值表中的元素应设置多少个存储单元；

(3)插值表中第k列第i行元素的计算公式。?; Aitken插值算法如下：

 inputx,(xi,yi),i=0,1,…,n

1k

L：

fori=k,k+1,…,ndo

 end

 ifk≠nthenk+1k,gotoL

 ifk=n,outputyn

;Aitken插值算法为二重循环。外循环为k循环，用于计算Aitken插值表中的第k列；内循环为i循环，用于计算Aitken插值表中的第k列中的第i个元素。;［例4］已知f(-1)=2,f(1)=1,f(2)=1,求f(x)的Aitken插值多项式。

解：设x0=-1,x1=1,x2=2;x;1.4存在惟一性定理;1.5插值余项;［例５］设f(x)=lnx,并假定已给出值表试近似计算ln(0.6)的值，并指出精度。

解：利用3次Lagrange插值公式,简单计算过程如下：

;;综合上述,我们有

真值：ln(0.6)=-0.510826，

近似值：p3(0.6)=-0.509975，

真误差：ln(0.6)-p3(0.6)=-0.000851， 估计的上界：|ln(0.6)-p3(0.6)|0.00391;［例６］给定 (x∈［-5,5］)。取等距节点xi=-5+i(i=0,1,…,10),试建立插值多项式L10(x),并作图形,观察L10(x)对f(x)的逼近效果。;图1-3例6的图形;1.6分段三次埃尔米特插值;已知xi,f(xi),f(xi)(i=0,1,…,n),求分段三次插值函数H(x)满足

H(xi)=f(xi),H(xi)=f(xi) (i=0