数值计算(分析)实验报告1.docVIP

南昌航空大学 数学与信息科学学院 实 验 报 告 课程名称： 实验名称： 实验类型： 验证性■综合性□设计性□ 实验室名称： 班级学号： 姓名： 任课教师（教师签名）： 成绩： 实验日期： 2012/3/1=，i=0,1,2，…，n，称为插值基函数，且满足条件: ,(i,j=0,1,…,n) 近似多项式的余项（误差）为： 其中。 2、编写一个用于牛顿均差插值公式计算函数值的程序，要求先输出均差表，再计算x点的函数值，并应用于下面的数据： x 20 21 22 23 24 y 1.30103 1.32222 1.34242 1.36173 1.38021 求x=21.4时的三次多项式的值。 拉格朗日插值多项结构紧凑，理论分析方便。但当节点有所增减时，整个公式都将改变，因此使用不甚方便。为了克服这一缺点，引入如下便于计算的形式： 其中，…为待定系数，由插值条件确定。 当插值节点等距分布时： （1）令=+kh（h=0,1，…，n），x=+th，0≤t≤1，则得到牛顿前插公式： 其中向前差分 插值余项为： （2）令=+kh（h=0,1，…，n），x=+th，-1≤t≤0，则得到牛顿后插值公式： 其中后向差分 插值余项为： 当插值节点非等距分布时，得到牛顿均差插值多项式： 其中均差 插值余项为： 程序代码： 题目一（拉格朗日插值）代码如下: function lag(x,y,xi) dxi=xi-x; n=length(x); L=zeros(size(y)); L(1)=prod(dxi(2:n))/prod(x(1)-x(2:n)); L(n)=prod(dxi(1:n-1))/prod(x(n)-x(1:n-1)); for j=2:n-1 num=prod(dxi(1:j-1))*prod(dxi(j+1:n)); den=prod(x(j)-x(1:j-1))*prod(x(j)-x(j+1:n)); L(j)=num/den; end yi=sum(y.*L); fprintf(x=%f,y=%f,xi,yi); 题目二（牛顿均差插值）代码如下: function [f,P]=junchachazhi(x,y,z) / x、y分别为已知数据和对应函数值，z为待求数据，f为均差表，P为函数值/ n=size(x,2); f=zeros(n-1,n-1); k=1; for i=1:n-1 /得出均差表 / for j=k:n-1 if i==1 f(j,i)=(y(j+1)-y(j))/(x(j+1)-x(j)); else f(j,i)=(f(j,i-1)-f(j-1,i-1))/(x(j+1)-x(j-k+1)); end end k=k+1; end P=y(1); /以下代码为计算多项式函数P(x)的值/ k=1; for i=1:3 /插值次数 此处为3次插值/ d=f(i,i); for j=1:k d=d*(z-x(j)); end k=k+1; P=P+d; end 实验过程中需要记录的实验数据表格 题目一 计算结果 （1） x=[0.30 0.42 0.50 0.58 0.66 0.72]; y=[1.04403 1.08462 1.11803 1.15603 1.19817 1.23223]; lag(x,y,0.46) x=0.460000,y=1.100724 （2） x=[0.30 0.42 0.50 0.58 0.66 0.72]; y=[1.04403 1.08462 1.11803 1.15603 1.19817 1.23223]; lag(x,y,0.55) x=0.550000,y=1.141271 （3） x=[0.30 0.42 0.50 0.58 0.66 0.72]; y=[1.04403 1.08462 1.11803 1.15603 1.19817 1.23223]; lag(x,y,0.60) x=0.600000,y=1.166194 题目二 计算结果 x=[20 21 22 23 24]; y=[1.30103 1.32222 1.34242