matlab语音识别系统(源代码)最新版

1、实用标准文案12.122.232.332.4VQ43.1.43.253.2.1mfcc53.2.2disteu53.2.3vqlbg.63.2.4test63.2.5testDB73.2.6train83.2.7melfb8.9.11文档大全附： GUI 程序代码 12一、设计任务及要求用 MATLAB 实现简单的语音识别功能；具体设计要求如下：用 MATLAB 实现简单的数字 19 的语音识别功能。二、语音识别的简单介绍基于 VQ 的说话人识别系统， 矢量量化起着双重作用。在训练阶段，把每一个说话者所提取的特征参数进行分类，产生不同码字所组成的码本。在识别(匹配)阶段，我们用VQ 方法计算平

2、均失真测度 (本系统在计算距离d 时，采用欧氏距离测度 )，从而判断说话人是谁。语音识别系统结构框图如图1 所示。1图 1 语音识别系统结构框图2.1 语者识别的概念语者识别就是根据说话人的语音信号来判别说话人的身份。语音是人的自然属性之一，由于说话人发音器官的生理差异以及后天形成的行为差异，每个人的语音都带有强烈的个人色彩，这就使得通过分析语音信号来识别说话人成为可能。用语音来鉴别说话人的身份有着许多独特的优点，如语音是人的固有的特征，不会丢失或遗忘； 语音信号的采集方便， 系统设备成本低； 利用电话网络还可实现远程客户服务等。因此，近几年来，说话人识别越来越多的受到人们的重视。与其他生物识

3、别技术如指纹识别、 手形识别等相比较，说话人识别不仅使用方便，而且属于非接触性，容易被用户接受，并且在已有的各种生物特征识别技术中，是唯一可以用作远程验证的识别技术。因此，说话人识别的应用前景非常广泛：今天，说话人识别技术已经关系到多学科的研究领域，不同领域中的进步都对说话人识别的发展做出了贡献。说话人识别技术是集声学、语言学、计算机、信息处理和人工智能等诸多领域的一项综合技术，应用需求将十分广阔。 在吃力语音信号的时候如何提取信号中关键的成分尤为重要。语音信号的特征参数的好坏直2接导致了辨别的准确性。2.2 特征参数的提取对于特征参数的选取，我们使用mfcc 的方法来提取。 MFCC 参数是

4、基于人的听觉特性利用人听觉的屏蔽效应，在Mel标度频率域提取出来的倒谱特征参数。MFCC 参数的提取过程如下：1. 对输入的语音信号进行分帧、 加窗，然后作离散傅立叶变换， 获得频谱分布信息。设语音信号的 DFT 为：N 1j 2 nkX a (k)x(n)e N ,0 k N 1（1 ）n 1其中式中 x(n) 为输入的语音信号， N 表示傅立叶变换的点数。2. 再求频谱幅度的平方，得到能量谱。3. 将能量谱通过一组 Mel 尺度的三角形滤波器组。我们定义一个有M 个滤波器的滤波器组（滤波器的个数和临界带的个数相近），采用的滤波器为三角滤波器，中心频率为f(m),m=1,2,3,，·

5、;M本系统取 M=100 。4. 计算每个滤波器组输出的对数能量。N 1S(m)ln(| Xa (k) |2 H m (k),（2）k 10 mM1其中 H m (k) 为三角滤波器的频率响应。5. 经过离散弦变换（ DCT）得到 MFCC 系数。3M1C(n)S (m)cos( n(m 0.5 / m),(3)m 00 n N1MFCC 系数个数通常取20 30 ，常常不用 0 阶倒谱系数，因为它反映的是频谱能量，故在一般识别系统中，将称为能量系数，并不作为倒谱系数，本系统选取 20 阶倒谱系数。2.3 用矢量量化聚类法生成码本我们将每个待识的说话人看作是一个信源，用一个码本来表征。码本是从

6、该说话人的训练序列中提取的MFCC 特征矢量聚类而生成。只要训练的序列足够长，可认为这个码本有效地包含了说话人的个人特征，而与讲话的内容无关。本系统采用基于分裂的LBG 的算法设计 VQ 码本， X k (k1,2, K ) 为训练序列， B 为码本。具体实现过程如下：1. 取提取出来的所有帧的特征矢量的型心 (均值 )作为第一个码字矢量 B1。2. 将当前的码本 Bm 根据以下规则分裂，形成 2m 个码字。 BmBm (1)（4 ）BmBm (1)其中 m 从 1变化到当前的码本的码字数，是分裂时的参数，本文 =0.01 。3. 根据得到的码本把所有的训练序列 (特征矢量 )进行分类，然后按

7、照下面两个公式计算训练矢量量化失真量的总和D n 以及相对失真 (n 为迭代次数，初始n=0 ， D 1 = ，B为当前的码书 )，若相对失真小于某一阈值，迭代结束，当前的码书就是设计好的 2m 个码字的码书，转。否则，转下一步。4量化失真量和：KD (n)min d ( X k , B)（5）k1相对失真：D (n 1)D n（6 ）|D n4. 重新计算各个区域的新型心，得到新的码书，转3。5. 重复 ， 和步，直到形成有 M 个码字的码书 (M 是所要求的码字数 )，其中 D0=10000 。2.4 VQ 的说话人识别设是未知的说话人的特征矢量 X1 ,K , XT ，共有T 帧是训练阶

8、段形成的码书，表示码书第 m 个码字，每一个码书有M 个码字。再计算测试者的平均量化失真 D ，并设置一个阈值，若D 小于此阈值，则是原训练者，反之则认为不是原训练者。D 1/ T min d(x ,B ) （ 7 ）j 1j m1 m M三、算法程序分析在具体的实现过程当中，采用了matlab软件来帮助完成这个项目。在matlab中主要由采集，分析，特征提取，比对几个重要部分。以下为在实际的操作中，具体用到得函数关系和作用一一列举在下面。3.1 函数关系主要有两类函数文件Train.m 和 Test.m5在 Train.m 调用 Vqlbg.m 获取训练录音的 vq 码本，而 Vqlbg.m

9、 调用mfcc.m获取单个录音的mel 倒谱系数，接着mfcc.m调用 Melfb.m-将能量谱通过一组 Mel 尺度的三角形滤波器组。在 Test.m 函数文件中调用 Disteu.m 计算训练录音 （提供 vq 码本）与测试录音（提供 mfcc ）mel 倒谱系数的距离， 即判断两声音是否为同一录音者提供。Disteu.m 调用 mfcc.m 获取单个录音的 mel 倒谱系数。mfcc.m 调用 Melfb.m-将能量谱通过一组Mel 尺度的三角形滤波器组。3.2 具体代码说明3.2.1 函数 mffc:function r = mfcc(s, fs)-m = 100;n = 256;l

10、= length(s);nbFrame = floor(l - n) / m) + 1;% 沿-方向取整for i = 1:nfor j = 1:nbFrameM(i, j) = s(j - 1) * m) + i);%对矩阵 M 赋值endendh = hamming(n);%加 hamming窗，以增加音框左端和右端的连续性6M2 = diag(h) * M;for i = 1:nbFrameframe(:,i) = fft(M2(:, i);%对信号进行快速傅里叶变换FFTendt = n / 2;tmax = l / fs;m = melfb(20, n, fs); %将上述线性频谱通

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

12、 ')%两个音频时间长度不相等endd = zeros(N, P);7if (N < P)% 在两个音频时间长度相等的前提下copies = zeros(1,P);for n = 1:Nd(n,:) = sum(x(:, n+copies) - y) .2, 1);endelsecopies = zeros(1,N);for p = 1:Pd(:,p) = sum(x - y(:, p+copies) .2, 1)'end% 成对欧氏距离的两个矩阵的列之间的距离endd = d.0.5;3.2.3 函数 vqlbg- 该函数利用矢量量化提取了音频的 vq 码本functi

13、on 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);8m,ind = min(z, , 2);t = 0;for j = 1:2ir(:, j) = mean(d(:, find(ind = j), 2);x = disteu(d(:, find(ind = j), r(:, j);for q = 1:length(x)t = t + x(q);endendif (dpr - t)/t) < e)bre

14、ak;elsedpr = t;endendend3.2.4 函数 testfunction 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);9v = mfcc(s, fs);% 得到测试人语音的mel 倒谱系数distmin = 4;% 阈值设置处% 就判断一次，因为模板里面只有一个文件d = disteu(v,code1);% 计算得到模板

15、和要判断的声音之间的“距离”dist = sum(min(d,2) / size(d,1);%变换得到一个距离的量%测试阈值数量级msgc = sprintf('与模板语音信号的差值为:%10f ', dist);disp(msgc);%此人匹配if dist <= distmin%一个阈值，小于阈值，则就是这个人。msg = sprintf('第 %d 位说话者与模板语音信号匹配,符合要求 !n',k);finalmsg = '此位说话者符合要求 !' % 界面显示语句，可随意设定disp(msg);end%此人不匹配if dist &g

16、t; distmin10msg=sprintf(' 第%d 位说话者与模板语音信号不匹配,不符合要求 !n', k);finalmsg= ' 此位说话者不符合要求 !' % 界面显示语句，可随意设定disp(msg);endend3.2.5 函数 testDB这个函数实际上是对数据库一个查询,根据测试者的声音,找相应的文件 ,并且给出是谁的提示function testmsg = testDB(testdir, n, code)nameList='1','2','3','4','5'

17、,'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;k1 = 0;11for l = 1:length(code)d = disteu(v, codel);dist = sum(min(d,2) / size(d,1);i

18、f dist < distmindistmin = dist;%这里和 test 函数里面一样但多了一个具体语者的识别k1 = l;endendmsg=nameListk1msgbox(msg);end3.2.6函数 train- 该函数就是对音频进行训练，也就是提取特征参数function code = train(traindir, n)k = 16;% number of centroids requiredfor i = 1:n% 对数据库中的代码形成码本file = sprintf('%ss%d.wav', traindir, i);disp(file);12s

19、, fs = wavread(file);v =mfcc(s,fs);% 计算 MFCC's 提取特征特征，返回值是Mel 倒谱系数，是一个log 的 dct 得到的codei= vqlbg(v,k);% 训练 VQ 码本通过矢量量化，得到原说话人的 VQ 码本end3.2.7函数 melfb- 确定矩阵的滤波器function m = melfb(p, n, fs)f0 = 700 / fs;fn2 = floor(n/2);lr = log(1 + 0.5/f0) / (p+1);% convert to fft bin numbers with 0 for DC term bl

20、= 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);13pm = pf - fp;r = fp(b2:b4) 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)

21、；四、演示分析我们的功能分为两部分 :对已经保存的 9 个数字的语音进行辨别和实时的判断说话人说的是否为一个数.在前者的实验过程中,先把 9 个数字的声音保存成wav 的格式 ,放在一个文件夹中 ,作为一个检测的数据库.然后对检测者实行识别,系统给出提示是哪个数字.在第二个功能中 ,实时的录取一段说话人的声音作为模板,提取 mfcc 特征参数 ,随后紧接着进行遇着识别,也就是让其他人再说相同的话,看是否是原说话者 .实验过程及具体功能如下：先打开 Matlab使 Current Directory为录音及程序所所在的文件夹再打开文件“ enter.m ”，点 run 运行，打开 enter 界

22、面，点击“进入”按钮进入系统。（注：文件包未封装完毕，目前只能通过此方式打开运行。）（如下图 figure1 ）14figure1在对数据库中已有的语者进行识别模块:15选择载入语音库语音个数；点击语音库录制模版进行已存语音信息的提取；点击录音 -test 进行现场录音；点击语者判断进行判断数字，并显示出来。在实时语者识别模块 :点击实时录制模板上的“录音-train ”按钮 ,是把新语者的声音以wav 格式存放在”实时模板”文件夹中, 接着点击“实时录制模板” ，把新的模板提取特征值。随后点击实时语者识别模板上的“录音-train ”按钮 ,是把语者的声音以wav 格式存放在”测试”文件夹中

23、,再点击“实时语者识别” ，在对测得的声音提取特征值的同时， 和实时模板进行比对， 然后得出是否是实时模板中的语者。另16外面板上的播放按钮都是播放相对应左边录取的声音。想要测量多次，只要接着录音，自动保存，然后程序比对音频就可以。退出只要点击菜单File/Exit ，退出程序。程序运行截图：（ fig.2 ）运行后系统界面五、心得体会17实验表明，该系统能较好地进行语音的识别，同时，基于矢量量化技术（）的语音识别系统具有分类准确，存储数据少， 实时响应速度快等综合性能好的特点矢量量化技术在语音识别的应用方面，尤其是在孤立词语音识别系统中得到很好的应用，特别是有限状态矢量量化技术，对于语音识别

24、更为有效。通过这次课程设计，我对语音识别有了更加形象化的认识，也强化了MATLAB 的应用，对将来的学习奠定了基础。附： GUI 程序代码functionvarargout = untitled2(varargin)% UNTITLED2 M-file for untitled2.fig% UNTITLED2, by itself, creates a new UNTITLED2 or raises the existing% singleton*.% H = UNTITLED2 returns the handle to a new UNTITLED2 or the handle to% t

25、he existing singleton*.18% UNTITLED2('CALLBACK',hObject,eventData,handles,.) calls the local% function named CALLBACK in UNTITLED2.M with the given input arguments.% UNTITLED2('Property','Value',.) creates a new UNTITLED2 or raises the%existing singleton*.Starting from the le

26、ft, property value pairs are%applied to the GUI before untitled2_OpeningFunction gets called.An% unrecognized property name or invalid value makes property application%stop.All inputs are passed to untitled2_OpeningFcn via varargin.%*See GUI Options on GUIDE's Tools menu.Choose "GUI allows

27、only one% instance to run (singleton)".% See also: GUIDE, GUIDATA, GUIHANDLES% Copyright 2002-2003 The MathWorks, Inc.% Edit the above text to modify the response to help untitled2% Last Modified by GUIDE v2.5 08-Jun-2010 23:58:57% Begin initialization code - DO NOT EDITgui_Singleton = 1;19gui_

28、State = struct('gui_Name','gui_Singleton''gui_OpeningFcn''gui_OutputFcn''gui_LayoutFcn''gui_Callback',mfilename,.,gui_Singleton,., untitled2_OpeningFcn,.,untitled2_OutputFcn,., , .);if nargin && ischar(varargin1)gui_State.gui_Callback = str2fun

29、c(varargin1);endif nargoutvarargout1:nargout = gui_mainfcn(gui_State, varargin:);elsegui_mainfcn(gui_State, varargin:);end% End initialization code - DO NOT EDIT% - Executes just before untitled2 is made visible.functionuntitled2_OpeningFcn(hObject, eventdata, handles, varargin)% This function has n

30、o output args, see OutputFcn.% hObjecthandle to figure% eventdatareserved - to be defined in a future version of MATLAB% handlesstructure with handles and user data (see GUIDATA)% varargincommand line arguments to untitled2 (see VARARGIN)20% Choose default command line output for untitled2 handles.o

31、utput = hObject;% Update handles structureguidata(hObject, handles);axes(findobj('tag' ,'axes13' );imshow('3.jpg' );axes(findobj('tag' ,'axes12' );imshow('1.jpg' );% UIWAIT makes untitled2 wait for user response (see UIRESUME)% uiwait(handles.figure1);

32、% - Outputs from this function are returned to the command line.functionvarargout = untitled2_OutputFcn(hObject, eventdata, handles)% varargoutcell array for returning output args (see VARARGOUT);% hObjecthandle to figure% eventdatareserved - to be defined in a future version of MATLAB% handlesstruc

33、ture with handles and user data (see GUIDATA)% Get default command line output from handles structure varargout1 = handles.output;% - Executes on button press in pushbutton1.functionpushbutton1_Callback(hObject, eventdata, handles)% hObjecthandle to pushbutton1 (see GCBO)% eventdatareserved - to be

34、defined in a future version of MATLAB21% handlesstructure with handles and user data (see GUIDATA)Channel_Str=get(handles.popupmenu3,'String');Channel_Number=str2double(Channel_Strget(handles.popupmenu3,'Value' );globalmoodle;moodle = train('模版 ' ,Channel_Number)%? y? ó&

35、#243;? ò ?DD ìáè? ±?% - Executes on button press in pushbutton2.functionpushbutton2_Callback(hObject, eventdata, handles)% hObjecthandle to pushbutton2 (see GCBO)% eventdatareserved - to be defined in a future version of MATLAB% handglobal data1; global moodle ;test( '测试

36、 ' ,1,moodle) %êê±ó? ò ? ì 2a% -functionOpen_Callback(hObject, eventdata, handles)% hObjecthandle to Open (see GCBO)% eventdatareserved - to be defined in a future version of MATLAB% handlesstructure with handles and user data (see GUIDATA)filename,pathname=uigetfil

37、e('' )file=get(handles.edits,filename,pathname)y,f,b=wavread(file);% -functionExit_Callback(hObject, eventdata, handles)22% hObjecthandle to Exit (see GCBO)% eventdatareserved - to be defined in a future version of MATLAB% handlesstructure with handles and user data (see GUIDATA)exit% -funct

38、ionAbout_Callback(hObject, eventdata, handles)% hObjecthandle to About (see GCBO)% eventdatareserved - to be defined in a future version of MATLAB% handlesstructure with handles and user data (see GUIDATA)H= ' 语者识别 'helpdlg(H,'help text' )% -functionFile_Callback(hObject, eventdata,

39、handles)% hObjecthandle to File (see GCBO)% eventdatareserved - to be defined in a future version of MATLAB% handlesstructure with handles and user data (see GUIDATA)% -functionEdit_Callback(hObject, eventdata, handles)% hObjecthandle to Edit (see GCBO)% eventdatareserved - to be defined in a future

40、 version of MATLAB% handlesstructure with handles and user data (see GUIDATA)% -23functionHelp_Callback(hObject, eventdata, handles)% hObjecthandle to Help (see GCBO)% eventdatareserved - to be defined in a future version of MATLAB% handlesstructure with handles and user data (see GUIDATA)% - Execut

41、es on button press in pushbutton7.functionpushbutton7_Callback(hObject, eventdata, handles)% hObjecthandle to pushbutton7 (see GCBO)% eventdatareserved - to be defined in a future version of MATLAB% handlesstructure with handles and user data (see GUIDATA)msg=' 请速度录音?'msgbox(msg)clearglobald

42、ata1;%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','Manual')block

43、size = get(AI,'SamplesPerTrigger');24Fs = 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');denoisewavwrite(data1,fname)msgbox(fname)% - Executes on button press in pushbutton8.functionpushbutton8_Callback(hObject, eventdata, handles)% hObjecthandle to pushbutton8 (see GCBO)% eventdatareserved - to be defined in a future version of MATLAB% handlesstructure with hand