無跡卡曼濾波器

**無跡卡曼濾波器**（）係一種卡曼濾波器，攞嚟解決種問題係擴展卡曼濾波器喺高度非線性系統模型下表現嘸佳噉嘅情況；呢種卡曼濾波器使到一種叫做無跡變換（UT）嘅確定性抽樣技巧，即喺平均值周圍揀一組最細嘅樣本點（喊做sigma點）。然之後，啲sigma點識透過非線性函數傳播（propagated），從中形成到個新估計值畀平均值同埋協方差。噉樣可以避免似喺EKF當中噉經由線性化個潛在嘅非線性模型嚟傳播啲協方差。係噉呢種濾波器嘅具體效能取決於點樣計UT啲轉換統計量同埋使到邊一組採樣到嘅sigma點。需要注意嘅係，啲新嘅UKF始終可以以一致（consistent）嘅方法整出。對於某啲系統，噉樣得到嘅UKF更加準確噉估計得到真實嘅平均值同協方差，個論斷可以透過蒙地卡羅抽樣或者啲後驗統計值個泰勒級數展開嚟驗證到。另外，種技術免走明式噉計算Jacobi矩陣嘅要求：對於啲複雜函數來講，一嘸係啲函數嘸微分得所以冇辦法計得出Jacobi矩陣、就係計算本身好難無論係使分析法（導數複雜）抑或係數值法（計算成本高）都好。 

## Sigma點
對於隨機向量\mathbf{x}=(x_1, \dots, x_L)，啲Sigma點係任何啲向量嘅集：

:  \{\mathbf{s}_0,\dots, \mathbf{s}_N \}=\bigl\{\begin{pmatrix} s_{0,1}& s_{0,2}&\ldots& s_{0,L} \end{pmatrix}, \dots, \begin{pmatrix} s_{N,1}& s_{N,2}&\ldots& s_{N,L} \end{pmatrix}\bigr\}

有：

* 一階權重W_0^a, \dots, W_N^a滿足

#  \sum_{j=0}^N W_j^a=1 
# 對所有啲i=1, \dots, L有： E[x_i]=\sum_{j=0}^N W_j^a s_{j,i} 

* 二階權重W_0^c, \dots, W_N^c滿足

#  \sum_{j=0}^N W_j^c=1 
# 對所有啲𠵿 (i,l) \in \{1,\dots, L\}^2: E[x_ix_l]=\sum_{j=0}^N W_j^c s_{j,i}s_{j,l} 

UKF演算法入便一種揀啲sigma點同埋權重畀\mathbf{x}_{k-1\mid k-1}嘅簡單方式係：

: \begin{align}
\mathbf{s}_0&=\hat \mathbf{x}_{k-1\mid k-1}\\
-1&

其中\hat \mathbf{x}_{k-1\mid k-1}係 \mathbf{x}_{k-1\mid k-1}個平均估計。枚向量\mathbf{A}_j係\mathbf{A}個第*j*棟，其中\mathbf{A}有\mathbf{P}_{k-1\mid k-1}=\mathbf{AA}^\textsf{T}。矩陣 \mathbf{A} 應使數值上高效又穩定嘅方法似Cholesky分解嚟計算。平均值嘅權重W_0可以任意噉揀。

第個流行嘅參數化（對上述普適化）方式係： 

: \begin{align}
\mathbf{s}_0&=\hat \mathbf{x}_{k-1\mid k-1}\\
W_0^a&= \frac{\alpha^2\kappa-L}{\alpha^2\kappa}\\
W_0^c&= W_0^a + 1-\alpha^2+\beta \\
\mathbf{s}_j&=\hat \mathbf{x}_{k-1\mid k-1} + \alpha\sqrt{\kappa} \mathbf{A}_j, \quad j=1, \dots, L\\
\mathbf{s}_{L+j}&=\hat \mathbf{x}_{k-1\mid k-1} - \alpha\sqrt{\kappa} \mathbf{A}_j, \quad j=1, \dots, L\\
W_j^a&=W_j^c=\frac{1}{2\alpha^2\kappa}, \quad j=1, \dots, 2L.
\end{align}


\alpha同\kappa啲sigma點嘅擴散。\beta同x嘅分佈有關。

啱嘅值取決於手頭係咩樣嘅問題，而典型建議係揀\alpha = 10^{-3}，\kappa = 1同埋\beta = 2。之不過，大啲嘅值畀\alpha（似\alpha = 1）可能捕捉返個分佈嘅擴散同埋啲可能有嘅非線性會好啲。若果個x真實分佈係高斯分佈，係噉揀\beta = 2係最啱嘅。

## 過程

### 預測
同EKF一樣，UKF預測可以獨立於UKF更新，可以佮線性（或者EKF）更新一齊用或者反過嚟。攞有估計值 \hat\mathbf{x}_{k-1\mid k-1}、\mathbf{P}_{k-1\mid k-1}畀平均值同協方差嘅話，可以攞到 N = 2L+1 粒sigma點，似上節所講。啲sigma點係透過轉換函數*f* （transition function）傳播嘅：

: \mathbf{x}_{j} = f\left(\mathbf{s}_{j}\right) \quad j = 0, \dots, 2L .

着傳播過嘅 sigma 點又着加權嚟產生啲預測嘅平均值同埋協方差。

: \begin{align}
   \hat{\mathbf{x}}_{k \mid k-1} &= \sum_{j=0}^{2L} W_j^a \mathbf{x}_j \\
         \mathbf{P}_{k \mid k-1} &= \sum_{j=0}^{2L} W_j^c \left(\mathbf{x}_j - \hat{\mathbf{x}}_{k \mid k-1}\right)\left(\mathbf{x}_j - \hat{\mathbf{x}}_{k \mid k-1}\right)^\textsf{T}+\mathbf{Q}_k
\end{align}

其中W_j^a同埋W_j^c分別係一階同埋二階權重喺啲原始sigma點上。矩陣 \mathbf{Q}_k 係轉換噪音\mathbf{w}_k嘅協方差。

### 更新
畀有預測估計\hat{\mathbf{x}}_{k \mid k-1}同\mathbf{P}_{k \mid k-1}，計得出新一組N = 2L+1粒sigma點\mathbf{s}_0, \dots, \mathbf{s}_{2L}有相應一階權重 W_0^a,\dots W_{2L}^a二階權重W_0^c,\dots, W_{2L}^c嘅。 啲sigma點通過h轉變：

:  \mathbf{z}_j=h(\mathbf{s}_j), \,\, j=0,1, \dots, 2L .

係噉轉變過嘅啲點嘅觀測平均值同埋協方差就噉計出：

: \begin{align}
  \hat{\mathbf{z}} &= \sum_{j=0}^{2L} W_j^a \mathbf{z}_j \\[6pt]
  \hat{\mathbf{S}}_k &= \sum_{j=0}^{2L} W_j^c (\mathbf{z}_j-\hat{\mathbf{z}})(\mathbf{z}_j-\hat{\mathbf{z}})^\textsf{T} + \mathbf{R}_k
\end{align}

其中\mathbf{R}_k係觀測噪音\mathbf{v}_k嘅協方差矩陣。另外仲要到交叉協方差（cross covariance ）矩陣：

: \begin{align}
  \mathbf{C_{sz}} &= \sum_{j=0}^{2L} W_j^c (\mathbf{s}_j-\hat\mathbf{x}_{k|k-1})(\mathbf{z}_j-\hat\mathbf{z})^\textsf{T}
\end{align}

其中\mathbf{s}_j係攞\hat{\mathbf{x}}_{k \mid k-1}、\mathbf{P}_{k \mid k-1}整出嘅啲未轉變過嘅sigma點。

個卡曼增益即係：

: \begin{align}
  \mathbf{K}_k=\mathbf{C_{sz}}\hat{\mathbf{S}}_k^{-1}.
\end{align}

更新後嘅估計值畀平均值同埋協方差即係：

: 
\begin{align}
\hat\mathbf{x}_{k\mid k}&=\hat\mathbf{x}_{k|k-1}+\mathbf{K}_k(\mathbf{z}_k-\hat\mathbf{z})\\
\mathbf{P}_{k\mid k}&=\mathbf{P}_{k\mid k-1}-\mathbf{K}_k\hat{\mathbf{S}}_k\mathbf{K}_k^\textsf{T}.
\end{align}


## 考