新人请教

yjj1008

load q.txt
sst =q;
n=length(q);
variance = std(sst)^2;
sst = (sst - mean(sst))/sqrt(variance);
dt = 1;
year = [0:n-1]+ 1962.0 ;
xlim = [1962,2011];   %
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]=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 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
我按照这个输入的，在设置路径加入了 wavelet.m文件
请大神帮我找出问题，谢谢

Wetter

索引超出矩阵纬度你没法理解么那个矩阵维度你又有没有检查一下？

yjj1008

Wetter 发表于 2019-4-6 14:16
索引超出矩阵纬度你没法理解么那个矩阵维度你又有没有检查一下？

我是刚学的，看到有现成的程序，我就直接套进去了，是否是那个程序有地方要修改还是我的文件，请指点一二，谢谢

Wetter

yjj1008 发表于 2019-4-6 14:56
我是刚学的，看到有现成的程序，我就直接套进去了，是否是那个程序有地方要修改还是我的文件，请指点一二 ...

直接套进去就能用么，你不改一下么

youlanqinghai

超出矩阵维度这个问题前几天刚遇到过。你可以检查一下加载进来的原始数据，添加一句转置的code，看看能不能解决