第6章节 插值函数.pptVIP

第6章 拉格朗日插值多项式 设函数y=f(x)在区间[a,b]上有定义，且已知在点a ≤x0x1…xn ≤b上的值y0,y1,…,yn,若存在一简单函数P(x),使 P（xi）=yi (i=0,1,2,…,n) （6.1.1）成立，就称P(x)为f(x)的插值函数，点x0，x1，…，xn称为插值节点，包含插值节点的区间[a,b]称为插值区间，求插值函数P(x)的方法称为插值法。若P(x)为插值多项式，即 P(x)=a0+a1x+…+anxn （6.1.2）其中ai为实数，就称P(x)为插值多项式，相应的插值法称为多项式插值。若P(x)为三角多项式，就称三角插值。从几何图形上看，插值法就是求曲线y=P(x),使其通过给定的n+1个点（xi,yi),i=0,1,…,n,并用它近似已知曲线y=f(x)。 6.1 插值多项式存在的唯一性 设P(x)是形如（6.1.2）的插值多项式，则满足条件（6.1.1）的插值多项式（6.1.2）存在唯一。6.2 线性插值与抛物插值 要求插值多项式P(x),可以通过方程组 ∑ aixi=yi的解a0,a1,…an得到,但这样计算复杂,且难于得到P(x)的简单表达式。为此我们先讨论n=1的情况。假定已知区间[XK,,XK+1]的端点处的函数值yk=f(xk), yk+1=f(xk+1),要求线性插值多项式L1(x),使它满足 L1(xk)=yk L1(xk+1)=yk+1y= L1(x)的几何意义就是通过两点[xk,yk]与[xk+1,yk+1]的直线。 L1(x)的表达式可由几何意义直接给出： L1(x)=(（xk+1-x）/ （xk+1-xk）)*yk+(（x-xk）/ （xk+1-xk）)*yk+1 由上面表达式可以看出，L1(x)是由两个线性函数 lk(x)= （x-x k+1）/ （xk-xk+1）和 lk+1(x)= （x-xk）/ （xk+1-xk）的组合得到的，其系数分别为yk及yk+1,即 L1(x)=yk lk(x)+yk+1 lk+1(x)显然，lk(x)及 lk+1(x)也是线性插值多项式，在节点xk及xk+1上满足条件 lk(xk)=1, lk(xk+1)=0 lk+1(xk)=0, lk+1(xk+1)=1我们称lk(x)及 lk+1(x)为线性插值基函数。 下面讨论n=2的情况。假定插值节点为xk-1,xk,xk+1,要求二次插值多项式L2（X），使它满足 L2（Xj）=Yj (j=k-1,k,k+1)我们知道y= L2（X）在几何上就是通过三点（xk-1,yk-1）, （xk,yk）, （xk+1,yk+1）的抛物线,为了求L2（X）的表达式，可采用基函数方法。此时基函数lk-1(x)， lk(x)及 lk+1(x)是二次函数，且在节点上满足条件 lk-1(xk-1)=1, lk-1(xj)=0 (j=k,k+1)lk(xk)=1, lk(xj)=0 (j=k-1,k+1)lk+1(xk+1)=1, lk+1(xj)=0 (j=k-1,k)满足上述条件的插值基函数可通过以下方式求。例如求lk-1(x)，因它有两个零点xk及xk+1,故可表为： lk-1(x)=A*(X-XK)（X-XK+1）其中A为待定系数，可由条件lk-1(xk-1)=1定出： A=1/（（XK-1-XK)（XK-1-XK+1））于是lk-1(x)=(X-XK)（X-XK+1）/（（XK-1-XK)（XK-1-XK+1））同理可得 lk(x)=(X-XK-1)（X-XK+1）/（（XK-XK-1)（XK-XK+1）） lk+1(x)=(X-XK-1)（X-XK）/（（XK+1-XK-1)（XK+1-XK））利用二次插值基函数lk-1(x)、 lk(x)和lk+1(x),可得二次插值多项式 L2(x)=yk-1lk-1(x)+yklk(x)+yk+1lk+1(x)显然，他满足条件L2(xj)=yj(j=k-1,k,k+1).将插值基函数代入可得二元插值多项式。 6.3 拉格朗日一元全节点插值多项式 将上述方法推广到一般情况，若一个n次插值