%%計算lymax
%?function?jeffcott（）
global?a;
Z=[];
d0=1e-8;
ti=0.1;
X1=[];
for?a=linspace（255026）
lsum=0;
x=0;y=1;z=0;
x1=0;y1=1;z1=d0;
for?i=1:1000
%????[T1Y1]=ode45（@Chen1[x;y;z;a;3;28]）;?
%????[T2Y2]=ode45（@Chen1[x1;y1;z1;a;3;28]）;
x0=[x;y;z];
x01=[x1;y1;z1];
Y1=solvechen（x0ti）;
Y2=solvechen（x01ti）;
???n1=length（Y1）;n2=length（Y2）;
???x=Y1（n11）;y=Y1（n12）;z=Y1（n13）;%????X=x;
???x1=Y2（n21）;y1=Y2（n22）;z1=Y2（n23）;
???d1=sqrt（（x-x1）^2+（y-y1）^2+（z-z1）^2）;
???x1=x+（d0/d1）*（x1-x）;
???y1=y+（d0/d1）*（y1-y）;
???z1=z+（d0/d1）*（z1-z）;
???if?i>500
???????lsum=lsum+log（d1/d0）;
???end
end
%?X1=[X1?X];
Z=[Z?lsum/（i-500）];
end
a=linspace（255026）;
figure;plot（aZ）;
title（‘Chen?系統最大lyapunov指數‘）
xlabel（‘parameter?a‘）ylabel（‘lyapunov?exponents‘）
%?figure;plot（aX1）;

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

?????文件????????826??2014-02-24?12:31??chen_lyapunov.m

?????文件????????371??2014-02-23?22:24??solvechen.m

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

?????????????????1197????????????????????2