%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');