基于MATLAB和HMM的说话人识别系统

本资源是基于MATLAB的说话人识别系统，是基于HMM模型 的，具有很高的识别效率，大家可以放心使用

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

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      36552  2001-10-21 16:38  hmm\6a.wav

     文件      35624  2001-10-21 16:38  hmm\6b.wav

     文件      36280  2001-10-21 16:39  hmm\7a.wav

     文件      33576  2001-10-21 16:39  hmm\7b.wav

     文件      21624  2001-10-21 16:39  hmm\8a.wav

     文件      21608  2001-10-21 16:39  hmm\8b.wav

     文件      25192  2001-10-21 16:40  hmm\9a.wav

     文件      27192  2001-10-21 16:40  hmm\9b.wav

     文件       1512  2010-03-01 21:01  hmm\baum.asv

     文件       1515  2010-03-01 21:21  hmm\baum.m

     文件       3363  2000-02-14 08:32  hmm\DISTEUSQ.M

     文件       2072  1998-05-12 12:46  hmm\ENframe.M

     文件       1990  2002-02-26 15:51  hmm\getparam.m

     文件      63856  2001-11-12 12:37  hmm\hmm.mat

     文件       1266  2010-03-01 20:11  hmm\inithmm.m

     文件       2001  1998-07-27 15:06  hmm\KMEANS.M

     文件        165  2001-11-13 11:40  hmm\main.m

     文件         83  2010-03-01 21:17  hmm\main1.m

     文件       3482  1998-05-12 15:22  hmm\MELBANKM.M

     文件        859  2010-02-27 11:11  hmm\mfcc.asv

     文件        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

     文件        390  2010-03-01 20:23  hmm\pre.m

     文件        245  2010-03-01 22:41  hmm\recog.asv

     文件        239  2010-03-01 22:46  hmm\recog.m

     文件        239  2010-03-02 08:42  hmm\recog1.m

     文件     630000  2001-11-12 12:38  hmm\samples.mat

     文件      29096  2001-10-21 16:33  hmm\source0.wav

     文件      33832  2001-10-21 16:34  hmm\source1.wav

............此处省略135个文件信息