卡爾曼濾波器在車輛跟蹤預(yù)測的應(yīng)用(上)

狀態(tài)預(yù)測

卡爾曼濾波首先考慮的問題是目標(biāo)狀態(tài)指么,這也叫系統(tǒng)狀態(tài)仙逻。以小車在高速公路上行駛為例,小車的狀態(tài)可以用位置和速度來表示涧尿。用數(shù)學(xué)符號表示即x_{t} =\left[ \begin{matrix}   p_{t}\\   v_{t}\\     \end{matrix}  \right] ,其中p_{t}表示位置,v_{t}表示速度檬贰。

小車在行駛姑廉,既可以做勻速直線運動,也可以做加速運動翁涤。假定u表示小車的加速度桥言,那么小車進(jìn)行勻速運動時,u=0t時刻葵礼,小車的位移和速度可以分別表示為:

p_{t}= p_{t-1}+v_{t-1}\times \Delta t+\frac{\Delta t^2}{2 } \times u_{t}

v_{t}= v_{t-1}+\Delta t\times u_{t}

因為當(dāng)前時刻的位置和速度分別是上一時刻的線性組合号阿,因此可以用矩陣的形式來描述上述兩個公式, 即:\left[ \begin{matrix} p_{t}\\ v_{t}\\ \end{matrix} \right] =\left[ \begin{matrix} 1&\Delta t\\ 0&1\\ \end{matrix} \right] \left[ \begin{matrix} p_{t-1}\\ v_{t-1}\\ \end{matrix} \right] +\left[ \begin{matrix} \frac{\Delta t^2 }{2} \\ \Delta t\\ \end{matrix} \right] u_{t}

F_{t} =\left[ \begin{matrix} 1&\Delta t\\ 0&1\\ \end{matrix} \right] ,B_{t} =\left[ \begin{matrix} \frac{\Delta t^2 }{2} \\ \Delta t\\ \end{matrix} \right],那么上面公式就轉(zhuǎn)換為:

\hat x_{t}^-= F_{t}\hat x_{t}^-+B_{t}u_{t}

劃重點鸳粉,此公式也就是卡爾曼濾波器的第一個公式扔涧,叫狀態(tài)預(yù)測公式。F_{t} 就是系統(tǒng)狀態(tài)轉(zhuǎn)移矩陣届谈,表示上一個時刻狀態(tài)怎么變換到下一個時刻的狀態(tài)枯夜,B_{t}是控制矩陣,表示控制量u是怎么作用于當(dāng)前狀態(tài)艰山。\hat x_{t}^- 就是我們的預(yù)測值湖雹。細(xì)心的可以發(fā)現(xiàn),這里的x帶有帽子曙搬,沒有寫錯摔吏,這說明了這里的x是系統(tǒng)狀態(tài)的估計值,這個值是通過上一個時刻的狀態(tài)(卡爾曼濾波修正后的狀態(tài)值)計算估計出來的纵装。既然x是估計值征讲,那么就存在噪聲,我們用協(xié)方差矩陣表示當(dāng)前估計值和真實值相比搂擦,帶來多大的不確定性稳诚,那么就有卡爾曼濾波的第二個公式,也就是不確定性在各個時刻的傳遞公式:

P_{t}^- =FP_{t-1}F^{T} +Q

P表示協(xié)方差矩陣,F是狀態(tài)轉(zhuǎn)移矩陣,考慮在估計下一時刻的狀態(tài)x瀑踢,本身存在的噪聲扳还,這個噪聲用Q來表示才避。這樣下一時刻的協(xié)方差矩陣P的估計值,就得到了氨距。綜上桑逝,卡爾曼濾波器的狀態(tài)預(yù)測兩個公式就說明清楚了。

觀測矩陣


  小車在行駛俏让,那么這個觀察矩陣是什么東西楞遏。這個跟實際的需求有關(guān),可以是位置首昔,可以是速度寡喝,也可以是速度和位置等。在這個例子中以小車的位置為觀察量勒奇,用數(shù)學(xué)符號z_{t}表示预鬓。因為預(yù)測值x_{t}和觀測值z_{t}之間是線性關(guān)系(卡爾曼濾波器處理線性問題),那么轉(zhuǎn)換關(guān)系就可以用矩陣H表示赊颠,H叫觀測矩陣格二。需要考慮的是估計值x_{t} 和觀測值z_{t}的維度不總是相同,那么觀測矩陣H是根據(jù)實際情況構(gòu)造的竣蹦。

  小車行駛例子,估計值x_{t}.shape = (2,1)顶猜,z_{t}.shape=(1,1),那么H.shape =(1,2),H=\left[ \begin{matrix}  1&0 \end{matrix}  \right] ,據(jù)此,得到卡爾曼濾波的第三公式:

z_{t}=Hx_{t}+v痘括。

v表示測量噪聲长窄,因為測量值不是100%準(zhǔn)確,即實際的測量值跟你測出來的測量值還是存在誤差纲菌,那么用R表示這個誤差抄淑,專業(yè)術(shù)語就是觀測噪聲的協(xié)方差矩陣。小車?yán)佑^測值z_{t}是一個標(biāo)量驰后,那么R就只代表方差了肆资。

狀態(tài)更新

  有了小車狀態(tài)的估計值和測量值,就可以用卡爾曼濾波器來修正測量值了灶芝。

?????? \hat x_{t}=\hat x_{t}^-+K_{t}(z_t-H\hat x_{t}^-)

???? K_t =P_t^{-}H^T{(HP_t^{-}H^{T}+R)}^{-1}

???? 其中(z_t-H\hat x_{t}^-)是測量值和估計值之間的殘差郑原,K_{t}是卡爾曼系數(shù),用來衡量預(yù)測狀態(tài)協(xié)方差P和觀測狀態(tài)協(xié)方差R的大小進(jìn)而相信預(yù)測模型多一點還是觀察模型多一點夜涕。相信預(yù)測模型多一點犯犁,則殘差的權(quán)重小一點,反之女器,殘差的權(quán)重大一點酸役。

噪聲協(xié)方差矩陣更新:

 卡爾曼濾波器的最后一個公式就是對最佳估計值的噪聲分布的更新了:

 P_{t}= (1-K_{t}H)P_{t}^-

以上,卡爾曼濾波器的所有的公式和參數(shù)就解釋完啦,下集將會給出自己的實戰(zhàn)例子涣澡。

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末贱呐,一起剝皮案震驚了整個濱河市,隨后出現(xiàn)的幾起案子入桂,更是在濱河造成了極大的恐慌奄薇,老刑警劉巖,帶你破解...
    沈念sama閱讀 207,248評論 6 481
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件抗愁,死亡現(xiàn)場離奇詭異馁蒂,居然都是意外死亡,警方通過查閱死者的電腦和手機蜘腌,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,681評論 2 381
  • 文/潘曉璐 我一進(jìn)店門沫屡,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人撮珠,你說我怎么就攤上這事谁鳍。” “怎么了劫瞳?”我有些...
    開封第一講書人閱讀 153,443評論 0 344
  • 文/不壞的土叔 我叫張陵,是天一觀的道長绷柒。 經(jīng)常有香客問我志于,道長,這世上最難降的妖魔是什么废睦? 我笑而不...
    開封第一講書人閱讀 55,475評論 1 279
  • 正文 為了忘掉前任伺绽,我火速辦了婚禮,結(jié)果婚禮上嗜湃,老公的妹妹穿的比我還像新娘奈应。我一直安慰自己,他們只是感情好购披,可當(dāng)我...
    茶點故事閱讀 64,458評論 5 374
  • 文/花漫 我一把揭開白布杖挣。 她就那樣靜靜地躺著,像睡著了一般刚陡。 火紅的嫁衣襯著肌膚如雪惩妇。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 49,185評論 1 284
  • 那天筐乳,我揣著相機與錄音歌殃,去河邊找鬼。 笑死蝙云,一個胖子當(dāng)著我的面吹牛氓皱,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播,決...
    沈念sama閱讀 38,451評論 3 401
  • 文/蒼蘭香墨 我猛地睜開眼波材,長吁一口氣:“原來是場噩夢啊……” “哼股淡!你這毒婦竟也來了?” 一聲冷哼從身側(cè)響起各聘,我...
    開封第一講書人閱讀 37,112評論 0 261
  • 序言:老撾萬榮一對情侶失蹤揣非,失蹤者是張志新(化名)和其女友劉穎,沒想到半個月后躲因,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體早敬,經(jīng)...
    沈念sama閱讀 43,609評論 1 300
  • 正文 獨居荒郊野嶺守林人離奇死亡,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 36,083評論 2 325
  • 正文 我和宋清朗相戀三年大脉,在試婚紗的時候發(fā)現(xiàn)自己被綠了搞监。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點故事閱讀 38,163評論 1 334
  • 序言:一個原本活蹦亂跳的男人離奇死亡镰矿,死狀恐怖琐驴,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情秤标,我是刑警寧澤绝淡,帶...
    沈念sama閱讀 33,803評論 4 323
  • 正文 年R本政府宣布,位于F島的核電站苍姜,受9級特大地震影響牢酵,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜衙猪,卻給世界環(huán)境...
    茶點故事閱讀 39,357評論 3 307
  • 文/蒙蒙 一馍乙、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧垫释,春花似錦丝格、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,357評論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至订咸,卻和暖如春琅束,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背算谈。 一陣腳步聲響...
    開封第一講書人閱讀 31,590評論 1 261
  • 我被黑心中介騙來泰國打工涩禀, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留,地道東北人然眼。 一個月前我還...
    沈念sama閱讀 45,636評論 2 355
  • 正文 我出身青樓艾船,卻偏偏與公主長得像,于是被迫代替她去往敵國和親。 傳聞我的和親對象是個殘疾皇子屿岂,可洞房花燭夜當(dāng)晚...
    茶點故事閱讀 42,925評論 2 344

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

  • 卡爾曼濾波在我當(dāng)學(xué)生的時候就用過践宴,但是當(dāng)年我似乎就是套公式,沒有理解其精髓爷怀,加之時間久了有點模糊阻肩,突然需要指導(dǎo)學(xué)生...
    Roger_羅杰閱讀 83,490評論 41 159
  • 經(jīng)過看各種博客和文章,讓我最清楚明白的运授,是xiahouzuoxin 的博客烤惊,之后又看了一些國外的文獻(xiàn)進(jìn)行自己的理解...
    marine0131閱讀 7,395評論 4 12
  • 預(yù)估器 我們希望可以最大限度地使用測量結(jié)果來估計移動物體的運動。所以吁朦,多個測量的累積可以讓我們檢測出不受噪聲影響的...
    JasonDing閱讀 6,744評論 4 10
  • 姓名:周小蓬 16019110037 轉(zhuǎn)載自:http://blog.csdn.net/MangZuo/artic...
    aeytifiw閱讀 3,156評論 1 13
  • 第五章四十五言 【文】 傅奕《道德經(jīng)古本篇》第五章 天地不仁柒室,以萬物為芻狗。圣人不仁逗宜,以百姓為芻狗雄右。天地之間其猶橐...
    孤鴻學(xué)士閱讀 816評論 5 7