粒子濾波器（英語：particle filter）是一種使用蒙特卡羅方法的遞迴濾波器，透過一組具有權重的隨機樣本（粒子）來表示隨機事件的後驗機率，從含有雜訊或不完整的觀測序列，估計出動態系統的狀態，粒子濾波器可以運用在任何狀態空間的模型上。粒子濾波器是卡爾曼濾波器的一般化方法，卡爾曼濾波器建立在線性的狀態空間和高斯分布的雜訊上；而粒子濾波器的狀態空間模型可以是非線性，且雜訊分布可以是任何型式。



歷史

啟發式演算法

從統計和機率的角度來看，粒子濾波器屬於遺傳算法分支過程的類別，並且是指平均場類型相互作用的粒子方法。 這些粒子方法的解釋取決於科學專業。 在進化計算中，平均場遺傳類型粒子方法通常用作啟發式和自然搜索算法（也稱為元啟發算法）。在計算物理學和分子化學中，它們用於解決費曼 - 卡茲 (Feynman-Kac) 路徑積分問題，或者用於計算玻爾茲曼 - 吉布斯 (Boltzmann-Gibbs) 度量，薛丁格算子的最高特徵值和基態。在生物學與遺傳學中，它們還代表了某些環境中一群個體或基因的進化。

平均場類型進化計算技術的起源可以追溯到 1950 年和 1954 年，這是艾倫・圖靈在遺傳類型突變選擇學習機上的開創性工作，以及紐澤西州普林斯頓高級研究所的尼爾斯・艾爾・巴里切利 (Nils Aall Barricelli) 的文章。統計學中的粒子濾波器最早可追溯到 1950 年代中期。哈默斯利 (Hammersley) 等人在 1954 年提出的「窮人的蒙地卡羅」法暗示了當今使用的遺傳類型粒子過濾方法。 1963 年，尼爾斯・艾爾・巴里切利 (Nils Aall Barricelli) 模擬了一種遺傳類型算法，以模仿個人玩簡單遊戲的能力。在進化計算文獻中，遺傳類型突變選擇算法通過 1970 年代初期約翰・霍蘭德 (John Holland) 的開創性工作而流行，特別是 1975 年出版的著作。

在生物學與遺傳學中，澳大利亞遺傳學家亞歷克斯・弗雷澤（Alex Fraser）在 1957 年還發表了一系列有關人工選擇生物的遺傳類型模擬的論文。 由生物學家進行的進化計算機模擬在 1960 年代初變得更加普遍，並且該方法在弗雷澤和伯內爾 (Burnell)（1970）和克羅斯比 (Crosby)（1973）的書中進行了描述。弗雷澤的模擬包括了現代突變選擇遺傳粒子算法的所有基本要素。

從數學觀點來看，給定一些部分和雜訊觀測結果的信號隨機狀態的條件分布，是由信號的隨機軌跡上的費曼 - 卡茲 (Feynman-Kac) 機率描述的，該機率由一系列概似勢函數加權。量子蒙地卡羅法，更具體地說是擴散蒙地卡羅法，也可以解釋為費曼 - 卡茲路徑積分的平均場遺傳類型粒子近似。量子蒙特卡羅方法的起源通常歸因於恩里科・費米 (Enrico Fermi) 和羅伯特・里克特邁耶 (Robert Richtmyer)，他們於 1948 年開發了一種中子鏈反應的均值場粒子解釋方法，但它是第一個類似啟發式和遺傳類型的粒子算法（又名「重採樣」或「重新配置」蒙地卡羅方法） ）是在 1984 年由 Jack H. Hetherington 提出的用於估計量子系統（在簡化的矩陣模型中）的基態能量的方法。還可以引用 1951 年出版的西奧多・愛德華・哈里斯 (Theodore E. Harris) 和赫爾曼・康恩 (Herman Kahn) 的開創性著作，該著作使用平均場，但類似於啟發式遺傳方法，用於估計粒子傳輸能量。在分子化學中，遺傳啟發式粒子方法（又名修剪和富集策略）的使用可以追溯到 1955 年馬歇爾（Marshall）的開創性工作。羅森布魯斯（N. Rosenbluth）和阿里安娜・W・羅森布盧斯 (Arianna. W. Rosenbluth)。

遺傳粒子算法在高級信號處理和貝葉斯推理中的使用是最近的。1993 年 1 月，北川源四郎開發了「蒙特卡羅濾波器」，該文在 1996 年出現了輕微修改。1993 年 4 月，戈登等人在他們的開創性著作中發表了遺傳類型算法在貝氏統計推斷中的應用。作者將其算法命名為「自舉濾波器」，並證明與其他濾波方法相比，其自舉算法不需要對該狀態空間或系統噪聲進行任何假設。獨立的是皮埃爾・德爾・莫拉爾 (Pierre Del Moral) 和希米爾孔・卡瓦略 (Himilcon Carvalho)，皮埃爾・德爾・莫拉爾 (Pierre Del Moral)，安德烈・莫寧 (André Monin) 和熱拉爾・薩魯特 (Gérard Salut) 在 1990 年代中期發表的有關粒子過濾器的論文。 1989-1992 年初，皮埃爾・德爾・莫拉爾 (P.Del Moral)，J.C. 諾耶 (J.C. Noyer)，G.Rigal 和 G.Salut 在 LAAS-CNRS 中也開發了信號處理技術，並通過 STCAN（服務技術）進行了一系列受限和分類研究 des Constructions and Armes Navales），IT 公司 DIGILOG 和 LAAS-CNRS（系統分析和體系結構實驗室）來解決 RADAR / SONAR 和 GPS 信號處理問題。

數學基礎

從 1950 年到 1996 年，所有有關粒子過濾器，遺傳算法的出版物，包括在計算物理和分子化學中引入的修剪和重採樣蒙地卡羅方法，都提出了適用於不同情況的類似於自然和啟發式算法，而沒有任何證據證明它
們的一致性，也沒有討論估計的偏差以及基於族譜和祖先樹的算法。

這些粒子算法的數學基礎和首次嚴格的分析是由皮埃爾・德爾・莫拉爾 (Pierre Del Moral) 在 1996 年提出的。本文還包含了概似函數和未標準化的條件機率測度的粒子近似的不偏性質的證明。本文介紹的概似函數的不偏粒子估計器如今已用於貝式統計推理中。

到 1990 年代末，Dan Crisan，Jessica Gaines 和 Terry Lyons 以及 Dan Crisan，Pierre Del Moral 和 Terry Lyons 也開發了具有不同種群大小的分支類型粒子方法。 P. Del Moral，A。Guionnet 和 L. Miclo 在
2000 年開發了該領域的進一步發展。第一個中央極限定理是由 Pierre Del Moral 和 Alice Guionnet 於 1999 年以及 Pierre Del Moral 和 Laurent Miclo 於 2000 年提出的。Pierre 於 1990 年代末提出了關於粒子過濾器
時間參數的第一個統一收斂結果。 Del Moral 和 Alice Guionnet。 2001 年，P。Del Moral 和 L. Miclo 首次對基於族譜樹的粒子濾波器平滑器進行了嚴格的分析。

關於 Feynman-Kac 粒子方法論的理論和相關的粒子過濾器算法已在 2000 年和 2004 年的書中得到發展。這些抽象的機率模型封裝了遺傳類型算法，粒子和自舉濾波器，交互的卡爾曼濾波器（又稱 Rao–Blackwellized
粒子濾波器），重要性採樣和重採樣樣式的粒子濾波器技術，包括基於族譜樹的方法和用於解決濾波和平滑問題的粒子後向方法。其他類別的粒子濾波方法包括基於族譜樹的模型，後向馬可夫粒子模型，自
適應平均場粒子模型，島型粒子模型和粒子馬可夫鏈蒙特卡羅方法。



濾波問題

目標

粒子濾波器的目標是通過觀測值估計狀態變數的事後機率。粒子濾波器是針對隱藏式馬可夫模型設計的，系統中既包含隱藏的變數，也包含可觀測的變數。可觀測變數（觀測過程）與隱藏變數（狀態過程）通過某種已知的函數關係相關聯。類似地，描述狀態變數演化的動態系統也是機率學上已知的。

一個通用粒子濾波器通過觀察測量過程估計隱藏狀態的事後分布。考察下圖所示的狀態空間：

　
: ${\begin{array}{cccccccccc}X_{0}&\to &X_{1}&\to &X_{2}&\to &X_{3}&\to &\cdots &{\text{signal}}\\\downarrow &&\downarrow &&\downarrow &&\downarrow &&\cdots &\\Y_{0}&&Y_{1}&&Y_{2}&&Y_{3}&&\cdots &{\text{observation}}\end{array}}$

濾波問題是要通過給定的觀測過程 $Y_{0},\cdots ,Y_{k},$ 在任意時刻 _k_ 的值，序列地估計隱藏狀態 $X_{k}$ 的值。

所有的對 $X_{k}$ 的貝氏估計服從事後分布 _p_(_x_k | _y_0,_y_1,…,_y_k) 。粒子濾波器方法使用經驗測量和通用粒子算法，提供了這些條件機率的近似值。與之對比，馬可夫鏈蒙地卡羅（MCMC）或重要性採樣方式會對整個事後機率 _p_(_x_0,_x_1,…,_x_k | _y_0,_y_1,…,_y_k) 進行建模。

信號 - 觀測模型

粒子濾波器能夠從一系列含有雜訊或不完整的觀測值中，估計出動態系統的內部狀態。在動態系統的分析中，需要兩個模型，一個用來描述狀態隨時間的變化 (系統模型)，另一個用來描述每個狀態下觀測到的雜訊 (觀測模型)，將這兩個模型都用機率來表示。

在許多情況下，每得到一個新的觀測值時，都必須對系統做出一次估計，利用遞迴濾波器，能夠有效地達到這樣的目的。遞迴濾波器會對得到的資料做連續處理，而非分批處理，因此不需要將完整的資料儲存起來，也不需要在得到新的觀測值時，將現有的資料重新做處理。遞迴濾波器包含兩個步驟：

預測：利用系統模型，由前一個狀態的資訊預測下一個狀態的機率密度函數。

更新：利用最新的觀測值，修改預測出的機率密度函數。

藉由貝氏推論 (Baysian inference)，我們可以描述出狀態空間的機率，並在得到新的觀測值時，對系統做出更新，因而達成上述目的。

近似貝式計算模型

在某些問題中，在信號的隨機狀態下，觀測值的條件分布可能不具有密度，或者不可能或太複雜而無法計算。在這種狀況下，我們需要求近似值。其中一個策略是藉由馬可夫鏈 ${\mathcal {X}}_{k}=\left(X_{k},Y_{k}\right)$ 來取代信號 $X_{k}$ ，採用虛擬觀測的形式

　
: ${\mathcal {Y}}_{k}=Y_{k}+\epsilon {\mathcal {V}}_{k}\quad {\mbox{for some parameter}}\quad \epsilon \in [0,1]$

對於具有已知機率密度函數的獨立序列中的某些序列，核心的概念是觀測此式子 (公式)

　
: ${\text{Law}}\left(X_{k}|{\mathcal {Y}}_{0}=y_{0},\cdots ,{\mathcal {Y}}_{k}=y_{k}\right)\approx _{\epsilon \downarrow 0}{\text{Law}}\left(X_{k}|Y_{0}=y_{0},\cdots ,Y_{k}=y_{k}\right)$

在部分的觀測 ${\mathcal {Y}}_{0}=y_{0},\cdots ,{\mathcal {Y}}_{k}=y_{k},$ 之下，粒子濾波器是和馬可夫鏈 ${\mathcal {X}}_{k}=\left(X_{k},Y_{k}\right)$ 有相關的，定義為根據演化出的 $\mathbb {R} ^{d_{x}+d_{y}}$ ，在一些明顯的濫用性符號 $p({\mathcal {Y}}_{k}|{\mathcal {X}}_{k})$ 下，帶有概似函數的粒子。這些機率性的技術跟近似貝式計算非常相關。在量子濾波器中，這些近似貝式量子濾波器的技術在 1998 年被皮埃爾・德爾・莫拉爾 (P. Del Moral), 讓・雅各德 (J. Jacod) 和菲力普・普羅特 (P. Protter) 所採用。且被皮埃爾・德爾・莫拉爾，阿諾・杜斯 (A. Doucet) 和 阿傑・賈斯拉 (A. Jasra) 更進一步發展。

非線性濾波方程式

藉由貝式定理在條件機率中可以得出

　
: $p(x_{0},\cdots ,x_{k}|y_{0},\cdots ,y_{k})={\frac {p(y_{0},\cdots ,y_{k}|x_{0},\cdots ,x_{k})p(x_{0},\cdots ,x_{k})}{p(y_{0},\cdots ,y_{k})}}$

當

　
: ${\begin{aligned}p(y_{0},\cdots ,y_{k})&=\int p(y_{0},\cdots ,y_{k}|x_{0},\cdots ,x_{k})p(x_{0},\cdots ,x_{k})dx_{0}\cdots dx_{k}\\p(y_{0},\cdots ,y_{k}|x_{0},\cdots ,x_{k})&=\prod _{l=0}^{k}p(y_{l}|x_{l})\\p(x_{0},\cdots ,x_{k})&=p_{0}(x_{0})\prod _{l=1}^{k}p(x_{l}|x_{l-1})\end{aligned}}$

粒子濾波器也是一個近似值，當有足夠的粒子時，結果會更加的準確。非線性濾波方程式可以藉由遞迴得出

　
:

依照慣例 $p(x_{0}|y_{0},\cdots ,y_{k-1})=p(x_{0})$ 在 k=0 時。非線性濾波方程式問題的關鍵在於依序計算這些條件機率分布。

費曼 - 卡茲公式

固定一個時間範圍和一個序列的觀測 $Y_{0}=y_{0},\cdots ,Y_{n}=y_{n}$ ，在 k=0....n，我們定義:

　
: $G_{k}(x_{k})=p(y_{k}|x_{k}).$

在這種表示法中，對於從原點 k = 0 到時間 k = n 的 軌跡集合上的任何有界函數 F，我們都可以用費曼 - 卡茲公式來表示

　
: ${\begin{aligned}\int F(x_{0},\cdots ,x_{n})p(x_{0},\cdots ,x_{n}|y_{0},\cdots ,y_{n})dx_{0}\cdots dx_{n}&={\frac {\int F(x_{0},\cdots ,x_{n})\left\{\prod \limits _{k=0}^{n}p(y_{k}|x_{k})\right\}p(x_{0},\cdots ,x_{n})dx_{0}\cdots dx_{n}}{\int \left\{\prod \limits _{k=0}^{n}p(y_{k}|x_{k})\right\}p(x_{0},\cdots ,x_{n})dx_{0}\cdots dx_{n}}}\\&={\frac {E\left(F(X_{0},\cdots ,X_{n})\prod \limits _{k=0}^{n}G_{k}(X_{k})\right)}{E\left(\prod \limits _{k=0}^{n}G_{k}(X_{k})\right)}}\end{aligned}}$

這些費曼 - 卡茲路徑積分模型出現在各種科學學科中，包括計算物理，生物學，資訊理論和計算機科學。他們的解釋取決於應用領域。舉例來說，如果我們選擇狀態空間某些子集的指標函數 $G_{n}(x_{n})=1_{A}(x_{n})$ ，它們代表馬可夫鏈在給定管中的條件分布。我們會得到:

　
: $E\left(F(X_{0},\cdots ,X_{n})|X_{0}\in A,\cdots ,X_{n}\in A\right)={\frac {E\left(F(X_{0},\cdots ,X_{n})\prod \limits _{k=0}^{n}G_{k}(X_{k})\right)}{E\left(\prod \limits _{k=0}^{n}G_{k}(X_{k})\right)}}$

和

　
: $P\left(X_{0}\in A,\cdots ,X_{n}\in A\right)=E\left(\prod \limits _{k=0}^{n}G_{k}(X_{k})\right)$

只要標準化常數嚴格為正。



非線性貝葉斯追蹤

在動態系統中，我們常常需要對物體做追蹤。假設有一動態系統，已知其狀態序列 $\{x_{k},k\in N\}$ 定義為

$$
x_{k}=f_{k}(x_{k-1},v_{k-1})
$$

其中 $f_{k}:R^{n_{x}}\times R^{n_{v}}\rightarrow R^{n_{x}}$ 是狀態轉移函數，可以是非線性的函數， $\{v_{k-1},k\in N\}$ 是一個獨立且相同分布 (i.i.d.) 的過程雜訊序列， $n_{x}$ 和 $n_{v}$ 分別代表狀態和過程雜訊向量的維度， $N$ 為自然數的集合。追蹤的目標是要遞迴地從觀測值 $y_{k}$ 估計出 $x_{k}$ ，而觀測值 $y_{k}$ 定義為

$$
y_{k}=h_{k}(x_{k},n_{k})
$$

其中 $h_{k}:R^{n_{x}}\times R^{n_{n}}\rightarrow R^{n_{y}}$ 是觀測函數，可以是非線性的函數， $\{n_{k},k\in N\}$ 是一個獨立且相同分布 (i.i.d.) 的觀測雜訊序列， $n_{y}$ 和 $n_{n}$ 分別代表觀測值和觀測雜訊向量的維度。

從貝氏定理的觀點來看，追蹤問題是要根據到時間 $k$ 為止的已知資訊 $y_{1:k}$ ，遞迴地計算出時間 $k$ 的狀態 $x_{k}$ 的可信度，這個可信度用機率密度函數 $p(x_{k}\mid y_{1:k})$ 來表示。假設初始機率密度函數 $p(x_{0}\mid y_{0})\equiv p(x_{0})$ (稱為先驗機率) 為已知，則利用 「預測」 和 「更新」 兩個步驟就能遞迴地計算出機率密度函數 $p(x_{k}\mid y_{1:k})$ 。

在此假設在時間 $k-1$ 的機率密度函數 $p(x_{k-1}\mid y_{1:k-1})$ 為已知。

預測

利用查普曼 - 科爾莫戈羅夫等式 (Chapman–Kolmogorov equation)，可以由狀態轉移函數與時間 $k-1$ 的機率密度函數 $p(x_{k-1}\mid y_{1:k-1})$ ，計算出時間 $k$ 的先驗機率 $p(x_{k}\mid y_{1:k-1})$

$$
{\begin{aligned}p(x_{k}\mid y_{1:k-1})&=\int p(x_{k},x_{k-1}\mid y_{1:k-1})d{x_{k-1}}\\&=\int p(x_{k}\mid x_{k-1},y_{1:k-1})p(x_{k-1}\mid y_{1:k-1})d{x_{k-1}}\\&=\int p(x_{k}\mid x_{k-1})p(x_{k-1}\mid y_{1:k-1})d{x_{k-1}}\\\end{aligned}}
$$

其中，由於狀態轉移模型被假設為一階馬可夫過程，時間 $k-1$ 的狀態只由時間 $k$ 決定，因此 $p(x_{k}\mid x_{k-1},y_{1:k-1})=p(x_{k}\mid x_{k-1})$ 。機率模型 $p(x_{k}\mid x_{k-1})$ 由狀態轉移函數 $x_{k}=f_{k}(x_{k-1},v_{k-1})$ 和 $v_{k-1}$ 的統計值得到。

更新

在時間 $k$ ，我們得到觀測值 $y_{k}$ ，因此可以利用貝氏定理，由先驗機率 $p(x_{k}\mid y_{1:k-1})$ 得到後驗機率 $p(x_{k}\mid y_{1:k})$ ，也就是考慮觀測值後得到的機率。

$$
p(x_{k}\mid y_{1:k})={p(y_{k}\mid x_{k})p(x_{k}\mid y_{1:k-1}) \over p(y_{k}\mid y_{1:k-1})}
$$

其中的歸一化常數為

$$
p(y_{k}\mid y_{1:k-1})=\int p(y_{k}\mid x_{k})p(x_{k}\mid y_{1:k-1})dx_{k}
$$

其中的概似函數 $p(y_{k}\mid x_{k})$ 由觀測函數 $y_{k}=h_{k}(x_{k},n_{k})$ 和 $n_{k}$ 的統計值得到。

上述 「預測」 與 「更新」 的遞迴關係，是貝葉斯最佳解的基本概念，然而公式中運用到的積分，對於一般非線性、非高斯的系統，難以得到解析解，因此需要利用蒙地卡羅方法來近似貝葉斯最佳解。



序列重要性重採樣 (Sequential Importance Resampling, SIR)

序列重要性採樣 (Sequential Importance Sampling, SIS)

SIR 是由 SIS 加上重採樣 (resampling) 改良而來，在 SIS 中，我們將後驗機率 $p(x_{k}\mid y_{1:k})$ 用 $N$ 個隨機採樣的樣本 (稱為粒子) 與各自的權重表示為

$$
{\{x_{k}^{(i)},w_{k}^{(i)}\}}_{i=1}^{N}
$$

其中的權重是根據重要性採樣的規則產生，且必須符合

$$
\sum _{i=1}^{N}w_{k}^{(i)}=1
$$

SIS 是將重要性採樣遞迴執行的一種方法，根據重要性採樣的理論，一個函數 $f$ 的期望值可以用加權平均來近似

$$
\int f(x_{k})p(x_{k}\mid y_{1:k})dx_{k}\approx \sum _{i=1}^{N}w_{k}^{(i)}f(x_{k}^{(i)})=1
$$

在每一次遞迴過程中，我們希望由前一次採樣的權重，計算出下一次採樣的權重。假設採樣的樣本分布可以表示為

$$
x^{(i)}\sim q(x),$ $i=1,...,N
$$

其中 $q(x)$ 稱為重要性密度 (importance density)。若樣本 $x_{k}^{(i)}$ 是由重要性密度 $q(x_{k}\mid y_{1:k})$ 抽取出來，則權重可表示為

$$
w_{k}^{(i)}\propto {p(x_{k}^{(i)}\mid y_{1:k}) \over q(x_{k}^{(i)}\mid y_{1:k})}
$$

我們將重要性密度分解為

$$
q(x_{k}\mid y_{1:k})=q(x_{k}\mid x_{k-1},y_{1:k})q(x_{k-1}\mid y_{1:k-1})
$$

再將後驗機率表示為

$$
{\begin{aligned}p(x_{k}\mid y_{1:k})&={p(y_{k}\mid x_{k},y_{1:k-1})p(x_{k}\mid y_{1:k-1}) \over p(y_{k}\mid y_{1:k-1})}\\&={p(y_{k}\mid x_{k},y_{1:k-1})p(x_{k}\mid x_{k-1},y_{1:k-1})p(x_{k-1}\mid y_{1:k-1}) \over p(y_{k}\mid y_{1:k-1})}\\&={p(y_{k}\mid x_{k})p(x_{k}\mid x_{k-1})p(x_{k-1}\mid y_{1:k-1}) \over p(y_{k}\mid y_{1:k-1})}\\&\propto p(y_{k}\mid x_{k})p(x_{k}\mid x_{k-1})p(x_{k-1}\mid y_{1:k-1})\\\end{aligned}}
$$

則權重的遞迴式可以表示為

$$
{\begin{aligned}w_{k}^{(i)}&\propto {p(y_{k}\mid x_{k}^{(i)})p(x_{k}^{(i)}\mid x_{k-1}^{(i)})p(x_{k-1}^{(i)}\mid y_{1:k-1}) \over q(x_{k}^{(i)}\mid x_{k-1}^{(i)},y_{1:k})q(x_{k-1}^{(i)}\mid y_{1:k-1})}\\&=w_{k-1}^{(i)}{p(y_{k}\mid x_{k}^{(i)})p(x_{k}^{(i)}\mid x_{k-1}^{(i)}) \over q(x_{k}^{(i)}\mid x_{k-1}^{(i)},y_{1:k})}\\\end{aligned}}
$$

重採樣（resampling）

SIS 可能會造成退化問題（degeneracy problem），也就是在經過幾次遞迴後，很多粒子的權重都變小到可以忽略不計，只剩少數粒子的權重較大，如此會浪費大量的計算量在幾乎沒有作用的粒子上，而使估計性能下降。由於在遞迴過程中，權重的變異數只會愈來愈大，因此退化問題無法被避免。

為了評估退化問題，我們定義有效粒子數為

$$
{\widehat {N_{eff}}}={1 \over \sum _{i=1}^{N}(w_{k}^{(i)})^{2}}
$$

在進行 SIS 時，若有效粒子數小於某一閾值，則對粒子做重採樣，即可減緩退化問題。重採樣的概念是去除權重過小的粒子，專注於權重較大的粒子。進行重採樣時，要由現有的粒子分布取樣，產生一組新的粒子 ${\{x_{k}^{(i)*}\}}_{i=1}^{N}$ ，由於產生出的新樣本為獨立且相同分布 (i.i.d.)，因此將權重重新設定為 $w_{k}^{(i)}={1 \over N}$ 。



參見


* 圖像處理

* 機器學習

* 機器人學

* 人工智慧

* 同時定位與地圖構建

* 滾動時域估計



參考文獻

1. M. S. Arulampalam, S. Maskell, N. Gordon and T. Clapp, "A tutorial on particle filters for online nonlinear/non-Gaussian Bayesian tracking." In IEEE Transactions on Signal Processing, vol. 50, no. 2, pp. 174-188, Feb 2002.
2. A. Doucet, N. de Freitas, N. Gordon, "An Introduction to Sequential Monte Carlo Methods." In A. Doucet, N. de Freitas, N. Gordon (eds.) Sequential Monte Carlo Methods in Practice. Statistics for Engineering and Information Science. Springer, New York, NY, 2001
3. A. Doucet, "On sequential Monte Carlo methods for Bayesian filtering." Dept. Eng., Univ. Cambridge, UK, Tech. Rep., 1998.