efv.m

来自「各种SVM分类算法」· M 代码 · 共 97 行

clear all
 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);%近似系数
   x=ca;
   y=enframe(x,128,64);%对信号x分帧处理，每帧有128个采样点，相邻帧之间有50%的叠加率
   [row,col]=size(y);%返回行，列，每一行代表一帧，每一列表示1到128个采样点
  %求短时能量
   A=abs(y).^2;
   energy=zeros(row,1);
   for j=1:row
     energy(j)=sum(A(j,:));
   end
  %求短时过零率
zcr=zeros(size(y,1),1);
delta=0.02;
for i=1:size(y,1)
    x=y(i,:);
    for j=1:length(x)-1
        if x(j)*x(j+1)<0 & abs(x(j)-x(j+1))>delta
            zcr(i)=zcr(i)+1;
        end
    end
end
EFV=zcr./energy;efv=EFV/abs(max(EFV));
musicEFV(k,1)=mean(efv);
end
musicEFV


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);%近似系数
   x=ca;
   y=enframe(x,128,64);%对信号x分帧处理，每帧有128个采样点，相邻帧之间有50%的叠加率
   [row,col]=size(y);%返回行，列，每一行代表一帧，每一列表示1到128个采样点
  %求短时能量
   A=abs(y).^2;
   energy=zeros(row,1);
   for j=1:row
     energy(j)=sum(A(j,:));
   end
  %求短时过零率
zcr=zeros(size(y,1),1);
delta=0.02;
for i=1:size(y,1)
    x=y(i,:);
    for j=1:length(x)-1
        if x(j)*x(j+1)<0 & abs(x(j)-x(j+1))>delta
            zcr(i)=zcr(i)+1;
        end
    end
end
EFV=zcr./energy;efv=EFV/abs(max(EFV));
speechEFV(k,1)=mean(efv);
end
speechEFV


clear all
 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);%近似系数
   x=ca;
   y=enframe(x,128,64);%对信号x分帧处理，每帧有128个采样点，相邻帧之间有50%的叠加率
   [row,col]=size(y);%返回行，列，每一行代表一帧，每一列表示1到128个采样点
  %求短时能量
   A=abs(y).^2;
   energy=zeros(row,1);
   for j=1:row
     energy(j)=sum(A(j,:));
   end
  %求短时过零率
zcr=zeros(size(y,1),1);
delta=0.02;
for i=1:size(y,1)
    x=y(i,:);
    for j=1:length(x)-1
        if x(j)*x(j+1)<0 & abs(x(j)-x(j+1))>delta
            zcr(i)=zcr(i)+1;
        end
    end
end
EFV=zcr./energy;efv=EFV/abs(max(EFV));
swmEFV(k,1)=mean(efv);
end
swmEFV