%%計算lymax
%?function?jeffcott（）
global?a;
Z=[];
d0=1e-8;
ti=0.1;
X1=[];
for?a=linspace（85012）
lsum=0;
x=0;y=1;z=0;
x1=0;y1=1;z1=d0;
for?i=1:1000
%????[T1Y1]=ode45（@Lorenz1[x;y;z;a;3;28]）;?
%????[T2Y2]=ode45（@Lorenz1[x1;y1;z1;a;3;28]）;
x0=[x;y;z];
x01=[x1;y1;z1];
Y1=solveLorenz（x0ti）;
Y2=solveLorenz（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（85012）;
figure;plot（aZ）;
title（‘Lorenz?系統最大lyapunov指數‘）
xlabel（‘parameter?a‘）ylabel（‘lyapunov?exponents‘）
%?figure;plot（aX1）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2015-12-17?09:19??Lyapunov_max_Lorenz\
?????文件?????????834??2015-12-17?09:18??Lyapunov_max_Lorenz\Lorenz_lyapunov.m
?????文件?????????373??2015-12-17?09:16??Lyapunov_max_Lorenz\solveLorenz.m