这个程序到底哪里除了问题（wave）

阿吉木柯

load jupinzong.mat;
x1 =AOzong;
x2 =Tzong;
x3 =Rzong; %x4 = math; x5 = matl;

n = 732;

variance1 = std(x1)^2;
x1 = (x1 - mean(x1))/sqrt(variance1);
variance2 = std(x2)^2;
x2 = (x2 - mean(x2))/sqrt(variance2);
variance3 = std(x3)^2;
x3 = (x3 - mean(x3))/sqrt(variance3);


dt = 1/12;
time = [0:n-1]*dt + 1951.0 ;
xlim = [1951,2011];
pad = 1;      
dj = 1/12;
s0 = 6*dt;   
j1 = 6.5/dj;  
lag1 = 0.72;  
mother = 'Morlet';


[wave1,period,scale,coi] = wavelet(x1,dt,pad,dj,s0,j1,mother);
[wave2,period,scale,coi] = wavelet(x2,dt,pad,dj,s0,j1,mother);
[wave3,period,scale,coi] = wavelet(x3,dt,pad,dj,s0,j1,mother);


m = length(scale);


power1 = (abs(wave1)).^2 ;  
power2 = (abs(wave2)).^2 ;
power3 = (abs(wave3)).^2 ;


global_wst1 = (sum(power1')/n);  global_wss1 = (sum(power1)/m);
global_wst2 = (sum(power2')/n);  global_wss2 = (sum(power2)/m);
global_wst3 = (sum(power3')/n);  global_wss3 = (sum(power3)/m);

for j = 1:m ;
  pow1(j,:) = power1(j,:)./global_wst1(j); psd1(j,:) = power1(j,:)./scale(j);
  pow2(j,:) = power2(j,:)./global_wst2(j); psd2(j,:) = power2(j,:)./scale(j);
  pow3(j,:) = power3(j,:)./global_wst3(j); psd3(j,:) = power3(j,:)./scale(j);
end;

spow1 = fivesmooth(psd1,m,n);
spow2 = fivesmooth(psd2,m,n);
spow3 = fivesmooth(psd3,m,n);



c12=wave1.*wave2;
c13=wave1.*wave3;
gct12 = (sum(real(c12)')/n); gct13 = (sum(real(c13)')/n);
gcs12 = (sum(real(c12))/m); gcs13 = (sum(real(c13))/m);
红色部分总是提示我错误，可是我真的不知道怎么改了。难道不能直接用wave么？之前不是已经算了么～

wlzhongouc

把提示的错误贴出来吧

bakoo

然后把whos结果和提示结果都贴出来

夜尽天明183

楼主其实已经贴出来了，红色标记的部分