还原图1.4(绘制Lagrange 10次插值曲线和原函数曲线),并加上图1.5中的分段线性差值曲线。(三条曲线绘制在同一张Figure里面,注意用不同的线形区分)
第一步 将function yy=Lagrange(x,y,xi)m=length(x);n=length(y);if m ~= n , error('向量x与y的长度必须一致'); end; s=0;for i=1 : n z=ones(1,length(xi)); for j=1 : n if j ~= i z=z.*(xi-x(j))./(x(i)-x(j)); end end s=s+z*y(i);endyy=s;保存为M文件。(文件→新建→M文件)第二步将x=[0.5610,0.56280,0.56401,0.56521];y=[0.82741,0.82659,0.82577,0.82495];xi=[0.5625,0.5635,0.5645];yi=Lagrange(x,y,xi)粘贴至(命令窗口)。