一類隨機蚊蟲種群模型及其漸近性態
2021-03-23 07:17:42胡良劍
安徽工程大學學報 2021年1期
胡 悅，胡良劍
(東華大學 理學院，上海 201620)
近年來，傳染病數學模型的研究得到了廣泛關注，這些研究為傳染病的預測和控制提供了理論依據。2020年8月，新加坡遇到史上最嚴重的登革熱疫情，我國廣西和廣東等省份也出現了相關病例。登革熱、流行性乙型腦炎、瘧疾、絲蟲病和黃熱病等危害性很強的傳染病是由一類蚊蟲傳播的自然疫源性疾病，稱為蚊媒傳染病。噴灑農藥等方式并不能有效地阻止這類傳染病，過量使用還會導致環境污染、生態破壞等后果。于是人們研究出一種新型的控制策略來阻斷這類蚊蟲傳播疾病，即利用沃爾巴克氏菌入侵這類蚊蟲。這種細菌的繁殖方式屬于母系遺傳，并且通過驅動細胞質不親和性這一特有的細菌繁殖機制(Cytoplasmic Incompatibility,簡稱CI機制)入侵蚊蟲種群，就可以獲得一種較為理想的病毒傳播特性。這種特性會導致未受沃爾巴克氏菌感染的雌性蚊蟲和受感染的雄性蚊蟲交配產生的胚胎死亡，因此蚊蟲的數量會大大降低，從而達到減少蚊媒傳染病傳播的目的。在此期間，受沃爾巴克氏菌感染的雌性蚊蟲數量則不受任何影響，這樣便保證了沃爾巴克氏菌入侵蚊蟲種群的優勢，使得蚊媒傳染病能夠獲得有效的控制，以有利于人類健康和社會穩定。
一些學者對沃爾巴克氏菌影響下的蚊蟲種群模型進行了深入的研究。Fine對庫蚊細胞質不相容的種群動態進行研究分析，指出了雄性的攜帶者和雌性的非攜帶者之間有不親和性，僅僅通過母體遺傳特征，蚊蟲就能夠在一個種群中傳播寄生蟲和傳染病。Huang等考慮了一個反應擴散模型，該模型是在下列沃爾巴克氏菌影響下的確定性蚊蟲種群模型的基礎上改進的：
(1)
a
和b
替換為a
(t
)和b
(t
),其他參數c
、m
保持不變的條件下進行的。Hu等分別考慮了環境條件在兩個亞體系和多個亞體系之間的切換，把原本固定的參數變成隨時間等因素變化的量，最后估計了沃爾巴克氏固定的閾值釋放水平，并定義了疾病傳播的警報期，以衡量蚊蟲傳播疾病的風險。
(2)
式中，B
(t
)為標準的布朗運動。由于噪聲項B
(t
)的影響, 式(2)的狀態成為擴散過程,因此研究式(1)的常微分方程理論不再適用于式(2),需要運用隨機分析理論來進行研究。
1 全局正解的存在唯一性
為了保證方程的現實意義，首先需要驗證系統(2)全局正解的存在唯一性。顯然模型的系數是局部李普希茲連續的，但不滿足線性增長條件，所以普通的存在唯一性定理就不能適用于此系統，研究將通過構造李雅普諾夫函數來進行證明。
構造李雅普諾夫函數
V
(x
,y
)=x
-lnx
+y
-lny
,根據伊藤公式得
由
得
EV
(x
(T
∧τ
),y
(T
∧τ
))≤V
(x
(0),y
(0))+QT
,
因此可以得到
k
→∞時,V
(x
(0),y
(0))+QT
≥∞，這與V
(x
(0),y
(0))+QT
<∞矛盾，因此τ
=∞a
.s
。定理得證。
2 滅絕性
在研究種群的動力學行為時，滅絕性和持久性是兩個重要的系統特性。本節先討論系統的滅絕性。在探究滅絕性的過程中，基本再生數是系統的一個閾值，當基本再生數小于1時，種群便會隨著時間的推移而滅絕。
證明
利用伊藤公式，
t
積分
t
→∞時，x
(t
)指數趨向于0。證畢。
3 持久性
本節研究系統的持久性。 在證明持久性之前，先介紹引理1，證明見文獻[20]。
引理
1 設N
(t
)是下列方程具有初值N
>0的解，dN
(t
)=N
(t
)[α
-βN
(t
)]dt
+γN
(t
)dB
(t
)，若2α
>γ
，則
并且
注：該引理給出了隨機Logistic單種群模型的基本特性，其他復雜的隨機種群模型的研究可借鑒該系統的性質作進一步探討。
即蚊蟲的總數是隨機平均持久的。
證明
由于dx
(t
)=x
[a
-c
(x
+y
)]dt
+σxdB
(t
)≤x
(a
-cx
)dt
+σxdB
(t
),令x
(t
)是方程dx
=x
(a
-cx
)dt
+σx
dB
(t
)的解,根據引理1有
顯然
(3)
令y
(t
)是方程dy
=y
(b
-my
)dt
(4)
的解,那么
由式(3)、式(4)得
dx
(t
)=x
[a
-c
(x
+y
)]dt
+σxdB
(t
)≥x
(a
-cx
-cy
)dt
+σxdB
(t
)。令x
(t
)是方程dx
(t
)=x
[a
-cx
-cy
)dt
+σx
dB
(t
)的解,那么根據文獻[18]，
(5)
從而可設
根據伊藤公式得到
t
積分
整理后有
以及
兩式相加有
其中,
從而
結合式(5)可知
那么，
根據大數定律，
結合比較定理可知
因此有
由
以及大數定律得到
定理得證。
4 數值例子
為了證明以上的理論結果，利用Euler-Maruyama的方法來進行數值驗證。研究取初值x
(0)=y
(0)=0.5，即受感染蚊蟲和未受感染的蚊蟲數量各占一半。在系統(2)中設0≤t
≤100，時間步長取1，a
=0.1,b
=0.2,c
=0.3,m
=0.4。σ
=0.9時動力學行為如圖1所示。σ
=0.05時動力學行為如圖2所示。圖1、圖2中的實線代表受感染蚊蟲數量，虛線代表未受感染蚊蟲數量。圖1中的參數滿足定理2的條件R=(2a
/σ
)<1，數值結果說明,噪聲干擾足夠大時，受感染蚊蟲會隨著時間趨于滅絕，而未受感染的蚊蟲則不會滅絕；圖2參數不滿足定理2的條件，可見噪聲干擾較小時，受感染的蚊蟲不會滅絕。系統的持久性如圖3所示。圖3中曲線代表蚊蟲總和，這里a
=0.7，σ
=0.4。由圖3可知，基本再生數R足夠大時，隨著時間的變化，蚊蟲總數長期平均值會穩定在一個正值的趨勢。
圖1 σ=0.9時動力學行為 圖2 σ=0.05時動力學行為
圖3 系統的持久性
5 結論
進一步，可以探討模型(2)數值方法的收斂性和穩定性。還可以考慮將模型的其他參數，如未受感染的蚊蟲出生率的隨機擾動，研究其全局正解的存在唯一性、滅絕性以及持久性等。
猜你喜歡
童話王國·奇妙邏輯推理(2024年5期)2024-06-19 16:03:38
工業設計(2022年8期)2022-09-09 07:43:20
體育科技文獻通報(2022年3期)2022-05-23 13:46:54
遼金歷史與考古(2021年0期)2021-07-29 01:06:54
軍民兩用技術與產品(2021年10期)2021-03-16 06:05:30
北京測繪(2020年12期)2020-12-29 01:33:58
中學生數理化·七年級數學人教版(2020年10期)2020-11-26 08:24:50
數學物理學報(2020年2期)2020-06-02 11:29:24
科技傳播(2019年22期)2020-01-14 03:06:54
民用飛機設計與研究(2019年4期)2019-05-21 07:21:24
