小波分析代码

复morlet小波分析，里面包含三个m文件

sst = sanHe;
%------------------------------------------------------ Computation
variance = std（sst）^2;
sst = （sst - mean（sst））/sqrt（variance） ;

n = length（sst）;
dt =1;
time = （1964:2011）;  
xlim = [19642011];  
pad = 1;      
dj = 0.1;   
s0 = 2*dt;    
j1 =fix（（log（n*dt/s0）/log（2））/dj）;  
%lag1 = 0.72;  
mother = ‘Morlet‘;

% Wavelet transform:
[waveperiodscalecoi] = wavelet（sstdtpaddjs0j1mother）;
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]）

Yticks = 2.^（fix（log2（min（period）））:fix（log2（max（period））））;
[ch]=contour（time（（period））realpart‘-‘）;  
clabel（ch）;
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）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1118  2013-01-05 09:08  小波分析\rain.m
     文件        2944  2012-12-28 10:37  小波分析\wave_bases.m
     文件        5460  2012-12-28 17:26  小波分析\wavelet.m
     目录           0  2013-01-08 17:40  小波分析\