非季节的差异

打开实时脚本

此示例显示了如何对时间序列进行非季节差异。时间序列是从1947年到2005年的季度美国GDP。

加载工具箱中包含的GDP数据集。

加载data_gdpy =数据;n =长度（y）;图图（y）xlim（[0，n]）标题（'我们。GDP'）

图包含一个轴对象。带有标题的Axes对象，USY GDP包含类型线的对象。

时间序列有明显的上升趋势。

采取该系列的第一区别以消除趋势，

$δ y_{t} = （（ 1 - l ） y_{t} = y_{t} - y_{t - 1} 。$

首先创建一个差异滞后操作员多项式对象，然后使用它来过滤观察到的串联。

d1 = lagop（{1，-1}，“滞后”，[0,1]）;dy =滤波器（d1，y）;图图（2：n，dy）xlim（[0，n]）标题（“第一个差异GDP系列”）

图包含一个轴对象。带有标题的轴对象首先差的GDP系列包含一个类型行的对象。

该系列在首次差异后仍然存在一些上升趋势。

在系列中取得第二个区别，

$δ^{2} y_{t} = （（ 1 - l ）^{2} y_{t} = y_{t} - 2 y_{t - 1} + y_{t - 2} 。$

d2 = d1*d1;ddy =滤波器（d2，y）;图图（3：n，ddy）xlim（[0，n]）标题（“第二差GDP系列”）

图包含一个轴对象。带有标题第二个GDP系列的轴对象包含类型行的对象。

第二差异的系列看起来更固定。

也可以看看