【正文】 歸一化到[1,1]（2）短時能量的作用： 可以較好地區(qū)分出濁音和靜音。不過，低門限被超過未必就是語音的開始，有可能是時間很短的噪聲引起的。于是我們可以用信號的幅度作為特征來區(qū)分靜音和語音。由于采樣頻率的差異，幀長和幀移所對應的實際采樣點數(shù)都不一樣。但是這種算法并不可靠，因為人的發(fā)音有濁音和清音之分，濁音為聲帶振動發(fā)出，幅度高，周期性明顯，而清音則不會有聲帶的震動，只是靠空氣在口腔中的摩擦和沖擊，其短時能量一般比較小。語音和噪聲的主要區(qū)別在它們的能量上,如圖43所示。因此在許多場合會將En用下式來代替：這樣就不會因為取平方而造成信號的小取樣值的大取樣值出現(xiàn)較大差異。圖43(II)和圖44(II)分別為數(shù)字6數(shù)字4的短時過零率。在本系統(tǒng)中經多次試驗取定δ=。但是一些突發(fā)性的噪聲當然也會引起過零率或者短時能量的數(shù)值變高，但是都不能維持足夠長的時間，比如開關門窗、物體的碰撞等等引起的噪聲。圖45 端點檢測結果雙門限端點檢測顧名思義需要兩級檢測，即短時能量檢測和短時過零率檢測。在過度段中，由于參數(shù)數(shù)值較小，還不能確定是否真的進入語音段，只有兩個參數(shù)的其中一個超越了高門限才被認為是進入語音段。門限值的設置還沒有一個通用可靠的方法，需要根據(jù)經驗和特定環(huán)境進行調整。在整個語音識別的過程中，MFCC參數(shù)為了方便建模及計算，是按照幀計算的。 n=0,1,…,N/21 M是濾波器的數(shù)目，通常取24；一幀語音信號的點數(shù)為N，通常取256，這樣計算FFT的方便。MFCC特征系數(shù)的計算過程如下：（1）預處理：首先確定N為每一幀語音采樣點，如：N=256，并對每一幀序列s(n)進行分幀、預加重和加窗處理；（2）離散功率譜的計算：先對預處理的每一幀進行離散傅里葉變換得到其頻譜，再通過模的平方計算離散功率譜S(n)，功率譜如圖52所示。不管是在測試和建立參考模塊還是在識別階段，語音的起點和終點都必須用端點檢測來確定。假設測試和參考模板分別用T和R表示，為了比較它們之間的相似度，可以計算它們之間的距離D[T,R]，距離越小則相似度越高。若N=M則可以直接計算，否則要考慮將T（n）和R（m）對齊。通常，規(guī)整函數(shù)被限制在一個平行四邊形的網格內，如圖61所示。圖61 DTW原理圖由于在模板匹配過程中限定了彎折的斜率，因此平行四邊形之外的格點對應的幀匹配距離是不需要計算的。在x軸上的每一幀不再需要與y軸上的每一幀進行比較，而只是與y軸上[ymin，ymax]間的幀進行比較，ymin和ymax的計算公式為：ymin=x/2，0≤x≤xb，2x+(M2N)，xb x≤Nymax=2x，0≤x≤xa，x/2+(MN/2)，xa x≤N如果出現(xiàn)xa xb的情況，則彎折匹配的三段為（1，xb），（xb+1，xa），（xa+1，N）。接下來進行動態(tài)規(guī)劃，為每個格點(i,j)都計算其三個可能的前續(xù)格點的累積距離1D，2D和3D。通過點擊語音識別組件，MATLAB將會自動讀取對應文件夾中的wav格式的語音，然后進行前期的端點檢測和預處理，通過計算非線性預測系數(shù)（MFCC），并保留為參考庫及測試庫。它讓我們用戶自己定制MATLAB與用戶的交互方法，使命令窗口不是唯一與MATLAB的交互方式。圖71 GUI打開方式（1）組件面板在GUI界面中放置組件， GUI組件布置完成后，可以單擊菜單欄中最后一個三角形按鈕選項來觀察GUIDE的設計結果。待所有按鈕的功能的回調函數(shù)都編寫完畢后，點擊保存，然后點擊菜單欄中的三角形按鈕，則會彈出用戶開發(fā)界面，在該用戶開發(fā)界面中，用戶可以單擊按鈕以實現(xiàn)該按鈕所擁有的功能。用戶可以在點擊錄音后，錄制字、詞，錄制完畢后，可以點擊播放進行回放自己所錄制的語音，通過保存（uiputfile函數(shù)），用戶就可以保存數(shù)個錄制的語音。圖73運行結果圖本設計采用端點檢測算法，再經過特征提取，最后運用DTW算法來實現(xiàn)語音的識別。由以上數(shù)據(jù)分析可知識別率為86%。在本次設計中遇到的難題有很多。從選課題題目到設計的順利完成，陳杰老師都一直給予我耐心的指導。,hObject,eventData,handles,...) calls the local% function named CALLBACK in with the given input arguments.%% GUI(39。,...) creates a new GUI or raises the% existing singleton*. Starting from the left, property value pairs are% applied to the GUI before GUI_OpeningFunction gets called. An% unrecognized property name or invalid value makes property application% stop. All inputs are passed to GUI_OpeningF via varargin.%% *See GUI Options on GUIDE39。, mfilename, ... 39。, GUI_OpeningF, ... 39。, [] , ... 39。amp。end% End initialization code DO NOT EDIT% Executes just before GUI is made visible.function GUI_OpeningF(hObject, eventdata, handles, varargin)% This function has no output args, see OutputF.% hObject handle to figure% eventdata reserved to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% varargin mand line arguments to GUI (see VARARGIN)% Choose default mand line output for GUI = hObject。% hObject handle to figure% eventdata reserved to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Get default mand line output from handles structurevarargout{1} = 。 a=wavrecord(t*fs,fs)。pause()。axes()。[filename pathname]=uiputfile({39。39。},39。str=strcat(pathname,filename)。wavwrite(a,16000,str)。global x2。All MATLAB Files(*.wav)39。Allfiles(*.*)39。MultiSelect39。file=[L,fname]。[x1,x2]=vad(S,fs)。% Executes on button press in pushbutton11.function 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 b。 =b。*.wav39。*.*39。168。save(str,39。guidata(hObject, handles)。% Executes on button press in pushbutton13.function pushbutton13_Callback(hObject, eventdata, handles)% hObject handle to pushbutton13 (see GCBO)% eventdata reserved to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)global b。guidata(hObject, handles)。[fname,L]=uigetfile({39。39。},39。[S,fs]=wavread(file)。%guidata(hObject, handles)。global x2。for i=1:100fname=sprintf(39。 [x1 x2]=vad(x,fs)。 test(i).mfcc=m2。,i)。 m1=mfcc(x)。catch breakendenddist=zeros(a1,b1)。正在匹配……39。,i