【视觉slam十四讲】【十二讲 建图】12.1 习题:证明两个正态分布的联合分布
习题证明:证明公式12.6。
其实就是证明两个正态分布的联合分布是什么,即附录A的A.6。
1. 两个高斯 PDF 的乘积
设两个正态分布:
p1(x)=12πσ12exp[−(x−μ1)22σ12] p_1(x) = \frac{1}{\sqrt{2\pi\sigma_1^2}} \exp\left[ -\frac{(x-\mu_1)^2}{2\sigma_1^2} \right] p1(x)=2πσ121exp[−2σ12(x−μ1)2]
p2(x)=12πσ22exp[−(x−μ2)22σ22] p_2(x) = \frac{1}{\sqrt{2\pi\sigma_2^2}} \exp\left[ -\frac{(x-\mu_2)^2}{2\sigma_2^2} \right] p2(x)=2πσ221exp[−2σ22(x−μ2)2]
它们的乘积(忽略归一化常数)为:
p1(x)p2(x)∝exp[−(x−μ1)22σ12−(x−μ2)22σ22] p_1(x) p_2(x) \propto \exp\left[ -\frac{(x-\mu_1)^2}{2\sigma_1^2} - \frac{(x-\mu_2)^2}{2\sigma_2^2} \right] p1(x)p2(x)∝exp[−2σ12(x−μ1)2−2σ22(x−μ2)2]
2. 合并指数项
指数部分为:
Q=−12[(x−μ1)2σ12+(x−μ2)2σ22] Q = -\frac{1}{2} \left[ \frac{(x-\mu_1)^2}{\sigma_1^2} + \frac{(x-\mu_2)^2}{\sigma_2^2} \right] Q=−21[σ12(x−μ1)2+σ22(x−μ2)2]
展开平方项:
(x−μ1)2σ12+(x−μ2)2σ22=x2−2μ1x+μ12σ12+x2−2μ2x+μ22σ22 \frac{(x-\mu_1)^2}{\sigma_1^2} + \frac{(x-\mu_2)^2}{\sigma_2^2} = \frac{x^2 - 2\mu_1 x + \mu_1^2}{\sigma_1^2} + \frac{x^2 - 2\mu_2 x + \mu_2^2}{\sigma_2^2} σ12(x−μ1)2+σ22(x−μ2)2=σ12x2−2μ1x+μ12+σ22x2−2μ2x+μ22
=(1σ12+1σ22)x2−2(μ1σ12+μ2σ22)x+(μ12σ12+μ22σ22) = \left( \frac{1}{\sigma_1^2} + \frac{1}{\sigma_2^2} \right) x^2 - 2\left( \frac{\mu_1}{\sigma_1^2} + \frac{\mu_2}{\sigma_2^2} \right) x + \left( \frac{\mu_1^2}{\sigma_1^2} + \frac{\mu_2^2}{\sigma_2^2} \right) =(σ121+σ221)x2−2(σ12μ1+σ22μ2)x+(σ12μ12+σ22μ22)
3. 写成高斯分布的形式
一个高斯分布 N(x∣μf,σf2)\mathcal{N}(x \mid \mu_f, \sigma_f^2)N(x∣μf,σf2) 的指数部分为(fff表示fusefusefuse的缩写):
−(x−μf)22σf2=−12σf2(x2−2μfx+μf2) -\frac{(x - \mu_f)^2}{2\sigma_f^2} = -\frac{1}{2\sigma_f^2} \left( x^2 - 2\mu_f x + \mu_f^2 \right) −2σf2(x−μf)2=−2σf21(x2−2μfx+μf2)
比较 x2x^2x2 和 xxx 的系数。
从我们展开的 QQQ 来看:
Q=−12[(1σ12+1σ22)x2−2(μ1σ12+μ2σ22)x+常数项] Q = -\frac12 \left[ \left( \frac{1}{\sigma_1^2} + \frac{1}{\sigma_2^2} \right) x^2 - 2\left( \frac{\mu_1}{\sigma_1^2} + \frac{\mu_2}{\sigma_2^2} \right) x + \text{常数项} \right] Q=−21[(σ121+σ221)x2−2(σ12μ1+σ22μ2)x+常数项]
与 −12σf2(x2−2μfx+μf2)-\frac{1}{2\sigma_f^2}(x^2 - 2\mu_f x + \mu_f^2)−2σf21(x2−2μfx+μf2) 比较:
- x2x^2x2 的系数:
1σ12+1σ22=1σf2 \frac{1}{\sigma_1^2} + \frac{1}{\sigma_2^2} = \frac{1}{\sigma_f^2} σ121+σ221=σf21
所以:
1σf2=1σ12+1σ22 \frac{1}{\sigma_f^2} = \frac{1}{\sigma_1^2} + \frac{1}{\sigma_2^2} σf21=σ121+σ221
这就是 精度相加 的公式。
- xxx 的系数:
μ1σ12+μ2σ22=μfσf2 \frac{\mu_1}{\sigma_1^2} + \frac{\mu_2}{\sigma_2^2} = \frac{\mu_f}{\sigma_f^2} σ12μ1+σ22μ2=σf2μf
所以:
μf=μ1σ12+μ2σ221σ12+1σ22 \mu_f = \frac{ \frac{\mu_1}{\sigma_1^2} + \frac{\mu_2}{\sigma_2^2} }{ \frac{1}{\sigma_1^2} + \frac{1}{\sigma_2^2} } μf=σ121+σ221σ12μ1+σ22μ2
4. 最终融合公式
融合后的方差:
σf2=11σ12+1σ22 \sigma_f^2 = \frac{1}{\frac{1}{\sigma_1^2} + \frac{1}{\sigma_2^2}} σf2=σ121+σ2211
σf2=σ12σ22σ12+σ22 \sigma_f^2 = \frac{\sigma_1^2 \sigma_2^2}{\sigma_1^2 + \sigma_2^2} σf2=σ12+σ22σ12σ22
融合后的均值:
μf=μ1σ12+μ2σ221σ12+1σ22 \mu_f = \frac{ \frac{\mu_1}{\sigma_1^2} + \frac{\mu_2}{\sigma_2^2} }{ \frac{1}{\sigma_1^2} + \frac{1}{\sigma_2^2} } μf=σ121+σ221σ12μ1+σ22μ2
μf=μ1σ22+μ2σ12σ12+σ22 \mu_f = \frac{\mu_1 \sigma_2^2 + \mu_2 \sigma_1^2}{\sigma_1^2 + \sigma_2^2} μf=σ12+σ22μ1σ22+μ2σ12
因此融合后的分布为:
pfuse(x)=N(x | μ1σ22+μ2σ12σ12+σ22, σ12σ22σ12+σ22) p_{\text{fuse}}(x) = \mathcal{N}\left( x \,\middle|\, \frac{\mu_1 \sigma_2^2 + \mu_2 \sigma_1^2}{\sigma_1^2 + \sigma_2^2},\; \frac{\sigma_1^2 \sigma_2^2}{\sigma_1^2 + \sigma_2^2} \right) pfuse(x)=N(xσ12+σ22μ1σ22+μ2σ12,σ12+σ22σ12σ22)
