楼主你好~我目前也遇到了这个问题。经过改错后已经可以运行：
问题原因：代码中没有处理singular问题
改正： 加几句代码，处理singular问题
把你第96行的function换成这个就好
function Px = calc_prob()
        Px = zeros(N, K);
        for k = 1:K
            Xshift = X-repmat(pMiu(k, : ), N, 1);%x-u
            lemda=1e-5;
            conv=pSigma(:, :, k)+lemda*diag(diag(ones(D)));%这里处理singular问题，为协方差矩阵加上一个很小lemda*I
            inv_pSigma = inv(conv);%协方差的逆
            tmp = sum((Xshift*inv_pSigma) .* Xshift, 2);%(X-U_k)sigma.*(X-U_k),tmp是个N*1的向量
            coef = (2*pi)^(-D/2) * sqrt(det(inv_pSigma));%前面的参数
            Px(:, k) = coef * exp(-0.5*tmp);%把数据点 x 带入到 Gaussian model 里得到的值
        end
    end
end

200。。请拿来啊~~

改成pinv
就行了
给我两百吧！

试过，
错误使用 svd
SVD 的输入不能包含 NaN 或 Inf。

出错 pinv (line 18)
[U,S,V] = svd(A,'econ');

出错 gmm/calc_prob (line 96)
            inv_pSigma = pinv(pSigma(:, :, k)+eye(D));

出错 gmm (line 40)
        Px = calc_prob();

悬赏还可以再增加，只要能解决问题！

我跑DSGE的过程中遇到了同样的问题，不知道是MATLAB版本的问题还是咋滴啊