微分方程（Blanchard Differential Equations 4th）中文版Section5.5

三维空间上的非线性系统

我们在第2.8节中看到，具有三个依赖变量的微分方程的解是三维空间中的曲线。这些曲线可以以非常复杂的方式相互环绕。在第3.8节中，我们研究了具有三个依赖变量的线性系统的行为。线性系统的行为可以通过特征值和特征向量来确定。然而，相较于平面系统，可能的行为列表要长得多。

在本节中，我们考虑两个三维非线性系统的例子。第一个例子是一个包含三个物种的食物链的种群模型。第二个例子，我们回到第2.8节讨论的洛伦兹方程。我们可以使用平衡点线性化和数值逼近（如欧拉方法）来获得关于这些系统的一些信息。正如我们在前面几节中强调的，对于三维及更高维的微分方程系统，几乎没有通用的工具。事实上，这些系统是数学研究中的活跃领域。

食物链模型

我们已经研究了在孤立环境中生活的种群以及在捕食者-猎物、合作和竞争系统中相互作用的物种对。这些系统仅仅触及了自然界中复杂互动的表面。一个可能的情况是形成一个由三个或更多物种组成的食物链。一个物种可以既是捕食者又是猎物。最近研究的一个例子涉及到冷杉树、驼鹿和狼。冷杉树被驼鹿吃掉，驼鹿（尤其是年轻和虚弱的个体）则被狼吃掉。一个自然的问题是，狼的种群变化是否会影响树木的种群。最近对位于罗瓦尔岛国家公园的树木/驼鹿/狼种群的研究表明，狼的种群变化可能会影响树木的种群。

我们构建一个包含三个物种的系统模型，每个物种都被下一个物种吃掉。为了方便起见，我们将这些物种称为树木、驼鹿和狼，与上面的例子类比。设定：

• $x(t)$ = 时间$t$时树木的种群，
• $y(t)$ = 时间$t$时驼鹿的种群，
• $z(t)$ = 时间$t$时狼的种群。

我们假设每个种群在孤立状态下可以用逻辑斯蒂增长方程来建模，而物种之间的互动效果与种群的乘积成正比。解的行为将取决于为增长率、承载能力和互动效应选择的参数。为了了解这种形式系统的解的行为，我们首先将所有参数设为1。（这对于树木、驼鹿和狼并非实际情况，但使接下来的计算更加简单。）我们的模型为：

$$\frac{dx}{dt}=x(1-x)-xy$$

$$\frac{dy}{dt}=y(1-y)+xy-yz$$

$$\frac{dz}{dt}=z(1-z)+yz$$

食物链模型

注意到，树木种群的增长受到驼鹿存在的抑制（即 $-xy$ 项），驼鹿种群的增长则受到树木的促进（即 $+xy$ 项）但受到狼的抑制（即 $-yz$ 项），而狼的种群增长受到驼鹿的促进（即 $+yz$ 项）。

我们可以通过将方程的右侧设为零并求解 $x$、$y$ 和 $z$ 来找到该系统的平衡点。平衡点为 $(0,0,0)$、$(1,0,0)$、$(0,1,0)$、$(0,0,1)$、$(1,0,1)$ 和 $(2/3,1/3,4/3)$。其中，只有点 $(2/3,1/3,4/3)$ 的三个坐标都不为零，因此三个物种可以在这些种群数量下共存于平衡状态。该系统在 $(x,y,z)$ 处的雅可比矩阵为：

$$\left(\begin{array}{ccc}1-2x-y& -x& 0\\ y& 1-2y+x-z& -y\\ 0& z& 1-2z+y\end{array}\right).$$

因此，在平衡点 $(2/3,1/3,4/3)$ 处，雅可比矩阵为：

( − 2 / 3 − 2 / 3 0 1 / 3 − 1 / 3 − 1 / 3 0 4 / 3 − 4 / 3 ) . \begin{pmatrix} -2/3 & -2/3 & 0 \\ 1/3 & -1/3 & -1/3 \\ 0 & 4/3 & -4/3 \end{pmatrix}.