中国人口增长预测数学建模_1.docx

中国人口增长预测数学建模高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开场后参赛队员不能以任何方式包括电话、电子邮件、网上咨询等与队外的任何人包括指导老师研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违背竞赛规则的,假如引用别人的成果或其他公开的资料包括网上查到的资料，必须根据规定的参考文献的表述方式在正文引用途和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违背竞赛规则的行为，我们将遭到严肃处理。我们参赛选择的题号是从A/B/C/D中选择一项填写：我们的参赛报名号为假如赛区设置报名号的话：所属学校请填写完好的全名：参赛队员(打印并签名)：1.2.3.指导老师或指导老师组负责人(打印并签名)：日期：年月日赛区评阅编号由赛区组委会评阅前进行编号：中国人口增长预测摘要中国乃泱泱人口大国，人口规模是城市规划和土地利用总体规划中一项重要的控制性指标，预测人口模型的合理性，不仅影响到将来地区经济和社会发展，而且会影响到地区生态环境可持续发展。因而，建立合理的模型，准确地预测将来人口的发展趋势，制定合理的人口规划和人口布局方案具有重大的理论意义和现实意义。对此，本文通过建立适当的模型，预测出了短期和中长期到2050年中国人口的变化趋势和走向，并给出了在这段时间内人口构造的详细预测数据和曲线走向，包括总人口数、年龄构造、出生率和死亡率等。在此模型中，为准确预测，我们用到了人口密度、生育率、死亡率、人口总数以及迁出率等影响人口的因数，并将我国人口整合为一个由城市男性、城市女性、城镇男性、城镇女性、乡村男性、乡村女性组成的1x6的矩阵。同时用人口密度、生育率、死亡率及迁出率作为参数并结合人口发展偏微分方程，再通过完善和改良，建立了一个一阶偏微分方程的模型。最后以此模型作为基础，进行人口数据的相关预测。对于求解一阶偏微分方程模型中的相关参数，我们首先用MATLAB和EXCEL等软件对题目所给的2001年到2005年的数据进行处理和适当挑选。在求解生育率时，通过用MATLAB的曲线拟合工具箱，经处理和比拟，最后选取了高斯分布作为建立求解生育率的模型，合理而准确；在求解死亡率时，用EXCEL软件作出了各年各年龄段的折线图，为使模型愈加准确和实际，通过观察我们把年龄分为三段，用分段函数概念，分别求解这三段的死亡率，其中在求解第三段衰老期时，使用了指数函数模型；在求解迁出率时，考虑城镇化进程对我国人口分布的影响，我们对复杂因素适当简化，建立了理想化的迁出率子模型。在所有参数求到之后，剩下的就是求解模型中的一阶偏微分方程，对此，我们对数据进行离散化处理，化偏微分方程为差分方程，运用计算机模拟的方法预测出相关数据，在这个经过中我们用MATLAB编程实现。把通过该模型预测出的结果与国内外专家、学者预测的数据作比照，其在定性趋势与定量分析上的结果基本一致，模型可信度较高。最后，在模型改良中迁出率用阻滞增长模型优化。利用MATLAB中CurveFittingTool拟合出城市化率增长模型，进而用相邻城市化率相减即得迁移率矩阵。此方法求出的迁移率更为准确，使预测结果更为准确。关键词：人口发展偏微分方程MATLAB曲线拟合阻滞增长模型方程离散化求解高斯分布一、问题重述中国是一个人口大国，人口问题始终是制约我国发展的关键因素之一。根据已有数据，运用数学建模的方法，对中国人口做出分析和预测是一个重要问题。近年来中国的人口发展出现了一些新的特点，例如，老龄化进程加速、出生人口性别比持续升高，以及乡村人口城镇化等因素，这些都影响着中国人口的增长。2007年初发布的（国家人口发展战略研究报告）(附录1)还做出了进一步的分析。关于中国人口问题已有多方面的研究，并积累了大量数据资料。附录2就是从（中国人口统计年鉴）上采集到的部分数据。试从中国的实际情况和人口增长的上述特点出发，参考附录2中的相关数据可以以搜索相关文献和补充新的数据，建立中国人口增长的数学模型，并由此对中国人口增长的中短期和长期趋势做出预测；十分要指出你们模型中的优点与缺乏之处。二、问题分析此题是一个人口发展预测的问题。人口发展与一般种群增长一样，是由自然增长率决定的。然而，人类个体是一种社会的个体，所以人口发展有本人的特点。想到人口的迁移，性别比例，城镇化等。同时，人口发展受政策的影响，例如计划生育；也要遭到人们意识的影响，像生育意识等。但是从社会层面上看，生育意识在整个社会上体现为妇女的生育形式，进而能够十分地去考虑。首先，数据的处理。在经过EXCEL分析和验证后，适当修正题中的个别有误数据后，利用有效数据进行建模求解，在此经过中，我们提取出死亡率、生育率等感念，且把人的一生按年龄分为青年期、衰老期等阶段。这是求解人口增长模型的必要经过和方法。其次，模型建立。和一般的预测模型一样，本模型也是个预测模型，所以考虑到用题目所给的五年的信息，来揣测今后几十年的人口的总数和构造情况。对此，我们选用差分方程模型和数据参数拟合等方法。同时，将死亡率与出生率分开分别计算和拟合，通过五年的实际数据拟合出相应函数的参数，再利用此函数进行评估和预测。最后，利用已有信息以及上述所求出的对应函数和方程，对中短期与长期进行估计和预测，进而得出人口构造、人口比例、人口数量等一系列的相关数据。三、符号讲明mx：表示最大的年龄；i=1,2,3,4,5,6其中1表示市男性，2表示市女性，3表示镇男性，4表示镇女性，5表示乡男性，6表示乡女性；A：表示婴儿性别比例矩阵；*：表示点乘；iP(x,t)：表示t时刻年龄为x的人口数量；ibir(x,t)：表示t时刻年龄为x的出生率；)(,ideaxt：表示t时刻年龄为x的死亡率；)(itk：表示t时刻婴儿的死亡率；itra(x,t)：表示t时刻年龄为x的人口迁出率；四、模型假设1.假设国内社会环境稳定，无异常大量死亡或出生情况发生，人口比例，人口总数不会出现突变状况；2.假设只存在乡向城镇迁出，不存在其他迁移方式，且不同年龄段迁移率一样；3.假设不考虑国家之间的迁入与迁出，把中国内部看为一个封闭的模型；4.对于90岁以上的人都根据90岁处理；5.假设只存在乡向城镇迁出，不存在其他迁移方式，且不同年龄段迁出率一样，根据0.6%均匀增长。五、模型建立1.首先为了求解方便，建立了人口数量的向量，死亡率向量，出生率向量，迁移率向量和新生婴儿死亡率向量下面的符号量见符号讲明：)(,Pxt=()(1,Pxt，)(2,Pxt，)(3,Pxt，)(4,Pxt，)(5,Pxt，)(6,PxtT(,)deaxt=(1(,)deaxt,2(,)deaxt,3(,)deaxt,4(,)deaxt,5(,)deaxt,6(,)deaxt)T(,)birxt=(1(,)birxt,2(,)birxt,3(,)birxt,4(,)birxt,5(,)birxt,6(,)birxtT(,)traxt=(1(,)traxt,2(,)traxt,3(,)traxt,4(,)traxt,5(,)traxt,6(,)traxt)T()kt=(1()kt,2()kt,3()kt,4()kt,5()kt,6()kt)T2.由已有的人口模型可知：()(),PxxttPxttx+?+?-+?+()(),PxttPxtt+?-?=tra(x,t)+dea(x,t)*P(x,t)使x?0，t?0，则可得：PPtx?+=?tra(x,t)+dea(x,t)*P(x,t)又由1994到2005年的城镇乡点的男女出生比例求平均值可得婴儿出生比例的矩阵A，即：00.528000000.47200000000.541000000.45900000000.545000000.455A?=?最后结合条件可得人口微分模型:00(,)(,)*(,)(0,)(1)*(,)*(,)(,0)()(,)0mxmPPtraxtdeaxtPxttxPtkAbirxtPxtdxPxPxPxt?+=-+?=-?=?=?1将其转化成差分方程可得：