我的一个朋友在用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