捕食者具有流行病的捕食扩散模型的分析_医学论文.docVIP

捕食者具有流行病的捕食扩散模型的分析_医学论文 捕食者具有流行病的捕食扩散模型的分析_医学论文 【摘要】 目的 建立并分析一类捕食者具有流行病的捕食扩散模型，重点研究该方程组解的定性性质。方法 利用线性化和特征值方法及Lyapunov稳定性理论。结果 得到了捕食者绝灭和疾病成为地方病的充分条件。结论 我们的结果表明，当捕食者的转化率大且接触率大时捕食者的疾病成为地方病；而当接触率小时该疾病最终消除。如果捕食者的接触率足够小时，连捕食者都最终消亡。 【关键词】 捕食模型 扩散 流行病 稳定性 Abstract： Objective To formulate and analyze a predator-prey diffusion model of predator with epidemic disease and to study its qualitative properties. Methods Linearization and eigenvalue and the theory of Lyapunov′s stability were used in the study. Result Predators′ extinction and disease were found to be the sufficient conditions of endemic disease. Conclusion Our results show that if the conversion efficiency and the contact rate of the predator are large, the predator′s disease becomes endemic. If the contact rate is small, the disease will be eliminated eventually. When the contact rate of the predators is sufficiently small, the predators themselves too tend to be extinguished. Key words： predator-prey model；diffusion；epidemic；stability 1 预备知识 种群动力学中捕食系统已被广泛地研究，目前已有一些工作开始考虑具有流行病的捕食模型。如文献［1-2］考虑了食饵染病的捕食-被捕食模型，文献［3］考虑了食饵有密度制约，疾病只在捕食者之间传播，而且染病的捕食者会因病死亡的生态-流行病模型，假设染病的捕食者不捕食饵，且假设这种传染病一旦染上不再康复，其模型为： dXdt=X(a-bX)-cXS dSdt=kcXS-d1S-βSI dIdt=βSI-d2I（1） 但是，只研究种群密度与时间之间的关系显然是不够的，应当重视种群在不同空间区域的扩散、迁移。为此，我们在系统(1)的基础上建立了一类捕食者具有流行病的捕食扩散模型： u1t-d1Δu1=u1(a1-bu1-pu2), (x,t)∈Ω×（0，+∞） u2t-d2Δu2=u2(-a2+kpu1-βu3),(x,t)∈Ω×（0，+∞） u3t-d3Δu3=u3(-a3+βu2),(x,t)∈Ω×（0，+∞） uiv=0,i=1,2,3,(x,t)∈Ω×（0，+∞） ui(x,0)=ηi(x),i=1,2,3,x∈Ω （2） 这里，u1表示食饵的密度，u2，u3分别表示易感类捕食者和染病类捕食者的密度，所有系数均为正常数，a1为食饵的内禀增长率，b为食饵的密度制约函数，p为捕食函数，k为转化函数，a2，a3分别表示易感类捕食者和染病类捕食者的死亡率。考虑到因病死亡因素，故a2＜a3，β为接触率，正常数di表示ui的扩散系数。Ω是Rn上的具有光滑边界的有界区域，齐次Neumann边界条件表示上述系统在边界上无种群迁移，初值ηi(x),(i=1,2,3)在Ω上非负且H连续，并在上满足相应的相容性条件（i=1，2，3）。本文主要研究系统（2）解的渐近行为，第2节给出解的局部渐近稳定性，第3节研究解的全局渐近稳定性。2 局部渐近稳定性 显然，系统(2)拥有平凡解E0(0,0,0)；半平凡解E1(ρ１,0,0),E2(ρ,bp(ρ1-ρ),0)和平衡解E3(ρ2,a3β,kpβ(ρ2-ρ))，其中ρ1=a1b，ρ2=ρ1-pa3bβ，ρ=a2kp。 既然密度是非负的，由初等知识易得 定理1 当ρ∈(0，ρ2)，E0，E1，E2，E3， 都出现，其中E3是系统（2）的惟一正平衡解；当ρ=ρ