2. AUV惯导基础

惯性导航是AUV水下自主航行的核心自主导航手段。它不依赖外部声学信标或卫星信号,仅依靠载体内部安装的陀螺仪和加速度计测量角速度与比力,通过航位推算连续输出姿态、速度和位置。理解惯导,必须同时抓住三件事:测的是什么(比力与角速度)、在哪个坐标系里算(惯性系、地球系、导航系、载体系)、怎么一步步算出来(姿态更新—速度更新—位置更新)。本节围绕这三条主线展开,为后续DVL与深度计的组合导航打下基础。

一、惯性导航的基本原理

惯性导航的物理基础是牛顿力学。加速度计测量的并不是载体相对地球的“运动加速度”,而是比力,即除引力之外所有外力(主要是支撑力、推力)产生的等效加速度。陀螺仪测量的是载体相对惯性空间的转动角速度。惯导解算的本质,就是在一个已知初始状态的基础上,用陀螺输出的角速度跟踪姿态变化,用加速度计输出的比力扣除重力影响后积分得到速度,再积分得到位置。

这里有一个关键认识:惯导是递推式系统,当前时刻的输出依赖上一时刻的状态。因此初始对准(确定初始姿态、初始速度和初始位置)的质量直接决定后续精度;同时,积分过程会把传感器零偏、噪声和标度因数误差随时间累积放大。典型MEMS惯导在无外部辅助时,位置误差会随时间快速发散,这也是AUV必须引入DVL测速和深度计测深进行组合的原因。

二、比力方程

比力方程是惯导速度解算的核心公式,它把加速度计测得的比力与导航系下的速度变化联系起来。在导航坐标系(n系)下,常用的比力方程形式为:

v̇ⁿ = fⁿ − (2ωᵢₑⁿ + ωₑₙⁿ) × vⁿ + gⁿ

其中:
  v̇ⁿ     导航系速度变化率
  fⁿ      比力在导航系的投影(fⁿ = Cᵦⁿ · fᵇ)
  ωᵢₑⁿ    地球自转角速度在导航系的投影
  ωₑₙⁿ    导航系相对地球系的转动角速度在导航系的投影
  vⁿ      导航系速度
  gⁿ      重力加速度在导航系的投影

逐项理解:fⁿ是加速度计比力经姿态矩阵Cᵦⁿ从载体系转换到导航系的结果;(2ωᵢₑⁿ + ωₑₙⁿ) × vⁿ是哥氏加速度与向心加速度修正项,来源于地球自转和载体在地球表面运动;gⁿ是重力项。对AUV而言,工作深度通常只有几十到几百米、航速较低(1~3 m/s),地球自转相关项量级很小,在短时中等精度应用中常被简化,但理解完整形式有助于认识误差来源。实际解算中,先由陀螺更新姿态矩阵Cᵦⁿ,再把fᵇ投影为fⁿ,代入上式积分得速度,再积分得位置。

三、四个坐标系的定义

惯导解算全程在多个坐标系之间转换,定义必须统一、符号必须一致,否则会出现方向性错误。

  • 惯性坐标系(i系):原点在地球中心,坐标轴相对恒星固定不随地球转动,是牛顿定律成立的参考基准,陀螺测量的角速度即相对i系。
  • 地球坐标系(e系):原点在地球中心,随地球自转,z轴沿地球自转轴指向北极,x轴指向本初子午线。用于描述载体相对地球的位置(经度、纬度、高度)。
  • 导航坐标系(n系):惯导解算的工作坐标系,AUV常用东北天(ENU)当地水平坐标系:x轴指东、y轴指北、z轴指天。速度、位置、姿态输出通常都在n系下表达。
  • 载体坐标系(b系):固连于AUV本体,通常取右—前—上(RFU)或前—右—下(FRD),原点在惯导安装基准点。加速度计和陀螺的原始输出都在b系下。
坐标系约定一旦确定,全文与图示必须保持一致。例如采用ENU导航系时,天向速度向上为正,深度计给出的深度向下为正,组合时需做符号转换,这是实际工程中最容易出错的地方之一。
惯导四坐标系关系与解算主链 坐标系定义与转换关系 惯性系 i 相对恒星固定 地球系 e 随地球自转 导航系 n(ENU) 东—北—天 载体系 b 右—前—上 自转 位置 姿态 Cᵦⁿ 陀螺测相对 i 系角速度 → 姿态矩阵 Cᵦⁿ;加速度计测 b 系比力 fᵇ → 投影到 n 系得 fⁿ 位置输出用 e 系经纬高,速度与姿态输出用 n 系,原始测量在 b 系 惯导解算主链(递推) 初始对准 姿态/速度/位置 初值 姿态更新 陀螺 ωᵇ → Cᵦⁿ 四元数/方向余弦 速度更新 比力方程积分 fⁿ − 哥氏 + gⁿ 位置更新 vⁿ 积分 经纬高/深度 递推反馈:上一时刻状态作为下一时刻初值,误差随时间累积 外部辅助:DVL 测速、深度计测深 → 抑制惯导积分发散(后续章节展开)

四、惯导解算流程概述

惯导解算是一个按固定周期(如100 Hz或200 Hz)循环执行的递推过程,主要分四步:

  1. 初始对准:AUV下水前或水面阶段,利用陀螺感知地球自转角速度和重力方向,确定初始姿态矩阵Cᵦⁿ,同时给定初始速度(通常为零或已知)和初始位置(GPS或装订值)。对准精度是后续误差的起点。
  2. 姿态更新:用陀螺输出的角速度ωᵇ,结合四元数微分方程或方向余弦矩阵微分方程,更新Cᵦⁿ。工程上常用四元数法避免奇异并减少计算量,再归一化后转换为方向余弦矩阵。
  3. 速度更新:将加速度计比力fᵇ经Cᵦⁿ投影为fⁿ,代入比力方程,扣除哥氏项和重力项后积分,得到导航系速度vⁿ。
  4. 位置更新:对vⁿ积分得到位置变化,结合地球模型换算为经度、纬度、高度(或深度)。

以一个教学示例说明量级:假设某MEMS陀螺零偏为10 °/h(约0.0485 °/s),在无辅助情况下,姿态误差约以10 °/h累积,1小时后航向误差可达约10 °;若AUV以1.5 m/s航行,航向误差引起的横向位置误差在1小时后可达数百米量级。这组数字属于合理假设的示例,用于说明为何纯惯导难以满足长航时AUV需求,而非实测结果。

实际实现中需注意三点:一是采样周期与积分算法要匹配,高频抖动需滤波;二是姿态矩阵必须正交归一化,否则长期递推会漂移;三是n系与b系的符号约定、深度正方向必须与DVL、深度计统一,组合导航前先做坐标一致性检查。

把上述内容串起来看:陀螺和加速度计在b系给出原始测量,姿态矩阵Cᵦⁿ完成b系到n系的桥梁作用,比力方程在n系完成速度解算,位置再由速度积分得到,而e系与i系则提供了地球自转与位置描述的物理背景。掌握这条主链,才能在后续引入DVL速度观测和深度计深度观测时,清楚每一个观测量该修正解算流程中的哪一环。