一类带负交叉扩散项的SIR传染病模型的空间Turing斑图

周 文,胡 伟,陈金琼,凯 歌

(1.安徽师范大学数学与统计学院,安徽芜湖 241002)

(2.北京工业大学机电学院,北京 100124)

1 引言

根据世界卫生组织的最新研究,传染病依旧是人类死亡的第一杀手.由于对传染病的研究不宜采用实验的形式,因此理论分析与数值模拟常被用于传染病的机理研究上.此时选择合适的传染病的动力学模型显得十分重要,常见的传染病模型有SI,SIR,SEIR等等.

在传染病的动力学研究中,许多学者做出有意义的结果,特别是传染病模型的空间斑图动力学[1−4].孙桂全,靳祯等人研究了一类带有时滞的SIR空间传染病模型的Turing失稳,通过数值模拟得出了条状与点状共存的斑图[1].王玮明等人研究了一类SI传染病模型,通过推导模型的振幅方程做出斑图选择,得到了不同类型的Turing斑图[2].马知恩和周义仓在研究中提出了这样的一个比例依赖型的模型[5]

由于R与S,I无关,所以作者利用稳定性理论讨论了下面一个微分方程组的平衡点的存在性与稳定性

考虑到疾病对易感者与感染者心理上的影响,本文在此基础上研究了一类带有负交叉扩散项的SI传染病模型

在这里S(t)和I(t)分别代表时间t时易感者和感染者的密度,γ和K分别代表了内禀增长率和环境的承载力,µ代表人口的自然死亡率,d代表了因病死亡率,c代表了染病者的恢复率,β是传染病系数,是传染率,α1和α2心理影响系数,并且α1和α2都是常数,∆代表空间Ω中的Laplace算子,Neumann边界条件表明模型(1.3)是自我封闭的,并且是零流量的,n代表着光滑边界∂Ω上的单位外法向量.

这里需要指出的是,在之前的文献[1–4]中,研究的传染病模型所带的交叉扩散系数都是正数,其生物学意义是人群总是从另一人群的高密度区域向低密度区域移动[6−7].而在实际生活中,人群从低密度区域向高密度区域移动的现象也是存在的.一方面,考虑到在生活中易感者有辨别染病者的能力并且会远离染病者,同时染病者也会远离易感者[6].另一方面,在疾病爆发初期,由于人们的心理因素,觉得人多的地方就是安全的地方,染病者反而会尽量接近易感者.所以在某种特定的情况下,这种染病者向易感者移动的现象在模型中则表现为交叉扩散系数为负数.据我们所知,带有负交叉扩散系数的传染病模型的Turing斑图在生物模型中很少被研究.因此本文将研究带有负交叉扩散系数的二维模型中的Turing斑图的生成问题.

2 Turing空间的确定

首先考虑模型(1.3)

易知系统(1.3)有很多平衡点,包括E0=(0,0),稳定节点E1=(K(1−,0),和正平衡点E∗=(S∗,I∗),从生物学上考虑,正平衡点更加有讨论的意义,其中

定理2.1不带有扩散项的系统(1.3)的平衡点(S∗,I∗)是局部渐近稳定的.

证 为了讨论方便,记

J的特征方程是

其中tr0=a11+a22,detJ=a11a22−a12a21.易知tr0<0和detJ>0当且仅当max{t1,t3}<β

通过Routh-Hurwitz,可知(S∗,I∗)是局部渐进稳定的.

现在考虑系统(1.3),并且对平衡点(S∗,I∗)进行线性化分析.如下面所示,在平衡点(S∗,I∗) 处作微扰:

这里pλ是在时间t上的扰动增长率;kx和ky是相应的振幅;i是虚数单位并且有i2=−1;k=是波数;S0和I0是两个正常数.把(2.2)式带入系统(1.3),并且省略所有的非线性项,可以得到特征方程

特征方程(2.3)的解为如下形式

其中

选择β作为分支参数.当Im(λk)6=0和Re(λk)=0在k=0时成立,系统出现Hopf分支,这样能得到Hopf分支曲线

利用稳定性定理[9−11],可知道当Im(λk)=0和Re(λk)=0在k=kT6=0时成立,Turing分支出现,且波数kT满足因此分支参数βT满足如下Turing分支曲线

根据Hopf和Turing分支曲线[12,13],能得到Hopf分支区域和Turing不稳定区域.

图1:模型(1.3)的分支图,其中d11=0.02,d22=5,d12=0.1,d21=−0.1

在图1中,可以看到系统(1.3)的分支图包含了Turing分支线和Hopf分支线,并且它们把γ−β参数空间分成了四个区域,区域D11被称为Turing空间,在这里发生Turing失稳,区域D12被称为Hopf空间,在这里发生Hopf失稳.

为了更好地理解参数对系统稳定性的影响作用,在图2,给出了随参数d21变化的色散关系图.线(3)对应着Turing临界值d21=−1.13,当d21=−0.9>−1.13时,Turing失稳发生;当d21=−1.2<−1.13时,Turing失稳消失,也就是说,此刻系统(1.3)的平衡解为稳定态.

图2:Re(λ)(特征值λ的实部)和k的关系,γ=0.2,K=1,µ=0.12,d=0.08,c=0.04,α1=0.4,α2=0.5,d11=0.02,d22=5,d12=0.1,β=0.1380和不同的d21:线(1):d21=−0.05;线(2):d21=−0.9;线(3):d21=−1.13;线(4):d21=−1.2

3 数值模拟

在这一部分,我们将通过Matlab对系统(1.3)进行一系列的数值模拟.所有的数值模拟均运用齐次Neumann边界条件.将空间区域离散为200×200个格子.对空间的离散采用有限差分法,设定空间步长为∆h=0.25,对时间的离散采用欧拉方法,取定时间步长为∆t=0.01.

首先设d11=0.02,d22=5,d12=0.1,γ=0.2,µ=0.12,β=0.1380,d=0.08,c=0.04,α1=0.4,α2=0.5.现在研究参数d21的不同值所产生的斑图.

图3:时间:(A)t=0;(B)t=20000;(C)t=40000;(D)t=800000

图4:时间:(A)t=0;(B)t=56000;(C)t=70000;(D)t=490000

在图3中,d21=−0.05,这时可见:(A)中颜色条数值基本不变,初值选取为平衡解加上一个随机扰动;(B)中出现类条状斑图;(C)中出现条状斑图;(D)条状斑图几乎占据了整个区域,且系统的动力学行为不再发生变化.

图4,图5分别是d21=−0.7和d21=−0.9时,染病者的时间演化图.由图4和图5可见:随着时间的演化,最终点状斑图和条状斑图共存.但图4中条状斑图占优;而当d21达到−0.9时,点状斑图会占优(图5(D)).由图6可见,当d21增至−1.1时,最终点状斑图几乎占满整个空间.

图5:时间:(A)t=0;(B)t=90000;(C)t=120000;(D)t=800000

图6:时间:(A)t=0;(B)t=520000;(C)t=600000;(D)t=800000

4 结论

本文研究了在Neumann边界条件下,负交叉扩散对带有非线性传染率的传染病模型的影响.具体表现为负交叉扩散可引起系统(1.3)在平衡点E∗处的Turing失稳,并由此得到了不同类型的斑图,包括点状斑图、条形斑图和点条混合斑图.