请问一个二维系统的稳定性,理论能够证明,但matlab仿真怎么不行
各位师兄师姐:这个的系统是
dx/dt=y;
dy/dt=-1.5(y^(5/3)+0.15x)^(1/5);
这里可以证明存在一个正定的Lyapunov 函数V,并且dV/dt<0,文献也画出图了,我怎么就画不出来呢?
我的程序是制定一个M文件,命名为gouzhao08.m
function dx=gouzhao08(t,x);
dx=[x(2);
-1.5*(x(2)^(5/3)+0.15*x(1))^(0.2)];
再调用
clear al1
=ode45(@gouzhao08,0:0.01:10,);
plot(t,x(:,2),'k')
得到是一个不收敛的图象,请师兄师姐帮我一下,谢谢! 从程序来看没有问题，会不会初始值有影响呢？ 不懂!
解出复数解是否合理?
页:
[1]