资源描述：
《分数阶混沌程序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+h^q1*35*(y0-x0)/(gamma(q1)*q1);y1(1)=y0+h^q2*(-7*x0-x0*z0+28*y0)/(gamma(q2)*q2);z1(1)=z0+h^q3*(x0*y0-3*z0)/(gamma(q3)*q3);x(1)=x0+35*h^q1*(y1(1)-x1(1)+q1*
2、(y0-x0))/gamma(q1+2);y(1)=y0+h^q2*(-7*x1(1)-x1(1)*z1(1)+28*y1(1)+q2*(-7*x0-x0*z0+28*y0))/gamma(q2+2);z(1)=z0+h^q3*(x1(1)*y1(1)-3*z1(1)+q3*(x0*y0-3*z0)/gamma(q3+2));forn=1:NM1=(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*
3、y0-3*z0);N1=((n+1)^q1-n^q1)*35*(y0-x0);N2=((n+1)^q2-n^q2)*(-7*x0-x0*z0+28*y0);N3=((n+1)^q3-n^q3)*(x0*y0-3*z0);forj=1:nM1=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)^(q
4、3+1)-2*(n-j+1)^(q3+1))*(x(j)*y(j)-3*z(j));N1=N1+((n-j+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+h^q1*N1/(gamma(q1)*q1);y1(n+1)=y0+h^q2*N2/(gamma(q2)*q2);z1(n+1)=z0+h^q3*N3/(gamma(q3)*q3);x(n+
5、1)=x0+h^q1*(35*(y1(n+1)-x1(n+1))+M1)/gamma(q1+2);y(n+1)=y0+h^q2*(-7*x1(n+1)-x1(n+1)*z1(n+1)+28*y1(n+1)+M2)/gamma(q2+2);z(n+1)=z0+h^q3*(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');sub
6、plot(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'),zlabel('z');gridonA=[max(x)max(y)max(z)];B=[min(x)min(y)min(z)];电路仿真
