单个音频熵值比较（同一个图形中显示）.m

各种SVM分类算法

clear
x=wavread('D:\study\师姐论文\邢峰师姐\邢峰论文程序\邢峰论文程序\audio files\music\3.wav');%读入原始信号
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)=sum(A(j,:));
 end
for j=1:row
   s(j)=0;
   s(j)=s(j)+energy(j).*log10(energy(j));
end
Err=-s;
plot(Err,':');hold on;xlabel('帧数');ylabel('熵m');
legend('music');
clear
x=wavread('D:\study\师姐论文\邢峰师姐\邢峰论文程序\邢峰论文程序\audio files\speech\3.wav');%读入原始信号
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)=sum(A(j,:));
 end
for j=1:row
   s(j)=0;
   s(j)=s(j)+energy(j).*log10(energy(j));
end
Err=-s;
plot(Err);hold on;xlabel('帧数');ylabel('熵s');
legend('music','speech');

gtext('err compare')