中文标题 分数阶模型及两株系SPFMV传播动力学分析：O株系与EA株系案例

《Scientific African》：Fractional-order model and analysis of a two-strain SPFMV transmission dynamics: A case of strain O and EA

　　

摘要翻译
本研究构建了甘薯羽状斑驳病毒(SPFMV)传播动力学的数学模型，重点关注两种株系：O株系与EA株系。模型采用分数阶导数，以刻画健康甘薯与感染甘薯及蚜虫媒介之间的复杂相互作用。将两株系框架与具有记忆依赖特征的动力学相结合，相较于经典整数阶模型可提供更高的真实性。关键结果表明，当分数阶(α)趋近于1时，健康甘薯与蚜虫的下降速度加快，从而导致感染显著增加。具体而言，在α=0.4时，20%的健康甘薯和10%的蚜虫因感染而损失；而在α=1时，损失分别升至100%和65%。EA株系的峰值感染率明显高于O株系，甘薯中的差异分别为78.6%和70%，蚜虫中的差异分别为83.8%和70%。研究结果强调，政策制定者应优先推广对两种株系均具有抗性的甘薯品种，并实施混作策略以减少蚜虫危害，从而提升甘薯可持续生产与粮食安全。

论文解读
该文发表于《Scientific African》，围绕甘薯羽状斑驳病毒(SPFMV)的传播过程展开建模分析。甘薯是重要粮食作物，但其生产常受病毒病严重制约，尤其SPFMV在全球分布广、经济危害大。已有研究表明，该病毒在东非及乌干达地区流行明显，且O株系与EA株系是主要致病株系之一。传统整数阶模型通常只考虑当前状态，难以反映植物—病毒—媒介系统中普遍存在的记忆效应与历史依赖。基于此，研究人员引入Caputo–Fabrizio分数阶导数(CF)，构建包含O与EA两株系及蚜虫媒介的两株系传播模型，以提升对SPFMV复杂动力学的刻画能力。研究不仅求得无病平衡点与地方流行平衡点，还分析了基本再生数(R₀)、平衡点稳定性及参数敏感性，最终指出EA株系的传播优势更强，且分数阶参数对疫情演化具有显著调控作用。这一工作对于理解甘薯病毒病传播规律、优化防控策略及提升粮食安全具有重要意义。

研究方法上，作者首先建立了由甘薯和蚜虫两类宿主组成的确定性分室模型，并以Caputo–Fabrizio分数阶框架将其推广至分数阶系统。随后采用下一代矩阵法计算基本再生数(R₀)，利用线性化分析与Lyapunov函数讨论无病平衡和地方病平衡的稳定性，并通过固定点理论证明分数阶模型解的存在唯一性。最后，使用Latin Hypercube Sampling(LHS)与Partial Rank Correlation Coefficient(PRCC)进行敏感性分析，并借助MATLAB数值模拟验证理论结果。

### 研究结果
**1. 模型构建与正性/有界性分析**
通过建立含健康甘薯、O株系感染甘薯、EA株系感染甘薯，以及相应健康蚜虫和两类感染蚜虫的六维模型，研究人员证明了在非负初值下各状态变量始终保持非负，并给出了种群总量的正不变区域，说明模型具有生物学合理性。

**2. 无病平衡点与基本再生数**
研究人员推导出无病平衡点E₀，并得到两类株系对应的阈值量R₀₁与R₀₂，总体基本再生数为R₀=max{R₀₁,R₀₂}。结果表明，当R₀₁<1且R₀₂<1时，无病平衡点局部渐近稳定；否则系统将发生流行。

**3. 地方病平衡点与稳定性**
当仅存在EA株系或仅存在O株系时，模型分别给出两个地方病平衡点E₁与E₂。利用Lyapunov函数分析表明，这些地方平衡点满足局部渐近稳定条件，说明病毒在相应条件下可持续存在于种群中。

**4. 分数阶模型的存在唯一性**
在Caputo–Fabrizio分数阶框架下，研究人员证明了系统解的存在性与唯一性，说明所建分数阶模型在数学上是良定的，可用于后续数值计算与动力学分析。

**5. 参数敏感性分析**
LHS/PRCC结果显示，甘薯种植率(Λ_p)、蚜虫补充率(Λ_a)以及各类感染率β₁、β₂、β_a1、β_a2对R₀具有显著正向影响；自然死亡率μ与病毒损失率τ₁、τ₂则表现为负相关。表明降低传播与补充参数、提高死亡或失毒速率，有助于控制SPFMV。

### 讨论与结论
数值模拟进一步表明，分数阶参数α对疾病进程影响显著。随着α接近1，系统行为更接近经典整数阶模型，健康甘薯和蚜虫下降更快，感染增长更迅速；当α较低时，记忆效应增强，可减缓传播。文中给出，α=0.4时健康甘薯和蚜虫分别损失约20%和10%，而α=1时分别损失约100%和65%。这说明历史状态对当前传播具有重要影响，分数阶建模能够更真实地反映植物病毒系统的非局部动力学。

研究还发现，EA株系在甘薯与蚜虫中的峰值感染均高于O株系，甘薯侧差异可达78.6%与70%，蚜虫侧可达83.8%与70%。该结果表明EA株系具有更强传播潜力或适应性，防控上应优先针对EA株系开展抗性品种选育和综合治理。讨论部分指出，模型虽能揭示记忆效应、阈值行为和株系优势，但仍受限于参数数据不足、空间异质性未纳入以及未考虑共感染等简化假设。

总之，研究人员通过构建并分析两株系SPFMV分数阶传播模型，证明了分数阶动力学可显著影响病毒传播速度与感染峰值，且EA株系为更强势传播株系。该研究为甘薯病毒病的阈值识别、媒介防控和抗病品种推广提供了定量依据，对可持续农业与粮食安全具有现实指导意义。