中國十大權(quán)威新聞媒體谷歌優(yōu)化教程
李代數(shù)與李群的關(guān)系
R ˙ R T \dot{R}R^{T} R˙RT 是一個反對稱矩陣,所以這個矩陣可以用一個1×3向量進行反對稱來表示
R ˙ R T = Φ ^ \dot{R}R^{T}=Φ^{\hat{}} R˙RT=Φ^ ,
根據(jù)十四講 4.8 的推導(dǎo),最后則有 R ( t ) ˙ = Φ ^ ? R ( t ) \dot{R(t)}=Φ^{\hat{}}·R(t) R(t)˙?=Φ^?R(t)
這個李代數(shù) Φ Φ Φ 反映了 R R R 的導(dǎo)數(shù)性質(zhì),所以李代數(shù)是李群 S O ( 3 ) SO(3) SO(3) 的正切空間,因為這里李群是9維的,所以切向量也是一個空間
李代數(shù)由一個方向和夾角構(gòu)成, Φ = θ a , ∣ ∣ a ? ∣ ∣ = 1 Φ=θa,||\vec{a}||=1 Φ=θa,∣∣a∣∣=1
指數(shù)映射,也是羅德里格斯公式
e x p ( Φ ∧ ) = e x p ( θ a ∧ ) = c o s θ I + ( 1 ? c o s θ ) a a T + s i n θ a ∧ = R exp(Φ^{\wedge})=exp(θa^{\wedge})=cosθI+(1-cosθ)aa^{T}+sinθa^{\wedge}=R exp(Φ∧)=exp(θa∧)=cosθI+(1?cosθ)aaT+sinθa∧=R
e x p ( Φ ∧ ) = R exp(Φ^{\wedge})=R exp(Φ∧)=R
意思就是李代數(shù) Φ Φ Φ 可以由角度 θ θ θ 和方向向量 a ? \vec{a} a 表示,通過羅德里格斯公式可以變換成對應(yīng)的旋轉(zhuǎn)矩陣 S O ( 3 ) SO(3) SO(3)
用李代數(shù)表示旋轉(zhuǎn)會有個問題,就是周期性,就是多個李代數(shù)可以對應(yīng)一個旋轉(zhuǎn)矩陣,如果固定旋轉(zhuǎn)角度在±π時就是唯一對應(yīng)的
由旋轉(zhuǎn)矩陣 R R R 求李代數(shù) Φ Φ Φ , θ = a r c c o s t r ( R ) ? 1 2 , R a = a θ=arccos\frac{tr(R)-1}{2},Ra=a θ=arccos2tr(R)?1?,Ra=a,由 l n ( R ) ∨ = Φ ln(R)^{\vee}=Φ ln(R)∨=Φ 表示
BCH一階線性近似表達
l n ( e x p ( ? 1 ∧ ) e x p ( ? 2 ∧ ) ) ∨ ln(exp(\phi^{\wedge}_{1})exp(\phi^{\wedge}_{2}))^{\vee} ln(exp(?1∧?)exp(?2∧?))∨
當(dāng)左邊的 ? 1 \phi_{1} ?1? 為小量時,此時相當(dāng)于左乘,則
= > J l ( ? 2 ) ? 1 ? 1 + ? 2 =>J_{l}(\phi_{2})^{-1}\phi_{1}+\phi_{2} =>Jl?(?2?)?1?1?+?2?
其實就是相當(dāng)于在 ? 2 \phi_{2} ?2? 的基礎(chǔ)上加上微量 J l ( ? 2 ) ? 1 ? 1 J_{l}(\phi_{2})^{-1}\phi_{1} Jl?(?2?)?1?1?
當(dāng)右邊的 ? 2 \phi_{2} ?2? 為小量時,此時相當(dāng)于右乘,則
= > J l ( ? 1 ) ? 1 ? 2 + ? 1 =>J_{l}(\phi_{1})^{-1}\phi_{2}+\phi_{1} =>Jl?(?1?)?1?2?+?1?
SLAM中我們構(gòu)建了與位姿相關(guān)的函數(shù)后,需要討論該函數(shù)對位姿的求導(dǎo),以估計當(dāng)前值
有兩種方法求導(dǎo)
1、用李代數(shù)表示位姿,根據(jù)李代數(shù)加法對李代數(shù)求導(dǎo)
2、對李群進行左乘或右乘進行擾動,對擾動求導(dǎo)
由于使用李代數(shù)求導(dǎo)要計算 雅可比 J J J ,這個形式比較復(fù)雜,工程中不用這個方法,都是用擾動模型,所以只看擾動模型
? ( R p ) ? φ = l i m φ → 0 e x p ( φ ∧ ) e x p ( ? ∧ ) p ? e x p ( ? ∧ ) p φ \frac{\partial(Rp)}{\partial\varphi}=lim_{\varphi\rightarrow0}\frac{exp(\varphi^{\wedge})exp(\phi^{\wedge})p-exp(\phi^{\wedge})p}{\varphi} ?φ?(Rp)?=limφ→0?φexp(φ∧)exp(?∧)p?exp(?∧)p?
e x p ( φ ∧ ) exp(\varphi^{\wedge}) exp(φ∧) 是微量,相當(dāng)于是對旋轉(zhuǎn)的導(dǎo)數(shù),則等于 ( I + φ ∧ ) (I+\varphi^{\wedge}) (I+φ∧)
具體推導(dǎo)如下
左乘擾動求導(dǎo)
這里是旋轉(zhuǎn) R R R 對向量 p p p 進行旋轉(zhuǎn),不停地左乘擾動來改變向量 p p p 的方向,左擾動 Δ R \Delta{R} ΔR 對應(yīng)的李代數(shù) φ \varphi φ , R p Rp Rp 的結(jié)果也是向量
? ( R p ) ? φ = l i m φ → 0 e x p ( φ ∧ ) e x p ( ? ∧ ) p ? e x p ( ? ∧ ) p φ \frac{\partial(Rp)}{\partial\varphi}=lim_{\varphi\rightarrow0}\frac{exp(\varphi^{\wedge})exp(\phi^{\wedge})p-exp(\phi^{\wedge})p}{\varphi} ?φ?(Rp)?=limφ→0?φexp(φ∧)exp(?∧)p?exp(?∧)p?
= l i m φ → 0 ( I + φ ∧ ) e x p ( ? ∧ ) p ? e x p ( ? ∧ ) p φ =lim_{\varphi\rightarrow0}\frac{(I+\varphi^{\wedge})exp(\phi^{\wedge})p-exp(\phi^{\wedge})p}{\varphi} =limφ→0?φ(I+φ∧)exp(?∧)p?exp(?∧)p?
乘進去相減,很明顯
= l i m φ → 0 φ ∧ R p φ =lim_{\varphi\rightarrow0}\frac{\varphi^{\wedge}Rp}{\varphi} =limφ→0?φφ∧Rp?
叉乘有一個性質(zhì), a ? × b ? = ? b ? × a ? \vec{a}×\vec=-\vec×\vec{a} a×b=?b×a
a ? × b ? = a ∧ ? b \vec{a}×\vec=a^{\wedge}·b a×b=a∧?b
? b ? × a ? = ? b ∧ ? a -\vec×\vec{a}=-b^{\wedge}·a ?b×a=?b∧?a
則上式等于
= l i m φ → 0 ? ( R p ) ∧ φ φ =lim_{\varphi\rightarrow0}\frac{-(Rp)^{\wedge}\varphi}{\varphi} =limφ→0?φ?(Rp)∧φ?
約掉 φ \varphi φ 則 l i m φ → 0 = ? ( R p ) ∧ lim_{\varphi\rightarrow0}=-(Rp)^{\wedge} limφ→0?=?(Rp)∧
這里省去了雅可比 J J J 的計算
對右乘也是一樣的方法
矩陣轉(zhuǎn)置的性質(zhì)
( A + B ) T = A T + B T (A+B)^{T}=A^{T}+B^{T} (A+B)T=AT+BT
( λ A ) T = λ A T (\lambda{A})^{T}=\lambda{A^{T}} (λA)T=λAT
( A B ) T = B T A T (AB)^{T}=B^{T}A^{T} (AB)T=BTAT
SO(3)的伴隨性質(zhì)
R T E x p ( ? ) R = E x p ( R T ? ) R^{T}Exp(\phi)R=Exp(R^{T}\phi) RTExp(?)R=Exp(RT?)
? \phi ? 為擾動量 Δ R \Delta{R} ΔR 對應(yīng)的李代數(shù)
對復(fù)合旋轉(zhuǎn)進行求導(dǎo)
? L o g ( R 1 R 2 ) ? R 1 \frac{\partial Log(R_{1}R_{2})}{\partial R_{1}} ?R1??Log(R1?R2?)?
這里是對兩個相乘的旋轉(zhuǎn)矩陣中的其中一個旋轉(zhuǎn)進行求導(dǎo),上面的左乘擾動例子是對矩陣相乘向量進行求導(dǎo)的,所以可以直接對矩陣進行擾動
但是這里不能直接對矩陣 S O ( 3 ) SO(3) SO(3) 進行擾動,不能直接說 R 1 R 2 R_{1}R_{2} R1?R2? 對 R 1 R_{1} R1? 或 R 2 R_{2} R2? 的導(dǎo)數(shù),這樣就變成矩陣對向量的求導(dǎo),因為擾動量是可以用向量表示的,前面的例子 R p Rp Rp 相乘后也是個向量,所以可以直接對擾動量進行求導(dǎo),但是這里兩個矩陣相乘還是矩陣,所以得用 L o g Log Log 將矩陣相乘結(jié)果變?yōu)槔畲鷶?shù),這樣才符合求導(dǎo)的定義
L o g ( R ) = l o g ( R ) ∨ Log(R)=log(R)^{\vee} Log(R)=log(R)∨,用 L o g Log Log 就是為了懶得寫 ∨ \vee ∨
對 R 1 R_{1} R1? 進行右擾動
? L o g ( R 1 R 2 ) ? R 1 = l i m φ → 0 L o g ( R 1 E x p ( ? ) R 2 ) ? L o g ( R 1 R 2 ) ? \frac{\partial Log(R_{1}R_{2})}{\partial R_{1}}=lim_{\varphi\rightarrow0}\frac{Log(R_{1}Exp(\phi)R_{2})-Log(R_{1}R_{2})}{\phi} ?R1??Log(R1?R2?)?=limφ→0??Log(R1?Exp(?)R2?)?Log(R1?R2?)?
推導(dǎo)明天再寫了