捕食者-猎物系统中可变环境承载力与Cosner功能响应的时空动态研究

《Franklin Open》：Spatio-temporal dynamics of a predator–prey system with variable carrying capacity and Cosner response function

研究人员针对生态系统中环境承载力(k)高度不可预测的问题，开展了捕食者-猎物系统动力学研究。通过建立包含Logistic增长率和Cosner功能响应的数学模型，分析了平衡点的稳定性条件，发现了Hopf分叉和鞍结分叉等非线性现象。研究揭示了捕食者处理时间(d1)和种群内竞争强度(d2)对系统稳定性的调控机制，为理解生态系统的空间格局形成提供了新见解。

　　

在生态系统中，捕食者与猎物的相互作用一直是研究的核心问题。传统的Lotka-Volterra模型虽然奠定了理论基础，但无法解释复杂多变的自然现象。特别是环境承载力(k)这一关键参数，传统模型通常假定其为常数，而实际上它会受到种群活动的影响而动态变化。此外，捕食者的功能响应形式也直接影响系统动力学，常见的Holling类型功能响应未能充分考虑捕食者的集群行为。这些问题限制了我们对真实生态系统的理解和预测能力。

针对这些挑战，研究人员开展了一项创新性研究，建立了包含可变环境承载力和Cosner功能响应的捕食者-猎物模型。该系统考虑了猎物种群对其自身环境承载力的影响，通过引入参数β来表征这种影响的方向和强度（β>0表示建设性影响，β<0表示破坏性影响）。同时采用Cosner等人提出的功能响应形式，更好地描述了捕食者集群觅食的行为特征。

研究采用了非线性动力学分析和数值模拟相结合的方法。通过稳定性分析和分叉理论，解析了系统平衡点的存在条件及其稳定性特征。运用相空间分析和数值计算，揭示了系统在不同参数条件下的动态行为。空间扩展模型通过扩散项考虑了种群的空间分布，采用图灵不稳定性理论分析了空间格局的形成机制。

研究首先建立了基础模型：
du/dt = ru(1-u/(k+βu)) - αuv²/(1+αhuv)
dv/dt = θαuv²/(1+αhuv) - d₁v - d₂v²

在"系统平衡点及其稳定性"部分，研究发现系统存在三种平衡点：灭绝平衡点E₀=(0,0)、无捕食者平衡点E₁=(u₁,0)和共存平衡点E*=(u*,v*)。通过Jacobian矩阵分析，确定了各平衡点的稳定性条件。特别发现当β变化时，系统会经历从稳定焦点到极限环的转变。

"分叉分析"部分揭示了丰富的非线性现象。Hopf分叉发生在迹条件Tr(J)=0而Det(J)>0时，导致稳定焦点失稳并产生周期振荡。鞍结分叉则出现在两个平衡点碰撞消失的情况下。数值模拟展示了当β=-0.35时系统无共存平衡点，而β=-0.3103828时出现唯一平衡点，β=0.05时则存在两个平衡点的有趣现象。

在"Bogdanov-Takens分叉"部分，研究发现当Jacobian矩阵出现二重零特征值时，系统会经历这种高阶分叉。数值计算确定了分叉点的精确位置，展示了分叉曲线如何将参数平面划分为不同动力学区域。

"Bautin分叉（广义Hopf）"部分探讨了Lyapunov系数l₁变号的情况。这种分叉标志着系统从超临界Hopf分叉（产生稳定极限环）向亚临界Hopf分叉（产生不稳定极限环）的转变，对理解系统的全局动力学至关重要。

研究结论指出，可变环境承载力和Cosner功能响应的引入显著丰富了捕食者-猎物系统的动力学行为。β参数的变化可以导致系统稳定性的根本改变，而捕食者处理时间h的增加则可能抑制振荡、促进系统稳定。这些发现为理解自然生态系统的复杂动态提供了新视角，特别是在解释种群空间分布格局的形成机制方面具有重要意义。

这项工作发表在《Franklin Open》上，不仅推动了理论生态学的发展，也为生态系统管理和物种保护提供了科学依据。未来研究可以进一步考虑环境噪声、时滞效应等因素，使模型更贴近真实的生态系统。