求助各位！帮我看一下程序！

micky_hkw

1. clear;close all;
   
2. fs=1024;%采样频率
   
3. N=1024;%采样点数
   
4. dt=1/fs;%采样间隔
   
5. t=(0:N-1)*dt;
   
6. f=50;%基波频率
   
7. w=2*pi*f;
   
8. s1=sin(w*t);
   
9. subplot(4,1,1);plot(s1);
   
10. axis([0 100 -1 1]);
    
11. s2=1/3*sin(3*w*t);
    
12. subplot(4,1,2);plot(s2);
    
13. axis([0 100 -1 1]);
    
14. s3=1/5*sin(5*w*t);
    
15. subplot(4,1,3);plot(s3);
    
16. axis([0 100 -1 1]);
    
17. s4=1/7*sin(7*w*t);
    
18. subplot(4,1,4);plot(s4);
    
19. axis([0 100 -1 1]);
    
20. s=s1+s2+s3+s4;
    
21. figure(2);
    
22. plot(s);xlabel('时间 t/s');ylabel('幅值 A');
    
23. axis([0 100 -1 1]);
    
24. [c,l]=wavedec(s,5,'sym6');
    
25. a5=wrcoef('a',c,l,'sym6',5);
    
26. a4=wrcoef('a',c,l,'sym6',4);
    
27. a3=wrcoef('a',c,l,'sym6',3);
    
28. a2=wrcoef('a',c,l,'sym6',2);
    
29. a1=wrcoef('a',c,l,'sym6',1);
    
30. figure(3);
    
31. subplot(5,1,1);plot(a5,'LineWidth',1);ylabel('a5');
    
32. axis([0 100 -1 1]);
    
33. subplot(5,1,2);plot(a4,'LineWidth',1);ylabel('a4');
    
34. axis([0 100 -1 1]);
    
35. subplot(5,1,3);plot(a3,'LineWidth',1);ylabel('a3');
    
36. axis([0 100 -1 1]);
    
37. subplot(5,1,4);plot(a2,'LineWidth',1);ylabel('a2');
    
38. axis([0 100 -1 1]);
    
39. subplot(5,1,5);plot(a1,'LineWidth',1);ylabel('a1');
    
40. axis([0 100 -1 1]);
    
41. xlabel('时间 t/s');
    
42. d5=wrcoef('d',c,l,'sym6',5);
    
43. d4=wrcoef('d',c,l,'sym6',4);
    
44. d3=wrcoef('d',c,l,'sym6',3);
    
45. d2=wrcoef('d',c,l,'sym6',2);
    
46. d1=wrcoef('d',c,l,'sym6',1);
    
47. figure(4);
    
48. subplot(5,1,1);plot(d5,'LineWidth',1);ylabel('d5');
    
49. axis([0 100 -1 1]);
    
50. subplot(5,1,2);plot(d4,'LineWidth',1);ylabel('d4');
    
51. axis([0 100 -1 1]);
    
52. subplot(5,1,3);plot(d3,'LineWidth',1);ylabel('d3');
    
53. axis([0 100 -1 1]);
    
54. subplot(5,1,4);plot(d2,'LineWidth',1);ylabel('d2');
    
55. axis([0 100 -1 1]);
    
56. subplot(5,1,5);plot(d1,'LineWidth',1);ylabel('d1');
    
57. axis([0 100 -1 1]);
    
58. xlabel('时间 t/s');

复制代码

以上是小波分析，然后我想要把a1～5和d1～5（每一层也就是一个波形）用FFT画出频谱图，该怎么做？

happy

在命令窗口输入doc fft自己看帮助吧，里边有例子

然后再看一下帖子
http://forum.vibunion.com/thread-10344-1-1.html

估计就没什么问题了

micky_hkw

要是纵坐标要用幅值表示该如何做？

yangzj

对fft的结果取abs