亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        多通道信號自校正分布式融合Kalman濾波器

        2014-03-05 18:29:11劉文強(qiáng)韓娜顧澤元陶貴麗
        現(xiàn)代電子技術(shù) 2014年3期

        劉文強(qiáng)+韓娜++顧澤元++陶貴麗

        摘 要: 對帶未知參數(shù)的多傳感器多通道自回歸滑動平均(ARMA)信號,采用多維遞推輔助變量(MRIV)方法得到自回歸模型參數(shù)估值,通過Gevers?Wouters算法辨識滑動平均模型參數(shù)估值,再用相關(guān)方法得到噪聲方差的估值。把所有的估值都代入到最優(yōu)分布式融合信息濾波器中得到自校正分布式融合Kalman信息濾波器。該濾波器具有漸近全局最優(yōu)性,一個多通道信號仿真例子驗(yàn)證了其有效性。

        關(guān)鍵詞: 多通道ARMA信號; 多段辨識方法; 多重遞推輔助變量法; 信息濾波器

        中圖分類號: TN911.7?34; O211.64 文獻(xiàn)標(biāo)識碼: A 文章編號: 1004?373X(2014)03?0061?04

        Kalman filter of self?tuning distributed fusion for multi?channel signal

        LIU Wen?qiang, HAN Na, GU Ze?yuan, TAO Gui?li

        (Heilongjiang University of Science and Technology, Harbin 150022, China)

        Abstract: For multisensor and multi?channel autoregressive moving average (ARMA) signal with unknown parameters, its autoregressive model parameter estimators are obtained by multi?dimensional recursive instrumental variable algorithm. The moving average model parameter estimators are identified through Gevers?Wouters algorithm, and the noise variances estimators are obtained with the correlation method. All of the estimators are substitute into optimal distributed fusion information filter to obtain self?tuning distributed fusion Kalman information filter, which has asymptotic global optimality. The effectiveness is verified through a multi?channel signal simulation example.

        Keywords: multi?channel ARMA signal; multi?stage identification method; multidimensional recursive auxiliary variable method; information filter

        0 引 言

        近年來,廣泛應(yīng)用于很多高科技領(lǐng)域的多傳感器信息融合Kalman濾波理論已經(jīng)引起特別的關(guān)注,它主要研究信息融合狀態(tài)估計(jì)問題或信號估計(jì)問題。最優(yōu)信息融合Kalman濾波包括集中式融合和分布式融合方法。其中集中式融合是全局最優(yōu)的,即集中式融合是把所有傳感器觀測方程合為一個高維數(shù)的觀測方程,但是此時系統(tǒng)的觀測方程維數(shù)有所增加,容錯性能差,計(jì)算負(fù)擔(dān)大。分布式融合有兩種方法:信息矩陣[1]和加權(quán)方差方法[2],其中加權(quán)方差分布式融合方法是全局次優(yōu)的,而基于信息矩陣的分布式融合同集中式融合,也能得到全局最優(yōu)的狀態(tài)估計(jì)。對于最優(yōu)信息融合Kalman信號濾波問題以往文獻(xiàn)有很多報(bào)道,但是,系統(tǒng)的模型參數(shù)和噪聲方差通常是完全或部分未知的。這種系統(tǒng)的濾波稱為自校正濾波[3]。文獻(xiàn)[4]針對單通道自回歸(AR)信號提出了自校正集中式濾波器,但未解決多通道信號辨識問題。文獻(xiàn)[5?6]對于多通道信號提出了自校正信息融合濾波器,但他們都不能解決自校正分布式信息濾波器問題。

        為了克服以往文獻(xiàn)的局限性,本文針對多傳感器多通道ARMA信號,當(dāng)模型參數(shù)和噪聲統(tǒng)計(jì)全部未知時,應(yīng)用多段信息融合辨識方法,其中利用MRIV方法求得AR模型參數(shù),用Gevers?Wouters算法辨識滑動平均(MA)模型參數(shù),用相關(guān)方法得到噪聲方差融合估值器,把未知參數(shù)帶到最優(yōu)分布式融合Kalman信息濾波器得到多通道信號自校正分布式Kalman信息濾波器,該濾波器是漸近全局最優(yōu)的。

        1 全局最優(yōu)的最優(yōu)分布式信息融合Kalman

        信號濾波器

        考慮多傳感器多通道ARMA信號:

        [A(q-1)s(t)=C(q-1)w(t)] (1)

        [yi(t)=s(t)+vi(t), i=1,2,…,L] (2)

        其中信號[s(t)∈Rm,]第[i]個傳感器的觀測[yi(t)∈Rm,][vi(t)∈][Rm]是觀測噪聲,[w(t)∈Rm]為輸入噪聲,[q-1]為單位滯后算子,[q-1s(t)=s(t-1)。][A(q-1)]與[C(q-1)]有形式[X(q-1)=X0+][X1q-1+…+Xnxq-nx,]其中[X0=Im,][Im]是[m×m]單位陣,[na]和[nc]分別是多項(xiàng)式[A(q-1)]和[C(q-1)]的階次,且[na=nc]。假設(shè)[(A(q-1),C(q-1))]左素,初始時刻[t0=-∞]。問題是當(dāng)[Ai,][Ci,][Qw]和[Qvi]均未知時,求漸近全局最優(yōu)自校正分布式融合信號濾波器[ss(t|t)]。

        假設(shè)1 [w(t)∈Rm]和[vi(t)∈Rm]分別是零均值、方差陣為[Qw]和[Qvi]的獨(dú)立白噪聲。

        假設(shè)2 各傳感器子系統(tǒng)是一致完全可觀和一致完全可控的。

        令[w(t)=w(t-1),]則[w(t)]有方差[Qw。]ARMA信號式(1)能夠改寫為[A(q-1)s(t)=C(q-1)w(t),]其中:[C(q-1)=C0+][C1q-1+…+Cncq-nc, ][C0=0,Ci=Ci-1,i=1,2,…,nc,nc=nc+1。]

        系統(tǒng)(1)和(2)有狀態(tài)空間模型:

        [x(t+1)=Φx(t)+Γw(t)] (3)

        [yi(t)=Hx(t)+vi(t), i=1,2,…,L] (4)

        [s(t)=Hx(t)] (5)

        [Φ=-A1?Imna-Ana+10…0Γ=C1C2?Cnc+1H=Im0…0] (6)

        其中[Aj=0(j>na),Cj=0(j>nc)]。集中式融合觀測方程為[y(t)=Η0x(t)+v(t),]定義[y(t)=yT1(t)…yTL(t)T,H0=][HT…HTT,][v(t)=vT1(t)…vTL(t)T,]融合觀測噪聲[v(t)]有方差[Qv=diagQv1…QvL。]

        引理1[7] 多傳感器系統(tǒng)式(3)和式(4)有基于信息矩陣的局部Kalman濾波器:

        [xi(t|t)=Ψi(t)xi(t-1|t-1)+Ki(t)yi(t)] (7)

        [Ki(t)=Pi(t|t)HTQ-1vi] (8)

        [Ψi(t)=Pi(t|t)Σ-1i(t|t-1)Φ] (9)

        [P-1i(t|t)=Σ-1i(t|t-1)+HTQ-1viH] (10)

        [Σi(t+1|t)=ΦPi(t|t)ΦΤ+ΓQwΓΤ] (11)

        引理2[7] 多傳感器系統(tǒng)式(3)和(4)在假設(shè)1和2下,全局最優(yōu)分布式融合Kalman信號濾波器為:

        [x(t|t)=Ψ(t)x(t-1|t-1)+P(t|t)i=1L[P-1i(t|t)xi(t|t)-Σ-1i(t|t-1)Φxi(t-1|t-1)]] (12)

        [P-1(t|t)=Σ-1(t|t-1)+i=1L[P-1i(t|t)-Σ-1i(t|t-1)]] (13)

        [Σ(t+1|t)=ΦP(t|t)ΦΤ+ΓQwΓΤ] (14)

        [Ψ(t)=P(t|t)Σ-1(t|t-1)Φ] (15)

        [s(t|t)=Hx(t|t)] (16)

        2 未知參數(shù)的信息融合辨識

        參數(shù)[Ai(i=1,2,…,na),][Ci(i=1,2,…,nc),][Qw]和[Qvi]都是未知的,信息融合多段辨識算法可以分為三段:

        第1段,用多維遞推輔助變量(MRIV)算法求自回歸(AR)模型參數(shù)的局部和融合估計(jì)

        由式和式(2)可得:

        [A(q-1)yi(t)=C(q-1)w(t)+A(q-1)vi(t)] (17)

        則式(17)可轉(zhuǎn)化為AR模型:

        [A(q-1)yi(t)=ri(t)] (18)

        [ri(t)=C(q-1)w(t)+A(q-1)vi(t)] (19)

        式(18)有最小二乘結(jié)構(gòu):

        [yi(t)=Θφi(t)+ri(t)] (20)

        [φTi(t)=[-yTi(t-1),…,-yTi(t-na)]],[Θ=[A1,…,Ana]] (21)

        引理3[8] 對于多傳感器多通道平穩(wěn)AR模型(20),[Θ]的多維遞推輔助變量(RIV)局部估值器[Θi(t)]為:

        [Θi(t)=Θi(t-1)+[yi(t)-Θi(t-1)φi(t)]φΤi(t)Pi(t-1)1+φΤi(t)Pi(t-1)φi(t)] (22)

        [Pi(t)=Pi(t-1)-Pi(t-1)φi(t)φTi(t)Pi(t-1)1+φTi(t)Pi(t-1)φi(t)] (23)

        [φi(t)=φi(t-na)] (24)

        且?guī)С踔礫Θi(t)=Θi0,][Pi(t0)=P0(k0)。][Θ]的信息融合估值定義為[Θf(t)=1Li=1LΘi(t)=A1(t),…,Ana(t)。]

        第2段,用Gevers?Wouters算法辨識[C(q-1)]和[Qw]

        定義滑動平均(MA)過程:

        [m(t)=C(q-1)w(t)] (25)

        它的相關(guān)函數(shù)[Rm(k)=E[m(t)mT(t-k)]]為[Rm(k)=α=kncCαQwCTα-k,k=0,1,2,…,nc。]計(jì)算式(25)兩邊的相關(guān)函數(shù)有[Rm(k)=Rrij(k),]由式(18)在時刻[t,]觀測過程[ri(t)]的估值定義為:

        [ri(t)=A(q-1)yi(t)] (26)

        觀測[ri(t)]的采樣相關(guān)函數(shù)估值定義為:

        [Rtrij(k)=1tu=1tri(u)rTj(u-k)] (27)

        它有遞推形式:

        [Rtrij(k)=Rt-1rij(k)+1tri(t)rTj(t-k)-Rt-1rij(k)] (28)

        初值為[R1rij(k)=1tj=1tri(1)rTj(1-k)]。把估值[ri(t)]代入到式(28)得到估值[Rtrij(k)],應(yīng)用帶死區(qū)的G?W算法[9],能夠得到局部的估值[Cuij(t)][(u=1,2,…,nc;i,j=1,2,…,L)]和[Qwij(t)。]基于所有傳感器的信息融合估值為:

        [Qw(t)=1L2i=1Lj=1LQwij(t),Cu(t)=1L2i=1Lj=1LCuij(t)] (29)

        第3段,基于相關(guān)方法辨識[Qvi]

        觀測[ri(t)]的自相關(guān)函數(shù)[Rrii(k)=E[ri(t)rTi(t-k)]]為:

        [Rrii(k)=α=kncCαQwCTα-k+β=knaAβQviATβ-k, k=0,1,2,…,nr,nr=max(Rrii(k)≠0)] (30)

        對于不同的[k],把采樣估值[Rtrii(k)],融合估值[Ai(t),][Ci(t)]和[Qw(t)]代入到式(30)解線性方程組就可以得到噪聲方差估值[Qvi(t)]。

        3 基于信息矩陣的自校正分布式融合Kalman

        信號濾波器

        將辨識得到的估值[Ai(t)]和[Ci(t)]代入到式(6)中得到[Φ(t)]和[Γ(t)],融合觀測噪聲[v(t)]有噪聲方差估值為[Qv(t)=diagQv1(t)…QvL(t)],把所有的估值代到引理1可得到自校正分布式融合Kalman信息濾波器:

        [xs(t|t)=Ψ(t)xs(t-1|t-1)+P(t|t)i=1L[P-1i(t|t)xsi(t|t)-Σ-1i(t|t-1)Φ(t)xsi(t-1|t-1)]] (31)

        [P-1(t|t)=Σ-1(t|t-1)+i=1L[P-1i(t|t)-Σ-1i(t|t-1)]] (32)

        [Σ(t+1|t)=Φ(t)P(t|t)ΦΤ(t)+Γ(t)Qw(t)ΓΤ(t)] (33)

        [Ψ(t)=P(t|t)Σ-1(t|t-1)Φ(t)] (34)

        [ss(t|t)=Hxs(t|t)] (35)

        4 仿真例子

        考慮多傳感器多通道ARMA信號(1)和(2),[i=1,2,3。]其中[A(q-1)=I2+A1q-1,][C(q-1)=I2+C1q-1,][w(t)]和[vi(t)]是零均值,方差各為[Qw]和[Qvi]的相互獨(dú)立白噪聲,假設(shè)[A1,C1,Qw]和[Qvi]未知,求自校正分布式融合Kalman信號濾波器。在仿真中?。?/p>

        [A1=-0.80-0.3-0.7,][C1=0.150-0.30.3,] [Qw=0.81001,][Qv1=0.1000.2,][Qv2=0.3000.49,][Qv3=0.4000.5]。

        基于信息融合多段辨識方法,得到[A1,C1,Qw]和[Qvi]的信息融合估值,辨識結(jié)果如圖1~圖6所示,其中真實(shí)值用直線表示,融合估計(jì)用實(shí)曲線表示,[D(k,r)]表示矩陣[D]的第[(k,r)]個元素。圖1所示為應(yīng)用引理1得到[A1]的融合估值曲線,用Gevers?Wouters算法得到[C1]和[Qw]的融合估值如圖2和圖3所示,應(yīng)用相關(guān)方法得到觀測噪聲方差[Qvi]的融合估值曲線如圖4~圖6所示。圖7是誤差曲線,說明最優(yōu)分布式融合Kalman信號濾波器[s(t|t)=s1(t|t)s2(t|t)T]和自校正分布式融合Kalman信號濾波器[ss(t|t)=ss1(t|t)ss2(t|t)T]之間的誤差趨于零。

        圖1 融合估值[A1(t)]的曲線

        圖2 融合估值[C1(t)]的曲線

        圖3 噪聲方差[Qw]融合估值曲線

        圖4 噪聲方差[Qv1]融合估值曲線

        圖5 噪聲方差[Qv2]融合估值曲線

        圖6 噪聲方差[Qv3]融合估值曲線

        圖7 最優(yōu)和自校正分布式融合Kalman信號濾波器的誤差曲線

        5 結(jié) 語

        本文對于多通道ARMA信號模型,提出了自校正分布式融合信息濾波器。該自校正分布式融合信息濾波器不需要求解Riccati方程,避免[mL×mL]矩陣[(H0Σ(t|t-1)HT0+Qv(t))]的求逆運(yùn)算,其中[m=m1+…+mL,]而僅需要計(jì)算[mi×mi]矩陣的逆,計(jì)算量較小,且該自校正分布式信息濾波器是漸近全局最優(yōu)的。

        參考文獻(xiàn)

        [1] CHANG K C, ZHI T, SAHA R K. Performance evaluation of track fusion with information matrix filter [J]. Aerospace and Electronic Systems, 2002, 38(2): 455?464.

        [2] DENG Z L, GAO Y, MAO L, et al. New approach to information fusion steady?state Kalman filtering [J]. Automatica, 2005, 41(10): 1695?1707.

        [3] MOIR T, GRIMBLE M J. Optimal self?tuning filtering, prediction, and Smoothing for multivariable processes [J]. IEEE Transactions on Automatic Control, 1984, 29(2): 128?137.

        [4] 劉文強(qiáng),陶貴麗,鄧自立.帶未知模型參數(shù)的自校正集中式融合信息濾波器[J].科學(xué)技術(shù)與工程,2010,10(33):9?14.

        [5] SUN X J, DENG Z L. Optimal and self?tuning weighted measurement fusion Wiener filter for the multisensor multichannel ARMA signals [J]. Signal Processing, 2009, 89(5): 738?752.

        [6] RAN C J, DENG Z L. Self?tuning distributed measurement fusion Kalman estimator for the multi?channel ARMA signal [J]. Signal Processing, 2011, 91(8): 2028?2041.

        [7] 鄧自立.信息融合濾波理論及其應(yīng)用[M].哈爾濱:哈爾濱工業(yè)大學(xué)出版社,2007.

        [8] 徐慧勤,鄧自立,張明波.多維和多重遞推輔助變量辨識算法[J].科學(xué)技術(shù)與工程,2010,10(2):366?370.

        [9] GEVERS M, WOUTERS W R E. An innovation approach to discrete?time stochastic realization problem [J]. Quarterly Journal of Automatic Control, 1978, 19(2): 90?110.

        [10] 趙大建,趙偉,張兆亮,等.基于FPGA的Kalman濾波器實(shí)現(xiàn)研究[J].現(xiàn)代電子技術(shù),2012,35(6):67?70.

        [Qw(t)=1L2i=1Lj=1LQwij(t),Cu(t)=1L2i=1Lj=1LCuij(t)] (29)

        第3段,基于相關(guān)方法辨識[Qvi]

        觀測[ri(t)]的自相關(guān)函數(shù)[Rrii(k)=E[ri(t)rTi(t-k)]]為:

        [Rrii(k)=α=kncCαQwCTα-k+β=knaAβQviATβ-k, k=0,1,2,…,nr,nr=max(Rrii(k)≠0)] (30)

        對于不同的[k],把采樣估值[Rtrii(k)],融合估值[Ai(t),][Ci(t)]和[Qw(t)]代入到式(30)解線性方程組就可以得到噪聲方差估值[Qvi(t)]。

        3 基于信息矩陣的自校正分布式融合Kalman

        信號濾波器

        將辨識得到的估值[Ai(t)]和[Ci(t)]代入到式(6)中得到[Φ(t)]和[Γ(t)],融合觀測噪聲[v(t)]有噪聲方差估值為[Qv(t)=diagQv1(t)…QvL(t)],把所有的估值代到引理1可得到自校正分布式融合Kalman信息濾波器:

        [xs(t|t)=Ψ(t)xs(t-1|t-1)+P(t|t)i=1L[P-1i(t|t)xsi(t|t)-Σ-1i(t|t-1)Φ(t)xsi(t-1|t-1)]] (31)

        [P-1(t|t)=Σ-1(t|t-1)+i=1L[P-1i(t|t)-Σ-1i(t|t-1)]] (32)

        [Σ(t+1|t)=Φ(t)P(t|t)ΦΤ(t)+Γ(t)Qw(t)ΓΤ(t)] (33)

        [Ψ(t)=P(t|t)Σ-1(t|t-1)Φ(t)] (34)

        [ss(t|t)=Hxs(t|t)] (35)

        4 仿真例子

        考慮多傳感器多通道ARMA信號(1)和(2),[i=1,2,3。]其中[A(q-1)=I2+A1q-1,][C(q-1)=I2+C1q-1,][w(t)]和[vi(t)]是零均值,方差各為[Qw]和[Qvi]的相互獨(dú)立白噪聲,假設(shè)[A1,C1,Qw]和[Qvi]未知,求自校正分布式融合Kalman信號濾波器。在仿真中?。?/p>

        [A1=-0.80-0.3-0.7,][C1=0.150-0.30.3,] [Qw=0.81001,][Qv1=0.1000.2,][Qv2=0.3000.49,][Qv3=0.4000.5]。

        基于信息融合多段辨識方法,得到[A1,C1,Qw]和[Qvi]的信息融合估值,辨識結(jié)果如圖1~圖6所示,其中真實(shí)值用直線表示,融合估計(jì)用實(shí)曲線表示,[D(k,r)]表示矩陣[D]的第[(k,r)]個元素。圖1所示為應(yīng)用引理1得到[A1]的融合估值曲線,用Gevers?Wouters算法得到[C1]和[Qw]的融合估值如圖2和圖3所示,應(yīng)用相關(guān)方法得到觀測噪聲方差[Qvi]的融合估值曲線如圖4~圖6所示。圖7是誤差曲線,說明最優(yōu)分布式融合Kalman信號濾波器[s(t|t)=s1(t|t)s2(t|t)T]和自校正分布式融合Kalman信號濾波器[ss(t|t)=ss1(t|t)ss2(t|t)T]之間的誤差趨于零。

        圖1 融合估值[A1(t)]的曲線

        圖2 融合估值[C1(t)]的曲線

        圖3 噪聲方差[Qw]融合估值曲線

        圖4 噪聲方差[Qv1]融合估值曲線

        圖5 噪聲方差[Qv2]融合估值曲線

        圖6 噪聲方差[Qv3]融合估值曲線

        圖7 最優(yōu)和自校正分布式融合Kalman信號濾波器的誤差曲線

        5 結(jié) 語

        本文對于多通道ARMA信號模型,提出了自校正分布式融合信息濾波器。該自校正分布式融合信息濾波器不需要求解Riccati方程,避免[mL×mL]矩陣[(H0Σ(t|t-1)HT0+Qv(t))]的求逆運(yùn)算,其中[m=m1+…+mL,]而僅需要計(jì)算[mi×mi]矩陣的逆,計(jì)算量較小,且該自校正分布式信息濾波器是漸近全局最優(yōu)的。

        參考文獻(xiàn)

        [1] CHANG K C, ZHI T, SAHA R K. Performance evaluation of track fusion with information matrix filter [J]. Aerospace and Electronic Systems, 2002, 38(2): 455?464.

        [2] DENG Z L, GAO Y, MAO L, et al. New approach to information fusion steady?state Kalman filtering [J]. Automatica, 2005, 41(10): 1695?1707.

        [3] MOIR T, GRIMBLE M J. Optimal self?tuning filtering, prediction, and Smoothing for multivariable processes [J]. IEEE Transactions on Automatic Control, 1984, 29(2): 128?137.

        [4] 劉文強(qiáng),陶貴麗,鄧自立.帶未知模型參數(shù)的自校正集中式融合信息濾波器[J].科學(xué)技術(shù)與工程,2010,10(33):9?14.

        [5] SUN X J, DENG Z L. Optimal and self?tuning weighted measurement fusion Wiener filter for the multisensor multichannel ARMA signals [J]. Signal Processing, 2009, 89(5): 738?752.

        [6] RAN C J, DENG Z L. Self?tuning distributed measurement fusion Kalman estimator for the multi?channel ARMA signal [J]. Signal Processing, 2011, 91(8): 2028?2041.

        [7] 鄧自立.信息融合濾波理論及其應(yīng)用[M].哈爾濱:哈爾濱工業(yè)大學(xué)出版社,2007.

        [8] 徐慧勤,鄧自立,張明波.多維和多重遞推輔助變量辨識算法[J].科學(xué)技術(shù)與工程,2010,10(2):366?370.

        [9] GEVERS M, WOUTERS W R E. An innovation approach to discrete?time stochastic realization problem [J]. Quarterly Journal of Automatic Control, 1978, 19(2): 90?110.

        [10] 趙大建,趙偉,張兆亮,等.基于FPGA的Kalman濾波器實(shí)現(xiàn)研究[J].現(xiàn)代電子技術(shù),2012,35(6):67?70.

        [Qw(t)=1L2i=1Lj=1LQwij(t),Cu(t)=1L2i=1Lj=1LCuij(t)] (29)

        第3段,基于相關(guān)方法辨識[Qvi]

        觀測[ri(t)]的自相關(guān)函數(shù)[Rrii(k)=E[ri(t)rTi(t-k)]]為:

        [Rrii(k)=α=kncCαQwCTα-k+β=knaAβQviATβ-k, k=0,1,2,…,nr,nr=max(Rrii(k)≠0)] (30)

        對于不同的[k],把采樣估值[Rtrii(k)],融合估值[Ai(t),][Ci(t)]和[Qw(t)]代入到式(30)解線性方程組就可以得到噪聲方差估值[Qvi(t)]。

        3 基于信息矩陣的自校正分布式融合Kalman

        信號濾波器

        將辨識得到的估值[Ai(t)]和[Ci(t)]代入到式(6)中得到[Φ(t)]和[Γ(t)],融合觀測噪聲[v(t)]有噪聲方差估值為[Qv(t)=diagQv1(t)…QvL(t)],把所有的估值代到引理1可得到自校正分布式融合Kalman信息濾波器:

        [xs(t|t)=Ψ(t)xs(t-1|t-1)+P(t|t)i=1L[P-1i(t|t)xsi(t|t)-Σ-1i(t|t-1)Φ(t)xsi(t-1|t-1)]] (31)

        [P-1(t|t)=Σ-1(t|t-1)+i=1L[P-1i(t|t)-Σ-1i(t|t-1)]] (32)

        [Σ(t+1|t)=Φ(t)P(t|t)ΦΤ(t)+Γ(t)Qw(t)ΓΤ(t)] (33)

        [Ψ(t)=P(t|t)Σ-1(t|t-1)Φ(t)] (34)

        [ss(t|t)=Hxs(t|t)] (35)

        4 仿真例子

        考慮多傳感器多通道ARMA信號(1)和(2),[i=1,2,3。]其中[A(q-1)=I2+A1q-1,][C(q-1)=I2+C1q-1,][w(t)]和[vi(t)]是零均值,方差各為[Qw]和[Qvi]的相互獨(dú)立白噪聲,假設(shè)[A1,C1,Qw]和[Qvi]未知,求自校正分布式融合Kalman信號濾波器。在仿真中取:

        [A1=-0.80-0.3-0.7,][C1=0.150-0.30.3,] [Qw=0.81001,][Qv1=0.1000.2,][Qv2=0.3000.49,][Qv3=0.4000.5]。

        基于信息融合多段辨識方法,得到[A1,C1,Qw]和[Qvi]的信息融合估值,辨識結(jié)果如圖1~圖6所示,其中真實(shí)值用直線表示,融合估計(jì)用實(shí)曲線表示,[D(k,r)]表示矩陣[D]的第[(k,r)]個元素。圖1所示為應(yīng)用引理1得到[A1]的融合估值曲線,用Gevers?Wouters算法得到[C1]和[Qw]的融合估值如圖2和圖3所示,應(yīng)用相關(guān)方法得到觀測噪聲方差[Qvi]的融合估值曲線如圖4~圖6所示。圖7是誤差曲線,說明最優(yōu)分布式融合Kalman信號濾波器[s(t|t)=s1(t|t)s2(t|t)T]和自校正分布式融合Kalman信號濾波器[ss(t|t)=ss1(t|t)ss2(t|t)T]之間的誤差趨于零。

        圖1 融合估值[A1(t)]的曲線

        圖2 融合估值[C1(t)]的曲線

        圖3 噪聲方差[Qw]融合估值曲線

        圖4 噪聲方差[Qv1]融合估值曲線

        圖5 噪聲方差[Qv2]融合估值曲線

        圖6 噪聲方差[Qv3]融合估值曲線

        圖7 最優(yōu)和自校正分布式融合Kalman信號濾波器的誤差曲線

        5 結(jié) 語

        本文對于多通道ARMA信號模型,提出了自校正分布式融合信息濾波器。該自校正分布式融合信息濾波器不需要求解Riccati方程,避免[mL×mL]矩陣[(H0Σ(t|t-1)HT0+Qv(t))]的求逆運(yùn)算,其中[m=m1+…+mL,]而僅需要計(jì)算[mi×mi]矩陣的逆,計(jì)算量較小,且該自校正分布式信息濾波器是漸近全局最優(yōu)的。

        參考文獻(xiàn)

        [1] CHANG K C, ZHI T, SAHA R K. Performance evaluation of track fusion with information matrix filter [J]. Aerospace and Electronic Systems, 2002, 38(2): 455?464.

        [2] DENG Z L, GAO Y, MAO L, et al. New approach to information fusion steady?state Kalman filtering [J]. Automatica, 2005, 41(10): 1695?1707.

        [3] MOIR T, GRIMBLE M J. Optimal self?tuning filtering, prediction, and Smoothing for multivariable processes [J]. IEEE Transactions on Automatic Control, 1984, 29(2): 128?137.

        [4] 劉文強(qiáng),陶貴麗,鄧自立.帶未知模型參數(shù)的自校正集中式融合信息濾波器[J].科學(xué)技術(shù)與工程,2010,10(33):9?14.

        [5] SUN X J, DENG Z L. Optimal and self?tuning weighted measurement fusion Wiener filter for the multisensor multichannel ARMA signals [J]. Signal Processing, 2009, 89(5): 738?752.

        [6] RAN C J, DENG Z L. Self?tuning distributed measurement fusion Kalman estimator for the multi?channel ARMA signal [J]. Signal Processing, 2011, 91(8): 2028?2041.

        [7] 鄧自立.信息融合濾波理論及其應(yīng)用[M].哈爾濱:哈爾濱工業(yè)大學(xué)出版社,2007.

        [8] 徐慧勤,鄧自立,張明波.多維和多重遞推輔助變量辨識算法[J].科學(xué)技術(shù)與工程,2010,10(2):366?370.

        [9] GEVERS M, WOUTERS W R E. An innovation approach to discrete?time stochastic realization problem [J]. Quarterly Journal of Automatic Control, 1978, 19(2): 90?110.

        [10] 趙大建,趙偉,張兆亮,等.基于FPGA的Kalman濾波器實(shí)現(xiàn)研究[J].現(xiàn)代電子技術(shù),2012,35(6):67?70.

        久久久久九九精品影院 | 日本大尺度吃奶呻吟视频| 无码精品一区二区免费AV| av永远在线免费观看| 男女啪啪视频高清视频| 国产99在线 | 亚洲| 国产精品视频牛仔裤一区| 宅男久久精品国产亚洲av麻豆| 日本少妇一区二区三区四区| 一区二区三区人妻少妇| 欧美成人片在线观看| 欧美日韩性视频| 日韩女同一区在线观看| 国产精品女同一区二区免费站| 97色伦综合在线欧美视频| 国产美熟女乱又伦av果冻传媒| 久久亚洲av午夜福利精品西区| 亚洲高清一区二区精品| 精品国产综合区久久久久久 | 久久青青草原国产毛片| 日本道精品一区二区三区| 欧美 亚洲 国产 日韩 综AⅤ| 素人系列免费在线观看| 日本一区二区三区光视频| 中国老熟女重囗味hdxx| 欧美精品在线一区| 久久精品国产白丝爆白浆| 精品高朝久久久久9999| 日本边添边摸边做边爱的网站| 亚洲欧美成人中文在线网站| 国产精品av免费网站| 精品人妻少妇av中文字幕| 成人区人妻精品一熟女| 精品91精品91精品国产片| 国产一区二区黑丝美胸| 久久精品丝袜高跟鞋| 亚州少妇无套内射激情视频| 国产亚洲精选美女久久久久 | 色综合中文综合网| 人人妻人人澡av| 99精品人妻少妇一区二区三区 |