EM算法

一、前言

1. Jesen不等式

在这里插入图片描述

2、二项式分布的最大似然估计

在这里插入图片描述

3、高斯分布的极大似然估计

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

4、高斯混合分布

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

二、EM算法

1. EM算法的提出

在这里插入图片描述
在这里插入图片描述
迭代法:

  1. 找出一个下界函数 r ( x ∣ θ ) r(x|\theta) r(xθ) r ( x ∣ θ ) < = l ( θ ) r(x|\theta) <= l(\theta) r(xθ)<=l(θ)),在 x 0 x_{0} x0处于 l ( θ ) l(\theta) l(θ)相切.
  2. 找出下界函数 r ( x ∣ θ ) r(x|\theta) r(xθ)的极大值处 x 1 x_{1} x1.由于 r ( x ∣ θ ) < = l ( θ ) r(x|\theta) <= l(\theta) r(xθ)<=l(θ),所以 l ( θ , x 1 ) > l = l ( θ , x 0 ) l(\theta,x_{1})>l=l(\theta,x_{0}) l(θ,x1)>l=l(θ,x0)
  3. x 1 x_{1} x1处找出下界函数 r ( x ∣ θ ) r(x|\theta) r(xθ),重复上述过程,直到在选 r ( x ∣ θ ) r(x|\theta) r(xθ)时, r ( x 1 ∣ θ ) = l ( θ , x 1 ) = r ( x m a x ∣ θ ) r(x_{1}|\theta) = l(\theta,x_{1})=r(x_{max}|\theta) r(x1θ)=l(θx1)=r(xmaxθ).

这种方法与初值的选择有关,不同初值,最后得到的极大值不同。
在这里插入图片描述

2. 找出下界函数

在这里插入图片描述

  • 下界函数可以选择不等式右边的
    r ( x , z ; θ ) = Σ i = 1 m Σ z ( i ) Q i ( z i ) ) l o g p ( x ( i ) , z ( i ) ; θ ) Q i ( z i ) r(x,z;\theta) = \Sigma_{i=1}^{m}\Sigma_{z^{(i)}}Q_{i}(z^{i}))log\frac{p(x^{(i)},z^{(i)};\theta)}{Q_{i}(z^{i})} r(x,z;θ)=Σi=1mΣz(i)Qi(zi))logQi(zi)p(x(i),z(i);θ)
  • 找出最近下界.
    在相切点处, r ( x , z ; θ ) = l ( θ , x ) r(x,z;\theta) = l(\theta,x) r(x,z;θ)=l(θ,x) p ( x ( i ) , z ( i ) ; θ ) Q i ( z i ) = C \frac{p(x^{(i)},z^{(i)};\theta)}{Q_{i}(z^{i})}=C Qi(zi)p(x(i),z(i);θ)=C

为什么此处二者相等?(Q是Z的一个分布)
在这里插入图片描述
在这里插入图片描述
EM算法的整体框架

  1. 求隐变量Z的分布 Q i ( z i ) Q_{i}(z^{i}) Qi(zi)
  2. 求下界函数取极大值时的待估参数 θ \theta θ
    在这里插入图片描述

3. EM算法的收敛性

在这里插入图片描述

三、用EM算法推导混合高斯分布

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
求均值
在这里插入图片描述
在这里插入图片描述
求方差
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值