中国人口增长预测
摘要
中国乃泱泱人口大国，人口规模是城市规划和土地利用总体规划中一项重要的控制性指标，预测人口模型的合理性，不仅影响到未来地区经济和社会发展，而且会影响到地区生态环境可持续发展。 因此，建立合理的模型，准确地预测未来人口的发展趋势，制定合理的人口规划和人口布局方案具有重大的理论意义和现实意义。
对此，本文通过建立适当的模型，预测出了短期和中长期（到2050年）中国人口的变化趋势和走向，并给出了在这段时间内人口结构的具体预测数据和曲线走向，包括总人口数、年龄结构、出生率和死亡率等。
在此模型中，为精确预测，我们用到了人口密度、生育率、死亡率、人口总数以及迁出率等影响人口的因数，并将我国人口整合为一个由城市男性、城市女性、城镇男性、城镇女性、乡村男性、乡村女性组成的1x6的矩阵。同时用人口密度、生育率、死亡率及迁出率作为参数并结合人口发展偏微分方程，再通过完善和改进，建立了一个一阶偏微分方程的模型。最后以此模型作为基础，进行人口数据的相关预测。
对于求解一阶偏微分方程模型中的相关参数，我们首先用MATLAB和EXCEL等软件对题目所给的2001年到2005年的数据进行处理和适当筛选。在求解生育率时，通过用MATLAB的曲线拟合工具箱，经处理和比较，最后选取了高斯分布作为建立求解生育率的模型，合理而精确；在求解死亡率时，用EXCEL软件作出了各年各年龄段的折线图，为使模型更加精确和实际，通过观察我们把年龄分为三段，用分段函数概念，分别求解这三段的死亡率，其中在求解第三段（衰老期）时，使用了指数函数模型；在求解迁出率时，考虑城镇化进程对我国人口分布的影响，我们对复杂因素适当简化，建立了理想化的迁出率子模型。
在所有参数求到之后，剩下的就是求解模型中的一阶偏微分方程，对此，我们对数据进行离散化处理，化偏微分方程为差分方程，运用计算机模拟的方法预测出相关数据，在这个过程中我们用MATLAB编程实现。把通过该模型预测出的结果与国内外专家、学者预测的数据作对比，其在定性趋势与定量分析上的结果基本一致，模型可信度较高。
最后，在模型改进中迁出率用阻滞增长模型优化。利用MATLAB中Curve Fitting Tool拟合出城市化率增长模型，进而用相邻城市化率相减即得迁移率矩阵。此方法求出的迁移率更为精确，使预测结果更为准确。
关键词： 人口发展偏微分方程 MATLAB曲线拟合 阻滞增长模型 方程离散化求解 高斯分布
一、问题重述
中国是一个人口大国，人口问题始终是制约我国发展的关键因素之一。根据已有数据，运用数学建模的方法，对中国人口做出分析和预测是一个重要问题。
近年来中国的人口发展出现了一些新的特点，例如，老龄化进程加速、出生人口性别比持续升高，以及乡村人口城镇化等因素，这些都影响着中国人口的增长。2007年初发布的《国家人口发展战略研究报告》(附录1) 还做出了进一步的分析。
关于中国人口问题已有多方面的研究，并积累了大量数据资料。附录2就是从《中国人口统计年鉴》上收集到的部分数据。
试从中国的实际情况和人口增长的上述特点出发，参考附录2中的相关数据（也可以搜索相关文献和补充新的数据），建立中国人口增长的数学模型，并由此对中国人口增长的中短期和长期趋势做出预测；特别要指出你们模型中的优点与不足之处。
二、问题分析
本题是一个人口发展预测的问题。人口发展与一般种群增长一样，是由自然增长率决定的。然而，人类个体是一种社会的个体，所以人口发展有自己的特点。想到人口的迁移，性别比例，城镇化等。同时，人口发展受政策的影响，例如计划生育；也要受到人们意识的影响，像生育意识等。但是从社会层面上看，生育意识在整个社会上体现为妇女的生育模式，进而可以特别地去考虑。
首先，数据的处理。在经过EXCEL分析和验证后，适当修正题中的个别有误数据后，利用有效数据进行建模求解，在此过程中，我们提取出死亡率、生育率等感念，且把人的一生按年龄分为青年期、衰老期等阶段。这是求解人口增长模型的必要过程和方法。
其次，模型建立。和一般的预测模型一样，本模型也是个预测模型，所以考虑到用题目所给的五年的信息，来推测今后几十年的人口的总数和结构情况。对此，我们选用差分方程模型和数据参数拟合等方法。同时，将死亡率与出生率分开分别计算和拟合，通过五年的实际数据拟合出相应函数的参数，再利用此函数进行评估和预测。
最后，利用已有信息以及上述所求出的对应函数和方程，对中短期与长期进行估计和预测，进而得出人口结构、人口比例、人口数量等一系列的相关数据。
三、符号说明
xm： 表示最大的年龄；
i=1,2,3,4,5,6 其中 1表示市男性，2表示市女性， 3表示镇男性，4表示镇女性， 5表示乡男性，6表示乡女性； A ：表示婴儿性别比例矩阵； * ：表示点乘；
Pi(x,t) ：表示t时刻年龄为x的人口数量；
biri(x,t) ：表示t时刻年龄为x的出生率；
deai?x,t? ：表示t时刻年龄为x的死亡率；
k?t? ：表示t时刻婴儿的死亡率；
itrai(x,t) ：表示t时刻年龄为x的人口迁出率；
四、模型假设
1. 假设国内社会环境稳定，无异常大量死亡或出生情况发生，人口比例，人口
总数不会出现突变状况； 2. 假设只存在乡向城镇迁出，不存在其他迁移方式，且不同年龄段迁移率相同； 3. 假设不考虑国家之间的迁入与迁出，把中国内部看为一个封闭的模型； 4. 对于90岁以上的人都按照90岁处理； 5. 假设只存在乡向城镇迁出，不存在其他迁移方式，且不同年龄段迁出率相同，
按照0.6%均匀增长。
五、模型建立
1.首先为了求解方便，建立了人口数量的向量，死亡率向量，出生率向量，迁移率向量和新生婴儿死亡率向量（以下的符号量见符号说明）：
TP?x,t?=(P1?x,t?，P2?x,t?，P3?x,t?，P4?x,t?，P5?x,t?，P6?x,t?）
dea(x,t)=( dea1(x,t), dea2(x,t), dea3(x,t), dea4(x,t), dea5(x,t), dea6(x,t)) T
bir(x,t)=(bir1(x,t),bir2(x,t),bir3(x,t),bir4(x,t),bir5(x,t),bir6(x,t) T
tra(x,t)=( tra1(x,t), tra2(x,t), tra3(x,t), tra4(x,t), tra5(x,t), tra6(x,t))T k(t)=( k1(t), k2(t), k3(t), k4(t), k5(t), k6(t)) T
2.由已有的人口模型可知：
P?x??x,t??t??P?x,t??t?P?x,t??t??P?x,t?+ =—[tra(x,t)+dea(x,t)]*P(x,t))
?x?t使?x ?0，?t ?0，则可得：
?P?P??—[tra(x,t)+dea(x,t)]*P(x,t)) ?t?x又由1994到2005年的城镇乡点的男女出生比例求平均值可得婴儿出生比例的矩阵A，即：
000??00.5280??00.4720000???0000.54100?A???
0000.45900???000000.545????000000.455???最后结合条件可得人口微分模型:
??P?P??t??x??[tra(x,t)?dea(x,t)]*P(x,t)?xm??P(0,t)?(1?k)*A?bir(x,t)*P(x,t)dx （1）
0??P(x,0)?P0(x)?P(xm,t)?0?将其转化成差分方程可得：
P?x??x,t??t??P?x,t??t?P?x,t??t??P?x,t????[tra(x,t)?dea(x,t)]*P(x,t)?x?tP(0,t)?(1?k)*A?bir(x,t)*P(x,t)?x0xm
P(x,0)?P0(x)P(xm,t)?0为了求解上述方程，则只需已知出生率，死亡率，迁移率的值，在模型求解中给出了其具体求解的方法。
六、模型求解
一、
死亡率的求解
由不同阶段导致死亡的因素不同，我们将其分成了3个阶段——婴儿期，青年期和老年期，然后对其分别求解。
1．衰老期死亡率求解
我们通过matlab曲线拟合工具箱进行回归分析和数据拟合，经过试验，可以得到指数函数是最理想的，用matlab分别作出城男、城女，镇男、镇女、乡男、乡女衰老期的死亡率散点图及拟合曲线如下：
城男性人口死亡率拟合曲线图 城女性人口死亡率拟合曲图
镇男性人口死亡率拟合曲线图 镇女性人口死亡率拟合曲图