clear all
 for i=1:20
   fname=sprintf('D:\\樣本\\speech with music\\%d.wav',i);
   x=wavread(fname);
   x=x(:,1);
   [c,l]=wavedec(x,3,'db4');
   ca=appcoef(c,l,'db4',3);%近似系數(shù)
   x=ca;
   enframex=enframe(x,128,64);%對信號x分幀處理，每幀有128個采樣點，相鄰幀之間有50%的疊加率
   [row,col]=size(enframex);%返回行，列，每一行代表一幀，每一列表示1到128個采樣點
  A=(abs(enframex).^2);
   energy=zeros(row,1);
   for j=1:row
     energy(j)=log10(sum(A(j,:)));
   end
   s=0;
     for j=1:row
     s=s+energy(j);
   end
   s=0.7*(s/row);
   swm(i)=0;
   for k=1:row
       if energy(k,1)<s
          swm(i)=swm(i)+1;
       end
   end
   swmshangbi(i,1)=swm(i)/row;
  end
 swmshangbi
  

 

clear
for i=1:20
   fname=sprintf('D:\\樣本\\speech\\%d.wav',i);
   x=wavread(fname);
   x=x(:,1);
   [c,l]=wavedec(x,3,'db4');
   ca=appcoef(c,l,'db4',3);%近似系數(shù)
   x=ca;
   enframex=enframe(x,128,64);%對信號x分幀處理，每幀有128個采樣點，相鄰幀之間有50%的疊加率
   [row,col]=size(enframex);%返回行，列，每一行代表一幀，每一列表示1到128個采樣點
   A=(abs(enframex).^2);
   energy=zeros(row,1);
   for j=1:row
     energy(j)=log10(sum(A(j,:)));
   end
   s=0;
     for j=1:row
     s=s+energy(j);
   end
   s=0.7*(s/row);
   speech(i)=0;
   for k=1:row
       if energy(k,1)<s
          speech(i)=speech(i)+1;
       end
   end
   speechshangbi(i,1)=speech(i)/row;
  end
 speechshangbi
 
   

clear
for i=1:20
   fname=sprintf('D:\\樣本\\music\\%d.wav',i);
   x=wavread(fname);
   x=x(:,1);
   [c,l]=wavedec(x,3,'db4');
   ca=appcoef(c,l,'db4',3);%近似系數(shù)
   x=ca;
   enframex=enframe(x,128,64);%對信號x分幀處理，每幀有128個采樣點，相鄰幀之間有50%的疊加率
   [row,col]=size(enframex);%返回行，列，每一行代表一幀，每一列表示1到128個采樣點
   A=(abs(enframex).^2);
   energy=zeros(row,1);
   for j=1:row
     energy(j)=log10(sum(A(j,:)));
   end
   s=0;
    for j=1:row
     s=s+energy(j);
   end
   s=0.7*(s/row);
   music(i)=0;
   for k=1:row
       if energy(k,1)<s
          music(i)=music(i)+1;
       end
   end
  musicshangbi(i,1)=music(i)/row;
  end
musicshangbi