matlab实验报告_2完整版.docVIP

MATLAB实验报告 1、在区间[-1,1]上分别取n=10、20用两组选中节点对龙格函数作多项式插值及三次样条插值，对每个n值，分别画出插值函数及f(x)的图形. 解： n=10 在matlab命令窗口中键入： x=-1:0.2:1; y=1./(1+25*x.^2); y1=interp1(x,y,pchip); y2=interp1(x,y,spline); plot(x,y1,o,x,y1,-,x,y2,*,x,y2,-.),grid legend(样本点,三次插值,三次样条插值) 回车得出： n=20 在matlab命令窗口中键入： x=-1:0.1:1; y=1./(1+25*x.^2); y1=interp1(x,y,pchip); y2=interp1(x,y,spline); plot(x,y1,o,x,y1,-,x,y2,*,x,y2,-.),grid legend(样本点,三次插值,三次样条插值)legend(样本点,三次插值,三次样条插值) 回车得出： 由结果可见，用两种方法画出的曲线在样本点之间取值并无太大差异，曲线亦基本上一致。 2、对于给定函数在区间[-1,1]上取,试求3次曲线拟合，试画出拟合曲线并打印出方程。 解：在matlab命令窗口中键入： x=-1:0.2:1;y=1./(25*x.^2+1); p=polyfit(x,y,3) 回车得出