数学建模 人口模型 人口预测

1、关于计划生育政策调整对人口数量、结构及其影响的研究【摘要】本文着重于讨论两个问题:1、从目前中国人口现状出发，对于中国未来人口数量进行预测。2、针对深圳市讨论单独二胎政策对未来人口数量、结构及其对教育、劳动力供给与就业、养老等方面的影响。对于问题1从中国的实际情况和人口增长的特点出发，针对中国未来人口的老龄化、出生人口性别比以及乡村人口城镇化等，提出了Logistic、灰色预测、等方法进行建模预测。首先，本文建立了Logistic阻滞增长模型，在最简单的假设下，依照中国人口的历史数据，运用线形最小二乘法对其进行拟合，对2014至2040年的人口数目进行了预测，得出在2040年时，中国人口有14

2、.32亿。在此模型中，由于并没有考虑人口的年龄、出生人数男女比例等因素，只是粗略的进行了预测，所以只对中短期人口做了预测，理论上很好，实用性不强，有一定的局限性。然后，为了减少人口的出生和死亡这些随机事件对预测的影响，本文建立了GM(1,1)灰色预测模型，对2014至2040年的人口数目进行了预测，同时还用2002至2013年的人口数据对模型进行了误差检验，结果表明，此模型的精度较高，适合中长期的预测，得出2040年时，中国人口有14.22亿。与阻滞增长模型相同，本模型也没有考虑年龄一类的因素，只是做出了人口总数的预测，没有进一步深入。对于问题2针对深圳市人口结构中非户籍人口比重大，流动人口多

3、这一特点，我们采用了灰色GM(1,1)模型，通过matlab对深圳市自2001至2010年的数据进行拟合，发现其人口变化近似呈线性增长，线性相关系数高达0.99，我们就此认定其为线性相关并给出线性方程。同理，针对其非户籍人口，我们进行matlab拟合发现，其为非线性相关，并得出相关函数。并做出了拟合函数X(t+1)二17255.8Xe0.0419775xt-16531.2。对于新政策的实施，我们做出了两个假设。在假设只有出生率改变的情况，人口呈现一次函数线性增加。并拟合出一次函数Y=17965.0Xe0.0327356xt-17372.5；在假设人口增长率增长20%时，做出了预测如果单独二胎政

4、策实施，到2021年，深圳市常住人口数将会到达1137.98千万人。关键词:GM(1,1)灰色模型Logistic阻滞增长模型线性拟合非线性拟合目录】TOC o 1-5 h z HYPERLINK l bookmark8 o Current Document 一、问题重述(4) HYPERLINK l bookmark10 o Current Document 二、符号定义与说明(4)三、模型假设(4)四、问题分析及模型建立及求解A、问题一:1、问题背景(5)2、问题分析(5)3、模型建立模型一:阻滞增长模型的建立(6)阻滞增长模型的求解(6)阻滞增长模型的分析(7)阻滞增长模型的优化(7)阻

5、滞增长模型优化后的分析(9)模型二:GM(11)灰色预测模型的建立(9)GM(11)灰色预测模型的求解(10)GM(11)灰色预测模型的分析(11)B、问题二:1、问题重述(11)2、问题假设(11)3、问题背景(12)4、灰色预测模型的建立(14)5、灰色预测模型的求解(14)6、模型的优化(新政策实施后的预测)(15)新政策下的建模(16)假设拟合成一次线性函数(16)假设按比例增长(17)数据分析及评价(17)五、模型总评价(18)六、参考文献(19)七、附录(19)问题重述】人口的数量和结构是影响经济社会发展的重要因素。从20世纪70年代后期以来，我国鼓励晚婚晚育，提倡一对夫妻生育一个

6、孩子。此计划生育是我国的一项基本国策。近年来中国的人口发展出现了一些新的特点，例如，老龄化进程加速、出生人口性别比持续升高，以及乡村人口城镇化等因素。党的十八届三中全会提出了开放单独二孩，今年以来许多省、市、自治区相继出台了具体的政策。收集一些典型的研究评论报告，根据每十年一次的全国人口普查数据，建立模型，对报告的假设和某些结论发表自己的独立见解，并针对深圳市或其他某个区域，讨论计划生育新政策（可综合考虑城镇化、延迟退休年龄、养老金统筹等政策因素，但只须选择某一方面作重点讨论）对未来人口数量、结构及其对教育、劳动力供给与就业、养老等方面的影响。符号定义与说明】符号符号定义单位或备注x(t)t年

7、的全国人口总数亿r(t)人口增长率函数(t)t年人口死亡率%0X总人口随时间变化的拟合函数z（D（k）白化背景值b灰作用量Y数据向量B数据矩阵u参数向量【模型假设】不考虑我国人口向国外搬迁，同时也不考虑外国人口向国内搬迁不考虑战争、灾害、疾病对人口数目的影响；假设一年内，各个地区，各个年龄段的死亡率不会发生变化；假设在一年内，处于生育年龄的妇女生育率不会发生变化；假设附件中所给数据真实可靠具有预测性；假设影响中国总人口数的主要因素是死亡率和出生率。问题分析】问题一:对中国未来的人口数量进行预测。一、问题背景中国是世界上人口最多的发展中国家，人口多，底子薄，人均耕地少，人均占有资源相对不足，是我

8、国的基本国情，人口问题一直是制约中国经济发展的首要因素。人口数量、质量和年龄分布直接影响一个地区的经济发展、资源配置、社会保障、社会稳定和城市活力。在我国现代化进程中，必须实现人口与经济、社会、资源、环境协调发展和可持续发展，进一步控制人口数量，提高人口质量，改善人口结构。对此，单纯的人口数量控制（如已实施多年的计划生育）不能体现人口规划的科学性。政府部门需要更详细、更系统的人口分析技术，为人口发展策略的制定提供指导和依据。长期以来，对人口年龄结构的研究仅限于粗线条的定性分析，只能预测年龄结构分布的大致范围，无法用于分析年龄结构的具体形态。随着对人口规划精准度要求的提高，通过数学方法来定量计算

9、各种人口指数的方法日益受到重视，这就是人口控制和预测。二、问题分析本题需要结合中国的实际情况和人口增长的特点来对中国人口增长的中短期和长期趋势做出预测。首先，我们从简单模型入手，利用已有年份的人口总量数据预测将来的人口总量的变化趋势，从总体上对人口发展做出预测。其次，把人口的增长特点考虑在内，利用动态模型并进行计算机模拟，得到符合中国实际情况的模型，包含了老龄化水平、性别比例、城镇化等更细致的结果。最后，我们对每个模型的预测结果进行对比，评判其各自的优点及缺点，并对政府部门提出一些建设性的意见。三、模型建立及求解1.31模型建立A.阻滞增长模型针对未来中国人口总数，我们建立简单的预测模型-阻滞

10、增长模型。（具体建立方法见附录1）我们可以得到以下等式：x(1)人口增长率函数：r（x）二r（1）xmt时刻人口数目函数：x（t）=一1+Ix0 xm-m1e_rt丿(2)dx/(3)单位时间内人口增量方程:一dt=rst,s=-xxmA模型求解首先，利用方程（3）以及2002年到2012年的数据用MATLAB软件对方程（3）进行线性最小二乘法拟合。得到r=0.0524,x=14.753m表1:实际数据和计算数据对比表年份200220032004200520062007200820092010实际人口（亿）12.8412.9212.9913.0713.1413.2113.2813.3413.4

11、0计算人口（亿）12.8412.9213.0013.0813.1613.2313.3013.3713.43年份201120122013实际人口（亿）13.4713.5413.60计算人口（亿）13.4913.5513.61图1:阻滞增长模型拟合中国人散点图（蓝色为计算数据，红色为实际数据）data1data213.813.713.613.5亿奴口人13.413.313.213.11312.91234567891011相对于2002年（为0）的年份图1模型分析：可以看出组织增长的中段与实际数据不太符合。模型优化:在上述模型中，我们把2002年到2013年的数据参与函数的拟合，目的是为了用它们作模

12、型的检验:我们用模型计算的数据和这一算时间的实际数据相比较，来检验模型是否是合适。经计算得至到20022013年人口计算与实际数据的相对误差如表2:年份200220032004200520062007200820092010相对误差（）00.040.130.130.170.190.190.240.27年份201120122013相对误差（）0.210.120.09表220022013年人口计算与实际数据的相对误差经过上面的分析求解，我们应当把2002年到2013年的及实际数据进行重新拟合，得到新的(4)/、14.537x(t)如下:x(t)=】(14.5371+(-1)e-0.0553t12.

13、84用（4）对未来人数进行预测，人口数据如表3年份201420152016201720182019202020212022预测人口总数（亿）13.6113.6613.7013.7413.7813.8213.8613.8913.92年份202320242025202620272028202920302031预测人口总数（亿）13.9613.9914.0214.0414.0714.0914.1214.1414.16年份203220332034203520362037203820392040预测人口总数（亿）14.2014.2214.2314.2514.2614.2814.2914.3114.32表

14、3阻滞增长模型预测20142040年的人口数据下图是优化后阻滞增长模型拟合中国人散点图。重新拟合后计算人口数和实际人口数对比图13.7.13.6-LLLLdata1data2-13.5-813.4-OO13.3-013.2“013.1-0013-12.9-12.8rrr2002200420062008201020122014年份图2优化后阻滞增长模型拟合中国人散点图亿/数口人优化后的模型分析:对于优化后的模型，在2002到2010年间坐到了高度符合。但是在2010年后的数据的拟合中还是出现一些小问题。1.3.2模型二:GM（11）灰色预测模型1.3.2.1模型建立由于人的出生和死亡是随机的，因

15、为我们利用灰色预测模型中的累加效果，尽量减少这种随机的影响，在此我们采用GM（l.l）灰色预测模型（具体参考附录）。为了使预测结果效果更佳，并不直接用总人口序列建模，而是先求出各年的净增人口序列，即20022013年各年的净增人口数据如表4，然后应用净增人口序列建模计算净增人口预测值。年份200220032004200520062007200820092010人口总数（亿）12.8412.9212.9913.0713.1413.2113.2813.3413.40净增人口数（万）827774761768692681673648641年份201120122013人口总数（亿）13.4713.541

16、3.60净增人口数（万）644669668表42002年2013年各年净增人口数据。1.3.2.2模型求解1设X(0)=(x(o)(1),x(0)(2),L,x(o)(n),X=(x(l),x(2),L,x(n)称X(0)(k)+ax(k)(k)=b为GM(1,1)模型的原始形式。则此题中原始样本序列为:q(0)二827774761768692681673648641644669668通过灰色GM(l.l)模型，我们得到未来t年的较2001年的累计人口净增量(记2002年t=1,2003年t=2,之后以此类推)q(1)(t)=32954.4一313254e-0.049579(t-1)(5)利用

17、q(0)(t)=q(1)(t)一q(t一1)以及(5)式，得到人口的净增量x(0)(t)后,x(t)=q(0)(t)一q(t一1)利用迭代可得出未来第t年的人口数量记2002年t=1，2003年t=2，之后以此类推)为x(t)=145660313254e-0.049579(t-i)(6)由(6)式预测从2014年到2040年的人口预测如下表5年份200520062007200820092010201120122013预测人口数(亿)13.0713.1413.2113.2813.3413.4013.4613.5113.57年份20142015201620172018201920202021202

18、2预测人口数(亿)13.6113.6313.6713.6913.7113.7413.7613.7913.81年份202320242025202620272028202920302031预测人口数(亿)13.8313.8613.8813.9013.9213.9513.9713.9814.01年份203220332034203520362037203820392040预测人口数(亿)14.0314.0614.0814.1114.1314.1514.1714.1914.22表5灰色GM(1.1)模型预测20052040年的人口数据1.3.2.3结果分析对我国2010年实际总人口数进行检验性预测，预测

19、值为13.40亿，实际值为13.4091亿，相对误差为0.06%。然后进行后误差检验,检验合格后的模型即可用于预测。对20052013年数据进行检验，得到19902005年人口计算与实际数据的相对误差(如表6)年份200520062007200820092010201120122013相对误差()00.040.030.020.0010.010.020.0030.002表620052013年人口计算与实际数据的相对误差S1=246.9858,S2=8.0308,C=S2/S1=0.0325v0.35，表明模型预测的精度很高。根据灰色系统理论，当发展系数ae(-2,2)且a-0.3时，所建GM(1

20、,1)模型可用于中长期预测。但是由于中国的人口发展特点并不能保证很长时间预测的准确性,比如对中国100年后的人口数量便难以预测。问题二计划生育新政策对深圳人口数量的影响一、问题重述人口问题有着悠久的研究历史,也有不少经典的理论和模型。这些理论和模型都依赖生育模式、生育率、死亡率和性别比等多个因素。这些因素与政策及人的观念、社会文化习俗有着紧密的关系,后者又受社会经济发展水平的影响。研究中用到的数据的置信水平也与调查统计有关。根据每十年一次的全国人口普查数据,建立模型,对报告的假设和某些结论发表自己的独立见解,并针对深圳市或其他某个区域,讨论计划生育新政策对未来人口数量、结构及其对教育、劳动力供

21、给与就业、养老等方面的影响。二、模型假设不考虑深圳市的迁入和迁出情况,只考虑常住人口不考虑战争、灾害、疾病等对人口的影响假设一年内,各个地区,各个年龄的死亡率不变化三、问题背景1、对于深圳的人口，做如下分析普查资料显示：全市常住人口中，014岁人口为10188万人，占984；1564岁人口为91564万人，占8840；65岁及以上人口为1828万人，占176。深圳人口年龄分布饼图，如图3图3深圳人口年龄分布饼图深圳仍属于年轻人城市，据初步统计测算，全市人口平均年龄为30岁左右；深圳还处于旺盛的“人口红利”期，1564岁人口占总人口的8840，比广东省的7636高出12个百分点，比上海市的812

22、5高出72个百分点，比北京市的8270高出57个百分点；按照国家65岁及以上人口占7以上即达到老年社会的标准，深圳离进入老龄化社会还较远。2、模型的建立与求解人民网强国论坛近日针对“单独家庭”做了一项生育意愿的调查。结果显示，在符合“单独两孩”条件者中，愿意生第二孩的超过一半，达50.73%。其中，50.41%的人表示新政一实施就生，49.59%的人表示过几年再生。与此同时，不愿生两孩的占1/3,为33.46%。不愿意生两孩的原因中，“经济成本太高”是主要因素，占56.16%。3、在深圳新闻网发起的“你会选择生二胎吗？”的调查中，截至昨晚7时，有76.08%的受访网友表示会生二胎，表示不会的网

23、友有23.92%。其中57.81%的受访网友都选择在“3040岁选择生二胎”。此外，有87.24%的网友表示生二胎是因为“觉得一个小孩太孤单了”，而在选择不会生二胎的网友中，60.18%的网友是因为“养不起”。4、通过收集的数据分析运用Excel画出1979年到2012年的常住人口的变化如图4。1979年-2012年深圳常住人口数目变化曲线1J-50goodC0IW803662喜份862年C0G2526旳2影2S2拓2662.2图41979-2012年深圳常住人口数目变化曲线通过这张图表可看出，深圳市的常住人口在2011年已经突破了1000万，当政策新出之后，势必会对这样的一种平稳发展的格局有

24、一定的冲击，或多或少，都会使这样的数量有所增长。在此，收集到一些出生率和死亡率的图表1979-2012深圳人口出生率和死亡率变化曲线图3020107V+/*4*I*f十+出生率-死亡奉_.丁P1|uiiii巧七丁|口|=巧|rB-1|200000000006666600066666666666000I_17=I_1t=4I_1I1Ii=1r(21年份图5深圳市人口出生率和死亡率变化曲线分析:从这张图可以看出，最近几年的深圳，出生率在不断提高，而死亡率在不断下降，从而导致了人口的激增。但是，新政策的实行，对于出生率有促进作用，死亡率却不会有明显变化，所以，人口将会呈现更快的增长，但是这样的增长是

25、在前期，随着时间推移累计效应释放后，波动在1.61.7四、模型建立首先，我们采用灰色GM(1,1)模型采用原因：灰色模型适用于小样本、贫信息、内在规律未充分外露的系统，按适当办法处理原始数据后得到规律性较强的生成函数本题给出的常住人口、非常住人口数据受到难以区分的多重因素影响，且数据量较小，适用于灰色模型.由于常住人口数量受历史影响较大，不易发生较大变化，且在数据处理中发现了较强的线性关系，我们之后采用了一元线性拟合来简化模型；而非常住人口受各方面因素影响较大，仍保持灰色模型不变。最后我们得出了与常住/非常住人口相关的人口结构变化规律。灰色模型的详细做法见问题一以及附录此预测模型是拟合参数模型

26、，通过原始数据累加生成，得到规律性较强的序列，用函数曲线拟合得到预测值.建立过程如下：设原始数据序列X(0)有n个观察值，X(0)二X(o)(l),X(o)(2),L,X(o)(n)，通过累加生成新序列X(0)二X(1)(1),X(1)(2),L,X(1)(n)，利用新生成的序列X拟合函数曲线.利用拟合出的函数求出新生序列X的预测值序列X.利用X(k)二X(k)-X(k-1)累减还原，得到灰色预测值序列(0)(1)(1)X二X(1),X(2),L,X(n+m)(共n+m个，m个未来预测值).将序列X(o)分0000为Y和Z，其中Y反映X(0)的确定性增长趋势，Z反映X(0)的平稳周期变化趋00

27、00势.对X(0)序列的确定增长趋势进行预测.五、模型求解整理得深圳市2001年2010年常住人口数，见下表7.表7：深圳市20012010年年末常住人口数根据上述数据建立含有10个观察值的原始数据序列X(o):X(0)=【724.57746.62778.27L954.28995.011037.2使用Matlab软件对X(0)进行一次累加，得到新数列X，见表8.表8：GM(1,1)算法拟合值及误差序号年份模型值残差相对误差级比偏差X(1)(2)2001724.5700X(1)(3)2002739.770.00920.123%-0.0121X2003771.490.00870.112%-0.00

28、04X(1)(5)2004804.56-0.00470.059%-0.0135X(1)(6)2005839.06-0.01340.1650%0.0089X(1)(7)2006875.03-0.00450.0518%0.0090X(1)(8)2007912.54-0.00020.0020%0.0043X(1)(9)2008951.660.00270.0287%0.0029X(1)(10)2009992.460.00260.0257%-0.0002X(1(11)20101035.010.00210.0204%-0.0005拟合函数：X(t+1)二17255.8Xe0.0419775xt-16531

29、.2六、模型优化A.假设新政策推出后，对于人口的变化很小，除了出生率改变，其他条件都不改变。猜想新政策实施后，人口按线性增长，即用简单的一次函数进行拟合。图6由图可见数据的线性关系很强，且在一段时间内仍保持线性增长趋势。由最小二乘法得拟合函数为X二35.215t-69759，相关系数为r二0.99388。常住人口拟合函数为Y=17965.0 xeo.o327356xt-17372.5。在这种情况下，深圳人口预测变化如下图:图7B.假设实行新政策后，出身率有所提升，并且，出生率提高了20%,这样的话，在未来的几年，在人口数量，人口结构分布都会发生一定的影响，深圳市的总体人口也将相应地出现一些变化

30、。当出生率提高20后，后10年的常住人口数表如下年份常住人口数（万人）20141075.1220151083.3420161092.6420171102.2320181111.2320191120.6720201128.7820211137.98表9变化曲线如下深圳市未来10年人数预测ooooooO6420864-I亠*亠*亠4亠-I亠1亠1亠4亠1亠1亠1亠4亠20-亠2062011820份17年20220图8七、数据分析在“单独二孩”政策落地之前，省卫计委组织了全国著名专家来广东调研，对于实施“单独二孩”后对我省人口会有那些影响做了评估。调查认为，广东放开“单独二孩”政策后对出生人口的数量

31、影响并不大，对总人口的影响也微乎其微。2010年广东全省独生子女人口总数是732.25万，其中户籍人口中法定婚龄的独生子女家庭24.35万（男满22周岁，女满20周岁）。包括已婚单独未生育家庭5.1万，单独生育一孩的家庭14.59万，单独生育二孩的3.41万，单独生育三孩的1.25万。也就是说，全省符合“单独二孩”政策的家庭不超过15万。如果实施“单独二孩”政策后，2014年广东全省预计新增人口1-1.5万人，2015预计新增人口7万人，2016年预计新增8.19万人，2017年预计新增7.8万人，2018年预计新增6.96万人，呈逐渐平稳趋势。所以，由人口增长模型可以分析，如果婴儿出生率有所

32、增加，那人口的结构也将变化。在人口构成方面，由于现在的深圳还未进入老龄化，并且新政策的实施，使得婴幼儿的数量增加，使人口构成方面呈现出成长型的格局。可以看到，在不久的将来，会是一个稳定型的人口城市。人口的增加同时也给公共资源带来了压力，而且，对教育事业也是一个不小的挑战。“单独二胎”放开有望在5年内新增750万新生儿，这些新生儿在20152019年将拉动上万亿元的市场消费。从政策影响的先后时间来看，首先是母婴医药，其次是奶粉、尿不湿等初生婴儿用品相关产业，再次是儿童服饰、家具、童车、玩具，最后是动漫、钢琴、文具等文化教育相关产业。同时，这样的情况也使得就业岗位有所增加，在短期内解决了人们的就业

33、问题。但是长远来看，当这些孩子长大，不得不面临就业更加艰难的局面。【模型总评价】本文首先从最简单的情况入手，不考虑人口的年龄层次和男女出生比例，采用Logistic模型和灰色预测模型对人口的总数进行预测。由于这两个模型对于种群数量是单调改变的情况预测准确度很高，而人口总数的变化并不总是单调的，所以只能对人口行短期或中短期预测，长期预测结果只可作为参考。在第二问中，本文也从相对简单的灰色预测模型入手，对深圳在单独二胎政策不实施的情况下，未来人口数量进行了拟合估计。在对单独二胎政策实施后，深圳的人口数量变化做了预测。但是由于时间匆忙，对于深圳实施单独二胎之后的对深圳人口数量，结构以及教育、劳动力供

34、给与就业、养老等方面的影响没有深入讨论。【参考文献】数学模型姜启源高等教育出版社中国统计年鉴2013中华人民共和国国家统计局深圳统计年鉴2013深圳统计局数学建模全国赛07年A题一等奖论文【附录】1.模型一阻滞增长模型分析人口增长对一定数量后增长率下降的主要原因，人们注意到，自然资源、环境条件等因素对人口的增长起着阻滞作用，并且随着人口的增长，阻滞的作用越来越大。所谓阻滞增长模型就是考虑到这个因素，对指数增长模型的基本假设进行修改后得到的。阻滞作用体现在对人口增长率R的影响上，是的R随着人口数量X的增长而下降。若将R表示为X的函数r(x)，则它应是减函数。于是有方程写作dx=r(x)x,x(0

35、)=x0(1)dt0对r(x)的一个简单的假定是，设r(x)为x的线性函数，即r(x)=r一sx(r0,s0)(2)00这里的r称固有增长率，表示人口很少时(理论上是x=o)的增长率。为了确定系数s的意义，引入自然资源和环境条件所能容纳的最大人口数量X，称人口容量。当X=X时,mmr人口不再增长，即人口增长率r(x)二0，代入(2)式得s二，于是(2)式为mxmr（x）二r（1一上）（3）0 xm式的另一种解释是,增长率r(x)与人口尚未实现部分的比例(x-x)/x成正比,mm比例系数为固有增长率厂。将(3)代入(1)得dxx4）=rx（1一）,x（0）=xdt0 x0mx方程(4)右端的因子

36、rx体现人口自身的增长趋势，因子(1-)则体现了资源和环境对0 xm人口增长的阻滞作用。显然，x越大，前一因子越大，后一因子越小，人口增长是两个因子共同作用的结果。2灰色预测模型运用GM(1,1)模型、灰色Verhulst模型、离散灰色模型三个模型对深圳人口数量进行预测研究。3.2.1GM(1,1)模型定义3.2.1设X(0)二(x(0)(1),x(0)(2),L,x(o)(n),X二(xd)(1),x(2),L,x(d(n)称X(0)(k)+ax(k)(k)=b为GM(1,1)模型的原始形式。其中G表示灰色(grey),M表示模型(Model),第一个1表示一阶方程，第二l表示1个变量。GM(1，1)模型首先对原始数据进行一阶累加生成，然后利用指数曲线拟合并预测，最后通过累减还原得到预测值。一般将原始数据序列记为X(0)，将一阶累加生成序列记为X。建GM(1,1)模型的步骤如下:(1)假定原始数据序列为X(0)二(x(0)(1),x(0)(2),L,x(0)(n)对原始数据序列进行一阶累加生成X(1)=(x(1)(1),x