clear;
dim=zeros(20,30);              %20次試驗，30個小波子帶構成的差均值矩陣
for p=1:20                     %20次試驗，20對音頻例子比較
fname1=sprintf('G:\\亂\\music\\%d.wav',p);
x1=wavread(fname1);
y1=enframe(x1,512,256);
fname2=sprintf('G:\\亂\\speech with music\\%d.wav',p);
x2=wavread(fname2);
y2=enframe(x2,512,256);
z=min(size(y1,1),size(y2,1));   %音頻例子分幀后幀數的較小者
d=zeros(z,30);                  %一對音頻例子間的距離矩陣

  for i=1:z
    m=y1(i,:);
    t1=wpdec(m,4,'db4');          %db4小波4層分解
    cfs10a=wpcoef(t1,[1,0]); cfs11a=wpcoef(t1,[1,1]); cfs20a=wpcoef(t1,[2,0]); cfs21a=wpcoef(t1,[2,1]);
    cfs22a=wpcoef(t1,[2,2]); cfs23a=wpcoef(t1,[2,3]); cfs30a=wpcoef(t1,[3,0]); cfs31a=wpcoef(t1,[3,1]);
    cfs32a=wpcoef(t1,[3,2]); cfs33a=wpcoef(t1,[3,3]); cfs34a=wpcoef(t1,[3,4]); cfs35a=wpcoef(t1,[3,5]);
    cfs36a=wpcoef(t1,[3,6]); cfs37a=wpcoef(t1,[3,7]); cfs40a=wpcoef(t1,[4,0]); cfs41a=wpcoef(t1,[4,1]);
    cfs42a=wpcoef(t1,[4,2]); cfs43a=wpcoef(t1,[4,3]); cfs44a=wpcoef(t1,[4,4]); cfs45a=wpcoef(t1,[4,5]);
    cfs46a=wpcoef(t1,[4,6]); cfs47a=wpcoef(t1,[4,7]); cfs48a=wpcoef(t1,[4,8]); cfs49a=wpcoef(t1,[4,9]);
    cfs410a=wpcoef(t1,[4,10]); cfs411a=wpcoef(t1,[4,11]); cfs412a=wpcoef(t1,[4,12]); 
    cfs413a=wpcoef(t1,[4,13]); cfs414a=wpcoef(t1,[4,14]); cfs415a=wpcoef(t1,[4,15]);

    n=y2(i,:);
    t2=wpdec(n,4,'db4');          %db4小波4層分解
    cfs10b=wpcoef(t2,[1,0]); cfs11b=wpcoef(t2,[1,1]); cfs20b=wpcoef(t2,[2,0]); cfs21b=wpcoef(t2,[2,1]);
    cfs22b=wpcoef(t2,[2,2]); cfs23b=wpcoef(t2,[2,3]); cfs30b=wpcoef(t2,[3,0]); cfs31b=wpcoef(t2,[3,1]);
    cfs32b=wpcoef(t2,[3,2]); cfs33b=wpcoef(t2,[3,3]); cfs34b=wpcoef(t2,[3,4]); cfs35b=wpcoef(t2,[3,5]);
    cfs36b=wpcoef(t2,[3,6]); cfs37b=wpcoef(t2,[3,7]); cfs40b=wpcoef(t2,[4,0]); cfs41b=wpcoef(t2,[4,1]);
    cfs42b=wpcoef(t2,[4,2]); cfs43b=wpcoef(t2,[4,3]); cfs44b=wpcoef(t2,[4,4]); cfs45b=wpcoef(t2,[4,5]);
    cfs46b=wpcoef(t2,[4,6]); cfs47b=wpcoef(t2,[4,7]); cfs48b=wpcoef(t2,[4,8]); cfs49b=wpcoef(t2,[4,9]);
    cfs410b=wpcoef(t2,[4,10]); cfs411b=wpcoef(t2,[4,11]); cfs412b=wpcoef(t2,[4,12]); 
    cfs413b=wpcoef(t2,[4,13]); cfs414b=wpcoef(t2,[4,14]); cfs415b=wpcoef(t2,[4,15]);
    
    d10=(cfs10a-cfs10b).^2; d11=(cfs11a-cfs11b).^2; d20=(cfs20a-cfs20b).^2; d21=(cfs21a-cfs21b).^2;
    d22=(cfs22a-cfs22b).^2; d23=(cfs23a-cfs23b).^2; d30=(cfs30a-cfs30b).^2; d31=(cfs31a-cfs31b).^2;
    d32=(cfs32a-cfs32b).^2; d33=(cfs33a-cfs33b).^2; d34=(cfs34a-cfs34b).^2; d35=(cfs35a-cfs35b).^2;
    d36=(cfs36a-cfs36b).^2; d37=(cfs37a-cfs37b).^2; d40=(cfs40a-cfs40b).^2; d41=(cfs41a-cfs41b).^2;
    d42=(cfs42a-cfs42b).^2; d43=(cfs43a-cfs43b).^2; d44=(cfs44a-cfs44b).^2; d45=(cfs45a-cfs45b).^2;
    d46=(cfs46a-cfs46b).^2; d47=(cfs47a-cfs47b).^2; d48=(cfs48a-cfs48b).^2; d49=(cfs49a-cfs49b).^2;
    d410=(cfs410a-cfs410b).^2; d411=(cfs411a-cfs411b).^2; d412=(cfs412a-cfs412b).^2; 
    d413=(cfs413a-cfs413b).^2; d414=(cfs414a-cfs414b).^2; d415=(cfs415a-cfs415b).^2;
    
%     d(i,:)=[mean(d10),mean(d11),mean(20),mean(d21),mean(d22),mean(d23),mean(d30),mean(d31),mean(d32),mean(d33),mean(d34),
%             mean(d35),mean(d36),mean(d37),mean(d40),mean(d41),mean(d42),mean(d43),mean(d44),mean(d45),mean(d46),
%             mean(d47),mean(d48),mean(d49),mean(d410),mean(d411),mean(d412),mean(d413),mean(d414),mean(d415)];
    md10=mean(d10); md11=mean(d11); md20=mean(d20); md21=mean(d21); md22=mean(d22); md23=mean(d23); md30=mean(d30);
    md31=mean(d31); md32=mean(d32); md33=mean(d33); md34=mean(d34); md35=mean(d35); md36=mean(d36); md37=mean(d37); 
    md40=mean(d40); md41=mean(d41); md42=mean(d42); md43=mean(d43); md44=mean(d44); md45=mean(d45); md46=mean(d46);
    md47=mean(d47); md48=mean(d48); md49=mean(d49); md410=mean(d410); md411=mean(d411); md412=mean(d412); md413=mean(d413);
    md414=mean(d414); md415=mean(d415);
    dsp1=[md10,md11];dsp2=[md20,md21,md22,md23];dsp3=[md30,md31,md32,md33,md34,md35,md36,md37];
    dsp4=[md40,md41,md42,md43,md44,md45,md46,md47,md48,md49,md410,md411,md412,md413,md414,md415];
    d(i,:)=[dsp1,dsp2,dsp3,dsp4];         
  end
dmean=mean(d);
dim(p,:)=dmean;
end
dimention=mean(dim);