学习笔记——RAIM中的误差投影

RAIM中的误差投影

一、从投影算子的数学本质说起

1.1 回顾两个核心投影矩阵

在最小二乘残差法中，我们有两个投影矩阵：

帽子矩阵 H（投影到列空间）：

\[\mathbf{H} = \mathbf{G}(\mathbf{G}^T\mathbf{G})^{-1}\mathbf{G}^T \]

残差投影矩阵 P（投影到左零空间）：

\[\mathbf{P} = \mathbf{I} - \mathbf{H} \]

关键性质：

• \(\mathbf{H}\) 和 \(\mathbf{P}\) 都是投影矩阵：\(\mathbf{H}^2 = \mathbf{H}\)，\(\mathbf{P}^2 = \mathbf{P}\)
• 它们正交：\(\mathbf{H}\mathbf{P} = \mathbf{0}\)
• 它们互补：\(\mathbf{H} + \mathbf{P} = \mathbf{I}\)

1.2 误差的分解

任何测量值 \(\mathbf{y}\) 都可以分解为：

\[\mathbf{y} = \underbrace{\mathbf{H}\mathbf{y}}_{\text{可解释部分（位置相关）}} + \underbrace{\mathbf{P}\mathbf{y}}_{\text{残差部分（位置无关）}} \]

现在，考虑测量值 \(\mathbf{y}\) 包含真实值 \(\mathbf{Gx}\) 和误差 \(\boldsymbol{\epsilon}\)：

\[\mathbf{y} = \mathbf{Gx} + \boldsymbol{\epsilon} \]

代入分解：

\[\mathbf{H}\mathbf{y} = \mathbf{H}(\mathbf{Gx} + \boldsymbol{\epsilon}) = \mathbf{Gx} + \mathbf{H}\boldsymbol{\epsilon} \]

\[\mathbf{P}\mathbf{y} = \mathbf{P}(\mathbf{Gx} + \boldsymbol{\epsilon}) = \mathbf{0} + \mathbf{P}\boldsymbol{\epsilon} \]

这就是核心发现：

• 位置信息 \(\mathbf{Gx}\) 完全在 \(\mathbf{H}\) 的像空间中
• 残差 \(\mathbf{w} = \mathbf{P}\boldsymbol{\epsilon}\) 只包含误差通过 \(\mathbf{P}\) 的投影
• 任何误差，无论来源，都被投影到了两个空间

---

二、误差的投影机制详解

2.1 误差向量的分解

设总误差向量 \(\boldsymbol{\epsilon}\)（n×1）包含：

• 测量噪声
• 多径误差
• 卫星钟差残差
• 轨道误差
• 电离层残差（单频）
• 等等

这些误差可以被分解为两个正交分量：

\[\boldsymbol{\epsilon} = \underbrace{\mathbf{H}\boldsymbol{\epsilon}}_{\text{影响位置的误差}} + \underbrace{\mathbf{P}\boldsymbol{\epsilon}}_{\text{不直接影响位置的误差}} \]

通俗理解：

• \(\mathbf{H}\boldsymbol{\epsilon}\)：误差中与卫星几何结构一致的部分，会污染定位结果
• \(\mathbf{P}\boldsymbol{\epsilon}\)：误差中与几何结构正交的部分，不影响定位，但会在残差中显现

2.2 一个具体的二维例子

为了直观，我们降到二维（n=3颗卫星，解2个参数）：

设 \(\mathbf{G}\) 矩阵（3×2）：

\[\mathbf{G} = \begin{bmatrix} 1 & 0 \\ 0 & 1 \\ 1 & 1 \end{bmatrix} \]

计算 \(\mathbf{H}\) 和 \(\mathbf{P}\)：

\[\mathbf{G}^T\mathbf{G} = \begin{bmatrix} 1 & 0 & 1 \\ 0 & 1 & 1 \end{bmatrix} \begin{bmatrix} 1 & 0 \\ 0 & 1 \\ 1 & 1 \end{bmatrix} = \begin{bmatrix} 2 & 1 \\ 1 & 2 \end{bmatrix} \]

\[(\mathbf{G}^T\mathbf{G})^{-1} = \frac{1}{3}\begin{bmatrix} 2 & -1 \\ -1 & 2 \end{bmatrix} \]

\[\mathbf{H} = \mathbf{G}(\mathbf{G}^T\mathbf{G})^{-1}\mathbf{G}^T = \frac{1}{3}\begin{bmatrix} 2 & -1 & 1 \\ -1 & 2 & 1 \\ 1 & 1 & 2 \end{bmatrix} \]

\[\mathbf{P} = \mathbf{I} - \mathbf{H} = \frac{1}{3}\begin{bmatrix} 1 & 1 & -1 \\ 1 & 1 & -1 \\ -1 & -1 & 1 \end{bmatrix} \]

验证性质：

• \(\mathbf{H}\mathbf{G} = \mathbf{G}\)：列空间投影不变性
• \(\mathbf{P}\mathbf{G} = \mathbf{0}\)：左零空间投影为零

现在看误差向量的投影：

假设误差向量 \(\boldsymbol{\epsilon} = [1, 0, 0]^T\)（只有卫星1有误差）：

\[\mathbf{H}\boldsymbol{\epsilon} = \frac{1}{3}[2, -1, 1]^T \]

\[\mathbf{P}\boldsymbol{\epsilon} = \frac{1}{3}[1, 1, -1]^T \]

这意味着：卫星1的1单位误差中：

• 2/3 被投影到位置解（污染定位）
• 1/3 进入残差（可被RAIM检测）

---