人口模型预测——数学建模作业.docVIP

上传就是为了分析数学得乐趣,请粘贴复制得时候也多思考哈。为了更多得学子们。

2014年数学建模论文

第二套

题目：人口增长模型得确定

专业、姓名：土木135

提交日期：2０1５/7/2晚上

题目:人口增长模型得确定

摘要

对美国人口数据得变化进行拟合，并进行未来人口预测，在第一个模型中,考虑到人口连续变化得规律，用微分方程得方法解出其数量随时间变化得方程，用ｍatｌaｂ里得ｃｆｔool工具箱求出参数，即人口净增长率r=0、02２22,对该模型与实际数据进行对比，并计算了从198０年后每隔１０年得人口数据,与实际对比，有很大出入。因此又改进出更为符合实际得阻滞增长模型,应用微分方程里得分离变量法与积分法解出其数量随时间变化得方程,求出参数人口增长率r=0、028５8与人口所能容纳最大值=258、9,与实际数据对比，拟合得很好,并预测出1980年后每隔10年得人口数据,与实际对比，比较符合。为了便于比较两个模型与实际数据得描述情况作对比,又做出了两个模型与实际数据得对比图,以及两个模型得误差图.

关键词:人口预测微分方程马尔萨斯人口增长模型阻滞增长模型

一、问题重述

1790-1980年间美国每隔10年得人口记录如下表所示。

表1人口记录表

年份

179

0

人口（?106）

3、9

5、3

7、2

９、6

12、９

17、１

23、２

31、4

38、6

50、2

年份

４

70

1980

人口（?106）

62、9

76、0

92、0

１0６、5

１２3、２

1３1、7

150、７

179、3

2０4、０

２２6、5

试用以上数据建立马尔萨斯（Malｔhus)人口指数增长模型,并对接下来得每隔十年预测五次人口数量，并查阅实际数据进行比对分析。

如果数据不相符,再对以上模型进行改进，寻找更为合适得模型进行预测。

二、问题分析

由于题目已经说明首先用马尔萨斯人口增长模型来刻划，列出人口增长指数增长方程并求解，并进行未来5０年内人口数据预测,但发现与实际数据有较大出入.考虑到实际得人口增长率就是受实际情况制约得，因此,使人口增长率为一变化得线性递减函数,列出人口增长微分方程，求出其方程解，并预测未来五十年内人口实际数据。

三、问题假设

１、假设所给得数据真实可靠;

2、各个年龄段得性别比例大致保持不变;

3、人口变化不受外界大得因素得影响;

４、马尔萨斯人口模型

（1)单位时间得人口增长率r为常数;

(2）将视为t得连续可微函数.

5、改进后得模型（阻滞增长模型)

（1）人口净增长率r为变化量。

四、变量说明

ｔ时刻得人口数量

初始时刻得人口数量

r人口净增长率

环境所能容纳得最大人口数量，即

五、模型建立

1。马尔萨斯人口增长模型

t=1790时得人口数为,在t到t＋Δt这一时间间隔内，人口得增长为

由于

则得到可建立含初始条件得微分方程=，=3、9(省略1０^6)

其解为

2.阻滞增长模型

假设人口增长生长率为人口得线性递减函数,即.

假设自然资源与环境条件所能承受得最大人口容量为,显然,当时，。所以。因此有

。于就是建立下列微分方程,

。把上式化为。分离常数并积分得到：。

六、模型求解

1、马尔萨斯模型求解

参数估计:r可以用实际数据得线性最小二乘法求解，对于,直接求解就是比较麻烦得,因此在两边取对数,即,记，＝ａ。则原方程化为（x)=３、９*exｐ（r＊（t－1790）).利用1７90—1900年得数据进行拟合，得到r=0、02142、所以也能求出方程程序见附录1。但本题还可以应用matlａｂ里得ｃftｏol工具箱求参数,在命令行中输入得到更精确得解：

Ｇｅnerａlｍｏｄel：

ｆ（x）=3、９*exp(r*（t—1７9０））

Coeｆfiｃｉｅｎｔs（wiｔh95%coｎｆidenｃｅboｕnds)：

ｒ＝０、02222(0、０２１6３,0、０2281）

得到如图所示结果,其中蓝线表示马尔萨斯人口模型预测人口数据，正方形黑点表示实际人口数据。

图１、马尔萨斯人口模型与实际人口数据

则每隔１0年预测人口为：，,,，，然而查阅相关年份美国实际人口数据，1９90年为2４8、7百万，2000年为２8１、４百万,201０年为307、０百万。对于20２0年与20３0年实际还没有统计，因为没有发生,但通过前三个数据就可以瞧出马尔萨斯模型预测人口与实际有很大出入,所以必须对该模型做出改进,得到更符合实际得预测模型。

２、阻滞增长模型求解

通过对求导得拐点在时,人口增长速度最大。在问题分析已经得到该模型得表达式,运用ｍatlab里得cfｔool工具箱拟