带宽均值7.txt

各种SVM分类算法

 clear all
for k=1:30
fname=sprintf('D:\\study\\音频采样\\speech\\%d.wav',k);
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)=sum(A(j,:));
 end
for j=1:row
   sum1(j)=0;
   for i=1:col
   sum1(j)=sum1(j)+i*A(j,i);
   end
end
centroid=sum1./energy;
for j=1:row
    sum2(j)=0;
    for i=1:col
        sum2(j)=sum2(j)+(i-centroid(j)).^2*A(j,i);
    end
end
 bandwidth=sum2./energy;
speechcentroidjuzhi(1,k)=mean(centroid);
speechbandwidthjuzhi(1,k)=mean(bandwidth);
end
speechcentroidjuzhi;
speechbw=speechbandwidthjuzhi'



for k=1:30
fname=sprintf('D:\\study\\音频采样\\music\\%d.wav',k);
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)=sum(A(j,:));
 end
for j=1:row
   sum1(j)=0;
   for i=1:col
   sum1(j)=sum1(j)+i*A(j,i);
   end
end
centroid=sum1./energy;
for j=1:row
    sum2(j)=0;
    for i=1:col
        sum2(j)=sum2(j)+(i-centroid(j)).^2*A(j,i);
    end
end
 bandwidth=sum2./energy;
musiccentroidjuzhi(1,k)=mean(centroid);
musicbandwidthjuzhi(1,k)=mean(bandwidth);
end
musiccentroidjuzhi;
musicbw=musicbandwidthjuzhi'



for k=1:30
fname=sprintf('D:\\study\\音频采样\\speechwithmusic\\%d.wav',k);
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)=sum(A(j,:));
 end
for j=1:row
   sum1(j)=0;
   for i=1:col
   sum1(j)=sum1(j)+i*A(j,i);
   end
end
centroid=sum1./energy;
for j=1:row
    sum2(j)=0;
    for i=1:col
        sum2(j)=sum2(j)+(i-centroid(j)).^2*A(j,i);
    end
end
 bandwidth=sum2./energy;
swmcentroidjuzhi(1,k)=mean(centroid);
swmbandwidthjuzhi(1,k)=mean(bandwidth);
end
swmcentroidjuzhi;
swmbw=swmbandwidthjuzhi'