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??2010-12-06?11:13??hmm聲音識別\
?????文件??????132344??2009-10-01?23:11??hmm聲音識別\a.wav
?????文件????????1527??2009-10-01?09:51??hmm聲音識別\baum.m
?????文件????????3363??2009-10-01?15:24??hmm聲音識別\DISTEUSQ.M
?????文件????????2072??2009-10-01?16:08??hmm聲音識別\ENframe.M
?????文件????????3030??2009-10-02?11:56??hmm聲音識別\frameDetect.asv
?????文件????????3073??2009-10-02?22:00??hmm聲音識別\frameDetect.m
?????文件????????3041??2009-10-02?22:00??hmm聲音識別\frameDetect1.m
?????文件????????1990??2002-02-26?15:51??hmm聲音識別\getparam.m
?????文件???????63856??2001-11-12?12:37??hmm聲音識別\hmm.mat
?????文件????????1266??2009-10-01?15:39??hmm聲音識別\inithmm.m
?????文件????????2001??2009-10-01?15:21??hmm聲音識別\KMEANS.M
?????文件????????1514??2010-12-08?21:35??hmm聲音識別\main.asv
?????文件????????1510??2010-08-03?19:16??hmm聲音識別\main.m
?????文件????????3482??2009-10-01?16:08??hmm聲音識別\MELBANKM.M
?????文件?????????859??2009-10-01?16:18??hmm聲音識別\mfcc.m
?????文件?????????357??2001-11-11?10:28??hmm聲音識別\mixture.m
?????文件?????????230??2001-11-10?21:48??hmm聲音識別\pdf.m
?????文件?????????574??2009-10-03?12:36??hmm聲音識別\recog.asv
?????文件?????????585??2009-10-03?13:08??hmm聲音識別\recog.m
?????文件??????630000??2001-11-12?12:38??hmm聲音識別\samples.mat
?????文件???????11151??2009-12-14?16:40??hmm聲音識別\soundrecognize.asv
?????文件???????11025??2009-12-23?10:35??hmm聲音識別\soundrecognize.fig
?????文件???????11688??2009-12-23?10:12??hmm聲音識別\soundrecognize.m
?????目錄???????????0??2009-11-03?11:10??hmm聲音識別\test\
?????文件???????44144??2005-06-22?20:44??hmm聲音識別\test\H0.wav
?????文件???????44144??2005-06-22?20:45??hmm聲音識別\test\H1.wav
?????文件???????44144??2005-06-22?20:48??hmm聲音識別\test\H2.wav
?????文件???????44144??2005-06-22?20:49??hmm聲音識別\test\H3.wav
?????文件???????44144??2005-06-22?20:49??hmm聲音識別\test\H4.wav
?????文件???????44144??2005-06-22?20:49??hmm聲音識別\test\H5.wav
............此處省略11個文件信息