小波分析

心灵之约

借用论坛ccnasyq的帖子来做小波分析，但是我运行时却出现了错误提示：??? Undefined function or method 'wavelet' for input arguments of type 'double'.

请教下这个错误怎么才能解决呢？
代码如下，.q文件为： q.txt (372 Bytes, 下载次数: 53)

2012-11-27 19:02 上传

点击文件名下载附件
小波分析源数据

1. clear;
   
2. load q.txt
   
3. sst =q;
   
4. n=length(q);
   
5. variance = std(sst)^2;
   
6. sst = (sst - mean(sst))/sqrt(variance);
   
7. dt = 1;
   
8. year = [0:n-1]+ 1962.0 ;
   
9. xlim = [1962,2011];   %
   
10. pad = 1;      
    
11. dj = 1/12;     
    
12. s0 = 1/2*dt;   
    
13. j1 = 6.5/dj;   
    
14. lag1 = 0.72;  
    
15. mother = 'Morlet';
    
16. [wave,period,scale,coi] = wavelet(sst,dt,pad,dj,s0,j1,mother);
    
17. power = (abs(wave)).^2 ;  %计算小波系数的模的平方
    
18. modulus=abs(wave);  %计算小波系数的模
    
19. variance1=sum(power')/n;%计算小波方差
    
20. %画小波系数实部等值线图
    
21. subplot(3,1,1)
    
22. levels = [0,0.5,1.0,1.5,2.0,2.5];
    
23. v = [0,0.5,1.0,1.5];
    
24. Yticks = 0:5:45;
    
25. [c,h]=contour(year,period,real(wave),levels,'k-');
    
26. clabel(c,h,v,'fontsize',5);
    
27. xlabel('年份/year')
    
28. ylabel('周期/年 period/year')
    
29. title('(a)')
    
30. set(gca,'XLim',xlim(:))
    
31. set(gca,'YLim',[0 50], ...
    
32.      'YDir','default', ...
    
33. 'YTick',Yticks(:), ...
    
34. 'YTickLabel',Yticks)
    
35. hold on
    
36. levels = [-0.5,-1.0,-1.5,-2.0,-2.5];
    
37. v = [-0.5,-1.0,-1.5];
    
38. [c,h] = contour(year,period,real(wave),levels,'r--');
    
39. clabel(c,h,v,'fontsize',5);
    
40. 
    

复制代码

维欻?

wavelet有问题呀，wavelet在我的matlab里面没有这个函数

心灵之约

维欻? 发表于 2012-11-27 20:49

                               
登录/注册后可看大图

wavelet有问题呀，wavelet在我的matlab里面没有这个函数

那该如何解决这个问题呢？请赐教。

维欻?

建议你去请教一下ccnasyq，还有去http://www.docin.com/p-524412199.html看看wavelet 工具箱介绍，可能有你需要的

维欻?

心灵之约 发表于 2012-11-27 21:38

                               
登录/注册后可看大图

那该如何解决这个问题呢？请赐教。

我不懂小波分析，看了你的程序感觉这儿有问题

心灵之约

维欻? 发表于 2012-11-28 08:39

                               
登录/注册后可看大图

建议你去请教一下ccnasyq，还有去http://www.docin.com/p-524412199.html看看wavelet 工具箱介绍，可能有你 ...

真是太感谢你的帮助了哈，ccnasyq刚才回答了，说我少了这个函数，要下载才行，我刚试了好像能实现了，感谢。

心灵之约

维欻? 发表于 2012-11-28 08:40

                               
登录/注册后可看大图

我不懂小波分析，看了你的程序感觉这儿有问题

确实是有问题，缺少wavelet这个函数包。

心灵之约

这个问题已经解决了，如果论坛有朋友用到小波分析的话，可以问我哈，我只要知道就会尽力提供帮忙。

爱的侍者

心灵之约 发表于 2012-11-28 08:44

                               
登录/注册后可看大图

真是太感谢你的帮助了哈，ccnasyq刚才回答了，说我少了这个函数，要下载才行，我刚试了好像能实现了，感谢 ...

可否把工具箱共享~楼主 谢谢

心灵之约

爱的侍者 发表于 2012-12-3 17:44

                               
登录/注册后可看大图

可否把工具箱共享~楼主 谢谢

1. http://bbs.06climate.com/forum.php?mod=viewthread&tid=7519&page=1#pid69639

复制代码

请到这个帖子里下载。