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

matlab 语音识别系统（源代码）最新版目录一、设计任务及要求⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯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 程序代码⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯12、设计任务及要求用MATLAB 实现简单的语音识别功能；具体设计要求如下：用MATLAB 实现简单的数字1~9 的语音识别功能语音识别的简单介绍基于VQ的说话人识别系统，矢量量化起着双重作用。

在训练阶段，把每一个说话者所提取的特征参数进行分类，产生不同码字所组成的码本。

在识别（匹配）阶段，我们用VQ方法计算平均失真测度（本系统在计算距离d 时，采用欧氏距离测度），从而判断说话人是谁。

语音识别系统结构框图如图1 所示。

图 1 语音识别系统结构框图2.1 语者识别的概念语者识别就是根据说话人的语音信号来判别说话人的身份。

语音是人的自然属性之一，由于说话人发音器官的生理差异以及后天形成的行为差异，每个人的语音都带有强烈的个人色彩，这就使得通过分析语音信号来识别说话人成为可能。

用语音来鉴别说话人的身份有着许多独特的优点，如语音是人的固有的特征，不会丢失或遗忘；语音信号的采集方便，系统设备成本低；利用电话网络还可实现远程客户服务等。

因此，近几年来，说话人识别越来越多的受到人们的重视。

与其他生物识别技术如指纹识别、手形识别等相比较，说话人识别不仅使用方便，而且属于非接触性，容易被用户接受，并且在已有的各种生物特征识别技术中，是唯一可以用作远程验证的识别技术。

因此，说话人识别的应用前景非常广泛：今天，说话人识别技术已经关系到多学科的研究领域，不同领域中的进步都对说话人识别的发展做出了贡献。

说话人识别技术是集声学、语言学、计算机、信息处理和人工智能等诸多领域的一项综合技术，应用需求将十分广阔。

在吃力语音信号的时候如何提取信号中关键的成分尤为重要。

语音信号的特征参数的好坏直接导致了辨别的准确性。

2.2 特征参数的提取对于特征参数的选取，我们使用mfcc 的方法来提取。

MFCC参数是基于人的听觉特性利用人听觉的屏蔽效应，在Mel 标度频率域提取出来的倒谱特征参数。

MFCC参数的提取过程如下：1. 对输入的语音信号进行分帧、加窗，然后作离散傅立叶变换，获得频谱分布信息。

设语音信号的DFT为：N 1j 2 nkX a (k) x(n)e N,0 k N 1(1)n1其中式中x(n) 为输入的语音信号，N 表示傅立叶变换的点数。

2. 再求频谱幅度的平方，得到能量谱。

3. 将能量谱通过一组Mel 尺度的三角形滤波器组。

我们定义一个有M个滤波器的滤波器组(滤波器的个数和临界带的个数相近) ，采用的滤波器为三角滤波器，中心频率为f(m),m=1,2,3, ···，M 本系统取M=100。

4. 计算每个滤波器组输出的对数能量。

N12S(m) ln( |X a(k) |2H m (k)), (2)0 m M 1其中H m (k) 为三角滤波器的频率响应5. 经过离散弦变换( DCT)得到MFCC系数。

M1C(n) S (m)cos( n(m 0.5 / m)),(3) m00 n N 1MFCC系数个数通常取20—30，常常不用0 阶倒谱系数，因为它反映的是频谱能量，故在一般识别系统中，将称为能量系数，并不作为倒谱系数，本系统选取20 阶倒谱系数。

2.3 用矢量量化聚类法生成码本我们将每个待识的说话人看作是一个信源，用一个码本来表征。

码本是从该说话人的训练序列中提取的MFCC特征矢量聚类而生成。

只要训练的序列足够长，可认为这个码本有效地包含了说话人的个人特征，而与讲话的内容无关。

本系统采用基于分裂的LBG的算法设计VQ码本， X k(k 1,2, , K )为训练序列，B为码本。

具体实现过程如下：1. 取提取出来的所有帧的特征矢量的型心 （ 均值） 作为第一个码字矢量 B1。

2. 将当前的码本 Bm 根据以下规则分裂，形成 2m 个码字。

B m B m (1 ) B m B m (1 )其中 m 从 1 变化到当前的码本的码字数，ε是分裂时的参数，本文ε=0.01 。

3. 根据得到的码本把所有的训练序列 （ 特征矢量 ）进行分类，然后按照下面 两个公式计算训练矢量量化失真量的总和 D [ n]以及相对失真 （n 为迭代次数，初始=∞， B 为当前的码书 ） ，若相对失真小于某一阈值ε，迭代结束，当2m 个码字的码书，转５。

否则，转下一步。

相对失真：D(n 1) D n|D n |4. 重新计算各个区域的新型心，得到新的码书，转 3 。

5. 重复２ ，３ 和４步，直到形成有 M 个码字的码书 （M 是所要求的码字数 ） 其中 D0=10000。

2.4 VQ 的说话人识别设是未知的说话人的特征矢量 {X 1,K ,X T }，共有 T 帧是训练阶段形成的码 书，表示码书第 m 个码字，每一个码书有 M 个码字。

再计算测试者的平均量化 失真 D ，并设置一个阈值，若 D 小于此阈值，则是原训练者，反之则认为不是原 训练者。

D 1/T j 1min[d （x j ,B m ）]（7）j 11mM三、 算法程序分析在具体的实现过程当中， 采用了 matlab 软件来帮助完成这个项目。

在 matlab 中主要由采集，分析，特征提取，比对几个重要部分。

以下为在实际的操作中， 具体用到得函数关系和作用一一列举在下面。

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 获取单个录音的 mel 倒谱系数。

mfcc.m 调用 Melfb.m--- 将能量谱通过一n=0， D [ 1]前的码书就是设计好的 量化失真量和：Kmin d （X k ,B ）k1D (n)5）6）组Mel 尺度的三角形滤波器组。

3.2 具体代码说明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)); % 对信号进行快速傅里叶变换FFTendt = n / 2;tmax = l / fs;m = melfb(20, n, fs); %将上述线性频谱通过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(' 不匹配！') % 两个音频时间长度不相等endd = zeros(N, P);if (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%% 成对欧氏距离的两个矩阵的列之间的距离 end d = d.^0.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:2^ir(:, j) = mean(d(:, find(ind == j)), 2); x = disteu(d(:, find(ind == j)), r(:, j)); for q = 1:length(x) t = t + x(q); end endif (((dpr - t)/t) < e) break; else dpr = t; end end 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); 得到测试人语音的 mel 倒谱系数 阈值设置处 就判断一次，因为模板里面只有一个文件 计算得到模板和要判断的声音之间的v = mfcc(s, fs); % distmin = 4; %%d = disteu(v, code{1}); %距离”dist = sum(min(d,[],2)) / size(d,1); %% 测试阈值数量级 msgc = sprintf(' 与模板语音信号的差值为 :%10f ', dist); disp(msgc); % 此人匹配3.2.5 函数 testDB 这个函数实际上是对数据库一个查询 ,根据测试者的声音 , 找相应的文件 ,并且 给出是谁的提示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 = mfcc(s, fs); distmin = inf; k1 = 0;for l = 1:length(code) d = disteu(v, code{l}); dist = sum(min(d,[],2)) / size(d,1);变换得到一个距离的量if dist <= distmin % msg = sprintf(' 求 !\n', k);finalmsg = ' 定 disp(msg); end% 此人不匹配一个阈值，小于阈值，则就是这个人 第 %d 位说话者与模板语音信号匹配 此位说话者符合要求 !'; % 界面显示语句，, 符合要if dist > distmin msg = sprintf(' 求 !\n', k);finalmsg = ' 意设定disp(msg); end end第%d 位说话者与模板语音信号不匹配 , 不符合要此位说话者不符合要求 !'; % 界面显示语句，可随数据库中每一个说话人的特征%找出文件的路径 对找到的文件取 mfcc 变换if dist < distmindistmin = dist;%%者的识别k1 = l;end这里和test 函数里面一样但多了一个具体语endmsg=nameList{k1}msgbox(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);[s, fs] = wavread(file); v = mfcc(s, fs); %code{i} = vqlbg(v, k); % 话人的VQ码本end 3.2.7 函数melfb--- 确定矩阵的滤波器function m = melfb(p, n, fs)计算MFCC's 提取特征特征，返回值是Mel 倒谱系数，训练VQ码本是一个log 的dct 得到的通过矢量量化，得到原说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 = 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)1+fp(1:b3)]; c = [b2:b41:b3] + 1;v = 2 * [1-pm(b2:b4) pm(1:b3)]; m = sparse(r, c, v, p, 1+fn2)四、演示分析我们的功能分为两部分: 对已经保存的9 个数字的语音进行辨别和实时的判断说话人说的是否为一个数. 在前者的实验过程中, 先把9 个数字的声音保存成wav的格式, 放在一个文件夹中, 作为一个检测的数据库. 然后对检测者实行识别, 系统给出提示是哪个数字.在第二个功能中,实时的录取一段说话人的声音作为模板, 提取mfcc 特征参数, 随后紧接着进行遇着识别, 也就是让其他人再说相同的话, 看是否是原说话者.实验过程及具体功能如下：先打开Matlab 使Current Directory 为录音及程序所所在的文件夹再打开文件“ enter .m”，点run 运行，打开enter 界面，点击“进入” 按钮进入系统。