張敏華
(陽光學(xué)院基礎(chǔ)教研部 福建福州 350015)
在客觀世界里普遍且大量存在時滯反應(yīng)現(xiàn)象與非局部空間作用,這可以構(gòu)建非局部時滯反應(yīng)擴散方程來進行描述。非局部時滯反應(yīng)擴散方程結(jié)合了時間與空間上的非局部現(xiàn)象,被廣泛應(yīng)用于生態(tài)學(xué)、傳染病學(xué)等多個領(lǐng)域,也成為了高等數(shù)學(xué)應(yīng)用研究的重點與難點內(nèi)容,相較于非局部反應(yīng)擴散方程更為復(fù)雜[1]。形如u(x,t)=u(x+ct)的行波解是非局部時滯反應(yīng)擴散方程一類重要的穩(wěn)態(tài)解,行波解能揭示非局部時滯反應(yīng)擴散方程重要性質(zhì),且由于其具有空間平移不變性的特征能有效表現(xiàn)與描述自然界內(nèi)的振動與有限速度傳播現(xiàn)象,如可以表現(xiàn)生物物種的入侵過程以及傳染疾病源的傳播規(guī)律等,可見非局部時滯反應(yīng)擴散方程的行波解具有很強的應(yīng)用實踐價值。目前伴隨著非局部時滯反應(yīng)方程應(yīng)用強度進一步深入,為了解非局部時滯反應(yīng)方程的重要性質(zhì)以及解決實際應(yīng)用問題,非局部時滯反應(yīng)方程的行波解展開研究是非常必要的。對此,在反應(yīng)擴散方程理論的基礎(chǔ)上,著重分析了具有時滯的反應(yīng)擴散的SIR 模型以及非擬單調(diào)時滯非局部擴散系統(tǒng)的行波解存在性。
自然界中擴散中普遍存在時間滯后和空間非局部作用,為進一步解決與描述這一現(xiàn)象,在上個世紀(jì)70年代起研究者將時滯與局部擴散相結(jié)合,如Levin將擴散反應(yīng)引入到單種群模型中得到了具有時滯擴散Logistic 模型[2]。但是隨著研究的深入發(fā)現(xiàn)提出的局部有限時滯模型具有明顯的局限性,模型中時滯與空間擴散是彼此獨立的。但是在實際現(xiàn)象中,在過去的時間里個體的空間內(nèi)位置并不是一成不變的,而是隨著時間變化個體在空間內(nèi)是不斷發(fā)生變化。針對這一問題Britton(1989)進行了全方位的考慮,利用空間加權(quán)平均的思想提出如下單物種種群模型,其中g(shù)是指定函數(shù),卷積項(g*u)等于
后來研究學(xué)者將這類具有結(jié)合了時空加權(quán)平均時滯項的反應(yīng)擴散方程統(tǒng)稱為非局部時滯反應(yīng)擴散方程。隨后Smith和Thime(1991)通過采用對具有年齡結(jié)構(gòu)的種群模型沿著種群特征積分的方法推導(dǎo)了一個具有成熟結(jié)構(gòu)的種群模型。隨后這種方法被So(2001)、Weng(2003)等人有效推廣至連續(xù)空間以及具有無限多個離散的斑塊環(huán)境下非局部時滯反應(yīng)擴散方程。此后,非局部時滯反應(yīng)擴散方程被應(yīng)用于多個領(lǐng)域,研究也愈發(fā)深入,包括Al-Omari與Gourley.S(2003)構(gòu)建的Loth-Volterra型兩物種競爭模型,Lian(2003)建立的具有全局交互效應(yīng)時滯反應(yīng)對流擴散模型[3-4]。
在非局部時滯反應(yīng)擴散方程行波解研究上,Xiao與Ruan(2004)通過線性鏈技術(shù)與幾何擾動理論對滿足的一類傳染病模型在且b大于a的條件下模型行波解存在性進行研究,研究結(jié)果顯示波速不小于時存在連接模型兩個平衡點的行波解5];J.Zhang(2007)同樣采用以上技術(shù)與理論方法對方程行波解的存在性進行了研究[6]。Xu和Zhao(2005)采用相平面以及譜理論等技術(shù)方法對一類傳染病模型的雙穩(wěn)波前解的存在性與全局漸近穩(wěn)定性等進行了研究。
對非局部時滯反應(yīng)擴散方程行波解的研究具有很強的實際應(yīng)用價值,近年來已成為研究的熱點問題。國內(nèi)眾多學(xué)者對這類方程的行波解也持以了很高的關(guān)注度,對行波解的存在性、穩(wěn)定性以及漸進性等展開了豐富的研究。馬萬彪等人(2013)通過采用慣性流行理論對具有小離散時滯反應(yīng)擴散方程的行波解的存在性進行了證明,結(jié)果顯示該方法能有效避免反應(yīng)項的假設(shè)設(shè)定問題[7];張衛(wèi)國等人(2016)對一類具有時滯Lotka-Volterra系統(tǒng)的行波解存在性通過簡化的方式尋找競爭系統(tǒng)模型的上下解[8];鄒霞(2018)則針對一類具有時空時滯的傳染病SIR模型,通過Schauder不動點定理以及求極限等方法來判斷行波解的存在性[9];李成林(2017)對具有Neumann邊界條件擴散反應(yīng)的三種物種時滯系統(tǒng)采用特征值法對系統(tǒng)的穩(wěn)定性進行了分析,并通過尋找上下解的方式證明在波速c足夠大時存在連接零、正兩個平衡點的行波解[10]。
對非局部時滯反應(yīng)擴散方程行波解研究是當(dāng)前的重點與難點,取得了一定的研究進展。但是縱觀目前的研究結(jié)果來看,國內(nèi)對于非局部時滯反應(yīng)擴散方程行波解的研究還處于起步階段,對此針對這一現(xiàn)狀有必要對非局部時滯反應(yīng)擴散方程行波解的存在性展開更深入的探討來豐富行波解的理論研究與增強應(yīng)用能力,具有重要的理論與實踐研究價值。
傳染病直接影響了人類的身體健康,傳染病在傳播過程中不僅存在時間滯后,也與攜帶者、易感染著等人口的空間位置是緊密相關(guān)的,具有典型的時間——空間關(guān)系,對此探討一類具有非局部擴散時滯反應(yīng)的SIR模型,模型如下所示:
該SIR模型還滿足以下條件:
(H1)當(dāng)S≥0時,f’(S)為有界正函數(shù);
對上述模型(2)行波解的存在性進行證明,即需尋找行波方程(3)的存在性,也就是尋找方程(3)滿足邊界條件的解,其中為大于0的常數(shù),代表了感染前易感人群的密度。
因此通過簡單的計算可得出以下引理。
引理1設(shè)定大于0,則有大于0的C*與λ*使得
引理2存在任意大于0的常數(shù)C,對任意大于的X 都有和
根據(jù)Schauder不動點定理可知(SX,IX)為算子的不動點,對于任意屬于(-X,X)范圍內(nèi)下的ξ,不動點(SX,IX)滿足下列關(guān)系式:
對式(6)進行(-X,X)區(qū)域內(nèi)的積分可得
對此可以得到下式(10):
同樣對式(7)進行(-X,X)積分,得到如(11)式
聯(lián)合式(10)至(12)式可得
根據(jù)前面的條件(H3)可以了解到式(13)是右側(cè)是有界的,故而存在一個大于0的常數(shù)C使得
因此存在兩個正常數(shù)L1與L2,對于屬于[-X,X]范圍內(nèi)的ξ和η有
根據(jù)式(6)可得到
又因核函數(shù)J 是Lipschitz 連續(xù),因此存在大于0 的Lipschitz 常數(shù) LJ,使得,其中 suppyJ=[-N,N],N 為大于0的常數(shù),故經(jīng)過計算與變換有
并根據(jù)式(16)可以得到如下式(19)
同理對式(7)作同樣處理最終可得到式(20)和(21),如下所示:
聯(lián)立式(7)、(14)和(21)可以得到
(二)行波解的存在性分析。本次證明中假設(shè)C>C*,且R0>1對此進行SIR模型的行波解的存在性分析,提出以下定理。
定理設(shè)定大于1,那么對于任意大于C*的C,SIR模型存在行波解滿 足并且有
證明過程如下所示:
將和帶入到式(3)第一個方程可得
在憑借文獻[12]的計算方法可知,當(dāng)y屬于suppyJ時,則有下式(23)
進一步簡化可得到式(24)
當(dāng)前對非局部時滯反應(yīng)擴散方程的行波解性質(zhì)研究成為了高等應(yīng)用數(shù)學(xué)的熱門研究話題。本文基于一類具有非局部擴散的時滯傳染病SIR模型,通過Schauder不動點定理、Fubini's定理以及極值法等方法對該模型行波解的存在性進行了分析與證明。通過研究發(fā)現(xiàn)本次研究模型存在滿足邊界條件的形態(tài)如的穩(wěn)態(tài)解,證明結(jié)果顯示能方法能有效解決該類SIR模型行波解的存在性問題。