文档介绍：-
. z.
大学数学与计算机科学学院
计算机上机实验报告
专业和班级
数学实验班
凌白顺
成绩
**
031202121
课程名称
数值计算方法
实验名-
. z.
大学数学与计算机科学学院
计算机上机实验报告
专业和班级
数学实验班
凌白顺
成绩
**
031202121
课程名称
数值计算方法
实验名称
插值算法的程序设计
实
验
目
的
和
要
求
实验目的
本实验主要涉及函数逼近问题中的插值与拟合问题。通过本实验可以使同学们更深入地理解这些数值方法的根本原理，并通过程序设计使同学们初步领会到数值计算方法与计算机程序设计的密切相关；
实验根本要求
将调试好的插值法程序解决一两组实验数据。
Matlab现成函数的使用。
分析数值结果
实
验
容
和
步
骤
实验的主要容
、用给出的拉格朗日插值算法程序解决一到两组函数逼近的问题。并给出问题和结果，用图形表示〔可以截图说明〕。
答：拉格朗日函数可编成M文件，
解决问题1：给定函数，画出函数图像；然后给定*=144,169,225；Y=12,13,15；运用拉格朗日函数得到插值函数，接着画出差值函数图像即可，并用拉格朗日法计算在*=175处的值与真实值进展比拟
-
. z.
L〔175〕=；而真实值为f〔175〕=；
解决问题2：给出函数y=ln(*);给出*=10,11,12,13;对应的Y=,,,;用拉格朗日插值计算ln()；并与真实值进展比拟：首先用已有的M文件得到插值函数并画出图像和真实图像进展比拟：
-
. z.
并且可以计算出L()=，真实值f()=;
设计牛顿插值算法的程序。将其与拉格朗日程序结果进展比拟得出结论。
设计Newton算法如下：function y=Newton(a,b,*)
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
-
. z.
y=b(1);
for k=2:n
T=1;
for r=1:k-1
T=T.*(*-a(r));
end
N=T*A(k,k);
y=y+N;
end
并用Newton法来解决前面的二道题：
并且计算的N〔175〕=；而真实值为f〔175〕=；
-
. z.
并且可以计算出N()=，真实值f()=;