微分方程 · 捕食者-猎物模型

捕食者-猎物模型

捕食者-猎物模型(洛特卡–沃尔泰拉方程)描述猎物和捕食者两个种群:猎物自己繁殖、被吃掉,捕食者靠吃猎物增长、否则就减少。两个数量绕着一个平衡点周期性地起落。

捕食者-猎物模型与本页画面

这里xx是猎物的数量,yy是捕食者的数量,单位都是千只。两者都按“每只每年”的比率变化:猎物以α\alpha繁殖,被吃掉的部分与捕食者数量成正比;捕食者以γ\gamma死亡,增长与猎物数量成正比。相平面上画的是点(x,y)(x, y)从起点出发走过的路径,下方的时间图逐年画出两个数量。

不捕捞、不设上限时的模型
x′=x(α−βy),y′=y(δx−γ)x' = x(\alpha - \beta y), \qquad y' = y(\delta x - \gamma)

猎物的增长率α−βy\alpha - \beta y在直线y=α/βy = \alpha/\beta上为零,捕食者的增长率δx−γ\delta x - \gamma在直线x=γ/δx = \gamma/\delta上为零。这两条虚线就是零增长线,它们交于平衡点(x∗,y∗)=(γ/δ,α/β)(x^*, y^*) = (\gamma/\delta, \alpha/\beta),把平面分成四块。在y∗y^*下方猎物增加,在上方猎物减少;在x∗x^*右侧捕食者增加,在左侧捕食者减少。所以路径沿逆时针方向绕圈。捕食者在停止增长的那一刻达到高峰,也就是猎物降回x∗x^*的时候,这比猎物自己的高峰晚。读数“捕食者的高峰晚”显示的就是这段时间。

绕圈本身并不保证回到原处,螺旋也在绕圈。让路径闭合的是一个沿每条解都不变的量:V=δx−γln⁡x+βy−αln⁡y,V′=(δx−γ)(α−βy)+(βy−α)(δx−γ)=0.\begin{aligned} V &= \delta x - \gamma \ln x + \beta y - \alpha \ln y, \\ V' &= (\delta x - \gamma)(\alpha - \beta y) + (\beta y - \alpha)(\delta x - \gamma) = 0. \end{aligned}VV高于最小值的每条等高线都是围着平衡点的闭合曲线,路径沿它回到起点,周而复始。加上上限KK后,拥挤会让猎物长得更慢,VV不再守恒;在上限那个精选里,路径不再闭合,而是一圈圈收向平衡点。

不管起伏多大,每一整圈的平均都正好是平衡点。把捕食者的方程除以yy,得到(ln⁡y)′=δx−γ(\ln y)' = \delta x - \gamma。一整圈后yy回到原值,左边积分为0,所以xx的平均是γ/δ=x∗\gamma/\delta = x^*。同样由猎物的方程得到yy的平均是α/β=y∗\alpha/\beta = y^*。起伏变大,改变的是最高最低和周期,不是平均。

两个种群都按同样的每只比率ε\varepsilon持续被捕捞,于是α\alpha变成α−ε\alpha - \varepsilon,γ\gamma变成γ+ε\gamma + \varepsilon。不设上限时,平衡点移到x∗=(γ+ε)/δx^* = (\gamma + \varepsilon)/\delta、y∗=(α−ε)/βy^* = (\alpha - \varepsilon)/\beta:平均下来猎物更多,捕食者更少。猎物的平均数正是让捕食者刚好收支相抵的数量,捕捞让捕食者死得更快,这个数量就随之升高。

值得注意

  • 在方程里,种群可以非常接近0,但永远到不了0。真实的种群少到这种程度可能就灭绝了。
  • 系数和单位只是示意,不对应任何真实物种。

操作

  • 起点不能放在坐标轴上:从0开始的种群会一直是0。

相关概念