急求！matlab高手们帮忙修改下小波分析程序

whmooncool

clear;
load nian.txt
sst = nian;
%------------------------------------------------------ Computation
variance = std(sst)^2;
sst = (sst - mean(sst))/sqrt(variance) ;
n = length(sst);
dt =1/1;
time = [0:length(sst)-1]*dt + 1960.0 ;  
xlim = [1960,2001];  
pad = 1;      
dj = 1/4;   
s0 = 2*dt;   
j1 =fix((log(n*dt/s0)/log(2))/dj);  
lag1 = 0.72;  
mother = 'Morlet';
% Wavelet transform:
[wave,period,scale,coi] = wavelet(sst,dt,pad,dj,s0,j1,mother);
power = (abs(wave)).^2 ;      
realpart=real(wave);
modulus=abs(wave);  
phase=atan2(imag(wave),real(wave));
%------------------------------------------------------ Plotting
% Plot wavelet coefficients of time series
subplot('position',[0.10 0.30 0.80 0.60])
levels = [1,2,4,8,16] ;
%level1 = [0,2,4,8];
v=[0,1,2,4];
Yticks = 2.^(fix(log2(min(period))):fix(log2(max(period))));
[c,h]=contour(time,log2(period),realpart,log2(levels),'-');  
%[c,h]=contour(time,log2(period),realpart,level1,'k-');
clabel(c,h,v);
xlabel('Time (year)')
ylabel('Period (years)')
%title('a) Morlet Wavelet Real Part of the Temperature Anomaly Series')
set(gca,'XLim',xlim(:))
set(gca,'YLim',log2([min(period),max(period)]), ...
    'YDir','default', ...
'YTick',log2(Yticks(:)), ...
'YTickLabel',Yticks)
hold on
levels = [0.0625,0.125,0.25,0.5,1] ;
%level1=[-8,-4,-2,0];
v=[-2,-4];
[c,h]=contour(time,log2(period),realpart,log2(levels),':');
%[c,h]=contour(time,log2(period),realpart,level1,'k:');
clabel(c,h,v);
hold off


我的数据是1960-2001的年降水资料，我是matlab的初学者，我找了个别人的小波分析的程序，但是不会修改。想做成如图图片。希望可以帮忙修改下

whmooncool

请问levels = [1,2,4,8,16] ;
levels = [0.0625,0.125,0.25,0.5,1] ;
这些数据都是怎么来的？

whmooncool

lag1 = 0.72;  
v=[-2,-4];还有power，realpart，modulus，phase有什么作用？

兰溪之水

@whmooncool建议你还是把小波分析的原理看看~

edwardli

兰溪之水 发表于 2012-2-27 14:30

                               
登录/注册后可看大图

@whmooncool建议你还是把小波分析的原理看看~

这个是Torrence的小波程序  首先建议你看他的A Practical Guide to Wavelet Analysis A_Practical_Guide_to_Wavelet_Analysis.pdf (1.06 MB, 下载次数: 139)

2012-3-2 11:21 上传

点击文件名下载附件 ，里面写的很清楚。
另外 个人建议是知道几个公式就行了
小波的形式：
小波系数的公式：

我用过的Morlet小波的形式：
关于本文献中的其他数据说明
只要你对号入座变量代表的含义就能改成你需要的语句了。
另外出图的语句，你可以参考matlab的Help或者百度。既然自己选择用这样的方法，那就基本掌握了它。对于以后也大有好处

毛真

A_Practical_Guide_to_Wavelet_Analysis.pdf,怎么是英文版的，看不懂呀！

陌小默的

学习了啊~~好厉害！

xiaoxiaoyu

我也过来学习，

njzqxt

一直想有这么段CODE

气象局

都说要看小波原理，可是好难哦，有朋友能指点一下到底怎么看啊？自己都没信心了！