计算方法第1章插值方法.pptVIP

第一章　插值方法 拉格朗日插值 牛顿插值 埃特金插值 存在唯一性定理 插值余项 分段三次埃米尔差值 三次样条插值 1.1 拉格朗日插值公式 拉格朗日插值公式的基本思想是:把Pn(x)的构造问题转化为n+1个插值基函数Li(x)（i=0,1,…,n)的构造。插值函数Pn(x)是这n＋1插值基函数的线性组合，其组合系数就是对应点上的函数值。 1.２牛顿插值公式 　牛顿插值公式的基本思想是:把Pn(x)设计为递推形式。 　Pn(x)＝Pn-1(x)+an(x-x0)…(x-xn-1) =a0+a1(x-x0)+a2(x-x0)(x-x1)+ …  + an(x-x0)…(x-xn-1) 满足Pn(xi)=f(xi)=Yi i=0,1,2,…,n 称为牛顿插值公式。 1.２牛顿插值公式 为了简化上式，引进记号： ?0(x)=1 ?i(x)=(x-xi-1) ?i-1(x) =(x-x0)(x-x1)…(x-xi-1)  i=1,2 …，n Pn(x)＝a0 ?0(x) +a1 ?1(x) + …+ an ?n(x) ?0(x) ,?1(x) ,… ,?n(x)称为牛顿插值以x0, x1,  …, xn为插值节点的基函数。 如果增加一个新节点，只需增加一个新项 an ?n(x) 1.２牛顿插值公式 其系数为各阶差商： 1.3埃特金插值公式 埃特金插值公式的基本思想是：平面上的两个点可以连成一条线，对应一个线性函数；把线性函数看作形式点，经线性组合，可以构成二次函数；把二次函数看作形式点，经线性组合，可以构成三次函数。因次埃特金插值公式也叫逐步线性插值。 过两点（X0,Y0),(X1,Y1）的一次线性插值函数为: 1.3埃特金插值公式 过两点（X0,Y0),(X2,Y2）的一次线性插值函数为: 1.3埃特金插值公式 对于一般情况：过两个形式点（Xn-1,P0,1,2…,n-1（x））和（Xn,P0,1,…n-2,n (x)）的线性插值函数为: 1.3埃特金插值公式 1.4 存在唯一性定理 定理1 有惟一的n次多项式Pn(x)，满足条件： Pn(xi)=yi （i=0,1,…,n) 1.5 插值余项 设Pn(x)是在点X0,X1,…，Xn处关于f(x)的插值多项式，当X?Xi时，f(x)与 Pn(x)的偏差Rn(x)＝ f(x)－Pn(x)，称为插值余项。 定理 2 若f(x)在包含插值节点X0,X1,…，Xn的区间［a,b]上n+1次可微分，则对任意x,x?[a,b],有与x有关的?(a?b)存在，使得： Rn(x)= f(x)－Pn(x)＝ 其中?(x)=(x-x0)(x-x1)…(x-xn) 输入插值节点xi(0,…n) 计算插值节点上的原函数yi(0,…n) 输入插值点的个数m h=(xn-x0)/(m-1) 计算插值点上的插值基函数l、插值函数ll(k)和原函数yy(k) for k=0 to m xx(k)=x0+i*h ll(k)=0 for i=0 to n l=1 for j=0 to n if ji l=l*(xx(k)-x(j))/(x(i)-x(j)) end if next j ll(k)=ll(k)+l*y(i) next I yy(k)=1/(1+xx(k)*xx(k)) next k 确定坐标 绘制图形 1.6 分段三次埃尔米特插值 1.6 分段三次埃尔米特插值 已知xi,f(xi),f′(xi)(i=0,1,2,…,n)，求三次插值函数H(x)满足： H(xi)=f(xi)， H′(xi)=f ′(xi) (i=0,1,2,…,n) 在任意的子区间[Xi,Xi+1], i?(0,1,2, ,…,n-1)上构造插值函数 H(x)=f(xi)h1(x)+f(xi)h2(x)+ f ′(xi)h3(x)+f ′(xi+1)h4(x) hk(x) (k=1,2,3,4)称为埃尔米特插值的插值基函数。 根据插值函数H(x)所满足的条件，可以