微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

MATLAB中方程的稳态和特征值

如何解决MATLAB中方程的稳态和特征值

enter image description here

我尝试使用以下代码解决上述给定问题。

syms y r K theta
eqn = r*y*(y/theta - 1)*(1-y/K) == 0;
S=solve(eqn,y)

derivative = diff(eqn)

eigenvalue1 = subs(derivative,y,theta)
eigenvalue2 = subs(derivative,K)
eigenvalue3 = subs(derivative,0)            

打印出的稳态为:0,theta,K。 因此,为了获得导数的特征值,我将这些值替换为y,但是我不确定这是否是正确的方法。特征值如下:

eigenvalue1 =
 
-r*(theta/K - 1) == 0
 
eigenvalue2 =
 
-r*(K/theta - 1) == 0
 
eigenvalue3 =
 
-r == 0 

这些对吗?如果是这样,我该如何回答它们是否稳定/不稳定的最后一个问题?

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。