原理
卡爾曼濾波器是一種基礎(chǔ)預(yù)測定位算法。原理非常簡單易懂。核心過程可以用一個(gè)圖說明:
![](https://file1.elecfans.com//web2/M00/94/A5/wKgaomTl9puAGdi7AABSGQezAxI620.jpg)
本質(zhì)上就是這兩個(gè)狀態(tài)過程的迭代,來逐步的準(zhǔn)確定位。
預(yù)測:當(dāng)前狀態(tài)環(huán)境下,對下一個(gè)時(shí)間段t的位置估計(jì)計(jì)算的值。
更新:更具傳感器獲取到比較準(zhǔn)確的位置信息后來更新當(dāng)前的預(yù)測問位置,也就是糾正預(yù)測的錯(cuò)誤。
你可能要問為什么有傳感器的數(shù)據(jù)了還要進(jìn)行更新?因?yàn)樵诂F(xiàn)實(shí)世界中傳感器是存在很多噪聲干擾的,所以也不能完全相信傳感器數(shù)據(jù)??柭惴ㄒ蕾囉诰€性計(jì)算,高斯分布,我們以一維定位來介紹算法的實(shí)現(xiàn)。
![](https://file1.elecfans.com//web2/M00/94/A5/wKgaomTl9puAaKTuAAA9CEpMGEI542.jpg)
![](https://file1.elecfans.com//web2/M00/94/A5/wKgaomTl9puAVRPGAAAqtHOQo0A955.png)
接下來我們開更新,預(yù)測后我們獲取到傳感器數(shù)據(jù),表示目前傳感器發(fā)現(xiàn)小車的位置應(yīng)該是在26這個(gè)位置,在這種情況下,我們肯定是覺得傳感器的準(zhǔn)確度比我之前的預(yù)測瞎猜要來的準(zhǔn)確。
所以方差自然會(huì)比較小,最終我們覺得真是的小車位置應(yīng)該是更靠近傳感器數(shù)據(jù)的,而且方差會(huì)縮小,以至于,想想也很清楚,我猜了一個(gè)預(yù)測值,現(xiàn)在有個(gè)專家告訴了我相對比較靠譜的數(shù)據(jù),那我對小車的位置的自信度肯定會(huì)上升啊。
最終小車的位置經(jīng)過這個(gè)時(shí)間段t的更新就是下圖紅色的高斯圖:
![](https://file1.elecfans.com//web2/M00/94/A5/wKgaomTl9puAcz2nAAB941nK9cw485.jpg)
就這樣不停的移動(dòng)更新,最終小車的位置就會(huì)越來越準(zhǔn)確。
一維模型下的Kalman公式:
預(yù)測
![](https://file1.elecfans.com//web2/M00/94/A5/wKgaomTl9pyAYSllAAAgzwT497s325.jpg)
更新
![](https://file1.elecfans.com//web2/M00/94/A5/wKgaomTl9pyAH5piAAAOhXX03Kw400.jpg)
![](https://file1.elecfans.com//web2/M00/94/A5/wKgaomTl9pyAXZZQAAAI67Vp6Is330.jpg)
參考代碼:
using namespace std;
double new_mean, new_var;
tuple<double, double> measurement_update(double mean1, double var1, double mean2, double var2)
{
new_mean = (var2 * mean1 + var1 * mean2) / (var1 + var2);
new_var = 1 / (1 / var1 + 1 / var2);
return make_tuple(new_mean, new_var);
}
tuple<double, double> state_prediction(double mean1, double var1, double mean2, double var2)
{
new_mean = mean1 + mean2;
new_var = var1 + var2;
return make_tuple(new_mean, new_var);
}
int main()
{
//Measurements and measurement variance
double measurements[5] = { 5, 6, 7, 9, 10 };
double measurement_sig = 4;
//Motions and motion variance
double motion[5] = { 1, 1, 2, 1, 1 };
double motion_sig = 2;
//Initial state
double mu = 0;
double sig = 1000;
for (int i = 0; i < sizeof(measurements) / sizeof(measurements[0]); i++) {
tie(mu, sig) = measurement_update(mu, sig, measurements[i], measurement_sig);
printf("update: [%f, %f] ", mu, sig);
tie(mu, sig) = state_prediction(mu, sig, motion[i], motion_sig);
printf("predict: [%f, %f] ", mu, sig);
}
return 0;
}
原文標(biāo)題:基礎(chǔ)卡爾曼濾波
文章出處:【微信公眾號:機(jī)器視覺智能檢測】歡迎添加關(guān)注!文章轉(zhuǎn)載請注明出處。
審核編輯:湯梓紅
-
濾波器
+關(guān)注
關(guān)注
161文章
7868瀏覽量
179021 -
算法
+關(guān)注
關(guān)注
23文章
4631瀏覽量
93406 -
原理
+關(guān)注
關(guān)注
4文章
550瀏覽量
44974
發(fā)布評論請先 登錄
相關(guān)推薦
卡爾曼濾波器的使用原理
卡爾曼濾波器是什么
卡爾曼濾波器參數(shù)分析與應(yīng)用方法研究
圖解卡爾曼濾波器
![圖解<b class='flag-5'>卡爾</b><b class='flag-5'>曼</b><b class='flag-5'>濾波器</b>](https://file.elecfans.com/web1/M00/45/E3/o4YBAFp6z72AYluhAAAgJNpnoxg334.jpg)
詳細(xì)了解卡爾曼濾波器的工作原理
![詳細(xì)了解<b class='flag-5'>卡爾</b><b class='flag-5'>曼</b><b class='flag-5'>濾波器</b>的<b class='flag-5'>工作原理</b>](https://file.elecfans.com/web1/M00/BA/D4/pIYBAF6c7Y2AP1Y-AABQ4FZh7jc501.png)
如何理解卡爾曼濾波器?卡爾曼濾波器狀態(tài)方程及測量方程
卡爾曼濾波器的特性及仿真
![<b class='flag-5'>卡爾</b><b class='flag-5'>曼</b><b class='flag-5'>濾波器</b>的特性及仿真](https://file1.elecfans.com/web2/M00/0A/C5/wKgZomcoQguAAeuwAAAAWwWCeQ0425.png)
評論