差商及插值多项式.ppt

  1. 1、本文档共20页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多

**********关于差商及插值多项式一般地,称k-1阶差商的一阶差商为f(x)关于点的k阶差商。例如,已知f(x)在的函数值为:可以求得第2页,共20页,2024年2月25日,星期天2.差商的性质性质1:k阶差商是由函数值的线性组合而成,即其中以k=2进行证明。由得到第3页,共20页,2024年2月25日,星期天由得到从而第4页,共20页,2024年2月25日,星期天由性质1立刻可得。性质2:差商具有对称性,即k阶差商f[x0,x1,…,xk-1,xk]中,任意调换xi,xj的次序,其值不变。再由数学归纳法可证得:第5页,共20页,2024年2月25日,星期天性质3:若f(x)为n次多项式,则f[x,x0]为关于x的n-1次多项式。证明:已知故类似的可以得到:也就是说,对多项式求一次差商,次数降低一次。由于是的根,所以第6页,共20页,2024年2月25日,星期天3.差商的计算为构造Newton插值多项式方便起见,计算差商时,采用列表的方式进行。第7页,共20页,2024年2月25日,星期天例2.2已知函数y=f(x)的如下离散数据(1,0)、(2,2)、(4,12)、(5,20)、(6,70),试求其各阶差商.解:列差商表计算xy一阶差商二阶差商三阶差商四阶差商1022412520670258501121051第8页,共20页,2024年2月25日,星期天二、Newton插值多项式对于区间[a,b]内的离散点及相应的函数值,计算如下差商:可以求得:第9页,共20页,2024年2月25日,星期天第10页,共20页,2024年2月25日,星期天依次类推得到:令:则可以将函数f(x)表示成:第11页,共20页,2024年2月25日,星期天由如上构造,容易验证因此Nn(x)满足插值条件,是一个n次插值多项式。并称为n次Newton插值多项式。如果f(x)≈Nn(x),则误差为:☆第12页,共20页,2024年2月25日,星期天关于Newton插值多项式,有以下几个特点:1Newton插值多项式与同次Lagrange插值多项式相同,因而误差相同因为Newton插值多项式与Lagrange插值多项式满足相同的插值条件,由插值多项式的存在唯一性知因此,Newton插值多项式与Lagrange插值多项式的误差相同。这样,由Nn(x)=Ln(x)得到第13页,共20页,2024年2月25日,星期天这个表达式给出了n+1阶差商与n+1阶导数之间的关系式。例3已知,试求其如下差商解:由差商与导数的关系式得到练习:提示:第14页,共20页,2024年2月25日,星期天2.Newton插值多项式具有递推式由可知所以,具有递推公式:第15页,共20页,2024年2月25日,星期天所以,具有递推公式:由此可知:当求出n次插值多项式后,再增加一个节点时,只需要增加一项的计算即可。由Newton插值多项式的结构可以看出,在构造Newton插值多项式时,必须首先计算各阶差商。3.Newton插值多项式的计算第16页,共20页,2024年2月25日,星期天例4已知f(x)的五组数据(1,0)、(2,2)、(3,12)、(4,42)、(5,116),求N4(x)。如果再增加一个节点(6,282),求出N5(x),并计算N4(1.5)、N5(1.5).解:先由前五组数据列差商表10223124425116210307441022240.5得到:如果,再增加一点(6,282),就在上表中增加一行计算差商。628216646810.1第17页,共20页,2024年2月25日,星期天由Newton公式的递推式得到:得到:练习题:已知离散数据(1,0)、(2,2)、(4,12)、(5,20)求三次Newton插值多项式,增加一点(6,70)后,再求出四次Newton插值多项式。第18页,共20页,2024年2月25日,星期天本节(§3)要点1.掌握差商及其性质,导数与差商的关系2.掌

您可能关注的文档

文档评论(0)

xiaolan118 + 关注
实名认证
内容提供者

你好,我好,大家好!

版权声明书
用户编号:7140162041000002

1亿VIP精品文档

相关文档