频谱估计算法的计算复杂度是多少?它们依赖什么?

信息处理 fft 频谱 功率谱密度 估计 频谱估计
2022-02-06 19:56:20

有很多频谱估计技术,每一种都有一些优点和缺点。算法如:

  • 音乐
  • 韦尔奇的方法
  • Yule-Walker AR 方法
  • 周期图
  • 修正协方差法
  • 多锥度法 (MTM)
  • 使用短时间的频谱图
  • 傅里叶变换
  • Burg 法 协方差法

和许多其他的。但是我不知道哪个更适合嵌入式系统实现,因为我不知道它们的计算复杂性,所以我对这种方法的计算复杂性的数学估计更感兴趣,特别是它们的复杂性如何受输入噪声等因素的影响,频谱占用等?

我找到了Spectrum Estimaton Methods,它极大地解释了这些方法,但没有从计算复杂性的角度比较它们。

1个回答

[*] 中讨论了各种非参数 PSD 估计方法的计算要求。复杂度取决于长度N您的数据记录。粗略地讲:

  1. 计算 DFT 幅度平方的最直接的周期图方法N样本的复杂度为NlogN.
  2. Bartlett 和 Welch 方法使用M重叠的长度窗口M每个并计算〜N/MDFT,所以复杂度是NMMlogM=NlogM.
  3. Yule-Walker 方法AR(p)反转一个p×pToeplitz 矩阵的复杂度为O(p2). Toeplitz 矩阵本身是由计算自相关值的估计形成的p不同的滞后N- 给出复杂度的长数据向量O(Np). 所以总复杂度是O(p2+Np).
  4. 音乐估计M×M来自数据的自协方差矩阵,通过平均外部产品,例如,N/M M- 原版的长块N-long 信号,然后计算特征分解M×M自协方差矩阵是O(M3). (M由用户选择。对于一个信号p频率成分,选择M>p.)

我认为复杂性不受输入噪声水平的影响,因为算法本身不会改变。

参考

[*] J. Proakis 数字信号处理,第 4 版 Ch。14 https://www.pearsonhighered.com/program/Proakis-Digital-Signal-Processing-4th-Edition/PGM258227.html