三种音频文件的熵值比较（同一个框框显示）new.m

各种SVM分类算法

clear all
 for i=1:20
   fname=sprintf('D:\\study\\师姐论文\\邢峰师姐\\邢峰论文程序\\邢峰论文程序\\audio files\\speech with music\\%d.wav',i);
   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个采样点
   A=abs(enframex).^2;
   for j=1:row
     energy(j,1)=sum(A(j,:));
   end
   for j=1:row
       se=sum(abs(energy(j,1)).^2);
   end
   for j=1:row
       p(j,1)=abs(energy(j,1)).^2./se;
   end
   for j=1:row
     s(j,1)=0;
     s(j,1)=s(j,1)+p(j,1).*log2(p(j,1));
   end
   Err=-s;
   subplot(311);plot(Err,':r');hold on;xlabel('帧数');ylabel('熵swm');
   legend('speech with music');
  end

for i=1:20
   fname=sprintf('D:\\study\\师姐论文\\邢峰师姐\\邢峰论文程序\\邢峰论文程序\\audio files\\speech\\%d.wav',i);
   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个采样点
   A=abs(enframex).^2;
   for j=1:row
     energy(j,1)=sum(A(j,:));
   end
   for j=1:row
       se=sum(abs(energy(j,1)).^2);
   end
   for j=1:row
       p(j,1)=abs(energy(j,1)).^2./se;
   end
   for j=1:row
     s(j,1)=0;
     s(j,1)=s(j,1)+p(j,1).*log2(p(j,1));
   end
   Err=-s;
   subplot(312);plot(Err);hold on;xlabel('帧数');ylabel('熵speech');
   legend('speech');
  end

for i=1:20
   fname=sprintf('D:\\study\\师姐论文\\邢峰师姐\\邢峰论文程序\\邢峰论文程序\\audio files\\music\\%d.wav',i);
   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个采样点
   A=abs(enframex).^2;
   for j=1:row
     energy(j,1)=sum(A(j,:));
   end
   for j=1:row
       se=sum(abs(energy(j,1)).^2);
   end
   for j=1:row
       p(j,1)=abs(energy(j,1)).^2./se;
   end
   for j=1:row
     s(j,1)=0;
     s(j,1)=s(j,1)+p(j,1).*log2(p(j,1));
   end
   Err=-s;
   subplot(313);plot(Err,'--g');hold on;xlabel('帧数');ylabel('熵music');
   legend('music');
  end