问题描述
我正在尝试编写一个MATLAB代码以使用fsolve来求解非线性ODE的根。在试图找到稳态下这两个方程的解时，我只得到一个解。是否有一种方法可以解决此问题，并获得所有可能的稳态解？这就是问题所在，下面是我到目前为止拥有的代码。
ODE dx / dt = [u1-x1-72x10 ^ 9exp（-10000 / x2）x1; 310-x2 + 1.44x10 ^ 13exp（-10000 / x2）x1 + u2 / 23.9]
给定稳态uss = [0.998,0.5]，在稳态xss = [x1，x2]处找到x
function F = root2d(x)
F(1) = 0.998 - x(1) - (((72*10^9).*exp(-10000/x(2))).*x(1));
F(2) = 310 - x(2) + (((1.44*10^13).*exp(-10000/x(2))).*x(1)) + 0.5/23.9;
fun=@root2d;
xo=[0,0];
x=fsolve(fun,xo)
x =
0.9973 310.1638
但是，我看到解决方案应该是xss = [0.5709,395.4047]。因此，我决定使用Wolfram Alph求解x，并获得了三个单独的xss解。解决方案1：xss = [0.00475972,508.669]解决方案2：xss = [0.5709,395.4047]解决方案3：xss = [0.9973,310.164]
解决方法
暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！
如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@）