多项式插值唯一性.ppt_第1页
多项式插值唯一性.ppt_第2页
多项式插值唯一性.ppt_第3页
多项式插值唯一性.ppt_第4页
多项式插值唯一性.ppt_第5页
已阅读5页,还剩13页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

多项式插值的存在唯一性拉格朗日插值,牛顿插值埃米特插值与三次样条数据拟合的线性模型两种典型的正交多项式,数值分析习题课III,若插值结点x0,x1,xn是(n+1)个互异点,则满足插值条件P(xk)=yk(k=0,1,n)的n次插值多项式P(x)=a0+a1x+anxn存在而且惟一。,多项式插值的存在唯一性定理,Laglarge插值公式,插值基,(k=0,1,2,n),2/18,插值误差余项,其中,线性插值误差:,二次插值误差:,思考:构造线性插值函数计算115的平方根近似值,估计近似值的误差并指出有效数位数。,3/18,已知x0,x1,xn处的值f(x0),f(x1),f(xn).,(j=0,1,n-1),(j=0,1,n-2),均差的定义,牛顿插值公式,(k=1,2,n),思考:证明一阶差商的对称性:fx0,x1=fx1,x0,进一步证明二阶差商的对称性。,4/18,牛顿插值余项,(j=0,1),三次Hermite插值,5/18,给定a,b的分划:a=x0x1xn=b.已知f(xj)=yj(j=0,1,n),如果,满足:(1)S(x)在xj,xj+1上为三次多项式;(2)S”(x)在区间a,b上连续;(3)S(xj)=yj(j=0,1,n).则称S(x)为三次样条插值函数.,三次样条的定义,6/18,(j=1,2,n-1),自然边界条件,三次样条一阶导数值:S(xj)=mj(j=0,1,n),三次样条二阶导数值:S”(xj)=Mj(j=0,1,n),j=1,n1,自然边界条件:M0=0,Mn=0,7/18,拟合函数:(x)=a00(x)+a11(x)+ann(x),数据拟合的线性模型,离散数据,超定方程组,超定方程组最小二乘解:,8/18,对连续函数f(x)的正交多项式平方逼近,其中,Ex1.设x0,x1,xn是互异的插值结点,l0(x)为对应于x0的拉格朗日插值基函数,试证明,9/18,Ex2.设x0,x1,x2,xn为互异的结点,求证Lagrange插值基函数满足下列恒等式,(1),(2),(k=1,n),证:(1)令,在插值结点处Pn(xj)=0(j=0,1,2,n),n次多项式Pn(x)有n+1个相异零点Pn(x)=0,10/18,所以,将f(x)=xk(kn)代入,得,(k=0,1,2,n),思考题:f(x)是(n+1)次多项式且最高次项系数为1,取互异的插值结点x0,x1,xn,构造插值多项式Pn(x),证明:f(x)=Pn(x)+(xx0)(xx1)(xxn),(2)取f(x)=xkf(n+1)(x)=0Rn(x)=0,11/18,Ex4.设x0x1x2,从函数表,出发,利用f(x)的二次拉格朗日插值多项式L2(x)推导出求f(x)的极值点x*的近似值计算公式.,Ex3.设P(x)是不超过n次的多项式,而n+1(x)=(xx0)(xx1)(xxn)证明存在常数Ak(k=0,1,n)使得,12/18,Ex5.设有数列:x1,x2,xn,(1).证明平方和数列为3阶等差数列,证明:(1)Sn=n2,2Sn=n2(n-1)2=2n-13Sn=(2n-1)-(2n-3)=2,故平方和数列为3阶等差数列.,(2).证明,则,(2)令g(n)=n(n+1)(2n+1)/6,13/18,同理,(k=1,2,n),显然,14/18,证明:Fx0,x1,xn=,Ex6.记n+1(x)=(xx0)(xx1)(xxn),(j=1,2,n),对比Lagrange插值和Newton插值中xn的系数,得,Fx0,x1,xn=,15/18,Ex7.2次埃尔米特插值的适定性问题,给定插值条件:f(x0)=y0,f(x1)=m1,f(x2)=y2,插值结点应满足什么条件能使插值问题有唯一解。,思考:构造带导数条件的二次插值多项式公式f(0)=y0,f(1)=y1,f(0)=m0;,16/18,解:设H(x)=a0+a1x+a2x2,H(x)=a1+2a2x,Ex8.如果xa,b,t-1,1,(1)证明联系两个区间的映射为,17/18,Ex9.一个量x被测量了n次,其结果是a1,a2,an.用最小二乘法解超定方程组x=a

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论