谱分析方法.doc

_第4章 谱分析方法1 绪论一 时间序列模型：通过分析自相关就获得描述与预测时间序列可能够用模型的第一印象。如这里与相关性较大，而与相关较弱，为什么？二分析时间序列的两种方法频谱法， 时间序列法-Box Jenkins方法三 时间序列模型的五个特征（最重要的）描述趋势有多种方法1 趋势 确定性趋势 随机趋势2 季节性： 是季节哑变量，定义为， ， 其它3 异常观测值异常观测值：在时间序列中，可能有一个或几个点，会对时间序列的建模与预测起到重要的作用。这样的数据点称为奇异观测值。4 条件异方差异常观测值倾向于成群出现，这个现象称为波动性集聚（vilatility clustering）条件异方差 5 非线性： 状态依赖机制转换特征2 谱分析一 时间序列分析的方法1 时序分析方法：也就是时序建模方法，ARMA等，也就是原序列的时间顺序不变。2 频谱建模方法：单变量频谱建模技术就是时间序列看作是有不同频率的正弦和余弦波组成。其基本思想是：把时间序列看作是互不相关的周期（频率）分量的叠加，通过研究和比较各分量的周期变化，以充分揭示时间序列的频域结构，掌握其主要波动特征。 做法：对某个时间序列剔除趋势和季节因素后的循环项（平稳）进行谱估计，根据估计出的普密度函数，找出序列中的主要频率分量，从而把握该序列的周期波动特征。优点：当频率分量的行为或内在机制互不相同时，谱分析可以避免时域方法带来的混淆（因为时域方法所衡量的只是各频率分量共同叠加后的结果）更精细的研究各种行为及因素。二 基本原理1 对于确定性函数，周期为，如果在可积，并在点连续，则有如下的傅立叶级数展开：其中 意义：周期函数通常可以分解为常数项与频率的正弦和余弦函数之和。为了方便，令从而得到复数形式的傅氏级数展开 其中为振幅，为频率因此，周期函数可以表示成不同频率及其对应的振幅的正弦和余弦函数之和，这就是的谱表示。在周期内消耗的能量等于每一个不同频率的三角函数分量所消耗的能量之和：单位时间上的能量消耗的功率为即的总功率等于各频率分量的功率之和。功率谱：功率依不同频率的分布，表示为称此函数为功率谱密度函数。三一般确定性非周期函数的谱表示构造以2T为周期的函数，满足则当在内可积，并在t点连续的情况下，进行傅立叶级数展开，在内：这里，当。如果积分存在，则对任何t，有的功率为这里，则是在频率f处的能量密度，或称为连续能量谱密度。而且非周期函数的功率为。3 平稳过程的频域分析一 因为无法保证实现的周期性和可积性，因而需采取相应的手段。对平稳过程的实现加以截取，构造新函数其傅立叶展开为这里是的能量谱密度函数。的总能量无限，但功率却可能有限。称为功率功率谱密度函数。如果功率谱密度函数的期望存在，则称为平稳过程的功率谱密度函数，简称功率谱或谱密度。称称为谱分布函数。二自协方差函数和谱密度的关系定理1：平稳过程的谱密度存在，则过程的自协方差函数有如下的傅立叶变换即如果的自协方差函数绝对可积，则的谱密度必定存在。方差：，因为自相关函数，即方差代表平稳过程的总平均功率。 标准化功率谱密度函数，他代表频率范围在内的分量对总功率的贡献率。即标准化谱密度是自相关函数的傅立叶变换。 谱密度的性质：1）2）对任何，。3）对实值过程，对于任何f，满足。积分谱具有下列性质：1） 。2） ，。3） 是f的单调非减函数，即当时，。定理2（维纳-辛钦定理）：自相关函数是某个平稳（连续）随机过程的自相关函数的充分必要条件是：存在一个函数，在上具有分布函数的性质（即，且是单调非减），使得对一切，可以表示成：分别称为自相关函数和自协方差函数的谱表示。 如果是绝对连续的分布函数，则对任何f都存在谱密度函数称随机过程具有纯连续谱。如过程。三 平稳时间序列的频域分析1 平稳时间序列是平稳随机过程的特殊情况。具有两种特殊性质：1） t只取整数值，因此自协方差和自相关函数只在整数点有定义。2） 谱密度只在范围内有定义。2 Wold定理：序列作为某个时间序列的自相关函数的充分必要条件是：存在一个上的单调非减函数，使得自相关函数自协方差称为非标准化积分谱。为标准化积分谱。当序列为实值，自相关函数为偶函数时，标准化谱密度函数：例：求白噪声的谱密度其自协方差函数则表明纯随机序列具有常值谱密度函数，意味着总功率或方差在上是均匀分布的，即每个频率成分对总功率或方差的贡献是一样的。一般情况下，纯连续谱的时间序列的谱密度可划分为三种类型：1） 谱密度从频率0到频率1/2递减，高谱密度值集中在相对低频处，表明序列异常周期波动为主。2） 在高频处显示高谱密度，说明序列以短周期波动为主，比白噪声还不规则的随机过程。3） 谱密度主要集中在某个特定频率附近，意味着序列的变动主要是由这个频率所确定的周期波动。 4 功率谱密度函数的估计方法一几种常用指标设为具有纯连续谱的平稳时间序列，其谱密度为，记为根据样本所得到的的估计量。1偏差当，则称此估计为无偏估计。当，为渐进无偏估计。2方差 ， 它描述了偏离均值的程度，方差越小，估计量越好。3均方误差：时， ，则是在均方意义下的一致估计。二估计方法：非参数方法和参数方法（一）非参数方法：周期图和窗谱估计1周期图估计：研究观测数据中可能隐藏的周期性模型为其中从而估计的基本思想（寻找隐含周期的基本思路）：用一组足够密的备选频率作选择，得到参数估计值，做出振幅的平方关于频率的图形，根据平方振幅足够大处（即图形中较明显的峰值）作对应的备选频率，便可识别出隐含周期频率，并同时得到相应的振幅。根据是否显著不同于零，可以判别其对应频率是否是隐含周期的频率。定义周期图： 其中；。因此在备选频率处，周期图为。而且等价 周期图恰好是从样本数据出发通过样本协方差函数的谱密度进行直接估计的统计量。它有以下特点：1）时，周期图方差趋向于谱密度的平方，即在均方误差意义下，不是一致估计；2）对于不同的频率f，是渐进不相关的，当增大时，特别不稳定，不平滑。3 窗谱估计：对周期图进行截断，即对 的两头的若干求和项去掉，从而得到谱密度的新估计量：，其中M为小于N-1的正整数，称这个M为截断点，而上面的估计式称为截断周期图。更一般的形式为其中为的权函数，称为时窗函数或时窗。为窗谱估计。（二）参数方法：ARMA和极大熵谱估计1AR和ARMA谱估计基于ARMA模型的谱估计方法是近代谱分析中比较感兴趣的一种方法。1）AR谱估计假设一个时间序列是由过程生成的，即其理论谱密度为其估计量为这种方法就是所谓的自回归谱估计或AR谱估计。2）ARMA谱估计其理论谱密度为其估计量为这种方法就是所谓的自回归谱估计或AR谱估计。3） 极大熵估计（现代谱分析）熵是一种不规则性的度量，现在用于表示不确定性的度量。可用来度量随机变量或随机过程的不确定程度。（1） 熵：称为谱熵。谱熵越大，则的随机性越强，所受人为干预越少，包含信息越多。（2） 极大熵准则，对于满足约束条件使谱熵达到最大的，记为（极大熵谱估计）。4 谱分析在经济中的应用须注意的问题1） 序列长度要大：至少100-200个数据。2） 数据需要预先处理：要求平稳序列、去掉长期趋势以及可能的季节性变化。3） 谱估计方法的选择。