噪声信号的 1/3 倍频程波段分析



我想对噪声信号进行1/3倍频程波段分析。我有一个代表声压(Pa(的时间信号。我有中心频率和下限和上限截止的阵列。最终目标是获得一个条形图,表示声级与 1/3 倍频程频段的关系。

根据我遵循的一个理论,建议首先做信号的FFT。然后在 1/3 倍频程频段重建信号。在每个构造中,计算压力的RMS值。就是这样。将这些RMS值转换为声级(dB(,并绘制与中心频率的关系图。

我面临的问题是如何使用IFFT函数重建信号,因为新信号的振幅点数量较少。因此,由于时间和振幅之间的维度不一致,重建基本上是不可能的。我被困在这一点上。

我在DSP方面的经验很少。因此,即使该方法与我上面尝试解释的方法不同,任何帮助也将不胜感激。提前谢谢。

要重建特定 1/3 倍频程频段的时域信号,请在执行 FFT 之前将频段外所有频率的幅度设置为零。 每个频段的 IFFT 大小与我们的原始 FFT 相同,因此您最终将获得每个频段输出相同的时间分辨率。

为每个频段做一个全尺寸的IFFT显然非常昂贵。 有一些方法可以优化这一点,但实际上可能没有理由首先重建时域信号

IFFT运算1不会改变信号的RMS值,因此只需使用频域分量计算即可,而完全忽略IFFT值。 由于RMS计算的工作方式,您甚至不需要记住频域幅度是复数 - 如果您以相同的方式处理实部和虚部,您将获得比例答案。 只是不要忘记在 RMS 除数中包含所有额外的零,因为每个波段都有不同数量的零。

1- 对于数学标准 IFFT 来说,这是正确的。 某些实现会将输出除以信号长度。 对于所有波段,这都是相同的,因此不会影响它们的相对值,并且易于校正。

最新更新