捕食系统的Volterra方程

由于捕捞能力系数ε的引入， 食用鱼的平均量有了增加， 而食肉鱼的平均量却有所下 降，ε越大，平衡点的移动 也越大。
根据P-P模型，我们可以导出以下结论： （1）食用鱼的平均量取决于参数r1与λ1 （2）食用鱼繁殖率r1的减小将导致食肉鱼平均量的 减小，食肉鱼捕食能力λ1的增大也会使自己的平均量减 小；反之，食肉鱼死亡率r2的降低或食饵对食肉鱼供养 效率λ2的提高都将导致食用鱼平均量的减少。 （3）捕鱼对食用鱼有利而对食肉鱼不利，多捕鱼（ 当然要在一定限度内，如ε<r1）能使食用鱼的平均数量增 加而使食肉鱼的平均数量减少。 P-P模型导出的结果虽非绝对直理，但在一定程度上是附 合客观实际的，有着广泛的应用前景。例如，当农作物发生病 虫害时，不要随随便便地使用杀虫剂，因为杀虫剂在杀死害虫 的同时也可能杀死这些害虫的天敌，（害虫与其天敌构成一个 双种群捕食系统），这样一来，使用杀虫剂的结果会适得其反 ，害虫更加猖獗了。
max
，则 0 max
( x1 ) ( x1) 。同样根据的性质知，当 x1<x1< x1 时
( x1 ) 。此时： ( x2 )
由 ( x2 ) 的性质， x2、x2 ，使 ( x1 ) ( x2 ) S 成立。 max S max 当x1= x1 或 x1时， ( x1 ) ， ( x2 ) ( x1 ) ( x1 ) 0 x2 x2 时才能成立。 仅当 时，由于 ( x1 ) ， ( x2 ) S max max x1或x1> x1 而当x1< ( x1 ) ( x1 ) 故 ( x1 ) ( x2 ) S 无解。
dx1 x1 (r1 1 x2 ) 将两方程相除消去时间t，得： dx2 x2 (r2 2 x1 ) 分离变量并两边积分得轨线方程：
r ( x1r2 e2 x1 )( x21 e1x2 ) S
r ( x1 ) ( x1r e x ) ( x2 ) ( x2 e x ) 令
当 S max max时，轨线为一封闭曲线（图3-21），即周期解。
证明具有周期解。) ( x1
m
x 只需证明：存在两点 x1及 x1 ， 1 < x1 当 x1<x1< x1时，方程（3.32）有两 x0 时，方程恰 2 个解，当x1= x1 或x1= x1 a 时，方 有一解，而在x1< x1或x1> x1
m
图3-20 (a)
( x1 )
P0
x2
图3-21 图3-20 (b)
程无解。
x1
0
x10
r2 2
x11
0 0
x1
x20
1
r1
x1
0
x1
x2
x1
x21
事实上，若 S max max，记
S
0 0 ( x1 ) 的性质， x1、x1 ， 1 x1 而 x1 x1 ，使得： x 由
x2 r2 2 x1 x2
平衡点P的两个坐标恰为 2 图3-22 食用鱼与食肉鱼在一个 周期中的平均值。 x1 0 x1 0 将其在一个周期长度为T的区间上积分，得 x2 0 x2 0 t T x (t T )
x
ln
1
0
x1 (t0 )
r2T 2
2、模型分析
Po(0,0)是平凡平衡点且明 显是不稳定，没必要研究
方程组（3.31）是非线性的，不易直接求解。容易看 出，该方程组共有两个平衡点，即： r2 r1 P 0,0 和 P , 所以x1、x2轴是方程组的 1 0 2 1 两条相轨线。 方程组还有两组平凡解：
1、模型建立 Volterra将鱼划分为两类。一类为食用鱼（食饵），数 量记为x1(t)，另一类为食肉鱼（捕食者），数量记为x2(t)， 并建立双房室系统模型。
对于食饵（Prey）系统 ： 大海中有食用鱼生存的足够资源，可假设食用鱼独立生 存将按增长率为r1的指数律增长（Malthus模型），既设：
年代 1914 1915 1916 1917 1918 他知道，捕获的各种鱼的比例近似地反映了地中海里各 种鱼类的比例。战争期间捕鱼量大幅下降，但捕获量的下降 百分比 11.9 21.4 22.1 21.2 36.4 为什么会导致鲨鱼、鳐鱼等食肉鱼比例的上升，即对捕食者 年代 1919 1920 1921 1922 1923 有利而不是对食饵有利呢？他百思不得其解，无法解释这一 现象，就去求教当时著名的意大利数学家V.Volterra，希望 百分比 27.3 16.0 15.9 14.8 10.7 他能建立一个数学模型研究这一问题。
(r1 ) x1 1 x1 x2 x1 r1 x1 1 x1 x2 x1 食用鱼的数量反而 x2 (r2 ) x2 2 x2 r2 x2 2 x1 x2 因捕捞它而增加， x1 x2 真的是这样？！ 平衡点P的位置移动到了： r2 r1 P , 2 1
dx1 r1 x1 dt
由于捕食者的存在，食用鱼数量因而减少，设减少的速 率与两者数量的乘积成正比（竞争项的统计筹算律），即：
dx1 1 x1 x2 dt 出
λ 1反映了捕食者掠取食饵的能力
对于捕食者（Predator）系统 ： 捕食者设其离开食饵独立存在时的死亡率为r2，即：
问题背景：
§3.8 捕食系统的Volterra方程
意大利生物学家D’Ancona曾致力于鱼类种群相互制约 关系的研究，在研究过程中他无意中发现了一些第一次世 界大战期间地中海沿岸港口捕获的几种鱼类占捕获总量百 分比的资料，从这些资料中他发现各种软骨掠肉鱼，如鲨 鱼、鳐鱼等我们称之为捕食者（或食肉鱼）的一些不是很 理想的鱼类占总渔获量的百分比。在 1914~1923年期间，意 大利阜姆港收购的鱼中食肉鱼所占的比例有明显的增加：
2 2 1
1 1 2
（3.32）
用微积分知识容易证明：
两者应具有类似的性质
(0) () 0
x1
x1
2
2
r2
r2
'( x1 ) 0
有： max
x2
r2 ' 0 2 r x1 2 '( x1 ) 0
2
同理：对 ( x2 )
x1 (t ) x1 (0)e r1t x1 (t ) 0 和 x2 (t ) 0 x2 (t ) x2 (0)e r2t
当x1(0)、x2(0)均不为零时， 0 ，应有x1(t)>0且x2(t)>0， t 相应的相轨线应保持在第一象限中。
求（3.31）的相轨线
0
t0
x1 (t )dt
等式左端为零，故可得： r1 1 t0 T r2 1 t0 T x1 (t )dt 同理： t0 x2 (t )dt 1 T 2 T t0
x1 0 x2 0
x x1 0 1 x2 0
解释D’Ancona发现的现象 引入捕捞能力系数ε，（0<ε<1），ε表示单位时间 内捕捞起来的鱼占总量的百分比。故Volterra方程应为：
1
r1
有： max
( x1 ) 与 ( x2 ) 的图形见图3-20
易知仅当 S max max时（3.32）才有解 记： 0 x1
当 S max max时，轨线退化为平衡点。
2
r2
, x2
0
1
r1
0 0 讨论平衡点 ( x1 , x2 ) 的性态。
dx2 r2 x2 dt 出
但食饵提供了食物，使生命得以延续。这一结果也要通过竞 争来实现，再次利用统计筹算律，得到： 方程组（3.31）反映了在没有 dx1 人工捕获的自然环境中食饵 2 x1 x2 与捕食者之间的相互制约关 dt 入 系。下面我们来分析该方程 综合以上分析，建立P-P模型（Volterra方程）的方程组： 组。 x1 x1 (r1 1 x2 ) （3.31） x2 x2 ( r2 2 x1 )
max S max ( x1 ) ( x1 )
得证。
Байду номын сангаас
确定闭曲线的走向
用直线
l1 : x1 l2 : x2
2 将第一象限划分成四个子区域 1
r1
r2
在每一子区域，1与 x2不变号，据此确定轨线的走向（图3-22） x
将Volterra方程中的第二个改写成：