bpf_extf.m

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

% ========================================================
% Filename : BPF_EXTF.M
% ========================================================
% Speech Feature Extractor ( Using BPF parameter )
% --------------------------------------------------------
% global variable : fft_point;FR_BPFB;fr_sz;win_sz;
%		    std_fr_num;ch_num
% --------------------------------------------------------
% call subprogram :
%   1.frbpfb.m	 : calculate BPF bank Frequency response
%   2.bpf.m	 : calculus BPF parameter of speech signal
% ========================================================
clear;
clear global;
global fft_point fr_sz win_sz std_fr_num ch_num FR_BPFB
ch1_freq=150;
ch_num=15;
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;
FR_BPFB=frbpfb(ch1_freq,sampling_rate);
% ========================================================
%    Speech Feature Exactor (Getting BPF Parameter)
% ========================================================
begin_time=clock;
for i=1:word_num
  processing_word_number=i
  sum=zeros(std_fr_num,ch_num+1);
  for j=1:person_num
    for k=1:serier_num
      eval(['load ','c:\spdb\c',int2str(i),'_',int2str(j),'_',int2str(k),'.pcm']);
      eval(['BPF_',int2str(i),'_',int2str(j),'_',int2str(k),'=bpf(c',int2str(i),'_',int2str(j),'_',int2str(k),');']);
      filename=['c:\spdb\bpf35pn\c',int2str(i),'_',int2str(j),'_',int2str(k),'.bpf'];
      fid=fopen(filename,'w');
      eval(['filedata=BPF_',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+BPF_',int2str(i),'_',int2str(j),'_',int2str(k),';']);
      eval(['clear BPF_',int2str(i),'_',int2str(j),'_',int2str(k)]);
    end
  end
  eval(['ave',int2str(i),'=sum./(person_num*serier_num);']);
  filename=['c:\spdb\bpf35pn\ave',int2str(i),'.bpf'];
  fid=fopen(filename,'w');
  eval(['filedata=ave',int2str(i),';']);
  fprintf(fid,'%g\n\r',filedata);
  fclose(fid);
  eval(['clear ave',int2str(i)]);
end
clear sum filename fid filedata i j k processing_word_num FR_BPFB
RunTime=ceil(etime(clock,begin_time)/60);
% ========================================================
RunTime
pgmover
% ========================================================