【发布时间】：2015-07-10 20:50:29
【问题描述】：
我正在尝试使用 EM 算法估计高斯混合的平均值、权重和协方差，但是对于向量的所有值，我都得到了“NaN”（不是数字）。这是我的代码：
function [nw,ns,nu]=est(wo,so,uo,w,s,u,o,J)
% o is the data points vector (size=nx2)
% J is the number of Gaussian curve
% maximum number of Gaussians=5
% wo is the old weight vector (size=1x5) and "w" is the actual weight vector (size=1x5)
% so is the old covariance vector (size=1x5) and "s" is the actual covariance value (size=1x5)
% uo is the old mean vector (size=5x2) and "u" is the actual mean vector (size=5x2)
o=unique(o,'rows');
t=length(o);
for i=1:5
for j=1:J
un=0;ud=0;s2n=0;
for T=1:t
ud=ud+probab(wo,o(T,:),uo,so,j,J);
end
for T=1:t
un=un+probab(wo,o(T,:),uo,so,j,J)*(o(T,:));
end
for T=1:t
s2n=s2n+probab(wo,o(T,:),uo,so,j,J)*(o(T,:)-uo(j,:))*(((o(T,:)-uo(j,:))'));
end
wtemp(j)=ud/t;
utemp(j,:)=un/ud;
stemp(j)=sqrt(s2n/(ud));
end
wo=wtemp;
uo=utemp;
so=stemp;
end
nw=wo;
ns=so;
nu=uo;
函数“概率”是：
function [pro]=probab(w,o,u,s,j,J)
pro=w(j).*(1/(2*pi*s(j))).*exp((-1/(2.*(s(j).^2)))*((o(1)-u(j,1)).^2+(o(2)-u(j,2)).^2))/(gaussiandistribution(w,u,s,o(1),o(2),J));
函数“高斯分布”是
function [z]=gaussiandistribution(w,u,s,X,Y,J)
z1=w(1).*(1/(2*pi*s(1))).*exp(-(((X-u(1,1)).^2)/(2.*s(1).^2)+((Y-u(1,2)).^2)/(2.*s(1).^2)));
z2=w(2).*(1/(2*pi*s(2))).*exp(-(((X-u(2,1)).^2)/(2.*s(2).^2)+((Y-u(2,2)).^2)/(2.*s(2).^2)));
z3=w(3).*(1/(2*pi*s(3))).*exp(-(((X-u(3,1)).^2)/(2.*s(3).^2)+((Y-u(3,2)).^2)/(2.*s(3).^2)));
z4=w(4).*(1/(2*pi*s(4))).*exp(-(((X-u(4,1)).^2)/(2.*s(4).^2)+((Y-u(4,2)).^2)/(2.*s(4).^2)));
z5=w(5).*(1/(2*pi*s(5))).*exp(-(((X-u(5,1)).^2)/(2.*s(5).^2)+((Y-u(5,2)).^2)/(2.*s(5).^2)));
if J==5
z=z1+z2+z3+z4+z5;
elseif J==4
z=z1+z2+z3+z4;
elseif J==3
z=z1+z2+z3;
elseif J==2
z=z1+z2;
elseif J==1
z=z1;
end
【问题讨论】：
-
"但它不起作用" 这很模糊。代码是否运行？不会报错吗？如果是这样，错误信息是什么？如果没有，它在哪里失败？没有人会为您完成所有这些代码，您需要自己付出一些努力。
-
通过说“它不工作”，我的意思是它给了我所有向量的“NaN”，并且为了你的信息，我正在研究它两个月，所以我不期待语法错误或其他东西，我我不指望有人通过所有的代码。不幸的是，即使我彻底阅读了有关该主题的一些论文，我也没有任何统计知识来找出确切的问题。也许这是一个错误的地方问...
-
OK，遍历代码，找到第一个
NaN出现的位置，告诉我们什么变量，代码的哪个阶段。这将帮助您解决问题。 -
实际上问题出现在权重向量接近零（或非常小的值）时。当权重向量接近零时，“gaussiandistribution”函数的输出将变为零，这导致“probab”函数的输出变为零。如您所见，probab 函数在更新均值、权重和协方差向量方面起着最大的作用。我怀疑“ud”变量接近零并除以零导致“NaN”。
标签： algorithm matlab statistics