美国人口增长预测模型

20162016 年数学建模论文年数学建模论文 第第 一一 套套 论文题目 论文题目 人口增长模型的确定人口增长模型的确定 组组 别 别 第第 35 组组 姓姓 名 名 耿晨耿晨 闫思娜闫思娜 王强王强 提交日期 提交日期 2016 年年 7 月月 4 日日 题目 美国人口增长预测模型题目 美国人口增长预测模型 摘摘 要要 本文根据近两个世纪美国每十年一次的人口统计数据 建立了指数增长模型 即 Malthus 模型 并通过 1790 1890 年的数据验证了它的准确性 但是 随着时间的推移 拟合函数与统计数据误差逐渐增大 所以 又建立了阻滞增长模型 即 Logistic 模型 这个模型的拟合函数与统计数据误差较小 并用该模型对美国未来几年的人口做出了 预测 总体来说 阻滞增长模型在预测准确度方面要明显优于原始的马尔萨斯人口指 数增长模型 关键词 关键词 指数增长模型 阻滞增长模型 人口预测 一 问题重述一 问题重述 1790 1980 年间美国每隔 10 年的人口记录如下表所示 表表 1 1 人口记录表 人口记录表 年份1790180018101820183018401850186018701880 人口 106 3 95 37 29 612 917 123 231 438 650 2 年份1890190019101920193019401950196019701980 人口 106 62 976 092 0106 5123 2131 7150 7179 3204 0226 5 1 试用以上数据建立马尔萨斯 Malthus 人口指数增长模型 并对接下来的每隔十 年预测五次人口数量 并查阅实际数据进行比对分析 2 如果数据不相符 再对以上模型进行改进 寻找更为合适的模型进行预测 并 对两次预测结果进行对比分析 3 查阅资料找出中国人口与表 1 同时期的人口数量 用以上建立的两个模型进行 人口预测与分析 二 问题分析二 问题分析 影响人口增长的因素很多 其中最主要的两个因素是出生率和死亡率 出生率受 到婴儿死亡率 对避孕的态度及措施效果 对堕胎的态度 怀孕期间的健康护理等因 素的影响 死亡率则受到卫生设施与公共卫生状况 战争 污染 医疗水平 饮食习 惯 心理压力和焦虑等因素的影响 此外 影响人口在一个地区增长的因素还有迁入 和迁出 生存空间的限制 水和食物 疾病等 在这些因素中 有些是常态的或者有 规律的 这些因素对人口的增长是恒定的 而有些因素是随机的 对人口的增长是没 有规律的 因此 当大范围 长时期研究人口增长问题时 对人口增长产生影响的随 机因素就不在考虑了 建立该模型的目的是要能通过模型预测美国后来每十年的人口数具体变化 并与 实际的数据进行对比 看误差的大小 在此基础上利用改进的模型对美国人口同时期 数量进行预测 并进行总结分析 三 问题假设三 问题假设 人口指数增长模型中采用以下基本假设 1 单位时间的人口总量增长与当时的人口呈正比 比例常数为 k 2 假设 t 时刻的人口为 N t 因为人口数一般是很大的 所以将 N t 近似地视 为连续 可微的函数 记初始时刻 t 0 的人口数为 N0 新生人口数百分率为 a 死 亡的百分率为 b 那么 经过 t 时间后 人口数量为 N t t 就是原来人口数量加 上 t 时间内新生人口数减去死亡人口数 四 变量说明四 变量说明 t0 数据的起始时间 即 1790 年 t 时间变量 r 人口固有增长率 N0 当时间 t 1790 时的人口数量 即 3 9 10 6 人 N t t 时刻人口数量 Nm 最大人口容量 五 模型建立五 模型建立 在在 Malthus 的人口指数增长模型中 根据假设我们可以得到 1 N ttN taN ttbN tt 上式进行变形 bN N aN ttkN t t 其中 tNN t N t 可见在一段时间内 人口的变化和人口的数量成正比 用瞬时变化率逼近平均变 化率得 2 00 dN kN dt N tN 式中 N0 初始时刻的人口数 N 人口数 用分离变量法得到上述方程的解为 3 0 0 kt kt N e N t e 即方程的解为 4 0 0 k t t NN e 六 模型求解六 模型求解 利用 MATLAB 数学工具 对数据中前十一年 1790 1890 的人口数拟合 4 式 现 对问题进行变形以便处理 对 4 式进行取对数 得 5 00 lnlnNk ttN 取 在 MATLAB 中输入程序 程序 1 0 0 t 可以得到 k 0 2808 lnN0 1 4107 即 N0 4 0988 所以 1790 年到 1890 年的拟合函数为线性函数 所以 ln0 28081 4107Nt 6 0 2808 4 0988 t Ne 绘制离散点和拟合函数的图像进行对比 在 MATLAB 中输入程序 程序 2 得图 1 如下 图图 1 1790 18901 1790 1890 年美国人口拟合曲线图年美国人口拟合曲线图 比较拟合的曲线图和散点图 图 1 可以发现与 19 世纪的人口增长情况相吻合 类似地 可以求出 1790 年至 1990 年的拟合函数为 7 0 2142 5 5918 t Ne 拟合过程如下 在 MATLAB 中输入程序见 程序 3 可以得到 k 0 2142 lnN0 1 7213即N0 5 5918 所以 1970 年到 1990 年的拟合函数为线性函数 所以 0 21421 7213yx 0 2142 5 5918 t Ne 绘制离散点和拟合函数的图像进行对比 在 MATLAB 中输入程序 程序 4 得图 2 如下 图图 2 1790 19902 1790 1990 年美国人口拟合曲线图年美国人口拟合曲线图 此时我们可以观察到虽然前期的数据能较好的吻合 但是随着时间的推移并不能 很好的反应后期的人口真实情况 函数估计值与统计数据的误差越来越大 这说明 用指数增长模型预测短期人口的数量可以得到较好的结果 但是从长期来看 任何地 区的人口数量都不可能无限制地增长 因此 指数增长模型不适合预测长时期人口的 增长情况 Malthus 模型中 我们只考虑了出生率和死亡率对人口的影响 而忽略了其他因素 如自然资源 生存环境等对人口的影响 然而这些因素对人口增长起着阻滞作用 并 且随着人口数量的增加 阻滞作用也会增大 因此需要将增长率 k 看作是人口数量的 函数 丹麦生物学家 Pierre Francois Verhulst 在指数增长模型的基础上建立了改进的 Malthus 模型 即 Logistic 模型 一般来说人口的增长率是变化的 当人口较少时 增长速度较快 增长率较大 当增加至一定的数量时 增长速度必然会减慢 增长率开始减小 因此增长率 K 应该 视为人口数量的函数 在修正的人口阻滞增长模型中 有以下假设 1 假设人口增长率 k t 是 t 时人口 N t 的函数 随着人口的增加 自然资源 环境条件等对人口增长的阻滞作用越来越明显 k N 应是 N 的减函数 一个简单的假设是 k N 是 N 的线性函数 k N r sN 其中 s 0 k 0 式中 r 称为 固有增长率 表示人口很少时的增长率 其中 Nm称为最大人口容量 m k s N 2 考虑自然资源和环境因素所能容纳的最大人口数量 Nm 当 N Nm时增长率 为 0 环境饱和 即 k Nm 0 考虑环境因素的限制因素 对于阻滞增长模型 在假设的前提下我们可以得到 8 1 m N k xr N 式中 人口增长率 k x Nm 最大人口容量 在上式的假设前提下 指数模型可以修改为 9 00 1 m dNN rN dtN N tN 上述方程称为阻滞增长模型 Logistic 模型 上述方程中我们可以看到 方程右 面的因子中 当 N 增大时 相应的将会减小 即人口增长由这两个因子控制 1 m N N 比较符合实际情况 方程 9 的解为 10 0 t 0 11 m r t m N N N e N 在 9 式中 当时间 t 无限大时 N 接近限制 Nm 这表明随时间推移人口数最 大达到饱和值 Nm 进一步考虑增长率变化问题 由 9 式可得 2 2 2 1 m d NdN N r dt dtN 令 0 则 这表明当人口数量达到最大人口数一半时 人口增长率 2 2 d N dt 2 m N N 达到最大 此后开始不断减小 不同时期的增长率见图三所示 见程序 5 图图 3 3 增长率随时间变化的关系图 增长率随时间变化的关系图 用 MATLAB 对函数式 10 进行模拟可得图 4 程序 6 图图 4 4 LogisticLogistic 模型曲线图与美国实际人口对比模型曲线图与美国实际人口对比 采用阻滞增长模型对美国人口进行预测的结果如图 4 所示 我们能够看到 阻滞 增长模型很好的处理了指数增长模型无限增长的弊端 对长远的人口情况估计更符合 实际情况 但是与统计值之相比较仍存在一定的误差 并且预测值小于统计值 这是 因为我们忽略了社会因素还有其他导致地区人口减少的因素 事实上 对于复杂生命 史和个体生长期较长的高等动植物简单的统计不足以反映出全部的问题 七 结果分析七 结果分析 通过对美国人口统计模型的建立和修改 最终采用更为实际的人口增长阻滞模型 并对中国同期人口进行统计 采用两种模型分别对中国人口进行预测 具体过程如下 首先调查统计同时期中国人口 人口数分布情况如表 2 所示 表表 2 2 中国人口调查记录表 中国人口调查记录表 1790 19801790 1980 年份1790180018101820183018401850186018701880 人口 106 301 5366 0340 0383 1409 0418 9430 0377 0357 7368 0 年份1890190019101920193019401950196019701980 人口 106 380 0400 0427 7472 0489 0518 8546 8667 1818 3981 2 表表 3 3 中国人口调查记录表 中国人口调查记录表 1991 20101991 2010 年份1991199219931994199519961997199819992000 人口 106 1158 21171 71185 21198 51211 21223 91236 31247 61257 91267 4 年份2001200220032004200520062007200820092010 人口 106 1276 31284 51292 31299 91307 61314 51321 31328 01334 51340 9 采用指数增长的马尔萨斯模型 首先对 1790 到 1980 年的人口数据作为依据 得 到相应的模型参数 在 MATLAB 中输入程序 程序 7 可以得到 k 0 0425 lnN0 5 7074 即 N0 301 0872 所以 1790 年到 1890 年的拟合函数为线性函数 所ln0 0425301 0872Nt 以 11 0 0425 301 0872 t Ne 绘制离散点和拟合函数的图像进行对比 在 MATLAB 中输入程序 程序 8 图图 5 5 马尔萨斯模型曲线图与中国实际人口对比 马尔萨斯模型曲线图与中国实际人口对比 由图可见 马尔萨斯不能很好的反应长期时间人口发展的情况 正如美国的情况 一样 因此人口的预测需要改进的马尔萨斯模型更为准确 采用阻滞增长的 Logistic 模型对中国人口增长情况的预测时 首先分析不同时期的 增长率变化情况 从而得出人口最大值 Nm 其实现程序如程序 9 在此基础上对 1790 2010 年人口进行模型与统计数据的对比 在 MATLAB 中输入 程序 10 得到口增长率变化图 6 图图 6 6 1790 20101790 2010 人口增长率变化情况人口增长率变化情况 针对增长率的变化和历史的角度可以发现 中国需要预测人口变化情况 必须得 从近二十年的人口变化来寻求规律 则需要重新对中国人口数据进行模拟及预测 对 1991 2010 年马尔萨斯拟合 在 MATLAB 中实现程序 11 得到图 7 图图 7 7 1991 20101991 2010 年改进的马尔萨斯人口增长率变化情况年改进的马尔萨斯人口增长率变化情况 参数确定后 我们采用该模型来预测中国未来几年的人口数情况 在 MATLAB 中键入 程序 12 得到图 8 图图 8 8 LogisticLogistic 模型曲线图与中国实际人口对比模型曲线图与中国实际人口对比 该模型最终很好的预测了人口的发展状况 实现了人口预测的目标 就统计值和 预测值之间的差正如前面所述 我们了解到高等动植物的生存期限会受多方面的因素 影响 因此我们可以说 在有条件的选则处理数据的时候尽可能选择对未来有相似生 存条件或者生活质量的数据 以保证在数据处理的后期能得到比较理想的模型 因此 我们可以说 本次模型建立是成功地实现了中国人口预测的母的 八 参考文献八 参考文献 1 张志涌 杨祖樱 MATLAB 教程 M 北京 北京航空航天大学出版社 2011 2 夏鸿鸣 魏艳华 王丙参 数学建模 成都 西南交通大学出版社 2014 3 周品 赵新芬 数学建模 北京 国防工业出版社 2009 4 王中鲜 MATLAB 建模与仿真应用 机械工业出版社 2010 九 附录九 附录 程序程序 1 1 x 0 1 10 y 3 9 5 3 7 2 9 6 12 9 17 1 23 2 31 4 38 6 50 2 62 9 y log y k polyfit x y 1 程序程序 2 2 x 0 1 10 y 3 9 5 3 7 2 9 6 12 9 17 1 23 2 31 4 38 6 50 2 62 9 y1 4 0988 exp 0 2808 x plot x y x y1 程序程序 3 x 0 1 19 y 3 9 5 3 7 2 9 6 12 9 17 1 23 2 31 4 38 6 50 2 62 9 76 0 92 0 106 5 123 2 131 7 150 7 179 3 204 0 226 5 y log y k polyfit x y 1 程序程序 4 4 x 0 1 19 y 3 9 5 3 7 2 9 6 12 9 17 1 23 2 31 4 38 6 50 2 62 9 76 0 92 0 106 5 123 2 131 7 150 7 179 3 204 0 226 5 y1 5 5918 exp 0 2142 x plot x y x y1 程序程序 5 5 A zeros 0 19 计算 1790 1980 年的人口增长率 A 为人口增长率的矩阵 for i 1 19 A i y i 1 y i y i end figure 2 A1 A 0 plot x A1 or 程序程序 6 6 q polyfit y A1 1 q 0 3389 0 0014 ym 242 0714 for i 21 25 A2 i 0 0014 i 0 3389 end A2 for i 1 20 y3 i 242 0714 1 61 0696 exp 0 3389 x i end y3 figure 3 plot x y ob x y3 or title 改进马尔萨斯 程序程序 7 7 x 0 1 19 y 301 5 366 0 340 0 383 1 409 0 418 9 430 0 377 0 357 7 368 0 380 0 400 0 427 7 472 0 489 0 518 8 546 8 667 1 818 3 981 2 y log y k polyfit x y 1 程序程序 8 8 x 0 1 19 y 301 5 366