卡爾曼濾波(2)數(shù)據(jù)融合,協(xié)方差矩陣,狀態(tài)空間

1數(shù)據(jù)融合
將多個(gè)不準(zhǔn)確的數(shù)值轉(zhuǎn)化成相對(duì)準(zhǔn)確的數(shù)值(自己的理解)癣猾。
同樣舉個(gè)例子:假設(shè)測(cè)量一個(gè)物體的質(zhì)量敛劝,測(cè)了兩次的質(zhì)量分別為Z1=30g,Z2=32g纷宇,兩組數(shù)據(jù)均有偏差夸盟,假設(shè)偏差分別為σ1=2g,σ2=4g像捶。滿(mǎn)足高斯正態(tài)分布上陕。
如何用倆次測(cè)量的不準(zhǔn)確的值來(lái)估計(jì)出物體真實(shí)的值,這便是數(shù)據(jù)融合拓春。

由上一張的公式
859D6902785396681800024919ADD59E.png

得到^Z= Z1 +K(Z2 - Z1),其中k就是卡爾曼增益(0<=K<=1)释簿。我們?cè)賮?lái)分析一下,當(dāng)K=0時(shí)硼莽,那么^Z = Z1,說(shuō)明信任第一次測(cè)量的值庶溶,同理,當(dāng)K=1時(shí)沉删,^Z = Z2渐尿,說(shuō)明信任第二次測(cè)量的值醉途。那么我們?nèi)绾吻笠粋€(gè)最接近真實(shí)值的^Z矾瑰? 這里就需要求出一個(gè)K使得方差最小,方差最小也就是標(biāo)準(zhǔn)差最邪妗(方差=標(biāo)準(zhǔn)差的平方)殴穴。我們用var(^Z)來(lái)表示方差。
IMG_3321.JPG

極值處货葬,方差最小采幌。我們便可以得到下面的式子。
IMG_3322.JPG

因此震桶,^Z = Z1 + K(Z2 - Z1) = 30.4g,
這便是數(shù)據(jù)融合的一個(gè)簡(jiǎn)單例子休傍。通過(guò)兩個(gè)不準(zhǔn)確的數(shù)值得到一個(gè)相對(duì)準(zhǔn)確的數(shù)值。

2協(xié)方差矩陣
協(xié)方差:介紹協(xié)方差矩陣之前蹲姐,要先對(duì)協(xié)方差要有一定的理解磨取。
首先來(lái)看它的定義式:Cov[X,Y]=E[(X?μ)(Y?ν)],其中μ=E(X),ν=E(Y)柴墩。而E()帶表了數(shù)學(xué)期望忙厌。而X,Y為隨機(jī)變量。這個(gè)式子是為了表示隨機(jī)變量與它的期望相差多少江咳。由此我們可以得到結(jié)論:
當(dāng)協(xié)方差為正值時(shí)逢净,兩個(gè)隨機(jī)變量正相關(guān);當(dāng)協(xié)方差為負(fù)值時(shí),兩個(gè)隨機(jī)變量負(fù)相關(guān)爹土,當(dāng)協(xié)方差為0的話(huà)甥雕,那么兩個(gè)隨機(jī)變量就不具有相關(guān)性。也就是說(shuō)胀茵,協(xié)防差可以表示變量之間的聯(lián)動(dòng)關(guān)系犀农。
協(xié)方差矩陣:將方差和協(xié)方差在同一個(gè)矩陣中表示出來(lái),從而判斷變量之間的相關(guān)性宰掉。
那么為什么要引入?yún)f(xié)方差矩陣呢呵哨?這里需要做一個(gè)解釋。
協(xié)方差只能處理二維問(wèn)題轨奄,對(duì)于三維以上數(shù)據(jù)孟害,就需要計(jì)算多個(gè)協(xié)方差,然后用矩陣將其組織起來(lái)挪拟,這就是協(xié)方差矩陣挨务。
下面列出一個(gè)n=3時(shí)的協(xié)防差矩陣,一看便知道怎么列出協(xié)方差矩陣玉组。


IMG_3323.JPG

狀態(tài)空間:用彈簧阻尼系統(tǒng)舉例谎柄。

IMG_3326.JPG

設(shè)x'為x的一階導(dǎo),x''為x的二階導(dǎo)惯雳,x為位移朝巫,x的一階導(dǎo)也就是速度,x的二階導(dǎo)也就是加速度石景,B為阻尼器的阻尼系數(shù)劈猿,阻尼系數(shù)與速度相乘就代表了阻尼器產(chǎn)生的力。由此我們可以列出表達(dá)式潮孽。
mx'' + Bx' + kx = F (下面均將F寫(xiě)成u揪荣,后面寫(xiě)成狀態(tài)空間表達(dá)形式方便觀察)
mx'' = u - Bx' - kx 我們?cè)O(shè)x1 = x,x2 = x'往史。那么x2' = x'' = 1/mu - B/mx2 - k/mx1.

定義z1為速度仗颈,z2為加速度(z1,z2為傳感器測(cè)量的數(shù)據(jù)椎例,x1挨决,x2為計(jì)算的數(shù)據(jù))z1=x=x1(位置)z2=x=x2(速度)


IMG_3327.JPG

上面的第一個(gè)式子就是狀態(tài)空間的常見(jiàn)表達(dá)形式。后面兩個(gè)式子便是將等式轉(zhuǎn)化成狀態(tài)空間的表達(dá)形式粟矿。而第三個(gè)矩陣也就是輸出矩陣凰棉,也叫觀測(cè)矩陣。

如果要用拓展卡爾曼陌粹,就要線性系統(tǒng)離散化撒犀。就是要找到x的k+1時(shí)刻的狀態(tài)量和x的k時(shí)刻狀態(tài)量的相對(duì)關(guān)系,以及z的k時(shí)刻的狀態(tài)關(guān)系和x的k時(shí)刻的狀態(tài)關(guān)系。(t)代表t時(shí)刻或舞。
x'(t)= Ax(t)+ Bu(t) z(t)= Cx(t)
接下來(lái)講下如何離散化荆姆。(下標(biāo)k為時(shí)刻)
X(k) = AX(k-1)+ BU(k-1)+ W(k-1)
Z(k) = H*X(k)+ V(k)
W表示X(t)和不確定性,也就是X(k)的不確定性映凳,X(k)也就是卡爾曼濾波的計(jì)算誤差胆筒。而當(dāng)Z的值(也就是觀測(cè)到的值)和狀態(tài)變量是一個(gè)值,H矩陣就可以等于C矩陣诈豌。W(k-1)也叫做過(guò)程噪音仆救,V(k)也叫做測(cè)量噪音。

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末矫渔,一起剝皮案震驚了整個(gè)濱河市彤蔽,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌庙洼,老刑警劉巖顿痪,帶你破解...
    沈念sama閱讀 217,826評(píng)論 6 506
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場(chǎng)離奇詭異油够,居然都是意外死亡蚁袭,警方通過(guò)查閱死者的電腦和手機(jī),發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,968評(píng)論 3 395
  • 文/潘曉璐 我一進(jìn)店門(mén)石咬,熙熙樓的掌柜王于貴愁眉苦臉地迎上來(lái)揩悄,“玉大人,你說(shuō)我怎么就攤上這事碌补÷彩” “怎么了棉饶?”我有些...
    開(kāi)封第一講書(shū)人閱讀 164,234評(píng)論 0 354
  • 文/不壞的土叔 我叫張陵厦章,是天一觀的道長(zhǎng)。 經(jīng)常有香客問(wèn)我照藻,道長(zhǎng)袜啃,這世上最難降的妖魔是什么? 我笑而不...
    開(kāi)封第一講書(shū)人閱讀 58,562評(píng)論 1 293
  • 正文 為了忘掉前任幸缕,我火速辦了婚禮群发,結(jié)果婚禮上,老公的妹妹穿的比我還像新娘发乔。我一直安慰自己熟妓,他們只是感情好,可當(dāng)我...
    茶點(diǎn)故事閱讀 67,611評(píng)論 6 392
  • 文/花漫 我一把揭開(kāi)白布栏尚。 她就那樣靜靜地躺著起愈,像睡著了一般。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上抬虽,一...
    開(kāi)封第一講書(shū)人閱讀 51,482評(píng)論 1 302
  • 那天官觅,我揣著相機(jī)與錄音,去河邊找鬼阐污。 笑死休涤,一個(gè)胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的笛辟。 我是一名探鬼主播功氨,決...
    沈念sama閱讀 40,271評(píng)論 3 418
  • 文/蒼蘭香墨 我猛地睜開(kāi)眼,長(zhǎng)吁一口氣:“原來(lái)是場(chǎng)噩夢(mèng)啊……” “哼手幢!你這毒婦竟也來(lái)了疑故?” 一聲冷哼從身側(cè)響起,我...
    開(kāi)封第一講書(shū)人閱讀 39,166評(píng)論 0 276
  • 序言:老撾萬(wàn)榮一對(duì)情侶失蹤弯菊,失蹤者是張志新(化名)和其女友劉穎纵势,沒(méi)想到半個(gè)月后,有當(dāng)?shù)厝嗽跇?shù)林里發(fā)現(xiàn)了一具尸體管钳,經(jīng)...
    沈念sama閱讀 45,608評(píng)論 1 314
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡钦铁,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 37,814評(píng)論 3 336
  • 正文 我和宋清朗相戀三年,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了才漆。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片牛曹。...
    茶點(diǎn)故事閱讀 39,926評(píng)論 1 348
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡,死狀恐怖醇滥,靈堂內(nèi)的尸體忽然破棺而出黎比,到底是詐尸還是另有隱情,我是刑警寧澤鸳玩,帶...
    沈念sama閱讀 35,644評(píng)論 5 346
  • 正文 年R本政府宣布阅虫,位于F島的核電站,受9級(jí)特大地震影響不跟,放射性物質(zhì)發(fā)生泄漏颓帝。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,249評(píng)論 3 329
  • 文/蒙蒙 一窝革、第九天 我趴在偏房一處隱蔽的房頂上張望购城。 院中可真熱鬧,春花似錦虐译、人聲如沸瘪板。這莊子的主人今日做“春日...
    開(kāi)封第一講書(shū)人閱讀 31,866評(píng)論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)侮攀。三九已至史侣,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間魏身,已是汗流浹背惊橱。 一陣腳步聲響...
    開(kāi)封第一講書(shū)人閱讀 32,991評(píng)論 1 269
  • 我被黑心中介騙來(lái)泰國(guó)打工, 沒(méi)想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留箭昵,地道東北人税朴。 一個(gè)月前我還...
    沈念sama閱讀 48,063評(píng)論 3 370
  • 正文 我出身青樓,卻偏偏與公主長(zhǎng)得像家制,于是被迫代替她去往敵國(guó)和親正林。 傳聞我的和親對(duì)象是個(gè)殘疾皇子,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 44,871評(píng)論 2 354

推薦閱讀更多精彩內(nèi)容