MATLAB中方程的稳态和特征值

问题描述

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

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 (将#修改为@）