1. 课程导入与频谱分析定位

在科研数据处理中,我们采集到的绝大多数信号都以“时间序列”的形式存在:示波器记录电压随时间的波动,加速度计输出振动幅值随时间的序列,脑电设备保存电极电位随时间的变化,光谱仪给出强度随波长(或频率)的分布。这些数据的第一种自然视角是时域视角——横轴是时间,纵轴是幅值,我们关心峰值何时出现、波形是否畸变、是否存在突变或漂移。时域视角回答的是“什么时候发生了什么”。

但很多科研问题并不止于“何时发生”。当我们需要判断一台旋转机械的振动主要来自哪个部件、一段语音信号的能量集中在哪个频带、一个传感器通道中是否混入了 50 Hz 工频干扰、某段脑电在 α 波段是否增强时,仅看时域波形往往难以直接得出结论。这时就需要切换到频域视角:把信号分解为不同频率成分的叠加,横轴变为频率,纵轴表示各频率成分的幅值、功率或相位。频域视角回答的是“信号由哪些频率成分构成、各自占多大比重”。

频谱分析正是实现这一视角转换的核心工具。它的本质是把时域信号通过傅里叶变换类方法映射到频域,得到频谱(幅值谱、功率谱、相位谱等),再据此识别主导频率、带宽、谐波结构、噪声底与干扰成分。需要强调的是,频谱分析不是“另一种画图方式”,而是一种特征提取与问题定位手段:它把难以直接比较的时域波形,转换为可量化、可对比、可判读的频率特征。

从数据采集到频谱解读的完整链路 ① 数据采集 传感器 / 采样率 fs 量程、抗混叠滤波 ② 预处理 去趋势 / 去均值 加窗、分段、去噪 ③ 频谱变换 FFT / PSD / STFT 频率分辨率 Δf ④ 特征提取 主峰频率、谐波 带宽、噪声底 ⑤ 解读 结论与 验证 时域视角:幅值 — 时间 t 幅值 波形叠加,难以直接分辨成分 频域视角:幅值 — 频率 f / Hz 幅值 f₁ 2f₁ f₂ 噪声 成分分离,主峰与谐波清晰可读 FFT

上图把整条链路拆成五个环节,并对比了同一信号在时域与频域中的呈现差异。理解这条链路,是后续所有章节的基础。下面逐环节说明关键要点。

采集环节决定频谱的“天花板”

频谱分析的质量上限在采集阶段就已确定。最关键的两个参数是采样率 fs 和采样时长 T。根据奈奎斯特采样定理,可无混叠分析的频率上限为 fs/2。例如以 fs = 1000 Hz 采集,则最高可分析到 500 Hz;若信号中存在 600 Hz 成分而抗混叠滤波不足,它会被折叠到 400 Hz 处,形成虚假谱峰——这类错误在频谱图上无法自行暴露,必须在采集端防范。

采样时长 T 则决定频率分辨率 Δf = 1/T。采集 1 s 数据,Δf = 1 Hz;采集 10 s,Δf = 0.1 Hz。若两个相邻谱峰相距 0.3 Hz,用 1 s 数据无法分辨,用 10 s 数据才能分开。因此“能不能看清”在采集时就已决定,后期无法凭空提高分辨率。

预处理不是可选项

原始序列常含直流偏置与缓慢漂移,直接做 FFT 会在 0 Hz 附近形成巨大谱峰,掩盖低频有效成分,因此通常先去均值、去线性趋势。若信号非整周期截断,频谱会出现频谱泄漏:单频成分的能量扩散到邻近频点,表现为谱峰变宽、旁瓣抬高。工程上常用汉宁窗(Hann)抑制泄漏,代价是主瓣变宽、幅值需按窗函数做增益修正。对非平稳信号(频率随时间变化),则需改用短时傅里叶变换(STFT)或小波分析,以“时间—频率”二维方式呈现。

从频谱到结论:看什么、怎么判

得到频谱后,判读通常围绕四类特征展开:主峰频率(能量最大处,对应主导振动或周期成分)、谐波结构(基频整数倍处的谱线,提示非线性或机械松动)、带宽与谱形(宽带抬升常对应随机噪声或冲击)、噪声底与信噪比(决定弱成分能否被可靠识别)。以一台转速 1500 r/min 的电机为例,其转频为 25 Hz;若频谱在 25 Hz 出现主峰、50 Hz 与 75 Hz 出现递减谐波,同时 100 Hz 处有独立窄峰,则可初步判断存在转频振动、谐波失真以及可能的 2 倍频电磁或结构共振成分,再结合工况与测点位置进一步验证。

判读频谱时务必同时记录采集参数(fs、T、窗函数、平均次数)。同一信号在不同参数下得到的谱形可能明显不同,脱离参数谈谱峰没有意义。

建立正确的分析习惯