- 积分
- 2304
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-3-13
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 黑今之马户 于 2014-11-27 14:51 编辑
部分源码在此
subplot('position',[0.1 0.10 0.65 0.38])
levels = [1,2,4,8,16] ;
v=[0,1,2,4];
Yticks = 2.^(fix(log2(min(period))):fix(log2(max(period))));
[c,h]=contour(time,log2(period),log2(power),log2(levels),'k-'); %*** or use 'contourfill'
%imagesc(time,log2(period),log2(power)); %*** uncomment for 'image' plot
clabel(c,h,v);
xlabel('年份')
ylabel('周期')
%title('b) NINO3 SST Wavelet Power Spectrum')
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] ;
v=[0,-2,-4];
[c,h]=contour(time,log2(period),log2(power),log2(levels),'k:');
clabel(c,h,v);
hold on
plot(time,log2(coi),'k')
hold off
%图片质量
print('-depsc','-tiff','-r300','CuteFred')
print -dbmp -r300 output |
|