牛顿插值法实验报告.pdf

姜改革太原科技大学数值分析实验报告

牛顿插值法
一、实验目的：学会牛顿插值法，并应用算法于实际问题。
二、实验内容：给定函数 f (x姜改革太原科技大学数值分析实验报告

牛顿插值法
一、实验目的：学会牛顿插值法，并应用算法于实际问题。
二、实验内容：给定函数 f (x)  x ，已知：
f ()  f ()  f () 
f ()  f () 
三、实验要求：
（1）用牛顿插值法求 4 次 Newton 插值多项式在 处的值，以此作为函数
的近似值  N () 。在 MATLAB 中用内部函数 ezplot 绘制出 4 次 Newton
插值多项式的函数图形。
（2）在 MATLAB 中用内部函数 ezplot 可直接绘制出以上函数的图形， 并与作
出的 4 次 Newton 插值多项式的图形进行比较。
四、实验过程：
1、编写主函数。打开 Editor 编辑器，输入 Newton 插值法主程序语句：
function [y,L]=newdscg(X,Y,x)
n=length(X); z=x; A=zeros(n,n);A(:,1)=Y';s=; p=;
for j=2:n
for i=j:n
A(i,j)=(A(i,j-1)- A(i-1,j-1))/(X(i)-X(i-j+1));
end
end
C=A(n,n);
for k=(n-1):-1:1
C=conv(C,poly(X(k)));
d=length(C);C(d)=C(d)+A(k,k);
end
y(k)= polyval(C, z);
L(k,:)=poly2sym(C);
1姜改革太原科技大学数值分析实验报告

%%%%%%%%%%%%%%%%%%
t=[2,,,,];
fx=sqrt(t);
wucha=fx-Y;
以文件名 保存。
2、运行程序。
（1）在 MATLAB 命令窗口输入：
>> X=[2,,,,]; Y
=[,,,,];
x=;[y,P]=newdscg(X,Y,x