轉(zhuǎn)移矩陣未知時(shí)帶乘性噪聲系統(tǒng)的自適應(yīng)濾波算法

褚東升,劉 祺,張 玲

(中國(guó)海洋大學(xué)工程學(xué)院,山東青島266100)

轉(zhuǎn)移矩陣未知時(shí)帶乘性噪聲系統(tǒng)的自適應(yīng)濾波算法

褚東升,劉 祺,張 玲

(中國(guó)海洋大學(xué)工程學(xué)院,山東青島266100)

針對(duì)帶乘性噪聲系統(tǒng)狀態(tài)轉(zhuǎn)移矩陣未知的情況,提出1種在線性最小方差意義下的系統(tǒng)參數(shù)和狀態(tài)聯(lián)合濾波算法。以迭代方式獲得模型參數(shù)和系統(tǒng)狀態(tài)的遞推估計(jì):首先,利用之前時(shí)刻的狀態(tài)估計(jì)值,根據(jù)投影定理,對(duì)系統(tǒng)未知參數(shù)即系統(tǒng)狀態(tài)轉(zhuǎn)移矩陣作出估計(jì);其次,利用已得到的系統(tǒng)參數(shù)估計(jì)值,獲取當(dāng)前時(shí)刻的狀態(tài)濾波。計(jì)算機(jī)仿真結(jié)果表明了算法的有效性。

乘性噪聲;轉(zhuǎn)移矩陣;自適應(yīng)濾波;線性最小方差估計(jì)

自Kalman濾波方法提出以來(lái),在信號(hào)估計(jì)領(lǐng)域取得了眾多的理論成果,并在空間技術(shù)、通訊、導(dǎo)航等許多領(lǐng)域得到了廣泛應(yīng)用[1-3]。但在使用Kalman濾波解決實(shí)際問(wèn)題時(shí)的1個(gè)很大限制是必須預(yù)先知道模型的參數(shù)或其統(tǒng)計(jì)性質(zhì),如系統(tǒng)初始值、各模型參數(shù)矩陣、噪聲的統(tǒng)計(jì)特性等。而實(shí)際中這些參數(shù)往往是未知的,或近似已知的,或部分已知的:這常常導(dǎo)致大的狀態(tài)估計(jì)誤差,甚至?xí)篂V波器發(fā)散[1]。因此,對(duì)自適應(yīng)Kalman濾波的研究十分必要。另外,在許多實(shí)際觀測(cè)過(guò)程中,如石油地震勘探、通訊工程、語(yǔ)音處理等,往往需要考慮系統(tǒng)的時(shí)變性、非線性畸變、能量衰減等各種復(fù)雜甚至不確定因素。這些因素在數(shù)學(xué)上可以近似歸結(jié)為乘性噪聲[4]。此類(lèi)估計(jì)問(wèn)題的難點(diǎn)在于,由于乘性噪聲的引入,系統(tǒng)量測(cè)方程不再是線性的。近年來(lái),這種帶乘性噪聲系統(tǒng)的濾波問(wèn)題已經(jīng)引起人們的極大關(guān)注[4-8]。

目前,針對(duì)不帶乘性噪聲隨機(jī)線性系統(tǒng)的自適應(yīng)濾波算法的理論與應(yīng)用研究已經(jīng)取得了較豐富的成果[9-14],其中,文獻(xiàn)[11]討論了當(dāng)系統(tǒng)模型矩陣、噪聲方差矩陣和初值向量中含有未知參向量時(shí)的線性估計(jì)。文獻(xiàn)[12]討論了狀態(tài)空間規(guī)范形式下轉(zhuǎn)移矩陣與控制矩陣含有未知參數(shù)時(shí)的參數(shù)估計(jì)和聯(lián)合濾波。文獻(xiàn)[13]給出了一種將未知參數(shù)作為擴(kuò)展?fàn)顟B(tài),借助非線性濾波對(duì)狀態(tài)空間模型進(jìn)行線性化處理,并利用擴(kuò)展Kalman濾波對(duì)模型轉(zhuǎn)移矩陣進(jìn)行估計(jì)的方法。但該方法只適用于轉(zhuǎn)移矩陣為單常數(shù)的情況,且是在不帶乘性噪聲的條件下進(jìn)行的。文獻(xiàn)[14]給出了控制矩陣和初值含未知參向量時(shí)情況下的估計(jì)公式,并給出了嚴(yán)格的證明。

在帶乘性噪聲系統(tǒng)自適應(yīng)濾波的研究上,文獻(xiàn)[6]提出了1種參數(shù)和狀態(tài)預(yù)報(bào)估計(jì)的兩段互耦式算法,并以石油勘探領(lǐng)域中的地震道模型為例,對(duì)其參數(shù)辨識(shí)方法進(jìn)行了詳細(xì)推導(dǎo)。文獻(xiàn)[17]對(duì)帶乘性噪聲系統(tǒng)在乘性噪聲及其統(tǒng)計(jì)參數(shù)未知的情況下,基于極大似然準(zhǔn)則,提出了1種分塊組合優(yōu)化估計(jì)算法。

本文將針對(duì)帶乘性噪聲系統(tǒng),研究狀態(tài)轉(zhuǎn)移矩陣的遞推估計(jì)方法及其系統(tǒng)狀態(tài)濾波。

1 問(wèn)題描述

對(duì)于帶乘性噪聲的離散隨機(jī)系統(tǒng),其狀態(tài)空間表達(dá)的一般形式是:

其中x(k)為n維狀態(tài)向量,z(k)為q維觀測(cè)向量, w(k)為h維動(dòng)態(tài)噪聲,v(k)為q維觀測(cè)噪聲。A∈Rn×n為未知狀態(tài)轉(zhuǎn)移矩陣,B∈Rn×n、C∈Rq×n分別為相應(yīng)維數(shù)的常值矩陣;m(k)為一維乘性噪聲。對(duì)任意整數(shù)k≥0,j≥0,滿(mǎn)足以下假設(shè)條件:

假設(shè)1 E{w(k)}=0,E{w(k)wT(j)}=Q(k)δkj;

假設(shè)2 E{v(k)}=0,E{v(k)vT(j)}=R(k)δkj;

假設(shè)3 E{m(k)}=￣m,E{[m(k)-￣m][m(j)-￣m]T}= σ(k)δkj;

假設(shè)4 E{x(0)}=0,E{x(0)xT(0)}=S(0);

假設(shè)5 {m(k)}、{w(k)}、{v(k)}及{x(0)}相互統(tǒng)計(jì)獨(dú)立。

2 一類(lèi)線性模型的線性最小方差估計(jì)算法

古典最小二乘法是最為常用的參數(shù)辨識(shí)方法,且在一定條件下,加權(quán)最小二乘估計(jì)即等價(jià)于線性最小方差估計(jì)[13]。下面利用這種等價(jià)關(guān)系,以引理的形式給出一類(lèi)線性模型的線性最小方差估計(jì)算法。

引理1[13]對(duì)于如下線性模型:

其中Z(k)為量測(cè)向量,G(k)為相應(yīng)維數(shù)的已知矩陣, J(k)為未知參向量,ε(k)為隨機(jī)噪聲向量且ε(k)相互獨(dú)立即E{ε(k)}=0,E{ε(k)εT(j)}=V(k)δkj。

有J(k)的線性最小方差無(wú)偏估計(jì)為:

本文將在以上引理的基礎(chǔ)上,對(duì)系統(tǒng)(1)和(2)中的未知轉(zhuǎn)移矩陣A進(jìn)行估計(jì),并利用所得結(jié)果計(jì)算系統(tǒng)狀態(tài)濾波值。

另外,在估計(jì)系統(tǒng)轉(zhuǎn)移矩陣時(shí),需要用到上一時(shí)刻系統(tǒng)狀態(tài)濾波值,因此,在這里以引理的形式不加證明的給出帶乘性噪聲系統(tǒng)的濾波遞推式。

引理2[4]滿(mǎn)足假設(shè)條件1～5的系統(tǒng)(1)和(2),有如下在線性最小方差意義下的狀態(tài)濾波遞推算法:

狀態(tài)濾波估計(jì)

^x(k)=^x(k/k-1)+K(k)(z(k)-￣mC^x(k/k-1))

增益矩陣

K(k)=￣m(k)P(k/k-1)CTR-1L(k)

新息方差陣

RL(k)=R(k)+￣m2CP(k/k-1)CT+σ(k)CS(k)CT

狀態(tài)相關(guān)矩陣

S(k)=A(k/k-)S(k-1)AT(k/k-1)+BQ(k-1)BT

濾波誤差方差陣為

P(k)=P(k/k-1)-￣mK(k)CP(k/k-1)

狀態(tài)一步預(yù)測(cè)

^x(k/k-1)=A(k/k-1)^x(k/k-1)

一步預(yù)測(cè)誤差方差矩陣

P(k/k-1)=A(k/k-1)P(k-1)AT(k/k-1)+BQ(k-1)BT

初始值^x(1/0)=0,P(1/0)=S(0)。

3 帶乘性噪聲系統(tǒng)狀態(tài)轉(zhuǎn)移矩陣未知時(shí)的自適應(yīng)濾波算法

下面利用之前給出的2個(gè)引理,給出線性最小方差意義下,狀態(tài)轉(zhuǎn)移矩陣的遞推估計(jì)算法。定理 考慮系統(tǒng)(1)、(2),若系統(tǒng)滿(mǎn)足假設(shè)條件1～5,則有系統(tǒng)狀態(tài)轉(zhuǎn)移矩陣在線性最小方差意義下的如下遞推估計(jì)算法。

證明 首先,對(duì)(1)式進(jìn)行改寫(xiě):

其中,Vec(A)是對(duì)矩陣A按列向量依次排成的向量,即矩陣向量化運(yùn)算[15]。

將(8)式代入(2)式得

其中,“?”即為矩陣的Kronecker乘積運(yùn)算[15]。

在(9)式中,令H(k)=mT(k)x(k-1)=m(k)x(k-1)

根據(jù)投影定理[16]有

其中^x(k-1)由引理2的狀態(tài)濾波遞推式得到。

將(10)式代入(9)式得

為證明方便,令

其中,^L(k)已知,Vec(A)為未知參向量,ξ(k)為噪聲向量。

將(12)式與(3)式類(lèi)比,由(4)式可得未知參向量Vec(A)的線性最小方差無(wú)偏估計(jì)如下:

為方便利用計(jì)算機(jī)計(jì)算,將(13)式化作像Kalman濾波那樣逐步遞推的形式[11]

取k=1,2,…,n,合并寫(xiě)作

此即(5)式。

用(20)式替換(17)式中相關(guān)項(xiàng),即得(7)式。

至此,定理證畢。

將由以上定理得出的轉(zhuǎn)移矩陣估計(jì)值代入引理2,即可對(duì)(1)和(2)式所描述的帶乘性噪聲系統(tǒng)進(jìn)行狀態(tài)濾波估計(jì),得到的濾波值用于下一時(shí)刻對(duì)狀態(tài)轉(zhuǎn)移矩陣的估計(jì)運(yùn)算,即兩者互為條件。如此遞推下去,即得到系統(tǒng)狀態(tài)濾波與狀態(tài)轉(zhuǎn)移矩陣的迭代濾波算法。

4 數(shù)值仿真

針對(duì)(1)和(2)式所描述的帶乘性噪聲系統(tǒng)模型,作者對(duì)本文提出的轉(zhuǎn)移矩陣未知情況下的狀態(tài)濾波算法進(jìn)行了大量的仿真研究,并且下面以狀態(tài)向量為二維的系統(tǒng)為例,給出仿真結(jié)果。各參數(shù)設(shè)置如下:

圖1 二維狀態(tài)轉(zhuǎn)移矩陣各元素估計(jì)值Fig.1 The estimate of the elements of the transfer matrix

圖2 狀態(tài)第一個(gè)分量的濾波Fig.2 The filtering of the first state component

圖3 狀態(tài)第二個(gè)分量的濾波Fig.3 The filtering of the second state component

在每一時(shí)刻,基于上一時(shí)刻的濾波結(jié)果,估計(jì)當(dāng)前時(shí)刻轉(zhuǎn)移矩陣,并利用所得估計(jì)值,對(duì)系統(tǒng)當(dāng)前時(shí)刻狀態(tài)進(jìn)行濾波估計(jì)。仿真效果如圖1～3所示。圖1中,對(duì)轉(zhuǎn)移矩陣的估計(jì)很快趨于穩(wěn)定,為便于更清晰的進(jìn)行觀察,在這里只取前一半時(shí)刻的估計(jì)結(jié)果。

5 結(jié)語(yǔ)

本文針對(duì)1種轉(zhuǎn)移矩陣未知情況下的帶乘性噪聲系統(tǒng),在其狀態(tài)濾波估計(jì)的基礎(chǔ)上,給出了轉(zhuǎn)移矩陣在線性最小方差意義下的估計(jì)算法,并將算法化作了便于計(jì)算機(jī)實(shí)現(xiàn)的遞推形式,實(shí)現(xiàn)了系統(tǒng)狀態(tài)濾波與狀態(tài)轉(zhuǎn)移矩陣之間的迭代估計(jì),從而形成了帶乘性噪聲系統(tǒng)在狀態(tài)轉(zhuǎn)移矩陣未知情況下的1種自適應(yīng)濾波算法。仿真結(jié)果表明了算法的有效性。

[1] Lewis F L,Xie L,Popa D.Optimal and robust estimation:with an introduction to stochastic control theory[M].NY:CRC Press,2008.

[2] Fan Wang,V Balakrishnan.Robust Kalman filters for linear timevarying systems with stochastic parametric uncertainties[J]. IEEE Trans on Signal Processing,2002,50(4):803-813.

[3] 付夢(mèng)印,鄧志紅,閻莉萍.卡爾曼濾波理論及其在導(dǎo)航系統(tǒng)中的應(yīng)用[M].北京:科學(xué)出版社,2010.

[4] Rajasekaran P K,Satyanarayana N,Srinath M D.Optimum linear estimation of stochastic signals in the presence of multiplicative noise[J].IEEE Trans Aerosp Electron Syst,1971,AES-7 462-468.

[5] Chow B S,Birkemeier W P.A recursive estimator that uses data in blocks with application to time invariant systems with multiplicative noise[C].USA:Proc IEEE Int Symp Circuits and Syst, 1987:98-101.

[6] 褚東升.帶乘性噪聲系統(tǒng)的估計(jì)理論與應(yīng)用[D].哈爾濱:哈爾濱工業(yè)大學(xué),1994.

[7] Primbs J A.Stochastic receding horizon control of constrained linear systems with state and control multiplicative noise[C].New York:2007 American Control Conference,2007:4470-4475.

[8] 褚東升,張征.帶乘性噪聲的一類(lèi)非線性系統(tǒng)的濾波算法[J].中國(guó)海洋大學(xué)學(xué)報(bào):自然科學(xué)版,2006,36(4):569-572.

[9] Shiuh-Ku Weng,Chung-Ming Kuo,Shu-Kang Tu.Video object tracking using adaptive Kalman filter[J].Journal of Visual Communication and Image Representation,2006(6):1190-1208.

[10] Dimitris G,Manolakis,Vinay K,et al.Statistical and Adaptive Signal Processing[M].Beijing:Tsinghua University Press, 2003.

[11] Nahi N E,Knobbe E J.Optimal liner recursive estimation with uncertain system parameters[J].IEEE Transactions on Automatic Control,1976(4):263-266.

[12] 丁峰,謝新民.系統(tǒng)參數(shù)和狀態(tài)聯(lián)合估計(jì)[J].控制與決策, 1994,9(3):223-225.

[13] 中國(guó)科學(xué)院數(shù)學(xué)研究所概率組.離散時(shí)間系統(tǒng)濾波的數(shù)學(xué)方法[M].北京:國(guó)防工業(yè)出版社,1975.

[14] 汪咬元.量測(cè)方程中合有參向量的遞推濾波[J].系統(tǒng)工程學(xué)報(bào),1995,72-84.

[15] 王松桂.線性模型的理論及其應(yīng)用[M].合肥:安徽教育出版社,1999.

[16] 鄧自立.最優(yōu)估計(jì)理論及其應(yīng)用[M].哈爾濱:工業(yè)大學(xué)出版社,2005.

[17] 張文林,褚東升,侯永海.帶乘性噪聲系統(tǒng)的極大似然最優(yōu)估計(jì)算法[J].青島海洋大學(xué)學(xué)報(bào),1999,29(2):332-338.

Abstract: In the sense of linear minimum-variance,this paper develops a corresponding recursive algorithm for stochastic linear systems with multiplicative noise whose transfer matrix is unknown.That is to say,estimations of parameters and states are obtained by the way of iteration.First,according to projection theorem,the system parameters are estimated by the using of the state estimation of previous moment.Then,based on the estimation of the system parameters,the optimal filtering of the current moment is obtained.Simulations are illustrated to show the validity of the algorithms.

Key words: multiplicative noise;transfer matrix;adaptive filtering;linear minimum-variance estimation

責(zé)任編輯 陳呈超

Adaptive Estimation Algorithm for Systems with Multiplicative Noise and Unknown Transfer Matrix

CHU Dong-Sheng,LIU Qi,ZHANG Ling
(College of Engineering,Ocean University of China,Qingdao 266100,China)

TP31

A

1672-5174(2011)06-115-05

國(guó)家自然科學(xué)基金項(xiàng)目(60704023)資助

2010-10-26;

2010-12-16

褚東升(1956-),男,教授。E-mail:chuds@yahoo.cn