如何用Matlab实现最大期望值算法 (EM)？

qqsweet

快郁闷死了，做这个完全是赶鸭子上架。。。我来求救了$郁闷$
请问大侠们，下面这4个公式如何用Matlab实现？能不能给我段Matlab code
聚元重重有赏$害羞$
叩谢各位了

[ 本帖最后由 qqsweet 于 2008-10-21 15:18 编辑 ]

qqsweet

$送花$ $送花$

aileute

lz确认下，第一第二个公式有没有写错，你这两个公式里面怎么，怎么求p的有Q，求Q的里面有p呢？

qqsweet

原帖由 aileute 于 2008-10-5 17:55 发表
lz确认下，第一第二个公式有没有写错，你这两个公式里面怎么，怎么求p的有Q，求Q的里面有p呢？

谢谢回复$送花$
公式是对的
好像是要对p给个初值，然后进行迭代。。。迭代我不会$郁闷$
第一个公式里有个f，f是概率密度函数，也需要求出来才可以算Q，f 跟mu 和 sigma 有关

nightblue

http://www.mathworks.de/matlabce ... amp;objectType=File

看看这个能不能用到

[ 本帖最后由 nightblue 于 2008-10-6 12:02 编辑 ]

aileute

那lz要先把理论搞清楚才可以写，搞清楚了理论，编码就变得很简单了。

qqsweet

原帖由 nightblue 于 2008-10-6 12:00 发表
http://www.mathworks.de/matlabce ... amp;objectType=File

看看这个能不能用到

非常感谢$送花$ 我已经下了，希望可以找到点感觉

qqsweet

原帖由 aileute 于 2008-10-6 18:49 发表
那lz要先把理论搞清楚才可以写，搞清楚了理论，编码就变得很简单了。

我就是理论还没完全理出头绪，又碰到这么一个很纠结的公式:mad:
基本原理就是用一个高斯分布函数去拟合一个波形，从而探测出波峰，有几个波峰。
Beispiel:
图1： 黑色的表示接受到的波形
     淡绿色是对接受的波形进行了预处理，即除噪
     红色的是我们设定初值的符合正态分布的波形，根据上面给出的
     公式计算出最符合绿色波的曲线，结果如图2。

公式1 中的Q是概率，f是概率密度函数，j是Gauss komponente,在这个例子中j=5
公式中的S表示采样数，在这个例子中 S=130
Ni是每个采样对应的振幅
譬如 Ni = [2 4 5 6 7 18 12 8 5 31 44 ... ]

情况大概就是这样$m19$

[ 本帖最后由 qqsweet 于 2008-10-21 15:19 编辑 ]

qqsweet

会处理这类问题的大侠请帮帮忙$握手$

recbio

如果是单纯的高斯分布拟合, 可以用LM方法, 具体的 可以看
http://www.dolc.de/forum/viewthr ... &extra=page%3D2

有问题欢迎讨论
:)