数学实验作业.doc

验证及改进马尔萨斯人口模型 （高迪，071111 摘要：人口规模的合理性不仅对未来地区经济和社会发展大有益处，而且对地区生态环境可持续发展也有重要意义。现有从1790年到1980年每隔十年的人口数目，欲结合这些数据利用数学软件Mathematica中的数据拟合及其他相关知识拟合Malthus（马尔萨斯）人口指数增长模型的参数，运用马尔萨斯（Malthus）人口模型，根据检验结果进一步讨论马尔萨斯人口模型的改进。 一、引言 人口数目每年甚至每时每刻都有增有减，然而仔细研究就会发现其实人口的增长是有一定的规律的，人口规模是城市规划和土地利用总体规划中一项重要的控制性指标。运用马尔萨斯人口模型，利用数学软件Mathematica可以进行合理分析并能做出较合理的预测，准确地预测未来人口的发展趋势，制定合理的人口规划和人口布局方案具有重大的理论意义和现实意义。 二、基础理论 此实验过程中用到的Mathematica命令有： 1）画图： ListPlot[data,PlotJoined-True];(数据连线) ListPlot[data,PlotStyle-PointSize[0.02]](0.02为点的大小)。 2）拟合函数： 注意，该拟合是非线性拟合，需要装入nonlineearfit.M软件包 使用NonlinearFit命令。 三、实验过程与结果 1、按照表中的数据，用软件画出数据连线图 选取1970年为初始年份，则相应的数据为：{0, 3.9}, {10, 5.3}, {20, 7.2}, {30, 9.6}, {40, 12.9}, {50, 17.1}, {60, 23.2}, {70, 31.4}, {80, 38.6}, {90, 50.2}, {100, 62.9}, {110, 76.0}, {120, 92.0}, {130, 106.5}, {140, 123.2}, {150, 131.7}, {160, 150.7}, {170, 179.3}, {180, 204.0}, {190, 226.5}。 用Mathmatica软件编程如下： a = {{0, 3.9}, {10, 5.3}, {20, 7.2}, {30, 9.6}, {40, 12.9}, {50, 17.1}, {60, 23.2}, {70, 31.4}, {80, 38.6}, {90, 50.2}, {100, 62.9}, {110, 76.0}, {120, 92.0}, {130, 106.5}, {140, 123.2}, {150, 131.7}, {160, 150.7}, {170, 179.3}, {180, 204.0}, {190, 226.5}}; P1=ListPlot[a, PlotJoined - True]; P2=ListPlot[a, PlotStyle - PointSize[0.02]]; Show[P1,P2]; 图形如下： 2、写出微分方程的解，带入初值，初值由表定 微分方程为 以1790年为初始年份，即t=0时，=3.9,将这组初始值代入上式得N(t)=3.9 3、用数据拟合模型中的参数k，写出人口增长方程 由于马尔萨斯人口模型的特点，因该选取前若干项，而不能选取全部年份数据，所以可以取{0, 3.9}, {10, 5.3}, {30, 9.6}, {50, 17.1},{70, 31.4}, {90, 50.2}这六组数据进行数据拟合。由于该拟合属于非线性拟合，需要用NonlinearFit命令。Mathmatica程序为： C:\\ProgramFiles\\WolframResearch\\Mathematica\\5.0\\AddOns\\StandardPackages\\Statistics\\NonlinearFit.m P3={{0,3.9},{10,5.3},{20,7.2},{30,9.6},{40,12.9},{50,17.1}}; NonlinearFit[P3,3.9*Exp[?*t],t,?] 所以=0.0297009，则 验证所得到的方程，分析是否符合现实数据 f[t_]=3.9*Exp[0.0297009*t] For[i=0,i20,i+=1,j=i*10; Print[f[,i*10,]=,a[[i+1,2]],...... f[,i*10,]=,f[j]];]//N f[ 0 ]= 3.9 ...... f[ 0 ]= 3.9 f[ 10 ]= 5.3 ...... f[ 10 ]=