分?jǐn)?shù)階混沌程序MATLAB

1、q1=0.9;q2=0.9;q3=0.9;h=0.01;n=4000;x0=2;y0=1;z0=3;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+hq1*35*(y0-x0)/(gamma(q1)*q1);y1(1)=y0+hq2*(-7*x0-x0*z0+28*y0)/(gamma(q2)*q2);z1(1)=z0+hq3*(x0*y0-3*z0)/(gamma(q3)*q3);x(1)=x0+35*hq1*(y1(1)-x1(1)+q1*(y0-x0)/gamma(q1+2)

2、;y(1)=y0+hq2*(-7*x1(1)-x1(1)*z1(1)+28*y1(1)+q2*(-7*x0-x0*z0+28*y0)/gamma(q2+2);z(1)=z0+hq3*(x1(1)*y1(1)-3*z1(1)+q3*(x0*y0-3*z0)/gamma(q3+2);for n=1:n m1=(n(q1+1)-(n-q1)*(n+1)q1)*35*(y0-x0); m2=(n(q2+1)-(n-q2)*(n+1)q2)*(-7*x0-x0*z0+28*y0); m3=(n(q3+1)-(n-q3)*(n+1)q3)*(x0*y0-3*z0); n1=(n+1)q1-nq1)*35*

3、(y0-x0); n2=(n+1)q2-nq2)*(-7*x0-x0*z0+28*y0); n3=(n+1)q3-nq3)*(x0*y0-3*z0);for j=1:n m1=m1+(n-j+2)(q1+1)+(n-j)(q1+1)-2*(n-j+1)(q1+1)*35*(y(j)-x(j);m2=m2+(n-j+2)(q2+1)+(n-j)(q2+1)-2*(n-j+1)(q2+1)*(-7*x(j)-x(j)*z(j)+28*y(j); m3=m3+(n-j+2)(q3+1)+(n-j)(q3+1)-2*(n-j+1)(q3+1)*(x(j)*y(j)-3*z(j); n1=n1+(n-j

4、+1)q1-(n-j)q1)*35*(y(j)-x(j); n2=n2+(n-j+1)q2-(n-j)q2)*(-7*x(j)-x(j)*z(j)+28*y(j); n3=n3+(n-j+1)q3-(n-j)q3)*(x(j)*y(j)-3*z(j);endx1(n+1)=x0+hq1*n1/(gamma(q1)*q1);y1(n+1)=y0+hq2*n2/(gamma(q2)*q2);z1(n+1)=z0+hq3*n3/(gamma(q3)*q3);x(n+1)=x0+hq1*(35*(y1(n+1)-x1(n+1)+m1)/gamma(q1+2);y(n+1)=y0+hq2*(-7*x1(

5、n+1)-x1(n+1)*z1(n+1)+28*y1(n+1)+m2)/gamma(q2+2);z(n+1)=z0+hq3*(x1(n+1)*y1(n+1)-3*z1(n+1)+m3)/gamma(q3+2);end%subplot(1,1,1),plot3(x,y,z);xlabel('x'),ylabel('y');zlabel('z');subplot(2,2,1),plot(x,y);xlabel('x'),ylabel('y');subplot(2,2,2),plot(x,z);xlabel('x'),ylabel('z');subplot(2,2,3),plot(y,z);ylabel('y'),ylabel('z');subplot(2,2,4),plot3(x,y,z);xlabel('x'),ylabel('y