【发布时间】:2012-01-06 17:27:35
【问题描述】:
我想使用基于模型的聚类对 1,225 个时间序列(每个 24 个周期)进行分类。我已经使用快速傅里叶变换分解了这些时间序列,并选择了能够解释样本中所有时间序列的时间序列方差的阈值百分比的谐波。我想对给定时间序列的每个变换元素的实部和虚部进行基于模型的聚类,因为它可能使我不必考虑跨时间序列周期的模型敬畏聚类中的时间自相关。我知道快速傅里叶变换的每个复数元素都独立于其他元素,但我不知道给定输出元素的输出的虚部和实部是否独立。我想知道,因为如果它们是,它将允许我保持 R 中 Mclust 包的默认假设,用于基于模型的聚类分析的变量具有多元高斯分布。
注意:不使用完整的 FFT,因为我已经丢弃了负频率的元素,并通过将频率 1 到奈奎斯特的频率乘以 2 将其从两侧转换为单侧频谱:How do I calculate amplitude and phase angle of fft() output from real-valued input?。
【问题讨论】: