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??小波分析\