clear all
for k=1:20
 fname=sprintf('D:\\study\\師姐論文\\邢峰師姐\\邢峰論文程序\\邢峰論文程序\\audio files\\speech\\%d.wav',k);
 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個(gè)采樣點(diǎn)，相鄰幀之間有50%的疊加率
[row,col]=size(enframex);%返回行，列，每一行代表一幀，每一列表示1到128個(gè)采樣點(diǎn)
 A=abs(enframex).^2;
 energy=zeros(row,1);
 for j=1:row
   energy(j)=sum(A(j,:));
 end
 num(k)=0;
 for j=1:row
     if energy(j)<0.1
     num(k)=num(k)+1;
     end
 end
 speechjingyinbi(1,k)=num(k)/row;
end
 speechjingyinbi


for k=1:20
 fname=sprintf('D:\\study\\師姐論文\\邢峰師姐\\邢峰論文程序\\邢峰論文程序\\audio files\\music\\%d.wav',k);
 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個(gè)采樣點(diǎn)，相鄰幀之間有50%的疊加率
[row,col]=size(enframex);%返回行，列，每一行代表一幀，每一列表示1到128個(gè)采樣點(diǎn)
 A=abs(enframex).^2;
 energy=zeros(row,1);
 for j=1:row
   energy(j)=sum(A(j,:));
 end
 num(k)=0;
 for j=1:row
     if energy(j)<0.1
     num(k)=num(k)+1;
     end
 end
 musicjingyinbi(1,k)=num(k)/row;
end
 musicjingyinbi

for k=1:20
 fname=sprintf('D:\\study\\師姐論文\\邢峰師姐\\邢峰論文程序\\邢峰論文程序\\audio files\\speech with music\\%d.wav',k);
 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個(gè)采樣點(diǎn)，相鄰幀之間有50%的疊加率
[row,col]=size(enframex);%返回行，列，每一行代表一幀，每一列表示1到128個(gè)采樣點(diǎn)
 A=abs(enframex).^2;
 energy=zeros(row,1);
 for j=1:row
   energy(j)=sum(A(j,:));
 end
 num(k)=0;
 for j=1:row
     if energy(j)<0.1
     num(k)=num(k)+1;
     end
 end
 swmjingyinbi(1,k)=num(k)/row;
end
 swmjingyinbi