改进Leslie–Gower比率依赖反应扩散捕食–猎物系统的扩散驱动失稳与Hopf分岔分析

《Journal of Low Frequency Noise, Vibration and Active Control》：Existence of diffusion-driven instability and Hopf bifurcation in a Spatio-temporal Predator–Prey system

　　

本研究对一类比率依赖（ratio-dependent）反应扩散（reaction–diffusion, RD）捕食–猎物模型展开系统动力学分析，重点考察非线性振荡行为与时空调制模式演化。研究人员以偏微分方程组为对象，通过局部与全局稳定性分析刻画平衡态对小扰动的响应，并建立长期动态稳定条件。研究特别关注Hopf分岔（Hopf bifurcation, HB）分析，辨识驱动自持振荡与周期动力学产生的临界参数。进一步，研究人员考察扩散驱动不稳定性（diffusion-driven instability, DDI）机制，揭示反应动力学与扩散耦合所诱发的空间分布振荡结构。数值模拟用于验证理论结果，并展示稳定态、周期态与空间结构化动态区间之间的转变。针对非标准有限差分（nonstandard finite difference, NSFD）格式，研究人员结合放大矩阵分析、网格加密检验并与经典显式有限差分方法对比，评估离散稳定性、一致性、保正性（positivity preservation）与小步长精度。结果表明，该模型能够刻画由非线性相互作用主导的振荡响应与复杂时空行为，为反应扩散捕食系统的稳定性、分岔与图灵型空间模式形成提供了统一数学框架。

研究背景与问题提出

在数学生态学中，捕食–猎物（predator–prey, PP）关系关乎生物多样性与生态系统稳定。早期Lotka–Volterra模型假设捕食率仅依赖猎物密度，忽视捕食者干扰与空间异质性。比率依赖（ratio-dependent, RD）机制在猎物密集时更符合实际，将捕食功能反应写为猎物与捕食者密度之比的形式。引入反应扩散（reaction–diffusion）项后，即便无扩散时稳定的均匀平衡态，也可能因物种扩散速率差异而产生扩散驱动不稳定性（diffusion-driven instability, DDI），形成图灵斑图（Turing patterns）。已有文献多聚焦单一机制，如时滞、阶段结构、交叉扩散或Turing–Hopf动力学，缺乏将时间稳定性、全局行为、Hopf分岔（Hopf bifurcation, HB）、DDI与数值保正性统一分析的框架。研究人员据此建立改进Leslie–Gower型RD捕食–猎物模型，系统回答“参数如何使系统从稳定共存走向振荡与空间结构化”这一问题。

研究内容与总体结论

研究人员构建带No-flux边界的一维RD系统：猎物种群p(x,t)具logistic自限增长与Holling II型比率依赖捕食项，捕食者q(x,t)依Leslie–Gower项增长；两物种均含扩散系数??、??。首先证明正解存在性（positivity），继而分析平衡点（灭绝态、仅猎物态、仅捕食者态及共存态F?(p?,q?)），以捕食转化效率s为关键分岔参数。局部与全局稳定性分析显示：s较小时系统收敛于共存平衡；s跨越临界值后出现HB，诱发周期振荡；进一步在扩散项中因??≠??产生DDI，使均匀态失稳并形成空间斑图。数值上设计保正NSFD格式，经放大矩阵（amplification-matrix）分析、网格细化测试与显式格式对比，确认离散稳定性、一致性与小步长精度。研究意义在于为非线性分布动力系统提供可迁移的稳定性—分岔—模式形成分析范式。

主要关键技术方法

研究人员采用五类方法：一，基于Jacobian（雅可比矩阵）的线性化与特征分析，判定平衡态局部稳定性与HB临界条件；二，正锥不变性与最大值原理证明解的非负性；三，以s为分岔参数的Hopf与DDI解析条件推导；四，非标准有限差分（NSFD）离散与放大矩阵稳定性分析；五，数值模拟与网格加密检验，对比经典显式差分，验证保正性与精度（无生物试剂与培养操作，属纯数学建模与计算实验）。

研究结果

1. 2.
   时空模型（Spatio-temporal Model）
   给出含扩散项的改进Leslie–Gower RD方程组及初边值条件，定义参数生物学含义：β为猎物种内竞争，m为捕食率，e?、e?为半饱和常数，??、??为迁移率，s为捕食者能量转化效率。通过引理与定理证明：在严格正初值下，p(x,t)、q(x,t)恒正。
2. 3.
   稳定性评估与Hopf分岔（Stability Assessment and Hopf-Bifurcation）
   3.1 平衡态分析（Analysis of Equilibrium State）
   枚举原点、猎物单独存活、捕食者单独存活及两类共存平衡；指出仅当F?(p?,q?)双分量正时具生态可行性，F?因非正被剔除。
   3.2 PDE局部稳定性（Local Stability of PDE）
   构造紧凑形式 ?z/?t = NΔz + H(z)，在F?处求Jacobian J，得h??、h??、h??=s、h??=?s。由Trace与Determinant符号变化识别HB：s增大使特征值穿越虚轴，均匀态由稳定变周期。
   3.3 扩散驱动不稳定性（Diffusion-driven Instability）
   在均匀平衡上叠加空间模态，分析扩散修正特征根。表格参数显示：s=0.28时Z?>0不触发DDI；s=0.31–0.53时Z?<0且扩散项放大扰动，系统进入图灵型不稳定；s=0.62又回归稳定，说明s在区间内存在“失稳窗口”。

数值与讨论部分总结

研究人员将连续模型离散为NSFD格式，放大矩阵分析表明该格式在合理步长下离散稳定且保持正性；网格细化结果收敛，与显式方法相比避免非负性破坏。讨论指出：模型虽源于生态学，但其“反应—输运”耦合、HB与DDI机制可迁移至分布式工程系统，如反应输运网络中的稳定性丧失与空间模态选择。统一框架的价值不在直接等同生态设备，而在提供可复用的非线性分布系统分析工具。

结论翻译

研究人员证实：改进Leslie–Gower比率依赖反应扩散模型在合适参数下，依次呈现稳定共存、Hopf诱导周期振荡与扩散驱动空间斑图三类体制；正解性可被严格保证，NSFD数值格式保真且保正；以s为控制参数可解析刻画从时间振荡到空间结构化的转变，该框架对非线性分布动力系统的稳定性与分岔研究具有通用参考价值。

论文解读已完成，接下来可以选择一个方向继续深化：

把NSFD格式与放大矩阵分析写成推导小节

对比本文与Turing–Hopf文献的方法差异

把NSFD格式与放大矩阵分析写成推导小节

对比本文与Turing–Hopf文献的方法差异