chmm.zip

chmm.zip

%   感谢亲亲使用此代码，此代码解决您的问题了吗~（@^_^@）~
%   没解决的话告诉亲亲一个好消息，登录淘宝店铺“大成软件工作室”，可以下载（????）1分钱成品代码（′▽‘〃）哦~
%   是的，亲亲真的没有看错，挠破头皮的问题真的1分钱就可以解决了（づ??????）づ
%   小的这就把传送门给您，记得要收藏好哦（づ￣3￣）づ╭?～
%   传送门：https://item.taobao.com/item.htm?spm=a1z10.1-c.w4004-15151018122.5.uwGoq5&id=538759553146
%   如果传送门失效，亲亲可以来店铺讨要，客服MM等亲亲来骚扰哦~（*/ω╲*）
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  2011-03-12 09:11  chmm\
     文件        1989  2016-09-22 15:44  chmm\baum.m
     文件        2681  2016-09-22 15:44  chmm\enframe.m
     文件        2464  2016-09-22 15:44  chmm\getparam.m
     文件       14700  2011-04-17 09:39  chmm\hmm.mat
     文件        1757  2016-09-22 15:44  chmm\inithmm.m
     文件        2948  2016-09-22 15:44  chmm\melbankm.m
     文件        1481  2016-09-22 15:44  chmm\mfcc.m
     文件     1444553  2010-11-30 19:35  chmm\mfcc.mat
     文件         831  2016-09-22 15:45  chmm\mixture.m
     文件         978  2016-09-22 15:45  chmm\mymfcc.m
     文件         704  2016-09-22 15:46  chmm\pdf.m
     文件         142  2010-11-21 22:18  chmm\read me.txt
     文件         375  2010-11-30 19:47  chmm\recog.asv
     文件         388  2011-02-28 15:57  chmm\recog.m
     文件    10475928  2011-04-17 09:31  chmm\samples.mat
     目录           0  2011-03-12 09:11  chmm\test\
     目录           0  2011-03-12 09:11  chmm\test1\
     目录           0  2011-03-12 09:11  chmm\test111\
     文件       28990  2010-05-25 21:35  chmm\test111\0.wav
     文件       26082  2010-05-25 21:36  chmm\test111\1.wav
     文件       23912  2010-05-25 21:36  chmm\test111\2.wav
     文件       33200  2010-05-25 21:37  chmm\test111\3.wav
     文件       45470  2010-05-25 21:38  chmm\test111\4.wav
     文件       46906  2010-04-28 21:03  chmm\test111\5.wav
     文件       46138  2010-04-28 21:03  chmm\test111\6.wav
     文件       46906  2010-04-28 21:03  chmm\test111\7.wav
     文件       41892  2010-05-25 21:39  chmm\test111\8.wav
     文件       43066  2010-04-28 21:04  chmm\test111\9.wav
     文件       46906  2009-08-27 21:33  chmm\test1\01.wav
     文件       64570  2009-08-27 21:35  chmm\test1\11.wav
............此处省略145个文件信息