xie_matrix.m

来自「基于有限高斯混合模型的EM算法的源程序代码」· M 代码 · 共 12 行

function xie = xie_matrix(Z,X,U,m,n,k)
temp1 = zeros(m,m);
temp2 = 0;
for j=1:k
    for i=1:n
        temp1 = temp1 + Z(i,j).*(X(i,:)-U(j,:))'*(X(i,:)-U(j,:));
        temp2 = temp2 + Z(i,j);
    end
    xie(:,:,j) = temp1/temp2;
    temp1 = zeros(m,m);
    temp2 = 0;
end