2. 参数辨识中的误差来源剖析

参数辨识得到的估计值几乎不可能与真实参数完全一致，理解误差从何而来、如何沿着辨识链路传播，是判断置信区间是否可信的前提。工程上通常把误差归为四类：测量噪声、模型结构误差、数值计算误差与数据不足。它们并非孤立存在，而是从数据采集开始，经过模型构建、算法求解，最终叠加到参数估计与置信区间上。本节先逐一剖析这四类误差的物理与数学本质，再通过一个直流电机参数辨识案例，追踪误差的传播路径。

一、测量噪声：随机误差的统计根源

测量噪声来自传感器精度、电磁干扰、量化误差与采样抖动，通常可建模为零均值随机序列。设真实输出为 y(k) = φᵀ(k)θ₀ + v(k)，其中 v(k) 为噪声，θ₀ 为真实参数。最小二乘估计 θ̂ = (ΦᵀΦ)⁻¹Φᵀy 的协方差为 σ²(ΦᵀΦ)⁻¹，可见噪声方差 σ² 直接放大参数估计的散布。噪声的影响有两个特点：一是可通过增加样本量或平均滤波降低，二是若噪声与回归量相关（如闭环辨识中的有色噪声），最小二乘会产生有偏估计，此时需改用工具变量法或增广最小二乘。

二、模型结构误差：系统性偏差的主因

模型结构误差源于所建模型与真实系统动态之间的失配，例如忽略执行器饱和、摩擦非线性、未建模的高频动态或时滞。这类误差不会因数据增多而消失，反而在数据越丰富时越明显。以电机为例，若只辨识一阶惯性模型而实际存在电气时间常数，则辨识出的机电时间常数会系统性偏大或偏小。结构误差的典型特征是残差呈现规律性（如与输入幅值相关），可通过残差自相关检验与交叉验证识别。工程上常采用“由简到繁”的模型阶次递增策略，结合AIC/BIC准则在拟合精度与复杂度间折中。

三、数值计算误差：算法与病态问题

数值误差包括浮点舍入、矩阵病态导致的求逆不稳定、迭代未收敛以及激励信号不充分引起的回归矩阵秩亏。当输入信号为阶跃或常值时，ΦᵀΦ 接近奇异，参数估计对扰动极其敏感，置信区间会被严重低估。实践中应使用伪逆或QR分解代替直接求逆，并对回归矩阵做条件数检查；若条件数超过 1e6，说明激励不足，需改用伪随机二进制序列（PRBS）等持续激励信号。此外，采样率过高会引入数值微分噪声，过低则丢失动态信息，需按系统带宽的5～10倍选取采样频率。

四、数据不足：信息量约束

数据不足既指样本数量少，也指激励频段覆盖不全。参数可辨识性要求输入信号在待辨识动态的频段内具有足够能量。例如辨识二阶系统需要输入包含接近自然频率的分量，否则阻尼比估计的方差极大。数据不足的直接后果是置信区间异常宽，甚至出现参数不可辨识（多个参数组合产生相同输出）。判断方法是观察Fisher信息矩阵的特征值分布，小特征值对应的参数方向即为弱可辨识方向。

五、工程案例：直流电机机电时间常数辨识

设某直流电机在电压阶跃激励下，转速响应近似一阶惯性环节 ω(s)/U(s) = K/(Ts+1)，待辨识参数为增益 K 与时间常数 T。采样频率取1 kHz，采样时长2 s，转速传感器分辨率0.5 r/min，叠加标准差约2 r/min的随机噪声。采用最小二乘拟合离散模型 ω(k) = a·ω(k-1) + b·u(k-1)，再由 T = -Ts/ln(a) 反算时间常数。

若仅用一段阶跃数据，回归矩阵条件数约为 3×10⁴，尚可接受；但若电机存在明显的静摩擦死区，低速段转速被“卡住”，模型结构误差使 a 的估计偏大，反算出的 T 偏大约15%。此时残差在低速段呈现系统性正偏，残差自相关函数在滞后1～3阶显著非零，提示结构失配。改用含摩擦补偿的二阶模型后，残差白性检验通过，T 的估计回落至合理范围，95%置信区间由原来的 ±0.08 s 收窄至 ±0.03 s。

进一步地，若把采样时长缩短到0.3 s，转速尚未进入稳态，ΦᵀΦ 的最小特征值骤降，T 的置信区间扩张到 ±0.25 s，几乎失去工程意义。这说明数据不足对置信区间的影响往往比噪声更致命——噪声只增加方差，而数据不足可能直接破坏可辨识性。

六、误差传播的定量视角

把四类误差统一到估计误差分解式 θ̂ - θ₀ = (ΦᵀΦ)⁻¹Φᵀv + (ΦᵀΦ)⁻¹ΦᵀΔη + δ_num，其中 Δη 为结构误差引起的等效输出偏差，δ_num 为数值误差。第一项随噪声方差线性增长、随信息矩阵增大而减小；第二项由结构失配决定，与样本量无关；第三项在病态条件下可能被急剧放大。置信区间若只按第一项计算，就会在存在结构误差时给出“过窄且偏离真值”的区间，这是工程中最危险的误判。因此，报告参数辨识结果时，应同时给出残差诊断、条件数与激励设计说明，使置信区间具备可解释性。