求助帖，这个小波分析怎么改啊

狂奔的蜗牛

>> 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;%计算小波方差
%画小波系数实部等值线图
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:30;
[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 30], ...
    '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
% 画小波方差图
fig2=figure;
%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
%画小波系数模
fig3=figure;
%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)
??? Error using ==> minus
Matrix dimensions must agree.
老是出错，菜鸟也不知道怎么改。另附上我的数据（第一幅图），想做个周期分析还有像第二幅图那样的结果的。不知道我的这些数据能不能做出来。求大家帮帮忙

kongfeng0824

这个是实部系数图。
还有一个方差图。
你的出现什么错误提示？

狂奔的蜗牛

kongfeng0824 发表于 2013-4-16 23:20

                               
登录/注册后可看大图

这个是实部系数图。
还有一个方差图。
你的出现什么错误提示？

??? Error using ==> minus
Matrix dimensions must agree.
老是出现这个错，我也不知道怎么改

狂奔的蜗牛

kongfeng0824 发表于 2013-4-16 23:20

                               
登录/注册后可看大图

这个是实部系数图。
还有一个方差图。
你的出现什么错误提示？

我的QQ754613382，能咨询咨询你吗？

kongfeng0824

狂奔的蜗牛 发表于 2013-4-17 08:50

                               
登录/注册后可看大图

??? Error using ==> minus
Matrix dimensions must agree.
老是出现这个错，我也不知道怎么改

你的数组维数不对。需要改成符合的维数。

狂奔的蜗牛

kongfeng0824 发表于 2013-4-17 10:57

                               
登录/注册后可看大图

你的数组维数不对。需要改成符合的维数。

数组维数改了之后，现在又出现??? Undefined command/function 'wavelet'.
这个错误了。这是什么错误啊

kongfeng0824

狂奔的蜗牛 发表于 2013-4-17 13:58

                               
登录/注册后可看大图

数组维数改了之后，现在又出现??? Undefined command/function 'wavelet'.
这个错误了。这是什么错误啊

wavelet是一个函数，你查一下是否是外调的函数？如果不是，就在你的这个程序之中。看看是不是哪里名称出现错误了。

阿吉木柯

狂奔的蜗牛 发表于 2013-4-17 13:58

                               
登录/注册后可看大图

数组维数改了之后，现在又出现??? Undefined command/function 'wavelet'.
这个错误了。这是什么错误啊

不知道你解决问题了没有，这个就是你要把wavelet等.m文件和你的程序文件放在一起。因为你的程序中要调用它们。

金陵蔡肉包

wavelet.m是第三方的程序，要add to path

zl.87625

一直用fortran，改天也试一下这个matlab 的