chenpredictmain.asv

来自「混沌多步预测算法 可实现较为精确的多步预测 实用性较强」· ASV 代码 · 共 36 行

%Chen's系统多步预报的主程序
clear all



roll=[5023	4647	9100	14100	15842	17800	20000	48104	32737	46789	41418	64608	68222	69613	73624	87750	92616	96674	112044	126382	138127	144547	150648	180512	179801	203523	218573	233704	257705	277089	304442	330408.3	376414.6	389147	416757.8	420063.98	429687.09	445290.53	467347.6	498337.4	538375.4	575297.6	652365.8]';%取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');