爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3749|回复: 1

matlab 中的for循环问题,小波分析

[复制链接]

新浪微博达人勋

发表于 2016-3-11 16:38:39 | 显示全部楼层 |阅读模式
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
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-3-11 16:54:42 | 显示全部楼层
我找到原因了,,还是感谢大家,
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表