📄 err均值.m
字号:
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)=sum(A(j,:));
end
se=0;
for j=1:row
se=se+energy(j);
end
for j=1:row
p(j)=energy(j)./se;
end
s=0;
for j=1:row
s=s+p(j).*log2(p(j));
end
Err=-s;
swmErrzhuzhi(i,1)=mean(Err);
end
swmErrzhuzhi
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)=sum(A(j,:));
end
se=0;
for j=1:row
se=se+energy(j);
end
for j=1:row
p(j)=energy(j)./se;
end
s=0;
for j=1:row
s=s+p(j).*log2(p(j));
end
Err=-s;
speechErrzhuzhi(i,1)=mean(Err);
end
speechErrzhuzhi
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)=sum(A(j,:));
end
se=0;
for j=1:row
se=se+energy(j);
end
for j=1:row
p(j)=energy(j)./se;
end
s=0;
for j=1:row
s=s+p(j).*log2(p(j));
end
Err=-s;
musicErrzhuzhi(i,1)=mean(Err);
end
musicErrzhuzhi
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -