人口预测MATLAB.docxVIP

天津商业大学宝德学院Tianjin University of Commerce Boustead College1982-2023中国人口分析与预测作 者: 李 紫 琳 何 龄 童 汪 晨 1982-2023中国人口分析与预测摘 要：本文以Logistic人口阻滞增长模型为基础建立了我国人口增长预测模型，并就1982-2012年实际人口与预测人口进行对比分析，最后预测了2013-2023我国全国总人口，从而为我国人口控制与管理提供一定的依据。关键词：Logistic模型；最小二乘法；人口增长；MATLAB软件1问题分析中国是一个人口大国,人口问题始终是制约我国经济发展的关键因素之一。而今全面建设小康社会时期是我国社会快速转型期,人口发展面临着前所未有的复杂局面,人口安全面临的风险依然存在。因此,如何准确地判断我国人口在未来若干年的发展趋势就显得非常重要。另外，我国人口发展经历了多个阶段,特别是自1979年以后，实施了计划生育政策使得中国人口的增长进入一个相对平稳的时期，所以本文选取1982-2012年的全国人口总人数作为依据，对中国未来的人口发展趋势作了一定的预测。2模型假设（1）自然资源和环境因素对人口的增长期阻滞作用，人口规模增大时，人口增长率降低；（2）自然资源和环境所容许的最大人口为常数Nm，并且人口总数的净相对增长率是是人口总数的线性递减函数，设为r（N）=r（1-N/Nm）,表示人口相对增长率随N（t）的增加而减少，其中r为固有增长率。当N（t）→Nm时，人口净相对增长率r（N）趋于零。3符号说明Nm：人口最大容量r：固有增长率t：年份N(t)：t时刻的人口总数N0：人口初始值T0：初始年份4模型建立与求解由上述假设，令 (1.15)方程（1.15）右端因子rN(t)体现人口自身的增长趋势，因子 则体现自然资源和环境对人口增长的阻滞作用。显然N（t）越大，前一因子越大，后一因子越小，表明人口的增长是两个因子共同作用的结果。用分离变量法可解得方程（1.15）的解为（1.16）5模型分析与讨论 对于Logistic模型的解（1.16），有如下结论：当r0时，随着t→+∞，必有N（t）→0；当r=0时，模型的解为常数N（t）=N0，称此常数解为平衡解；当r0时，不论N0取何值，均有N（t）→Nm，即此时任何解都趋于平衡解，此解是稳定的。称0N0Nm时模型的解曲线为Logistic曲线，如图1.16所示。图1.16 NNmNm/2Ot为了考察何时人口增长速度最快，对求导，得到 不难看出，当0N时，单调增；当N时，单调减；当N=时，取最大值，此时人口增长速最大，即人口增长速度达到拐点。表1.1年份实际人口/亿198210.1654198310.3008198410.4357198510.5851198610.7507198710.9300198811.1026198911.2704199011.4333199111.5823199211.7171199311.8517199411.9850199512.1121199612.2389199712.3626199812.4761199912.5786200012.6743200112.7627200212.8453200312.9227200412.9988200513.0756200613.1448200713.2129200813.2802200913.3450201013.4091201113.4735201213.5404数据来源：国家统计局注：1981年及以前人口数据为户籍统计数；1982、1990、2000、2010年数据为当年人口普查数据推算数；其余年份数据为年度人口抽样调查推算数据。6 模型的参数估计与检验采用表（1.1）中的数据来模拟我国人口数，用Logistic模型作预测，需要估计（1.16）中的三个参数Nm，r和N0 ，用MATLAB工具箱进行拟合如图1.18，图1.18计算出参数Nm=14.43 ， r =0.0622 ，N0 =10.08 ，代入（1.16）式得 =（1.17）令t=1982,1983，…，2023，分别代入（1.17）式即可算得2013—2023年的预测人口数。如表1.19和图1.20，表1,21和图1.22。表1.19年份实际人口/亿预测人口/亿绝对误差相对误差198210.165410.08000.08540. 84%198310.300810.26660.03420. 33%198410.435710.4488-0.01270. 12%198510.585110.6252-0.04010. 38%198610.750710.7969-0.04620.43%198710.930010.9634-0.03