急求！怎么在代码里改小波分析的时间尺度的？！！

狂奔的蜗牛

我根据http://bbs.06climate.com/forum.php?mod=viewthread&tid=11368做的一个小波分析，代码如下：
clear;
load q.txt
sst =q;
n=length(q);
variance = std(sst)^2;
sst = (sst - mean(sst))/sqrt(variance);
dt = 1;
year = [0:n-1]+ 1644.0 ;
xlim = [1644,1911];   %
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;%计算小波方差
%画小波系数实部等值线图
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:45;
[c,h]=contourf(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 50], ...
    '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
% 画小波方差图
subplot(3,1,2)
plot(period,variance1,'k-')
hold on;
levels= [1,5,10,15,20,25,30,35,40,45];
title('(b)')
set(gca,'XLim',[1,50], ...        
    'XTick',levels,...
   'XTickLabel',levels)
xlabel('周期/a')
ylabel('方差 variance')
hold on
%画小波系数模
subplot(3,1,3)
levels = [0,0.5,1.0,1.5,2.0,2.5];
v = [0,0.5,1.0,1.5];
Yticks = [0:5:30];
[c,h]=contour(year,period,abs(wave),levels,'k-');
clabel(c,h,v,'fontsize',5);
title('(c)')
xlabel('年份/year')
ylabel('周期/年 period/year')
set(gca,'XLim',xlim(:))
set(gca,'YLim',[0 30], ...
    'YDir','default', ...
'YTick',Yticks(:), ...
'YTickLabel',Yticks)
结果图是这样的：

由于该代码的时间分析尺度是50年的，我的数据是270年的，所以我想把尺度改成100年的，我自己尝试了好几遍都改不上来。
我把上面红色部分改成：
set(gca,'YLim',[0 100], ...
    'YDir','default', ...
'YTick',Yticks(:), ...
'YTickLabel',Yticks)
出来的图如下：

真心求大神帮忙改改，跪谢了！！！

狂奔的蜗牛

有没有高手帮帮忙啊！！

bbs

如果你红色字体改成[0 100]，那“Yticks = 0:5:45;“也要把45改到100.。。。应该是这样的，我没试过，现在也是初学阶段。。。。

lyqx

我也要高手回答啊

风儿飘飘

正如楼上所说的，修改yticks的范围即可，与你的小波实部范围相匹配

蜗牛快跑

请问楼主的问题解决了吗？我修改Yticks = 0:20:240;和set(gca,'YLim',[0 240], 得到的图还是只显示0-40年周期部分.请问如何修改？谢谢！

fhzf1234

请问楼主的问题解决了吗？

杜佳泽

有没有运行时提示“未定义函数或变量 'wavelet'”的，我是初学者，请大神们指正！！！

zoiujing456

过来学习学习