主题:代朋友提个问题 挺急的 多谢
我的一个朋友在用matlab做毕设,不大顺利,我不懂这个
看她挺着急,我来代她提个问题
她做的题目是 基于小波包的语音特征提取
她做的程序在附件里
其他的,我也不是很懂
希望高手帮忙看看
用记事本打开是这样的:
% 归一化mel滤波器组系数
bank=melbankm(24,256,8000,0,0.5,'m');
bank=full(bank);
bank=bank/max(bank(:));
% DCT系数,12*24
for k=1:12
n=0:23;
dctcoef(k,:)=cos((2*n+1)*k*pi/(2*24));
end
% 归一化倒谱提升窗口
w = 1 + 6 * sin(pi * (1:12) ./ 12);
w = w/max(w);
% the current sound x in the voicebox
% load signal
x=wavread('63.wav');
% 预加重滤波器
xx=double(x);
xx=filter([1 -0.9375],1,xx);
% 语音信号分帧
xx=enframe(xx,256,80);
% 小波包分解得WPXkm
% decompoese xx at depth6 with db2 wavelet packets
wpt=wpdec(xx,6,'db2');
% read packet (a,b) coefficient
for a=1:6
for b=2^(a-1)
cfsab=wpcoef(wpt,[a,b]);
end
end
% 计算对数能量
% the number of the coefficients of packet(a,b)
for k=1:24
for m=1:size(cfsab,1)
cfskm=cfsab;
WPXkm=cfskm;
Sk=log(cumsum((WPXkm)^2/size(cfsab,1)));
end
end
% 经过DCT变换得到新参数
for i=1:24
Ci=cumsum(S(k).cos(pi*i*(k-1/2)/24));
end
看她挺着急,我来代她提个问题
她做的题目是 基于小波包的语音特征提取
她做的程序在附件里
其他的,我也不是很懂
希望高手帮忙看看
用记事本打开是这样的:
% 归一化mel滤波器组系数
bank=melbankm(24,256,8000,0,0.5,'m');
bank=full(bank);
bank=bank/max(bank(:));
% DCT系数,12*24
for k=1:12
n=0:23;
dctcoef(k,:)=cos((2*n+1)*k*pi/(2*24));
end
% 归一化倒谱提升窗口
w = 1 + 6 * sin(pi * (1:12) ./ 12);
w = w/max(w);
% the current sound x in the voicebox
% load signal
x=wavread('63.wav');
% 预加重滤波器
xx=double(x);
xx=filter([1 -0.9375],1,xx);
% 语音信号分帧
xx=enframe(xx,256,80);
% 小波包分解得WPXkm
% decompoese xx at depth6 with db2 wavelet packets
wpt=wpdec(xx,6,'db2');
% read packet (a,b) coefficient
for a=1:6
for b=2^(a-1)
cfsab=wpcoef(wpt,[a,b]);
end
end
% 计算对数能量
% the number of the coefficients of packet(a,b)
for k=1:24
for m=1:size(cfsab,1)
cfskm=cfsab;
WPXkm=cfskm;
Sk=log(cumsum((WPXkm)^2/size(cfsab,1)));
end
end
% 经过DCT变换得到新参数
for i=1:24
Ci=cumsum(S(k).cos(pi*i*(k-1/2)/24));
end