福州大学数值计算实验报告一.docVIP

福州大学数学与计算机科学学院 计算机上机实验报告 专业和班级 数学实验班 姓名 凌白顺 成绩 学号 031201321 课程名称 数值计算方法 实验名称 插值算法的程序设计 实 验 目 的 和 要 求 实验目的 本实验主要涉及函数逼近问题中的插值与拟合问题。通过本实验可以使同学们更深入地理解这些数值方法的基本原理，并通过程序设计使同学们初步领会到数值计算方法与计算机程序设计的密切相关； 实验基本要求 将调试好的插值法程序解决一两组实验数据。 Matlab现成函数的使用。 分析数值结果 实 验 内 容 和 步 骤 实验的主要内容 、用给出的拉格朗日插值算法程序解决一到两组函数逼近的问题。并给出问题和结果，用图形表示（可以截图说明）。 答：已知拉格朗日函数可编成M文件， 解决问题1：给定函数，画出函数图像；然后给定X=144,169,225；Y=12,13,15；运用拉格朗日函数得到插值函数，接着画出差值函数图像即可，并用拉格朗日法计算在x=175处的值与真实值进行比较 L（175）=13.2302；而真实值为f（175）=13.2288； 解决问题2：给出函数y=ln(x);给出X=10,11,12,13;对应的Y=2.3026,2.3979,2.4849,2.5649;用拉格朗日插值计算ln(11.75)；并与真实值进行比较：首先用已有的M文件得到插值函数并画出图像和真实图像进行比较： 并且可以计算出L(11.75)=2.4639，真实值f(11.75)=2.4639; 设计牛顿插值算法的程序。将其与拉格朗日程序结果进行比较得出结论。 设计Newton算法如下：function y=Newton(a,b,x) n=length(a); A=zeros(n); A(:,1)=b; for k=2:n for r=k:n A(r,k)=(A(r,k-1)-A(r-1,k-1))/(a(r)-a(r-k+1)); end end y=b(1); for k=2:n T=1; for r=1:k-1 T=T.*(x-a(r)); end N=T*A(k,k); y=y+N; end 并用Newton法来解决前面的二道题： 并且计算的N（175）=13.2302；而真实值为f（175）=13.2288； 并且可以计算出N(11.75)=2.4639，真实值f(11.75)=2.4639; 利用三次样条插值函数spline()解决一两组问题。 给定数据 x = 0:10; y = sin(x); xx = 0:.25:10; yy = spline(x,y,xx); plot(x,y,o,xx,yy) 研 究 与 探 讨 通过用拉格朗日插值、牛顿插值以及三次样条插值来处理数据，并且对数据进行分析，然后可以估算出一些未知数据。通过对三种插值方法进行编程以及解决一些实际问题时发现各种方法都各有利弊，首先最基本的是拉格朗日插值方法，对于拉格朗日插值方法来说编程是非常简单的，但是一旦增加一组数据时，程序就要相应做出改变，而对于现实中的问题来说，增加或减少数据这是常见的事，这使得拉格朗日的实用性不强，而对于牛顿法来说，编程相对与拉格朗日方法来说更加难一点，插值效果上面通过前面的实验可以看出，牛顿法的插值效果和拉格朗日插值的效果是一样的，但是牛顿法可以有效的改善拉格朗日法弊端，如果增加或者减少数据时，程序不用发生改变，这使得牛顿法在日常使用中更加具有实用性。但是对于拉格朗日插值和牛顿插值来说，当插值的图像是一段段的折线形成的，其光滑性不好，精确度也不是很好，相对于拉格朗日法个牛顿法，三次样条插值是一种更加精确的插值方法，图像也更加光滑，同时用法也比较简单，三次样条插值也是数学实验中常常使用的一种方法。