mfcc.m

滤波器的设计以及零极点对应频率响应分析

%fid=fopen('F:\pcm_error\feature\feibiaozhun\F16A271.mfc','r');
%mfcc=fread(fid,[39,inf],'float');
fid=fopen('F:\Align\mono\feibiaozhun_mfc.scp.bak','r');
dr=fread(fid,'char');
di=char(dr);
l=length(di);
j=1;
i=1;
k=di(46);
k=1;
while i<=l
   
  if di(i) == di(46)
     dirm(k)={di(j:i,:)'};
     i=i+1;
     j=i;
     k=k+1;
       
  else
      i=i+1;
  end
end
fid1=fopen(dirm{1},'r');
mfcc=fread(fid,[39,inf],'float');