Matlab 寻找非线性常微分方程的多重解？

我试图编写一个MATLAB代码，使用非线性ODE的

fsolve

求解根。当我试图在稳态下找到这两个方程的解时，我只得到了一个解。有没有一种方法可以获得这个问题所有可能的稳态解？这就是问题所在，下面是我到目前为止的代码

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.5709395.4047]。所以我决定用wolfram alph来解x，我得到了三个单独的xss解。解决方案1:xss=[0.00475972508.669]解决方案2:xss=[0.5709395.4047]解决方案3:xss=[0.9973310.164]

你可能想问它

fsolve

只给你一个根，如果你想要另一个，使用不同的开始猜测你可能想问它

fsolve

只给你一个根，使用不同的起始猜测