www.gusucode.com > SAE RBM 程序MATLAB源码代码实现的一个关于sae的例子 > minFunc/logistic/mylogsumexp.m
function lse = mylogsumexp(b) % does logsumexp across columns B = max(b,[],2); if issparse(b) lse = log(sum(exp(b-repmat(B,[1 size(b,2)])),2))+B; else lse = log(sum(exp(b-repmatC(B,[1 size(b,2)])),2))+B; end end