gao.m

来自「em算法是一种估计最优参数的方法 又名最大期望算法」· M 代码 · 共 6 行

function G = Gao(X,U,xie,m,n,k)
for i = 1:n
    for j=1:k
        G(i,j) = exp((-1/2).*(X(i,:)- U(j,:))*pinv(xie(:,:,j))*(X(i,:)-U(j,:)).')*1/((2*pi).^(m/2)*sqrt(det(xie(:,:,j))));
    end 
end