4. 数据采集与预处理
模型参数辨识的精度上限,往往在数据采集阶段就已经被决定。辨识算法只能从数据中提取信息,如果传感器量程不合适、采样率不足、标定存在系统偏差,或者原始信号中混入了未处理的异常值,后续再精细的优化也只能拟合噪声。因此,数据采集与预处理不是“准备工作”,而是辨识链路中与模型结构同等重要的一环。本节围绕传感器选型与标定、采集系统配置、噪声与异常值处理、滤波与重采样、去趋势与数据分段、数据质量评估六个方面展开,重点说明每一步的实际操作与判断依据。
传感器选型与标定
传感器选型的核心是让被测物理量的实际变化范围落在传感器有效量程的中段,而不是接近满量程或接近零。以某电机转速辨识试验为例,若目标转速范围为 0~3000 r/min,选用量程 0~5000 r/min 的编码器,则常用工作区间位于量程的 60% 附近,分辨率与线性度都能得到较好利用;若选用 0~100000 r/min 的传感器,虽然量程覆盖,但低转速段的分辨率会严重不足。选型时还需关注灵敏度、分辨率、频响带宽、非线性度、重复性和温度漂移等指标。频响带宽应至少覆盖被测信号最高关注频率的 5~10 倍,否则高频成分会被传感器自身衰减,导致辨识出的动态参数偏小。
标定是建立传感器输出与被测物理量之间确定关系的过程。静态标定通常采用多点法:在量程内均匀选取 5~10 个标准输入点,记录传感器输出,用最小二乘拟合直线或多项式。例如某压力传感器标定数据中,标准压力为 0、2、4、6、8、10 MPa,对应输出为 0.002、1.998、4.003、5.997、8.001、10.002 V,拟合得到 y = 0.9998x + 0.0012,斜率接近 1 说明灵敏度准确,截距 0.0012 V 对应约 0.0012 MPa 的零点偏置,应在采集软件中作为偏置补偿扣除。动态标定则用于确定传感器幅频和相频特性,对高频辨识尤为重要。标定结果应记录标定日期、标准器编号、环境温度和拟合残差,作为数据可信度的依据。
数据采集系统配置
采集系统配置的关键参数包括采样率、分辨率、量程、耦合方式和触发方式。采样率的选择依据奈奎斯特采样定理,同时考虑抗混叠滤波器的截止频率。若关注信号最高频率为 f_max,采样率 f_s 应满足 f_s ≥ 2.5~5 × f_max,工程上常取 5~10 倍。例如辨识某机械系统 0~50 Hz 的动态特性,采样率取 500 Hz 或 1000 Hz 较为稳妥。分辨率方面,16 位 ADC 在 ±10 V 量程下的理论分辨率为 20 V / 65536 ≈ 0.305 mV,若信号有效变化仅 50 mV,则量化台阶约占信号幅值的 0.6%,通常可接受;若信号更小,应缩小量程或改用更高位数 ADC。
多通道采集时,必须保证各通道同步采样,否则通道间相位差会直接污染辨识结果。对于振动、电流、电压等多物理量同步采集,应使用同一时钟源触发所有通道,避免使用软件轮询方式依次读取。触发方式上,自由运行适合连续监测,边沿触发适合捕捉瞬态过程,预触发功能可以保留触发点之前的数据,对辨识冲击响应特别有用。采集前应进行通道自检:短接输入端观察零点漂移,接入标准信号观察增益是否正确,确认无通道接反或接地环路。
噪声与异常值处理
噪声分为随机噪声和系统性干扰。随机噪声表现为高频小幅波动,可通过统计方法抑制;系统性干扰如 50 Hz 工频串扰、接地环路、电磁耦合,则需要在硬件层面解决,例如使用屏蔽线、单点接地、隔离放大器。异常值处理首先要识别,常用方法包括 3σ 准则、四分位距(IQR)准则和基于滑动窗口的局部离群检测。以 3σ 准则为例,计算序列均值 μ 和标准差 σ,将 |x_i - μ| > 3σ 的点标记为异常。但 3σ 对非高斯分布和趋势数据不稳健,此时 IQR 更合适:Q1 和 Q3 为第一、第三四分位数,IQR = Q3 - Q1,超出 [Q1 - 1.5IQR, Q3 + 1.5IQR] 的点视为异常。
异常值的处理方式有剔除、替换和保留标记三种。对于孤立尖峰,可用相邻点的线性插值替换;对于连续异常段,应剔除该段并记录缺失区间,而不是强行插值。需要特别注意的是,异常值可能携带真实信息,例如传感器饱和、机械碰撞或电路故障,直接删除可能掩盖系统问题。因此,异常值处理应保留原始数据副本,并在处理日志中记录异常点位置、判定依据和处理方式。
滤波与重采样
滤波的目的是在保留有效信号的前提下抑制噪声。低通滤波是最常用的手段,截止频率应略高于关注信号最高频率,避免过度滤波导致相位滞后和幅值衰减。以二阶巴特沃斯低通滤波器为例,若关注频率为 50 Hz,截止频率取 80~100 Hz,可在 50 Hz 处保持较平坦的幅频响应,同时对 200 Hz 以上噪声提供约 20~30 dB 的衰减。对于存在明显工频干扰的信号,可叠加 50 Hz 陷波滤波器,但陷波会引入相位失真,使用前应评估对辨识结果的影响。
重采样用于统一不同采样率的数据或降低数据量。降采样前必须先进行抗混叠低通滤波,否则高频噪声会折叠到低频段,污染辨识结果。例如将 1000 Hz 数据降采样到 200 Hz,应先将信号低通滤波到 100 Hz 以下,再每 5 个点取 1 个点。升采样则通过插值实现,常用线性插值、样条插值和 sinc 插值,其中 sinc 插值理论上最精确但计算量大。重采样后应检查时间轴是否均匀,避免因插值引入额外的时间偏差。
去趋势与数据分段
去趋势用于消除信号中的缓慢漂移或直流偏置,使数据围绕零均值波动。常用方法包括多项式拟合去趋势、滑动平均去趋势和差分去趋势。例如某温度信号在 10 min 内从 25.0 ℃ 缓慢漂移到 25.8 ℃,若直接用于辨识,漂移会被误认为低频动态。用一次多项式拟合该漂移并减去,可保留真实的温度波动成分。去趋势的阶次不宜过高,否则会削弱真实的低频动态信息。
数据分段是将长序列划分为若干适合辨识的短段。分段依据可以是时间等长、工况变化点或激励事件。例如辨识系统阶跃响应时,以阶跃发生时刻为起点,截取 5~10 倍时间常数的数据段;辨识频响函数时,按激励信号的周期分段。分段时应保证每段包含完整的动态过程,避免截断在瞬态中间。段与段之间可设置重叠,以提高数据利用率,但重叠率过高会导致段间相关性增强,影响统计独立性。
数据质量评估方法
数据质量评估应贯穿采集与预处理全过程,常用指标包括信噪比、有效数据率、缺失率、异常值比例、频谱一致性、重复性误差和残差白度。信噪比可通过信号功率与噪声功率之比估算,工程上要求辨识用数据信噪比一般不低于 20 dB。有效数据率指剔除异常和缺失后可用于辨识的数据占比,低于 70% 时应重新采集。频谱一致性用于检查同一工况多次采集的频谱是否一致,若差异显著,说明系统或采集条件不稳定。
残差白度检验是辨识后验证数据质量的重要手段:若模型残差仍存在明显自相关或周期性成分,说明数据中仍有未被模型解释的结构,可能是噪声处理不充分或模型结构不完整。实际工作中,建议在预处理前后分别计算上述指标,形成数据质量报告,作为辨识结果可信度的支撑材料。