数学实验实验报告二答案.doc

文档介绍：实验二实验项目名称:插值与拟合实验内容:1、。Newscript源程序:functiony=lagr(x0,y0,x)n=length(x0);m=length(x);fori=1:mz=x(i);s=0;fork=1:np=1;forj=1:nifj~=kp=p*(z-x0(j))/(x0(k)-x0(j));endends=p*y0(k)+s;endy(i)=s;、完成机床加工问题中数据增加密度,利用三种插值方法计算加工所需数据,并画出曲线,比较分析三种插值结果。(参考教材P77)源程序:x0=[0,3,5,7,9,11,12,13,14,15];y0=[0,,,,,,,,,];x=0::15;y1=lagr(x0,y0,x);y2=interp1(x0,y0,x);y3=interp1(x0,y0,x,'spline');[x',y1',y2',y3'];subplot(3,1,1)plot(x0,y0,'k+',x,y1,'r')grid,title('lagrange')subplot(3,1,2)plot(x0,y0,'k+',x,y2,'r')grid,title('piecewiselinear')subplot(3,1,3)plot(x0,y0,'k+',x,y3,'r')grid,title('spline')实验结果:3、完成教材P80习题1中的函数y=sinx的三种插值结果的分析和比较。源程序:n=10;m=100;x=linspace(0,2*pi,100);y=sin(x);z=x*0;x0=linspace(0,2*pi,6);y0=sin(x0);y1=lagr(x0,y0,x);y2=interp1(x0,y0,x);y3=interp1(x0,y0,x,'spline');plot(x,z,x,y,'k',x,y1,x,y2,x,y3,'y')实验结果:4、完成教材P81习题2。源程序:x0=[0,1,4,9,16];y0=sqrt(x0);x=5;y=s