西南大学A彭作祥 陈嵩箐 王霞 李敏.docx

1、西南大学A彭作祥 陈嵩箐 王霞 李敏2007高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）：

2、A 我们的参赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）： 西南大学 参赛队员 (打印并签名) ：1. 陈嵩箐 2. 王 霞 3. 李 敏 指导教师或指导教师组负责人 (打印并签名)： 彭作祥 日期：2007年9 月 24 日赛区评阅编号（由赛区组委会评阅前进行编号）：教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：中国人口预测摘 要 人口的增长与育龄妇女的生育率成正比与死亡率成反比，同时受到

3、人口的性别比率，年龄结构，社会的经济水平以及人文历史各个方面因素的影响。本文以2005年中国人口数据为基础，建立人口模型，并根据模型预测未来50年中国人口的发展趋势。首先对一些参数作了适当的假设，我们预计中国的人口仍然有一个缓慢的增长过程。大约于2016年左右达到最大值13.4亿左右，从2017年起中国人口会有一个短暂的平衡时期，到2022年左右，人口将开始下降。最后到2055年中国的人口将会下降到10亿左右。为了实现我国人口的战略目标，并且能在促进人的全面发展的基础上，实现人口发展自身的协调和可持续，实现人口与经济社会资源环境之间的协调和良性互动，我们对模型进行了优化。在从提高人口素质，合理

4、布局人口以及缓解人口老龄化趋势和降低人口以外死亡率等各个方面考虑以后，我们重新对变化趋势作出假设，得到了一个相对理想的模型。预测结果表明，从2005年到2055年期间，总人口是先有一段持续上升的过程，2010年人口总数将达到13.54亿，至2022年达到峰值14.381亿，2023年至2039年，2040年左右人口总数将保持在一个相对稳定的状态。从2023年至2039年，人口总数将保持在一个相对稳定的水平。从2040年左右人口将出现回落，至2050 年大约将减少到 13.6亿。 关键词 总和生育率， 死亡率， 性别比， 人口预测模型.一 问题重述中国是一个人口大国，人口问题始终是制约我国发展的

5、关键因素之一。认识人口数量的变化规律，建立人口模型并做出准确的预测，是有效控制人口增长和保持人口稳定的前提。因此，根据已有数据，运用数学建模的方法，对中国人口做出分析和预测是一个重要问题。近年来我国的人口发展出现了一些新的特点，例如：老龄化进程加速，出生人口性别比持续升高，以及乡村人口城镇化等等，这些都影响着中国人口的增长。对此，2007年初发布的国家人口发展战略研究报告还做出了进一步的分析。请参考附录2中的相关数据，建立一个明确完整的人口增长模型，并指出求解方法。然后根据上述模型，对中国未来中短期和长期的人口增长趋势做出预测。最后对模型做出评价并加以改进和推广。二问题分析影响人口增长的因素多

6、种多样，必须从诸多因素中选择一些最重要的因素并从这些方面来描述人口增长。首先育龄妇女的生育率和人口的死亡率是直接影响人口增长的因素；其次新生儿的性别比率也是影响人口增长的重要因素。以上三种因素是我们重点讨论的对象。至于人口的老龄化趋势和乡村人口城镇化等因素，由于本身影响效应较小且数据有限，暂不纳入模型中加以考虑。以某年人口数据为起始，并以妇女生育率，人口的死亡率和新生儿的性别比率为参数，我们得到了以时间（单位年）为自变量的初步模型。由于所给数据未必是1%抽样调查的结果，不能从中得到全国总人口。我们从国家统计局下载到全国总人口如下：年份2001年2002年2003年2004年2005年全国总人口

7、（万人）127627128453129227129988130628 对所给数据分析发现，2003年的生育率异常，都不及其它年份的四分之一，大多只有十分之一，如下表。于是，用2002年与2004年的平均值代替2003年的相应值。城市生育率2001年2002年2003年2004年2005年15岁0.4500.0100.0216岁00.490.0400.2517岁0.1400.080.471.0718岁0.730.50.182.443.0719岁2.41.970.76.789.8320岁4.0710.531.8623.8329.0121岁23.4234.473.7339.7948.1722岁52.

8、7952.76.1860.5664.9523岁85.0779.368.9495.7784.1724岁99.8101.911.37117.7896.4525岁128.84111.111.51123.8596.5226岁137.75117.2410.65123.6692.4227岁112.3101.739.91108.1283.128岁86.4284.167.5386.5469.4829岁73.2274.035.5766.385530岁47.353.744.5149.3945.3431岁38.5133.383.2936.7735.4332岁23.931.162.3426.5527.0933岁27.

9、9814.821.921.5121.2834岁18.0918.031.3417.5417.535岁11.1511.150.9411.1712.4636岁7.528.550.77.88.9837岁7.276.290.526.266.3638岁4.294.670.395.114.5139岁3.82.490.214.333.3740岁2.271.620.192.952.741岁01.040.231.611.7642岁0.950.620.120.91.2343岁0.570.580.0800.944岁00.950.0600.6445岁00.510.0300.7646岁00.10.050.130.947岁

10、0.3600.030.130.548岁0.360.230.010.050.5249岁0.360.420.010.140.74三假设说明1 按国家统计局的统计，2001年全国总人口数为127627万，2005年全国总人口数为130628万；2 以每1岁为1个年龄组，1年为1个时间段；3 假设在预测年限内，不发生战争及其他重大自然灾害，不考虑大的人口迁移因素；4 假设抽样调查中的数据能完全反映总体的水平；5 假设每年生育子女的妇女只生一个子女，不考虑多胞胎的情形；6 假设对异常数据的处理不影响整体的水平。四符号说明：第年岁妇女的生育率；：第年男性人口数比女性人口数；：第年岁男性的死亡率；：第年岁女

11、性的死亡率；：第年岁男性人口的存活人数；：第年岁女性人口的存活人数；：第年岁男性占总人口的比重；：第年岁女性占总人口的比重；：第年人口总数。五模型的建立与求解1.模型的建立：为了简化模型并降低数据处理的难度，我们将全国总人口分成几大部分分别加以讨论。由计算知：年龄大于90岁的人口数仅占总人口的0.001，所以在模型中只考虑089岁的人，以预测全国城市男性人口的发展趋势为例，我们建立如下的模型。设第年的全国城市男性人口总数为：，以此类推，我们得到了下面的人口滚动模型： 同理可得城市女性，城镇男性，女性和乡村男性女性分别的人口发展的模型。将以上每年六种不同类别的人口数相加就得到了全国总人口的发展情

12、况的模型。2.模型的求解：首先，我们用人口总数，以2001年为起始年，利用2001年的全国人口数量，取生育率、死亡率、出生性别比为20012005年的平均值来预测未来十年人口发展的趋势，并把预测结果和实际情况相比较来检验模型的合理性。（下图）此时预测出2005年的总人口为131776万，与实际总人口130628万的误差为0.88% 。这说明本模型能较好的预测出人口总数的发展趋势。但是从长远来看，存在一定的误差，有待改进。再以2005年为起始年，利用05年人口总数。仍然假设妇女生育率，人口死亡率和出生性别比是2001年到2005年的平均值。（下图）2005-2015近十年内，人口逐步增长，在20

13、16年左右达到人口峰值为。到2055年人口数量下降到10亿以内。分别以最大、最小生育率来预测未来人口，并与平均生育率预测结果作比较得到下图：.上图表明：生育率水平提高时，峰值增大，人口峰值出现时间延迟。计算2005年总和生育率得2005年的总和生育率为1.3385，确实明显低于更替水平。由以上分析，模型有待改进。3.模型改进：由上面的分析，可以发现参数对结果影响很大。于是我们希望通过优化参数设置得到最优的模型。根据生育更替水平的标准，总和生育率越接近2.1，出生和死亡将逐渐趋于均衡，在没有大量人口迁移的情况下，总人口将趋于稳定状态。由上面的计算知：2005年的总和生育率为1.3385，低于更替

14、水平。于是我们对生育率水平作相应的调整，尽量使得总的生育率接近根替水平。此外，出生性别比在极大程度上，对人口的增长有重大的影响。专家分析:“一般来说，每出生100名女婴，其男婴出生数置于103107之间。”此分析明确认定了出生性别比的通常值域为103107之间。从此，出生性别比值下限不低于103、上限不超过107的值域一直被国际社会公认为通常理论值，其他值域则被视为异常。即正常情况下，出生性别比保持在103107是最好的。我们从以下方面加以改进。死亡率 考虑到以后科技水平和生活水平会提高，环境的改善，死亡率会有所下降，故死亡率取2001年到2005年的最小值。总和生育率 以05年数据为基准，适

15、当调高各个年龄阶段育龄妇女的生育率，使得2005年的总和生育率为2.1左右。然后取性别比为1.07，预测未来50年的结果。采用同样的方法，适当调整各个年龄阶段育龄妇女的生育率，使得2005年的总和生育率为1.8左右，其他参数在相同的条件下作出预测，得到下图。从图上可以看出，总和生育率为2.1时，人口总数会逐步上升，到2039年左右达到峰值15.5亿，之后会稳定一段时间一直到2055年以后。总和生育率为1.8时，人口总数也会逐步增加，到2020年左右达到峰值14.4亿，之后会稳定一段时间，到2039年左右开始下降，到2055年降到大约12.5亿。在给定性别比率和死亡率的条件下，考率不同的总和生育

16、率对老年人口发展的影响，得到下图：对比上图的四条曲线，当总和生育率在1.8左右时老年人口将会相对平稳。在相同条件下，考率不同的总和生育率对总人口发展的影响，得到下图：. 对比上图的四条曲线，当总和生育率在1.8左右时总人口将会相对平稳。故选定总和生育率1.8为最终模型中的参数。性别比 用05年数据为基准，死亡率给定，总和生育率取1.8时，性别比分别取1.07和1.03来预测未来50年总人口数的发展，得到下图：比较发现，性别比取1.03时，总人口的变化更加平稳。所以取定性别比为1.03。从上图中可以看出，总人口是先持续上升，至2022年左右达到峰值14.4亿。2022年2039年总人口数有所缓和

17、。2040年总人数开始明显下降，至2055年总人口数大约将降至13.3亿。长期的预测：利用以上所确立的模型作长期预测（100年）。可以清晰地看到2055以后，总和生育率为2.1时，总人口持续上升。总和生育率为1.8时，总人口将继续下降，往后推迟缓慢下降。按照更替水平，总和生育率为2.1时，应该是处于稳定状态，但由于我国的人口老龄化在逐渐加重，而生育率在更替水平时，总人口就会上升，如果生育率一直处于1.8，总人口就会逐步下降。因此，长期的预测应该采取中短期以前总和生育率靠近2.1，中短期之后靠近1.8，这样容易达到稳定。七模型的评价与推广模型的优点：1. 本模型能准确地预测未来各年龄段的人口数量

18、情况和变化趋势，对国家的人口调整工作有很大的帮助，推广容易和实用性强。2. 模型思路清晰，容易理解和灵活应用。3. 本模型作了三个重要的参数设定，逐步优化模型，便于解释一些人口由于受到诸多复杂因素影响而出现波动的状况，较为准确地作出人口预报。 模型的缺点： 未考虑战争、重大灾害和大量人口迁移的情况，若要考虑这些因素，需对模型加以调整。推广： 本论文的模型先固定参数，再采用参数变化解决问题的办法，在其他领域也有许多类似的方法。如：工程预算，市场评估，国民生产总值预算等，尤其对预测和估算方面的问题更是合适。参考文献1 姜启源，数学模型（第三版），北京：高等教育出版社，20032 吴建国，数学建模案

19、例精编，北京：中国水利水电出版社，20053 傅鹂，数学实验，北京：科学出版社，20004 中华人民共和国国家统计局, 2007-9-215 第五次人口普查公报， 2007-9-216 中国人口， 2007-9-227 中国网， 2007-9-22 附录程序1： 计算各类各年各岁人口数。如0年0岁城市男人数，l1是基准年男性各年龄人口（0到89岁），l2是基准年女性各年龄人口（0到89岁），F是各岁妇女生育率（15到49岁），. s是出生的性别比，D1是男性各年龄人口死亡率（0到89岁），D2是女性各年龄人口死亡率（0到89岁）。function L1,L2=fL12(l1,l2,F,s,D1

20、,D2) L1(:,1)=l1; %用矩阵L1表示各年各年龄的男性人口。同一列，同一年；同一行，同一岁L2(:,1)=l2; %用矩阵L2表示各年各年龄的女性人口。同一列，同一年；同一行，同一岁N =100; %预测 N 年for j=2:N+1 %一年一年地预测 for i=2:89+1 %一岁到89岁 L1(i,j)=L1(i-1,j-1)*(1-D1(i-1); %男性i-1岁的人口是上一年i-2岁的人口死亡后剩余的 L2(i,j)=L2(i-1,j-1)*(1-D2(i-1); %女性i-1岁的人口是上一年i-2岁的人口死亡后剩余的EndL1(1,j)=sum(L2(15:49,j).

21、*F)*s(j)/(1+s(j); %0岁的人口是当年15到49岁生育的总和*男性比率L2(1,j)=sum(L2(15:49,j).*F)/(1+s(j); %0岁的人口是当年15到49岁生育的总和*女性比率End程序2： 计算各年各类人口增长率,N指预测N年。function nt=fnt(fx,Hxt1,Hxt2,D1,D2,N) lao=1/1000; %把大于等于90岁的人口假设为总人口的1千分之一dlao=0.2; %假设大于等于90岁的人口死亡率为dlaoforj=1:N+1 nt(j)=sum(fx.*Hxt2(15:49,j)-sum(D1(1:90).*Hxt1(:,j)-

22、sum(D2(1:90).*Hxt2(:,j)-lao*dlao;end程序3： 预测50或100年各类各年各岁的人口数。用05年为基准，05年的生育率的2.1/1.334或1.8/1.334倍，性别比都取107或103，死亡率都取历年最小值。N=100; sex0=1.07; %设定性别比 sex=linspace(sex0,sex0,N+1);zhs=2.1; %设定总和生育率cL1,cL2=fL12(cnan05,cnv05,cs5*0.001*zhs/1.334,sex,cdn*0.001,cdnv*0.001);tL1,tL2=fL12(tnan05,tnv05,ts5*0.001*

23、zhs/1.334,sex,tdn*0.001,tdnv*0.001);vL1,vL2=fL12(vnan05,vnv05,vs5*0.001*zhs/1.334,sex,vdn*0.001,vdnv*0.001);L05g=cL1+cL2+tL1+tL2+vL1+vL2; %基准年开始N+1年的0到89岁人口L05g1=cL1+tL1+vL1; %基准年开始N+1年的0到89岁男性人口 L05g2=cL2+tL2+vL2; %基准年开始N+1年的0到89岁女性人口程序4： 计算各年各岁人口所占比例。(如0年0岁的城市男占0年城市男女总人口的比例)，N指预测N年。function Hxt1,H

24、xt2=fhxt(L1,L2,N) lao=1/1000; %把大于等于90岁的人口假设为总人口的1千分之一for j=1:N+1 for i=1:90 Hxt1(i,j)=L1(i,j)/(sum(L1(:,j)+sum(L2(:,j)*(1+lao); %j-1年i-1岁男性除以j-1年男女总人口 Hxt2(i,j)=L2(i,j)/(sum(L1(:,j)+sum(L2(:,j)*(1+lao); %j-1年i-1岁女性除以j-1年男女总人口 endend附表一： 用05年为基准，05年的生育率的1.8/1.334倍，性别比都取1.03，死亡率都取历年最小值，预测100年的结果2005年

25、2006年2007年2008年2009年2010年2011年老年人1.7E+081.76E+081.84E+081.91E+082E+082.09E+082.18E+08总人数1.31E+091.32E+091.32E+091.33E+091.34E+091.35E+091.36E+092012年2013年2014年2015年2016年2017年2018年老年人2.29E+082.4E+082.52E+082.63E+082.73E+082.84E+082.92E+08总人数1.37E+091.38E+091.39E+091.4E+091.41E+091.42E+091.42E+092019年

26、2020年2021年2022年2023年2024年2025年老年人2.96E+083E+083.01E+083.11E+083.27E+083.39E+083.51E+08总人数1.43E+091.43E+091.44E+091.44E+091.44E+091.44E+091.44E+092026年2027年2028年2029年2030年2031年2032年老年人3.63E+083.71E+083.84E+083.94E+084.06E+084.15E+084.22E+08总人数1.44E+091.44E+091.44E+091.44E+091.44E+091.44E+091.44E+0920

27、33年2034年2035年2036年2037年2038年2039年老年人4.28E+084.33E+084.35E+084.36E+084.35E+084.34E+084.32E+08总人数1.44E+091.43E+091.43E+091.43E+091.43E+091.43E+091.43E+092040年2041年2042年2043年2044年2045年2046年老年人4.29E+084.26E+084.24E+084.19E+084.13E+084.07E+084.04E+08总人数1.42E+091.42E+091.41E+091.41E+091.4E+091.4E+091.39E+

28、092047年2048年2049年2050年2051年2052年2053年老年人4.04E+084.05E+084.09E+084.15E+084.19E+084.19E+084.16E+08总人数1.38E+091.38E+091.37E+091.36E+091.36E+091.35E+091.34E+092054年2055年2056年2057年2058年2059年2060年老年人4.12E+084.09E+084.05E+084.01E+083.96E+083.89E+083.83E+08总人数1.33E+091.33E+091.32E+091.31E+091.31E+091.3E+091.29E+092061年2062年2063年2064年2065年