matalab对0-9数字识别

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