请问一个二维系统的稳定性,理论能够证明,但matlab仿真怎么不行

xiangwei27 · 发表于 2009-9-8 01:13

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

各位师兄师姐:
这个的系统是
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
[t,x]=ode45(@gouzhao08,0:0.01:10,[5 -3]);
plot(t,x(:,2),'k')
得到是一个不收敛的图象,请师兄师姐帮我一下,谢谢!

无水1324 · 发表于 2009-9-11 16:52

从程序来看没有问题，会不会初始值有影响呢？

ChaChing · 发表于 2009-9-11 18:45

不懂!
解出复数解是否合理?

账号		自动登录	找回密码
密码			我要加入

[动力学和稳定性] 请问一个二维系统的稳定性,理论能够证明,但matlab仿真怎么不行

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。