4. 记录数据预处理实战

飞控故障定位的准确性,很大程度上取决于进入分析环节的数据是否“干净、同轴、同尺度”。原始飞行数据记录器(FDR/QAR)或遥测链路落盘的数据,往往来自多个采集通道、不同采样率、不同时间基准,甚至不同物理单位。若直接把这些数据送入频谱分析、残差比较或故障树推理,很容易把“数据问题”误判为“飞控问题”。本节围绕五个关键动作展开:时间对齐、单位换算、野值剔除、缺失值处理与数据质量评估,并给出可复现的处理流程与判断依据。

一、时间对齐:先统一时间轴,再谈相关性

时间对齐的目标是让不同参数在同一时间坐标下可比较。常见问题有三类:一是采样率不同,例如姿态角100 Hz、大气数据25 Hz、发动机参数10 Hz;二是起始时刻不同,各通道上电或记录触发存在毫秒级到秒级偏移;三是时钟漂移,长航时记录中晶振误差会累积。工程上通常以主飞控计算机的帧计数或GPS秒脉冲(PPS)作为基准时间轴,将其他通道重采样到统一时间栅格。

对齐时先做“粗对齐”:用GPS时间或记录器起始标志把各通道拉到同一绝对时间;再做“精对齐”:选取一个在多个通道中都能观察到的公共事件(如起落架收放、舵面阶跃、发动机转速突变),计算互相关峰值对应的时延,对通道做平移修正。对于采样率差异,推荐使用线性插值或零阶保持:连续量(高度、速度、姿态)用线性插值,离散状态量(开关量、档位)用零阶保持,避免插值出“中间态”。

对齐后应保留一份“时间轴映射表”,记录每个原始通道的采样率、起始偏移、插值方式和修正量,便于后续复现与审计。

二、单位换算:把物理量拉回同一量纲

单位不统一是故障定位中最隐蔽的陷阱之一。例如俯仰角可能以弧度记录,而判据以度为单位;空速可能以节(kt)记录,而模型以米/秒(m/s)计算;高度可能以英尺(ft)记录,而气压高度以米(m)参与比较。换算必须依据记录格式说明(如ARINC 429标签、IRIG 106 Chapter 10参数表)中的LSB(最低有效位)和偏移量进行,而不是凭经验猜测。

典型换算关系包括:角度 rad→deg 乘以 57.2958;速度 kt→m/s 乘以 0.514444;高度 ft→m 乘以 0.3048;温度 ℃→K 加 273.15;角速率 deg/s→rad/s 乘以 0.0174533。对于带偏移的工程量,使用 物理量 = 原始码值 × LSB + 偏移。换算完成后,建议对关键参数做“量级检查”:俯仰角应在±30°量级、空速在0~300 m/s量级、高度在0~12000 m量级,若出现数量级异常,优先怀疑单位或LSB错误。

三、野值剔除:区分“真故障”与“假跳变”

野值(outlier)是单点或短时的大幅跳变,可能来自电磁干扰、总线误码、传感器瞬时失效。剔除野值的原则是:不能把真实的故障响应误删。因此推荐“先检测、后判定、再处理”的三步法。检测可用三种互补方法:一是物理范围门限,如俯仰角超过±90°直接判为无效;二是变化率门限,如角速率超过±200°/s且持续小于2个采样周期;三是统计方法,如对滑动窗口计算中位数绝对偏差(MAD),当 |x - median| > k × MAD(k取3~5)时标记为可疑。

标记后不要立即删除,而是结合上下文判定:若该点前后数据连续、且与相关通道(如俯仰角与俯仰角速率、升降舵偏度)不一致,则判为野值;若多个相关通道同时出现同向跳变,则更可能是真实故障或真实机动,应保留并进入故障分析。处理方式优先选择线性插值或前值保持,并在数据质量标记中记录被修正的样本位置与数量。

四、缺失值处理:短缺口插值,长缺口标记

缺失值通常表现为记录中断、参数未激活或校验失败。处理策略按缺口长度区分:短缺口(一般小于0.5 s,约几十个采样点)可用线性插值或样条插值补齐;中等缺口(0.5~5 s)建议用前值保持并附加“低置信”标记;长缺口(大于5 s)不应强行插值,应保留为NaN并标记为“不可用”,避免用虚构数据支撑故障结论。对于离散状态量,缺失期间应保持前一有效状态,并记录状态保持时长。

缺失值处理完成后,应统计每个关键参数的缺失率。若某参数缺失率超过5%,该参数不宜作为故障定位的主证据,只能作为辅助参考。

五、数据质量评估:用指标决定“能不能用”

数据质量评估不是主观判断,而应形成可量化的指标集。建议至少包含:时间对齐残差(互相关峰对应的时延,单位ms)、单位换算校验通过率、野值比例、缺失率、采样率一致性、关键参数物理范围合规率。综合这些指标可给出“可用、降级可用、不可用”三级结论。例如:野值比例小于0.1%、缺失率小于1%、时间对齐残差小于20 ms,可判为“可用”;野值比例0.1%~1%或缺失率1%~5%,判为“降级可用”,需在报告中注明;超过上述范围则判为“不可用”,应先修复数据再分析。

飞行记录数据预处理链路与质量判定 原始记录数据 姿态 100 Hz 大气数据 25 Hz 发动机参数 10 Hz 开关量 5 Hz ① 时间对齐 GPS/PPS 基准轴 互相关求时延 线性插值 / 零阶保持 残差 < 20 ms ② 单位换算 rad→deg ×57.2958 kt→m/s ×0.514444 ft→m ×0.3048 码值×LSB+偏移 ③ 野值剔除 物理范围门限 变化率门限 MAD 统计检测 k=3~5 相关通道一致性判定 ④ 缺失值处理 <0.5 s 线性插值 0.5~5 s 前值保持 >5 s 保留 NaN 缺失率 >5% 降权 ⑤ 数据质量评估指标 • 时间对齐残差(ms) • 单位换算校验通过率 • 野值比例 / 缺失率 • 采样率一致性 • 物理范围合规率 • 关键参数置信标记 质量判定结论 可用 野值<0.1% 缺失<1% 残差<20ms 降级可用 野值0.1~1% 或 缺失1~5% 不可用 超出上述范围,先修复数据 注:图中阈值与参数为教学示例,实际项目应依据记录格式说明与机型规范确定

六、实战示例:一次俯仰角异常的处理过程

假设某次飞行记录中,俯仰角通道在 t=125.4 s 附近出现单点跳变至 85°,而俯仰角速率、升降舵偏度和法向过载均无对应变化。处理过程如下:首先检查单位,确认俯仰角以弧度记录,换算后 85° 对应约 1.483 rad,量级异常;其次用 MAD 检测,该点偏离窗口中位数超过 6 倍 MAD,标记为可疑;再检查相关通道,角速率与舵面均无同向响应,判定为野值;最后用前后各 5 个采样点线性插值替换,并在质量标记中记录“t=125.4 s,俯仰角,野值修正,1 点”。若该点同时伴随角速率和舵面同向跳变,则应保留并转入故障分析,而不是简单剔除。

完成全部预处理后,建议输出一份数据质量摘要表,字段包括:参数名、原始采样率、对齐后采样率、单位换算方式、野值数量与比例、缺失数量与比例、质量等级。该表既是分析报告的一部分,也是后续复现和审计的依据。只有数据质量达到“可用”或明确标注“降级可用”时,飞控故障定位结论才具备可信度。