分数阶李指数的程序

大江东去

最近在搞分数阶混沌，不太了解软件，这是一个matlab程序，不知道怎么运行，求大神给看看如何出来分数阶李指数的结果。
程序代码：

1. n=3500;x=zeros(n+1,4);
   
2. x(1,1)=0;x(1,2)=5;x(1,3)=8;x(1,4)=10;q1=0.95;q2=0.95;q3=0.95;a=35;b=3;c=28;h=0.01;
   
3. m=1;r=1;w=1;xx(1,2)=5;xx(1,3)=8;xx(1,4)=10;kk1=1.8;kk2=0;kk3=0;tau1=0.45;tau2=0;tau3=0;
   
4. for i=2:n+1
   
5.     x(i,1)=x(i-1,1)+h;
   
6.     k1=0;
   
7.     for j=0:i-2
   
8.         k1=k1+1/gamma(q1)*h^q1/q1*[(i-j-1)^q1-(i-j-2)^q1]*(a*[x(j+1,3)-x(j+1,2)]+kk1*(x(j+1,2)-xx(j+1,2)));
   
9.     end
   
10.     k1=k1+x(1,2);
    
11.     k2=0;
    
12.     for j=0:i-2
    
13.         k2=k2+1/gamma(q2)*h^q2/q2*[(i-j-1)^q2-(i-j-2)^q2]*([(c-a)*x(j+1,2)-x(j+1,4)*x(j+1,2)+c*x(j+1,3)]+kk2*(x(j+1,3)-xx(j+1,3)));
    
14.     end
    
15.     k2=k2+x(1,3);
    
16.     k3=0;
    
17.     for j=0:i-2
    
18.         k3=k3+1/gamma(q3)*h^q3/q3*[(i-j-1)^q3-(i-j-2)^q3]*([x(j+1,2)*x(j+1,3)-b*x(j+1,4)]+kk3*(x(j+1,4)-xx(j+1,4)));
    
19.     end
    
20.     k3=k3+x(1,4);
    
21.     k4=0;
    
22.     for j=0:i-2
    
23.         if j==0
    
24.             z=(i-2)^(q1+1)-(i-2-q1)*(i-1)^q1;
    
25.         else
    
26.             z=(i-j)^(q1+1)+(i-j-2)^(q1+1)-2*(i-j-1)^(q1+1);
    
27.         end
    
28.         k4=k4+z*(a*[x(j+1,3)-x(j+1,2)]+kk1*(x(j+1,2)-xx(j+1,2)));
    
29.     end
    
30.     k5=0;
    
31.     for j=0:i-2
    
32.         if j==0
    
33.             z=(i-2)^(q2+1)-(i-2-q2)*(i-1)^q2;
    
34.         else
    
35.             z=(i-j)^(q2+1)+(i-j-2)^(q2+1)-2*(i-j-1)^(q2+1);
    
36.         end
    
37.         k5=k5+z*([(c-a)*x(j+1,2)-x(j+1,4)*x(j+1,2)+c*x(j+1,3)]+kk2*(x(j+1,3)-xx(j+1,3)));
    
38.     end
    
39.     k6=0;
    
40.     for j=0:i-2
    
41.         if j==0
    
42.             z=(i-2)^(q3+1)-(i-2-q3)*(i-1)^q3;
    
43.         else
    
44.             z=(i-j)^(q3+1)+(i-j-2)^(q3+1)-2*(i-j-1)^(q3+1);
    
45.         end
    
46.         k6=k6+z*([x(j+1,2)*x(j+1,3)-b*x(j+1,4)]+kk3*(x(j+1,4)-xx(j+1,4)));
    
47.     end
    
48.     x(i,2)=x(1,2)+h^q1/gamma(q1+2)*[a*(k2-k1)+k4];
    
49.     x(i,3)=x(1,3)+h^q2/gamma(q2+2)*[(c-a)*k1-k1*k3+c*k2+k5];
    
50.     x(i,4)=x(1,4)+h^q3/gamma(q3+2)*(k1*k2-b*k3+k6);
    
51.     if x(i,1)<=tau1
    
52.         xx(i,2)=x(i,2);
    
53.         m=m+1;
    
54.     else
    
55.         xx(i,2)=x(i-m,2);
    
56.     end
    
57.     if x(i,1)<=tau2
    
58.         xx(i,3)=x(i,3);
    
59.         r=r+1;
    
60.     else
    
61.         xx(i,3)=x(i-r,3);
    
62.     end
    
63.     if x(i,1)<=tau3
    
64.         xx(i,4)=x(i,4);
    
65.         w=w+1;
    
66.     else
    
67.         xx(i,4)=x(i-w,4);
    
68.     end
    
69. end
    
70. Y=x(500:3500,2:4);
    
71. l=lyapunovdelay(1e-2,Y)
    

复制代码

一生有你

这个lyapunovdelay(1e-2,Y)函数怎么没有定义

lihaitao123

楼主的分数阶le指数的程序有没有完整的呢？384331310@qq.com

大江东去

做不出来啊