【導(dǎo)讀】第一章語音知識(shí)基礎(chǔ)---------------------------------------------------------------6. 第一節(jié)語音識(shí)別的基本內(nèi)容-------------------------------------------6. 第二章HMM的理論基礎(chǔ)--------------------------------------------------------10. 第三章HMM算法實(shí)現(xiàn)的問題----------------------------------------------21. 第四章語音識(shí)別系統(tǒng)的設(shè)計(jì)---------------------------------------------------32. 第一節(jié)語音識(shí)別系統(tǒng)的開發(fā)環(huán)境-----------------------------------32. 計(jì)特征的動(dòng)態(tài)特性。做語音識(shí)別的識(shí)別率是很高的，其中用到了HMM的三個(gè)算法。對(duì)話的夢(mèng)想逐漸接近現(xiàn)實(shí)。人類的語言在人類的的智能組。意商是通過語言的方式有效的完成。作為人與人之間交流。算機(jī)交流的媒介。也是人機(jī)交互最重要的一步。之后，出現(xiàn)了諸多實(shí)用化的研究方向。

　　

【正文】 denom = denom + sum(tmp(:))。 57 end (l).weight(j) = nom/denom。 end fprintf(39。\n39。) end 五、 HMM 參數(shù)的初始化 在進(jìn)行 BaumWelch 迭代前，首先要對(duì) HMM 的參數(shù)進(jìn)行初始化。函數(shù) 就是完成這一工作的。初始化的過程中，首先將初始概率 賦值為只有第一個(gè)元素為 1,其他元素為 0 的數(shù)組，其長度為該 HMM 的狀態(tài)數(shù) N。接著將轉(zhuǎn)移概率賦值為如下形式的矩陣： 最后是對(duì)各狀態(tài)的混合高斯函數(shù)的均值、方差和權(quán)系數(shù)進(jìn)行初始化。這里采用了按狀態(tài)對(duì)觀察序列的參數(shù)平均分段，然后將所有觀察序 列中屬于一個(gè)段的參數(shù)組成一個(gè)大的矩陣，調(diào)用 58 Voice Box 工具箱中的聚類函數(shù) kmeans，由聚類結(jié)果，計(jì)算得到各 pdf 的均值、方差和權(quán)系數(shù)。 function hmm = inithmm(samples, M) K = length(samples)。 %語音樣本 N = length(M)。 %狀態(tài)數(shù) = N。 = M。 % 初始概率矩陣 = zeros(N,1)。 (1) = 1。 % 轉(zhuǎn)移概 率矩陣 =zeros(N,N)。 for i=1:N1 (i,i) = 。 (i,i+1) = 。 end (N,N) = 1。 % 概率密度函數(shù)的初始聚類 59 %平均分段 for k = 1:K T = size(samples(k).data,1)。 samples(k).segment=floor([1:T/N:T T+1])。 end %對(duì)屬于每個(gè)狀態(tài)的向量進(jìn)行 K 均值聚類，得 到連續(xù)混合正態(tài)分布 for i = 1:N %把相同聚類和相同狀態(tài)的向量組合到一個(gè)向量中 vector = []。 for k = 1:K seg1 = samples(k).segment(i)。 seg2 = samples(k).segment(i+1)1。 vector = [vector 。 samples(k).data(seg1:seg2,:)]。 end mix(i) = getmix(vector, M(i))。 end = mix。 function mix = getmix(vector, M) [mean esq nn] = kmeans(vector,M)。 60 % 計(jì)算每個(gè)聚類的標(biāo)準(zhǔn)差，對(duì)角陣，只保存對(duì)角線上的元素 for j = 1:M ind = find(j==nn)。 tmp = vector(ind,:)。 var(j,:) = std(tmp)。 end % 計(jì)算每個(gè)聚類中的元素?cái)?shù)，歸一化為各 pdf 的權(quán)重 weight = zeros(M,1)。 for j = 1:M weight(j) = sum(find(j==nn))。 end weight = weight/sum(weight)。 % 保存結(jié)果 = M。 = mean。 % M*SIZE = var.^2。 % M*SIZE = weight。 % M*1 六、訓(xùn)練主程序 上面的函數(shù) 只是實(shí)現(xiàn)了一次迭代，實(shí)際應(yīng)用中應(yīng)該是多次迭代才能得到 結(jié)果，同時(shí)還應(yīng)該給出一個(gè)結(jié)束迭代的 61 條件。本例中是計(jì)算所有觀察序列的輸出概率，對(duì)其進(jìn)行累加，得到總和輸出概率。當(dāng)此概率的相對(duì)變化小到一定數(shù)值的時(shí)候，結(jié)束迭代。另外，設(shè)定最大迭代次數(shù)為一個(gè)常數(shù)，當(dāng)?shù)螖?shù)超過此常數(shù)時(shí)，也停止迭代。下面是完整的訓(xùn)練程序 ： function [hmm, pout] = train(samples, M) %輸入： % samples – 樣本結(jié)構(gòu) % M 為每個(gè)狀態(tài)指定 pdf 個(gè)數(shù)，如： [3 3 3 3] %202。228。179。246。: % hmm 訓(xùn)練完成后的 hmm K = length(samples)。 % 計(jì)算語音參數(shù) disp(39。正在計(jì)算語音參數(shù) 39。)。 for k = 1:K if isfield(samples(k),39。data39。) amp。 ~isempty(samples(k).data) continue。 else samples(k).data = mfcc(samples(k).wave)。 end 62 end hmm = inithmm(samples, M)。 for loop = 1:40 fprintf(39。\n 第 %d 遍訓(xùn)練 \n\n39。,loop) hmm = baum(hmm, samples)。 %計(jì)算輸出總概率 pout(loop)=0。 for k = 1:K pout(loop) = pout(loop) + viterbi(hmm, samples(k).data)。 end fprintf(39。總和輸出概率 (log)=%d\n39。, pout(loop)) %比較兩個(gè) HMM 距離 if loop1 if abs((pout(loop)pout(loop1))/pout(loop)) 5e3 fprintf(39。收斂 !\n39。)。 return end end 63 end disp(39。迭代 40 次仍不收斂 , 退出 39。)。 程序中，輸入?yún)?shù)有兩個(gè)，結(jié)構(gòu)數(shù)組 samples 內(nèi)包含了觀察序列的信息，每個(gè) samples(k)都包含兩個(gè)成員 samples(k).wave和 samples(k).data，分別為該觀察序列的原始語音和參數(shù)。其中成員 data 可以在調(diào)用之前計(jì)算，也可以由 train 程序內(nèi)計(jì)算。數(shù)組 M 包含了各狀態(tài)對(duì)應(yīng)的高斯混合數(shù)。 七、測(cè)試程序 利用訓(xùn)練函數(shù) 和識(shí)別函數(shù) ，就可以對(duì)觀察序列進(jìn)行訓(xùn)練和識(shí)別了。下面是一段訓(xùn)練程序的腳本文件。 for i=1:11 sample=[]。 for k=1:10 file=sprintf(39。c:\\Program Files\\MATLAB71\\work\\cdhmm\\samples1\\%d\\%39。,i,k)。 disp(file)。 samplesz{i}{k}=wavread(file)。 sample(k).wave=samplesz{i}{k}。 sample(k).data=[]。 64 end hmm{i}=train(sample,[3 3 3 3])。 end save(39。hmm239。, 39。hmm39。)。 識(shí)別 結(jié)果 hmm 為一個(gè) cell 數(shù) 組，每個(gè)元素為一個(gè) hmm 結(jié)構(gòu)，下面一段識(shí)別程序。首先用函數(shù) wavread 將其讀入，再用函數(shù) vad 進(jìn)行端點(diǎn)檢測(cè)，計(jì)算出 MFCC 參數(shù)后，交由識(shí)別函數(shù) 計(jì)算得到其對(duì)數(shù)形式的輸出概率，最后用 max 函數(shù)找到識(shí)別結(jié)果。 下面是識(shí)別程序 的代碼： re=zeros(1,11)。 for i=1:11 for k=1:10 fname = sprintf(39。c:\\Program Files\\MATLAB71\\work\\cdhmm\\samples1\\%d\\%39。,i,k)。 % disp(fname)。 x = wavread(fname)。 [x1 x2] =Evad(x)。 m = mfcc(x)。 m = m(x12:x22,:)。 for j=1:11 pout(j)=viterbi(hmm{j}, m)。 end 65 [d,n] = max(pout)。 % fprintf(39。第 %d180。識(shí)別為類型 %d\n39。,i,n) fprintf(39。 詞 匯 類 別 ： 186。%d。 第 %d 識(shí) 別 為 類型 %d\n39。,i,k,n) if i==n re(i)=re(i)+1。 end end end rate=sum(re)/66。 fprintf(39。系統(tǒng)識(shí)別率： %d\n39。,rate)。 第 三節(jié) 實(shí)驗(yàn)結(jié)果 本語音識(shí)別系統(tǒng)所用的語音庫是十個(gè)足球術(shù)語的語音，它們分別為鏟球、沖撞、點(diǎn)球、防守、紅牌、加速、進(jìn)球、起動(dòng)、手球、暫停。 一、訓(xùn)練結(jié)果 第一次： 正在計(jì)算語音參數(shù) 第 1 遍訓(xùn)練 66 計(jì)算樣本參數(shù) ... 1 2 3 4 5 6 7 8 9 10 重估 轉(zhuǎn)移概率矩陣 A... 重估混合高斯的參數(shù) ... 1,1 1,2 1,3 2,1 2,2 2,3 3,1 3,2 3,3 4,1 4,2 4,3 總和輸出概率 (log)=+004 第 2 遍訓(xùn)練 計(jì)算樣本參數(shù) ... 1 2 3 4 5 6 7 8 9 10 重估轉(zhuǎn)移概率矩陣 A... 重估混合高斯的參數(shù) ... 1,1 1,2 1,3 2,1 2,2 2,3 3,1 3,2 3,3 4,1 4,2 4,3 總和輸出概率 (log)=+004 收斂 ! 67 第二次 ： 正在計(jì)算語音參數(shù) 第 1 遍訓(xùn)練 計(jì)算樣本參數(shù) ... 1 2 3 4 5 6 7 8 9 10 重估轉(zhuǎn)移概率矩陣 A... 重估混合高斯的參數(shù) ... 1,1 1,2 1,3 2,1 2,2 2,3 3,1 3,2 3,3 4,1 4,2 4,3 總和輸出概率 (log)=+004 第 2 遍訓(xùn)練 計(jì)算樣本參數(shù) ... 1 2 3 4 5 6 7 8 9 10 重估轉(zhuǎn)移概率矩陣 A... 重估混合高斯的參數(shù) ... 1,1 1,2 1,3 2,1 2,2 2,3 68 3,1 3,2 3,3 4,1 4,2 4,3 總和輸出概率 (log)=+004 收斂 ! 第三次： 正在計(jì)算語音參數(shù) 第 1 遍訓(xùn)練 計(jì)算樣本參數(shù) ... 1 2 3 4 5 6 7 8 9 10 重估轉(zhuǎn)移概率矩陣 A... 重估混合高斯的參數(shù) ... 1,1 1,2 1,3 2,1 2,2 2,3 3,1 3,2 3,3 4,1 4,2 4,3 總和輸出概率 (log)=+004 第 2 遍訓(xùn)練 計(jì)算樣本參數(shù) ... 1 2 3 4 5 6 7 8 9 10 69 重估轉(zhuǎn)移概率矩陣 A... 重估混合高 斯的參數(shù) ... 1,1 1,2 1,3 2,1 2,2 2,3 3,1 3,2 3,3 4,1 4,2 4,3 總和輸出概率 (log)=+004 收斂 ! 第四次： 正在計(jì)算語音參數(shù) 第 1 遍訓(xùn)練 計(jì)算樣本參數(shù) ... 1 2 3 4 5 6 7 8 9 10 重估轉(zhuǎn)移概率矩陣 A... 重估混合高斯的參數(shù) ... 1,1 1,2 1,3 2,1 2,2 2,3 3,1 3,2 3,3 4,1 4,2 4,3 總和輸出概率 (log)=+004 70 第 2 遍訓(xùn)練 計(jì)算 樣本參數(shù) ... 1 2 3 4 5 6 7 8 9 10 重估轉(zhuǎn)移概率矩陣 A... 重估混合高斯的參數(shù) ... 1,1 1,2 1,3 2,1 2,2 2,3 3,1 3,2 3,3 4,1 4,2 4,3 總和輸出概率 (log)=+004 收斂 ! 第五次： 正在計(jì)算語音參數(shù) 第 1 遍訓(xùn)練 重估轉(zhuǎn)移概率矩陣 A... 重估混合高斯的參數(shù) ... 1,1 1,2 1,3 2,1 2,2 2,3 3,1 3,2 3,3 71 4,1 4,2 4,3 總和輸出概率 (log)=+004 第 2 遍訓(xùn)練 計(jì)算樣本參數(shù) ... 1 2 3 4 5 6 7 8 9 10 重估轉(zhuǎn)移概率矩陣 A... 重估混合高斯的參數(shù)