function?[num?sbjg]=shibiesb（file）
%識別原始的語音文件
%filename=‘D:\work\語音\1107\語音文件\t222‘;
filename=file;
[qszznum]=readlab（[filename?‘.lab‘]）;
[dfpdttct]=readhtk（[filename?‘.wav‘]）;
pd=zeros（43000）;%定義一個臨時的矩陣pd用于存儲4個語音波形序列序列長度不超過3000
%通過.lab和.wav文件將語音文件中的聲音信號讀入pd中
for?i=1:4
????qs1（i）=round（qs（i）/fp/10000000）;
????zz1（i）=round（zz（i）/fp/10000000）;
????cb（i）=zz1（i）-qs1（i）+1;
????pd（i1）=cb（i）;
????for?j1=qs1（i）:zz1（i）
????????pd（ij1-qs1（i）+2）=d（j1）;
????end????????????????
end
%分別對4個數字進行識別??????
for?sb=1:4
????[cepsfreqrespfbfbreconfreqrecon]?=?	mfcc（pd（sb2:pd（sb1））?8000?160）;?%計算每個數字的mfcc參數???
????[a1?a2]=size（ceps）;
????if?（a2>14）
????????datacell（::sb）=ceps（:1:15）;%截取前13*15部分
????????
????%讀取0到9十個數字的HMM特征參數并計算概率
????for?i=1:10
????????fid?