4、信号预处理：去趋势、去均值、滤波（低通/高通/带通）、异常值处理、重采样与数据清洗

振动信号在采集过程中不可避免地会混入直流偏置、趋势项、环境噪声、瞬态干扰以及采样率不匹配等问题。若未经预处理直接进入特征提取与故障诊断环节，轻则降低模型灵敏度，重则导致误判漏判。本章围绕去趋势、去均值、滤波、异常值处理、重采样与数据清洗六个核心环节，系统讲解其原理、操作流程、工程案例与常见误区，帮助学员建立一套稳健的信号预处理流水线。

4.1 去趋势与去均值：消除基线漂移与直流分量

核心原理：振动传感器（如压电式加速度计）在低频段容易受到温度漂移、前置放大器零点漂移以及积分电路饱和的影响，导致信号中混入缓慢变化的趋势项。趋势项在频域表现为极低频分量，会掩盖真实的低频故障特征（如转频及其谐波）。去均值则是消除信号的直流分量，使信号围绕零值波动，便于后续FFT分析和特征统计。

判断方法与操作流程：

1. 趋势项检测：对原始信号进行最小二乘线性拟合，若拟合斜率显著（R²>0.3），则判定存在线性趋势；若残差仍呈现缓慢波动，需考虑高阶多项式拟合（2~3阶）。
2. 去趋势操作：使用detrend函数（Python/Matlab均支持），指定type='linear'或type='constant'；对于非线性趋势，可先进行高通滤波（截止频率设为最低关注频率的1/3）。
3. 去均值：计算信号均值并减去，公式为 x_clean = x - mean(x)。注意：去均值应在去趋势之后进行，否则趋势项会残留部分直流分量。

4.2 滤波处理：低通、高通、带通与带阻

核心原理：滤波的本质是频率选择。低通滤波保留低于截止频率的成分，用于去除高频噪声；高通滤波保留高于截止频率的成分，用于去除低频漂移；带通滤波保留特定频段，常用于提取轴承故障特征频率（如BPFO、BPFI）；带阻滤波用于抑制特定干扰（如50Hz工频）。

滤波器类型对比：

滤波器类型	典型应用	截止频率设置	注意事项
低通（Butterworth）	去除高频噪声，保留低频振动	设为最高关注频率的1.5~2倍	阶数过高会引起相位滞后
高通（Chebyshev I）	去除趋势项和直流偏置	设为最低关注频率的1/3	通带纹波可能影响幅值精度
带通（FIR）	提取特定故障特征频段	中心频率±10%~20%带宽	需保证过渡带陡峭
带阻（IIR）	抑制工频及谐波干扰	阻带宽度设为干扰频率±2Hz	注意阻带引起的相位突变

操作流程：

1. 根据故障特征频率范围确定滤波器类型与截止频率。
2. 使用scipy.signal.butter设计滤波器系数，filtfilt实现零相位滤波（避免相位偏移）。
3. 验证滤波效果：对比滤波前后频谱，确认目标频段保留、干扰频段衰减≥40dB。

4.3 异常值处理：识别与修复数据中的离群点

核心原理：异常值（Outliers）可能来自传感器瞬时掉电、电磁干扰脉冲、电缆接触不良或数据采集卡溢出。这些离群点会严重扭曲均方根值、峰值因子等时域指标，并在频谱中产生虚假的宽带噪声。

判断方法：

• 3σ准则：计算信号均值μ和标准差σ，超出μ±3σ范围的点视为异常值。适用于信号近似正态分布的场景。
• Hampel滤波：基于滑动窗口的中值偏差（MAD），对局部突变更鲁棒。窗口长度通常设为信号长度的1%~5%。
• 差分阈值法：计算相邻采样点差值，若差值超过设定阈值（如最大正常差值的5倍），则判定为异常。

处理策略：

1. 剔除：直接删除异常点，但需保证删除后信号连续，否则需插值。
2. 替换：用窗口内中值或线性插值替换异常值，保持信号长度不变。
3. 标记：在数据清洗阶段标记异常段，供后续特征提取时排除。

4.4 重采样与数据清洗：统一时间基准与提升数据质量

核心原理：多通道振动采集系统可能因采样率不一致、触发延迟或时钟漂移导致数据时间轴不对齐。重采样通过插值算法将信号统一到目标采样率，保证后续多通道联合分析（如全息谱、阵列信号处理）的准确性。数据清洗则是对整体数据质量进行系统性校验与修复。

重采样操作流程：

1. 确定目标采样率：通常取各通道采样率的最小公倍数，或统一到最高采样率。
2. 抗混叠滤波：重采样前先进行低通滤波，截止频率设为目标采样率的0.4倍（满足奈奎斯特准则）。
3. 插值重采样：使用scipy.signal.resample（FFT法）或resample_poly（多相滤波法），后者计算效率更高。

数据清洗清单：

检查项	判定标准	处理方式
数据完整性	缺失率 < 0.1%	线性插值或删除该段
重复时间戳	无重复	保留首条，删除后续
时间戳单调性	严格递增	排序并标记异常段
信号幅值范围	在传感器量程内	超限点按异常值处理
通道间延迟	延迟 < 1个采样周期	互相关法校正

实战场景：某旋转机械试验台使用4个加速度计，其中2个通道采样率为51.2kHz，另2个为25.6kHz。为进行全息谱分析，需将所有通道统一至51.2kHz。先对25.6kHz通道进行抗混叠滤波（截止频率=20kHz），再使用多相滤波器进行2倍插值，重采样后各通道频谱一致性误差小于0.5dB。

4.5 预处理流程的工程化落地与质量验证

核心原理：预处理不是孤立的步骤，而是需要形成标准化流水线，并建立质量验证机制。建议将预处理流程封装为可复用的函数或类，输入原始信号，输出清洗后的信号及质量指标。

推荐流程：

1. 原始信号导入 → 2. 去趋势（线性/多项式） → 3. 去均值 → 4. 带通滤波（根据故障类型） → 5. 异常值检测与修复 → 6. 重采样（如需） → 7. 数据清洗校验 → 8. 输出信号与质量报告。

质量验证指标：

• 信噪比（SNR）：预处理后SNR应提升至少10dB。
• 趋势项残留：对预处理后信号再次拟合，斜率应接近零。
• 频谱保真度：对比预处理前后已知故障特征频率的幅值变化，衰减应小于3dB。
• 相位失真：使用零相位滤波（filtfilt）确保相位无偏移。