快速傅里叶变换(FFT)计算频谱分布,避免遗漏特征频谱的处理方式

为了“尽量不漏掉特征频谱”,关键不是只把 FFT 算出来,而是把采样→预处理→谱估计→时频补充→特征提取/验证这一整条链路做得更鲁棒。下面按实战常用程度给出一套“防漏特征”的处理方式清单。

一、从源头避免“看不见”:采样与频率分辨率设计

1)采样率与抗混叠(anti-aliasing)

  • 采样率要保证:感兴趣最高频率 fmaxf_{max}​ 远小于 fs/2f_s/2,并留出余量(尤其有冲击/尖峰时,高频能量会很宽)。
  • 必须有模拟抗混叠低通(或设备内置同等能力),否则高频折叠到低频,容易把“假特征”当真特征,也会掩盖真特征。

2)频率分辨率 Δf

  • Δf=fs/N=1/T\Delta f = f_s/N = 1/T(T 是采样时长)。要分开很近的谱线,就要更长的记录长度 T。
  • 常见“漏特征”原因之一:分辨率不够导致两个特征峰粘在一起,或峰被噪声底淹没。

3)多档位采集(宽带扫一遍 + 重点放大)

  • 先用较高采样率做宽带普查(不漏大范围),再对可疑频段做更长时长/更高分辨率的“精扫”。

二、预处理:减少泄漏与低频污染,提升可检出性

4)去直流/去趋势(de-mean / detrend)

  • 直流分量、慢漂移会抬高低频噪声底,掩盖低频特征或制造“假低频峰”。

5)合理窗函数,控制谱泄漏

  • 特征频率不一定正好落在 FFT 的频点上,会发生谱泄漏scalloping loss(栅栏效应)
  • 经验:
    • Hann/Hamming:通用,泄漏较小,幅值损失可接受;
    • Blackman / Blackman-Harris:更强抑制旁瓣,适合强弱峰共存(高动态范围);
    • Flat-top:幅值测量更准,但主瓣更宽(分辨率变差),适合“要精确幅值”而不是“分开近峰”。

6)幅值校正

  • 加窗后要做coherent gain / ENBW相关校正,否则峰值幅度会系统性偏小,弱峰更容易“看起来不存在”。

三、谱估计更稳:避免“偶然性”与“噪声底掩盖”

7)Welch 分段平均(强烈推荐)

  • 把信号分段、加窗、重叠(如 50%~75%)、求功率谱再平均:
    • 优点:显著降低方差、稳定噪声底,弱峰更容易冒出来;
    • 代价:频率分辨率略降(每段更短)。

8)多锥窗 Multi-taper(更“防漏”)

  • 对窄带谱线+噪声场景特别好:用多组正交窗得到多份谱再平均,泄漏更可控、估计方差更低。
  • 在“强峰旁边找弱峰”时非常有价值。

9)时域同步平均 / 转速同步(若存在周期性)

  • 机械系统里很多特征与转速、工频、啮合频相关:
    • 若转速稳定:做Time Synchronous Averaging (TSA),能把随机噪声压下去;
    • 若转速波动:做Order Tracking(阶次分析),把频率轴换成“阶次”,否则峰会被“拉宽”看不见。

四、对准频点:减少“峰值落在两个 bin 之间”导致的漏检

10)增加记录长度(最有效) vs 零填充(zero-padding)

  • 增加记录长度:真正提高分辨率(更不漏)。
  • 零填充:不提高真实分辨率,但能让谱看起来更平滑、便于峰值定位(尤其配合峰值插值)。

11)峰值插值(bin 间插值)

  • 对检测谱线频率很有用:在峰附近做抛物线/高斯插值,可显著降低“栅栏效应”带来的低估与漏检风险。

12)Zoom-FFT / CZT(Chirp-Z Transform)

  • 如果你只关心某一小段频带:用 CZT 或 Zoom-FFT 直接对该频带做高分辨率分析,比“全频段超大 N”更高效、更不漏细节。

五、非平稳/瞬态特征:单次 FFT 最容易漏掉的部分

很多故障特征不是“全程稳定一条线”,而是间歇出现、频率漂移、冲击调制。这时“整段 FFT”会把能量摊平,导致看不到峰。

13)STFT(短时傅里叶)

  • 用滑动窗口看随时间变化的谱(谱图)。
  • 对“门机动作过程”“启停瞬态”“偶发冲击”特别有效。

14)小波/连续小波 CWT

  • 对冲击类、非平稳信号更敏感,时频分辨率自适应;常用来避免漏掉“短促但关键”的异常能量带。

15)包络解调(Envelope / Demodulation)

  • 轴承、齿轮等常见“高频共振被低频调制”:直接看原始 FFT 可能只看到一坨宽带。
  • 做带通(抓共振带)→希尔伯特包络→对包络做 FFT,经常能把“隐藏的调制频率”拉出来。

六、让“弱特征”不被强特征淹没:动态范围与去干扰

16)对数坐标/分段显示

  • PSD 用 dB 显示;或者同一张图做“全频 + 局部放大”,不然强峰会把弱峰压扁到看不见。

17)去已知干扰(工频/谐波/固定噪声源)

  • 例如 50/60 Hz 及其谐波、PWM 载波等:可用陷波/同步抑制,避免它们的泄漏旁瓣淹没邻近弱峰。

18)预白化/高通(谨慎用)

  • 对严重 1/f 噪声或低频漂移,可适度高通提高可检出性,但要确保不会把你真正关心的低频特征削掉。

七、实战建议:一套“防漏”默认配置(可直接套用)

如果你现在就要一个稳妥的默认流程(尤其适合振动/电流这类工业信号):

  1. 去直流 + detrend
  2. Hann 或 Blackman-Harris 加窗(强弱峰共存选后者)
  3. Welch:50%~75% 重叠,段长按“想要的 Δf\Delta f”定
  4. dB 显示 + 局部频段放大
  5. 峰值检测 + 峰值插值
  6. 对动作过程/非平稳:补一张 STFT 谱图
  7. 怀疑调制:做包络谱(带通→包络→FFT)
  8. 转速波动:上 Order Tracking

发表评论

您的邮箱地址不会被公开。 必填项已用 * 标注