我有一个原始周期图,需要用自回归过程建模。但是该模型不正确,因为它不能很好地拟合原始周期图。
我怀疑我在参加 IFFT 时做错了什么。
- 标准化原始周期图I(f)使得(归一化频率上的积分总和为 1)
- 创建一个2N−1×1两侧周期图N标准化周期图样本
- 采取I(f)−−−−√最后一步的信号
- 从 [0,2pi] 生成随机相位值,并将信号点与相位中的相量相乘
- 执行IFFT(Scipy IFFT,0频率内容-正面内容-负面内容)
- 自相关生成的 IFFT 信号
- 使用自相关方法求解 AR 参数
- 最后拟合AR模型并生成Periodogram
I(f)=1|A^(f)|2=1|a[0]+a[1]e−j2πf1+a[2]e−j2πf2+....+a[p]e−j2πfp|2
在哪里p是 AR 模型顺序
抱歉,我不能展示整个情节是有原因的。我只展示了正面频谱的结束
我究竟做错了什么 ?