GPS從入門到放棄（十三） --- 接收機自主完好性監測（RAIM）

GPS從入門到放棄（十三） — 接收機自主完好性監測（RAIM）

接收機自主完好性監測（RAIM: Receiver Autonomous Integrity Monitoring）是根據用戶接收機的冗餘觀測值監測用戶定位結果的完好性，其目的是在導航過程中檢測出發生故障的衛星，並保障導航定位精度。

爲了能進行接收機自主完好性監測，必須有冗餘的觀測量。一般來說，需要可見衛星數5顆以上纔可進行完好性檢測；需要有6顆以上纔可能辨識出故障衛星。

RAIM 的增強版本爲 RAIM-FDE（FDE: Fault Detection Exclusion），即故障檢測與排除技術，因爲需要排除故障衛星，所以必須使用最少6顆可見衛星。

RAIM 算法對於安全性有嚴格要求的應用非常重要，如民航、航空之類。

RAIM 有不同的實現算法：僞距殘差判決法、僞距比較法、校驗向量法以及最大解分離法。

僞距殘差判決法

在定位方程解算中我們講過了定位方程組
$\boldsymbol{G} \left[ \begin{array}{c} \Delta x\\ \Delta y\\ \Delta z\\ c\cdot\Delta \delta_{t} \end{array} \right] = \boldsymbol{b}$

其中 $\boldsymbol{b}$ 爲僞距殘差。在用最小二乘法定位結束後，如果參與定位的衛星數量大於4顆，若定位結果正確，則僞距殘差應該是比較小的；若定位結果有錯誤，則僞距殘差會比較大。於是我們可以根據定位後僞距殘差的大小來判決定位結果是否正確。這種方法就叫僞距殘差判決法。

假設兩次定位之間的間隔時間很短（一般情況下也是事實），那麼接收機在兩次定位之間位置和鐘差變化並不大，於是幾何矩陣 $\boldsymbol{G}$ 可以認爲是不變的。設前一次定位的僞距殘差爲 $\boldsymbol{b}$，本次定位的僞距殘差爲 $\boldsymbol{z}$，則：

$\boldsymbol{z} = \boldsymbol{b} - \boldsymbol{G} \left[ \begin{array}{c} \Delta x\\ \Delta y\\ \Delta z\\ c\cdot\Delta \delta_{t} \end{array} \right] = \boldsymbol{b}-\boldsymbol{G}(\boldsymbol{G}^T\boldsymbol{G})^{-1}\boldsymbol{G}^T\boldsymbol{b} = [\boldsymbol{I} - \boldsymbol{G}(\boldsymbol{G}^T\boldsymbol{G})^{-1}\boldsymbol{G}^T]\boldsymbol{b}$

令 $\boldsymbol{S}=\boldsymbol{I} - \boldsymbol{G}(\boldsymbol{G}^T\boldsymbol{G})^{-1}\boldsymbol{G}^T$，則 $\boldsymbol{z = Sb}$，可知我們無須計算最小二乘解就可以由 $\boldsymbol{b}$ 算得 $\boldsymbol{z}$，然後檢查 $\boldsymbol{z}$ 的各個分量值是否滿足判決條件。

一般來說，會設定某個門限值，若 $\boldsymbol{z}$ 的某個分量值高於門限值，則認爲某個測量值有問題。

也可以計算僞距殘差的平方和：
$SSE = \boldsymbol{z}^T\boldsymbol{z}$

它是一個標量，符合自由度爲 $n-4$ 的 $\chi^2$ 分佈，其中 $n$ 爲觀測量個數。可以根據 $SSE$ 的概率分佈來計算出給定虛警率下的檢測門限值，然後根據這個門限值去判決是否存在錯誤觀測量。

僞距比較法

僞距比較法是直接利用僞距來判決。步驟如下：

1. 在觀測數量大於4的情況下，對 $n$ 個觀測量進行分組，一組爲4個，另一組爲 $n-4$ 個。
2. 對4個的那一組進行最小二乘法解算出接收機位置。
3. 對 $n-4$ 個的那一組計算衛星位置。
4. 根據解算出的接收機位置和衛星位置估算僞距。
5. 將估算的僞距與實際觀測的僞距進行比較。若比較結果都很接近，說明觀測結果沒問題，第2步中解算的位置也沒問題。否則說明有問題。
6. 對有問題的情況進行處理。得看具體情況。若只有一個比較結果差異很大，則很大概率是那一個觀測值有問題；若很多比較結果都差異很大，則很大概率是最小二乘解有問題。若有必要，可以重新分組，重新計算比較判定。

校驗向量法

此方法是構造出一個校驗矩陣 $\boldsymbol{P}$，然後計算出校驗向量 $\boldsymbol{p} = \boldsymbol{Pb}$。判決方法是用 $\boldsymbol{p}$ 向量的模來判決是否出現了觀測量錯誤，然後用 $\boldsymbol{p}$ 向量的元素和原點之間的連線的斜率來判決哪一個僞距觀測值出錯的概率最大。

$\boldsymbol{P}$ 矩陣的選取並不唯一，可以通過對幾何矩陣 $\boldsymbol{G}$ 做 $\boldsymbol{QR}$ 分解來得到，即若 $\boldsymbol{G = QR}$，則 $\boldsymbol{Q}$ 的轉置矩陣的底部 $n-4$ 行就是 $\boldsymbol{P}$ 矩陣。

最大解分離法

其基本思想是將 $N (N \ge 5)$ 個僞距觀測量分成 $N$ 組，每一組不包含第 i 個僞距觀測量，$i=1,2,\cdots,N$。則每組包含 $N-1$ 個觀測量，可以用最小二乘法求得一個解，於是可以得到 $N$ 個解。計算這 $N$ 個解之間的兩兩距離，若距離都很接近，說明全部觀測量沒有錯誤；反之，則至少有一個觀測量有錯誤。通過求出兩兩距離的最大值來跟一個預設的門限值比較來判定是否有觀測量出錯。

這個方法的缺點是計算量很大，而且出錯時不容易排除錯誤。