cp_extf.m

DSK6713 Speech.rar語音訊號處理應用範例程序

% ========================================================
% Filename : CP_EXTF.M
% ========================================================
% Speech Feature Extractor ( Using Cepstrum parameter )
% --------------------------------------------------------
% global variable : CBF;fft_point;fr_sz;win_sz;
%		    std_fr_num;ch_num
% --------------------------------------------------------
% call subprogram :
%   1.cbf.m	 : calculate critical bank Frequency response
%   2.cp.m	 : calculus cepstrum parameter of speech signal
% ========================================================
clear;
clear global;
global fft_point fr_sz win_sz std_fr_num order CBF
order=16;
word_num=30;
person_num=20;
serier_num=3;
sampling_rate=20000;
fft_point=2048;
fr_sz=200;
win_sz=400;
std_fr_num=35;
CBF=cbf(sampling_rate);
% ========================================================
%    Speech Feature Exactor (Getting CP Parameter)
% ========================================================
begin_time=clock;
for i=1:word_num
  processing_word_number=i
  sum=zeros(std_fr_num,order);
  for j=1:person_num
    for k=1:serier_num
      eval(['load ','c:\spdb\c',int2str(i),'_',int2str(j),'_',int2str(k),'.pcm']);
      eval(['CP_',int2str(i),'_',int2str(j),'_',int2str(k),'=cp(c',int2str(i),'_',int2str(j),'_',int2str(k),');']);
      filename=['c:\spdb\cp\c',int2str(i),'_',int2str(j),'_',int2str(k),'.cp'];
      fid=fopen(filename,'w');
      eval(['filedata=CP_',int2str(i),'_',int2str(j),'_',int2str(k),';']);
      fprintf(fid,'%g\n\r',filedata);
      fclose(fid);
      eval(['clear c',int2str(i),'_',int2str(j),'_',int2str(k)]);
      eval(['sum=sum+CP_',int2str(i),'_',int2str(j),'_',int2str(k),';']);
%      eval(['clear CP_',int2str(i),'_',int2str(j),'_',int2str(k)]);
    end
  end
  eval(['cp_ave',int2str(i),'=sum./(person_num*serier_num);']);
  filename=['c:\spdb\cp\cp_ave',int2str(i),'.cp'];
  fid=fopen(filename,'w');
  eval(['filedata=cp_ave',int2str(i),';']);
  fprintf(fid,'%g\n\r',filedata);
  fclose(fid);
  eval(['clear cp_ave',int2str(i)]);
end
clear sum filename fid filedata i j k processing_word_num CBF
RunTime=ceil(etime(clock,begin_time)/60);
% ========================================================
RunTime
pgmover
% ========================================================