用matlab编写一个用EM算法估计参数的程序

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

请问谁有EM算法估参的程序呢？要matlab的，高手帮帮忙哇。。。！！！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏2 回帖

关键词：MATLAB atlab matla EM算法 Lab matlab 程序

为何木有人回应。。。。？？

以前写过，不过换电脑没保存
去Mathworks的用户中心搜，应该有很多
http://www.mathworks.cn/matlabcentral/?s_tid=gn_mlc

windblood 发表于 2014-5-12 21:45
以前写过，不过换电脑没保存
去Mathworks的用户中心搜，应该有很多
http://www.mathworks.cn/matlabcentr ...

哦，谢谢你哈，我试试

using R it is worth reading this
http://epub.ub.uni-muenchen.de/11700/1/BA_Ernst.pdf

% 目标函数：高斯混合模N[x(i),mu,sm]=0.8N[x(i),mu1,sm1]+0.2N[x(i),mu2,sm2]
%                                           mu1=2,sm1=1;    mu2=1,sm2=1
clc
clear all
N=30;                               %变量x数目为30个
for i=1:N                            %构造目标函数
     if rand(1)>0.2    x(i)=randn(1)+2;
     else              x(i)=randn(1)+1;
     end
end
a(1,1)=0.8;mu(1,1)=3;sm(1,1)=2;      %给定初始值
a(2,1)=0.2;mu(2,1)=-2;sm(2,1)=7;
for n=1:30                          %EM迭代次数为30次
    for j=1:2                        %高斯分量2个
       for i=1:N
        p1(j)=a(j,n)*normpdf(x(i),mu(j,n),sm(j,n));
        P=0;  %清零
        for s=1:2
            P=a(s,n)*normpdf(x(i),mu(s,n),sm(s,n))+P;  %累加
        end
        p(j,i)=p1(j)/P;
       end
       temp1=0; temp2=0;temp3=0;
       for i=1:N
           temp1=p(j,i)+temp1;  %累加
           temp2=p(j,i)*x(i)+temp2;  %累加
           temp3=p(j,i)*((x(i)-(temp2/temp1)).^2)+temp3;  %累加
       end
       a(j,n+1)=temp1/N;  %得到后面的项
       mu(j,n+1)=temp2/temp1;  
       sm(j,n+1)=sqrt(temp3/temp1);   
    end
end

matlab-007 发表于 2015-2-25 21:17
% 目标函数：高斯混合模N[x(i),mu,sm]=0.8N[x(i),mu1,sm1]+0.2N[x(i),mu2,sm2]
%                       ...

这么就才看到，不好意思。
我的论文已经弄完了，不过，还是谢谢您。

空了的城堡 发表于 2015-3-31 15:24
这么就才看到，不好意思。
我的论文已经弄完了，不过，还是谢谢您。

你好  我论文也想用这个  能把你的程序发我参考一下吗 谢谢！ly892256813@sina.com

也给我来一份呗
Lpeach@163.com
QQQ

租金在做毕设，能不能发我一份，liaohuan5203@163.com