4. 信号与数据预处理
在航空航天电子系统中,传感器输出的原始信号几乎不可能直接用于故障检测与诊断。无论是惯性测量单元(IMU)、大气数据计算机、GPS接收机,还是发动机振动传感器、温度/压力变送器,其原始数据都混杂着电子噪声、工频干扰、量化误差、野值以及不同步的采样时钟。如果把这些“脏数据”直接送入诊断算法,轻则虚警率飙升,重则漏检真实故障。因此,信号与数据预处理是故障诊断链路中决定成败的第一道工序。本节围绕噪声抑制、滤波、归一化、异常值处理、时间同步与多源数据对齐六个环节,说明它们各自解决什么问题、实际怎么做,以及如何判断处理是否到位。
一、噪声抑制与滤波:先分清“噪声”和“信号”
噪声抑制的目标不是把信号变得“光滑好看”,而是在保留故障特征频带的前提下压低无用成分。工程上首先要区分两类噪声:一类是宽带随机噪声,如热噪声、散粒噪声,通常用低通滤波或滑动平均抑制;另一类是窄带干扰,如50 Hz工频及其谐波、开关电源的固定频率纹波,用陷波器或自适应滤波更有效。
以某型无人机IMU的俯仰角速度信号为例,采样率 fs = 200 Hz,有用姿态运动带宽约 0~10 Hz,但信号中叠加了约 50 Hz 的电机振动干扰。若直接使用一阶低通滤波,截止频率设为 10 Hz,其离散化实现为:
alpha = dt / (RC + dt) # dt = 1/fs = 0.005 s
y[n] = alpha * x[n] + (1 - alpha) * y[n-1]
# 取 RC = 1/(2*pi*10) ≈ 0.0159 s
# alpha ≈ 0.005 / (0.0159 + 0.005) ≈ 0.239
这样50 Hz干扰会被显著衰减,但一阶滤波在10 Hz附近仍有明显相位滞后,对姿态控制不利。更实用的做法是采用二阶巴特沃斯低通,截止频率10 Hz,通带内平坦、相位滞后可控;若50 Hz干扰特别强,再级联一个中心频率50 Hz、带宽2 Hz的IIR陷波器。判断滤波是否合适,不能只看时域曲线是否平滑,而应观察频谱:有用频带内幅值衰减应小于1 dB,干扰频点衰减应大于20 dB,同时阶跃响应超调不宜过大。
二、归一化:让不同量纲的信号可比
航空航天电子系统中,同一诊断模型往往同时接收温度(℃)、电压(V)、转速(r/min)、振动加速度(g)等量纲差异巨大的信号。若不做归一化,数值范围大的通道会主导距离计算或梯度更新。常用方法有两种:最小-最大归一化把数据映射到[0,1],适合边界已知且无极端野值的场景;Z-score标准化把数据变换为均值0、标准差1,适合近似高斯分布且存在少量离群点的场景。
例如某发动机排气温度通道,正常巡航段样本均值 μ = 620 ℃,标准差 σ = 18 ℃。某时刻测得 685 ℃,则Z-score为 (685-620)/18 ≈ 3.61,已经超过3σ,应触发异常值复核。这里要注意:归一化参数必须用训练集或健康基线数据统计得到,并固化到诊断系统中;若每帧都用当前数据重新计算均值和方差,故障引起的均值漂移会被“归一化掉”,反而掩盖故障。
三、异常值处理:先识别,再决定剔除还是修正
异常值来源包括传感器瞬时失效、通信误码、电磁脉冲干扰等。处理原则是:不能无条件删除,也不能无条件保留。工程上常用“3σ准则+中值滤波+物理限幅”三级策略。以某气压高度通道为例,采样率10 Hz,连续5帧高度跳变超过300 m且超出飞机最大爬升率对应范围,则判为野值。
识别后有三种处理方式:一是直接剔除该帧,适用于孤立野值;二是用前值保持或线性插值替代,适用于短时丢包;三是用中值滤波输出替代,适用于连续多帧受扰。若某通道连续超过0.5 s均为异常,则不应继续“修补”,而应将该通道标记为失效并切换到余度通道,同时向诊断系统上报“传感器故障”而非“数据异常”。这一区分非常关键:数据预处理只能处理“数据问题”,不能替代故障诊断本身。
四、时间同步与多源数据对齐:把不同时钟的数据放到同一时间轴
多源数据对齐是航空航天电子诊断中最容易被低估的环节。IMU可能以200 Hz输出,GPS以10 Hz输出,发动机控制器以50 Hz输出,且各设备时钟存在晶振漂移。若直接把不同时刻的数据拼成一个特征向量,诊断模型看到的其实是“时间错位”的伪相关。
实际做法分三步。第一步,统一时间基准:所有数据打上同一系统时间戳,通常以PPS(秒脉冲)或IEEE 1588精密时间协议对齐到主时钟。第二步,重采样到统一频率:对低频信号采用线性插值或样条插值升采样,对高频信号采用抗混叠滤波后抽取降采样。第三步,估计并补偿固定时延:例如GPS位置相对IMU存在约80 ms传输延迟,需在融合前把GPS数据时间戳前移80 ms。
判断对齐是否成功,可用一个简单检验:飞机匀速平飞时,IMU积分得到的高度变化率应与气压高度变化率、GPS高度变化率趋势一致;若三者之间存在固定时间差,则说明时延补偿未做好。另一个检验是互相关函数峰值位置,理想情况下应位于零延迟附近。
五、工程落地中的顺序与常见陷阱
上述五个环节并非随意排列。推荐顺序是:先滤波去噪,再异常值处理,再归一化,最后时间同步与多源对齐。原因在于:若先归一化再滤波,野值会被归一化放大;若先对齐再滤波,不同采样率信号在插值过程中会把噪声扩散到相邻时刻。另一个常见陷阱是“用故障数据统计归一化参数”,这会让故障样本被压缩到正常范围内,直接导致漏检。因此归一化参数必须来自健康基线,并在系统运行中定期用健康数据更新,但更新频率不宜过高。
最后需要强调:预处理的目标是“保真地保留故障特征”,而不是“让数据看起来漂亮”。过度滤波会抹掉轴承早期故障的高频冲击成分,过度平滑会掩盖传感器间歇性失效的跳变特征。工程上应通过频谱对比、残差分析和诊断算法闭环验证来确认预处理参数是否合适。只有预处理链路稳定可靠,后续的特征提取、故障分类与健康评估才有意义。