概述
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))/gamma(q3+2);
最后
以上就是狂野毛巾为你收集整理的分数阶 matlab,分数阶统一混沌系统matlab程序.doc的全部内容,希望文章能够帮你解决分数阶 matlab,分数阶统一混沌系统matlab程序.doc所遇到的程序开发问题。
如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。
发表评论 取消回复