3. 采样定理与工程化采样策略
在动态信号分析中,采样率的选择是决定测量成败的第一道关口。无论后续的FFT、包络解调还是阶次分析多么精密,一旦采样环节出现混叠或分辨率不足,数据中携带的故障特征就已经被不可逆地破坏。本节从奈奎斯特准则出发,逐步展开过采样、欠采样风险、采样率与频率分辨率的权衡,并给出针对旋转机械与瞬态事件的采样参数设计方法。
奈奎斯特准则与混叠的物理含义
奈奎斯特采样定理指出:若连续信号中最高频率分量为 fmax,则采样率 fs 必须满足 fs > 2fmax,才能由离散序列无失真地重建原信号。2fmax 称为奈奎斯特率,fs/2 称为奈奎斯特频率或折叠频率。当输入信号中含有高于 fs/2 的频率成分时,这些成分不会消失,而是被“折叠”到 0~fs/2 区间内,形成混叠频率:
f_alias = | f_signal − k · f_s |,k 取使结果落入 [0, f_s/2] 的整数
例如采样率 fs = 1000 Hz,奈奎斯特频率为 500 Hz。若信号中存在 800 Hz 成分,则 k=1 时混叠频率为 |800−1000| = 200 Hz。在频谱上,这个 200 Hz 的虚假峰值与真实低频成分无法区分,极易被误判为故障特征频率。因此,抗混叠滤波器必须在ADC之前将高于 fs/2 的成分衰减到噪声底以下,这是硬件层面不可省略的环节。
过采样与欠采样:两种策略的工程取舍
过采样是指实际采样率显著高于奈奎斯特率,通常取 fs = (2.5~10)×fmax。过采样的收益包括:降低对抗混叠滤波器过渡带陡度的要求、提高信噪比(每倍频程约改善3 dB)、为后续数字滤波和抽取留出余量。代价是数据量增大、存储与传输压力上升。在旋转机械监测中,常用过采样配合数字抽取,先以高采样率采集,再按分析需求降采样。
欠采样则指采样率低于奈奎斯特率,通常被视为错误操作。但在特定场景下,欠采样可被有意利用:当只关心某个窄带高频载波附近的调制信息时,通过带通采样将高频段搬移到低频,从而用较低的采样率获取该频段信息。带通采样的条件是信号带宽 B 满足 fs > 2B,且载波频率与采样率满足特定整数关系。这种技术对滤波器要求极高,工程中需谨慎使用。
采样率与频率分辨率的权衡
频率分辨率 Δf 由采样率 fs 和采样点数 N 共同决定:Δf = fs/N。同时,采样时长 T = N/fs,因此 Δf = 1/T。这意味着:要提高频率分辨率,必须延长采样时间,而不是单纯提高采样率。提高采样率只会扩展分析带宽,不会改善分辨率。
这一关系对旋转机械诊断至关重要。假设某齿轮箱输入轴转频为 25 Hz,需要分辨边带间隔 1 Hz 以识别调制现象,则采样时长至少需要 1 秒。若同时希望覆盖 5 kHz 的轴承高频段,采样率至少取 12.8 kHz(2.56倍过采样)。此时若采样 1 秒,N = 12800 点,Δf = 1 Hz,满足要求。若只采样 0.2 秒,Δf = 5 Hz,边带将无法分辨。
旋转机械的采样参数设计
旋转机械的振动信号通常包含丰富的频率成分:轴转频及其谐波、齿轮啮合频率及其边带、轴承故障特征频率(BPFO、BPFI、BSF、FTF)以及结构共振。设计采样参数时,需要先明确分析目标频段。若关注轴承早期故障,需覆盖 5~10 kHz 的高频共振区;若关注齿轮调制边带,则需保证足够的频率分辨率。
以某电机驱动齿轮箱为例,输入轴转频 fr = 25 Hz,齿轮啮合频率 fm = 500 Hz,轴承故障特征频率最高约 3.2 kHz。分析目标包括:啮合频率周围的调制边带(间隔约 25 Hz)和轴承高频共振(中心约 8 kHz)。
第一步确定分析带宽:轴承共振中心 8 kHz,考虑调制边带,最高分析频率取 10 kHz。第二步确定采样率:fs ≥ 2.56 × 10 kHz = 25.6 kHz,取 25.6 kHz 或 51.2 kHz。第三步确定频率分辨率:要分辨 25 Hz 的边带,Δf 至少取 25/4 ≈ 6.25 Hz,工程上取 Δf = 1 Hz 更稳妥。第四步确定采样时长:T = 1/Δf = 1 秒。第五步确定采样点数:N = fs × T = 25600 点。若设备内存有限,可分段采集并做平均,但每段仍需满足 T ≥ 1 秒。
瞬态事件的采样策略
瞬态事件(如轴承局部剥落产生的冲击、齿轮断齿、碰摩)持续时间短、频带宽,对采样率的要求与稳态信号不同。冲击脉冲的上升时间 τ 决定了其频谱宽度,近似有 fmax ≈ 1/(πτ)。工程上通常要求采样率至少为冲击带宽的 5~10 倍,即 fs ≥ 10/τ。
例如,某轴承滚道剥落产生的冲击上升时间约 0.5 ms,则 fmax ≈ 1/(π×0.0005) ≈ 637 Hz,但实际冲击能量常延伸至 10 kHz 以上,因此 fs 取 100 kHz 以上才能完整捕捉冲击波形。同时,瞬态采集需要预触发功能:在触发条件满足前持续缓存数据,触发后保留触发点之前的一段数据,确保冲击起始沿不被丢失。触发方式可选择幅值触发、包络触发或转速触发。
工程化采样参数设计清单
综合以上分析,针对旋转机械与瞬态事件的采样参数设计可归纳为以下可执行步骤:
| 步骤 | 操作内容 | 关键公式/判据 | 示例值 |
|---|---|---|---|
| 1 | 确定分析目标频段 | 根据故障类型选择 f_max | 轴承共振 10 kHz |
| 2 | 确定采样率 | f_s ≥ 2.56 × f_max | 25.6 kHz |
| 3 | 确定频率分辨率 | Δf ≤ 边带间隔 / 4 | 1 Hz |
| 4 | 确定采样时长 | T = 1 / Δf | 1 s |
| 5 | 确定采样点数 | N = f_s × T | 25600 |
| 6 | 配置抗混叠滤波 | 截止频率 ≤ f_s/2 | 12.8 kHz |
| 7 | 瞬态事件附加 | f_s ≥ 10/τ,启用预触发 | τ=0.5ms → 100 kHz |
采样参数设计是动态信号分析中最基础也最容易被忽视的环节。奈奎斯特准则给出了理论下限,而工程实践中需要综合考虑抗混叠滤波器特性、频率分辨率需求、数据存储能力和分析目标频段。对于旋转机械,重点在于覆盖轴承高频段的同时保证边带分辨率;对于瞬态事件,重点在于足够高的采样率捕捉冲击波形并利用预触发保留完整事件。掌握这些原则,才能在后续的频谱分析和故障诊断中获得可信的数据基础。