这是一个求解分数阶混沌系统的Matlab实现程序

这是一个求解分数阶混沌系统的Matlab实现程序，希望有帮助

function fractional_chaos（xyz）
   clearclf
   hold on
   h=0.02;t=0:0.02:90;
   N=length（t）;
   x（1）=-15;y（1）=-18;z（1）=15;
   a=35;b=3;c=27;q=0.9;
   p=min（[2 1+q]）; 
      for n=2:N
        for j=1:n-1
          if j==1alpha（（n-1）j）=（n-1）^（q+1）-（n-q-1）*n^q;                 
          else alpha（（n-1）j）=（n-j+2）^（q+1）+（n-j）^（q+1）-2*（（n-j+1）^（q+1））;              
          end
          
          beta（（n-1）j）=（（h^q）/q）*（（n-j+1）^q-（n-j）^q）; 
          
     
          u（j）=a*（y（j）-x（j））;
          v（j）=（c-a-z（j））*x（j）+c*y（j）;
          w（j）=x（j）*y（j）-b*z（j）;
          
          C1（j）=alpha（（n-1）j）*u（j）;          
          C2（j）=alpha（（n-1）j）*v（j）;
          C3（j）=alpha（（n-1）j）*w（j）;
          
          xA=sum（C1）;
          yA=sum（C2）;
          zA=sum（C3）; 
          
          D1（j）=beta（（n-1）j）*u（j）;          
          D2（j）=beta（（n-1）j）*v（j）;
          D3（j）=beta（（n-1）j）*w（j）;
          
          xB=sum（D1）;
          yB=sum（D2）;
          zB=sum（D3）;          
        end           
          xp（n）=x（1）+1/gamma（q）*xB;
          yp（n）=y（1）+1/gamma（q）*yB;
          zp（n）=z（1）+1/gamma（q）*zB;
          
          x（n）=x（1）+h^q/gamma（2+q）*a*（yp（n）-xp（n））+h^q/gamma（2+q）*xA;
          y（n）=y（1）+h^q/gamma（2+q）*（（c-a）*xp（n）-xp（n）*zp（n）+c*yp（n））+h^q/gamma（2+q）*yA;
          z（n）=z（1）+h^q/gamma（2+q）*（xp（n）*yp（n）-b*zp（n））+h^q/gamma（2+q）*zA;
       end
     subplot（221）plot（xy）;
      xlabel（‘x‘）;ylabel（‘y‘）;
     
     subplot（222）plot（xz）;
     xlabel（‘x‘）;ylabel（‘z‘）;
     
     subplot（223）plot（yz）;
     xlabel（‘y‘）;ylabel（‘z‘）;
     
     subplot（224）plot3（xyz）;
     xlabel（‘x‘）;ylabel（‘y‘）;zlabel（‘z‘）;
     axis（[-20 20 -25 25 5 40]）
     hold off 
             
   

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1776  2007-09-07 21:41  fractional_chaos.m

----------- ---------  ---------- -----  ----

                 1994                    2