人口增长模型的确定

1、题目：人口增长模型的确定摘要人口问题己成为当前世界上最普遍关注的问题之一，人口增长规律的发现以及人口增长的预测问题对一个国家制定长远的发展规划有着非常重要的意义。本文分别使用了马尔萨斯人口指数增长模型和阻滞增长模型，以美国1790-1980年间每隔10年的人口数量为依据，对接下來的每隔十年进行了预测五次人口数量。通过对比我们可以发现阻滞增长模型在预测准确度方面要明显优于原始的马尔萨斯人口指数增长模型。关键词：人口增长；马尔萨斯人口指数增长模型；阻滞增长模型；人口预测50*问丿重述1.1问题背景1790-1980年间美国每隔10年的人口记录如下表所示。表1人口记录表年份1790180018101

2、820183018401S50186018701880人II(xlO6)3.95.37.29.612.917.123.231.438.650.2年份1890190019101920193019101950196019701980人II(xlO6)62.976.092.0106.5123.2131.7150.7179.3204.0226.51.2问题提出我们需要解决以下问题：试用以上数据建立马尔萨斯(Malthus)人口指数增长模型，并对接下來的每隔十年预测五次人口数量，并查阅实际数据进行比对分析。如果数据不相符，再对以上模型进行改进，寻找更为合适的模型进行预测，并对两次预测结果进行对比分析。查

3、阅资料找出中国人口与表1同时期的人口数量，用以上建立的两个模型进行人口预测与分析。二、问题分析首先，我们运用Matlab软件绘制出1790到1980年的美国人口数据图，如图1。250200150100(7801800182018401860188019001920194019601980图11790到1980年的美国人口数据图从图表中我们可以清晰地看到人口数在1790-1980年是呈增长趋势的，而且我们很容易发现上述图表和我们学过指数函数的图表有很大的相似性，所以我们很自然想到建立指数模型。因此我们首先建立马尔萨斯模型，马尔萨斯生物总数增长定律指出：在孤立的生物群体中，生物总数N的变化率与生物

4、总数成正比。三、问题假设为简化问题，我们做出如下假设：在模型中预期的时间内，人口不会因发生大的自然灾害，突发事件或战争而受到大的影响；所给出的数据具有代表性，能够反映普遍情况；一段时间内我国人口死亡率不发生大的波动；在查阅的资料与文献中，所得数据可信；假设人口净增长率为常数。四、变量说明在此，对本文所使用的符号进行定义。表2变量说明符号符号说明N(0)起始年人口容纳量N(t)t年后人口容纳量t年份r增长率五、模型建立5.1问题一：马尔萨斯(Malthus)人口指数增长模型设：t表示年份(起始年份t=0),r表示人口增长率，N表示t年后的人口数量。当考察一个国家或一个很大地区的人口时，N(t)是

5、很大的整数。为了利用微积分这一数学工具，将N(t)视为连续、可微函数。记初始时刻(t=0)的人口为N(0),人口增长率为r,r是单位时间内N的增量与N(t)的比例系数。根据r是常数的基本假设，于是N满足如下的微分方程：dN(t)/dt=r*N(t)(5-1)由这个线性常系数微分方程容易解出：N(t)=N(0)ert(5-2)表明人口将按指数规律无限增长(r0)o将以t年为单位，上式表明，人口以e为公比的等比数列增长。因为这时1-表示年增长率，通常11,所以可用近似关系e“l+r可得出N(t)=N(O)(l十Qt(5-3)(5-3)式即人口增长模型。5.2问题二：改进模型-阻滞增长模型(Logi

6、stic模型)自从英国人口学家和政治经济学家托马斯罗伯特马尔萨斯1798年发表人口学原理后，马尔萨斯(Malthus)人口指数增长模型在世界上引起了轩然大波，并在后來的人口预测中扮演着重要的角色。但是随着时间的发展，由于现代社会与自然环境的改变，马尔萨斯人口指数增长模型在预测未來人口时，误差可能会比较大。上述模型对较早时期的统计数据吻合得较好，但也存在问题，即人口是呈指数规律无止境地增长，此时人口的自然增长率随人口的增长而增长，这不可能。由于社会的快速发展，自然环境遭受严重破坏，人口的高速增长等一系列原因，人口的增长率不能按照马尔萨斯所假设为一个常数I不改变。一般说來，当人口较少时增长得越來越

7、快，即增长率在变大；人口增长到一定数量以后，增长就会慢下來，即增长率变小这是因为，自然资源、环境条件等因素不允许人口无限制地增长，它们对人口的增长起着阻滞作用，而且随着人口的增加,阻滞作用越來越大。而且人口最终会饱和，趋于某一个常数Xoo,我们假设人口的净增长率为i(l-X(t)/Xoe),即人口的净增长率随着人口的增长而不断减小，当tTS时，净增长率趋于零。按照这个假设，得到：色=口1-型)山耳(5-4)?Vo)=X。这便是荷兰数学家Veihulst于19世纪中叶提出的阻滞增长模型(Logistic模型)。在MATLAB命令窗口键入dsolve(，Dx=r*x*(1-x/c)，,，x(179

8、0)=39)输出：ans=c/(l+l/39*exp(-r*t)*exp(1790*r)*(10*c-39)其中C=Xx因此，人口的变化规律为：一1+(甥心-I)。”S)5.3问三模型建立经调查，1790-2010年间中国每隔10年的人口记录如下表所示。表7中国人口记录表年份179018001S10182018301840185018601870人11(xlO6)323.5341.6360.8381.0409.0412.0412.0377.0358.0年份188018901900191019201930191019501960人1(x10)368.0380.0400.0423.0472.048

9、9.0518.8546.8667.1年份19701980199020002010人1(x10)818.3981.21135.21264.11311.0我们分别应用马尔萨斯人口指数增长模型和Logistic人口阻滞增长模型來对中国人口进行预测。六、模型求解6.1问题一模型求解在应用预测模型的过程中考虑到，若要提高预测结果的准确性，就必须增加预测方案的数量，对比各方案的预测值和误差，选取误差最低的一组预测方案。特别是马尔萨斯模型中，人口增长率T是一定时期内人口增加的综合结果，在预测中它的取值直接关系到预测结果的精度，因此在进行不同阶段的人口预测时根据实际情况对人口增长率r加以分类和处理才能得到理想

10、的预测结果。本文根据1790-1980年计算美国常住人口每年的增长率，按照人口增长率的大小设置了高中低三个方案，以此加强预测结果的对比，提高预测的准确度。表3美国每10年自然增长率序号123156789r0.3590.3580.3330.3430.3260.3570.3530.2290.301101112131415161718190.2530.2080.2110.1580.1570.0690.1440.1900.1380.110通过表3可以确定自然增长率高，中，低三个方案。通过数据分析可得，上述表格为10年的累计增长率，而自然增长率强调一年，所以可近似除以10求得，高方案中自然增长率为0.0

11、33,中方案中自然增长率为0.029,低方案中自然增长率为0.013o依据人口增长率的大小分为高、中、低、三个预测方案，将预测值与实际值进行拟合比较。图21=0.033时马尔萨斯模型曲线拟合图31=0.029马尔萨斯模型曲线拟合图41=0.013马尔萨斯模型曲线拟合根据上述分析，及曲线拟合可知，取中方案即1=0.029时，马尔萨斯模型更符合实际情况。因此本文自然增长率取1=0.029來预测美国人口数量并与实际情况对比。由预测公式预测1790-1980年的人口数量，由指数增长模型可得各个年份的真实值与预测值之间的差别如下表：表41790年-1980年美国人口真实值与预测值年份1790180018

12、101820183018401850186018701880真实人II(xlO6)3.95.37.29.612.917.123.231.43&650.2误差(xlO6)00.10.290.410.660.812.550.20.95.1年份1890190019101920193019401950196019701980真实人II(xlO6)62.976.092.0106.5123.2131.7150.7179.3204.0226.5误差(xlO6)14.528.450.3190.2152.3227.3323.3465.6664.7800.5通过调查得知1990-2010年人口数量统计如下表表51

13、990年-2010年美国人口真实值与预测值时间（年）199020002010人口248.7281.4308.7误差(xlO6)93912971792.4图5美国人口真实值与预测值曲线拟合通过上图可以发现，1790-1870指数增长模型确实拟合的比较好，但从1870年开始往后发现误差越來越大，可知指数增长模型只适合于短期的人口预测。为了生存以及人类的发展，人们自然会釆取有效措施來控制人口的过度增长，自然资源、环境资源的条件也限制了人口数量的过度增长。因此为了使人口预报模型适合长期的发展趋势，更好地符合实际情况，必须修改指数增长模型中关于人口增长率是常数这个基本假设了，这时必将导致更适合人类发展的

14、规律的新数学模型的产生。6.2问题二模型求解利用MATLAB软件中的“ciirvefit”命令和式(5-5)来拟合所给的人口统计数据，从而确定出(5-5)中的待定参数I和Xooo查阅资料可得1-的初值取为小于1的数，比如取a=200,0.1时，得到a=311.95570.0280,yl=267.1959,即(5-5)中的r=0.0280.x=311.9557,2010年美国的人口预计为267.1959百万人。这个结果还比较合理，当t趋于无穷时，静增长率趋于零，人口数趋F311.9557百7J人，即极限人口Xx=311.9557百万百万。拟合效果见图5。根据该题己给数据可作如下图形:300图6L

15、ogistic模型拟合曲线从图6可以看出，在前一段吻合得比较图，但在最上面，若拟合曲线更接近原始数据，对将來人口的预测应该更好。因此略加修改将拟合准则改为：TOC o 1-5 h z/I21ininE=工（/（匚）一兀才+w工（/匕）一兀）2（5-6）；=1i=/hi其中W为右端儿个点的误差权重，在此处应该取为大于1的数，这样会使右边的拟合误差减小，相应的，其他点的误差会有所增加。我们要使这些误差的增减恰当，可以通过调整歹和门的具体取值，比较他们取各种不同值时的拟合效果，从而确定出一个合适的数值。1）先取n二17,w二1.5,运行上述程序，得到结果a二324.0666,0.0276;xl=27

16、2.7996.2）再取n二16,w=2,运行上述程序，得到结果a二345.1439,0.0270;xl=280.0539.我们把两种情况的拟合曲线画在同一个坐标系中，很容易作出比较，见图6。第二种情形后半段的变化趋势与原始数据更吻合，因此，对将來人口的预测应该更好。300图7Logistic模型优化拟合曲线经过修改，得到了一个较满意的结果，人口增长率1=0.0270,极限人口Xm=345.1439（白万），并预测1990年2010年美国人口。通过调查得知1990-2010年人口数量统计如下表表61990年-2010年Logistic模型美国人口真实值与预测值时间（年）199020002010人

17、口248.7281.4308.7预测值249.6253.5280.0误差0.927.928.66.3问6.3.1马尔萨斯（Malthus）人口指数増长模型求解中国人口问题参照问题一，我们来求解中国人口问题。表8中国每10年自然增长率序号123456789r0.0560.0560.0560.0730.0070-0.085-0.050.028101112131115161718190.0330.0530.0580.1160.0360.0610.0540.220.2270.199同样，取高、中、低三种自然增长率方案，高方案中自然增长率为0.022,中方案中自然增长率为0.005,低方案中自然增长率为

18、0.001o图8r二0.022时马尔萨斯模型曲线拟合图9r二0.005时马尔萨斯模型曲线拟合图10r二0001时马尔萨斯模型曲线拟合根据上述分析，及曲线拟合可知，取中方案即1=0.005时，马尔萨斯模型更符合实际情况。因此本文自然增长率取1=0.005來预测中国人口数量并与实际情况对比。由预测公式预测1790-1980年的人口数量，由指数增长模型可得各个年份的真实值与预测值之间的差别如下表：表91790年4980年中国人口真实值与预测值年份1790180018101820183018401850186018701880真实人1(xlO6)323.5341.6360.83814094124123

19、77358368误差(xlO6)061139240364548782110714871925年份1890190019101920193019401950196019701980真实人1(xlO6)380400423472489518.8546.8667.1818.3981.2误差(xlO6)2471314439825004639179449973124101543919228通过调查得知1990-2010年人口数量统计如下表表101990年-2010年中国人口真实值与预测值时间(年)199020002010人口1135.21264.11341.0误差(xlO6)239872996537480图

20、11中国人口真实值与预测值曲线拟合6.3.2Matlab中cftool()工具箱求解应用Matlab中cftool()C具箱來进行图像拟合。Data.|Fitting.Exclude.|Plotting.Analysis.图12cftool工具箱拟合经过不同模型的应用比对，我们发现应用3阶高斯分布可以达到较好的拟合效果,误差相对较小。此时，al=1002(663.7,1339)bl=1992(1978,2006)cl=34.03(9.35,58.71)a2=307.6(58.96,556.3)b2=1932(1907,1958)c2=36.64(13.75,59.53)a3=403.7(391

21、.6,415.8)b3=1837(182&1846)c3=93.28(68.17,11&4)七、结果分析综合做出假设的两种模型与原始数据所描述的图形如下：可以看出，当世界人口总数不大时，生存空间，资源等极充裕，人口总数指数的增长是可能的，但当人口总数非常大时，指数增长的线性模型则不会反映这样的现实。阻滞增长模型从一定程度上克服了指数增长模型的不足，而抬数增长模型在做人口的短期预测时因为其形式的相对简单性也常被采用，但是依据实际情况，根据己给数据和查资料所得数据通过图像拟和可得二次函数模型与其最为匹配。事实上，人口的预测是一个相当复杂的问题，影响人口增长的因素除了人口基数与可利用资源量外，还和医药卫生条件的改善、人们生育观念的变化等因素有关，特别在做中短期预测时，我们希望得到满足一定预测精度的结果，比如在刚刚经历过战争或是由于在特定的历史条件下采纳了特殊的人口政策等，这些因素本身以及由此引起的年龄结构变化就变得相当重要，也要予以考虑。八、参考文献张志涌，杨祖櫻.MATLAB教程M.北京：北京航空航天人学出版社，2011.李晓梅.人II预测模型研究及应用M四川：西南财经大学岀版社，2011.姜启源等数学模型(第三版)M高等教育出版社2003.08P10何春.马尔萨斯人【I模型在广州市人口预测中的应用J.广东工业大学