虽然比较简单但还是记录一下吧,光看容易忘记
我们首先找到条件分布 p(xa∣xb) 的表达式。根据概率的乘法规则,我们可以从联合分布 p(x)=p(xa,xb) 出发,通过固定 xb 为观测值并归一化结果表达式来获得一个有效的概率分布。
为了更高效地得到解,可以考虑指数中的二次形式(由 (3.27) 给出),然后在计算结束时重新引入归一化系数。如果我们利用分块形式 (3.49),(3.50),和 (3.53),我们得到:
−21(x−μ)TΣ−1(x−μ)=
−21(xa−μa)TΛaa(xa−μa)−21(xa−μa)TΛab(xb−μb)
−21(xb−μb)TΛba(xa−μa)−21(xb−μb)TΛbb(xb−μb).(3.54)
我们看到,作为 xa 的函数,这仍然是一个二次形式,因此对应的条件分布 p(xa∣xb) 将是高斯分布。因为这个分布完全由其均值和协方差表征,我们的目标将是通过检查 (3.54) 来识别 p(xa∣xb) 的均值和协方差的表达式。
前提知识
-
多元高斯分布:
p(x)=(2π)D/2∣Σ∣1/21exp(−21(x−μ)TΣ−1(x−μ))
其中 x 是一个 D-维向量,μ 是均值向量,Σ 是协方差矩阵。
-
分块矩阵:
Σ=(ΣaaΣbaΣabΣbb)
其中 Σaa 和 Σbb 分别是子向量 xa 和 xb 的协方差矩阵,而 Σab 和 Σba 是它们之间的协方差矩阵。
-
精度矩阵:
Λ=Σ−1=(ΛaaΛbaΛabΛbb)
推导过程
1. 联合分布的二次形式
联合分布的指数部分可以写为:
−21(x−μ)TΣ−1(x−μ)
2. 分块表示
将 x 和 μ 分块表示:
x=(xa xb),μ=(μa μb)
3. 替换并展开
将上述分块矩阵代入指数部分:
−21(xa−μaxb−μb)T(ΛaaΛbaΛabΛbb)(xa−μaxb−μb)
4. 展开矩阵乘法
−21((xa−μa)TΛaa(xa−μa)+(xa−μa)TΛab(xb−μb)
+(xb−μb)TΛba(xa−μa)+(xb−μb)TΛbb(xb−μb))
5. 最终结果
整理后得到:
−21(x−μ)TΣ−1(x−μ)=
−21(xa−μa)TΛaa(xa−μa)−21(xa−μa)TΛab(xb−μb)
−21(xb−μb)TΛba(xa−μa)−21(xb−μb)TΛbb(xb−μb).
这就是公式 (3.54) 的推导过程。通过这种方式,我们可以进一步分析条件分布 p(xa∣xb) 的均值和协方差。