时间序列分析matlab_平稳时间序列分析之模型识别

时间序列识别过程利用自相关系数ACF和偏自相关系数PACF，其中ACF用于衡量序列中较早的数据值是否与后面时间值有某种关系，PACF用于捕获变量和变量滞后量之间的关系。

计算出样本自相关系数和偏相关系数的值之后，就要根据它们表现出来的性质，确定相应的p，q的数值，从而选择适当阶数的 ARMA(p,q)模型。因此，模型识别过程也称为模型定阶过程。

1、ARMA模型定阶的基本原则为：

由于样本的随机性，样本的自相关系数不会呈现出理论截尾的完美情况，本应截尾的样本自相关系数或偏自相关系数人会呈现小值振荡。同时，由于平稳时间序列通常都具有短期相关性，随着延迟阶数k的增大，样本自相关系数或偏自相关系数都会衰减值零值附近作小指波动。

这种现象促使我们必须思考：

(1)当样本自相关系或偏自相关系数在延迟若干阶之后衰减为小值波动时，什么情况下该看做相关系数截尾？

(2)什么情况下该看做相关系数在延迟若干阶之后衰减到零值附近作拖尾波动？

这实际上没有绝对的标准，很大程度上依靠分析人员的主观经验。但样本自相关系数和偏自相关系数的近似分布可以帮助分许人员做出尽量合理的判断。

如果样本自相关系数或偏相关系数在最初的d阶明显超过2倍标准差范围，而后几乎95%的自相关关系都落在2倍标准差范围以内，而且由非零自相关系数衰减为小值波动的过程非常突然，这时，通常视为自相关系数截尾，截尾阶数为d。

如果有超过5%的样本系数落入2倍标准差范围之外，或者有显著非零的自相关系数衰减为小值波动的过程比较缓慢或者非常连续，这时，通常视为自相关系数不截尾。

2、案例分析

选择合适的ARMA模型拟合1880—1985年全球气表平均温度改变值差分序列。

①读入数据，绘制时序图

a