求大神指导 matlab小波分析出图问题

shany

用论坛上的程序，自己改了数字后，运行结果如图，为什么纵坐标 上部缺失？是程序哪里弄错了吗 ？新手小白求指导，用的是55年的温度数据


程序如下：

sst=wendu;
n=length(wendu);
variance = std(sst)^2;
sst = (sst - mean(sst))/sqrt(variance);
dt = 1;
year = [0:n-1]+ 1961.0 ;
xlim = [1961,2015];   %
pad = 1;      
dj = 1/12;     
s0 = 1/2*dt;   
j1 = 6.5/dj;   
lag1 = 0.72;  
mother = 'Morlet';
[wave,period,scale,coi] = wavelet(sst,dt,pad,dj,s0,j1,mother);
power = (abs(wave)).^2 ;  %计算小波系数的模的平方
modulus=abs(wave);  %计算小波系数的模
variance1=sum(power')/n;%计算小波方差
%画小波系数实部等值线图
fig1=figure;
%subplot(3,1,1)
levels = [0,0.5,1.0,1.5,2.0,2.5];
v = [0,0.5,1.0,1.5];
Yticks = 0:5:55;
[c,h]=contour(year,period,real(wave),levels,'k-');
clabel(c,h,v,'fontsize',5);
xlabel('年份/year')
ylabel('周期/年 period/year')
title('(a)')
set(gca,'XLim',xlim(:))
set(gca,'YLim',[0 55], ...
    'YDir','default', ...
'YTick',Yticks(:), ...
'YTickLabel',Yticks)
hold on
levels = [-0.5,-1.0,-1.5,-2.0,-2.5];
v = [-0.5,-1.0,-1.5];
[c,h] = contour(year,period,real(wave),levels,'r--');
clabel(c,h,v,'fontsize',5);
hold on

nunu18

你这个时间尺度是55年，纵坐标肯定缺啊，而且你再分析50年以上的尺度也没什么意义，你调整纵坐标就行了

晴天

这个小波图的周期是怎么看的？

su27peng

大家一起研究研究啊！

shany

nunu18 发表于 2018-3-8 15:30
你这个时间尺度是55年，纵坐标肯定缺啊，而且你再分析50年以上的尺度也没什么意义，你调整纵坐标就行了

好的，谢谢您解答

shany

nunu18 发表于 2018-3-8 15:30
你这个时间尺度是55年，纵坐标肯定缺啊，而且你再分析50年以上的尺度也没什么意义，你调整纵坐标就行了

您好，还有一个就是小波方差图，横坐标后的也确实，这怎么调整啊 ，想模拟一个周

nunu18

你是用matlab画的图，我对matlab很陌生，不知道怎么改，但是论坛有很多小波的程序，也有高人指出过小波尺度最终怎么确定，具体的算法忘记了好像是年限的一半儿，自己查查