时间序列分析报告报告材料方法第06章谱分析报告报告材料

1、第六章谱分析Spectral A nalysis到目前为止，t时刻变量Yt的数值一般都表示成为一系列随机扰动的函数形式，一般的模型形式为：Ytj t jj 0我们研究的重点在于，这个结构对不同时点t和 上的变量Yt和Y的协方差具有什么样 的启示。这种方法被称为在时间域(time domain) 上分析时间序列¥的性质。在本章中，我们讨论如何利用型如 cos( t)和sin( t)的周期函数的加权组合来描述时间 序列Yt数值的方法，这里表示特定的频率，表示形式为：Yt0 ( )cos( t)d 0 ( )sin( t)d上述分析的目的在于判断不同频率的周期在解释时间序列Yt 性质时所发

2、挥的重要程度如何。如此方法被称为频域分析(frequencydomain analysis)或者谱分析(spectralanalysis)。我们将要看到，时域分析和频域分析之间不是相互排斥的，任何协方差平稳过程既有时域表示，也有频域表示，由一种表示可以描述的任何数据性质，都可以利用另一种表示来加以体现。对某些性质来说，时域表示可能简单一些；而对另外一些性质，可能频域表 示更为简单。§6.1母体谱我们首先介绍母体谱，然后讨论它的性质。母体谱及性质假设Yt是一个具有均值的协方差平稳过程，第j个自协方差为：j cov(Yt ,Yt j) E(Yt)(Yt j )假设这些自协方差函数是绝对可

3、加的，则自协方差生成函数为：gY jzjj这里z表示复变量。将上述函数除以2，并将复数z表示成为指数虚数形式则得到的结果(表达式)称为变量Y的母体谱：丄 d2 j j给定任何特定的z exp( i ) ， i1isy( ) 2 gY(e )注意到谱是 的函数:值和自协方差j的序列 j，原则上都可以计算sy()的数值。利用De Moivre 定理，我们可以将e i j cos( j) isin( j)因此，谱函数可以等价地表示成为：1sy( ) 2 jcos( j) isin(j)j表示成为:注意到对于协方差平稳过程而言，有:j j ,因此上述谱函数化简为:sY ()jcos( j) cos(

4、j) isin( j) isin( j)1 12- 0cos(0) isin(0)尸利用三角函数的奇偶性，可以得到:Sy ()02 j cos( j)j i的实值、对称、连续函数。由于对任意期函数，如果我们知道了0,内的所有Sy(假设随机过程Yt(L) t服从MA()过程:Yt这里:(L)j0 jL，I j I ，E( t s) 02, s t0, s t假设自协方差序列 j 是绝对可加的，则可以证明上述谱函数 Sy()存在，并且是2 k，有：Sy(2 k) Sy()，因此 Sy()是周)的值，我们可以获得任意时的Sy()值。§6.2不同过程下母体谱的计算根据前面关于 MA()过程自

5、协方差生成函数的推导:2 1gY(z)(Z) (Z )因此得到MA()过程的母体谱为:Sy( )2- 2 (e i ) (ei )例如，对白噪声过程而言,(z)1，这时它的母体谱函数是常数:Sy()F面我们考虑MA(1)过程,显然，当0时，谱函数Sy()在0,内是的单调递减函数；当0时，谱函数Yttt i此时：(Z)1Z，则母体谱为：Sy()122(1e i )(1ei)12 “ii2、(1ee)2可以化简成为Sy()12122 cos()2Sy()在0,内是的单调递增函数。对AR(1)过程而言，有:Yt 1 t这时只要|Sy()| 1，则有：2(1 e i )(1 1/(1 JR 2Z)，因

6、此谱函数为：2ii2、(1 e e )丄_2 12 cos()该谱函数的性质为：当 0时，谱函数Sy()在0,内是 的单调递增函数；当11Z11ZqZq(11Z)(12Z)(1qZ)pZP (11Z)(12Z) (1 pZ)时，谱函数Sy()在0,内是 的单调递减函数。一般地，对ARMA(p, q)过程而言：Yt c 弭12论2pYtpt1 t 12 t 2则母体谱函数为：2(1i1ei22eqe iq )Sy ()2(1i1e22epe ip )(1i1e22eqeiq )(1i1e22epeip )如果移动平均和自回归算子多项式可以进行下述因式分解:22Z2 2Z则母体谱函数可以表示为:S

7、y ()q21j 12 j cos()p1j 12 j cos()从母体谱函数中计算自协方差如果我们知道了自协方差序列 j，原则上我们就可以计算出任意的谱函数Sy()的数值。反过来也是对的：如果对所有在0,内的，已知谱函数Sy()的数值,则对任意给定的整数 k，我们也能够计算k阶自协方差k。这意味着母体谱函数 Sy()和自协方差序列 j 包含着相同的信息。其中任何一个都无法为我们提供另外一个无法给出 的推断。下面的命题为从谱函数计算自协方差提供了一个有用的公式：命题6.1 假设 j 是绝对可加的自协方差序列，则母体谱函数与自协方差之间的关系为：ksY( )ei kd上述公式也可以等价地表示为:

8、sY ( ) cos( k)d利用上述谱公式，可以实现谱函数与自协方差函数之间的转换。解释母体谱函数o，计算公式为:假设k 0，则利用命题6.1可以得到时间序列的方差，即Sy( )d根据定积分的几何意义，上式说明母体谱函数在区间,内的面积就是 °，也就是过程的方差。更一般的，由于谱函数 Sy()是非负的，对任意 10,，如果我们能够计算：11 Sy( )d这个积分结果也是一个正的数值，可以解释为Yt的方差中与频率的绝对值小于!的成分相关的部分。注意到谱函数也是对称的，因此也可以表示为：1 1Sy( )d 20 Sy( )d这个积分表示频率小于1的随机成分对Yt方差的贡献。但是，频率小

9、于1的随机成分对斗方差的贡献意味着什么？为了探索这个问题，我们考虑更为特殊一些的时间序列模型：MYt jcos( jt) j sin( jt)j 1这里j和j是零均值的随机变量， 这意味着对所有时间t，有EYt 0。进一步假设序列 j；1和 j；1是序列不相关和相互不相关的:j, j kE( j k) j， E( j k)0, j kE( j k)0，对所有的j和k2j0,这时Yt的方差是：M2 2 2 2 2E(Yt ) E( j)cos ( jt) E( j)sin ( jt)j 1M2 2 2j cos ( jt) sin ( jt)j 1M2jj 1j的周期成分对Yt的方差的贡献部分是

10、2。如果频率是有顺序的：0，则Yt的方差中由频率小于或者等于j的周期因此，对这个过程来说，具有频率形成的部分是:这种情形下Yt的k阶自协方差为：ME(YtYtQE( 2)cos( jt)cos j (t k) E(2)sin( jt)sin j (t k)j 1M 2jcos( jt)cos j (t k) sin( jt)sin j (t k) j 1M2 cos( jk)j 1因为过程Yt的均值和自协方差函数都不是时间的函数，因此这个过程是协方差平稳 过程。但是，可以验证此时的自协方差序列 kk 0不是绝对可加的。虽然在上述过程中，我们已经过程的方差分解为频率低于某种程度的周期成分的贡献，

11、 我们能够这样做的原因在于这个过程是比较特殊的。对于一般的情形，著名的谱表示定理(the spectral representation theorem)说明：任何协方差平稳过程都可以表示成为不同频率周期成分的和形式。对任意给定的固定频率0,，我们定义随机变量()和()，并假设可以将一个具有绝对可加自协方差的协方差平稳过程表示为:Yt0 ( )cos( t) ( )sin( t)d这里需要对随机变量()和()的相关性给出更为具体的假设，但是上述公式便是谱表示定理的一般形式。§6.2 样本周期图 Sample Periodogram对一个具有绝对可加自协方差的协方差平稳过程Yt，我们已

12、经定义在频率处的谱函数值为：1 1sy( ) gY(ei )je i J， j E(Yt)(丫 j )22 j注意到母体谱是利用 jj 0表示的，而 jj 0表示的是母体的二阶矩性质。给定由yy2.,yT表示的T个样本，我们可以利用下述公式计算直到(T 1)阶的样本自协方差:J)1T(yty)(yt j y), J 0,1, ,t 1-1 T，yytI t 11, 2, T 1，我们可以获得母体谱密度对应的样本情形，我们称其为样本周期图:T 1?e i J对于给定的1 sY( ) F2 j t 1 _样本周期图也可以表示成为如下形式：1T 1口 2?J cos( j)2j 1我们可以证明样本周

13、期图下的面积等于样本方差:sY()类似地,j 1?y( )d样本周期图也是关于原点对称的，因此也有:?o 2o sY( )d更为重要的是，谱表示定理在样本情形也有类似的表示。我们将要说明，对于平稳过程的任意一个容量为T的观测值序列yy2, ,yi，存在频率1, 2, , m和系数？，?1,?2, , ?M， ?, ?2, ?M使得t期的y值可以表示成为：MYt? ?j cos j(t 1)?jsin j(t 1)j 1其中：当 j k 时，?j cos j(t 1)与?k cos k(t 1)不相关；当 j k 时，?j sin j(t 1)与?ksin k(t 1)不相关；对于所有的j和k，

14、 ?J cos j(t 1)与？kSin k(t 1)不相关。y的样本方差是t 1 ：1(yt y)2，该方差中可以归因于频率为J的周期成分的部分由样本周期图Sy( j)给出。我们对样本容量是奇数的情形展开讨论上述谱表示模式。这时yt可以表示成为由M (T 1)/2个不同频率构成的周期函数，频率1, 2, , m如下：2M"T因此最高频率为:2M2(T1)2T我们考虑yt基于常数项、正弦函数和余弦函数的线性回归:Myt jcos j(t 1) j sin j(t 1) utj i将这个回归方程表示成为下述方式:yt 说 ut其中：xt 1, cos建 1), sin建 1),L ,c

15、os M(t 1),sin m (t 1),这是 一个具有(2M1) T个解释变量的回归方程，因此解释变量与观测值是一样多的。我们将证明解释变量之间是线性无关的，这意味着yt基于Xt回归的OLS估计具有惟一解。该回归方程的 系数具有显著的统计意义：(？2?2)/2表示yt中可以归因于频率j的周期成分的那部分。这就是说，任意观测到的序列,yT，它都可以利用上述周期函数形式表示，并且不同频率的周期成分对方差的贡献都可以在样本周期图中找到。命题6.2 假设样本容量是奇数，定义M (T 1)/2，并设定j 2 i/T ,j 1,2, ,M，假设解释变量为：xt 1, cos & 1), sin

16、建 1),L ,cos m (t 1),sin m (t 1)',1,1 ,2 ,2丄，M , M 则有：TXtXtt 1T 00 (T/2)It 1进一步，假设yy2,y是任意T个实数，则下述推断成立:(a) 过程yt可以表示为:Myt? ?j cos j(t 1) ?sin j(t 1)j 1这里：2 T ? 2 T? y , ?j 二yt cos j(t 1), ?二yt sin j(t 1)T t 1T t 1(b) yt的样本方差可以表示为：TMT(yty)2 2(?2?2)I t 12 j 1样本方差可以归因于频率为j的周期成分的部分为(？2?2)/2。(c) yt的样本方

17、差中可以归因于频率为j的周期成分的部分还可以表示为:?2)*?y( j)其中?y( j)是样本周期图在频率j处的值。T上述结果说明，XtXt是对角矩阵，这意味着包含在向量Xt中的向量之间是相互正交t 1的。这个命题断言：任何奇数个观测到的时间序列Yi，Y2 , , Yt可以表示成为一个常数加上具有(T 1)/2个不同频率的(T 1)个周期成分的加权和。当 T是偶数整数的时候，类似的 结果也是成立的。因此，这个命题给出了类似谱表示定理的有限样本的类似情况。这个命题进一步表明了样本周期图的特征是将y的方差按部分分解为不同频率的周期成分的贡献。注意到解释y的方差的频率j都落在区间0,中。为什么不使用

18、负的频率0 ?假设数据确实是由上述过程的一种特殊情形生成的：Ytcos( t) sin( t)这里0代表某个特殊的负频率，和 是零均值的随机变量，利用三角函数的奇偶性，可以将Yt表示为：Ytcos( t) sin( t)因此，利用上述式子无法从数据中识别数据是从正发频率还是负的频率生成的。这时 一种简单的方式是假设数据是从具有正的频率中生成的。为什么只考虑作为最大的频率呢？假设数据真的是从频率的周期函数中生成的，例如 3 /2 :Ytcos(3 /2)tsin(3 /2)t这时正弦和余弦函数的周期性质表明，上式可以表示成为：Ytcos(/2) tsi n(/2)t因此，根据以前的讨论，具有频率

19、3 /2的周期在观测值上等价于具有频率/2的周期。注意到频率和周期之间的关系，频率对应的周期为2 / 。由于我们考虑的最高频率为，因此我们所观测到的能够自己重复的最短阶段是2 /2。如果 3 /2，则周期是每4/3阶段重复自己。但是，如果数据是整数阶段观测的，因此数据可以观测的时间间隔仍然是每4个阶段观测到，这对应着周期频率是/2。例如，函数cos( /2)t和函数cos(3 /2)t在整数的时间间隔上，它们的观测值是一致的。命题6.2也为计算在频率2 i/T(j 1,2, ,M )上的样本周期图的数值提供了方法。定义:sy(j)工(？2?2)8 (?j?j)这里：?. zj T t 1因此可

20、以得到:TYt cos j (t 1)，2 TYt sin j (t1)?y( j) ATYt cost 1j(t21)TYt sint 12j(t 1)§5.3估计总本谱Estimati ng the Population Spectrum上面我们介绍了母体谱的意义和性质，下面我们面对的问题是：获得了观测样本yi, y2, , yM以后，如何估计母体谱函数 Sy ( ) ?样本周期图的大样本性质一个显然的方法是利用样本周期图Sy()去估计母体谱函数 SY()。但是，这种方法具有显著的限制。假设对于无限移动平均过程而言：Ytj t jj 0这里系数 jj o是绝对可加的， Jt 是具有均值E( t)0和方差var( J 2的