2. 磁场与磁测量基础

三轴磁力计的本质，是把空间某一点的磁场矢量分解到三个互相垂直的敏感轴上，再通过标定还原出真实的磁场方向和大小。因此，在讨论传感器结构和标定算法之前，必须先建立对“磁场是什么、地磁场长什么样、如何用数学描述它”的统一认识。本节围绕地磁场特性、磁场矢量表示、磁感应强度单位以及坐标系定义四个基础问题展开，为后续的零偏、灵敏度、非正交和硬磁/软磁标定打下概念地基。

一、地磁场的基本特性

地球本身相当于一个巨大的磁偶极子，其磁力线从南半球附近发出、进入北半球附近，形成近似偶极子的空间分布。对于地面附近的磁测量应用，我们通常把地磁场视为一个缓慢变化的准静态矢量场，它具有以下几个对工程实现非常关键的特性：

• 强度微弱：地表总场强大约在 25 μT 到 65 μT 之间，仅为普通冰箱贴表面磁场的千分之一量级，因此对传感器噪声和干扰非常敏感。
• 方向稳定：在局部区域内，地磁场的水平分量始终指向磁北，这是电子罗盘、航向解算的物理基础。
• 随位置变化：不同经纬度、不同海拔下，总场强和磁倾角都会变化，全球范围内需要用球谐模型（如 WMM、IGRF）描述。
• 随时间缓慢漂移：地磁极本身在缓慢移动，长期应用中需要引入年变率修正。
• 易受环境干扰：铁磁材料、电机、扬声器、电流回路都会产生局部附加磁场，其强度往往远大于地磁信号。

理解这些特性后可以得出一个工程结论：磁力计测量的从来不是“地磁场本身”，而是地磁场与所有环境干扰磁场的叠加。标定的目标，正是把传感器自身的误差和可建模的干扰从测量值中剥离出来。

二、磁场矢量表示与磁感应强度单位

磁场在空间中是一个矢量，记作 B，在直角坐标系中可写成三分量形式：

B = [Bx, By, Bz]^T
|B| = sqrt(Bx^2 + By^2 + Bz^2)

其中 |B| 为总场强，方向由三个分量共同决定。在导航和姿态解算中，常用的还有水平分量、磁倾角等派生量：

水平分量  Bh = sqrt(Bx^2 + By^2)
磁倾角    I  = atan2(Bz, Bh)
磁偏角    D  = atan2(By, Bx)   // 相对磁北的方位

磁感应强度的国际单位是特斯拉（T），地磁测量中常用微特（μT）或纳特（nT）。历史上还广泛使用高斯（G）和毫高斯（mG），二者换算关系必须牢记：

三、坐标系定义：从地理系到传感器系

磁测量中至少涉及三套坐标系，混淆它们是初学者最常见的错误来源：

1. 地理坐标系（NED / ENU）：以当地水平面为基准，N 指北、E 指东、D 指地（或 U 指天）。地磁场在此系中的方向由磁偏角和磁倾角确定。
2. 载体坐标系（Body）：固连在设备上，通常定义 X 向前、Y 向右、Z 向下，随设备姿态变化。
3. 传感器坐标系（Sensor）：由三轴磁力计三个敏感轴的实际安装方向决定，理想情况下与载体系重合，实际存在安装误差和非正交误差。

三者的关系可以概括为：地理系中的地磁矢量，经过姿态旋转矩阵 R 变换到载体系，再经过传感器误差矩阵 M 和零偏 b 变换，最终得到磁力计读数：

B_meas = M * R * B_geo + b + n
其中：
  B_geo  : 地理系下的地磁矢量
  R      : 姿态旋转矩阵（载体相对地理系）
  M      : 3x3 综合误差矩阵（灵敏度、非正交、软磁）
  b      : 零偏与硬磁偏置
  n      : 噪声

这个模型是后续所有标定方法的出发点。标定要估计的，正是 M 和 b 这两组参数。

四、典型场景与常见误区

假设一台手持设备在室外空旷处测量，读数为 Bx=22.1 μT、By=-8.4 μT、Bz=41.3 μT。计算总场强约为 47.6 μT，处于地磁正常范围；水平分量约 23.6 μT，磁倾角约 60°，符合中纬度地区特征。若把设备放到铁桌上，总场强突然跳到 120 μT，说明引入了强硬磁干扰，此时任何航向解算都不可信。

另一个高频误区是忽略坐标系一致性：用地理系公式去解释传感器系数据，或在不同坐标系之间直接比较分量大小。正确做法是先明确每个矢量的参考系，再通过旋转矩阵统一到同一坐标系下运算。

把上述内容串起来看：地磁场是一个微弱、方向稳定但随位置和时间变化的矢量场；磁力计读数是在传感器坐标系下对该矢量叠加误差后的观测结果；单位换算和坐标系定义则是所有计算的前提。只有先厘清这三层关系，后续的椭球拟合、最小二乘标定和航向解算才有正确的物理意义。