请问这个系统怎么得到分维数,怎么判断是不是混沌啊?

lyxxl

我是新手，感觉自己太菜了，请高手指点一下，怎么判断分形维或者其他混沌特征啊,以下是M文件和做出的图相
t0=0;
tf=5;
tspan=[t0,tf]
Y0=[35,0,1256,5,0,628,0,0]
[t,YY]=ode45('DYdt',tspan,Y0)
X=YY(:,1);
Y=YY(:,4);
Z=YY(:,7);
C=YY(:,8);
plot(X,Y,'b','linewidth',2);
hold on

function Yd=Dydt(t,Y)
Yd=[-(-36*0.1*0.5*50*(1+0.2*Y(2))+36*0.1*0.5*Y(1)-2*25*50*0.2*Y(3)-0.001^2*25*50*0.2*Y(3))/0.001^2/(2*3000+25);
    Y(3);
    -4*3.14^2*100^2*Y(2);
   (36*3.14*0.1*0.001*0.5*10*(1+0.1*Y(5))-36*3.14*0.1*0.001*0.5*Y(4)+2*3.14*0.001*25*10*0.1*Y(6)+3.14*0.001^3*25*10*0.1*Y(6)+12*(-3.14/6*0.001^3*3000*9.8+1.61*(0.1*25)^0.5*0.001^2*(50*(1+0.2*Y(2))-Y(1))*(10*0.2*Y(3)/Y(4))^0.5))/3.14/0.001^3/(2*3000+25);
   Y(6);
   -4*3.14^2*50^2*Y(5);
   Y(1);
   Y(4)]

gghhjj

http://forum.vibunion.com/forum/?thread-15918-1-1.html