matlab语音识别系统

1、智能仪器课程设计题 目：MATLAB实现语音识别功能班 级：学 号：姓 名：同组人员：任课教师：完成时间：2012/11/3目录一、设计任务及要求 1二、语音识别的简单介绍2.1语者识别的概念22.2 特征参数的提取 32.3用矢量量化聚类法生成码本 32.4VQ的说话人识别 4三、算法程序分析3.1函数关系.43.2 代码说明 53.2.1 函数mfcc 53.2.2 函数disteu 53.2.3 函数 vqlbg .63.2.4 函数test 63.2.5 函数testDB 73.2.6 函数 train 83.2.7 函数melfb 8四、演示分析 .9五、心得体会 .11附：GUI程

2、序代码12>设计任务及要求用MATLAB实现简单的语音识别功能；具体设计要求如下：用MATLAB实现简单的数字19的语音识别功能 、语音识别的简单介绍基于VQ的说话人识别系统，矢量量化起着双重作用。在训练阶段，把每一个 说话者所提取的特征参数进行分类， 产生不同码字所组成的码本。在识别（匹配） 阶段，我们用VQ方法计算平均失真测度（本系统在计算距离d时，采用欧氏距离 测度），从而判断说话人是谁。语音识别系统结构框图如图1所示。图1语音识别系统结构框图2.1语者识别的概念语者识别就是根据说话人的语音信号来判别说话人的身份。语音是人的自然属性之一，由于说话人发音器官的生理差异以及后天形成的行

3、为差异，每个人的语音都带有强烈的个人色彩，这就使得通过分析语音信号来识别说话人成为可 能。用语音来鉴别说话人的身份有着许多独特的优点，如语音是人的固有的特征，不会丢失或遗忘；语音信号的采集方便，系统设备成本低；利用电话网络还可实 现远程客户服务等。因此，近几年来，说话人识别越来越多的受到人们的重视。 和其他生物识别技术如指纹识别、手形识别等相比较，说话人识别不仅使用方便， 而且属于非接触性，容易被用户接受，并且在已有的各种生物特征识别技术中， 是唯一可以用作远程验证的识别技术。因此，说话人识别的使用前景非常广泛： 今天，说话人识别技术已经关系到多学科的研究领域， 不同领域中的进步都对说 话人识

4、别的发展做出了贡献。说话人识别技术是集声学、语言学、计算机、信息 处理和人工智能等诸多领域的一项综合技术，使用需求将十分广阔。在吃力语音 信号的时候如何提取信号中关键的成分尤为重要。 语音信号的特征参数的好坏直 接导致了辨别的准确性。2.2特征参数的提取对于特征参数的选取，我们使用 mfcc的方法来提取。MFC（参数是基于人的 听觉特性利用人听觉的屏蔽效应，在 Mel标度频率域提取出来的倒谱特征参数。MFC（参数的提取过程如下：1. 对输入的语音信号进行分帧、加窗，然后作离散傅立叶变换，获得频谱 分布信息。设语音信号的DFT为：N 4j 2nkXa（k）» x（n）e P,0k EN

5、1（ 1）n =1其中式中x(n)为输入的语音信号，N表示傅立叶变换的点数。2. 再求频谱幅度的平方，得到能量谱。3. 将能量谱通过一组Mel尺度的三角形滤波器组。我们定义一个有M个滤波器的滤波器组(滤波器的个数和临界带的个数相近)， 采用的滤波器为三角滤波器，中心频率为f(m),m=1,2,3,，M本系统取M=10Q4. 计算每个滤波器组输出的对数能量。N丄2S(m) =ln(:JXa(k) | Hm(k),0 _ m - M -1其中Hm(k)为三角滤波器的频率响应。5. 经过离散弦变换(DCT得到MFC系数。M丄C(n) S(m)cos(二n(m -0.5/m),(3)m z00三n三N

6、 -1MFCC系数个数通常取2030,常常不用0阶倒谱系数，因为它反映的是频谱能 量，故在一般识别系统中，将称为能量系数，并不作为倒谱系数，本系统选取 20阶倒谱系数。2.3用矢量量化聚类法生成码本我们将每个待识的说话人看作是一个信源，用一个码本来表征。码本是从该 说话人的训练序列中提取的 MFC(特征矢量聚类而生成。只要训练的序列足够长， 可认为这个码本有效地包含了说话人的个人特征，而和讲话的内容无关。本系统采用基于分裂的LBG的算法设计VQ码本，Xk(k =1,2，,K)为训练序 列，B为码本。具体实现过程如下：1. 取提取出来的所有帧的特征矢量的型心(均值)作为第一个码字矢量B1。2.

7、将当前的码本Bn根据以下规则分裂，形成2m个码字。盅二時1 J(4)BBm(1 - ；)其中m从1变化到当前的码本的码字数，&是分裂时的参数，本文& =0.01。3. 根据得到的码本把所有的训练序列(特征矢量)进行分类，然后按照下面两个公式计算训练矢量量化失真量的总和Dn以及相对失真(n为迭代次数，初始n=0， D4=x，B为当前的码书)，若相对失真小于某一阈值£，迭代结束，当前 的码书就是设计好的2m个码字的码书，转5。否则，转下一步。量化失真量和：D八 min d(Xk,B)(5)k 4相对失真：Dn(6)4. 重新计算各个区域的新型心，得到新的码书，转3。5重复

8、2,3和4步，直到形成有M个码字的码书(M是所要求的码字数),其中 D0=1000O2.4 VQ的说话人识别设是未知的说话人的特征矢量Xi,|(,Xt，共有T帧是训练阶段形成的码书, 表示码书第m个码字，每一个码书有M个码字。再计算测试者的平均量化失真 D, 并设置一个阈值，若D小于此阈值，则是原训练者，反之则认为不是原训练者。D =1/V mind(Xj,Bm)l (7)j 生1 ：m -JM三、算法程序分析在具体的实现过程当中，采用了 matlab软件来帮助完成这个项目。在matlab 中主要由采集，分析，特征提取，比对几个重要部分。以下为在实际的操作中， 具体用到得函数关系和作用一一列举

9、在下面。3.1函数关系主要有两类函数文件Train.m和Test.m在Train.m 调用Vqlbg.m获取训练录音的vq码本，而Vqlbg.m调用mfcc.m 获取单个录音的mel倒谱系数，接着mfcc.m调用Melfb.m-将能量谱通过一组 Mel尺度的三角形滤波器组。在Test.m函数文件中调用Disteu.m计算训练录音(提供vq码本)和测试 录音(提供mfcc) mel倒谱系数的距离，即判断两声音是否为同一录音者提供。 Disteu.m 调用mfcc.m获取单个录音的 meI倒谱系数。mfcc.m调用Melfb.m- 将能量谱通过一组Mel尺度的三角形滤波器组。3.2具体代码说明3.

10、2.1 函数 mffc:fun ctio n r = mfcc(s, fs)m = 100;n = 256;沿-%方向取整对矩阵M赋值l = len gth(s);n bFrame = floor(l - n) / m) + 1; % for i = 1:nfor j = 1:n bFrameM(i, j) = s(j - 1) * m) + i); %endendh = hamming(n); % 加 hamming 窗，以增加音框左端和右端的连续性M2 = diag(h) * M;for i = 1:nbFrameframe(:,i) = fft(M2(:, i); % 对信号进行快速傅里

11、叶变换 FFTendt = n / 2;tmax = l / fs;m= melfb(20, n, fs); %将上述线性频谱通过 Mel 频率滤波器组得到 Mel 频 谱, 下面在将其转化成对数频谱n2 = 1 + floor(n / 2);z = m * abs(frame(1: n2, :).A2;r = dct(log(z); %将上述对数频谱，经过离散余弦变换(DCT)变换到倒谱域，即可得到Mel倒谱系数(MFCC参数)3.2.2 函数 disteu- 计算测试者和模板码本的距离function d = disteu(x, y)M, N = size(x); % 音频x赋值给【M N

12、】M2, P = size(y); % 音频 y 赋值给【 M2，P】if (M = M2)error(' 不匹配！ ') % 两个音频时间长度不相等endd = zeros(N, P);if (N < P)% 在两个音频时间长度相等的前提下copies = zeros(1,P);for n = 1:Nd(n,:) = sum(x(:, n+copies) - y) .A2, 1);endelsecopies = zeros(1,N);for p = 1:Pd(:,p) = sum(x - y(:, p+copies) .A2, 1)'end% 成对欧氏距离的两个

13、矩阵的列之间的距离endd = d.A0.5;3.2.3 函数 vqlbg- 该函数利用矢量量化提取了音频的 vq 码本function r = vqlbg(d,k)e = .01;r = mean(d, 2);dpr = 10000;for i = 1:log2(k)r = r*(1+e), r*(1-e);while (1 = 1)z = disteu(d, r);m,ind = min(z, , 2);t = 0;for j = 1:29r(:, j) = mean(d(:, find(ind = j), 2); x = disteu(d(:, find(ind = j), r(:, j

14、);for q = 1:length(x)t = t + x(q);endendif (dpr - t)/t) < e) break;elsedpr = t;endend end3.2.4 函数 test function finalmsg = test(testdir, n, code)for k = 1:n% read test sound file of each speakerfile = sprintf('%ss%d.wav', testdir, k);s, fs = wavread(file);v = mfcc(s, fs); % distmin = 4;%

15、d = disteu(v, code1); % 距离”dist = sum(min(d,2) / size(d,1); %变换得到一个距离的量测试阈值数量级得到测试人语音的 mel 倒谱系数 阈值设置处 就判断一次，因为模板里面只有一个文件 计算得到模板和要判断的声音之间的msgc = sprintf(' disp(msgc);% 此人匹配if dist <= distmin % msg = sprintf(' 求 !n', k);finalmsg = ' 定disp(msg);end% 此人不匹配if dist > distmin msg = sp

16、rintf(' 求 !n', k);finalmsg第%d位说话者和模板语音信号不匹配，不符合要和模板语音信号的差值为 :%10f ', dist);一个阈值，小于阈值，则就是这个人。第 %d 位说话者和模板语音信号匹配 , 符合要此位说话者符合要求 !' % 界面显示语句，可随意设= ' 此位说话者不符合要求 !'%界面显示语句，可随意设定disp(msg);endend3.2.5 函数 testDB这个函数实际上是对数据库一个查询 ,根据测试者的声音 ,找相应的文件 ,并且给 出是谁的提示function testmsg = testDB(t

17、estdir, n, code) nameList='1','2','3','4','5','6','7','8','9' ;% 这个是我们要识别的 9个数for k = 1:n% 数据库中每一个说话人的特征file = sprintf('%ss%d.wav', testdir, k);%找出文件的路径s, fs = wavread(file);v = mfcc(s, fs);%对找到的文件取 mfcc 变换distmin = inf;

18、k1 = 0;for l = 1:length(code)d = disteu(v, codel);dist = sum(min(d,2) / size(d,1);if dist < distmindistmin = dist;%这里和 test 函数里面一样 但多了一个具体语者的识别k1 = l;endendmsg=nameListk1 msgbox(msg);end3.2.6 函数 train- 该函数就是对音频进行训练，也就是提取特征参数 function code = train(traindir, n)k = 16;% number of centroids requiredf

19、or i = 1:n% 对数据库中的代码形成码本file = sprintf('%ss%d.wav', traindir, i); disp(file);计算 MFCC's 提取特征特征，返回值是Mel 倒谱系数，是一个 log 的 dct 得到的训练VQ码本通过矢量量化，得到原说s, fs = wavread(file); v = mfcc(s, fs);%codei = vqlbg(v, k); %话人的VQ码本end3.2.7 函数 melfb- 确定矩阵的滤波器 function m = melfb(p, n, fs) f0 = 700 / fs;fn2 = f

20、loor(n/2);lr = log(1 + 0.5/f0) / (p+1);% convert to fft bin numbers with 0 for DC termbl = n * (f0 * (exp(0 1 p p+1 * lr) - 1); 直接转换为 FFT 的数字模型 b1 = floor(bl(1) + 1;b2 = ceil(bl(2);b3 = floor(bl(3);b4 = min(fn2, ceil(bl(4) - 1;pf = log(1 + (b1:b4)/n/f0) / lr;fp = floor(pf);pm = pf - fp;r = fp(b2:b4)

21、 1+fp(1:b3);c = b2:b4 1:b3 + 1;v = 2 * 1-pm(b2:b4) pm(1:b3);m = sparse(r, c, v, p, 1+fn2) ；四、演示分析我们的功能分为两部分 : 对已经保存的 9 个数字的语音进行辨别和实时的 判断说话人说的是否为一个数 . 在前者的实验过程中 ,先把 9 个数字的声音保存 成wav的格式,放在一个文件夹中，作为一个检测的数据库然后对检测者实行识 别 , 系统给出提示是哪个数字 .在第二个功能中 , 实时的录取一段说话人的声音作为模板 , 提取 mfcc 特征 参数, 随后紧接着进行遇着识别 , 也就是让其他人再说相同的

22、话 , 看是否是原说话 者.实验过程及具体功能如下 ：先打开 Matlab 使 Current Directory 为录音及程序所所在的文件夹再打开文件“ enter .m”，点run运行，打开enter界面，点击“进入”按 钮进入系统。(注：文件包未封装完毕，目前只能通过此方式打开运行。 )(如下 图 figure1 )Q enterI 口 SpEaker Recognition SystemGroup 10退岀figurel在对数据库中已有的语者进行识别模块选择载入语音库语音个数；点击语音库录制模版进行已存语音信息的提取; 点击录音-test进行现场录音；点击语者判断进行判断数字，并显示出

23、来。 在实时语者识别模块：点击实时录制模板上的“录音-train ”按钮,是把新语者的声音以wav格 式存放在”实时模板”文件夹中，接着点击“实时录制模板”，把新的模板提取 特征值。随后点击实时语者识别模板上的“录音 -train ”按钮,是把语者的声音 以wav格式存放在”测试”文件夹中，再点击“实时语者识别”，在对测得的声音 提取特征值的同时，和实时模板进行比对，然后得出是否是实时模板中的语者 另外面板上的播放按钮都是播放相对应左边录取的声音。想要测量多次，只要接着录音，自动保存，然后程序比对音频就可以 退出只要点击菜单File/Exit ，退出程序。程序运行截图：（fig.2 ）运行后系

24、统界面I UntHIcdI E五、心得体会实验表明，该系统能较好地进行语音的识别，同时，基于矢量量化技术 （V Q）的语音识别系统具有分类准确，存储数据少，实时响应速度快等综合性能好 的特点.矢量量化技术在语音识别的使用方面，尤其是在孤立词语音识别系统中得到很好的使用，特别是有限状态矢量量化技术，对于语音识别更为有效。通过这次课程设计，我对语音识别有了更加形象化的认识，也强化了 MATLAB 的使用，对将来的学习奠定了基础。附： GUI 程序代码function pushbutton1_Callback(hObject, eventdata, handles)% hObject handle

25、to pushbutton1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)Channel_Str=get(handles.popupmenu3, 'String' );Channel_Number=str2double(Channel_Strget(handles.popupmenu3, 'Value' ) );global moodle

26、;moodle = train('模版' ,Channel_Number) %? y? oo ?o ?DDide ?士？% - Executes on button press in pushbutton2.function pushbutton2_Callback(hObject, eventdata, handles)% hObject handle to pushbutton2 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handglobal data1;glo

27、bal moodle ;test('测试' ,1,moodle) %ee±o ?o ? i 2a% function Open_Callback(hObject, eventdata, handles)% hObject handle to Open (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) filename,pathname=uiget

28、file('' )file=get(handles.edits,filename,pathname)y,f,b=wavread(file);% function Exit_Callback(hObject, eventdata, handles)% hObject handle to Exit (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) exit%

29、function About_Callback(hObject, eventdata, handles)% hObject handle to About (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)H= ' 语者识别 ' helpdlg(H, 'help text' )function File_Callback(hObject, ev

30、entdata, handles)% hObject handle to File (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% function Edit_Callback(hObject, eventdata, handles)% hObject handle to Edit (see GCBO)% eventdata reserved - to be defin

31、ed in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)%functionHelp_Callback(hObject, eventdata, handles)% hObject handle to Help (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUID

32、ATA)% - Executes on button press in pushbutton7.functionpushbutton7_Callback(hObject, eventdata, handles)% hObject handle to pushbutton7 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA) msg= ' 请速度录音 ?'msg

33、box(msg)clearglobal data1;%global dataDN1;AI = analoginput('winsound' );chan = addchannel(AI,1:2);duration = 3; %1 second acquisitionset(AI, 'SampleRate' ,8000)ActualRate = get(AI,'SampleRate');set(AI, 'SamplesPerTrigger' ,duration*ActualRate) set(AI, 'TriggerType

34、', 'Manual' )blocksize = get(AI,'SamplesPerTrigger' );Fs = ActualRate;start(AI) trigger(AI) data1,time,abstime,events = getdata(AI);fname=sprintf( 'E:Matlab 语音识别系统 实时模版 s1.wav' ) %dataDN1=wden(data1,'heursure','s','one',5,'sym8');denoise wa

35、vwrite(data1,fname)msgbox(fname)% - Executes on button press in pushbutton8.function pushbutton8_Callback(hObject, eventdata, handles)% hObject handle to pushbutton8 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDAT

36、A) global data1;%global dataDN1;sound(data1)%sound(dataDN1) axes(handles.axes1) %set to plot at axes1 plot(data1);%plot(dataDN1);xlabel( ' 训练采样序列 ' ),ylabel( ' 信号幅 ' );%xlabel('?2 e ? u Doa D'),ylabel('sym8D?2 ？卩?o 6 ?D?o? u ');grid on;clear% - Executes on button pres

37、s in pushbutton9.function pushbutton9_Callback(hObject, eventdata, handles)% hObject handle to pushbutton9 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA) msg= ' 请速度录音 ?'msgbox(msg)clearglobal data2;%glo

38、bal dataDN2;AI = analoginput('winsound' );chan = addchannel(AI,1:2); duration = 3;%1 second acquisitionset(AI, 'SampleRate' ,8000)ActualRate = get(AI,'SampleRate' );set(AI, 'SamplesPerTrigger',duration*ActualRate)set(AI, 'TriggerType' , 'Manual' ) bloc

39、ksize = get(AI,'SamplesPerTrigger');Fs = ActualRate; start(AI) trigger(AI)data2,time,abstime,events = getdata(AI);fname=sprintf( 'E:Matlab 语音识别系统 测试 s1.wav')%dataDN1=wden(data1,'heursure','s','one',5,'sym8');denoisewavwrite(data2,fname) msgbox(fname)%

40、- Executes on button press in pushbutton10. function pushbutton10_Callback(hObject, eventdata, handles)% hObject handle to pushbutton10 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA) global data2;%global dataDN

41、2;sound(data2)%sound(dataDN2)axes(handles.axes2) %set to plot at axes1 plot(data2);%plot(dataDN2);xlabel( ' 测试采样序列 ' ),ylabel( ' 信号幅 ' );%xlabel('2a e ?2 e ?u Dod D'),ylabel('sym8D?2 - ? ?o?D?o? u ');%grid on;clear% - Executes on button press in pushbutton11.function

42、pushbutton11_Callback(hObject, eventdata, handles)% hObject handle to pushbutton11 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA) global moodle ;testDB( ' 测试 ' ,1,moodle)% - Executes on button press in

43、pushbutton12.function pushbutton12_Callback(hObject, eventdata, handles)% hObject handle to pushbutton12 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)global moodle;moodle = train(' 实时模板 ' ,1)% - Execut

44、es on selection change in popupmenu3.function popupmenu3_Callback(hObject, eventdata, handles)% hObject handle to popupmenu3 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)as% Hints:contents = get(hObject,'String')returns popupmenu3 contentscell array% contentsget(hObject,'Value') returns selected item from popupmenu3str=get(handles.popupmenu3,'String' );val=str2num(strget(handles.popupmenu3,'Value' );switch valcase 1case 2case 3case