马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
在论坛上下了这个代码,可是有些地方没看懂,希望有人帮忙解释下,不胜感激。代码如下:- %Chen's系统多步预报的主程序
- clear all
- fid=fopen('chen.txt','r');
- a=fread(fid);
- b=char(a');
- sj=str2num(b);%数据读取完毕
- fclose(fid);
- roll=sj(:,2);%取X数据
- oldroll=roll;
- AverRoll=roll-mean(roll);
- roll=AverRoll;
- mtbp=70;
- taow=mtbp;
- deltaT=0.001;
- tao = 10;
- taowdivtao=round(taow/tao);
- d=8;%taowdivtao+1;
- prestep=5000;
- N = length(roll)-prestep;
- MaxStep = 2000;
- Step = 1 : 1 : MaxStep;
- [PredictedData] = FunctionChaosPredict(roll(1:N),N,mtbp,deltaT,tao,d,MaxStep);%调用AOLMM进行多步预报
- subplot(2,1,1);plot(roll(1:9000+MaxStep)+mean(roll));
- hold on;
- subplot(2,1,1);
- plot(Step + N,PredictedData+mean(oldroll),'r:');
- hold on;
- ylabel('Chen''s System');
- legend('Original','Predict');
- axis([0 N+MaxStep -40 40]);
- subplot(2,1,2);
- plot(Step + N,PredictedData+mean(oldroll),'r:');
- axis([N N+MaxStep -40 40]);
- hold on;
- subplot(2,1,2);
- plot(Step + N,oldroll(Step + N),'b-');
- ylabel('Chen''s System');
- legend('Predict','Original');
复制代码 不懂得地方:‘’prestep=5000;‘’这个值是怎么确定的,我比较了一下,是减去后5000个数据,想请教下这个值如确定
|