微分方程组建模
在本节中,我们讨论两个非常不同的现象的模型-捕食者-猎物系统中两个种群的进化以及质量弹簧系统的运动。最初,这些模型似乎完全不同,但是从正确的角度来看,它们具有许多相似之处。
捕食者-猎物系统的再探
考虑两种捕食者-猎物模型的微分方程组,第1.1节中简要讨论过。回顾一下, R ( t ) R(t) R(t) 表示时间 t t t 时猎物的种群(以千为单位,或以百万为单位,或其他单位), F ( t ) F(t) F(t) 表示捕食者的种群。我们假设 R ( t ) R(t) R(t) 和 F ( t ) F(t) F(t) 都是非负的。可能描述这两个物种种群变化的微分方程组是:
d R d t = 2 R − 1.2 R F \frac{dR}{dt} = 2R - 1.2RF dtdR=2R−1.2RF
d F d t = − F + 0.9 R F \frac{dF}{dt} = -F + 0.9RF dtdF=−F+0.9RF
在关于 d R d t \frac{dR}{dt} dtdR 的方程中, 2 R 2R 2R 项表示在没有捕食者的情况下猎物的指数增长,而 − 1.2 R F -1.2RF −1.2RF 项对应于捕食者与猎物相互作用对猎物的负面影响。在关于 d F d t \frac{dF}{dt} dtdF 的方程中, − F -F −F 项对应于假设捕食者在没有猎物的情况下会死亡,而 0.9 R F 0.9RF 0.9RF 项对应于捕食者与猎物相互作用对捕食者的正面影响。系数 2 、 − 1.2 、 − 1 2、-1.2、-1 2、−1.2、−1 和 0.9 0.9 0.9 依赖于所涉及的物种。不同系数的类似系统在练习中有所考虑。(我们选择这些参数值仅仅是为了方便。)方程中的 R F RF RF 项使得系统难以求解。我们无法推导出一般解的显式公式,但有一些初始条件可以得到简单的解。例如,假设 R = 0 R = 0 R=0 和 F = 0 F = 0 F=0。那么两个方程的右侧在所有 t t t 时刻都为零( d R d t = d F d t = 0 \frac{dR}{dt} = \frac{dF}{dt} = 0 dtdR=dtdF=0),因此常数函数 R ( t ) = 0 R(t) = 0 R(t)=0 和 F ( t ) = 0 F(t) = 0 F(t)=0 形成了该系统的解。类似一阶方程,我们将这样的常数函数对称为系统的平衡解。这一对平衡解完全有意义:如果捕食者和猎物的种群都消失了,我们当然不期望这些种群在之后会增长。
下面尝试寻求方程组的其他常数解。我们将系统重新表示为
d R d t = ( 2 − 1.2 F ) R \frac{dR}{dt} = (2 - 1.2F)R dt