数学建模人口增长模型.docx

1、数学建模人口增长模型数学建模人口增长模型摘要 :人口的增长是当前世界上引起普遍关注的问题作为世界上人口最多 的国家，我国的人口问题是十分突出的由于人口基数大尽管我国已经实行了 20 多年的计划生育政策人口的增长依然很快，巨大人口压力会给我国的社会 政治 经济医疗就业等带来了一系列的问题。 因此研究和解决人口问题在我国显得尤为 重要。我们经常在报刊上看见关于人口增长预报，说到本世纪，或下世纪中叶， 全世界的人口将达到多少亿。 你可能注意到不同报刊对同一时间人口的预报在数 字商场有较大的区别，这显然是由于用了不同的人口整张模型计算出来的结果。人类社会进入 20 世纪以来，在科学和技术和生产力飞速发

2、展的同时世界人 口也以空前的规模增长。 人口每增加十亿的时间， 有一百年缩短为十几年。 我们 赖以生存的地球已经携带着他的 60 亿子民踏入下一个世纪。长期以来，人类的繁殖一直在自然地进行着， 只是由于人口数量的迅速膨胀 和环境质量的急剧恶化， 人们才猛然醒悟， 开始研究人类和自然的关系、 人口数 量的变化规律以及如何惊醒人口控制等问题。本论文中有两个模型：（1）：中国人口的指数增长模型，并用该模型进行预测，与实际人口数据 进行比较。（2）：中国人口的 Logistic 图形，标出中国人口的实际统计数据进行比较。 而且利用 MATLAB 图形 ，标出中国人口的实际统计数据， 并画出两种模型的预

3、测曲线关键字： 人口预测；Malthus模型；Logistic模型；MATLAB软件、问题背景及重述1.1问题的背景中国是一个人口大国，人口问题始终是制约我国发展的关键因素之一。我国自1973年全面推行计划生育以来，生育率迅速下降，取得了举世瞩目的成就，但 全面建设小康社会仍面临着人口的形势和严峻挑战。 随着我国经济的发展、国家人口政策的实施，未来我国人口高峰期到底有多少人口， 专家学者们的预测结果 不一。因此，根据已有数据，运用数学建模的方法，对中国人口做出分析和预测 是一个重要问题。1.2问题的重述下表列出了中国19821998年的人口统计数据，取1982年为起始年（t=0），1982年的

4、人口 101654万人，人口自然增长率为14%。，以36亿作为我国人口的容 纳量，试建立一个较好的人口数学模型并给出相应的算法和程序， 并与实际人口进行比较时间（年）198219831984198519861987人口（万人）101654103008104357105851107507109300时间（年）198819891990199119921993人口（万人）111026112704114333115823117171118517时间（年）19941995199619971998人口（万人）119850121121122389123626124810问题分析对于人口增长的问题，其影响因素

5、有很多，比如：人口基数，出生率，死亡 率，人口男女比例，人口年龄结构的组成，人口的迁入率和迁出率，人口的生育率和生育模式，国家的医疗发展情况，国家的政治策略等众多的因素。如果把这 些因素都要考虑进去，则该问题根本无从下手。因此，应该根据中国人口自身发 展的特点，选取相应的能够体现我国人口发展特点的模型。人口发展模型有连续形式和离散形式，因为题目所给的数据是每个年份的具 体数据，可以将这些数据视为连续的。根据表格中的数据，我们使用 MATLAB编程（附录1）画出散点图。1.15 -1.1 -1.05 -图1中国19821998年的人口数据图从图中我们可以看到人口数在 1982-1998年是呈增长

6、趋势的，且增长趋势类似于指数型增长，因此，我们可以先建立一个指数增长模型（ Malthus模型）但是，由于地球上的资源是有限的， 它只能提供一定数量的生命生存所需的条件， 因此人口不可能无限制增加。 随着人口数量的增加， 自然资源， 环境条件等对人 口再增长的限制作用将越来越显著。于是我们假设在人口较少时， 可以把人口增长率看成常数， 但随着人口的增 加，我们应该把人口增长率视为一个随着人口增加而减小的量， 从而我们可以将 模型一（Malthus模型）优化为一个阻滞增长模型（Logistic模型）。三、 模型假设（一） Malthus模型假设我国人口的增长符合人口指数增长的规律，即满足 Mal

7、thus模型的两个前提：第一， 食物是人类生存所必需的；第二， 两性间的情欲是必然的，而且几乎会保持现状。从这两个“人类本性的固定法则” 出发，可以得出一个最基本的经济 比例：食物或生活资料的增长与人口的增殖之间的关系。马尔萨斯说， 人口的增殖比生活资料增长的要快，人口是按几何级数增长的，而生活资料则 只按算术级数增长。但是，马尔萨斯并不认为这两个级数就是人口规律的反映，他提出，保 持两个级数平衡的唯一出路就是抑制人口的增长。 他把所谓支配人类命运的 永恒的人口自然法则， 归纳成以下三个定理。三个定理：第一点是人口的制约原理，说明人口与生活资料之间必然存在某种正常的比例， 即“人口的增长，必然

8、要受到生活资料的限制”；第二点是人口的增殖原理，即“生活资料增加， 人口也常随着增加” ；第三点是马尔萨斯人口原理的核心， 称之为人口的均衡原理， 即“占 优势的人口繁殖力为贫困和罪恶所抑制， 因而使现实的人口得以与生活资料 保持平衡”。这个原理与前两个原理是紧密相连的，它说明人口与生活资料 之间最终将实现均衡， 但是这种均衡不是自然实现的，而是种种“抑制” 的产物。所以，Malthus模型假设条件如下：1.设P(t)表示t时刻的人口数，且P(t)连续可微。2人口的增长率 r 是常数(增长率=出生率-死亡率)。3.人口数量的变化是封闭，即人口数量的增加与减少只取决于人口中个体 的生育与死亡，且

9、每一个都具有同样的的生育能力与死亡率。(二) Logistic 模型由于地球上的资源有限， 当人口数量发展到一定阶段后， 会产生一系列问题， 如食物短缺、居住和交通拥挤等。另外，随着人口密度的增加，疾病将会增 多，死亡率会上升，因此，人口的增长率不会是 Malthus所假设的是一个常数不改变，而是会随着人口数量增加而减少。假设增长率 r表示P(t)的函数r(p),且r(p)为p的减函数。1.设r(p)为p的线性函数，r(p)=r-kp。2.自然资源与环境条件所能容纳的最大人口数为 Pm,即当P=Pm时，增长率 r(p)=0。四、变量说明符号表示意乂P人口数量t年份r人口自然增长率pm人口最大容

10、纳量p起始年人口五、模型建立与模型求解5.1 Malthus模型由分离变量法解得模型的解 为p(t) p0 ? ert对该模型两边同时取对数得到一次线性拟合函数y lnp,a ln p0, y rt a取表中 1982到 1998年的数据在 MATLAB 中 M 文件(附录 2)进行线性最小 乘拟合可得出 :f =0.013141 t - 14.5121所以可知 r=0.013141,p(t)=101654*exp(0.013141*(t-1982)用 MATLAB 进行指数拟合得到下图可以看出拟合曲线基本吻合，但是随着时间t的增加其误差逐渐加大，所以需要 对其修正。5.2 Logistic

11、模型由假设二可知，记p(t)是第t年的人口数量，人口增长率r(p)是p的线性函dpdt数，r(p)=r-kp。最大人口容纳量为Pm。即当P=Pm时，增长率r(p)=0。所以，r(p)?p rp(1 )Pm (5.2.1)P（t。） p。同样利用分离变量法求得其解p(t)Pmdpp根据（5.2.1）式作出dt 的曲线图（图1）以及由（5.2.2式作出p-t曲线图（图2）从上述曲线图以及表达式中，我们可以总结出如下规律:lim P(t) pmt ，它表明不管人口初始状态是什么样，人口总数最终都将趋于最大人口容纳量。dp dp当p(t)pm时，dt 0 ；当p(t)0。它表明当人口数量超过最大人口容

12、纳量时，人口数量将减少，当人口数量小于最大人口容纳量时， 人口数量将增 加。dp p 二人口变化率dt在 2时取到最大值，即人口总数达到极限值一半之前是加速生长的，经过此点后，增长率会逐渐减小至 0。采用非线性最小二乘估计法对参数 r和pm进行估计，通过使用matlab编写程序(附录4)可得:r =0.01137, pm =3.7465e+04用MATLAB拟合图像如下图3六、模型检验及结果分析经过前面模型建立的工作，已建立出Malthus模型和Logisic模型。现在根据所建 立的模型预测相关年份的人口数量，并与实际人口数量相比较以检验模型的优劣 性。Malthus模型与Logistic模型

13、对我国人数据的拟合结果年份实际人口计算人口计算人口/万P1P219821016541016501016501983103008103000102500198410435710436010335019851058511057401042019861075071071401050601987109300108560105920198811102610999010679019891127041114501076601990114333112920108530199111582311442010941019921171711159301102901993118517117460111180199411

14、98501190201120701995121121120590112960199612238912219011386019971236261238001147701998124810125440115670对表中数据进行分析可知：对于短期预测，这两个模型基本一致，但使用模型一 精选文档更简单；对于中长期预测，模型二要强于模型一。七、 模型评价与推广一、 优点：首先我们采用图表结合法， 比较直观地表达出题中所给的信息， 并据此得出 了人口增长的基本规律。根据所给出的数据，对其进行分析得出了人口增长率与人口总数的线性关 系，从而建立了人口阻滞增长模型， 对未来人口数的预测作出了较为准确的判断。模

15、型一是依据英国神父T Malthus的发现建立了指数型增长模型，经过我 们实际数据的检验，发现其人口早期的增长情况与 Malthus模型的预测基本相符， 然而随着时间的增加， 该模型的预测结果明显出现了不合理性。 其原因就是我们 将人口增长率视为常数，因此需要对 r进行修正。所以，我们将r表示为p的减 函数，从而推导建立了模型二二、 缺点：本文对模型一中的参数只做了线性估计，所以其计算结果与实际误差较大模型二中仅考虑了 r 与 p 的关系是线性的，没有考虑非线性关系八、 参考文献1.司守奎，孙兆亮，孙玺菁，周刚，仲维杰，康淑瑰 .数学建模算法与应用 （第二版） .国防工业出版社， 2016年2

16、.姜启源，谢金星，叶俊数学模型（第四版）高等教育出版社北京.20113.储昌木，沈长春.数学建模及其应用西南交通大学出版社成都.2015年4胡守信，李柏年.基于MATLAB的数学实验M.北京科学出版社.2004年6月5扬启帆，康旭升，等.数学模型M.北京：高等教育出版社.2006年5月6于学军.中国人口科学 2000年第 2期，时间： 2000-4-6，中国人口信息网.附录：1)symsx yx0=1982:1:1998y0=101654 103008 104357 105851 107507 109300 111026 112704114333 115823 117171 118517 11

17、9850 121121 122389 123626 124810xlabel(x)ylabel(y)plot(x0,y0,*)2)t=1982:1:1998;y=log(101654) log(103008) log(104357) log(105851) log(107507) log(109300)log(111026) log(112704) log(114333) log(115823) log(117171) log(118517) log(119850) log(121121) log(122389) log(123626) log(124810)p1=polyfit(t,y,1)

18、;f=poly2str(p1,t)3)syms x y px0=1982:1:1998y0=101654 103008 104357 105851 107507 109300 111026 112704114333 115823 117171 118517 119850 121121 122389 123626 124810xlabel(x) ylabel(y)plot(x0,y0,*)hold ont=1982:1:1998p=101654.*exp(0.013141.*(t-1982)plot(x,p,r,LineWidth,0.5);lege nd(原始数据散点图,指数拟合曲线);grid on;4)clc, cleara=textread(data1.txt);p=a(2:2:6,:);p=nonzeros(p);t=1982:1:1998;t0=t(1); p0=p(1);fun=(cs,td)cs(1)./(1+(cs(1)/p0-1)*exp(-cs(2)*(td-t0);cs=lsqcurvefit(fun,rand(2,1),t(2:end),p(2:end),zeros(2,1); r=cs(2),pm=cs(1)