clear all
for k=1:20
  fname=sprintf('D:\\樣本\\speech\\%d.wav',k);
  x=wavread(fname);
  x=x(:,1);
  [c,l]=wavedec(x,3,'db4');
  ca=appcoef(c,l,'db4',3);%近似系數
  x=ca;
  enframex=enframe(x,128,64);%對信號x分幀處理，每幀有128個采樣點，相鄰幀之間有50%的疊加率
  [row,col]=size(enframex);%返回行，列，每一行代表一幀，每一列表示1到128個采樣點
  c=melcepst(x,22050);
  [row,col]=size(c);
  mfcc=zeros(1,12);
  for i=1:col
      for j=1:row  
        mfcc(1,i)=mean(c(j,i));
      end
  end
mfccspeech(k,:)=mfcc;
end
mfccspeech


for k=1:20
  fname=sprintf('D:\\樣本\\music\\%d.wav',k);
  x=wavread(fname);
  x=x(:,1);
  [c,l]=wavedec(x,3,'db4');
  ca=appcoef(c,l,'db4',3);%近似系數
  x=ca;
  enframex=enframe(x,128,64);%對信號x分幀處理，每幀有128個采樣點，相鄰幀之間有50%的疊加率
  [row,col]=size(enframex);%返回行，列，每一行代表一幀，每一列表示1到128個采樣點
  c=melcepst(x,22050);
  [row,col]=size(c);
  mfcc=zeros(1,12);
  for i=1:col
      for j=1:row  
        mfcc(1,i)=mean(c(j,i));
      end
  end
mfccmusic(k,:)=mfcc;
end
mfccmusic


for k=1:20
  fname=sprintf('D:\\樣本\\speech with music\\%d.wav',k);
  x=wavread(fname);
  x=x(:,1);
  [c,l]=wavedec(x,3,'db4');
  ca=appcoef(c,l,'db4',3);%近似系數
  x=ca;
  enframex=enframe(x,128,64);%對信號x分幀處理，每幀有128個采樣點，相鄰幀之間有50%的疊加率
  [row,col]=size(enframex);%返回行，列，每一行代表一幀，每一列表示1到128個采樣點
  c=melcepst(x,22050);
  [row,col]=size(c);
  mfcc=zeros(1,12);
  for i=1:col
      for j=1:row  
        mfcc(1,i)=mean(c(j,i));
      end
  end
mfccspeechwithmusic(k,:)=mfcc;
end
mfccspeechwithmusic