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

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????10580??2004-11-16?20:19??whisper\0_01.wav

?????文件??????14444??2004-11-16?20:19??whisper\0_02.wav

?????文件??????14484??2004-11-16?20:19??whisper\0_03.wav

?????文件??????14554??2004-11-16?20:19??whisper\0_04.wav

?????文件??????13956??2004-11-16?20:19??whisper\0_05.wav

?????文件??????14376??2004-11-16?20:19??whisper\0_06.wav

?????文件??????14182??2004-11-16?20:19??whisper\0_07.wav

?????文件??????15242??2004-11-16?20:19??whisper\0_08.wav

?????文件??????15080??2004-11-16?20:19??whisper\0_09.wav

?????文件??????12766??2009-05-19?09:45??whisper\0_10.wav

?????文件??????12286??2009-05-19?09:45??whisper\0_11.wav

?????文件???????8926??2009-05-19?09:45??whisper\0_12.wav

?????文件???????9566??2009-05-19?09:45??whisper\0_13.wav

?????文件??????12126??2009-05-19?09:45??whisper\0_14.wav

?????文件??????11326??2009-05-19?09:45??whisper\0_15.wav

?????文件??????11006??2009-05-19?09:45??whisper\0_16.wav

?????文件???????9886??2009-05-19?09:45??whisper\0_17.wav

?????文件??????12286??2009-05-19?09:45??whisper\0_18.wav

?????文件???????8926??2009-05-19?09:45??whisper\0_19.wav

?????文件???????8926??2009-05-19?09:45??whisper\0_20.wav

?????文件??????10846??2009-05-19?09:45??whisper\0_21.wav

?????文件??????11966??2009-05-19?09:45??whisper\0_22.wav

?????文件???????8126??2009-05-19?09:45??whisper\0_23.wav

?????文件??????11006??2009-05-19?09:45??whisper\0_24.wav

?????文件???????8606??2009-05-19?09:45??whisper\0_25.wav

?????文件??????13458??2004-11-16?20:20??whisper\1_01.wav

?????文件??????14444??2004-11-16?20:20??whisper\1_02.wav

?????文件??????15574??2004-11-16?20:20??whisper\1_03.wav

?????文件??????15316??2004-11-16?20:20??whisper\1_04.wav

?????文件??????15664??2004-11-16?20:20??whisper\1_05.wav

............此處省略644個文件信息