人口增长模型的确定.docx
- 文档编号:23729761
- 上传时间:2023-05-20
- 格式:DOCX
- 页数:14
- 大小:37.70KB
人口增长模型的确定.docx
《人口增长模型的确定.docx》由会员分享,可在线阅读,更多相关《人口增长模型的确定.docx(14页珍藏版)》请在冰豆网上搜索。
人口增长模型的确定
Documentnumber:
NOCG-YUNOO-BUYTT-UU986-1986UT
人口增长模型的确定
题目:
人口增长模型的确定
摘要
人口问题已成为当前世界上最普遍关注的问题之一,人口增长规律的发现以及人口增长的预测问题对一个国家制定长远的发展规划有着非常重要的意义。
本文分别使用了马尔萨斯人口指数增长模型和阻滞增长模型,以美国1790-1980年间每隔10年的人口数量为依据,对接下来的每隔十年进行了预测五次人口数量。
通过对比我们可以发现阻滞增长模型在预测准确度方面要明显优于原始的马尔萨斯人口指数增长模型。
关键词:
人口增长;马尔萨斯人口指数增长模型;阻滞增长模型;人口预测
一、问题重述
问题背景
1790-1980年间美国每隔10年的人口记录如下表所示。
表1人口记录表
年份
1790
1800
1810
1820
1830
1840
1850
1860
1870
1880
人口(106)
年份
1890
1900
1910
1920
1930
1940
1950
1960
1970
1980
人口(106)
问题提出
我们需要解决以下问题:
1.试用以上数据建立马尔萨斯(Malthus)人口指数增长模型,并对接下来的每隔十年预测五次人口数量,并查阅实际数据进行比对分析。
2.如果数据不相符,再对以上模型进行改进,寻找更为合适的模型进行预测,并对两次预测结果进行对比分析。
3.查阅资料找出中国人口与表1同时期的人口数量,用以上建立的两个模型进行人口预测与分析。
二、问题分析
首先,我们运用Matlab软件绘制出1790到1980年的美国人口数据图,如图1。
图11790到1980年的美国人口数据图
从图表中我们可以清晰地看到人口数在1790—1980年是呈增长趋势的,而且我们很容易发现上述图表和我们学过指数函数的图表有很大的相似性,所以我们很自然想到建立指数模型。
因此我们首先建立马尔萨斯模型,马尔萨斯生物总数增长定律指出:
在孤立的生物群体中,生物总数N的变化率与生物总数成正比。
三、问题假设
为简化问题,我们做出如下假设:
(1)在模型中预期的时间内,人口不会因发生大的自然灾害,突发事件或战争而受到大的影响;
(2)所给出的数据具有代表性,能够反映普遍情况;
(3)一段时间内我国人口死亡率不发生大的波动;
(4)在查阅的资料与文献中,所得数据可信;
(5)假设人口净增长率为常数。
四、变量说明
在此,对本文所使用的符号进行定义。
表2变量说明
符号
符号说明
N(0)
起始年人口容纳量
N(t)
t年后人口容纳量
t
年份
r
增长率
五、模型建立
问题一:
马尔萨斯(Malthus)人口指数增长模型
设:
t表示年份(起始年份t=0),r表示人口增长率,N(t)表示t年后的人口数量。
当考察一个国家或一个很大地区的人口时,N(t)是很大的整数。
为了利用微积分这一数学工具,将N(t)视为连续、可微函数。
记初始时刻(t=0)的人口为N(0),人口增长率为r,r是单位时间内N(t)的增量与N(t)的比例系数。
根据r是常数的基本假设,于是N(t)满足如下的微分方程:
dN(t)/dt=r*N(t)(5-1)
由这个线性常系数微分方程容易解出:
N(t)=N(0)ert(5-2)
表明人口将按指数规律无限增长(r>0)。
将以t年为单位,上式表明,人口以er为公比的等比数列增长。
因为这时r表示年增长率,通常r<<1,所以可用近似关系er≈1+r可得出
N(t)=N(0)(1+r)t(5-3)
(5-3)式即人口增长模型。
问题二:
改进模型-阻滞增长模型(Logistic模型)
自从英国人口学家和政治经济学家托马斯·罗伯特·马尔萨斯1798年发表《人口学原理》后,马尔萨斯(Malthus)人口指数增长模型在世界上引起了轩然大波,并在后来的人口预测中扮演着重要的角色。
但是随着时间的发展,由于现代社会与自然环境的改变,马尔萨斯人口指数增长模型在预测未来人口时,误差可能会比较大。
上述模型对较早时期的统计数据吻合得较好,但也存在问题,即人口是呈指数规律无止境地增长,此时人口的自然增长率随人口的增长而增长,这不可能。
由于社会的快速发展,自然环境遭受严重破坏,人口的高速增长等一系列原因,人口的增长率不能按照马尔萨斯所假设为一个常数r不改变。
一般说来,当人口较少时增长得越来越快,即增长率在变大;人口增长到一定数量以后,增长就会慢下来,即增长率变小这是因为,自然资源、环境条件等因素不允许人口无限制地增长,它们对人口的增长起着阻滞作用,而且随着人口的增加,阻滞作用越来越大。
而且人口最终会饱和,趋于某一个常数x,我们假设人口的净增长率为r(1-x(t)/x),即人口的净增长率随着人口的增长而不断减小,当t时,净增长率趋于零。
按照这个假设,得到:
(5-4)
这便是荷兰数学家Verhulst于19世纪中叶提出的阻滞增长模型(Logistic模型)。
在MATLAB命令窗口键入
dsolve(‘Dx=r*x*(1-x/c)’,’x(1790)=’)
输出:
ans=c/(1+1/39*exp(-r*t)*exp(1790*r)*(10*c-39))
其中c=x
因此,人口的变化规律为:
(5-5)
问题三模型建立
经调查,1790-2010年间中国每隔10年的人口记录如下表所示。
表7中国人口记录表
年份
1790
1800
1810
1820
1830
1840
1850
1860
1870
人口(106)
年份
1880
1890
1900
1910
1920
1930
1940
1950
1960
人口(106)
年份
1970
1980
1990
2000
2010
人口(106)
我们分别应用马尔萨斯人口指数增长模型和Logistic人口阻滞增长模型来对中国人口进行预测。
六、模型求解
问题一模型求解
在应用预测模型的过程中考虑到,若要提高预测结果的准确性,就必须增加预测方案的数量,对比各方案的预测值和误差,选取误差最低的一组预测方案。
特别是马尔萨斯模型中,人口增长率r是一定时期内人口增加的综合结果,在预测中它的取值直接关系到预测结果的精度,因此在进行不同阶段的人口预测时根据实际情况对人口增长率r加以分类和处理才能得到理想的预测结果。
本文根据1790-1980年计算美国常住人口每年的增长率,按照人口增长率r的大小设置了高中低三个方案,以此加强预测结果的对比,提高预测的准确度。
表3美国每10年自然增长率
序号
1
2
3
4
5
6
7
8
9
r
10
11
12
13
14
15
16
17
18
19
通过表3可以确定自然增长率高,中,低三个方案。
通过数据分析可得,上述表格为10年的累计增长率,而自然增长率强调一年,所以可近似除以10求得,高方案中自然增长率为,中方案中自然增长率为,低方案中自然增长率为。
依据人口增长率的大小分为高、中、低、三个预测方案,将预测值与实际值进行拟合比较。
图2r=时马尔萨斯模型曲线拟合
图3r=马尔萨斯模型曲线拟合
图4r=马尔萨斯模型曲线拟合
根据上述分析,及曲线拟合可知,取中方案即r=时,马尔萨斯模型更符合实际情况。
因此本文自然增长率取r=来预测美国人口数量并与实际情况对比。
由预测公式预测1790-1980年的人口数量,由指数增长模型可得各个年份的真实值与预测值之间的差别如下表:
表41790年-1980年美国人口真实值与预测值
年份
1790
1800
1810
1820
1830
1840
1850
1860
1870
1880
真实人口(106)
误差(106)
0
年份
1890
1900
1910
1920
1930
1940
1950
1960
1970
1980
真实人口(106)
误差(106)
通过调查得知1990-2010年人口数量统计如下表
表51990年-2010年美国人口真实值与预测值
时间(年)
1990
2000
2010
人口
误差(106)
939
1297
图5美国人口真实值与预测值曲线拟合
通过上图可以发现,1790-1870指数增长模型确实拟合的比较好,但从1870年开始往后发现误差越来越大,可知指数增长模型只适合于短期的人口预测。
为了生存以及人类的发展,人们自然会采取有效措施来控制人口的过度增长,自然资源、环境资源的条件也限制了人口数量的过度增长。
因此为了使人口预报模型适合长期的发展趋势,更好地符合实际情况,必须修改指数增长模型中关于人口增长率是常数这个基本假设了,这时必将导致更适合人类发展的规律的新数学模型的产生。
问题二模型求解
利用MATLAB软件中的“curvefit”命令和式(5-5)来拟合所给的人口统计数据,从而确定出(5-5)中的待定参数r和x。
查阅资料可得r的初值取为小于1的数,比如取a=[200,]时,得到a=[],y1=,即(5-5)中的r=,x=,2010年美国的人口预计为百万人。
这个结果还比较合理,当t趋于无穷时,静增长率趋于零,人口数趋于百万人,即极限人口x=百万百万。
拟合效果见图5。
根据该题已给数据可作如下图形:
图6Logistic模型拟合曲线
从图6可以看出,在前一段吻合得比较图,但在最上面,若拟合曲线更接近原始数据,对将来人口的预测应该更好。
因此略加修改将拟合准则改为:
(5-6)
其中w为右端几个点的误差权重,在此处应该取为大于1的数,这样会使右边的拟合误差减小,相应的,其他点的误差会有所增加。
我们要使这些误差的增减恰当,可以通过调整w和n的具体取值,比较他们取各种不同值时的拟合效果,从而确定出一个合适的数值。
1)先取n=17,w=,运行上述程序,得到结果a=[,];
x1=.
2)再取n=16,w=2,运行上述程序,得到结果a=[,];x1=.
我们把两种情况的拟合曲线画在同一个坐标系中,很容易作出比较,见图6。
第二种情形后半段的变化趋势与原始数据更吻合,因此,对将来人口的预测应该更好。
图7Logistic模型优化拟合曲线
经过修改,得到了一个较满意的结果,人口增长率r=,极限人口xm=(百万),并预测1990年--2010年美国人口。
通过调查得知1990-2010年人口数量统计如下表
表61990年-2010年Logistic模型美国人口真实值与预测值
时间(年)
1990
2000
2010
人口
预测值
误差
问题三模型求解
马尔萨斯(Malthus)人口指数增长模型求解中国人口问题
参照问题一,我们来求解中国人口问题。
表8中国每10年自然增长率
序号
1
2
3
4
5
6
7
8
9
r
0
10
11
12
13
14
15
16
17
18
19
同样,取高、中、低三种自然增长率方案,高方案中自然增长率为,中方案中自然增长率为,低方案中自然增长率为。
图8r=时马尔萨斯模型曲线拟合
图9r=时马尔萨斯模型曲线拟合
图10r=时马尔萨斯模型曲线拟合
根据上述分析,及曲线拟合可知,取中方案即r=时,马尔萨斯模型更符合实际情况。
因此本文自然增长率取r=来预测中国人口数量并与实际情况对比。
由预测公式预测1790-1980年的人口数量,由指数增长模型可得各个年份的真实值与预测值之间的差别如下表:
表91790年-1980年中国人口真实值与预测值
年份
1790
1800
1810
1820
1830
1840
1850
1860
1870
1880
真实人口(106)
381
409
412
412
377
358
368
误差(106)
0
61
139
240
364
548
782
1107
1487
1925
年份
1890
1900
1910
1920
1930
1940
1950
1960
1970
1980
真实人口(106)
380
400
423
472
489
误差(106)
2471
3144
3982
5004
6391
7944
9973
12410
15439
19228
通过调查得知1990-2010年人口数量统计如下表
表101990年-2010年中国人口真实值与预测值
时间(年)
1990
2000
2010
人口
误差(106)
23987
29965
37480
图11中国人口真实值与预测值曲线拟合
Matlab中cftool()工具箱求解
应用Matlab中cftool()工具箱来进行图像拟合。
图12cftool工具箱拟合
经过不同模型的应用比对,我们发现应用3阶高斯分布可以达到较好的拟合效果,误差相对较小。
此时,a1=1002,1339)
b1=1992(1978,2006)
c1=,
a2=,
b2=1932(1907,1958)
c2=,
a3=,
b3=1837(1828,1846)
c3=,
七、结果分析
综合做出假设的两种模型与原始数据所描述的图形如下:
可以看出,当世界人口总数不大时,生存空间,资源等极充裕,人口总数指数的增长是可能的,但当人口总数非常大时,指数增长的线性模型则不会反映这样的现实。
阻滞增长模型从一定程度上克服了指数增长模型的不足,而指数增长模型在做人口的短期预测时因为其形式的相对简单性也常被采用,但是依据实际情况,根据已给数据和查资料所得数据通过图像拟和可得二次函数模型与其最为匹配。
事实上,人口的预测是一个相当复杂的问题,影响人口增长的因素除了人口基数与可利用资源量外,还和医药卫生条件的改善、人们生育观念的变化等因素有关,特别在做中短期预测时,我们希望得到满足一定预测精度的结果,比如在刚刚经历过战争或是由于在特定的历史条件下采纳了特殊的人口政策等,这些因素本身以及由此引起的年龄结构变化就变得相当重要,也要予以考虑。
八、参考文献
[1]张志涌,杨祖樱.MATLAB教程[M].北京:
北京航空航天大学出版社,2011.
[2]李晓梅.人口预测模型研究及应用[M]四川:
西南财经大学出版社,2011.
[3]姜启源等《数学模型》(第三版)[M].高等教育出版社P10
[4]何春.马尔萨斯人口模型在广州市人口预测中的应用[J].广东工业大学学报,2010,(03):
31-34.
[5]孙鹏.两种人口预测方法在孝感市人口预测中的应用[J].时代农机,2017
(2):
99-101.
[6]张金明,李骞.基于马尔萨斯模型的北京市人口预测[J].特区经济,2013(7):
63-65.
[7]金鑫,高丽娟,金辉,等.基于马尔萨斯模型对甘肃省人口的预测[J].齐齐哈尔大学学报(自然科学版),2012,28(5):
90-92.
九、附录
程序1指数增长模型:
functionzhishuzengzhang
clearall;
t=1790:
10:
2010;r=;
N=[].*10^6;
N
(1)=*10^6;Nm=3*10^7;r=;
N1=N
(1)*((1+r).^(t-1790));
N1/10^6
plot(t,N,'k*',t,N1,'B-')
legend('原始数据','指数模型')
程序2Logistic模型:
functionf=fun3(a,t)
f=a
(1)./(1+(a
(1)/*exp(-(t-1790)*a
(2)));
x=1790:
10:
1990;
y=[92];
plot(x,y,'*',x,y);
a0=[,1];
a=curvefit('fun3',a0,x,y)
xi=1790:
10:
2020;
yi=fun3(a,xi);
holdon
plot(xi,yi);
x1=2010;
y1=fun3(a,x1)
holdoff
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 人口增长 模型 的确