bpf_vq.m

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

% ========================================================
% Filename : BPF_VQ.M
% ========================================================
% Clustering speech feature ( Using BPF parameter )
% --------------------------------------------------------
% main variable : std_fr_num;ch_num;knn_group_num
% --------------------------------------------------------
% call subprogram :
%   1.knn.m      : calculate 2-mean of speech feature
% ========================================================
clear;
clear global;
begin_time=clock;
ch_num=15;
order=ch_num+1;
word_num=30;
std_fr_num=35;
knn_group_order=6;
knn_group_num=2^(knn_group_order);
% ========================================================
% Load BPF Parameter
% ========================================================
filedata=zeros(std_fr_num,order);
for i=1:word_num
  filename=['c:\spdb\bpf\ave',int2str(i),'.bpf'];
  fid=fopen(filename,'r');
  filedata=fscanf(fid,'%g\n\r',[35 inf]);
  fclose(fid);
  V1((std_fr_num*(i-1)+1):(std_fr_num*i),:)=filedata;
end
BPF_VQ_CodeBook=zeros(knn_group_num,order);
for i=1:knn_group_order
  for j=1:2^(i-1)
    k=2^(i-1)+j-1;
    eval(['[V',int2str(2*k),' V',int2str(2*k+1),']=knn(V',int2str(k),');']);
  end
end
for i=1:knn_group_num
  eval(['BPF_VQ_CodeBook(',int2str(i),',:)=mean(V',int2str(knn_group_num+i-1),');']);
end
filename=['c:\spdb\bpf\bpf35p.vq'];
fid=fopen(filename,'w');
fprintf(fid,'%g\n\r',BPF_VQ_CodeBook);
fclose(fid);
clear i j k filedata filename fid
save BPF35pvq.mat BPF_VQ_CodeBook
RunTime=ceil(etime(clock,begin_time)/60);
% --------------------------------------------------------
figure(1);
mesh(BPF_VQ_CodeBook);
title('BPF Parameter Vector Quantize Codebook ( using KNN rule for the patterns" centres )');
xlabel('kth BPF Parameter');
ylabel('Group number');
zlabel('BPF Value');
axis([0 16 0 64 floor(min(min(BPF_VQ_CodeBook))) ceil(max(max(BPF_VQ_CodeBook)))]);
grid;
% --------------------------------------------------------
RunTime
%pgmover
% ========================================================