matlab - 径向功率谱分析

Question

首先 - 我需要根据波数 |k| 绘制这个径向平均光谱一维。但是我很早就计算了|k| 使用两个矩阵，所以我有一个表示|k|的矩阵。所以我也必须对这个矩阵进行径向平均？

第二 - 如果我在使用 FFT2 函数之前从我的数据中删除一个平均值，我预计只会看到值幅度的差异，但实际上我注意到径向平均频谱的形状存在差异。

我正在尝试使用它来分析磁场数据。谢谢你的帮助！！

Answer 1

您的第二个问题的解释来自 FFT 的加性属性：

FFT(c) = FFT(a+b) = FFT(a) + FFT(b)

如果b是一个偏移量 - 即平均值 - 您从 signal 中删除c，那么 的频谱c-b将等于您的原始频谱减去常数的频谱b，但常数的 FFT 会导致 sinc 摆动！因此，通过消除信号的平均值，您可以消除零频率的 sinc 毛刺。