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

-

大学数学与计算机科学学院

计算机上机实验报告

凌白顺

专业和班级数学实验班成绩

**031202121

课程名称数值计算方法实验名称插值算法的程序设计

实验目的

实本实验主要涉及函数逼近问题中的插值与拟合问题。通过本

验实验可以使同学们更深入地理解这些数值方法的根本原理，并通

目过程序设计使同学们初步领会到数值计算方法与计算机程序设计

的的密切相关；

和实验根本要求

要1.将调试好的插值法程序解决一两组实验数据。

求2.Matlab现成函数的使用。

3.分析数值结果

实验的主要容

实

（1）、用给出的拉格朗日插值算法程序解决一到两组函数逼近的

验

问题。并给出问题和结果，用图形表示〔可以截图说明〕。

答：拉格朗日函数可编成M文件，

容

解决问题1：给定函数yx，画出函数图像；然后给定*=144,169,225；

和

Y=12,13,15；运用拉格朗日函数得到插值函数，接着画出差值函数图像即可，

步

并用拉格朗日法计算在*=175处的值与真实值进展比拟

骤

.z.

-

L〔175〕=13.2302；而真实值为f〔175〕=13.2288；

解决问题2：给出函数y=ln(*);给出*=10,11,12,13;对应的

Y=2.3026,2.3979,2.4849,2.5649;用拉格朗日插值计算ln(11.75)；并与真实

值进展比拟：首先用已有的M文件得到插值函数并画出图像和真实图像进

展比拟：

.z.

-

并且可以计算出L(11.75)=2.4639，真实值f(11.75)=2.4639;

(2)设计牛顿插值算法的程序。将其与拉格朗日程序结果进展比拟

得出结论。

设计Newton算法如下：functiony=Newton(a,b,*)

n=length(a);

A=zeros(n);

A(:,1)=b;

fork=2:n

forr=k:n

A(r,k)=(A(r,k-1)-A(r-1,k-1))/(a(r)-a(r-k+1));

end

end

.z.

-

y=b(1);

fork=2:n

T=1;

forr=1:k-1

T=T.*(*-a(r));

end

N=T*A(k,k);

y=y+N;

end

并用Newton法来解决前面的二道题：

并且计算的N〔175〕=13.2302；而真实值为f〔175〕=13.2288；