趙丹丹 劉靜娜 賀康建

摘要：高斯粒子濾波是一種免重采樣的粒子濾波，不會(huì)出現(xiàn)粒子退化，但其重要性密度函數(shù)由于沒(méi)有考慮到最新量測(cè)信息，使得濾波性能明顯下降，且該算法沒(méi)有較高的實(shí)時(shí)性。針對(duì)這個(gè)問(wèn)題提出一種基于CKF的高斯粒子濾波算法—CKGPF算法。該算法利用CKF算法構(gòu)造高斯粒子濾波的重要性密度函數(shù)，且在時(shí)間更新階段借助CKF算法來(lái)完成只對(duì)高斯分布參數(shù)的更新。仿真結(jié)果表明，CKGPF算法相比于標(biāo)準(zhǔn)GPF算法不僅提高了濾波精度，而且還具有較好的實(shí)時(shí)性。

關(guān)鍵詞：高斯粒子濾波；重要性密度函數(shù)；實(shí)時(shí)性；容積卡爾曼濾波

中圖分類(lèi)號(hào)：TP391 文獻(xiàn)標(biāo)識(shí)碼：ADOI：10.3969/j.issn.10036199.2017.01.017

粒子濾波[1]是一種基于蒙特卡洛思想的非線(xiàn)性、非高斯系統(tǒng)的濾波方法，已在定位、目標(biāo)跟蹤、無(wú)線(xiàn)通信、目標(biāo)識(shí)別等領(lǐng)域得到了深入研究[2-4]和廣泛應(yīng)用。但在標(biāo)準(zhǔn)的粒子濾波算法中，一般把先驗(yàn)概率密度函數(shù)作為重要性密度函數(shù)，這種方法沒(méi)有把最新的量測(cè)值考慮進(jìn)去，因而使得從重要性概率密度函數(shù)采樣得到的樣本與從真實(shí)后驗(yàn)概率密度函數(shù)采樣得到的樣本之間存在很大的偏差，特別是當(dāng)似然分布比較陡峭或是似然函數(shù)位于狀態(tài)轉(zhuǎn)移概率密度的尾部時(shí)，出現(xiàn)的偏差就更為顯著。因此，為了在算法中將最新的量測(cè)信息考慮進(jìn)去，通常是將不同卡爾曼濾波和粒子濾波結(jié)合來(lái)構(gòu)造重要性密度函數(shù)[5-7]，但文獻(xiàn)[8]已給出證明盲目將兩種算法結(jié)合來(lái)構(gòu)造新的重要性密度函數(shù)并不一定能獲得較好的濾波精度，且會(huì)增加計(jì)算復(fù)雜度，降低算法的實(shí)時(shí)性，而且粒子濾波本身就存在實(shí)時(shí)性較差的問(wèn)題。隨著有限集多目標(biāo)跟蹤的出現(xiàn)， 粒子濾波也進(jìn)入了隨機(jī)集粒子濾波[9]新的發(fā)展階段。文獻(xiàn)[10]更是實(shí)現(xiàn)了隨機(jī)集粒子濾波對(duì)多目標(biāo)跟蹤的應(yīng)用。高斯粒子濾波[11]是針對(duì)解決粒子濾波重采樣中樣本枯竭問(wèn)題的一種改進(jìn)的粒子濾波，具有更高的濾波精度。但GPF同樣也是把先驗(yàn)概率密度函數(shù)作為重要性密度函數(shù)，沒(méi)有考慮最新的量測(cè)信息，更新階段對(duì)所有的粒子進(jìn)行更新，使得高斯粒子濾波沒(méi)有較高的實(shí)時(shí)性。文獻(xiàn)[12]提出了一種基于GHF的高斯粒子濾波算法，該算法采用GHF構(gòu)造高斯粒子濾波的重要性密度函數(shù)但沒(méi)有考慮算法實(shí)時(shí)性問(wèn)題。本文提出選用濾波精度明顯高于EKF和UKF的CKF來(lái)構(gòu)造GPF的重要性密度函數(shù)，同時(shí)在GPF的時(shí)間更新階段通過(guò)CKF來(lái)完成只對(duì)高斯分布函數(shù)參數(shù)的更新，仿真結(jié)果表明，本文算法不僅提高了濾波精度，還具有更好的實(shí)時(shí)性。

2CKF算法

CKF算法是以非線(xiàn)性高斯濾波用以框架為基礎(chǔ)，利用Cubature規(guī)則選取2n個(gè)等權(quán)值的Cubature點(diǎn)，將這些Cubature點(diǎn)經(jīng)過(guò)非線(xiàn)性的系統(tǒng)方程進(jìn)行變換，然后進(jìn)行加權(quán)求和來(lái)近似系統(tǒng)狀態(tài)的后驗(yàn)均值和協(xié)方差。相比于EKF和UKF算法，CKF算法具有嚴(yán)密的理論推導(dǎo)，簡(jiǎn)單易于實(shí)現(xiàn)且具有較高的濾波精度。

CKF具體算法和步驟如下[13]：

（1）時(shí)間更新

6結(jié)論

本文提出的基于CKF的高斯粒子濾波算法，通過(guò)CKF構(gòu)造重要性密度函數(shù)，考慮了最新的量測(cè)信息，增加了有效粒子數(shù)，使得濾波精度明顯提高，同時(shí)本文算法在時(shí)間更新階段利用CKF來(lái)完成只對(duì)高斯分布參數(shù)的更新，減少了算法的計(jì)算量，使得算法具有更高的實(shí)時(shí)性。仿真結(jié)果表明，本文算法無(wú)論是濾波精度還是實(shí)時(shí)性，都有了一定的改善。

參考文獻(xiàn)

[1]GORDON N J，SALMOND D J，SMITH A F M.Novel approach to nonlinear/non-Gaussian Bayesian state estimation[J].IEE Proceedings F Radar and Signal Processing， 1993，140（2）：107-113.

[2]ZHOU S K，CHELLAPPA R，MOGHADDAM B.Visual tracking and recognition using appearance adaptive models in particle filters[J].IEEE Trans.on Image Processing，2004，13（11）：1492-1506.

[3]SABUNCE M，DEMIREKLER M.IMMPF for increased performance for maneuvering target[J]. IEEE on Signal Processing and Communications Applications Conf. Antalya：IEEE，2011：1044-1047.

[4]BANANI S A，VAUGHAN R G.Blind channel estimation and discrete speed tracking in wireless systems using independent component analysis with particle filtering[J].IET Communications，2012，6（2）：224-234.

[5]李良群，姬紅兵，羅軍輝.迭代擴(kuò)展卡爾曼粒子濾波器[J].西安電子科技大學(xué)學(xué)報(bào)：自然科學(xué)版，2007，34（2）：233-238.

[6]QI C，BONDON P.A new unscented particle filter[C].IEEE International Conference on Acoustics，Speech and Signal Processing，2008：3417-3420.

[7]袁澤劍，鄭南寧，賈新春.高斯厄米特粒子濾波器[J].電子學(xué)報(bào)，2003，31（7）：970-973.

[8]LI T C，CORCHADO J M，BAJO J，et al.Effectiveness of Bayesian filters： an information fusion perspective[J]. Information Sciences， 2016， 329： 670-689.

[9]李天成，范紅旗， 孫樹(shù)棟.粒子濾波理論、方法及其在多目標(biāo)跟蹤中的應(yīng)用[J].自動(dòng)化學(xué)報(bào)，2015，41（12）： 1981-2002.

[10]UNEY M，CLARK D E，JULIER S J.Distributed fusion of PHD filters via exponential mixture densities[J]. IEEE Journal of Selected Topics in Signal Processing， 2013， 7（3）： 521-531.

[11]KOTECHA J H，DJURIC P M.Gaussian particle filtering[J]. IEEE Transactions on Signal Processing， 2003， 51（10）：2592-2601.

[12]鄭作虎， 王首勇.一種基于GHF 的高斯粒子濾波算法[J].控制與決策，2014，29（9）：1698-1702.

[13]ARASARATNAM I，HAYKIN S.Cubature Kalman Filters[J]. IEEETransactions on Automatic Control ，2009，54（6）：1254-1269.