请问各位大侠如何绘制小波的频谱（例如morlet小波）

快乐人生662

请给出源程序，不胜感激

c2019

关注，借贴同问

happy

本程序来自网络，版权归原作者所有，由于原作者无法考究，故此处没有标明
如原作者看到本贴，请和管理员联系添加相关说明或做其他处理

1. %%======================%%
   
2. %%小波包分析程序
   
3. %%=======================%%
   
4. clc;clear all;close all;
   
5. % %-----------------输入实际信号---------------
   
6. N=8192;
   
7. fs=20480;
   
8. [DATAfile DATApath]=uigetfile('.txt','输入信号');%显示一个取文件的窗口
   
9. FILENAME=[DATApath,DATAfile]
   
10. DATA=load(FILENAME);
    
11. 
    
12. Data=DATA(1:N);
    
13. DataLen=length(Data);%数据长度
    
14. t=[0:(DataLen-1)]'/fs;%时间数组
    
15. % -------------------小波包分解---------------
    
16. n=3;%分解层数
    
17. rr=wpdec(Data,n,'db15');%小波包分解
    
18. 
    
19. % ---------对特征频带波形进行分析-----
    
20. for i=1:2^n
    
21.      Datar(i,:)=wprcoef(rr,[n,i-1]);%求取每个频带的小波包系数
    
22. end
    
23. Datar=[Datar(1,:);Datar(2,:);Datar(4,:);Datar(3,:);Datar(7,:);Datar(8,:);Datar(6,:);Datar(5,:)];%数据调整
    
24. 
    
25. f=(0:N/2-1)*fs/N;
    
26. for i=1:2^n
    
27.     DataFFT(i,:)=abs(fft(Datar(i,:)))/(N/2);  %%fft
    
28. end
    
29. % %-----------小波包分解计算能量-------
    
30. for i=0:(2^n-1)
    
31.     rcfs=wprcoef(rr,[n i]);
    
32.     r(i+1,:)=rcfs;
    
33. end
    
34. r_check=[r(1,:);r(2,:);r(4,:);r(3,:);r(7,:);r(8,:);r(6,:);r(5,:)];
    
35. for i=1:2^n
    
36.     xx(i)=sum(r_check(i,:).^2)/N;%求能量
    
37. end
    
38. temp=sum(xx);
    
39. for i=1:2^n
    
40.     xx(i)=xx(i)/temp;%%能量归一化
    
41. end
    
42. %---------画小波分频带图---------
    
43. % figure()
    
44. % for i=1:2^n/2
    
45. %     subplot(4,1,i);
    
46. %     plot(t,Datar(i,:));ylabel(i);axis([min(t) max(t) min(Datar(i,:)) max(Datar(i,:))]);
    
47. %     if i==1
    
48. %        title('小波包分解')
    
49. %     end
    
50. % end
    
51. % xlabel('time');
    
52. % grid on
    
53. % figure()
    
54. % for i=2^n/2+1:2^n
    
55. %     subplot(4,1,i-2^n/2);
    
56. %     plot(t,Datar(i,:));ylabel(i);axis([min(t) max(t) min(Datar(i,:)) max(Datar(i,:))]);
    
57. %         if i== 2^n/2+1
    
58. %            title('小波包分解')
    
59. %         end
    
60. % end
    
61. % xlabel('time');
    
62. % grid on
    
63. %%-----------对每个频带求取频谱--------
    
64. 
    
65. figure()
    
66. set(gcf,'Name','ff1')
    
67. for i=1:2^n/2
    
68.     subplot(4,1,i);
    
69.     plot(f,DataFFT(i,1:N/2));ylabel(i);axis([0 10240 0 max(DataFFT(i,1:N/2))]);grid on
    
70.     if i==1
    
71.        title('小波包各频带频谱图')
    
72.     end
    
73. end
    
74. xlabel('frequency')
    
75. 
    
76. figure()
    
77. set(gcf,'Name','ff2')
    
78. for i=2^n/2+1:2^n
    
79.     subplot(4,1,i-2^n/2);
    
80.     plot(f,DataFFT(i,1:N/2));ylabel(i);axis([0 10240 0 max(DataFFT(i,1:N/2))]);grid on
    
81.     if i==2^n/2+1
    
82.        title('小波包各频带频谱图')
    
83.     end
    
84. end
    
85. xlabel('frequency')
    
86. 
    
87. %-------------画出能量谱-----------
    
88. figure()
    
89. set(gcf,'Name','能量谱')
    
90. width=0.5;
    
91. bar(xx,width);
    
92. colormap(cool);
    
93. title('小波包能量谱图');
    
94. xlabel('频 带');
    
95. ylabel('能量归一化');
    
96. axis([0 9 0 1]);
    
97. grid on

复制代码

快乐人生662

回复 happy 的帖子

万分感谢

快乐人生662

回复 happy 的帖子

你好能否再为我解答几个问题？

小波多分辨分析方面的问题
http://forum.vibunion.com/forum- ... fromuid-151591.html