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,边带将无法分辨。

采样参数设计:从信号到频谱的约束链 原始信号 f_max(最高频率分量) B(信号带宽) T(采样时长需求) 抗混叠滤波器 截止频率 ≤ f_s/2 过渡带陡度要求 过采样可放宽 ADC 采样 f_s ≥ 2.56 × f_max N = f_s × T Δf = f_s / N = 1/T 频谱输出 0 ~ f_s/2 Δf 分辨率 无混叠 混叠风险:当 f_s 不足时 0 f_s/2(奈奎斯特频率) f_s 真实 f=800Hz 折叠 混叠 f=200Hz 旋转机械采样设计 • 转频 f_r → 需分辨边带 Δf ≤ 0.1×f_r • 轴承高频段 → f_s ≥ 2.56×f_max • 采样时长 T ≥ 1/Δf • 示例:f_r=25Hz,Δf=1Hz,T≥1s 瞬态事件采样设计 • 冲击持续时间 τ → f_s ≥ 10/τ • 预触发缓冲保留事件前数据 • 示例:τ=1ms → f_s ≥ 10kHz • 配合峰值保持或高速连续采集

旋转机械的采样参数设计

旋转机械的振动信号通常包含丰富的频率成分:轴转频及其谐波、齿轮啮合频率及其边带、轴承故障特征频率(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
工程提示:采样率并非越高越好。过高的采样率会带来数据量膨胀、存储成本上升和后续处理负担。在满足分析带宽和分辨率的前提下,选择适度的过采样倍数(2.56~5倍)即可。对于长期在线监测,可采用“高采样率触发采集 + 低采样率连续监测”的双模式策略。
常见误区:提高采样率不能改善频率分辨率。若发现边带无法分辨,应延长采样时间而非提高采样率。另外,抗混叠滤波器必须位于ADC之前,数字滤波无法消除已经发生的混叠。

采样参数设计是动态信号分析中最基础也最容易被忽视的环节。奈奎斯特准则给出了理论下限,而工程实践中需要综合考虑抗混叠滤波器特性、频率分辨率需求、数据存储能力和分析目标频段。对于旋转机械,重点在于覆盖轴承高频段的同时保证边带分辨率;对于瞬态事件,重点在于足够高的采样率捕捉冲击波形并利用预触发保留完整事件。掌握这些原则,才能在后续的频谱分析和故障诊断中获得可信的数据基础。