求助：如何在matlab里面把小波实部彩色图变成有虚实的等值线

quanmaomao

代码是从论坛下载的：
t= x;%读取时间
zf= y; %读取气温
figure(1);
subplot(1,1,1);
index=1960:2013;
t=1:54;
c1=cwt(zf,[1:1:54],'morl','plot');%（morl是Morlet小波，mexh为墨西哥帽小波）
contourf(index,t,c1,25);
colormap(jet);
xlabel('年份');ylabel('尺度');
title('MORLET时频分布图');
colorbar('EASTOUTSIDE');
print(gcf,'-djpeg','MORLET时频分布图');
figure(2);%绘制小波时频分析图
subplot(1,1,1);
contourf(index,t,c1,25);
colormap(jet);
xlabel('年份');ylabel('尺度');
title('MORLET小波实部时频分布等高线');
colorbar('EASTOUTSIDE');
print(gcf,'-djpeg','MORLET时频Coef时频分布等高线');

figure(3);%绘制立体图
subplot(1,1,1);
meshc(c1);
surfc(c1);
shading interp;colormap(jet);
xlabel('年份');ylabel('尺度');zlabel('MORLET小波变换系数Coef');
title('MORLET小波变换系数时频分布立体图');
colorbar('EASTOUTSIDE');
print(gcf,'-djpeg','MORLET时频Coef时频分布立体图');

figure(4);
y=var(c1);%小波方差
plot(y);
title('小波方差图');
ylabel('Morlet小波方差');xlabel('年份');
grid;

figure(5);
subplot(1,1,1);
b=(abs(c1)).^2;%小波系数摸平方
[c,h]=contourf(index,t,b,25);
colormap(jet);
ylabel('尺度');xlabel('年份');
title('小波系数模平方');

quanmaomao

有没有大神帮帮忙，因为不会写代码，所以举步维艰

二爷名声在外

1.contour
2.正负分开画

quanmaomao

二爷名声在外 发表于 2017-3-29 10:02
1.contour
2.正负分开画

二爷，能详细说明一下不？？是在哪一部分添加什么代码呢？