计算方法-实验一牛顿插值法.doc

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

计算方法课程设计报告 实验一 牛顿K次插值多项式 姓名:黄仁化 学号:031010151551017 班级:计算机科学与技术2004班 日期 :二○○六年六月十日 一、实验目的: 1、掌握牛顿插值法的基本思路和步骤。 培养编程与上机调试能力。 二、牛顿插值法基本思路与计算步骤: 给定插值点序列(。构造牛顿插值多项式。输入要计算的函数点并计算的值,利用牛顿插值公式,当增加一个节点时,只需在后面多计算一项,而前面的计算仍有用;另一方面的各项系数恰好又是各阶均差,而各阶均差可用均差公式来计算。 牛顿插值法计算步骤: 输入值及(;要计算的函数点。 对给定的由 计算的值。 3.输出。 三:程序流程图: 四:程序清单: function[c, d]=newpoly(x, y) %牛顿插值的MATLAB实现 %这里 x为n个节点的横坐标所组成的向量,y为纵坐标所组成的向量。 %c为所求的牛顿插值多项式的系数构成的向量。 n=length(x);%取x的个数。 d=zeros(n, n); d(: , 1)=y; for j=2 : n for k=j : n d(k, j)=(d(k, j-1) - d(k-1, j-1)) / (x(k)-x(k-j+1)); end end c =d(n, n); for k=(n-1) : - 1 : 1 c =conv(c, poly(x(k))); m=length(c); c(m)=c(m)+d(k, k); end 五、测试数据与结果: 测试数据:(第三章习题第三题第2题) f(x)=lnx的数值如表所示, 构造牛顿插值X 0.4 0.5 0.6 0.7 0.8 lnx -0.916291 -0.693147 -0.510826 -0.357765 -0.223144 解Y0=-0.916291, y1=-0.693147, y2=-0.510826, y3=-0.357765, y4=-0.223144 newpoly([0.4,0.5,0.6,0.7,0.8],[ -0.916291, -0.693147, -0.510826, -0.357765, -0.223144]) 计算结果如下: ans = -0.3096 2.6083 -5.4861 5.6921 -2.4744 由此看出所求的牛顿多项式为: P(x)= -0.3096x4+2.6083x3-5.4861x2+5.6921x-2.4744 P(0.53)= -0.6347。 六、小结: 本实验通过MATLAB编程实现求解牛顿K次插值多项式,能加深对牛顿插值法的基本思路和步骤的理解。 1 姓名:黄仁化          学号:031010151551017        实验一 牛顿K次插值多项式

文档评论(0)

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

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档