• 大小: 1KB
    文件类型: .zip
    金币: 1
    下载: 0 次
    发布日期: 2021-06-02
  • 语言: Matlab
  • 标签: lyapunov  matlab  

资源简介

chen系统最大lyapunov指数。求解方法为定义法。两条相轨线的步长为初设距离d0的基础上加上相对分量。求解时运行chen_lyapunov.m。可移植性强,只需要改变变量数目即可。

资源截图

代码片段和文件信息

%%计算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  liyapunupfu\chen_lyapunov.m
     文件         371  2014-02-23 22:24  liyapunupfu\solvechen.m
     目录           0  2018-04-09 23:40  liyapunupfu\

评论

共有 条评论