1. USBL安装偏角标定任务与坐标系梳理

标定任务边界与坐标系梳理

本例的对象是一台刚性固定在船体上的USBL换能器(下称“换能器”)。标定的目标不是测声速、也不是测换能器内部阵元相位,而是求出一组安装偏角:换能器坐标系相对船体坐标系的三个欧拉角。只要这组角确定,后续就能把USBL测得的换能器坐标系下的目标方位,正确旋转到船体坐标系,再叠加姿态与航向进入大地坐标系。

先明确任务边界,避免把标定问题扩大化:

  • 待估量:3个安装偏角(航向偏角 Δψ、俯仰偏角 Δθ、横滚偏角 Δφ),外加换能器在船体坐标系下的安装位置矢量(若仅做方位标定,位置可暂按设计值固定)。
  • 已知量:船体坐标系定义、换能器机械安装图、GNSS/惯导提供的姿态与航向、声学观测(斜距、方位角、俯仰角)。
  • 不估量:声速剖面误差、换能器阵元相位误差、时延偏差、船体挠曲——这些属于其他标定或补偿环节,本例只把它们列为误差源,不纳入待估清单。

三个坐标系的定义与转换链

标定能否收敛,首先取决于坐标系定义是否唯一、旋转顺序是否固定。本例采用如下约定(右手系,角度单位:度):

  • 船体坐标系 b:原点取船体基准点(如惯导安装中心),xb 指向船艏,yb 指向右舷,zb 指向船底(下正)。
  • 换能器坐标系 s:原点取换能器声学中心,xs 为换能器零方位参考方向,ys 为右舷方向,zs 向下。理想安装时 s 与 b 三轴平行,实际存在小角度偏差。
  • 大地坐标系 n:北-东-地(NED),由惯导姿态与航向给出 b 到 n 的旋转。

旋转顺序统一为 Z-Y-X(航向→俯仰→横滚),即 Rs→b = Rz(Δψ)·Ry(Δθ)·Rx(Δφ)。这个顺序一旦写进标定模型,后续所有观测方程都必须一致,否则三个角会互相“吸收”误差,导致标定结果不可复现。

坐标系转换链:s → b → n,待估偏角位于 s→b 环节 换能器坐标系 s 原点:声学中心 x_s:零方位参考 y_s:右舷 z_s:向下 观测:斜距 r、方位 α、俯仰 β (声学测量在此坐标系表达) 船体坐标系 b 原点:惯导安装中心 x_b:船艏 y_b:右舷 z_b:向下 安装位置矢量 t_s^b (本例按设计值固定,不参与估计) 大地坐标系 n(NED) 北-东-地 由惯导姿态 + 航向给出 R_b→n(已知) 目标大地坐标 / 真值比对 (标定结果的检验基准) R_s→b 待估 R_b→n 已知 待估偏角清单(R_s→b = R_z(Δψ)·R_y(Δθ)·R_x(Δφ),Z-Y-X 顺序) Δψ 航向偏角 绕 z_s 轴,换能器零方位 相对船艏的偏转 典型量级:0.1°~3° 主要激励:多航向航迹 Δθ 俯仰偏角 绕 y_s 轴,换能器声轴 相对水平面的抬/低头 典型量级:0.1°~2° 主要激励:多距离段 Δφ 横滚偏角 绕 x_s 轴,换能器阵 左右倾斜 典型量级:0.1°~2° 主要激励:横摇/多航向

观测方程与待估量的关系

把声学观测写成向量形式,标定模型就清楚了。换能器测得目标在 s 系下的单位方向 us(由 α、β 构造),则目标在 b 系下的方向为 ub = Rs→b·us,再经 Rb→n 进入大地系。标定时通常用“目标大地坐标已知”或“往返一致性”构造残差:

u_s = [cosβ·cosα, cosβ·sinα, sinβ]ᵀ        # 换能器系单位方向
u_b = Rz(Δψ)·Ry(Δθ)·Rx(Δφ) · u_s            # 转到船体系
u_n = R_b→n · u_b                            # 转到大地系
残差 = u_n - u_n_ref                          # 与真值/参考方向比对

三个偏角对残差的敏感方向不同:Δψ 主要改变水平方位,Δθ 主要改变垂直分量,Δφ 在目标位于换能器正前方时几乎不可观测、在侧方目标上才明显。这正是后续航迹设计要“多航向、多距离段”的根本原因——本例只需把这一可观测性要求写进任务边界,不必展开航迹。

坐标系定义一旦确定,应写入标定任务书并冻结:原点、轴向、旋转顺序、角度正方向、单位。任何一项在标定中途变更,都会使已采集数据作废。

标定任务边界清单(可直接填写)

项目本例取值/约定说明
标定对象船体刚性安装USBL换能器不含船体挠曲、不含阵元相位
待估量Δψ、Δθ、Δφ 三个安装偏角Z-Y-X 顺序,单位:度
固定量安装位置矢量 t_s^b按机械安装图取值,不参与估计
已知输入R_b→n、声学斜距/方位/俯仰由惯导与USBL提供
参考基准目标大地坐标或往返一致性用于构造残差
不纳入估计声速、时延、阵元相位、船体挠曲列为误差源,另行处理
验收判据残差RMS与偏角重复性阈值按任务精度要求设定

常见定义陷阱

  • z 轴朝上还是朝下:本例取“下正”。若惯导输出为“上正”,必须在进入标定前统一翻转,否则 Δθ 符号会整体反号。
  • 航向正方向:取顺时针(从上方看)为正,与航海习惯一致;若软件内部用逆时针,需在 Rz 中体现。
  • 旋转顺序不可交换:小角度下 Z-Y-X 与 X-Y-Z 差异很小,但会系统性混入 Δφ 与 Δψ 的耦合,导致重复标定结果漂移。
  • 零方位参考:xs 必须与换能器出厂零方位标记一致,不能凭安装外观臆定。
本例只输出“坐标系定义 + 待估偏角清单 + 观测方程形式”,不涉及具体数值标定结果。任何偏角数值都必须由实际采集数据解算得到,不得以设计值或经验值冒充标定结果。