3. 信号采样与频谱分析前提

在频谱分析中，我们最终看到的谱线位置、幅值乃至“是否存在某个频率成分”，往往并不完全由原始物理信号决定，而是被采样环节提前“锁定”了。采样是模拟信号进入数字分析世界的第一道关口，一旦这道关口设置不当，后续再精细的FFT、加窗或平均都无法恢复已经丢失或混淆的信息。本节围绕采样定理、奈奎斯特频率、混叠现象与抗混叠滤波四个相互关联的概念展开，并通过一个可复算的采样案例，说明采样参数如何决定频谱结果。

采样定理与奈奎斯特频率

对连续信号 x(t) 以采样间隔 Δt 进行理想采样，得到序列 x[n] = x(nΔt)，采样频率 f_s = 1/Δt。采样定理指出：若信号中最高频率成分不超过 f_max，则当 f_s > 2f_max 时，可由采样序列无失真地重建原信号。这里的 2f_max 是临界值，而 f_s/2 被称为奈奎斯特频率（折叠频率）。换言之，采样后频谱以 f_s 为周期重复，只有位于 [0, f_s/2] 区间内的频率才能被唯一表示。

工程上不能取等号。若 f_s = 2f_max，相位与幅值的重建对采样时刻极其敏感，实际系统总留有裕量。常见做法是让 f_s 达到感兴趣最高频率的 2.5～5 倍，并在采样前用模拟低通滤波器把 f_s/2 以上的能量压制到可接受水平。

混叠：高频“伪装”成低频

当信号中含有超过 f_s/2 的频率成分时，该成分不会消失，而是被折叠到 [0, f_s/2] 区间内，形成一个虚假的低频峰，这就是混叠。设真实频率为 f₀，采样频率为 f_s，则其混叠后的表观频率为：

f_alias = | f0 - k·fs |，取使结果落在 [0, fs/2] 的整数 k

例如 f_s = 1000 Hz，信号含 800 Hz 成分，则 k = 1 时表观频率为 |800 − 1000| = 200 Hz。频谱上会在 200 Hz 处出现一个本不存在的峰，而 800 Hz 处什么也没有。更麻烦的是，混叠后的频率与真实低频成分完全重合，仅凭频谱无法区分，因此混叠必须在采样前预防，而不能在采样后修正。

抗混叠滤波器的位置与指标

抗混叠滤波器必须是模拟低通滤波器，且位于采样器（ADC）之前。数字滤波器工作在采样之后，此时混叠已经发生，无能为力。设计时需关注三个指标：截止频率 f_c（通常取略低于 f_s/2）、过渡带宽度、阻带衰减。理想情况下希望在 f_s/2 处已有足够衰减；实际滤波器有过渡带，因此常采用“过采样 + 数字抽取”策略：先用较高频率采样，使模拟滤波器过渡带相对宽松，再用数字滤波器精确限带并降采样。

采样案例：同一信号、两种采样率

设一个教学示例信号由 120 Hz 与 800 Hz 两个正弦成分叠加，幅值均为 1。用两种采样率各采集 1 秒并做 FFT，观察结果差异（以下为按公式推算的示例结果，用于说明原理，非实测数据）。

方案 B 中，800 Hz 折叠为 |800 − 1000| = 200 Hz，与 120 Hz 成分并存，频谱上看到 120 Hz 和 200 Hz 两个峰。若分析者不知道原始信号含 800 Hz，很可能把 200 Hz 误判为真实成分，进而得出错误的物理结论。这正是采样参数“决定”频谱结果的典型体现。

实操中的参数选择与检查

面对一个待分析的信号，可按以下思路确定采样参数：先估计或实测信号中最高有效频率 f_max（包括可能存在的噪声与干扰高频）；据此选 f_s ≥ 2.5f_max，若关心幅值精度可提高到 5 倍；确认 ADC 前端模拟低通截止频率低于 f_s/2，并核对过渡带与阻带衰减是否足以压制 f_s/2 以上的能量。采集后还应做一次自检：改变采样率重新采集，若某谱峰位置随 f_s 改变而移动，则该峰很可能是混叠产物。

此外，采样长度决定频率分辨率 Δf = f_s/N（N 为采样点数），它与采样率、分析带宽共同构成一组需要权衡的参数。提高 f_s 扩展了可分析带宽，但在 N 固定时降低了分辨率；反之降低 f_s 可提高分辨率，却压缩了无混叠带宽。理解这一制约关系，才能在具体分析任务中做出合理取舍，而不是盲目追求高采样率或长采样时间。