多项式拟合与插值的计算方法.docVIP

安徽财经大学统计与数学模型分析实验中心 《数学软件》实验报告 实验名称：曲线拟合与插值 使用软件： Matlab 实 验 目 的 熟练掌握多项式拟合与插值的计算方法 实 验 内 容（具体题目及程序） 1.某种合金中的主要成分为A,B两种金属，经过试验发现：这两种金属成分之和x与合金的膨胀系数y有如下关系，建立描述这种关系的数学表达式. x 37 37.5 38 38.5 39 39.5 40 40.5 41 41.5 42 42.5 43 y 3.4 3 3 2.27 2.1 1.83 1.53 1.7 1.8 1.9 2.35 2.54 2.9 程序如下： x=37:0.5:43; y=[3.4,3,3,2.27,2.1,1.83,1.53,1.7,1.8,1.9,2.35,2.54,2.9]; plot(x,y,*); p=polyfit(x,y,2) y1=0.1660*x.^2-13.3866*x+271.6231 plot(x,y,*,x,y1,r*-)) R=sum((polyval(p,x)-y).^2) 2.对 y=cosx的数据进行插值，比较各种插值方法 程序如下： x=[-2*pi:0.5*pi:2*pi]; y=cos(x); x1=-2*pi:0.3*pi:2*pi; y_nearest=interp1(x,y, x1,nearest); y_linear= interp1(x,y,x1); y_spline= interp1(x,y,x1, spline ); y_cubic= interp1(x,y,x1, cubic ); plot(x,y,o,x1,y_nearest,-,x1,y_linear, r* , x1,y_spline,k:,x1,y_cubic,k -); legend (original data,nearest,linear,spline,cubic) 3.在某处测得海洋不同深处水温 如下： 深度（m） 446 714 950 1422 1634 水温（°C） 7.04 4.28 3.40 2.54 2.13 试求在深度为500米，1000米，1500米处的水温。 h=[446,714,950,1422,1634]; w=[7.04,4.28,3.40,2.54,2.13]; hi=[500,1000,1500]; p=spline(h,w); ci=ppval(p,hi) plot(h,w, or) ,hold on,plot(hi,ci,*) legend (original data,spline) 实 验 结 果 分 析 1. 首先作出散点图: 发现： 有点像抛物线，故选二次函数拟合. p = 0.1660 -13.3866 271.6231 所以 将两个图形进行比较： 误差平方和： R=sum((polyval(p,x)-y).^2)= 0.2523 2. 3. 思考题： 何时应采取多项式拟合？ 通常，在解决实际问题时先将已知数据的散点图画出，然后设计拟合的曲线类型，最后根据某种准则选定最佳的曲线. 成绩 1