重建周期图中的振铃/振荡

信息处理 功率谱密度 自相关 自回归模型 周期图
2022-02-08 00:45:20

我有一个原始周期图,需要用自回归过程建模。但是该模型不正确,因为它不能很好地拟合原始周期图。

我怀疑我在参加 IFFT 时做错了什么。

  • 标准化原始周期图I(f)使得(归一化频率上的积分总和为 1)
  • 创建一个2N1×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]ej2πf1+a[2]ej2πf2+....+a[p]ej2πfp|2

在哪里p是 AR 模型顺序

在此处输入图像描述

抱歉,我不能展示整个情节是有原因的。我只展示了正面频谱的结束

我究竟做错了什么 ?

0个回答
没有发现任何回复~