为什么小波功率谱图下边部分是空白？

黑今之马户

部分源码在此
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

一水天下

时间间隔1年的数据应该不能分析出1年以下的周期了吧，我猜的哈哈

消魂张

很希望斑竹能贴出wavelet power spectrum的实现部分代码，学习一下

造字仓颉

好东西。希望能够好好学习一下。谢谢楼主。

fhzf1234

遇到同样问题，请问解决了么

孙玻玻

我认为不是你画图代码的问题，应该是前面有一个j1=（数字）/dj 的问题，那个数字你要根据你的研究时段来更改一下