数学建模的时间序列建模(软件应用)

实验名称:时间序列建模

实验目的:

学会通过各种手段检验序列的平稳性;学会根据自相关系数和偏自相关系数来初步判断ARMA模型的阶数p和q,学会利用最小二乘法等方法对ARMA模型进行估计,学会利用信息准则对估计的ARMA模型进行诊断,以及掌握利用ARMA模型进行预测。掌握在实证研究中如何运用Eviews软件进行ARMA模型的识别、诊断、估计和预测和相关具体操作。

实验内容:

(1)根据时序图判断序列的平稳性;

(2)观察相关图,初步确定移动平均阶数q和自回归阶数p;

3)运用经典B-J方法对某企业201个连续生产数据建立合适的ARMA()模型,并能够利用此模型进行短期预测。

(1)数据录入

   打开Eviews软件,选择“File”菜单中的“New--Workfile”选项,在“Workfile structure type”栏选择“Unstructured /Undated”,在“Date range”栏中输入数据个数201,点击ok,这样就建立了一个工作文件。点击File/Import,找到相应的Excel数据集,打开数据集导入数据。

(2)绘制序列时序图

    双击序列production,点击view/Graph/line,则出现序列时序图,时序图看出201个连续生产的数据是平稳的,这个判断比较粗糙,需要用统计方法进一步验证。

(3)绘制序列相关图

双击序列production,点击view/Correlogram,我们对原始数据序列做相关图,因此在“Correlogram of”对话框中选择“Level”即表示对原始序列做相关,在滞后阶数中选择14(

),点击ok。

    从相关图看出,自相关系数迅速衰减为0,说明序列平稳,但最后一列白噪声检验的Q统计量和相应的伴随概率表明序列存在相关性,因此序列为平稳非白噪声序列。我们可以对序列采用B-J方法建模研究。

(4)ADF检验序列的平稳性

通过时序图和相关图判断序列是平稳的,我们通过统计检验来进一步证实这个结论,双击序列production,点击view/unit root test,出现对话框,我们对序列本身进行检验,序列不存在明显的趋势,所以选择对常数项,不带趋势的模型进行检验,其他采用默认设置,点击ok,出现图2-7的检验结果,表明拒绝存在一个单位根的原假设,序列平稳。

(5)模型定阶

     由图2-5看出,偏自相关系数在k=3后很快趋于0即3阶截尾,尝试拟合AR(3);自相关系数在k=1处显著不为0,当k=2时在2倍标准差的置信带边缘,可以考虑拟合MA(1)或MA(2);同时可以考虑ARMA(3,1)模型等。

    在序列工作文件窗口点击View/Descriptive Statistics/Histogram and States对原序列做描述统计分析见下图,可见序列均值非0,我们通常对0均值平稳序列做建模分析,所以需要在原序列基础上生成一个新的0均值序列。点击主菜单Quick/Generate Series,在对话框中输入赋值语句Series x=production-84.11940,点击ok则生成新序列x,这个序列是0均值的平稳非白噪声序列,相当于在原序列基础上作了个整体平移,所以统计特性没有发生根本改变。我们对序列x进行分析。

2、模型参数估计

(1)尝试AR模型。经过模型识别所确定的阶数,可以初步建立AR (3),可用菜单或命令两种方式分别建立。在主菜单选择Quick/Estimate Equation,出现图2-10的方程定义对话框,在方程定义空白区键入x  ar(1)  ar(2)  ar(3)  ,其中ar(i)(i=1,2…)表示自回归系数;估计方法选择最小二乘估计(LS)。

即对应的模型为:

(2)尝试MA模型。按上面介绍方法,方程定义空白区键入x  ma(1)(其中ma(j),j=1,2…代表移动平均系数)或在主窗口输入ls  x  ma(1)  ma(2) 。模型输出结果见下图。

即对应的模型为:

(3)尝试ARMA模型

由模型定阶发现,p可能等于3,q可能等于2或1,经过进一步筛选,逐步剔除不显著的滞后项或移动平均项,最后得到如下ARMA(2,1)模型,参数估计结果见下图

ARMA(2,1)模型拟合图如下显示

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值