3. 特征工程与健康指标构建
在高端装备健康评估中，传感器采集到的原始振动、电流、温度、压力等信号，往往包含大量与退化无关的成分：工况波动、负载变化、环境噪声、装配差异。若直接把原始波形送入模型，模型学到的很可能是“今天转速多少”而不是“轴承还剩多少寿命”。特征工程的任务，就是把原始信号转换成对退化敏感、对工况鲁棒的数值描述，并进一步压缩成单一或少数几个健康指标（Health Indicator, HI），让退化轨迹可观察、可比较、可报警。
一、时域特征：最直接但需区分量纲
时域特征直接对采样序列 x[1..N] 计算，物理意义清晰、计算量小，适合在线监测。常用量包括：
- 有量纲指标：均值、峰值、峰峰值、均方根值（RMS）、方差。RMS 反映能量水平，对磨损引起的冲击能量上升较敏感，但受负载影响大。
- 无量纲指标：峭度（Kurtosis）、偏度、峰值因子（Crest Factor）、波形因子、裕度因子。它们对幅值变化不敏感，更适合变工况场景。峭度对早期冲击性故障尤其敏感，但接近失效时反而可能下降。
以某型轴承试验为例，采样率 20 kHz，每 10 分钟取 1 秒数据。正常段 RMS 约 0.35 g，峭度约 3.0；进入早期点蚀后，RMS 缓慢升至 0.5 g，而峭度先升至 5～6，随后在严重剥落阶段回落至 3.5 左右。这说明单一指标难以覆盖全寿命，需要组合使用。
二、频域特征：定位故障频率成分
频域特征通过 FFT 或功率谱密度（PSD）获得，核心价值在于把故障与特定频率对应起来。对旋转机械，需先计算转频 fr 及其倍频，再计算轴承外圈、内圈、滚动体、保持架的故障特征频率（BPFO、BPFI、BSF、FTF）。常用频域特征包括：
- 各阶转频幅值及其谐波能量占比；
- 故障特征频率处的幅值、边频带宽度；
- 谱重心、谱熵、频带能量比（如 0～1 kHz、1～5 kHz、5～10 kHz 分段能量）。
例如，某齿轮箱在 1500 r/min（fr=25 Hz）下运行，若在 125 Hz（5 倍频）附近出现明显边频，且边频间隔等于 25 Hz，通常提示该级齿轮存在局部损伤。频域特征的优势是物理可解释，缺点是对转速波动和分辨率敏感，需配合阶次分析或重采样。
三、时频域特征：兼顾非平稳与瞬态冲击
高端装备常处于变转速、变负载工况，单纯时域或频域都不够。时频域方法把信号映射到“时间—频率”平面，常用手段包括短时傅里叶变换（STFT）、小波变换（WT）、经验模态分解（EMD）及其变体。可提取的特征有：
- 小波包各频带能量比；
- IMF 分量的能量、样本熵、排列熵；
- 时频图的纹理特征、瞬时频率均值与方差。
样本熵和排列熵对冲击成分的出现非常敏感，适合刻画退化过程中的复杂度变化。实际中常把时频特征与工况参数（转速、负载）一起归一化，避免把工况变化误判为退化。
四、特征选择与降维：从几十维到可用维度
提取出几十甚至上百个特征后，并非越多越好。冗余特征会放大噪声、增加过拟合风险。工程上通常分两步：
- 过滤式筛选：计算每个特征与退化时间或剩余寿命的单调性、趋势性、相关性，剔除方差过小或与工况强相关的特征。常用指标包括 Spearman 相关系数、单调性指标 Monotonicity = |正负导数符号一致的比例之差|。
- 降维融合：对保留特征做主成分分析（PCA）、核主成分分析或自编码器压缩。PCA 前需标准化，否则量纲大的 RMS 会主导主成分。若希望保留物理可解释性，可改用加权融合而非黑箱降维。
例如，从 42 个候选特征中，先按单调性 > 0.6、与转速相关性 < 0.3 筛出 9 个，再用 PCA 取前 2 个主成分，累计方差贡献率约 87%，即可作为后续 HI 的输入。
五、健康指标构造准则与验证
一个合格的 HI 应满足三条基本准则：
- 单调性：随退化进程总体单调上升或下降，避免频繁反复。可用单调性指标量化，越接近 1 越好。
- 趋势性：与时间或寿命之间存在稳定趋势关系，便于外推预测。
- 鲁棒性：在相同退化阶段、不同工况或噪声水平下，HI 波动小，不因负载变化而误报。
构造方法上，常见有：加权融合（对筛选后的特征赋权求和）、马氏距离（以健康基线为参考，计算当前特征向量到基线的距离）、自编码器重构误差（健康数据训练，退化时重构误差增大）。加权融合可解释性强，适合工程落地；马氏距离能考虑特征间相关性；自编码器适合高维非线性场景，但需要足够健康数据。
验证时，建议用全寿命数据回放：画出 HI 随时间的曲线，检查是否单调、是否在失效前明显抬升、是否在早期工况波动段保持平稳。若 HI 在健康段就大幅波动，说明鲁棒性不足，应回退检查特征是否混入工况量，或调整归一化方式。只有通过这三项验证的 HI，才适合送入退化建模与剩余寿命预测环节。