分数阶 统一混沌系统 matlab 代码 .txt

function dy=united-fra-chaos

q1=0.9;q2=0.9;q3=0.8;

h=0.01;N=2000;

a=1;

x0=2;y0=1;z0=3;

%x0=-3.5;y0=4.2;z0=2.5;

M1=0;M2=0;M3=0;

x(N+1)=[0];y(N+1)=[0];z(N+1)=[0];

x1(N+1)=[0];y1(N+1)=[0];z1(N+1)=[0];

 

x1(1)=x0+h^q1*(25*a+10)*(y0-x0)/(gamma(q1)*q1);

y1(1)=y0+h^q2*((28-35*a)*x0-x0*z0+(29*a-1)*y0)/(gamma(q2)*q2);

z1(1)=z0+h^q3*(x0*y0-(8+a)*z0/3)/(gamma(q3)*q3);

 

x(1)=x0+h^q1*((25*a+10)*(y1(1)-x1(1))+q1*(25*a+10) *(y0-x0))/gamma(q1+2);