数值分析15切比雪夫插值.ppt

《数值分析15切比雪夫插值.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、切比雪夫插值节点带导数条件的插值函数分段插值函数二元函数插值简介《数值分析》15取插值结点:a≤x0＜x1＜······＜xn≤b满足Ln(xk)=f(xk)的n次多项式插值余项其中,选取:x0,x1,······,xn,使结论:切比雪夫多项式Tn+1(x)的全部零点。拉格朗日插值余项2/18n+1阶切比雪夫多项式:Tn+1=cos(n+1)cos=x代入得Tn+1(x)=cos((n+1)arccosx)即(k=0,1,···,n)取f(x)∈C[–1,1],令x=cos,则有[–1,1][0,]将g()=f(c

2、os)展开成余弦级数——切比雪夫结点3/18例1.函数取等距插值结点:-5,-4,-3,-2,-1,0,1,2,3,4,5x∈[-5,5]11(x)=(x+5)(x+4)(x+3)(x+2)(x+1)x(x-1)(x-2)(x-3)(x-4)(x-5)11(x)4/18-4.9491-4.5482-3.7787-2.7032-1.40870.00001.40872.70323.77874.54824.9491在[-5,5]区间上,取11个切比雪夫结点(k=10,9,8,···,1,0)11(x)=(x–x0)(x–x1

3、)(x–x2)······(x–x10)5/1811(x)插值函数L10(x)取切比雪夫结点插值插值函数L10(x)取等距结点插值6/18已知节点x0和x1处的函数值及导数值求三次插值函数H(x)=a0+a1x+a2x2+a3x3满足插值条件(j=0,1)三次Hermite插值问题xx0x1H(x)y0y1H’(x)m0m17/18例2.已知插值条件:求3次插值函数.解:设得a0=0,a1=0,列出方程组求解,得a2=3,a3=–2所以,有H(x)=3x2–2x3=(3–2x)x2x01H(x)01H’(x)008/18利用

4、基函数表示Hermite插值x0x110000100xx0x100100001x9/18两点Hermite插值的误差估计式证明:由插值条件知R(x0)=R’(x0)=0,R(x1)=R’(x1)=0构造辅助函数利用f(x)–H(x)=C(x)(x–x0)2(x–x1)2取x异于x0和x1,设10/18反复应用Roll定理,得F(4)(t)有一个零点设为ξ显然,F(t)有三个零点x0,x,x1,由Roll定理知,存在F’(t)的两个零点t0,t1满足x0

5、异零点.11/18分段线性插值插值节点满足:x0

6、(x1,y1)=z1,P(x2,y1)=z2,P(x2,y2)=z3,P(x1,y2)=z4P(x,y)=ax+by+cxy+dl1(u，v)=(1–u)(1–v)l2(u，v)=u(1–v)l3(u，v)=uvl4(u，v)=(1–u)v其中P(x,y)=z1(1–u)(1–v)+z2u(1–v)+z3uv+z4(1–u)v16/18[u,v]=meshgrid(0:0.1:1)；L1=(1-u).*(1-v);surf(u,v,L1)figureL2=u.*(1-v);surf(u,v,L2)figureL3=u.*v;su

7、rf(u,v,Lu3)figureL4=(1-u).*v;surf(u,v,L4)x=asinφcosθy=asinφsinθz=bcosφ三角形区域线性插值插值条件:z1=P(x1,y1)z2=P(x2,y2)z3=P(x3,y3)(x1,y1)(x3,y3)(x2,y2)拉格朗日方法P(x,y)=l1(x,y)z1+l2(x,y)z2+l3(x,y)z3P(x,y)=ax+by+c(x，y)(x1，y1)(x2，y2)(x3，y3)l1(x,y)100l2(x,y)010l3(x,y)00117/18l1(x,y)的图形

8、是空间三角形分片线性插值18/18P(x,y)=l1(x,y)z1+l2(x,y)z2+l3(x,y)z3图形是空间三角形