mixture.m

来自「MATLAB扩展编程by 孙永康.ZIP.有需要的朋友请来下载.」· M 代码 · 共 19 行

function prob = mixture(mix, x)
%计算输出概率
%输入:
%  mix  -- 混合高斯结构
%  x    -- 输入向量, SIZE*1
%输出:
%  prob -- 输出概率

prob = 0;
for j = 1:mix.M
	m = mix.mean(j,:);
	v = mix.var (j,:);
	w = mix.weight(j);
	prob = prob + w * pdf(m, v, x);
end

% 加上realmin, 以防止viterbi.m中计算log(prob)时溢出
if prob==0, prob=realmin; end