常微分方程建模R包ecode（三）——探寻平衡点

在建立常微分方程后，我们常常会问的一个问题是： 这个常微分方程系统中是否存在平衡点？如果存在，这个平衡点是否稳定？ 本节将展示如何通过ecode包来回答这个问题。

我们首先利用ecode包建立一个Lotka–Volterra竞争模型（详见第一节），
$$\begin{gathered} \frac{dx}{dt}=(r_1-a_{11}x-a_{12}y)x,(1) \\ \frac{dy}{dt}=(r_2-a_{21}x-a_{22}y)y,(2) \end{gathered}$$
其中，$x$代表物种1的种群个体数，$x$代表物种2的种群个体数，$r_1,r_2$为种群增长率，$a_{11},a_{12},a_{21},a_{22}$为两物种之间的竞争系数。

library(ecode)

eq1 <- function(x, y, r1 = 1, a11 = 1, a12 = 2) (r1 - a11 * x - a12 * y) * x
eq2 <- function(x, y, r2 = 1, a21 = 2, a22 = 1) (r2 - a21 * x - a22 * y) * y
x <- eode(dxdt = eq1, dydt = eq2)

一、数学原理

对任意一个常微分方程组，
$$\begin{gathered} \frac{d{x}_1}{dt}=f_1(x_1,x_2,...,x_n), \\ \frac{d{x}_2}{dt}=f_2(x_1,x \end{gathered}$$