关于水文小波分析，请大家翻译，帮忙一下，谢谢

wendeping

clear
sst=[-3.246319444
-3.066319444
-0.786319444
-0.296319444
-0.536319444
2.963680556
5.053680556
3.223680556
3.933680556
3.113680556
1.823680556
-1.616319444
-2.386319444
-1.456319444
-1.556319444
0.233680556
2.293680556
2.323680556
3.343680556
1.173680556
2.543680556
1.153680556
-0.836319444
-2.746319444
-3.546319444
-3.426319444
-1.966319444
-0.096319444
4.783680556
4.353680556
6.073680556
6.063680556
4.433680556
-0.056319444
0.413680556
-0.726319444
-1.526319444
-1.296319444
-0.086319444
0.423680556
2.873680556
2.793680556
5.943680556
3.213680556
3.413680556
1.723680556
-2.886319444
-4.146319444
-4.996319444
-5.366319444
-3.266319444
-2.076319444
1.303680556
2.573680556
3.913680556
3.743680556
3.203680556
0.633680556
-2.186319444
-3.506319444
-3.926319444
-1.466319444
-1.496319444
-1.286319444
1.253680556
4.563680556
3.433680556
3.763680556
4.773680556
1.583680556
-0.376319444
-3.356319444
-3.876319444
-4.056319444
-1.026319444
-0.046319444
1.513680556
3.323680556
3.023680556
0.543680556
-1.926319444
-3.356319444
-3.596319444
-3.626319444
-4.586319444
-4.366319444
-2.046319444
-2.246319444
-1.666319444
1.243680556
3.883680556
4.913680556
3.643680556
-0.176319444
-3.576319444
-5.106319444
-4.926319444
-4.146319444
-4.036319444
-0.366319444
-0.516319444
2.003680556
2.823680556
3.543680556
4.363680556
0.233680556
0.953680556
-2.936319444
-4.796319444
-4.746319444
-0.916319444
-0.916319444
1.163680556
2.193680556
3.163680556
4.023680556
1.793680556
-2.856319444
-3.806319444
-4.586319444
-4.586319444
-2.856319444
-1.056319444
1.233680556
3.413680556
5.343680556
7.143680556
5.663680556
4.103680556
1.383680556
-2.546319444
-3.206319444
-3.246319444
-3.846319444
-3.536319444
-3.336319444
-2.966319444
1.683680556
2.793680556
1.043680556
-0.936319444
-1.826319444
-2.226319444
-3.986319444
];
n=length(sst);
mother = 'Morlet';
coef=cwt(sst,1:n,'cmorl');
[x,y]=meshgrid(1957:2011,1:n);
z=zeros(n,n);
z=griddata(1957:2011,1:n,coef,x,y,'v4');
[c,h] =contourf(x,y,z);
clabel(c,h)
hold on
[wave,period,scale,coi] = wavelet(sst,1,1,1,n-1,n-1,mother);
x1=1957:0.5:2011;
y1=interp1(1957:2011,coi,x1,'spline');
plot(x1,y1,'k')
set(gca,'YDir','reverse');
set(gca,'xtick',1957:5:2011,'xticklabel',1957:5:2011)
set(gca,'ytick',0:5,'yticklabel',0:5:40)

wendeping

想做时频分析图和小波方差曲线，

rs06

群里有呢，自己找找