shangchuan

用matlab实现混沌系统的李雅普诺夫指数图 分数阶程序。

global x y z
a=1; %parameter a
q1=0.9;q2=0.9;q3=0.9; %the fractional order
x0=3;y0=1;z0=5;
h=0.01;N=1000;
x（N+1）=[0];y（N+1）=[0];z（N+1）=[0];
x1（N+1）=[0];y1（N+1）=[0];z1（N+1）=[0];
M1=0;M2=0;M3=0;N1=0;N2=0;N3=0;
 
x1（1）=x0+h^q1*（25*a+10）*（y0-x0）/（q1*gamma（q1））;
y1（1）=y0+h^q2*（（28-35*a-z0）*x0+（29*a-1）*y0）/（q2*gamma（q2））;
z1（1）=z0+h^q3*（x0*y0-z0*（8+a）/3）/（q3*gamma（q3））;
 
x1=x0+h^q1*（25*a+10）*（y1（1）-x1（1）+q1*（y0-x0））/gamma（q1+2）;
y1=y0+h^q2*（（28-35*a-z1（1））*x1（1）+（29*a-1）*y1（1）+q2*（（28-35*a-z0）*x0+（29*a-1）*y0））/gamma（q2+2）;
z1=z0+h^q3*（x1（1）*y1（1）-z1（1）*（8+a）/3+q3*（x0*y0-z0*（8+a）/3））/gamma（q3+2）;
 
for n=1:N
    M1=（n^（q1+1）-（n-q1）*（n+1）^q1）*（y0-x0）;
    M2=（n^（q2+1）-（n-q2）*（n+1）^q2）*（（28-35*a-z0）*x0+（29*a-1）*y0）;
    M3=（n^（q3+1）-（n-q3）*（n+1）^q3）*（x0*y0-z0*（8+a）/3）;
    
 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1814  2014-10-30 11:15  上传\BIANLIANG.m
     文件         826  2014-02-24 12:31  上传\chen_lyapunov.m
     文件         480  2014-10-30 11:43  上传\HANSHU.m
     文件         370  2014-11-06 09:06  上传\solvechen.m
     文件        2496  2014-09-27 18:36  上传\Unified.m