径向功率谱分析

Posted 2023-02-25

【中文标题】径向功率谱分析

【英文标题】：Radially Power Spectrum Analysis

【发布时间】：2013-08-04 15:17:30

【问题描述】：

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

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

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

【问题讨论】：

您可能想在物理（或数学？）stackexchange 网站上发布此内容。我不太确定第一点，但更多信息可能会有所帮助，例如 |k|不是只是半径的函数吗？

一个更普遍的问题是：您期望（或观察到）图像中的哪种对称性？是给你|k|的操作吗？生成具有相同对称性的矩阵？

【参考方案1】：

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

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

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