a=35;
b=3;
c=28;
x（1）=0;
y（1）=1.001;
z（1）=0;
h=0.0001;
for?i=1:100000
????K1=a*（y（i）-x（i））;
????L1=（c-a）*x（i）-x（i）*z（i）+c*y（i）;
????M1=x（i）*y（i）-b*z（i）;

????K2=a*（（y（i）+h/2*L1）-（x（i）+h/2*K1））;
????L2=（c-a）*（x（i）+h/2*K1）-（x（i）+h/2*K1）*（z（i）+h/2*M1）+c*（y（i）+h/2*L1）;
????M2=（x（i）+h/2*K1）*（y（i）+h/2*L1）-b*（z（i）+h/2*M1）;

????K3=a*（（y（i）+h/2*L2）-（x（i）+h/2*K2））;
????L3=（c-a）*（x（i）+h/2*K2）-（x（i）+h/2*K2）*（z（i）+h/2*M2）+c*（y（i）+h/2*L2）;
????M3=（x（i）+h/2*K2）*（y（i）+h/2*L2）-b*（z（i）+h/2*M2）;

????K4=a*（（y（i）+h*L3）-（x（i）+h*K3））;
????L4=（c-a）*（x（i）+h*K3）-（x（i）+h*K3）*（z（i）+h*M3）+c*（y（i）+h*L3）;
????M4=（x（i）+h*K3）*（y（i）+h*L3）-b*（z（i）+h*M3）;

????x（i+1）=x（i）+h/6*（K1+2*K2+2*K3+K4）;
????y（i+1）=y（i）+h/6*（L1+2*L2+2*L3+L4）;
????z（i+1）=z（i）+h/6*（M1+2*M2+2*M3+M4）;
end
%for?j=1:250000
%???x（j）=x（10000+j）;
%???y（j）=y（10000+j）;
%???z（j）=z（10000+j）;
%end
figure（1）;
plot3（xyz）;
title（‘Chen等式的三維相空間吸引子‘）;
xlabel（‘X（t）‘）;
ylabel（‘Y（t）‘）;
zlabel（‘Z（t）‘）;
grid;
figure（2）;
plot（x）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1087??2008-04-16?12:32??chaos\chen.m

?????文件?????????53??2012-06-02?18:19??chaos\logistic.asv

?????文件????????177??2008-04-14?21:56??chaos\logistic.m

?????文件???????1427??2008-04-15?23:06??chaos\lorenz.m

?????文件????????949??2008-04-15?23:18??chaos\rossler.m

?????目錄??????????0??2012-06-02?16:04??chaos

-----------?---------??----------?-----??----

?????????????????3693????????????????????6