文档介绍：《数值分析》实验报告班级：
姓名：
学号：
指导老师：
实验基本要求
一、上机前的准备工作1、复****和掌握与本次实验有关的教学内容。
2、根据本次实验要求，在纸上编写算法及上机的程序，并经过人工模拟运
行检验，减少不必要0675
算法2functiony2(i)formatlongy2=zeros(1,21);y2(21)=;fori=21:-1:2y2(i-1)=/(i-1)-*y2(i);endy2(1:21)ans=**********.，算法1不稳定。
实验二、插值法
一、实验目的1、理解插值的基本概念，掌握各种插值方法，包括拉格朗日插值和牛顿插值等，注意其不同特点；2、通过实验进一步理解并掌握各种插值的基本算法。
二、实验原理
插值法是函数逼近的一种重要方法，它是数值积分、微分方程数值解等数值计算的基础与工具，其中多项式插值是最常用和最基本的方法。拉格朗日插值多项式的优点是表达式简单明确，形式对称，便于记忆，它的缺点是如果想要增加插值节点，公式必须整个改变，这就增加了计算工作量。而牛顿插值多项式对此做了改进，当增加一个节点时只需在原牛顿插值多项式基础上增加一项，此时原有的项无需改变，从而达到节省计算次数、节约存储单元、应用较少节点达到应有精度的目的。
三、实验任务1、已知函数表
xi
yi
。
2、已知函数表
用牛顿插值多项式求N3()
和N4()。
四、实验源程序及结果1、二次拉格朗日插值function[y,R]=lagranzi(X,Y,x,M)formatlongx=;M=2;X=[,,,];Y=[,,,];n=length(X);m=length(x);fori=1:mz=x(i);s=;fork=1:np=;q1=;c1=;forj=1:nifj~=kp=p*(z-X(j))/(X(k)-X(j));endq1=abs(q1*(z-X(j)));c1=c1*j;ends=p*Y(k)+s;endy(i)=s;endR=M.*q1./c1;ans=、牛顿插值多项式求N3()和N4()1）N3()：
function[y,R]=newcz(X,Y,x,M)M=3;X=[,,,,];Y=[,,,,];x1=;
n=length(X);m=length(x1);fort=1:mz=x1(t);A=zeros(n,n);A(:,1)=Y';s=;p=;q1=;c1=;forj=2:nfori=j:nA(i,j)=(A(i,j-1)-A(i-1,j-1))./(X(i)-X(i-j+1));endq1=abs(q1.*(z-X(j-1)));c1=c1.*j;endC=A(n,n);q1=abs(q1.*(z-X(n)));fork=(n-1):-1:1