function?hmm?=?baum（hmm?samples）

mix??=?hmm.mix;			%高斯混合
N????=?length（mix）;		%HMM狀態數
K????=?length（samples）;	%語音樣本數
SIZE?=?size（samples（1）.data2）;?%參數階數

%?計算前向?后向概率矩陣?考慮多觀察序列和下溢問題
disp（‘計算樣本參數...‘）;
for?k?=?1:K
????fprintf（‘%d?‘k）
	param（k）?=?getparam（hmm?samples（k）.data）;
end
fprintf（‘\n‘）

%?重估轉移概率矩陣A:?trans
disp（‘重估轉移概率矩陣A...‘）
for?i?=?1:N-1
	denom?=?0;
	for?k?=?1:K
		tmp???=?param（k）.ksai（:i:）;
		denom?=?denom?+?sum（tmp（:））;
	end
	for?j?=?i:i+1
		nom?=?0;
		for?k?=?1:K
			tmp?=?param（k）.ksai（:ij）;
			nom?=?nom???+?sum（tmp（:））;
		end
		hmm.trans（ij）?=?nom?/?denom;
	end
end

%?重估混合高斯的參數
disp（‘重估混合高斯的參數...‘）
for?l?=?1:N
for?j?=?1:hmm.M（l）
	fprintf（‘%d%d?‘lj）
	%?計算各pdf的均值和方差
	nommean?=?zeros（1SIZE）;?
	nomvar??=?zeros（1SIZE）;?
	denom???=?0;
	for?k?=?1:K
		T?=?size（samples（k）.data1）;
		for?t?=?1:T
			x	????=?samples（k）.data（t:）;
			nommean?=?nommean?+?param（k）.gama（tlj）?*?x;
			nomvar??=?nomvar??+?param（k）.gama（tlj）?*?（x-mix（l）.mean（j:））.^2;
			denom???=?denom???+?param（k）.gama（tlj）;
		end
	end
%?	hmm.mix（l）.mean（j:）?=?nommean?/?denom;
	hmm.mix（l）.var?（j:）?=?nomvar??/?denom;

	%?計算各pdf的權
	nom???=?0;
	denom?=?0;
	for?k?=?1:K
		tmp?=?param（k）.gama（:lj）;????nom???=?nom???+?sum（tmp（:））;
		tmp?=?param（k）.gama（:l:）;????denom?=?denom?+?sum（tmp（:））;
	end
	hmm.mix（l）.weight（j）?=?nom/denom;
end
fprintf（‘\n‘）
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2012-12-14?20:08??hmm\
?????文件????????1513??2012-12-14?20:38??hmm\baum.m
?????文件????????2072??2012-12-14?20:08??hmm\enframe.m
?????文件????????1990??2002-02-28?11:05??hmm\getparam.m
?????文件???????63856??2001-11-12?12:37??hmm\hmm.mat
?????文件????????1266??2009-10-19?09:33??hmm\inithmm.asv
?????文件????????1266??2009-10-19?09:41??hmm\inithmm.m
?????文件?????????180??2012-12-14?20:01??hmm\main.m
?????文件????????3482??2012-12-14?20:06??hmm\melbankm.m
?????文件?????????859??2001-11-09?09:30??hmm\mfcc.m
?????文件?????????357??2001-11-11?10:28??hmm\mixture.m
?????文件?????????230??2001-11-10?21:48??hmm\pdf.m
?????文件?????????264??2009-10-19?08:52??hmm\recog.asv
?????文件?????????252??2009-10-18?17:25??hmm\recog.m
?????文件??????630000??2001-11-12?12:38??hmm\samples.mat
?????文件?????????839??2001-11-13?10:12??hmm\train.m
?????文件????????1860??2009-10-18?17:08??hmm\vad.m
?????文件????????1026??2001-11-11?09:21??hmm\viterbi.m