时间序列的算法ARMA算法的参数设定的问题

matlab的程序

autocorr(Series) %画出自相关图,图中上下两条横线分别表示自相关系数的上下界,超出边界的部分表示存在相关关系。
[a,b] = autocorr(Series) %a 为各阶的相关系数,b 为滞后阶数
parcorr(Series) %画出偏自相关图
[c,d] = parcorr(Series) %c 为各阶的偏自相关系数,d 为滞后阶数

子相关系数 偏相关系数 模型定阶
拖尾 p阶截尾 AR(p)模型
q阶截尾 拖尾 MA(q)模型
拖尾 拖尾 ARMA(p,q)模型

其实,观察自相关图与偏相关图最主要的目的还是确定序列的ARMA(p,q)模型的具体形式。首先,需要明确这样几对概念:
第一,自回归过程与移动平均过程。自回归由序列的滞后变量的线性组合以及白声噪(符合0均值固定方差的随机干扰项)相加而成,移动平均过程为白声噪的线性组合构成;
第二,拖尾和截尾。这一对概念从图表上很容易看出,前者指AC或者PAC呈几何衰减(指数式衰减或者正弦式衰减),后者指AC或者PAC在某一阶之前明显不为0,之后突然接近或者等于0.其实,从字面上也很好理解,拖尾就是拖拖拉拉,截尾就是抽刀断水。
其次是对ARMA模型的分解。
AR(p)模型,从自相关函数ACF来看,在自回归方程的基础上可以很简单地构造自相关系数,最后发现自相关系数等于w^k(w为自回归系数),对于平稳时间序列(注意这一前提条件,如果放开这一条件图形将会很难识别),|w|<1,所以当w>0时,ACF呈现为指数式衰减至0。当w<0时,ACF则正负交替呈指数衰减至0,整体表现则是正弦式衰减;从偏相关函数PACF来看,这就相当明显了,因为PACF与自回归方程的形式完全一样,只是自回归方程只有滞后p期,而PACF则有更多的滞后项。于是乎,很明显,当k<=p,偏相关系数不等于0,当k>p,偏相关系数等于0,明显呈现出截尾现象。
MA(q)模型,从自相关函数ACF来看,在移动平均方程的基础上也可以很简单地构造自相关系数,这时候的自相关函数为分段函数,当k<=q,偏相关系数不等于0,当k>q,偏相关系数等于0,明显呈现出截尾现象;从偏相关函数PACF来看,任何一个可逆的MA(q)过程都可以转换成一个无限阶、系数按几何衰减的AR过程(将白噪声替换为序列的滞后形式即可),呈现拖尾现象。与AR(p)不同的是,当v>0(v为移动平均系数)时,PACF呈现为交替式正弦衰减。当v<0时,PACF则呈指数衰减至0。ARMA(p,q)模型则是两者的结合,实际判别p、q值时还是比较依赖经验的。

下表给出了不同 ARMA 模型的自相关函数和偏自相关函数:
这里写图片描述

这里写图片描述

这里写图片描述

参考的网址:
https://wenku.baidu.com/view/025c8e3350e2524de4187e2e.html
http://www.ilovematlab.cn/thread-41758-1-1.html
http://bbs.pinggu.org/thread-3096618-1-1.html

  • 1
    点赞
  • 9
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值