数值计算实验报告.docVIP

北方名族大学 信息与计算科学学院 数值分析课程设计 指导老师：杨录峰 学生姓名：李俊 2010级信计2班 学号 2012年6月 牛顿插值法计算函数值 算法分析 本题在算法上需要解决的问题主要是：求出第二问中的Newton插值多项式。如此，则第三、四问则迎刃而解。计算牛顿插值多项式的算法如下： 求Newton插值多项式，算法组织如下： Newton插值多项式的表达式如下： 其中每一项的系数ci的表达式如下： 根据以上公式，计算的步骤如下： 算法实现 输入数据 计算参数 直接采用牛顿插值公式计算 4.利用表达式进行插值 利用插值余项公式可知 说明每增加一个结点，Newton插值多项式只增加 一项，克服了 Lagrange插值的缺点。实际计算时可 利用差商表。 下面建立差值表 一阶插商 二阶插商 三阶插商 四阶插商 四 源代码 %牛顿插值法 Function yy=newton(x,y,z) x=a:(b-a)/n:b; %插值节点 y=f(x); plot(x,y,b) %用蓝色线作被插函数图象 hold on z=a:(b-a)/(2*n):b; n=length(x); for j=2:n for i=n:-1:j y(i)=(y(i)-y(i-1))/(x(i)-x(i-j+1));%计算差商 end end u=y(n); m=length(z); for j=1:m for i=n-1:-1:1 u=y(i)+u*(z(j)-x(i)); %计算牛顿插值多项式的值 v(j)=u; end u=y(n); end plot(z,v,r) %用红色线作牛顿插值多项式图象 3、不论n为多少，是不会改变的，其值存储在矩阵yy中； 为了使整个结果直观，实验的最终结果还用图形进行的重现，牛顿插值结果存在数组变量Nn。 图5- 与原始值的对比图 根据第3问中得到的数据可以很容易的得到它们的值如下表所示： n值 E(Nn) 5 0.43269 20 58.2781 当n=20时，使用Newton插值多项式出现龙格现象，其最大误差达到58.2781可见，n越大，Newton插值越可能偏离被插值函数. 六．结果分析 牛顿法具有非常强的承袭性，这就保证了前面计算的数据不会造成浪费，而且在后面的计算中，下一步的结果是综合了前面两个数据，因此结果更加准确，精度更加高，这个公式的计算也是比较容易的，特别是其系数可以依据牛顿插值中的定理4来快速计算。但是，它也有缺点，那就是随着数据的增多，以前微小的误差就累计了一起然后就成了巨大的误差，就出现了书中1.7的龙格现象。 七．心得体会 通过半学期的的数值计算学习，收获虽然不小，但是更多的是经验与教训，比如matlab的编程，总感觉心有余而力不足，最开始简单程序就编得很吃力，导致学习没有太大激情，下课又不强加练习，以至于后来面对复杂的程序时无法下手，从而慢慢觉得数值计算就是难，就是无法学！但是虽然课程不简单，但是老师的讲解却是面面俱到，每一种方法，每一个程序，老师都细细的评讲，耐心的询问我们有何不明白，然后还不厌其烦的为我们一一讲解。使得我对这课程慢慢拾起信心，慢慢接受它，慢慢学好它。所以，我非常感谢杨老师的谆谆教导。 1