zoukankan      html  css  js  c++  java
  • SPSS时间序列:频谱分析

    一、频谱分析(分析-预测-频谱分析)

        “频谱图”过程用于标识时间序列中的周期行为。它不需要分析一个时间点与下一个时间点之间的变异,只要按不同频率的周期性成分分析整体序列的变异。平滑序列在低频率具有更强的周期性成分;而随机变异(“白噪声”)将成分强度分布到所有频率。不能使用该过程分析包含缺失数据的序列。

    1、示例。建造新住房的比率是一个国家/地区经济的重要晴雨表。有关住房的数据开始时通常会表现出一个较强的季节性成分。但在估计当前数字时,分析人员需要注意数据中是否呈现了较长的周期。

    2、统计量。正弦和余弦变换、周期图值和每个频率或周期成分的谱密度估计。在选择双变量分析时:交叉周期图的实部和虚部、余谱密度、正交谱、增益、平方一致和每个频率或周期成分的相位谱。

    3、图。对于单变量和双变量分析:周期图和频谱密度。对于双变量分析:平方一致性、正交谱、交叉振幅、余谱密度、相位谱和增益。

    4、数据。变量应为数值型。

    5、假设。变量不应包含任何内嵌的缺失数据。要分析的时间序列应该是平稳的,任何

    非零均值应该从序列中删除。

    平稳. 要用ARIMA 模型进行拟合的时间序列所必须满足的条件。纯的MA 序列是平稳

    的,但AR 和ARMA 序列可能不是。平稳序列的均值和方差不随时间改变。

    二、频谱图(分析-预测-频谱分析)

    1、选择其中一个“频谱窗口”选项来选择如何平滑周期图,以便获得谱密度估计值。可用的平滑选项有“Tukey-Hamming”、“Tukey”、“Parzen”、“Bartlett”、“Daniell(单元)”和“无”。

    1.1、Tukey-Hamming. 权重为Wk = .54Dp(2 pi fk) + .23Dp (2 pi fk + pi/p) + .23Dp (2pi fk - pi/p),k = 0, ..., p,其中p 是一半跨度的整数部分,Dp 是阶数p 的Dirichlet 内核。

    1.2、Tukey. 权重为Wk = 0.5Dp(2 pi fk) + 0.25Dp (2 pi fk + pi/p) + 0.25Dp(2 pi fk -pi/p),k = 0, ..., p,其中p 是一半跨度的整数部分,Dp 是阶数p 的Dirichlet 内核。

    1.3、Parzen. 权重为Wk = 1/p(2 + cos(2 pi fk)) (F[p/2] (2 pi fk))**2,k= 0, ... p,其中p 是一半跨度的整数部分,而F[p/2] 是阶数p/2 的Fejer 内核。

    1.4、Bartlett. 谱窗口的形状,窗口上半部分的权重按如下公式计算:Wk = Fp(2*pi*fk),k = 0, ... p,其中p 是半跨度的整数部分,Fp 是阶数p 的Fejer 内核。下半部分与上半部分对称。

    1.5、Daniell(单元). 所有权重均等于1 的频谱窗口形状。

    1.6、无. 无平滑。如果选择了此选项,则频谱密度估计与周期图相同。

    2、跨度. 一个连续值范围,在该范围上将执行平滑。通常使用奇数。较大的跨度对谱密度图进行的平滑比较小的跨度程度大。

    3、变量中心化. 调整序列以使在计算谱之前其均值为0,并且移去可能与序列均值关联的较大项。

    4、图。周期图和谱密度对单变量分析和双变量分析均可用。其他所有选项仅对双变量分析可用。

    4.1、周期图. 针对频率或周期绘制的未平滑谱振幅图(绘制在对数刻度中)。低频率变动是平滑序列的特征。均匀地分布在所有频率上的变动则表示“白噪音”。

    4.2、平方一致性. 两个序列的增益的乘积。

    4.3、正交谱. 交叉周期图的虚部,是两个时间序列的异相频率成分的相关性的测量。成分的异相为pi/2 弧度。

    4.4、交叉振幅. 余谱密度平方和正交谱平方之和的平方根。

    4.5、谱密度. 已进行平滑而移去了不规则变动的周期图。

    4.6、余谱密度. 交叉周期图的实部,是两个时间序列的同相频率分量的相关性的测量。

    4.7、相位谱. 一个序列的每个频率成分提前或延迟另一个序列的程度的测量。

    4.8、增益. 用一个序列的谱密度除以跨振幅的商。这两个序列都有自己的获得值。

  • 相关阅读:
    CJOJ 1308 【HNOI 2002 】营业额统计 / CodeVS 1296 营业额统计(STL,二分)
    CJOJ 1131 机器分配 / Luogu 2066 机器分配 (动态规划)
    CJOJ 2482 【POI2000】促销活动(STL优先队列,大根堆,小根堆)
    CJOJ 2484 函数最小值 / Luogu 2085 函数最小值(STL优先队列,堆)
    POJ 2296 Map Labeler / ZOJ 2493 Map Labeler / HIT 2369 Map Labeler / UVAlive 2973 Map Labeler(2-sat 二分)
    Luogu 1111 修复公路(最小生成树)
    POJ 3683 Priest John's Busiest Day / OpenJ_Bailian 3788 Priest John's Busiest Day(2-sat问题)
    POJ 3207 Ikki's Story IV
    洛谷 P1456Monkey King
    洛谷 P1231教辅的组成
  • 原文地址:https://www.cnblogs.com/China3S/p/3863766.html
Copyright © 2011-2022 走看看