一類隨機戒煙模型的滅絕性與均值持久性

戴祥軍，徐松金，劉 聲，李 通

?

一類隨機戒煙模型的滅絕性與均值持久性

戴祥軍，徐松金，劉 聲，李 通

（銅仁學院 大數據學院，貴州 銅仁 554300 ）

主要探討了一類隨機戒煙模型的動力學性質。首先討論了該系統存在唯一的全局正解，接下來獲得了該系統各類人群的滅絕性與均值持久性的一個充分條件，最后應用數值模擬對理論結果進行了驗證。

隨機干擾； 戒煙模型； 滅絕性； 持久性

0． 引言

眾所周知，吸煙對人類身體健康有著潛在的危害。2012年，世界衛生組織在報告上說，全世界每年因吸煙造成的各類疾病而死亡的人數高達500萬，也就是說香煙每小時平均殺死570人，超過了艾滋病引起的死亡人數。如果繼續按照這種情況下去，世界健康組織預測到2030年因吸煙造成的各類疾病而死亡的人數將高達1000萬[1]。最近幾年中國每年有100多萬人因吸煙造成各類疾病而死亡。雖然如此，仍然有大量的公民對此熟視無睹繼續保持著鋌而走險的抽煙習慣。因此我們每個人都應該了解吸煙的危害進而變成自覺戒煙的擁護者，生物數學工作者也貢獻自己的力量[2-10]。

在1997年，Castillo-Garsow etal[2]提出了一類簡單的戒煙模型。他們根據吸煙的數量多少把他們分成三類：P表示潛在吸煙者，S表示吸煙者，Q表示戒煙者。表示如下：

1． 準備工作

本節給出一些將要用到的記號和引理[11-14]。

為了方便起見，我們給出下列記號

2． 全局正解的存在唯一性

模型中所有的變量都是對人口密度的描述，那么它應該都是非負的。因此為了進一步對戒煙系統的研究，我們首先來討論模型(3)中的解的存在性與唯一性。

定義一個函數

應用伊藤公式計算可得

其中

那么

3． 滅絕性與均值持久性

本節我們主要是探討系統的滅絕性（即吸煙人群趨向滅絕，吸煙將會得到控制）和持久性（即吸煙人群持久存在，吸煙將會盛行）。

證明：對方程（3）我們應用伊藤公式可得，

對方程（3）兩邊同時積分后可得

通過計算可知

其中

把(7)式代入不等式(6)可得

其中

如果條件（a）成立，由不等式(8)，我們兩邊同時取上極限可得

同樣，由不等式(6)可知

如果條件（b）成立，則有

解方程(3)的第三個方程可得

解之得

另一方面，把(7)式代入到(6)式得

對上式兩邊同時取下極限可得

4． 數值模擬

圖2 當時隨機系統（1.2）的解

圖3 當時隨機系統（1.2）的解

5． 結論

[1] Zeb A, Zaman G, Momani S. Square dynamics of a giving up smoking model[J]. Appl. math. model, 2013, 37: 5326-5334.

[2] Castillo-Garsow C, Jordan-Salivia G, Herrera A R. Mathematical models for the dynamics of tobacco use, recovery, and relapse. Technical Report Series BU-1505-M, Cornell University, Ithaca, New York, 1997.

[3] Zaman G. Qualitative behavior of giving up smoking models[J]. Bulletin of the Malaysian Mathematical Sciences Society, 2011, 34（2）：403-415.

[4] Sharomi O, Gumel A B. Curtailing smoking dynamics: a mathematical modeling approach[J]. Applied Mathematics and Computation, 2008, 195（2）：475-499.

[5] Erturk V S, Zaman G, Momanl S. A numeric-analytic method for approximating a giving up smoking model containing fractional derivatives[J]. Computers Mathematics with Applications, 2008, 64（10）：3065- 3074.

[6] 霍海峰，羅琪．一類具有階段結構的戒煙模型的全局漸近穩定性[J]．蘭州理工大學學報，2015，41（5）：148-151．

[7] Qamar D, Muhammad O, Takasar H. Qualitative behavior of smoking models[J]. Advances in Difference Equations, 2016（1）：1-12.

[8] McCluskey C C, vanden Driessche P. Global analysis of two tuberculosis models[J]. Journal of Dynamics and Differential Equations, 2004, 16（1）：139-166.

[9] CAI L, LI X, GHOSH M. Stability analysis of an HIV/AIDS epidemic model with treatment[J]. Journal of Computational and Applied Mathematics, 2009, 229（1）：313-323.

[10] Lahrouz A, OMARI L, Kiouach D. Deterministic and stochastic stability of a mathematical model of smoking[J]. Statistics Probability Letters, 2011, 81（8）：1276-1284.

[11] 王克．隨機生物數學[M]．北京：科學出版社，2010．

[12] Zhao Y N, Jiang D Q , O’REGAN D. The extinction and persistence of the stochastic SIS epidemic model with vaccination[J]. Physica A, 2013, 392（1）：4916-4927.

[13] Ji C Y, Jiang D Q, SHI N Z. Two-group SIR Epidemic Model with Stochastic Perturbation[J]. Acta Mathematica Sinica, English Series, 2012, 28（12）：2545-2560.

[14] Higham D J. An algorithmic introduction to numerical simulation of stochastic differential equations[J]. SIAM Rev, 2001, 43：525-546．

Extinction and Persistence in the Mean of a Stochastic Giving up Smoking Model

DAI Xiangjun, XUSongjin, Liu Sheng, LI Tong

( School of Data Science, Tongren University, Tongren 554300, Guizhou, China)

In this paper, the dynamics of a stochastic giving up smoking model are considered. First, we show that there is a unique global positive solution of this system. And then, sufficient conditions for the persistence and extinction of this system are established. The last, the results are illustrated by computer simulations.

stochastic perturbation, giving up smoking model, persistence, extinction

2018-09-04

貴州省創新群體重大研究項目（黔教合KY字[2016]051）；貴州省一流大學（一期）課程培育項目（2017YLDXXM005）；貴州省普通高等學校創新人才團隊項目（黔教合人才團隊字[2015]64）；大學生創新訓練項目（201710665018）。

戴祥軍（1989-），男，湖南邵陽人，講師，碩士，研究方向：隨機生物數學。

O175

A

1673-9639 (2018) 12-0111-06

（責任編輯 毛 志）（責任校對 印有家）