matlab语音识别系统

《matlab语音识别系统》由会员分享，可在线阅读，更多相关《matlab语音识别系统（11页珍藏版）》请在装配图网上搜索。

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

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

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

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

5、其中式中x(n)为输入的语音信号，N表示傅立叶变换的点数。2. 再求频谱幅度的平方，得到能量谱。3. 将能量谱通过一组Me 1尺度的三角形滤波器组。我们定义一个有M个滤波器的滤波器组(滤波器的个数和临界带的个数相近)， 采用的滤波器为三角滤波器，中心频率为f (m),m=l,2,3,，M 本系统取M=100。4. 计算每个滤波器组输出的对数能量。S(m)= 1n(艺 | Xa (k) |2 Hm(k),(2)k=10 mM1其中H (k)为三角滤波器的频率响应。m5. 经过离散弦变换(DCT)得到MFCC系数。C(n) = E S (m)cos5 n(m 0.5 / m),(3)m=00n N

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

7、的码本Bm根据以下规则分裂，形成2m个码字。严二 B （1+8） m m（4）B-二 B （1 -8）mm其中m从1变化到当前的码本的码字数，是分裂时的参数，本文=0.01。3. 根据得到的码本把所有的训练序列（特征矢量）进行分类，然后按照下面两个公式计算训练矢量量化失真量的总和Dn以及相对失真（n为迭代次数，初始 n=0, D-1 = -, B为当前的码书），若相对失真小于某一阈值，迭代结束，当前 的码书就是设计好的2m个码字的码书，转5。否则，转下一步。量化失真量和：D(n)= mmd(X ,B)(5)kk=1相对失真：| D(n-1) Dn |Dn4. 重新计算各个区域的新型心，得到新的

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

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

10、码说明3.2.1 函数 mffc:function r = mfcc(s, fs)m = 100;n = 256;l = 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 窗，以增加音框左端和右端的连续性M2 = diag(h) * M;for i = 1:nbFrameframe(:, i) = fft(M2(:, i);%对信号进行快速傅里叶变换F

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

12、ze(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+eopies) - 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-该函数利用矢量量

13、化提取了音频的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: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

14、(dpr - t)/t) e) break;else dpr = 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);v = 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 distminmsg = sprintf(第%小位说话者与模板语音信号不匹配，不符合要 求!n, k);finalmsg =此位说话者不符合要求!； 界面显示语句，可随意 设定disp(msg);endend3.2.5 函数 testDB 这个函数实际上是对数据库一个查询,根据测试者的声音,找相应的文件,并且给 出是谁的提示

16、function testmsg = testDB(testdir, 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 = distmin k1 = 0;mfcc(s, fs); inf;%对找到的文件取mfcc变换for l = 1:length(code)d = disteu(v, codel);dist = sum(min(d,2) / size(d

17、,1);if dist distmindis tmin = dis t;%这里和t es t函数里面一样 但多了一个具体语者的识别k1 = l;endendmsg=nameListk1 msgbox(msg);end3.2.6 函数 train-该函数就是对音频进行训练，也就是提取特征参数 function code = train（traindir, n）k = l6;for i = l:n% number of centroids required% 对数据库中的代码形成码本file = sprintf(%ss%d.wav, traindir, i); disp(file);s， fs =

18、 wavread(file);v = mfcc(s， fs);% 计算 MFCCs 提取特征特征，返回值是Mel倒谱系数，是一个log的det得到的codei = vqlbg(v， k); 话人的VQ码本end%训练VQ码本 通过矢量量化，得到原说四、演示分析我们的功能分为两部分:对已经保存的 9 个数字的语音进行辨别和实时的 判断说话人说的是否为一个数.在前者的实验过程中,先把 9 个数字的声音保存 成wav的格式，放在一个文件夹中，作为一个检测的数据库.然后对检测者实行识 别,系统给出提示是哪个数字.在第二个功能中，实时的录取一段说话人的声音作为模板，提取 mfcc 特征 参数，随后紧接着

19、进行遇着识别，也就是让其他人再说相同的话，看是否是原说话 者.实验过程及具体功能如下：先打开Ma tlab使Curren t Direc to ry为录音及程序所所在的文件夹再打开文件“enter.m”，点run运行，打开enter界面，点击“进入按 钮进入系统。（注：文件包未封装完毕，目前只能通过此方式打开运行。）（如下 图 figurel）进入退出enterSpeaker Recognition SystemGoup 10figure1在对数据库中已有的语者进行识别模块:选择载入语音库语音个数；点击语音库录制模版进行已存语音信息的提取点击录音-test进行现场录音；点击语者判断进行判断数字

20、，并显示出来。录b-train在实时语者识别模块:檣放-tr日in实时录制模板播放-馆吐实时语者识别实时谙者识别点击实时录制模板上的“录音-train按钮，是把新语者的声音以wav格 式存放在”实时模板”文件夹中, 接着点击“实时录制模板”，把新的模板提取 特征值。随后点击实时语者识别模板上的“录音-t rain ”按钮，是把语者的声音 以 wav 格式存放在”测试”文件夹中，再点击“实时语者识别”，在对测得的声音 提取特征值的同时，和实时模板进行比对，然后得出是否是实时模板中的语者。 另外面板上的播放按钮都是播放相对应左边录取的声音。想要测量多次，只要接着录音，自动保存，然后程序比对音频就可以。 退出只要点击菜单File/Exit,退出程序。