久久国产成人av_抖音国产毛片_a片网站免费观看_A片无码播放手机在线观看,色五月在线观看,亚洲精品m在线观看,女人自慰的免费网址,悠悠在线观看精品视频,一级日本片免费的,亚洲精品久,国产精品成人久久久久久久

分享

超全,!SLAM中的濾波算法大全

 InfoRich 2023-02-07 發(fā)布于江蘇

快速獲得最新干貨

本文轉(zhuǎn)載自自動駕駛與AI

概率論基礎(chǔ)

邊緣概率

邊緣概率是相對于聯(lián)合概率而言的的,,雖然你有兩個變量(x,y)但是你可以只考慮x或者y的分部,好像另外一個不存在一樣,,寫作 P(x) 或者 P(y) ,。

Image

離散概率和為1,即:

Image

為了簡化符號,,在可能時通常省略隨機變量的明確表示,,而是使用常見的縮寫 P(x) 代替 P(X=x) 。

聯(lián)合概率和獨立

兩個隨機變量X和Y的聯(lián)合概率表示兩件事都發(fā)生的概率,,由下式給出:

Image

如果X和Y相互獨立,,有 P(x,y)=P(x)P(y).

獨立與條件獨立

Image

左邊是獨立的定義,,右邊是條件獨立的定義??梢钥闯觯簵l件獨立是指在某個條件已知的前提下獨立,。

注意:獨立與條件獨立不能互推。

https://zhuanlan.zhihu.com/p/58593725

條件獨立的理解

如果 A,B 關(guān)于 C 條件獨立:C 的發(fā)生,,使本來不獨立的 A 和 B 變得獨立起來,,也即事件 C 的發(fā)生,解開了 A 和 B 的依賴關(guān)系,;

條件概率

假定已經(jīng)知道Y=y,,想知道基于以上事實條件X=x的概率。這樣的概率表示為

Image

稱為條件概率,。

如果p(y)>0,則條件概率定義為 .

Image

通過統(tǒng)計方法分析事件發(fā)生關(guān)系的占比,,以做到已知結(jié)果倒推原因,。

如果X和Y相互獨立,則有 .

Image

聯(lián)合概率,、邊緣概率與條件概率之間的關(guān)系

Image

將概率轉(zhuǎn)化為面積來理解:

聯(lián)合概率P(X=a,Y=b)滿足X=a且Y=b的面積 邊緣概率P(X=a)不考慮Y的取值,,所有滿足X=a的區(qū)域的總面積 條件概率P(X=a|Y=b)在Y=b的前提下,滿足X=a的面積(比例)

Image
三種概率的關(guān)系

全概率定理

從條件概率和概率測量公理得出的一個有趣事實經(jīng)常被稱為全概率定理

Image

貝葉斯定理

https://www.cnblogs.com/HuZihu/p/9368355.html

該定理將條件概率p(x|y)與其“逆”概率p(ylx)聯(lián)系起來,。如此處所闡述的,,準則要求p(y)>0:

Image

如果x是一個希望由y推測出來的數(shù)值,其中y稱為數(shù)據(jù)(data)(結(jié)果),,也就是傳感器測量值,。P(x) 是先驗概率(prior probability), P(y|x) 是條件概率(conditional probability)(原因),, P(x}y) 是后驗概率(posterior probability),。

SLAM中很多理論,諸如濾波方法,、BA和圖優(yōu)化等,,都是基于貝葉斯定理的思想而來的。

貝葉斯濾波

貝葉斯濾波基于馬爾可夫假設(shè),。

馬爾可夫假設(shè)

馬爾可夫假設(shè):系統(tǒng)當前時刻的狀態(tài)只與上一個時刻有關(guān),,與之前任意時刻的狀態(tài)都沒有關(guān)系。當然,,馬爾可夫假設(shè)的成立是有條件的,,它要求狀態(tài)變量具有完整性,換句話說,,系統(tǒng)的全部信息都包含在了 中,,不存在其它隨機變量對系統(tǒng)狀態(tài)產(chǎn)生影響,。

公式展示

Image
Image

公式推導

證明來自《概率機器人》p23~25頁,用到了馬爾可夫假設(shè),、全概率定理,、條件獨立、貝葉斯定理等,。

Image

卡爾曼濾波(KF和EKF)

https://zhuanlan.zhihu.com/p/558449409

EKF局限性

所以EKF面臨的一個重要問題是,,當一個高斯分布經(jīng)過非線性變換后,如何用另一個高斯分布近似它,?按照它現(xiàn)在的做法,,存在以下的局限性:(注意是濾波器自己的局限性,還沒談在SLAM問題里的局限性),。

  1. 即使是高斯分布,,經(jīng)過一個非線性變換后也不是高斯分布。EKF只計算均值與協(xié)方差,,是在用高斯近似這個非線性變換后的結(jié)果,。(實際中這個近似可能很差)。
  2. 系統(tǒng)本身線性化過程中,,丟掉了高階項,。
  3. 線性化的工作點往往不是輸入狀態(tài)真實的均值,而是一個估計的均值,。于是,,在這個工作點下計算的,也不是最好的,。
  4. 在估計非線性輸出的均值時,,EKF算的是的形式。這個結(jié)果幾乎不會是輸出分布的真正期望值,。協(xié)方差也是同理,。

解決辦法

那么,怎么克服以上的缺點呢,?途徑很多,,主要看我們想不想維持EKF的假設(shè)。如果我們只是希望維持高斯分布假設(shè),,可以這樣子改:

  1. 為了克服第2條,,可以使用ESKF。它的線性化是總是在0附近,,并且雅可比矩陣的形式簡單,,計算迅速。因此線性化比EKF更準確,。雖然,,ESKF主要用于有imu的系統(tǒng),。
  2. 為了克服第3條工作點的問題,我們以EKF估計的結(jié)果為工作點,,重新計算一遍EKF,,直到這個工作點變化夠小。是為迭代EKF(Iterated EKF, IEKF),。
  3. 為了克服第4條,,我們除了計算 , f()再計算其他幾個精心挑選的采樣點,,然后用這幾個點估計輸出的高斯分布,。是為Sigma Point KF(SPKF,或UKF),。
  4. 如果我們不要高斯分布假設(shè),,憑什么要用高斯去近似一個長得根本不高斯的分布呢?于是問題變?yōu)?,丟掉高斯假設(shè)后,,怎么描述輸出函數(shù)的分布就成了一個問題。一種比較暴力的方式是:用足夠多的采樣點,,來表達輸出的分布。這種蒙特卡洛的方式,,也就是粒子濾波的思路,。
  5. 如果再進一步,可以丟棄濾波器思路,,說:為什么要用前一個時刻的值來估計下一個時刻呢,?我們可以把所有狀態(tài)看成變量,把運動方程和觀測方程看成變量間的約束,,構(gòu)造誤差函數(shù),,然后最小化這個誤差的二次型。這樣就會得到非線性優(yōu)化的方法,,在SLAM里就走向圖優(yōu)化那條路上去了,。不過,非線性優(yōu)化也需要對誤差函數(shù)不斷地求梯度,,并根據(jù)梯度方向迭代,,因而局部線性化是不可避免的。

無跡卡爾曼濾波(UKF)

https://blog.csdn.net/gangdanerya/article/details/105215446

該算法的核心思想是:采用無損變換(Unscented Transform,,UT),,利用一組Sigma采樣點來描述隨機變量的高斯分布,然后通過非線性函數(shù)的傳遞,,再利用加權(quán)統(tǒng)計線性回歸技術(shù)來近似非x線性函數(shù)的后驗均值和方差,。

相比于EKF,,UKF的估計精度能夠達到泰勒級數(shù)展開的二階精度。其本質(zhì)是:近似非線性函數(shù)的均值和方差遠比近似非線性函數(shù)本身更容易,。

UT變換

y=f(x)是一個非線性變換, x為 n 維隨機變量,,UT變換根據(jù)一定的采樣策略,獲得一組Sigma采樣點,,并設(shè)定均值權(quán)值W和方差權(quán)值,,來近似非線性函數(shù)的后驗均值和方差。利用選取的Sigma采樣點集進行非線性函數(shù)傳遞可得:

Image

其中,,y_i為Sigma采樣經(jīng)過非線性函數(shù)傳遞后對應的點,。根據(jù)加權(quán)統(tǒng)計線性回歸技術(shù),可以近似得到y(tǒng)的統(tǒng)計特性:

Image

Sigma點采樣策略

下面介紹兩種經(jīng)常使用的采樣策略:比例采樣和比例修正對稱采樣,。

Image
Image
Image
Image

根據(jù)Sigma點采樣策略不同,,相應的Sigma點以及均值權(quán)值和方差權(quán)值也不盡相同,因此UT變換的估計精度也會有差異,。但總體來說,,其估計精度能夠達到泰勒級數(shù)展開的二階精度。

使用參數(shù)beta對高斯表示的附加的分布信息進行編碼,。如果分布是精確的高斯分布,,則beta=2是最佳的選擇。

Image

其中11行,,參考 Will:EKF公式推導 中的證明的協(xié)方差的其他形式,。

UKF采樣原則

為保證隨機變量x經(jīng)過采樣之后得到的Sigma采樣點仍具有原變量的必要特性,所以采樣點的選取應滿足:

Image

式中{}內(nèi)的符號分別為Sigma采樣點,,均值權(quán)值和方差權(quán)值組成的集合,;L為采樣點個數(shù),P為隨機變量x的密度函數(shù),,g[~]確定x的相關(guān)信息,。如果密度函數(shù)P(x)只有一、二階矩時,,可以寫成:

Image
Image

UKF的特點

UKF相比于EKF的精度更高一些,,其精度相當于二階泰勒展開,但速度會略慢一點,。UKF另一個巨大優(yōu)勢是不需要計算雅克比矩陣,,而有些時候雅克比矩陣也確實的我們無法獲得的。

另外UKF與PF(粒子濾波)也有相似之處,,只是無跡變換中選擇的粒子是明確的,,而粒子濾波中的粒子是隨機的。隨機的好處是可以用于任意分布,但也具有其局限性,。因此對于分布近似為高斯分布的,,采用UKF將更有效。

信息濾波(IF)

卡爾曼濾波(KF)的對偶濾波算法就是信息濾波(Information Filter,,IF),,形式和KF差不多。

其中,,用信息向量 信息矩陣 表示KF中的 .轉(zhuǎn)換關(guān)系為:

Image

Error-state卡爾曼濾波(ESKF)

ESKF(error state Kalman filter)是Kalman濾波的一種特殊形式,,采用ESKF的原因是由于對誤差的線性化要比直接對函數(shù)進行線性化更加符合實際情況。如下圖:

Image

狀態(tài)誤差卡爾曼(ESKF)的應用,,它是卡爾曼濾波器的變種中應用最為廣泛的一種,,與EKF一樣,它也是一種針對時變系統(tǒng)的非線性濾波器,。但是與EKF不同的是,,它的線性化是總是在0附近,因此線性化更準確,。絕大部分的場景,,ESKF就足夠使用了。如果對于濾波有更高的要求,,可以選擇UKF,,甚至PF。

ESKF數(shù)學模型

Error-state Kalman Filter(ESKF)就是一種傳感器融合的算法,,它的基礎(chǔ)仍然是卡爾曼濾波,。它的核心思想是把系統(tǒng)的狀態(tài)分為三類:

  • true state:實際狀態(tài),系統(tǒng)實際的運行狀態(tài)
  • nominal state:名義狀態(tài),,描述了運動狀態(tài)的主要趨勢,主導成分,。(large-signal,,非線性)
  • error state:誤差狀態(tài),實際狀態(tài)與名義狀態(tài)之間的差值(small-signal,,近似線性,,適合線性高斯濾波)

基于以上狀態(tài)分類,我們可以將關(guān)心的true-state,,分為兩部分,,分別進行估計,即nominal stateerror state,,然后再進行二者疊加,。

Image

ESKF的全過程

  • 對IMU數(shù)據(jù)進行積分,獲得名義狀態(tài)X,,注意這個X并沒有考慮噪聲,,所以必然引入了累計誤差,;
  • 利用KF算法估計Error State,包括狀態(tài)更新和測量更新,,這個過程是考慮了噪聲的,,而且由于這個誤差狀態(tài)的方程式近似線性的,直接使用KF就可以,;
  • 將Error State加到Nominal State上,,獲取“真值”;
  • 重置Error State,,等待下一次更新,;

ESKF的優(yōu)勢

Error State中的參數(shù)數(shù)量與運動自由度是相等的,避免了過參數(shù)化引起的協(xié)方差矩陣奇異的風險,;
Error State總是接近于0,,Kalman Filter工作在原點附近。因此,,遠離奇異值,、萬向節(jié)鎖,并且保證了線性化的合理性和有效性,;
Error State總是很小,,因此二階項都可以忽略,因此雅可比矩陣的計算會很簡單,,很迅速,;
Error State的變化平緩,因此KF修正的頻率不需要太高,。

Image

公式推導

https://zhuanlan.zhihu.com/p/557496107

迭代卡爾曼濾波(IKF/IEKF)

多狀態(tài)約束下的卡爾曼濾波(MSCKF)

粒子濾波(PF)

https://blog.csdn.net/weixin_38145317/article/details/82852580

總結(jié)和比較

卡爾曼濾波和其他方法的比較

  1. 卡爾曼濾波是遞歸的線性高斯系統(tǒng)最優(yōu)估計,。
  2. EKF將NLNG系統(tǒng)在工作點附近近似為LG進行處理。
  3. IEKF對工作點進行迭代,。
  4. UKF沒有線性化近似,,而是把sigma point進行非線性變換后再用高斯近似。
  5. PF去掉高斯假設(shè),,以粒子作為采樣點來描述分布,。
  6. 優(yōu)化方式同時考慮所有幀間約束,迭代線性化求解,。

濾波和優(yōu)化方法的區(qū)別

Image
EKF的局限性

濾波方法(EKF):

(1)濾波方法具有一階馬爾可夫假設(shè),,k時刻的狀態(tài)只與k-1時刻相關(guān)。

(2)可以粗略認為EKF僅是優(yōu)化中的一次迭代,,除非是使用IKF/IEKF

(3)計算量小,,實時性較好,容易積累誤差。

優(yōu)化方法(BA):

(1)沒有馬爾可夫假設(shè),,更傾向于使用所有的歷史數(shù)據(jù),。稱為全體時間上的SLAM(Full-SLAM),計算量變大

(2)沒有迭代次數(shù)的限制

(3)累計誤差少,。最理想的情況是回環(huán)檢測之后,,可以消除整個軌跡的誤差。

(4)計算量大,。不過目前基于雅各比和海森矩陣稀疏性的研究,,使得計算資源的消耗大大降低。

下文采自:Sky Shaw:多傳感器融合 | 各類濾波器方法整理 https://zhuanlan.zhihu.com/p/542440266

卡爾曼濾波器:從k-1時刻后驗推k時刻先驗,,從k時刻先驗推k時刻后驗(馬爾可夫假設(shè)),;
擴展卡爾曼濾波器:對卡爾曼濾波器進行修正,針對不是線性的情況,,采用一階泰勒展開近似線性(馬爾可夫假設(shè)),;
BA優(yōu)化:把一路上的所有坐標點(像素坐標對應的空間點等)與位姿整體放在一起作為自變量進行非線性優(yōu)化
PoseGraph優(yōu)化:先通過一路遞推方式算出的各點位姿,通過數(shù)學方式計算得到一個位姿的變換A,,再通過單獨拿出兩張圖像來算出一個位姿變換B,,爭取讓B=A
增量因子圖優(yōu)化:保留中間結(jié)果,每加入一個點,,對不需要重新計算的就直接用之前的中間結(jié)果,,需要重新計算的再去計算,從而避免冗余計算,。iSAM是增量的處理后端優(yōu)化,。由于機器人是運動的,不同的邊和節(jié)點會被不斷加入圖中,。每加入一個點,,普通圖優(yōu)化是對整個圖進行優(yōu)化,所以比較麻煩和耗時,。iSAM的話相當于是保留中間結(jié)果,,每加入一個點,對不需要重新計算的就直接用之前的中間結(jié)果,,對需要重新計算的再去計算。以這種方式加速計算,,避免冗余計算,。

從全文的整理來看,特別是在IKF章節(jié),,我們可以發(fā)現(xiàn)濾波算法其實就是將Sliding Window的大小設(shè)置為1的優(yōu)化算法,,不論是優(yōu)化算法還是濾波算法都是期望求解出問題概率模型的最大似然估計(MLE),本質(zhì)上濾波就是基于馬爾可夫假設(shè)將優(yōu)化問題的建模進行了“特征化”的處理。再進一步分析,,正因為濾波器進行了范圍上維度的“簡化”,、模型的近似處理,所以其計算消耗較于優(yōu)化算法更低,,但由此引發(fā)的代價就是精度上的損失,。另一個需要考慮的方面是,濾波器是由先驗信息+運動模型+觀測信息三個方面點順序執(zhí)行,,以實現(xiàn)位姿狀及其協(xié)方差的估計和更新,,也正因為濾波器的框架如此,若是先驗(上一時刻)狀態(tài)出現(xiàn)了問題,,比如位姿跟丟,、計算錯誤等,那么在該時刻之后的狀態(tài)都會出現(xiàn)問題以致糾正不回來了,,而基于優(yōu)化方法的位姿求解則會考慮更多時刻(關(guān)鍵幀)下的狀態(tài)信息和觀測信息,,即使有某一時刻的狀態(tài)量和協(xié)方差是outlier,系統(tǒng)也有一定的能力維穩(wěn),。

如果在處理器算力充足且精度為第一需求的前提下,,那么在位姿估算問題處理上是首推優(yōu)化算法,但若是效率是第一前提條件,,那么就需要根據(jù)實際的應用情況和機器人問題模型選擇合適的濾波器了,。

參考

Sky Shaw:多傳感器融合 | 各類濾波器方法整理 https://zhuanlan.zhihu.com/p/542440266

SLAM中的EKF,UKF,,PF原理簡介 - 半閑居士 - 博客園 https://www.cnblogs.com/gaoxiang12/p/5560360.html

SLAM 中的 Kalman Filter 推導
https://www./38267

學機器人SLAM/3D視覺,,就上cvlife.net 

點擊領(lǐng)取學習資料 → 機器人SLAM學習資料大禮包

獨家重磅課程官網(wǎng):cvlife.net

全國最大的機器人SLAM開發(fā)者社區(qū)

技術(shù)交流群

—   版權(quán)聲明  —

    本站是提供個人知識管理的網(wǎng)絡存儲空間,,所有內(nèi)容均由用戶發(fā)布,,不代表本站觀點。請注意甄別內(nèi)容中的聯(lián)系方式,、誘導購買等信息,,謹防詐騙。如發(fā)現(xiàn)有害或侵權(quán)內(nèi)容,,請點擊一鍵舉報,。
    轉(zhuǎn)藏 分享 獻花(0

    0條評論

    發(fā)表

    請遵守用戶 評論公約

    類似文章 更多