[分享]MFCC,即倒谱系数，也是语音的一个特征之一

happy

1. % MFCC,即倒谱系数，也是语音的一个特征之一，能够充分利用人耳的特性
   
2. % 具体的求法就是预处理，然后加窗，进行傅立叶变换，再求出功率普
   
3. % 然后得出其自然对数，最后进行dct变换
   
4. % 代码是用matlab写的，用到的一个语音工具箱，可以从网络上下载。
   
5. 
   
6. function mfcoef=mfcccoef(filename);
   
7. x=wavread(filename);
   
8. 
   
9. %归一化mel filter cof
   
10. bank=melbankm(24,256,8000,0,0.5,'m');
    
11. bank=full(bank);
    
12. bank=bank/max(bank(:));
    
13. 
    
14. 
    
15. % dct 系数
    
16. for k=1:12
    
17. n=0:23;
    
18. dctcoef(k,:)=cos((2*n+1)*k*pi/(2*24));
    
19. end
    
20. 
    
21. %归一化倒普窗口
    
22. w=1+6*sin(pi*[1:12]./12);
    
23. w=w/max(w);
    
24. 
    
25. % 一阶高通滤波器 1-0.9375z(-1)
    
26. 
    
27. xx=double(x);
    
28. xx=filter([1 -0.9375],1,xx);
    
29. 
    
30. % 分帧
    
31. y=enframe(xx,256,80);
    
32. [NumFrame,FrameLen]=size(y);
    
33. 
    
34. 
    
35. for i=1:NumFrame
    
36. y=xx(i,:);
    
37. s=y'.*hamming(256); %hamming 窗
    
38. t=abs(fft(s)); %傅立叶变换
    
39. t=t.^2; %功率普
    
40. c1=dctcoef*log(bank*t(1:129)+eps);
    
41. %c1=dct(log(bank*t(1:129)+eps));
    
42. c2=c1.*w';
    
43. m(i,:)=c2';
    
44. end
    
45. NumFrame
    
46. mfcoef=m;

复制代码

snort

为啥将3处用圆的图标遮住了

maoyaoyong

这两天在做一个问题，需要预处理，用的是加窗傅立叶变换，这方面我没学过，工作又催得紧，能否帮我编个程序，处理一下，以后我会慢慢学的。谢谢！下面是要处理的东西：
clear;
load f.mat;
rand('state',8);
A = rand(2,3)-0.5;
X=A*[f(1,:);f(2,:);f(3,:)];
for t=1:7000
    nor=norm(X(:,t));
    Y(:,t)=abs(X(:,t)/nor);
end
我解释一下，现在我要对Y进行分段傅立叶变换，每段长L=2048，每段乘上一个Hanning窗，连续短重叠的长度为d=614。
还要对Y竟进行相应的逆傅立叶变换，也帮我写出来吧！


谢谢！！！我9号就要用了，急求助！！！！
邮箱：maoyaoyong@163.com
谢谢！！！！！！

[ 本帖最后由 maoyaoyong 于 2006-11-5 14:39 编辑 ]

suffer

原帖由 maoyaoyong 于 2006-11-5 14:34 发表
这两天在做一个问题，需要预处理，用的是加窗傅立叶变换，这方面我没学过，工作又催得紧，能否帮我编个程序，处理一下，以后我会慢慢学的。谢谢！下面是要处理的东西：
clear;
load f.mat;
rand('state',8);
...

这个你可以参考 《MATLAB6.x信号处理》这本书，邹鲲写的
这本书你可以在本站ftp的/ebook/mathtools/Matlab/书籍及附盘/ 目录下找到

suffer

原帖由 snort 于 2006-7-17 13:42 发表

不太懂你的意思