数学建模人口增长模型

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

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

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

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

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

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

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

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

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

10、lthus 模型由假设一，t时刻到t十A t时刻人口增量为P(t t) - P(t) = r P(t) tdp于是可得dt二 rpp(t0) = p0由分离变量法解得模型的解 为p(t): poert对该模型两边同时取对数得到一次线性拟合函数y= lnpa= lnp,y= rt a取表中1982到1998年的数据在MATLAB M文件(附录2)进行线性最小二乘 拟合可得出：f =0.013141 t -14.5121所以可知 r=0.013141,p=101654*exp(0.013141*(t-1982)用MATLA进行指数拟合得到下图图2可以看出拟合曲线基本吻合，但是随着时间t的增加其误差

11、逐渐加大，所以需要 对其修正。5.2 Logistic 模型由假设二可知，记p(t)是第t年的人口数量，人口增长率r(p)是p的线性 函数，r(p)=r-kp o最大人口容纳量为Pm即当P=Pnfl寸，增长率r(p)=0。所以,dpp、r(p) p = rp(1 )dtpm (5.2.1)p(t0) = p0同样利用分离变量法求得其解Dmp(t)二一np.pm1 (1) exp(- r (t -1。)p0(5.2.2)dpp根据(5.2.1 )式作出出 的曲线图(图1)以及由(5.2.2 )式作出p-t曲线 图(图2)图2 p-t曲线图图1dpp曲线图 dt从上述曲线图以及表达式中，我们可以总

12、结出如下规律:lim Rd); pmtT00,它表明不管人口初始状态是什么样，人口总数最终都将趋于最大人口容纳量。dpdp当p(t)>pm时，出<0;当p(t)<pm时，出>0。它表明当人口数量超过最大人口容纳量时，人口数量将减少，当人口数量小于最大人口容纳量时，人口数量将增加。dp p = pm人口变化率出在 2时取到最大值，即人口总数达到极限值一半之前是加速 生长的，经过此点后，增长率会逐渐减小至 00采用非线性最小二乘估计法对参数 r和pm进行估计，通过使用matlab编写程序（附录4）可得：r =0.01137 , pm =3.7465e+04用MATLA版合图

13、像如下图3六、模型检验及结果分析经过前面模型建立的工作，已建立出Malthus模型和Logisic模型。现在根据所 建立的模型预测相关年份的人口数量，并与实际人口数量相比较以检验模型的优 劣性。Malthus模型与Logistic模型对我国人数据的拟合结果年份 实际人口 计算人口计算人口/万P1P2198210165410165010165019831030081030001025001984104357104360103350198510585110574010420198610750710714010506019871093001085601059201988111026109990106

14、7901989112704111450107660199011433311292010853019911158231144201094101992117171115930110290199311851711746011118019941198501190201120701995121121120590112960199612238912219011386019971236261238001147701998124810125440115670对表中数据进行分析可知：对于短期预测，这两个模型基本一致，但使用模型一 更简单；对于中长期预测，模型二要强于模型一。七、模型评价与推广一、优点：首先我们采

15、用图表结合法，比较直观地表达出题中所给的信息， 并据此得出 了人口增长的基本规律。根据所给出的数据，对其进行分析得出了人口增长率与人口总数的线性关系，从而建立了人口阻滞增长模型，对未来人口数的预测作出了较为准确的判断。模型一是依据英国神父T- Malthus的发现建立了指数型增长模型，经过我们实际数据的检验，发现其人口早期的增长情况与Malthus 模型的预测基本相符， 然而随着时间的增加，该模型的预测结果明显出现了不合理性。其原因就是我们将人口增长率视为常数，因此需要对r 进行修正。所以，我们将r 表示为 p的减函数，从而推导建立了模型二二、缺点：本文对模型一中的参数只做了线性估计，所以其计

16、算结果与实际误差较大模型二中仅考虑了r 与 p 的关系是线性的，没有考虑非线性关系八、参考文献1 . 司守奎，孙兆亮，孙玺菁，周刚，仲维杰，康淑瑰. 数学建模算法与应用（第二版）. 国防工业出版社，2016年2 .姜启源，谢金星，叶俊.数学模型（第四版）.高等教育出版社北京.2011年3 .储昌木，沈长春.数学建模及其应用.西南交通大学出版社成都.2015年4胡守信，李柏年.基于MATLAB勺数学实验M.北京科学出版社.2004年6月5 扬启帆，康旭升，等. 数学模型 M. 北京：高等教育出版社.2006 年 5月6 于学军. 中国人口科学2000 年第 2 期，时间：2000-4-6 ，中国

17、人口信息网 .附录：1) syms xyx0=1982:1:1998y0=101654 103008 104357 105851 107507 109300 111026 112704114333 115823 117171 118517 119850 121121 122389 123626124810xlabel( 'x' )ylabel( 'y' )plot(x0,y0, '*' )2) t=1982:1:1998;y=log(101654) log(103008) log(104357) log(105851) log(107507)lo

18、g(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);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 123626124810xlabel('x')ylabel('y')plot(x0,y0,'*')hold ont=1982:1:1998p=101654.*exp(0.013141.*(t-1982)plot(x,p,'r','LineWidth',0.