濾波問題

在随机过程理論中的濾波問題(Filtering problem)是指針對信号处理及相關領域中,許多狀態估測問題的數學模型。大致概念是從不完整的、可能包括雜訊的觀測值中,建立有關系統真實值的「最佳估測」。最佳非線性濾波問題(甚至也包括非平稳过程問題)由(1959年、1960年)找到解答,在的研究及的研究中也有提到,Zakai建立了濾波器在條件機率未歸一情況下的簡化動態模型,稱為。不過一般情形下的解是無限維的。

目前已針對一些近似以及一些特定條件有深入的研究。例如在高斯隨機變數的假設下,最佳解是線性濾波器,也稱為维纳滤波及卡尔曼滤波。更一般的情形下,其解為無限維度,為了在有限記憶體的電腦中計算,需要進行有限維度的近似,有限維的近似型非線性濾波器比較會以启发為基礎,例如擴展型卡爾曼濾波器或是假定密度濾波器(Assumed Density Filters),也有更方法論導向的作法,例如Projection Filters,其中有些子系列恰好和假定密度濾波器相同。

一般來說,若可以適用分離原理,這些濾波器也可以成為最优控制問題解的一部份。例如在LQG控制最佳控制問題中,其估測部份的解就是卡爾曼濾波。

數學表示
考慮概率空間 (Ω, Σ, P),並且假設在n維度欧几里得空间 Rn的系統,其在時間t的(隨機)狀態Yt為随机变量 Yt : Ω → Rn,可以由以下形式伊藤清隨機微分方程的解來求得

:\mathrm{d} Y_{t} = b(t, Y_{t}) \, \mathrm{d} t + \sigma (t, Y_{t}) \, \mathrm{d} B_{t},

其中B是標準p維布朗运动,b : [0, +∞) × Rn → Rn為漂移場(drift field),且σ : [0, +∞) × Rn → Rn×p是擴散場(diffusion field)。假設Rm內在每一個時間的觀測Ht(其中mn可能不同)由下式決定

:H_{t} = c(t, Y_{t}) + \gamma (t, Y_{t}) \cdot \mbox{noise}.

配合隨機微分方程的伊藤表示法,令

: Z_{t} = \int_{0}^{t} H_{s} \, \mathrm{d} s,

因此可以得到有關觀測Zt的隨機積分表示式:

:\mathrm{d} Z_{t} = c(t, Y_{t}) \, \mathrm{d} t + \gamma (t, Y_{t}) \, \mathrm{d} W_{t},

其中W表示標準r維的布朗运动,和B和初始條件Y0無關,c : [0, +∞) × Rn → Rn,且 γ : [0, +∞) × Rn → Rn×r

可以在所有tx,以及特定常數C的情形下,使下式成立:

:\big| c (t, x) \big| + \big| \gamma (t, x) \big| \leq C \big( 1 + | x | \big)

濾波問題如下:給定在0 ≤ s ≤ t時間內的觀測量Zs for 0 ≤ s ≤ t,依上述觀測值,針對系統真實狀態Yt的最佳估測Ŷt是什麼?

因為「依上述觀測量為基礎」,表示Ŷt是根據Zs觀測量中Σ-代数下的可測函數。令K = K(Zt) 是所有數值為Rn,平方可積分,而且Gt可量測隨機函數Y的集合:

:K = K(Z, t) = L^{2} (\Omega, G_{t}, \mathbf{P}; \mathbf{R}^{n}).

因為要求是「最佳估測」,表示Ŷt會讓YtK集合內所有候選估測值之間的均方差有最小值:

:\mathbf{E} \left[ \big| Y_{t} - \hat{Y}_{t} \big|^{2} \right] = \inf_{Y \in K} \mathbf{E} \left[ \big| Y_{t} - Y \big|^{2} \right]. \qquad \mbox{(M)}

基本結論:正交投影
候選估測值的空間K(Zt)是希尔伯特空间,根據希尔伯特空间的理論,可以推得最小值問題(M)的解Ŷt可以表示為下式

:\hat{Y}_{t} = P_{K(Z, t)} \big( Y_{t} \big),

其中PK(Z,t)表示將L2(Ω, Σ, PRn)映射到线性子空间 K(Zt) = L2(Ω, GtPRn)的正交投影。而且,有關其条件期望,可知道若F是Σ中的次σ代數,則正交投影
:P_{K} : L^{2} (\Omega, \Sigma, \mathbf{P}; \mathbf{R}^{n}) \to L^{2} (\Omega, F, \mathbf{P}; \mathbf{R}^{n})

也就是條件期望運算子E[·|F],也就是說

:P_{K} (X) = \mathbf{E} \big[ X \big | F \big].

因此

:\hat{Y}_{t} = P_{K(Z, t)} \big( Y_{t} \big) = \mathbf{E} \big[ Y_{t} \big | G_{t} \big].

這個基本結果是濾波理論中,廣義Fujisaki-Kallianpur-Kunita方程的基礎。

相關條目

  • 和濾波問題有緊密關係。
  • 濾波器
  • 信號處理中的
  • 卡尔曼滤波是濾波問題及平滑問題中最著名的解
  • 平滑

參考資料
*

  • (See Section 6.1)

评论 (0)

  • 还没有评论,来抢沙发吧。