数模第一次作业-(1).docVIP

2016年数学建模论文 第 套 论文 题目： 专业、姓名： 专业、姓名： 专业、姓名： 提交日期： 2016.6.27 题目：人口增长模型的确定 摘 要 对美国人口数据的变化进行拟合，并进行未来人口预测，在第一个模型中，考虑到人口连续变化的规律，用微分方程的方法解出其数量随时间变化的方程，先求对数用matlab里线性拟合求出参数，即人口净增长率r=0.0214,对该模型与实际数据进行对比，并计算了从1980年后每隔10年的人口数据，与实际对比，有很大出入。因此又改进出更为符合实际的阻滞增长模型，应用微分方程里的分离变量法和积分法解出其数量随时间变化的方程，求出参数人口增长率r=0.0268和人口所能容纳最大值=285.89,与实际数据对比，拟合得很好，并预测出1980年后每隔10年的人口数据，与实际对比，比较符合。为了便于比较两个模型与实际数据的描述情况作对比，又做出了两个模型与实际数据的对比图，并计算了误差。 关键词：人口预测 微分方程 马尔萨斯人口增长模型 阻滞增长模型 一、问题重述 1790-1980年间美国每隔10年的人口记录如下表所示： 年份 1790 1800 1810 1820 1830 1840 1850 1860 1870 1880 人口(?106) 3.9 5.3 7.2 9.6 12.9 17.1 23.2 31.4 38.6 50.2 年份 1890 1900 1910 1920 1930 1940 1950 1960 1970 1980 人口(?106) 62.9 76.0 92.0 106.5 123.2 131.7 150.7 179.3 204.0 226.5 表1 人口记录表 试用以上数据建立马尔萨斯(Malthus)人口指数增长模型，并对接下来的每隔十年预测五次人口数量，并查阅实际数据进行比对分析。 如果数据不相符，再对以上模型进行改进，寻找更为合适的模型进行预测。 二、问题分析 由于题目已经说明首先用马尔萨斯人口增长模型来刻划，列出人口增长指数增长方程并求解，并进行未来50年内人口数据预测，但发现与实际数据有较大出入。考虑到实际的人口增长率是受实际情况制约的，因此，使人口增长率为一变化的线性递减函数，列出人口增长微分方程，求出其方程解，并预测未来五十年内人口实际数据。 三、问题假设 1.假设所给的数据真实可靠; 2.各个年龄段的性别比例大致保持不变; 3.人口变化不受外界大的因素的影响； 4.马尔萨斯人口模型 （1）单位时间的人口增长率r为常数； （2）将视为t的连续可微函数。 5.改进后的模型（阻滞增长模型） (1)人口净增长率r为变化量。 四、变量说明 t时刻的人口数量 初始时刻的人口数量 人口净增长率 环境所能容纳的最大人口数量 表2：变量说明 五、模型建立 1．马尔萨斯人口增长模型 t=1790时的人口数为,在t到t+Δt这一时间间隔内，人口的增长为 由于 则得到可建立含初始条件的微分方程 2．阻滞增长模型 由韦尔候斯特假定，马尔萨斯模型应该改为 上式就是逻辑模型，该方程分离变量，其解为 六、模型求解 1.马尔萨斯模型求解 对每年的人口数取对数，用线性拟合求出和r，计算误差和对以后每隔十年进行人口预测。 程序结果为 p = 0.0214 -36.6198 N0 = 1.2480e-16 r = 0.0214 RM = 1.6975e+04 Malthus1 =405.5324 502.4002 622.4063 771.0778 955.261 图1.马尔萨斯人口模型与实际人口数据 得到马尔萨斯模型为 2.阻滞增长模型求解 我们假定美国人口上限为400，最大增长率为3%。 我们用lsqcurvefit进行参数拟合并得到人口上限和最大增长率。 结果为： x = 285.8973 0.0286 RM = 658.8203 Verhulst1 = 230.9164 242.5097 252.0172 259.6666 265.7272 图2 ；阻滞人口模型与实际人口数据 所以我们得到阻滞人口模型表达式为： 七、结果分析 1.马尔萨斯模型结果分析 从1990年每隔10年一直到2030年预测人口。结果为:405.5324 502.4002 622.4063 771.0778 955.2618。然而查阅相关年份美国实际人口数据，1990年为248.7百万，2000年为281.4百万，2010年为307.0百万。对于20