请问这个系统怎么得到分维数,怎么判断是不是混沌啊?
我是新手，感觉自己太菜了，请高手指点一下，怎么判断分形维或者其他混沌特征啊,以下是M文件和做出的图相t0=0;
tf=5;
tspan=
Y0=
=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)] http://forum.vibunion.com/forum/?thread-15918-1-1.html
页:
[1]