分数阶统一混沌系统matlab程序.docVIP

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); y(1)=y0+h^q2*((28-35*a)*x1(1)-x1(1)*z1(1)+(29*a-1)*y1(1)+q2*((28-35*a)*x0-x0*z0+(29*a-1)*y0))/gamma(q2+2); z(1)=z0+h^q3*(x1(1)*y1(1)-(8+a)*z1(1)/3+q3*(x0*y0-(8+a)*z0/3)