3. 水下INS误差传播机理

纯惯性导航的误差不是突然出现的，而是由初始对准残差、惯性器件零偏与标度因数误差、以及导航解算中的耦合关系逐步“放大”和“调制”出来的。要理解水下INS为何在长时间潜航后位置误差会发散，必须回到姿态、速度、位置误差的微分方程，看清误差在三个回路之间如何传递，以及舒勒振荡、傅科效应如何改变误差的形态与方向。

一、误差状态与坐标系约定

本节统一采用当地水平指北地理坐标系（n系：东-北-天，E-N-U）作为导航解算参考，载体坐标系（b系：右-前-上）与n系之间的姿态用欧拉角表示。定义姿态误差角 φ = [φ_E, φ_N, φ_U] 为计算n系相对真实n系的小角度旋转，速度误差 δv = [δv_E, δv_N, δv_U]，位置误差用纬度误差 δL、经度误差 δλ 和高度误差 δh 描述。所有误差均定义为“计算值减真实值”，全文与图示保持这一符号约定。

二、姿态、速度、位置误差微分方程

在静基座或低速水下航行条件下，忽略高度通道，姿态误差方程可写为：

φ̇_E = -δv_N / R + φ_N·ω_ie·sinL - φ_U·ω_ie·cosL - ε_E
φ̇_N =  δv_E / R - φ_E·ω_ie·sinL - ε_N
φ̇_U =  δv_E·tanL / R + φ_E·ω_ie·cosL - ε_U

其中 R 为地球曲率半径，ω_ie 为地球自转角速度，L 为纬度，ε 为等效陀螺漂移在n系的分量。可以看到，姿态误差一方面被速度误差驱动（δv/R 项），另一方面被地球自转与纬度调制（ω_ie·sinL、ω_ie·cosL 项），后者正是傅科效应的数学来源。

速度误差方程的核心是比力投影误差与哥氏项：

δv̇_E = -f_U·φ_N + f_N·φ_U + 2ω_ie·sinL·δv_N + δv_E·v_N·tanL/R + ∇_E
δv̇_N =  f_U·φ_E - f_E·φ_U - 2ω_ie·sinL·δv_E - δv_E·v_E·tanL/R + ∇_N

式中 f 为比力，∇ 为加速度计零偏在n系的分量。关键点在于：姿态误差 φ 通过比力 f 耦合进速度误差，而速度误差又通过 δv/R 反馈到姿态误差，形成闭合回路。位置误差方程则是速度误差的积分：

δL̇ = δv_N / R
δλ̇ = δv_E / (R·cosL)

三者构成一个典型的“姿态—速度—位置”串联加反馈结构，任何一环的常值偏差都会在回路中转化为周期性或发散性误差。

三、舒勒振荡与傅科效应的作用

舒勒振荡是惯性导航中最基本的误差抑制机制。当平台或数学平台满足舒勒调谐条件时，姿态误差与速度误差形成周期约84.4分钟的等幅振荡，而不是单调发散。这意味着单纯由初始姿态误差引起的速度误差不会无限增长，而是被“锁”在一个有界振荡中。但要注意：舒勒回路只约束姿态和速度，位置误差是速度误差的积分，因此位置误差仍会随时间线性累积——这正是纯惯性导航位置发散的根源。

傅科效应则来自地球自转的调制作用。在误差方程中，ω_ie·sinL 和 ω_ie·cosL 项使舒勒振荡的平面绕当地垂线缓慢旋转，旋转角速率约为 ω_ie·sinL。在纬度30°处，该速率约为每小时7.5°，一天旋转约180°。其工程含义是：误差振荡的方向不是固定的，而是随时间缓慢改变，导致东向和北向误差分量交替主导，位置误差轨迹呈现螺旋或椭圆形态，而不是简单的直线漂移。

四、仿真示例：纯惯性误差发散过程

下面给出一个教学用仿真设定（参数为合理假设，非实测数据）：陀螺常值漂移 ε = 0.01°/h，加速度计零偏 ∇ = 50 μg，初始姿态误差 φ = [0.02°, 0.02°, 0.1°]，初始速度误差为零，航行纬度30°，航向正北，航速5 kn。用上述误差方程做数值积分，可观察到以下典型现象：

• 0～1小时：姿态误差以84.4分钟周期振荡，幅值约0.02°～0.05°；速度误差振荡幅值约0.1～0.3 m/s；位置误差尚小，约数十米量级。
• 1～6小时：舒勒振荡持续，但加速度计零偏经二次积分后使位置误差呈抛物线增长，6小时末位置误差可达约1～2 km量级（取决于零偏方向与航向关系）。
• 6～24小时：傅科效应使误差振荡平面旋转，东向与北向误差交替增大，位置误差轨迹呈缓慢展开的螺旋；24小时末位置误差可达数十公里量级，纯惯性导航已无法满足水下隐蔽航行的定位需求。

五、工程上的观察与判断

在实际水下INS调试中，可通过以下方式验证误差传播机理是否符合预期：静基座下记录姿态与速度输出，做频谱分析应能看到约84.4分钟的舒勒周期峰值；改变初始纬度设定，观察振荡平面的旋转速率是否随 sinL 变化；人为注入已知陀螺漂移，检查位置误差增长斜率是否与理论二次积分趋势一致。若实测位置误差增长远快于理论估算，应优先排查初始对准残差、标度因数误差或温度引起的零偏变化，而不是简单归因于“惯导漂移大”。

理解误差传播机理的最终目的，是为后续的阻尼、重调和组合导航设计提供依据：舒勒回路决定了阻尼网络的设计频段，傅科效应决定了误差方向随时间的变化规律，而位置误差的发散速率则直接决定了外部参考的更新周期与精度指标。