时间序列ARIMA期末论文-ARIMA模型在总人口预测中的应用.docx

ARIMA模型在总人口预测中的应用
【摘要】人口发展与社会经济的发展是密不可分的,研究我国总人口的发展,对我国人口数进行分析和预测,有利于及时控制人口的增长调节人口平衡,利于政府及时了解发展趋势并做出反应对策使我国人口发展步入健康的轨道。本文利用时间序列建模原理和思路，并结合R3.2.1软件对1962年——2014年我国年底总人口数据做分析和预测。找到对原始数据有着较好的拟合度和较高的预测精度的模型。利用此模型可对我国年底总人口进行合理的预测。
【关键词】ARIMA建模总人口人口预测

目录

一、引言 3
1.1研究背景 3
1.2研究现状 4
二、模型建立 5
2.1模型识别 5
2.2模型的参数估计 8
2.3模型的诊断 10
2.模型的预测 12
三、 模型的优缺点及推广 13
3.1模型的优缺点 13
3.2模型的推广 13
结束语.............................................................14
【参考文献】........................................................ 15
附录...............................................................16

一、 引言
1.1研究背景
我国是世界上人口最多的国家，自1980年开始，年末中国大陆总人口就已经超过了10亿，并一直保持约占世界总人口的五分之一，亚洲人口的三分之一。中国人口的发展同中国社会的发展一样经过了漫长而曲折的道路。在世纪的进程中，目前我国进入了一个全新的时
代，要想在21世纪——这个充满竞争与挑战的时代中变的富强、屹立于世界民族之林，实现我们的中国梦，这全取决于人。能否顺利解决人口现状等问题，是我国乃自世界共同面临的问题，由于地球的资源是有限的，它不可能无限制的容纳人口，当人口过多，会由于经济跟不上，工作岗位欠缺，医疗等水平不足，从而导致整个社会处于一种动荡之中；然而如果人口过少，又会由于人员不足，导致各方面人力资源不足，无法正常完成各项必须社会活动，这也会极大地限制一个国家的发展，因此，对人口的研究是具有相当的意义的。
我国由于幅员广阔，民族众多，各民族发展水平不一，同时作为世界第一人口大国，我国的耕地面积却相对不足，因此我国每年都需要从国外大量进口粮食，由于过分依赖于进口这对我国的发展影响巨大，为此甚至有国外反华势力叫嚣只要断绝给中国供粮，三五年之内中国必定大乱。当然那只是敌对势力的一厢情愿与恶意诋毁，但我们自己却必须认识到在由于人口的问题而导致的一系列问题，关于人口问题我国必须重视,并根据其趋势做出反应对策。因此，认真分析我国当前人口现状，从中发现其变化的趋势，并对未来总人口进行短期预测，及时采取必要的政治及经济措施来解决人口发展问题，对树立未来的发展目标很有必要。总之，人口是构成社会的主体，在我国社会主义现代化建设中，人口问题始终是极为重要的问题，而人口问题的本质是发展问题。人口发展与社会经济的发展也是密不可分的。基于此，我们利用时间序列中的ARMA模型对我国人口进行预测，对人口的控制起到指导作用，有利于政府采取必要的政治及经济措施来进行调控。所以，对其进行分析和测试是非常有意义的工作。

1.2研究现状
在对人口问题的研究上，国内外学者做了相当多的工作。在国内程等利用自限模型对我国的人口增长进行了预测，认为中国在2010年-2019年人口数依次会缓慢增加，2016年突破14亿大关，且未来15年人口净增加量不会超过1亿；蒋慧基于多元统计模型对广西人口增长进行了分析，得出了人口增长的综合因子，并提出了稳定人口增长的建议；丁明等运用相空间重构神经网络模型对我国人口增长进行预测，很好的解决了非线性的问题，为我国人口增长预测提供了一种新的方法；王保等，利用Logistic模型进行人口预测，并检验了2005年—2007年的数据误差，取得了理想的效果。在国外，Rosen
利用Malthusian模型对人口进行研究，也取得了不错的成绩。本文基于时间序列在研究时间相关问题上的优势，以1949年-2014年的年末总人口数据，利用时间序列知识建模，找到适合人口增长的过程的时间序列模型（模型识别），然后利用参数估计估计出模型的参数（参数估计），再对模型进行诊断，判断模型的好坏（模型诊断），最后利用已经建立的模型对未来的给定的时间进行预测（预测）。

二、 模型建立
2.1模型识别
首先需要对数据进行预先的处理，观察其时间序列图是否为平稳序列，可以用R软件完成（具体程序见于附录1），可得到时间序列图为：

图1：人口时间序列图
由图一可以看出，年底总人口数随着时间增加的同时也在逐年上涨，有着明显的上升趋势。因此可得出这列数据是不平稳的、方差也是不平稳的结论。并且数据大致是呈线性变化的，因此可以考虑做差分变换。
先对数据做一阶差分变换后再观察序列是否平稳，可运用R软件编程得其变换后的序列及变换后序列的自相关图和偏自相关图(具体程序见附录2)：

图二：一阶差分序列图及其ACF、PACF图
由图二可知，对原始数据在进行一阶差分之后的时间序列图显示并不算平稳，并且一阶差分后的ACF图仍呈现出近似直线下降趋势，因此可以考虑再做一次差分运算。通过编程可以得到（具体程序见附录3）：

图三：对数二阶差分序列图及其ACF、PACF图
再观察其二阶差分后的时间序列图，基本上趋于平稳，而其自相关图(ACF)和偏自相关图（PACF）呈现出明显的拖尾形式，且PACF图在滞后6阶比较显著，ACF在滞后1,3,6阶相对较显著，由此可以认为原序列基本上可以用ARIMA(1,2,1)，ARIMA(1,2,2)，ARIMA(1,2,3)，进行拟合。
对ARIMA(1,2,1)模型序列满足：
Yt+Yt-2-2Yt-1=∅Yt-1+Yt-3-2Yt-2+et-i=16θiet-i ①
由①式可以得到：
Yt=2+∅Yt-1+1-2∅Yt-2-∅Yt-3+et-θet-1 ②
同理可以得到 ARIMA(1,2,2)模型：
Yt=2+∅Yt-1+1-2∅Yt-2-∅Yt-3+et-i=12θiet-i ③
ARIMA(1,2,3) 的模型为：
Yt=2+∅Yt-1+1-2∅Yt-2-∅Yt-3+et-i=13θiet-i ④

2.2模型的参数估计
根据2.1节，已经找到几个可能用于拟合的模型，分别为ARIMA(1,2,1)，ARIMA(1,2,2)，ARIMA(1,2,3)，那么接下来就应该对其进行参数估计，在对模型进行参数估计时有多种方法可以选择，这里选择用最小二乘估计，通过R编程序可以得到如下结果（具体程序见于附录4）：
Call:
arima(x = x, order = c(1, 2, 1), method = "CSS")
Coefficients:
ar1 ma1
-0.4794 0.5150
s.e. 0.1033 0.1901

sigma^2 estimated as 21752: part log likelihood = -327.05

Call:
arima(x = x, order = c(1, 2, 2), method = "CSS")
Coefficients:
ar1 ma1 ma2
-0.5216 0.8442 -0.3147
s.e. 0.0009 0.0912 0.1014
sigma^2 estimated as 16939: part log likelihood = -320.67

Call:
arima(x = x, order = c(1, 2, 3), method = "CSS")
Coefficients:
ar1 ma1 ma2 ma3
-0.5203 0.7646 -0.1111 0.3405
s.e. 0.0008 0.1327 0.1435 0.1278
sigma^2 estimated as 15191: part log likelihood = -317.89
由上面的输出可以确定在②式中的系数，并且由估计值的误差项se判断该系数是否显著（通过判定在0是否在区间[∅-2se,∅+2se]内，若在则系数不显著，如不在则显著）。
从而代入②可以确定ARIMA(1,2,1)模型为：
Yt=1.5206Yt-1-0.0412Yt-2-0.4794Yt-2+et-0.515et-1 ⑤
类似代入③式得ARIMA(1,2,2)表达式为：
Yt=1.4784Yt-1+0.0432Yt-2-0.5216Yt-2+et+0.8442et-1-0.3147et-2 ⑥
代入④式得ARIMA(1,2,3)模型的表达式：
Yt=1.4797Yt-1+0.0406Yt-2-0.5203Yt-2+et+0.7646et-1-0.3405et-3 ⑦

2.3模型的诊断
在对模型完成了识别和参数估计之后，需要对模型进行诊断，诊断模型是否具有对原时间序列数据的很好的拟合效果。主要进行拟合模型的残差分析和分析过度参数化；对模型进行过度参数分析主要看在进行差分时是否出现过度差分的状况，而对残差进行分析主要需要做以下几个方面的工作：
⑴，检验残差是否是随机的，一个模型如果能很好的拟合，那么拟合后的残差基本上是随机的，残差是应该围绕在某条平行于x=0这一条直线上下波动的，并且波动的幅度不会很大，这可以用做残差序列图观察得到；
⑵，检验残差是否呈正态性，一个模型如果能很好的拟合，那么其残差应该是呈正态性的，这里用残差QQ图和S-W正态性检验（原假设为：H0:数据是呈正态性的）进行；
⑶，判断残差之间是否是相互独立的，一个模型如果能很好的拟合，那么其残差之间相对是比较独立的，这里主要用残差的自相关序列图和L-B检验（原假设为：H0:原数据的残差之间是不相关的）进行。
诊断ARIMA(1,2,1)模型，用R软件编程序（具体程序见于附录5）输出为：

图4：ARIMA(1,2,1)模型残差的序列图、ACF图和QQ图
Shapiro-Wilk normality test

W = 0.93417, p-value = 0.00592

Box-Ljung test
X-squared = 23.846, df = 23, p-value = 0.4124
由图4中的残差序列图可以看出残差基本基本上是围绕x=0这条直线上下波动的，因此可以认为ARIMA(1,2,1)模型拟合满足残差是随机的条件；又由图4中的QQ可以看出残差基本上是集中在一条直线上的，由S-W检验的的p=0.00592<0.05，检验也可以认为残差是非正态的；又在残差的自相关图中，只有滞后二阶是是显著的，因此可以认为残差之间基本上也是不相关的，特别由B-L检验的p=0.4121>0.05，因此没有充分的理由拒绝原假设，应该认为残差是相互独立的。到此就已经对模型行进行了诊断，由于残差正态性不足，说明用ARIMA(1,2,1)模型拟合原数据不是十分合适。
同理可对ARIMA(1,2,2)和ARIMA(1,2,3),进行诊断，这里图形和数据检验就不再一一呈现在论文中。在诊断中发现模型ARIMA(1,2,2)用于拟合是满足条件的，而ARIMA(1,2,3)模型也由于残差正态性不足而不适合用于拟合合。因此在对人口模型进行预测时采用ARIMA(1,2,1)模型。并且在进行参数冗余分析时发现，对于该时间序列用ARIMA(1,2,2)模型拟合后，对ARIMA(p,d,q)中p,d,q任意一个变小都不能再满足条件，因此用ARIMA(1,2,2)模型是适合且简化的，因此后文选用ARIMA(1,2,2)模型进行预测。

2.4模型预测
由2.1,2.2,2.3三节已经完成了对模型的识别、参数估计和模型诊断，现在就需要运用该模型对原时间序列趋势进行预测。在实际情况中，过去的状态已经是即成的事实，我们关心的更多的是未来是什么状态，从而就可以判断在目前的情况未来的情况，进而可以进行相应的措施予以应对。这里假设对未来5年进行预测，编程[7]得到（具体程序见附录6）：
Time Series:
Start = 2015
End = 2019
Frequency = 1
[1] 137487.2 138208.4 138936.0 139666.2 140397.4

Time Series:
Start = 2015
End = 2019
Frequency = 1
[1] 162.5725 305.4282 518.2338 796.2256 1129.1609
由上输出，可知由模型ARIMA (3,2,6)预测的2015年-2019年年末总人口数分别为：137487.2,138208.4 ,138936.0,139666.2,140397.4（单位：万人）。并且还可以对人口数进行区间估计，在此就不列出。

三、 模型优缺点及推广
3.1模型优缺点
优点：ARIMA模型在对许多时间序列都适用，并且在建模过程中有多种方法可以选择，并且在模型的诊断中能过对模型的好坏进行评价，对预测起到很好的作用。
缺点：ARIMA模型对数据进行预测时，对短期预测效果还是比较好的，但随着时间的延长，其预测误差就比较大了。
3.2模型的推广
ARIMA模型对大多以时间为变量的趋势都有比较好的拟合效果，可广泛运用于气温、股票、降水量、商品价格等的研究之中。
结束语
人口问题是人类社会伴随始终的问题，在新中国成立之初，由于对人口问题的认识不足，导致我国大量的人员过剩，从而影响到后面几十年的发展，给教育、医疗社会带来了巨大的压力，即便后来在意识到问题的严重性后，在计划生育政策下，人口得到了有效的控制，然而由于人口基数过大，人们传统观念的根深蒂固，我国人口还是处于不断增长的时期，最近几年又由于伴随着人口老龄化情况加剧，社会经济压力又有了很大的压力，这对我国的经济在一定程度上有着不良影响，慢慢政府在对计划生育问题上又有了一定放松，在一定条件下允许第二胎。这一系列的重大举措反应了当代人对人口问题的重视。本文对未来几年利用模型进行了预测，对政策是一定的补充说明。
论文结束之际，在此向各位老师表示感谢，对在相关领域做出贡献的学者表示感谢。

参考文献：
[1]程华，高孝成，幺焕民《我国人口增长预测问题的自限模型》，[J],哈尔滨师范大学自然科学学报，2009年，第25卷第3期：34-36
[2]蒋慧萍，《基于多元统计模型的广西人口增长分析》，[J],企业科技与发展，2010年，第10期：190-192
[3]丁明磊，杨晓娜，曹连海《相位空间重构神经网络模型在我国人口增长预测中的应用》，[J],华北水利水电学院学报，2008年6月，第29卷第3期：95-98
[4]王保学，蔡果兰，《Logistic模型的 参数估计及人口预测》[J]，北京工商大学学报（自然科学版），2009年11月，第27卷第6期：75-79
[5] Rosen F,《The principle of population as political theory: Godwin's Of population and the Malthusian controversy》,[J],Journal Of The History Of Ideas [J Hist Ideas] 1970 Jan-Mar; Vol. 31 (1), pp. 33-48.
[6]潘红宇等译，《时间序列分析及应用》[M]，机械工业出版社，2011年1月
[7]薛毅，陈丽萍，《统计建模与R软件（下册）》[M]，清华大学出版社

附录
1、做时间序列的序列图
plot(x)

2、做序列的一阶差分时间按序列图及处理后的ACF图,PACF图程序：
mydata<-read.table("C:/Users/Administrator/Desktop/shuju.txt",header=T)
x<-ts(mydata,frequency =1, start =1962)
par(mfrow=c(3,1))
plot(x)
plot(diff(x,difference=1)) #做差分后的时间序列图
acf(diff(x,difference=1),ci.type='ma') #做差分后的自相关图
pacf(diff(x,difference=1)) #做差分后的偏自相关图

3、做序列的二阶差分时间按序列图及处理后的ACF图,PACF图程序：
mydata<-read.table("C:/Users/Administrator/Desktop/shuju.txt",header=T)
x<-ts(mydata,frequency =1, start =1962)
par(mfrow=c(3,1))
plot(diff(x,difference=2))#做差分后的时间序列图
acf(diff(x,difference=2),ci.type='ma') #做差分后的自相关图
pacf(diff(x,difference=2)) #做差分后的偏自相关图

4、估计模型中的参数程序：
mydata<-read.table("C:/Users/Administrator/Desktop/shuju.txt",header=T)
x<-ts(mydata,frequency =1, start =1962)
arima(x,order=c(1,2,1),method="CSS")#CSS代表条件最小二乘估计
arima(x,order=c(1,2,2),method="CSS")#CSS代表条件最小二乘估计
arima(x,order=c(1,2,3),method="CSS")#CSS代表条件最小二乘估计

5、对模型的残差进行检验程序：
library(TSA)
mydata<-read.table("C:/Users/Administrator/Desktop/shuju.txt",header=T)
x<-ts(mydata,frequency =1, start =1962)
m<-matrix(c(1:3,2),2,2)
layout(m);layout.show(3);
model=arima(x,order=c(1,2,1)) #当需要改诊断其他模型时就改变此处的p,d,q值即可
plot(rstandard(model,infl = iflSR),ylab ='Standardized Residuals', type='o');abline(h=0);#画残差图
qqnorm(residuals(model))#画残差的QQ图
qqline(residuals(model))#画残差图QQ图
shapiro.test(residuals(model))#对残差做正态性检验
acf(residuals(model))#对残差做自相关图
LB.test(model,lag=25)#求残差Ljung-Box统计量

6、对未来5年年末人口数的预测程序：
library(TSA)
mydata<-read.table("C:/Users/Administrator/Desktop/shuju.txt",header=T)
x<-ts(mydata,frequency =1, start =1962)
model=arima(x,order=c(1,2,2))#确定模型
predict(model,c(5))#对未来五年在model模型下的预测

1. 基于C8051F单片机直流电动机反馈控制系统的设计与研究
2. 基于单片机的嵌入式Web服务器的研究
3. MOTOROLA单片机MC68HC（8）05PV8/A内嵌EEPROM的工艺和制程方法及对良率的影响研究
4. 基于模糊控制的电阻钎焊单片机温度控制系统的研制
5. 基于MCS-51系列单片机的通用控制模块的研究
6. 基于单片机实现的供暖系统最佳启停自校正（STR）调节器
7. 单片机控制的二级倒立摆系统的研究
8. 基于增强型51系列单片机的TCP/IP协议栈的实现
9. 基于单片机的蓄电池自动监测系统
10. 基于32位嵌入式单片机系统的图像采集与处理技术的研究
11. 基于单片机的作物营养诊断专家系统的研究
12. 基于单片机的交流伺服电机运动控制系统研究与开发
13. 基于单片机的泵管内壁硬度测试仪的研制
14. 基于单片机的自动找平控制系统研究
15. 基于C8051F040单片机的嵌入式系统开发
16. 基于单片机的液压动力系统状态监测仪开发
17. 模糊Smith智能控制方法的研究及其单片机实现
18. 一种基于单片机的轴快流CO〈,2〉激光器的手持控制面板的研制
19. 基于双单片机冲床数控系统的研究
20. 基于CYGNAL单片机的在线间歇式浊度仪的研制
21. 基于单片机的喷油泵试验台控制器的研制
22. 基于单片机的软起动器的研究和设计
23. 基于单片机控制的高速快走丝电火花线切割机床短循环走丝方式研究
24. 基于单片机的机电产品控制系统开发
25. 基于PIC单片机的智能手机充电器
26. 基于单片机的实时内核设计及其应用研究
27. 基于单片机的远程抄表系统的设计与研究
28. 基于单片机的烟气二氧化硫浓度检测仪的研制
29. 基于微型光谱仪的单片机系统
30. 单片机系统软件构件开发的技术研究
31. 基于单片机的液体点滴速度自动检测仪的研制
32. 基于单片机系统的多功能温度测量仪的研制
33. 基于PIC单片机的电能采集终端的设计和应用
34. 基于单片机的光纤光栅解调仪的研制
35. 气压式线性摩擦焊机单片机控制系统的研制
36. 基于单片机的数字磁通门传感器
37. 基于单片机的旋转变压器-数字转换器的研究
38. 基于单片机的光纤Bragg光栅解调系统的研究
39. 单片机控制的便携式多功能乳腺治疗仪的研制
40. 基于C8051F020单片机的多生理信号检测仪
41. 基于单片机的电机运动控制系统设计
42. Pico专用单片机核的可测性设计研究
43. 基于MCS-51单片机的热量计
44. 基于双单片机的智能遥测微型气象站
45. MCS-51单片机构建机器人的实践研究
46. 基于单片机的轮轨力检测
47. 基于单片机的GPS定位仪的研究与实现
48. 基于单片机的电液伺服控制系统
49. 用于单片机系统的MMC卡文件系统研制
50. 基于单片机的时控和计数系统性能优化的研究
51. 基于单片机和CPLD的粗光栅位移测量系统研究
52. 单片机控制的后备式方波UPS
53. 提升高职学生单片机应用能力的探究
54. 基于单片机控制的自动低频减载装置研究
55. 基于单片机控制的水下焊接电源的研究
56. 基于单片机的多通道数据采集系统
57. 基于uPSD3234单片机的氚表面污染测量仪的研制
58. 基于单片机的红外测油仪的研究
59. 96系列单片机仿真器研究与设计
60. 基于单片机的单晶金刚石刀具刃磨设备的数控改造
61. 基于单片机的温度智能控制系统的设计与实现
62. 基于MSP430单片机的电梯门机控制器的研制
63. 基于单片机的气体测漏仪的研究
64. 基于三菱M16C/6N系列单片机的CAN/USB协议转换器
65. 基于单片机和DSP的变压器油色谱在线监测技术研究
66. 基于单片机的膛壁温度报警系统设计
67. 基于AVR单片机的低压无功补偿控制器的设计
68. 基于单片机船舶电力推进电机监测系统
69. 基于单片机网络的振动信号的采集系统
70. 基于单片机的大容量数据存储技术的应用研究
71. 基于单片机的叠图机研究与教学方法实践
72. 基于单片机嵌入式Web服务器技术的研究及实现
73. 基于AT89S52单片机的通用数据采集系统
74. 基于单片机的多道脉冲幅度分析仪研究
75. 机器人旋转电弧传感角焊缝跟踪单片机控制系统
76. 基于单片机的控制系统在PLC虚拟教学实验中的应用研究
77. 基于单片机系统的网络通信研究与应用
78. 基于PIC16F877单片机的莫尔斯码自动译码系统设计与研究
79. 基于单片机的模糊控制器在工业电阻炉上的应用研究
80. 基于双单片机冲床数控系统的研究与开发
81. 基于Cygnal单片机的μC/OS-Ⅱ的研究
82. 基于单片机的一体化智能差示扫描量热仪系统研究
83. 基于TCP/IP协议的单片机与Internet互联的研究与实现
84. 变频调速液压电梯单片机控制器的研究
85. 基于单片机γ-免疫计数器自动换样功能的研究与实现
86. 基于单片机的倒立摆控制系统设计与实现
87. 单片机嵌入式以太网防盗报警系统
88. 基于51单片机的嵌入式Internet系统的设计与实现
89. 单片机监测系统在挤压机上的应用
90. MSP430单片机在智能水表系统上的研究与应用
91. 基于单片机的嵌入式系统中TCP/IP协议栈的实现与应用
92. 单片机在高楼恒压供水系统中的应用
93. 基于ATmega16单片机的流量控制器的开发
94. 基于MSP430单片机的远程抄表系统及智能网络水表的设计
95. 基于MSP430单片机具有数据存储与回放功能的嵌入式电子血压计的设计
96. 基于单片机的氨分解率检测系统的研究与开发
97. 锅炉的单片机控制系统
98. 基于单片机控制的电磁振动式播种控制系统的设计
99. 基于单片机技术的WDR-01型聚氨酯导热系数测试仪的研制
100. 一种RISC结构8位单片机的设计与实现
101. 基于单片机的公寓用电智能管理系统设计
102. 基于单片机的温度测控系统在温室大棚中的设计与实现
103. 基于MSP430单片机的数字化超声电源的研制
104. 基于ADμC841单片机的防爆软起动综合控制器的研究
105. 基于单片机控制的井下低爆综合保护系统的设计
106. 基于单片机的空调器故障诊断系统的设计研究
107. 单片机实现的寻呼机编码器
108. 单片机实现的鲁棒MRACS及其在液压系统中的应用研究
109. 自适应控制的单片机实现方法及基上隅角瓦斯积聚处理中的应用研究
110. 基于单片机的锅炉智能控制器的设计与研究
111. 超精密机床床身隔振的单片机主动控制
112. PIC单片机在空调中的应用
113. 单片机控制力矩加载控制系统的研究
项目论证，项目可行性研究报告，可行性研究报告，项目推广，项目研究报告，项目设计，项目建议书，项目可研报告，本文档支持完整下载，支持任意编辑！选择我们，选择成功！
项目论证，项目可行性研究报告，可行性研究报告，项目推广，项目研究报告，项目设计，项目建议书，项目可研报告，本文档支持完整下载，支持任意编辑！选择我们，选择成功！
单片机论文，毕业设计，毕业论文，单片机设计，硕士论文，研究生论文，单片机研究论文，单片机设计论文，优秀毕业论文，毕业论文设计，毕业过关论文，毕业设计，毕业设计说明，毕业论文，单片机论文，基于单片机论文，毕业论文终稿，毕业论文初稿，本文档支持完整下载，支持任意编辑！本文档全网独一无二，放心使用，下载这篇文档，定会成功！

19