3. 偏差数据采集与实验设计

数字孪生模型校准的质量,首先取决于偏差数据本身是否“可辨识”。如果传感器布置不合理、采样频率过低、激励信号无法激发关键动态,那么后续再精细的优化算法也只能在低信息量数据上做无效拟合。本节围绕“面向校准的数据采集与实验设计”展开,重点说明传感器布置、采样频率选择、激励信号设计、数据预处理与异常值剔除五个环节的工程做法。

一、传感器布置:让偏差在空间上可观测

传感器布置的核心目标是使模型关键状态量与输出量的偏差在测量空间中可观测、可分离。工程上通常遵循三条原则:第一,覆盖模型敏感区域,即对模型输出影响大的位置优先布点;第二,避免冗余,两个传感器若长期高度相关,则其中一个对参数辨识贡献有限;第三,兼顾物理可达性与安装约束。

以某三自由度机械臂数字孪生校准为例,模型参数包括关节摩擦系数、连杆等效惯量和关节刚度。若只在末端布置一个位置传感器,则摩擦与刚度偏差在末端轨迹上高度耦合,难以分离。合理做法是在每个关节处布置编码器测关节角,在末端布置视觉标记测末端位姿,并在关节驱动电流回路中采集电流信号以间接反映摩擦力矩。这样,关节角序列对惯量与刚度敏感,电流序列对摩擦敏感,末端位姿对综合误差敏感,三类信号形成互补。

定量上可用敏感度矩阵辅助布点。设模型输出 y 对参数 θ 的敏感度为 S = ∂y/∂θ,若某传感器位置的 S 列向量接近线性相关,则说明该位置对参数分离贡献低。实际布点时可计算候选位置的敏感度矩阵条件数,条件数越小,参数可辨识性越好。

二、采样频率选择:兼顾动态信息与噪声抑制

采样频率不是越高越好。过低会丢失动态信息并产生混叠,过高则引入更多高频噪声、增加存储与计算负担,还可能因差分运算放大噪声。工程上通常依据信号最高关注频率 fmax 选择采样频率 fs,满足 fs ≥ (5~10)fmax。对于需要微分求速度、加速度的场合,可提高到 10~20 倍。

以关节伺服系统为例,若关注带宽约 50 Hz,则采样频率取 500 Hz~1 kHz 较为合适。若同时采集电流与位置,建议使用同一硬件时钟触发,避免通道间时间戳漂移。对于温度、环境湿度等慢变量,采样频率可降至 1~10 Hz,但需与快变量数据在时间轴上对齐,通常采用线性插值或零阶保持后重采样。

三、激励信号设计:激发参数敏感的动态

激励信号应满足持续激励条件,即输入信号在关注频段内具有足够能量,使待辨识参数对应的动态被充分激发。常用激励包括阶跃、正弦扫频、伪随机二进制序列(PRBS)和多正弦组合信号。

阶跃信号适合辨识静态增益和低频动态,但频谱能量集中在低频,对高频动态激发不足。正弦扫频可在设定频段内均匀激励,适合辨识频响特性。PRBS 信号实现简单、频带较宽,适合在线辨识。多正弦信号可针对特定频点加强激励,适合已知关注频段的场景。

以机械臂关节摩擦辨识为例,若仅做缓慢点到点运动,关节速度长期处于低速段,摩擦模型中的 Stribeck 效应无法被激发。可设计一组包含低速匀速段、中速往复段和换向段的激励轨迹,使关节速度覆盖 0.01~1 rad/s,并在换向点附近加密采样,从而让库仑摩擦、粘滞摩擦和 Stribeck 参数均获得足够激励。

面向校准的偏差数据采集链路 传感器布置 关节编码器 → 关节角 电流传感器 → 力矩 末端视觉 → 位姿 温度传感器 → 热漂移 敏感度矩阵条件数最小化 避免冗余与耦合 采样频率选择 快变量:500 Hz~1 kHz 慢变量:1~10 Hz fs ≥ (5~10) fmax 统一时钟触发 防混叠、防时间漂移 重采样对齐时间轴 激励信号设计 阶跃:低频增益 扫频:频响特性 PRBS:宽频在线辨识 多正弦:指定频点加强 持续激励条件 覆盖关注频段与工况 数据预处理 去均值、去趋势 低通滤波 时间对齐与重采样 单位统一与坐标变换 保留动态、抑制噪声 避免相位失真 异常值剔除方法 3σ 准则:适用于近似正态分布残差 MAD 法:对离群点更稳健 滑动窗口残差检测:适合在线剔除 物理约束校验:超限即标记 剔除后需记录并回溯原因,不可静默丢弃 校准数据集 输入-输出对 + 质量标记 注:各环节参数为工程示例,实际取值需结合对象带宽、噪声水平与校准目标确定

四、数据预处理:把原始信号变成可校准数据

原始采集数据通常包含直流偏置、趋势项、高频噪声和通道间时间偏差,直接用于校准会引入系统性误差。预处理一般包括以下步骤。

去均值与去趋势。传感器零漂和温漂会表现为缓慢趋势项。对位置信号可先减去初始值,再用多项式或高通滤波去除低频趋势。对电流信号,需注意去均值不能去掉真实的平均力矩分量,应结合工况判断。

低通滤波。为抑制高频噪声,可设计截止频率略高于关注带宽的低通滤波器。例如关注 50 Hz 以内动态,可选截止频率 80~100 Hz 的二阶巴特沃斯滤波器。滤波会引入相位滞后,若后续需要计算速度、加速度,建议采用零相位滤波(如 filtfilt)避免相位失真。

时间对齐与重采样。多传感器数据往往来自不同采样率。可统一重采样到最高采样率,或统一到校准算法要求的采样率。重采样时优先使用抗混叠插值,避免直接抽取造成混叠。

单位与坐标统一。编码器脉冲需换算为弧度,电流需结合力矩常数换算为力矩,视觉位姿需统一到同一坐标系。坐标变换顺序和方向约定必须在校准前明确,并在全文及图示中保持一致。

五、异常值剔除:保护校准不被坏点带偏

异常值来源包括传感器瞬时故障、通信丢包、电磁干扰和机械碰撞。常用剔除方法有:

  • 3σ 准则:计算残差均值与标准差,超出 ±3σ 的点标记为异常。适用于近似正态分布且样本量较大的场景。
  • MAD 法:用中位数绝对偏差替代标准差,对离群点更稳健。当数据中存在少量极端值时,MAD 比 3σ 更可靠。
  • 滑动窗口残差检测:在窗口内建立局部模型,预测值与实测值偏差超过阈值即标记。适合在线采集时实时剔除。
  • 物理约束校验:关节角超出机械限位、电流超过额定值、速度超过设计上限等,可直接判定为无效数据。
异常值剔除后必须记录剔除位置、数量与原因,不能静默丢弃。若某通道异常值比例超过 5%,应优先排查传感器安装、接地或通信问题,而不是继续加大剔除力度。

六、实验设计闭环示例

以某关节摩擦参数校准为例,可形成如下实验设计闭环:

环节示例设置判断依据
传感器布置关节编码器 + 电流传感器 + 末端视觉敏感度矩阵条件数 < 20
采样频率位置/电流 1 kHz,温度 1 Hzfs ≥ 10×关注带宽
激励信号低速匀速 + 中速往复 + 换向段速度覆盖 0.01~1 rad/s
预处理去趋势 + 80 Hz 零相位低通 + 重采样相位滞后 < 2°
异常值剔除MAD + 物理限位校验剔除比例 < 3%

完成上述步骤后,得到的是带有质量标记的校准数据集。该数据集应包含输入激励、输出响应、时间戳、通道标识和有效标记,供后续参数辨识与偏差修正使用。若数据质量不达标,应回到传感器布置或激励设计环节调整,而不是直接进入优化算法。

实践建议:在正式采集前先做一次短时预采集,检查各通道信噪比、时间同步和量程设置。预采集可显著减少正式实验中的返工。