小波分析源程序：内置强大小波工具箱的使用

kongfeng0824

1. load D:\Wavelet-analysis\小波程序\qiuxun.txt;
   
2. s=qiuxun(1:50);index=1961:2010;
   
3. t=1:60;
   
4. figure(1);
   
5. plot(index,s);
   
6. title('降水量距平过程');
   
7. 
   
8. figure(2);
   
9. w=cwt(s,[1:1:60],'morl');%小波变换（marr为墨西哥帽小波）
   
10. %plot(index,w);
    
11. contourf(index,t,w,-3000:500:-100,':')
    
12. hold on;
    
13. contourf(index,t,w,0:500:3000)
    
14. hold off;
    
15. ylabel('尺度a');xlabel('年份b');
    
16. title('降水量距平序列的小波变换实部的时频分布');
    
17. 
    
18. figure(3)
    
19. b=(abs(w)).^2;%小波系数摸平方
    
20. [c,h]=contourf(index,t,b,100);
    
21. ylabel('尺度a');xlabel('年份b');
    
22. title('小波系数模平方');
    
23. 
    
24. y=var(w);%小波方差
    
25. figure(4)
    
26. plot(y);title('降水量过程小波方差图');
    
27. ylabel('Var(a)');xlabel('a');
    
28. grid;
    
29. 
    
30. figure(5)
    
31. w1=cwt(s,13,'morl');
    
32. w2=cwt(s,22,'morl');
    
33. w3=cwt(s,56,'morl');
    
34. %w4=cwt(s,65,'morl');
    
35. %w5=cwt(s,78,'morl');
    
36. plot(index,w1,'g');
    
37. hold on;
    
38. plot(index,w2,'r');
    
39. hold on;
    
40. plot(index,w3,'b');
    
41. hold on;
    
42. %plot(index,w4,'r:');
    
43. hold on;
    
44. %plot(index,w5,'b:');
    
45. title('不同尺度下实部变化过程')
    
46. 
    
47. 
    
48. 
    
49. 
    
50. 
    
51. 
    
52. 附上测试数据供大家参考讨论。最近我也在不断地探索当中。由于是跨专业，对大气方面也是门外汉。很多都是刚刚学习的。请大家指点！！！！！
    

复制代码

尽头的尽头

不错~~~支持一个

zhouyuanshen

太棒了，谢谢楼主分享

zhouyuanshen

试了，很给力的说。现在学习怎么看程序作出的图

kongfeng0824

zhouyuanshen 发表于 2013-4-5 21:14

                               
登录/注册后可看大图

试了，很给力的说。现在学习怎么看程序作出的图

有了结果别忘了发帖子分享一下你的经验呵呵

zhu

谢谢分享！！！

zhu

无法下载（自动回复：请不要使用迅雷等下载工具，点我查看下载帮助）·············

zhu

1106272162@qq.com  给传一个呗，无法下载（自动回复：请不要使用迅雷等下载工具，点我查看下载帮助）

zhu

用Morlet做的话，程序怎么改？

kongfeng0824

zhu 发表于 2013-4-10 20:42

                               
登录/注册后可看大图

用Morlet做的话，程序怎么改？

程序中有一个选择调用函数的地方第9行