《07年全国数学建模优秀论文》由会员分享，可在线阅读，更多相关《07年全国数学建模优秀论文（15页珍藏版）》请在金锄头文库上搜索。
1、华中科技大学 2007CUMCM 获奖论文 文峰 郭宪 刘飞，指导教师：高亮 1 基于差分方程的人口预测模型 基于差分方程的人口预测模型 摘要摘要 本文针对我国人口增长中出现的新特点，建立了两个符合实际情况的预测模型。 模型一：基于 Logistic 模型，建立了含市、镇、乡人口相互流动关系的微分方程模 型， 求得全国总人口数在短期内将持续增长， 到 2010 年、 2020 年分别为 13.59 亿和 14.44 亿，具有较好的中短期预测效果。 模型二：从年龄转移与总和生育率出发，建立了离散型人口发展模型。 1. 针对性别比例，引入女性比例转移矩阵，利用计算机进行随机模拟，建立起动态 的女性
2、比例转移矩阵； 2. 针对人口迁移，引入人口迁移率矩阵，将迁移率标准化后利用平均迁移率实现了 对迁移人数的预测； 3. 针对死亡率，利用分段加权法估计其随时间的变化，得到了较好的预测结果； 4. 针对老龄化和出生高峰， 将其转化为育龄妇女占总人口的比例， 实现了量化预测。 综合考虑上述因素，利用 MATLAB 编程求解，全国总人口将呈现短期增加、中期 平稳、 长期缓慢下降的趋势。 其中总人口在 2010 年、 2020 年将分别达到 13.63 亿和 14.57 亿人；在 2039 年将达到人口高峰 14.65 亿人，性别比逐步逼近 46.7，69 岁以上人口 将超过 4 亿人，占总人口的 2
3、7.48。 模型一需要的原始数据少，操作简单，适合于中短期预测，但长期预测效果不佳； 模型二综合考虑了各因素，对中短期和长期均有较好的预测效果，但所需数据量大，操 作较为复杂。 关键字关键字 Logistic 模型 人口发展模型 转移矩阵 计算机模拟 迁移率 华中科技大学 2007CUMCM 获奖论文 文峰 郭宪 刘飞，指导教师：高亮 2 1. 问题重述问题重述 人口预测是国家工作中的重点，关系着国家的发展方向和命运。我国是一个人口大 国，人口问题始终是制约我国发展的关键因素之一。根据已有数据，运用数学建模的方 法，对我国人口做出分析和预测是一个重要问题。 近年来我国的人口发展出现了一些新的特
4、点，例如，老龄化进程加速、出生人口性 别比持续升高，以及乡村人口城镇化等因素，这些都影响着我国人口的增长。2007 年初 发布的国家人口发展战略研究报告对此做出了进一步的分析。 从我国的实际情况和人口增长的上述特点出发，参考相关数据，建立我国人口增长 的数学模型，并由此对我国人口增长的中短期和长期趋势做出预测；特别注意指出模型 中的优点与不足之处。 2. 问题分析问题分析 一个社会（国家、省市、地区）人口的变化和随时间的发展过程，是由很多因素决 定的，社会制度、自然环境、生活水平、科学文化水平、战争、自然灾害和移民等等， 都能严重地影响社会人口的发展过程。然而，婴儿的出生、人口的死亡、居民的迁
5、移却 是决定该社会人口变化的直接原因，近年来我国人口发展出现的一些新特点，如老龄化 进程加速、出生人口性别比持续升高，以及乡村人口城镇化等因素，都直接或间接地通 过这三个现象表现出来。综合考虑这些因素成为构建符合我国国情的人口增长模型关 键。 建立模型对人口发展过程进行定量预测，就是根据现有的人口统计资料和原始数 据，从当前实际的人口状况出发，并对未来的人口发展过程，提出合理的控制要求和假 定，应用科学的方法，预测出未来几年、几十年甚至上百年的人口发展趋势，包括人口 总数、人口的性别、年龄和城乡结构，人口出生、死亡和自然增长率的变化以及在未来 的人口构成中劳力和抚养水平及老化水平等。 3. 模
6、型假设模型假设 针对本问题，建立如下合理的假设： 1. 题中所给数据能反映我国人口变化的基本情况； 2. 一些重大事件，如战争、自然灾害等对人口预测的影响暂不考虑； 3. 所给数据都是年末数据， 也即下一年年初数据， 如 2001 年总人口实质上也表示 2002 年初的总人口； 4. 今年所统计的i岁的人口在下一年年初均为 1 i + 岁； 5. 生育模式不随时间变化。 华中科技大学 2007CUMCM 获奖论文 文峰 郭宪 刘飞，指导教师：高亮 3 4. 符号符号说明说明 符号 说明 单位 ( ) N t 在t年初的总人口数 万人 m N 最大人口容量 万人 r 人口固有增长率 各地区间人口
7、迁移系数 ( ) i x t 第t年年初有i岁的人口数 万人 ( ) i t 第t年i岁人的死亡率 ( ) i g t 第t年i岁人口的迁入量 万人 ( ) i k t 第t年女性人口在i岁人口中所占比例 ( ) i h t 第t年生育模式 ( ) t 第t年总和生育率 1, 2 a a 育龄期范围， 115, 249 aa = 岁 ( ) i t 第t年i岁人的后向死亡率 ( ) i t 第t年i岁人的前向死亡率 ( ) i t 第t年i岁人的迁移率 5. 模型模型建立与求解建立与求解 5.1. 数据预处理数据预处理 题中所给 5 年我国人口 1%调查数据是对人口的抽样调查数据，由于数据的不
8、完备 性，并不能由它来估计当时的全国总人口数。但基于抽样调查的等概率性，可以认为它 所反应的市、镇、乡三个地区的人口比例及男女比例是与实际较为接近的。从中国人 口统计年鉴 2006 1 可以得到 20012005 年具体的全国总人口数。进而可以得到各部分 人口数。所得数据见表 1。 例： 当年抽样调查中市男性总人口 2001年市实际男性总人口=2001年全国总人口 当年抽样调查总人口 表 1 5 年间各部分人口数 年份 2001 2002 2003 2004 2005 城市男 15466 16803 16779 16677 18149 城市女 15420 16804 16843 16889 1
9、8092 镇男 8394 8213 9961 10060 11276 镇女 8154 7904 9705 9901 11117 乡男 41271 40566 39222 39360 36661 乡女 38923 38164 36718 37100 35461 总人口 127627 128453 129227 129988 130756 华中科技大学 2007CUMCM 获奖论文 文峰 郭宪 刘飞，指导教师：高亮 4 5.2. 模型一：基于人口迁移的模型一：基于人口迁移的 Logistic 阻滞增长阻滞增长模型模型 5.2.1. 模型模型建立：建立： 考虑单一地区只受自然资源和环境条件限制的 L
10、ogistic 阻滞增长模型 2 ： ( )( ) ( )(1) m dN tN t rN t dtN = 其中r 为固有增长率。 ( ), m N tN 分别表示第t年的人口数和此条件限制下的最大人口 容量。 因子 ( ) rN t 体现了人口自身的增长趋势， 因子(1( ) m N tN 则体现了资源和环境对 人口增长的阻滞作用，人口增长是两个因子共同作用的结果。 现在考虑市、镇、乡三个地区的总人口变化规律，假定在一段时期内，这三个地区 各自均处于稳定状态，出生率和死亡率都没有太大变化。然而居民在三个地区间的迁移 会对总人口产生一定的影响。对原 Logistic 模型进行修改，得到多地区条
11、件下的改进模 型： 1122 112131 122 3 221 221232 213 33 12 331323 312 ( )( )( )( ) ( )(1) ( ) ( )( )( ) ( )(1) ( )( ) ( )( ) ( )(1) mmm mmm mmm dN tN tN tN t rN t dtNNN N t dN tN tN t r N t dtNNN dN tN t N tN t r N t dtNNN =+ =+ =+ 其中 ,1,2,3. p q = 分别表示市、镇、乡这三个地区， qp 是迁移系数，正数时表示迁 入，负数时表示迁出。 / qpm N 体现了居民从地区q迁
12、到地区 p时对地区 p人口增长的影 响。 为求解此方程，不失一般性的取 1 t = ，则原微分方程组化为差分方程组： 122 11112131 122 3 21 22221232 213 3 12 33331323 312 ( )( )( ) (1)( )( )(1) ( ) ( )( ) (1)( )( )(1) ( ) ( )( ) (1)( )( )(1) mmm mmm mmm N tN tN t N tN trN t NNN N t N tN t N tN tr N t NNN N t N tN t N tN tr N t NNN +=+ +=+ +=+ 5.2.2. 模型模型求解：
13、求解： 1. 求解r 考虑到育龄妇女生育率信息的残缺，用年龄为i的育龄期妇女占总人口的比例 ( ) i p t 与对应每个年龄生育率的乘积和来表示年t的出生率。 2 1 ( )( ) a ii i a rp tt = = 取 2001 t = ,可得到市、镇、乡的固有增长率： 123 0.00907,0.0107,0.01233 rrr = 华中科技大学 2007CUMCM 获奖论文 文峰 郭宪 刘飞，指导教师：高亮 5 2. 求解 m N 由经验公式 0 ( ) (1) m NN tr r = ， 其中 0 0.029 r = ， 为世界人口固有增长率， 可得 2001 年各地区最大人口容量
14、（单位：万人） ： 123 139488,26228,44944 NNN = 。 3. 求解 为求解 ，把 2002 年和 2001 年的数据都带入阻滞增长模型，联合以下三个迁入迁 出方程可解迁移系数 。 （见表） 121212 21 232323 32 313131 13 mm mm mm rr NN rr NN rr NN = = = 表 2 迁移系数 市 1 p = 镇 2 p = 乡 3 p = 市 1 q = 0 21.13720 0.1108025 镇 2 q = 14.55182 0 3.735552 乡 3 q = 0.4674880 22.89316 0 4. 求解 (1)
15、p Nt + 由所求得的系数 , m r N ，把 2001 年作为初始年，用 MATLAB 求解可得以后各年 全国总人口数： 表 3 模型一对全国总人口数的预测 年代 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 预测人数（亿人） 12.94 13.03 13.12 13.22 13.32 13.41 13.50 13.59 13.68 13.77 实际人数（亿人） 12.92 13.00 13.08 误差（） 0.13 0.26 0.39 年代 2013 2014 2015 2016 2017 2018 2019 2020 2021 2022 预测人数（亿人） 13.86 13.94 14.03 14.11 14.19 14.28 14.26 14.44 14.52 14.61 年代 2023 2024 2025 2026 2027 2028 2030 2035 2040 2045 预测人数（亿人） 14.69 14.78 14.86 14.95 15.04 15.12 15.30 15.76 16.24 16.82 5.2.3. 模型分析模型分析： 1. 结果分析 由表 3 可以看出，我国人口在短时间内仍会继续增长，增长水平为每年 1000 万人 左右，到 2010 年和 2020 年，全国总人
