数学模型课程设计二.docVIP

课程设计名称： 设计二：数据拟合 指导教师： 课程设计时数： 课程设计设备：安装了Matlab、C++软件的计算机 课程设计日期： 实验地点： 课程设计目的： 1. 了解最小二乘拟合的原理，掌握用MATLAB作最小二乘拟合的方法； 2. 学会利用曲线拟合的方法建立数学模型。 课程设计准备： 在开始本实验之前，请回顾相关内容； 需要一台准备安装Windows XP Professional操作系统和装有数学软件的计算机。 课程设计内容及要求 要求：设计过程必须包括问题的简要叙述、问题分析、实验程序及注释、实验数据及结果分析和实验结论几个主要部分。 1. 用切削机床进行金属品加工时，为了适当地调整机床，需要测定刀具的磨损速度，在一定的时间测量刀具的厚度，得数据如表所示，请选用合适的函数来描述切削时间与刀具厚度的关系。 切削时间t/h 0 1 2 3 4 5 6 7 刀具厚度y/cm 30.0 29.1 29.8 28.1 28.0 27.7 27.5 27.2 切削时间t/h 8 9 10 11 12 13 14 15 刀具厚度y/cm 27.0 26.8 26.5 26.3 26.1 25.7 25.3 24.8 在M文件下输入如下代码： t=0:15;%切削时间 y=[30.0,29.1,29.8,28.1,28.0,27.7,27.5,27.2,27.0,26.8,26.5,26.3,26.1,25.7,25.3,24.8];%输入需要拟合的数据 A=polyfit(t,y,1)%用一次多项式拟合数据 z=polyval(A,t);%根据拟合的曲线估计A的值 plot(t,y,k+,t,z,r)%执行绘图命令 hold on B=polyfit(t,y,2)%用二次多项式拟合数据 z=polyval(B,t);%根据拟合的曲线估计Ｂ的值 plot(t,y,k+,t,z,b) hold on C=polyfit(t,y,3)%用三次多项式拟合数据 z=polyval(C,t);%根据拟合的曲线估计Ｃ的值 plot(t,y,k+,t,z,g) hold on D=polyfit(t,y,4)%用四次多项式拟合数据 z=polyval(D,t);%根据拟合的曲线估计Ｄ的值 plot(t,y,k+,t,z,y)%执行绘图命令，在同一坐标系中比较拟合程度 得出一到四次多项式的参数分别为： A = -0.3099 29.5676 B = 0.0060 -0.3994 29.7766 C = -0.0021 0.0524 -0.6691 30.0583 D = -0.0002 0.0031 0.0035 -0.5178 29.9777 一次多项式为：y = -0.3099*t + 29.5676 二次多项式为：y = 0.0060*t^2 - 0.3994*t + 29.7766 三次多项式为：y = -0.0021*t^3 + 0.0524*t^2 - 0.6691*t + 30.0583 四次多项式为：y = -0.0002*t^4 + 0.0031*t^3 + 0.0035*t^2 - 0.5178*t + 29.9777 从图像上看三次多项式能尽可能多的将给出的离散点拟合在曲线上，所以y = -0.0021*t^3 + 0.0524*t^2 - 0.6691*t + 30.0583能合适的描述切削时间与刀具厚度的关系。 2. Malthus人口指数增长模型 年 1790 1800 1810 1820 1830 1840 1850 1860 人口（亿） 3.9 5.3 7.2 9.6 12.9 17.1 23.2 31.4 年 1870 1880 1890 1900 1910 1920 1930 1940 人口（亿） 38.6 50.2 62.9 76.0 92.0 106.5 123.2 131.7 年 1950 1960 1970 1980 1990 2000 人口（亿） 150.7 179.3 204.0 226.5 251.4 281.4 用以上数据检验马尔萨斯人口指数增长模型，根据检验结果进一步讨论马尔萨斯人口模型的改进。 在M文件下输入如下代码： %马尔萨斯人口增长模型：N(t)=N0*e^r(t-t0) t= 1790:10:2000; %输入时间 y=[3.9,5.