- 积分
- 496
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-1-11
- 最后登录
- 1970-1-1
|
MATLAB
问题截图: |
- |
问题概况: |
采用小波分析通过MATLAB编程,过程中for循环运行第二遍的时候受阻,在删除高亮标记的语句后就能正常运行,主要就是小波工具箱使用的那一句语句 ,求解释,难道使用了小波工具箱就不能用循环么?MATLAB显示的错误也请求大神解释一下!! |
我看过提问的智慧: |
看过 |
自己思考时长(天): |
1 |
系统平台: |
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
错误如下:
Error using cell/strcat (line 23)
Inputs must be cell arrays or strings.
Error in fei1 (line 19)
infile=strcat(sta(m),'-',period(mm),'-',scen(mmm),'.txt');
程序语句:
sta={'wz','ljd','mg'};
period={'shice','jizhun','weilai'};
scen={'A2','B2'};
for m=1:3
for mm=1:3
for mmm=1:2
infile=strcat(sta(m),'-',period(mm),'-',scen(mmm),'.txt');
x=infile{1};
sst=load (x,'-ascii');
if mm==1||mm==2
qishi=1961;
chidu=74;
else
qishi=2010;
chidu=79;
end
mmmm=m+mm+mmm-2;
variance = std(sst)^2;
sst = (sst - mean(sst))/sqrt(variance) ;
n = length(sst);
dt = 1;
time = [0:length(sst)-1]*dt + qishi ; % construct time array
xlim = [qishi,qishi+n-1]; % plotting range
pad = 1; % pad the time series with zeroes (recommended)
dj = 1/12; %填补模式为1。本时间序列长度为47,不足26,故选择0填补
s0 = 1/2*dt; %可分辨的最小尺度为0.5年
j1 = 6.5/dj; % 为26中的6,即以2的幂指数,根据时间序列长度确定
lag1 = 0.72; % lag-1 autocorrelation for red noise background
mother = 'Morlet';
% Wavelet transform:
[wave,period,scale,coi] = wavelet(sst,dt,pad,dj,s0,j1,mother);
%power = (abs(wave)).^2 ; % compute wavelet power spectrum计算小波系数的模的平方
%modulus=abs(wave); %计算小波系数的模
%variance1=sum(power')/n;%计算小波方差
end
end
end
|
|