4. 惯性导航误差源(一)

惯性导航系统依靠陀螺仪感知载体相对惯性空间的角运动,再通过积分把角速度变成姿态,把姿态变成加速度的投影方向,最终得到速度和位置。这条链路中,陀螺仪的任何微小误差都会随着积分不断累积,因此理解陀螺仪的误差源,是分析水下导航精度衰减的起点。本节聚焦三类最基础、也最常被建模的陀螺误差:零偏、标度因数误差与随机噪声,从物理机理出发给出数学描述,并说明它们在工程上如何被辨识和使用。

一、零偏:角速度为零时的“虚假输出”

零偏(Bias)指陀螺仪在输入角速度为零时输出的非零角速度,单位通常为 °/h 或 rad/s。它的物理来源包括:MEMS 陀螺中质量块支撑结构的残余应力与封装热应力、光纤陀螺中光纤环的 Shupe 效应与光源波长漂移、激光陀螺的锁区效应,以及电路零点漂移和 A/D 量化偏置。零偏最直接的影响是姿态误差随时间线性增长:若零偏为 b,则经过时间 t 后姿态误差约为 b·t。例如 0.1 °/h 的常值零偏,1 小时后产生约 0.1° 的姿态误差,10 小时后约 1°,而姿态误差又会通过重力投影把误差耦合进水平加速度通道,形成随时间的振荡发散。

工程上常把零偏拆成三部分描述:常值零偏 b₀(上电重复性)、零偏稳定性( Allan 方差在 1 小时或 10 小时内的波动)以及零偏温度系数。数学上可写为:

ω_meas(t) = ω_true(t) + b₀ + b_T(T) + b_r(t) + K·ω_true(t) + n(t)
其中:
  b₀        常值零偏(° /h)
  b_T(T)    温度相关零偏,常用 b_T = k_T·(T - T₀) 近似
  b_r(t)    零偏随机游走(慢变)
  K         标度因数误差(无量纲)
  n(t)      角度随机噪声

需要注意,零偏并非完全不可观测。在静基座或已知航向的初始对准阶段,可以通过长时间平均陀螺输出估计水平零偏;水下载体若配有 DVL 或 USBL,也可在组合导航滤波中把零偏作为状态量在线估计。但零偏的可观测性依赖于载体机动,直线匀速航行时航向零偏往往难以与航向误差区分,这是水下组合导航中需要特别警惕的问题。

二、标度因数误差:输入与输出之间的比例失配

标度因数(Scale Factor)是陀螺输出与真实角速度之间的比例系数。理想情况下输出等于输入,但实际器件因敏感结构几何尺寸偏差、检测电路增益误差、温度引起的材料参数变化,会使比例系数偏离标称值。标度因数误差通常用相对误差 K 表示,即输出 = (1+K)·真实角速度。它与零偏的关键区别在于:零偏是加性误差,标度因数误差是乘性误差,只有存在角运动时才显现,且误差大小与输入角速度成正比。

对水下载体而言,这带来一个容易被忽视的后果:转弯、俯仰、横滚等机动越剧烈,标度因数误差造成的姿态误差越大;而长时间直航时它几乎不显现,导致标定不充分的系统在机动后出现“突然变差”的现象。此外,三轴陀螺之间的非正交误差(交叉耦合)与标度因数误差常被合并为标定矩阵处理:

[ωx]   [1+Kx   mxy   mxz] [ωx_true]
[ωy] = [myx  1+Ky   myz]·[ωy_true]
[ωz]   [mzx   mzy  1+Kz] [ωz_true]

其中对角项为各轴标度因数误差,非对角项为轴间交叉耦合。工程标定的典型做法是使用转台在多个角速率点(如 ±1、±5、±30、±100 °/s)正反旋转,记录输出并做最小二乘拟合,得到每个轴的增益与零偏;再通过三轴依次转动辨识交叉耦合项。标定后的残差通常可降到 10⁻⁴ 量级,但温度变化会使标定矩阵漂移,因此高精度系统会建立温度分段标定表。

三、随机噪声:无法用常值修正的部分

即使零偏和标度因数都被标定,陀螺输出仍会存在随机波动。这类误差不能用确定性的常值补偿,只能用统计特性描述。最常用的两个指标是角度随机游走(ARW,单位 °/√h)和零偏不稳定性(Bias Instability,单位 °/h)。ARW 来源于白噪声角速率积分,其对应的角度误差标准差随时间按 √t 增长;零偏不稳定性来源于 1/f 噪声和慢变过程,表现为 Allan 方差曲线上的最低点。

Allan 方差是辨识这些噪声类型的标准工具:把陀螺静止采样的数据按不同簇时长 τ 分组,计算相邻簇平均值的差,得到方差 σ²(τ)。在双对数坐标上,斜率 −1/2 段对应角度随机游走,斜率 0 的平台对应零偏不稳定性,斜率 +1/2 段对应角速率随机游走。工程上常据此为卡尔曼滤波器设定过程噪声协方差 Q:ARW 决定姿态积分中的白噪声强度,零偏不稳定性决定零偏状态随机游走的驱动噪声。若 Q 设置过小,滤波器会过度信任陀螺,导致姿态误差累积;设置过大则过度依赖外部观测,在 DVL 失锁时反而更差。

陀螺仪误差源 → 数学描述 → 姿态误差表现 零偏 b 加性误差 零输入仍有输出 单位:°/h 标度因数误差 K 乘性误差 误差 ∝ 输入角速度 含轴间交叉耦合 随机噪声 统计特性描述 ARW / 零偏不稳定性 单位:°/√h、°/h ω_meas = ω_true + b 姿态误差 ≈ b·t 随时间线性增长 可通过对准/组合估计 ω_meas = (1+K)·ω_true 机动时误差放大 直航时几乎不显现 转台多点标定 σ(θ) ∝ ARW·√t Allan 方差辨识 斜率 −1/2:ARW 平台段:零偏不稳定性 姿态缓慢漂移 经重力投影耦合 进入水平速度通道 形成舒勒振荡 机动后姿态跳变 转弯/俯仰后误差累积 需温度分段标定 组合导航中在线修正 滤波参数设定 ARW → 姿态过程噪声 零偏不稳定性 → 零偏 随机游走驱动噪声 零偏为加性、标度因数为乘性、随机噪声为统计性;三者叠加构成陀螺总误差,需分别标定与建模

四、从误差模型到工程使用

把三类误差放在一起看,可以得到一个实用的认识:零偏决定“静止时漂多快”,标度因数误差决定“动起来后偏多少”,随机噪声决定“滤波时该信多少”。在实际水下导航任务中,这三者并不是孤立处理的。以一台典型光纤陀螺为例,若其零偏稳定性为 0.05 °/h、ARW 为 0.02 °/√h、标度因数误差为 50 ppm,则 2 小时直航后零偏贡献约 0.1° 姿态误差,ARW 贡献约 0.028°,而一次 30 °/s 的转向若标度因数未修正,可带来约 0.0015° 的瞬时误差——单次不大,但频繁机动会累积。因此工程上通常先做实验室转台标定消除常值零偏与标度因数,再用 Allan 方差分析确定随机噪声参数,最后在组合导航滤波中把残余零偏作为状态量在线估计。

辨识零偏时,务必区分“上电零偏”和“运行中零偏漂移”。前者可通过初始对准平均消除,后者只能靠组合导航观测或定期旋转调制抑制。水下载体长时间无法上浮时,后者往往是航向误差的主要来源。

还需要强调参照系约定的一致性:上述零偏与标度因数均在陀螺敏感轴坐标系下定义,姿态误差则在导航坐标系(如当地水平指北系)下观察。做误差传播分析时,必须通过方向余弦矩阵把敏感轴误差投影到导航系,否则正负号和量级都会出错。这一投影关系也是后续分析加速度计误差和组合导航观测方程的基础。