Dynamics of an infectious diseases with media/psychology induced non-smooth incidence

  • Received: 01 May 2012 Accepted: 29 June 2018 Published: 01 January 2013
  • MSC : 92C50, 92C60, 54B20.

  • This paper proposes and analyzes a mathematical model on aninfectious disease system with a piecewise smooth incidence rateconcerning media/psychological effect. The proposed models extendthe classic models with media coverage by including a piecewisesmooth incidence rate to represent that the reduction factor becauseof media coverage depends on both the number of cases and the rateof changes in case number. On the basis of properties of Lambert Wfunction the implicitly defined model has been converted into apiecewise smooth system with explicit definition, and the globaldynamic behavior is theoretically examined. The disease-free isglobally asymptotically stable when a certain threshold is less thanunity, while the endemic equilibrium is globally asymptoticallystable for otherwise. The media/psychological impact although doesnot affect the epidemic threshold, delays the epidemic peak andresults in a lower size of outbreak (or equilibrium level ofinfected individuals).

    Citation: Yanni Xiao, Tingting Zhao, Sanyi Tang. Dynamics of an infectious diseases with media/psychology induced non-smooth incidence[J]. Mathematical Biosciences and Engineering, 2013, 10(2): 445-461. doi: 10.3934/mbe.2013.10.445

    Related Papers:

    [1] Jawhar Hbil, Mohamed Zaway . Some results on the convergence of Hessian operator and $ m-$subharmonic functions. AIMS Mathematics, 2022, 7(5): 9023-9038. doi: 10.3934/math.2022502
    [2] M. Emin Özdemir, Saad I. Butt, Bahtiyar Bayraktar, Jamshed Nasir . Several integral inequalities for (α, s,m)-convex functions. AIMS Mathematics, 2020, 5(4): 3906-3921. doi: 10.3934/math.2020253
    [3] Ghulam Farid, Saira Bano Akbar, Shafiq Ur Rehman, Josip Pečarić . Boundedness of fractional integral operators containing Mittag-Leffler functions via (s,m)-convexity. AIMS Mathematics, 2020, 5(2): 966-978. doi: 10.3934/math.2020067
    [4] Siamak Khalehoghli, Hamidreza Rahimi, Madjid Eshaghi Gordji . Fixed point theorems in R-metric spaces with applications. AIMS Mathematics, 2020, 5(4): 3125-3137. doi: 10.3934/math.2020201
    [5] Ehsan Movahednia, Young Cho, Choonkil Park, Siriluk Paokanta . On approximate solution of lattice functional equations in Banach f-algebras. AIMS Mathematics, 2020, 5(6): 5458-5469. doi: 10.3934/math.2020350
    [6] Khadija Gherairi, Zayd Hajjej, Haiyan Li, Hedi Regeiba . Some properties of $ n $-quasi-$ (m, q) $-isometric operators on a Banach space. AIMS Mathematics, 2023, 8(12): 31246-31257. doi: 10.3934/math.20231599
    [7] Khadija Gherairi, Zayd Hajjej, Haiyan Li, Hedi Regeiba . $ n $-quasi-$ A $-$ (m, q) $-isometry on a Banach space. AIMS Mathematics, 2023, 8(12): 28308-28321. doi: 10.3934/math.20231448
    [8] Murali Ramdoss, Divyakumari Pachaiyappan, Inho Hwang, Choonkil Park . Stability of an n-variable mixed type functional equation in probabilistic modular spaces. AIMS Mathematics, 2020, 5(6): 5903-5915. doi: 10.3934/math.2020378
    [9] Haikun Liu, Yongqiang Fu . On the variable exponential fractional Sobolev space Ws(·),p(·). AIMS Mathematics, 2020, 5(6): 6261-6276. doi: 10.3934/math.2020403
    [10] Jae-Hyeong Bae, Won-Gil Park . Approximate solution of a multivariable Cauchy-Jensen functional equation. AIMS Mathematics, 2024, 9(3): 7084-7094. doi: 10.3934/math.2024345
  • This paper proposes and analyzes a mathematical model on aninfectious disease system with a piecewise smooth incidence rateconcerning media/psychological effect. The proposed models extendthe classic models with media coverage by including a piecewisesmooth incidence rate to represent that the reduction factor becauseof media coverage depends on both the number of cases and the rateof changes in case number. On the basis of properties of Lambert Wfunction the implicitly defined model has been converted into apiecewise smooth system with explicit definition, and the globaldynamic behavior is theoretically examined. The disease-free isglobally asymptotically stable when a certain threshold is less thanunity, while the endemic equilibrium is globally asymptoticallystable for otherwise. The media/psychological impact although doesnot affect the epidemic threshold, delays the epidemic peak andresults in a lower size of outbreak (or equilibrium level ofinfected individuals).


    Theory of $ m $-subharmonic functions was recently developed by many mathematicians such as Li [20], Błocki [9], Dinew and Kołodziej [14,15], Lu [21,22], Sadullaev and Abdullaev [30], Nguyen [23,24], Åhag, Czyż and Hed [3,4] and many others. The notion of $ m $-subharmonicity appears naturally in generalization of subharmonicity and plurisubharmonicity. For the similarities and the differences between these notions, we refer the readers to the paper [15].

    A bounded domain $ \Omega\subset \mathbb{C}^n $ is called $ m $-hyperconvex if there exists an $ m $-subharmonic function $ \rho\colon \Omega \to (-\infty, 0) $ such that the closure of the set $ \{z\in \Omega \colon \rho(z) < c\} $ is compact in $ \Omega $ for every $ c\in (-\infty, 0) $. In what follows we will always assume that $ \Omega $ is an $ m $-hyperconvex domain. Denote by $ SH_m(\Omega) $ the set of all $ m $-subharmonic functions in $ \Omega $. Let the cones $ \mathcal{E}_{0, m}, \mathcal{E}_{p, m}, \mathcal{F}_m $ be defined in the similar way as in [21,25]:

    $ E0,m={uSHm(Ω)L(Ω):limzΩu(z)=0 and ΩHm(u)<},Ep,m={uSHm(Ω): {uj}E0,m,uju,supjΩ(uj)pHm(uj)<},Fm={uSHm(Ω): {uj}E0,m,uju and supjΩHm(uj)<}.
    $

    For the properties and applications of these classes, see [1,21,22,25,26,27].

    We use the notation $ \delta \mathcal{K} = \mathcal{K} = \mathcal{K} $ for $ \mathcal{K} $ be one of the classes $ \mathcal{E}_{0, m}, \mathcal{E}_{p, m}, \mathcal{F}_m $. Define

    $ ||u||p,m=infu=u1u2u1,u2Ep,m{(Ω(u1u2)pHm(u1+u2))1m+p},
    $
    (1.1)

    with the convention that $ (-u_1-u_2)^p = 1 $ if $ p = 0 $. For the reason why this quasi-norm is effective, please see [2,13,16,22,29]. It was proved in [25] that $ (\delta \mathcal{E}_{p, m}, ||\cdot||_{p, m}) $ is a quasi-Banach space for $ p > 0, p\neq 1 $ and it is a Banach space if $ p = 1 $. Moreover in [17] it was proved that $ (\delta \mathcal{F}_m, ||\cdot||_{0, m}) $ is a Banach space. The authors in [12] show that $ (\delta \mathcal{E}_{p, m}, ||\cdot||_{p, m}) $ can not be a Banach space. These facts are counterparts of [5,6,10,18] in $ m $-subharmonic setting.

    In Section 2, we shall show that $ \mathcal{E}_{0, m} $ and $ \delta \mathcal{E}_{0, m} $ are closed neither in $ (\delta \mathcal{E}_{p, m}, ||\cdot||_{p, m}) $ nor in $ (\delta \mathcal{F}_m, ||\cdot||_{0, m}) $. Moreover we prove that the inclusions $ \mathcal{E}_{0, m}\subseteq \mathcal{F}_m, \delta \mathcal{E}_{0, m}\subseteq \delta \mathcal{F}_m $ are proper in the space $ (\delta \mathcal{F}_m, ||\cdot||_{0, m}) $.

    In Section 3, we prove that the convergence in $ \delta \mathcal{E}_{p, m} $ implies the convergence in $ m $-capacity (Theorem 3). But the convergence in $ m $-capacity is not a sufficient condition for the convergence in $ \delta \mathcal{E}_{p, m} $ (Example 3). Similar results in plurisubharmonic setting have been proved by Czyż in [11].

    In plurisubharmonic case, the following proposition was proved in (see [11]). Let $ \mathbb{B} = \mathbb{B} (0, 1)\subset \mathbb{C}^n $ be the unit ball in $ \mathbb{C}^n $. Then the cones $ \mathcal{E}_{0, m}(\mathbb{B}) $ and $ \delta \mathcal{E}_{0, m}(\mathbb{B}) $ are not closed respectively in $ (\delta \mathcal{F}_m(\mathbb{B}), ||\cdot||_{0, m}) $ and $ (\delta \mathcal{E}_{p, m}(\mathbb{B}), ||\cdot||_{p, m}) $.

    Proof. We define

    $ v(z) = {ln|z|if m=n,1|z|22nmif 1m<n.
    $

    We obtain that $ H_m(v): = dd^c(v) \wedge\beta^{n-m} = c(n, m)\delta_0 $, where $ c(n, m) $ is a constant depending only on $ n $ and $ m $, $ \delta_0 $ is the Dirac measure at the origin $ 0 $ (see [28]). For each $ j \in \mathbb{N} $, define the function $ v_j \colon \mathbb{B} \to \mathbb{R}\cup \{-\infty\} $ by

    $ v_j(z) = \max(a_jv(z),-b_j), $

    where $ a_j = \frac{1}{2^j}, b_j = \frac{1}{j} $.

    We can see that $ v_j\in \mathcal{E}_{0, m}(\mathbb{B}) $, for each $ j $. Therefore, the function $ u_k: = \sum_{j = 1}^k v_j $ belongs to $ \mathcal{E}_{0, m}(\mathbb{B}) $. For $ k > l $ we can compute

    $ ||ukul||m0,m=||kj=l+1vj||m=BHm(kj=l+1vj)=c(n,m)(kj=l+1aj)m,
    $
    (2.1)

    and

    $ ||ukul||p+mp,m=||kj=l+1vj||p+mp+m=ep,m(kj=l+1vj)=B(kj=l+1vj)pHm(kj=l+1vj)=c(n,m)kj1,,jm=l+1[kr=l+1vr(max(tj1,,tjm))]paj1ajmc(n,m)kj1,,jm=l+1[uk(max(tj1,,tjm))]paj1ajmc(n,m)[kj=l+1(uk(tj))pmaj]m,
    $

    where

    $ t_j = {(1+bjaj)m2(mn),if 1m<n,ebjaj,if m=n.
    $

    The last inequality is a consequence of the fact that $ v_j $ is increasing function for each $ j $. Since

    $ v_l(t_j) = {1l,if 1lj,2jj2l,if l>j,
    $

    we have

    $ -u_k(t_j) = \sum\limits_{l = 1}^j \frac{1}{l}+\frac{2^j}{j}\sum\limits_{l = j+1}^k\frac{1}{2^l}\leq j+1. $

    Hence

    $ ||ukul||p+mp,mc(n,m)(kj=l+1(j+1)pm2j)m.
    $
    (2.2)

    Let $ u\colon \mathbb{B}\to \mathbb{R}\cup\{-\infty\} $ be defined by $ u = \lim_{k\to\infty}u_k $. Observe that $ u $ is the limit of a decreasing sequence of $ m $-subharmonic functions and $ u(z) > -\infty $ on the boundary of the ball $ B(0, \frac{1}{2}) $. Hence $ u $ is $ m $-subharmonic. Moreover $ u\notin \mathcal{E}_{0, m}(\mathbb{B}) $ since it is not bounded on $ \mathbb{B} $, its value is not bounded below at the origin. Equality (2.1) shows that $ \{u_k\} $ is a Cauchy sequence in the space $ \delta \mathcal{F}_m(\mathbb{B}) $. Thus the cone $ \mathcal{E}_{0, m}(\mathbb{B}) $ and the space $ \delta \mathcal{E}_{0, m}(\mathbb{B}) $ are not closed in $ (\delta \mathcal{F}_m(\mathbb{B}), ||\cdot||_{0, m}) $.

    The series $ \sum_{j = 1}^\infty \frac{(j+1)^{\frac{p}{m}}}{2^j} $ is convergent by the ratio test. Therefore $ \{u_k\} $ is a Cauchy sequence in $ \delta \mathcal{E}_{p, m} $ by (2.2). We have proved that the cone $ \mathcal{E}_{0, m}(\mathbb{B}) $ and the space $ \delta \mathcal{E}_{0, m}(\mathbb{B}) $ are not closed in $ (\delta \mathcal{E}_{p, m}(\mathbb{B}), ||\cdot||_{p, m}) $.

    The following proposition shows that the closure of the cone $ \mathcal{E}_{0, m} $ (resp. $ \delta \mathcal{E}_{0, m} $) is strictly smaller than $ \mathcal{F}_m $ (resp. $ \delta \mathcal{F}_m $) in the space $ (\delta \mathcal{F}_m, ||\cdot||_{0, m}) $. We have $ \overline{ \mathcal{E}_{0, m}}\subsetneq \mathcal{F}_m $ and $ \overline{\delta \mathcal{E}_{0, m}}\subsetneq \delta \mathcal{F}_m $ in the space $ (\delta \mathcal{F}_m, ||\cdot||_{0, m}) $.

    Proof. The definition of the $ m $-Lelong number of a function $ v\in SH_m(\Omega) $ at $ a\in \Omega $ is the following

    $ \nu_{m,a}(v) = \lim\limits_{r\to 0^+}\int_{|z-a|\leq r} dd^cv \wedge \left[dd^c(-|z-a|^{2-\frac{2n}{m}})\right]^{m-1} \wedge\beta^{n-m} $

    It is easy to see that $ m $-Lelong number is a linear functional on $ \delta \mathcal{F}_m $. Moreover, as in [7, Remark 1], for a function $ \varphi\in \mathcal{F}_m $ then

    $ \nu_{m,a}(\varphi)\leq (H_m(\varphi)(\{a\}))^{\frac{1}{m}}\leq (H_m(\varphi)(\Omega))^\frac{1}{m}. $

    Hence, for any representation $ u = u_1-u_2 $ of $ u\in \delta \mathcal{F}_m $ we have

    $ |\nu_{m,a}(u)|\leq (H_m(u_1+u_2)(\Omega))^\frac{1}{m}. $

    This implies that $ m $-Lelong number is a bounded functional on the space $ \delta \mathcal{F}_m $. We have shown that $ m $-Lelong number is continuous on the Banach space $ (\delta \mathcal{F}_m, ||\cdot||_{0, m}) $. We recall the definition of $ m $-Green function with pole at $ a $

    $ g_{m,\Omega,a}(z) = \sup\{v\in SH_m^-(\Omega)\colon u(z)+|z-a|^{2-\frac{2n}{m}}\leq O(1) \text{ as } z\to a\}. $

    The readers can find more properties of $ m $-Green function in [31]. Assume that $ \overline{ \mathcal{E}_{0, m}} = \mathcal{F}_m $. Then there exists a sequence $ \{u_j\} $ in $ \mathcal{E}_{0, m} $ that converges to $ g_{m, \Omega, a} $ in the space $ \delta \mathcal{F}_m $ as $ j\to \infty $. The $ m $-Lelong number of all $ u_j $ at $ a $ vanishes since $ u_j $ is bounded, but the $ m $-Lelong number of $ g_{m, \Omega, a} $ at $ a $ is 1. Hence we get a contradiction. Thus, $ \overline{ \mathcal{E}_{0, m}}\subsetneq \mathcal{F}_m $. By the same argument, if $ \overline{\delta \mathcal{E}_{0, m}} = \delta \mathcal{F}_m $, then there exists a sequence $ \{u_j\} $ in $ \mathcal{E}_{0, m} $ that converges to $ g_{m, \Omega, a} $ in the space $ \delta \mathcal{F}_m $ as $ j\to \infty $, but this is impossible since $ \nu_{m, a}(u_j) = 0 $.

    We are going to recall a Błocki type inequality (see [8]) for the class $ \mathcal{E}_{p, m} $. Similar results for the class $ \mathcal{F}_m $ were proved by Hung and Phu in [19, Proposition 5.3] (see also [1]) and for locally bounded functions were proved by Wan and Wang [31]. Assume that $ v\in \mathcal{E}_{p, m} $ and $ h\in SH_m $ is such that $ -1\leq h\leq 0 $. Then

    $ \int_\Omega (-v)^{m+p}H_m(h)\leq m! \int_\Omega (-v)^pH_m(v). $

    Proof. See the proof of [19, Proposition 5.3].

    Recall that the relative $ m $-capacity of a Borel set $ E\subset\Omega $ with respect to $ \Omega $ is defined by

    $ cap_{m,\Omega}(E) = \sup\{\int_E H_m(u)\colon u\in SH_m(\Omega), -1\leq u\leq 0\}. $

    We are going to recall the convergence in $ m $-capacity. We say that a sequence $ \{u_j\}\subset SH_m(\Omega) $ converges to $ u\in SH_m(\Omega) $ in $ m $-capacity if for any $ \epsilon > 0 $ and $ K\Subset\Omega $ then we have

    $ \lim\limits_{j\to\infty} cap_{m,\Omega}(K\cap \{|u_j-u| \gt \epsilon\}) = 0. $

    Let $ \{u_j\}\subset \delta \mathcal{E}_{p, m} $ be a sequence that converges to a function $ u\in \delta \mathcal{E}_{p, m} $ as $ j $ tends to $ \infty $. Then $ \{u_j\} $ converges to $ u $ in $ m $-capacity.

    Proof. Replacing $ u_j $ by $ u_j-u $, we can assume that $ u = 0 $. By the definition of $ \delta \mathcal{E}_{p, m} $, there exist functions $ v_j, w_j\in \mathcal{E}_{p, m} $ such that $ u_j = v_j-w_j $ and $ e_p(v_j+w_j)\to 0 $ as $ j\to \infty $. By [25],

    $ \max(e_{p,m}(v_i),e_{p,m}(w_j))\leq e_{p,m}(v_j+w_j), $

    which implies that $ e_{p, m}(v_j), e_{p, m}(w_j) $ tend to $ 0 $ as $ j\to \infty $. Given $ \epsilon > 0 $ and $ K\Subset \Omega $. For a function $ \varphi\in SH_m(\Omega) $, $ -1\leq \varphi\leq 0 $, we have

    $ {|vj|>ϵ}KHm(φ)1ϵp+mΩ(vj)p+mHm(φ)m!ϵp+mep,m(vj).
    $
    (3.1)

    The last inequality comes from Lemma 3. Hence, by taking the supremum over all functions $ \varphi $ in inequality (3.1), we get

    $ capm,Ω({|vj|>ϵ}K)m!ϵm+pep,m(vj).
    $
    (3.2)

    Similarly,

    $ capm,Ω({|wj|>ϵ}K)m!ϵm+pep,m(wj).
    $
    (3.3)

    From (3.2), (3.3) we obtain

    $ capm,Ω({|uj|>ϵ}K)capm,Ω({|vj|>ϵ2}K)+capm,Ω({|wj|>ϵ2}K)m!2m+pϵm+p(ep,m(vj)+ep,m(wj))0 as j.
    $

    Hence the sequence $ \{u_j\} $ tends to $ 0 $ in $ m $-capacity and the proof is finished.

    A similar result for the space $ \delta \mathcal{F}_m $ is proved in [17]. But the convergence in $ m $-capacity is not a sufficient condition for the convergence in the space $ \delta \mathcal{E}_{p, m} $. The following example shows that convergence in $ m $-capacity is strictly weaker than convergences in both $ \delta \mathcal{E}_{p, m} $ and $ \delta \mathcal{F}_m $. The case $ m = n $ has been showed in [11, Example 3.3]. Let $ v(z) $ be the function defined in the unit ball in $ \mathbb{C}^n $ as in the proof of Proposition 2. We define

    $ u_j(z) = \max(j^{\frac{p}{m}}v(z), -\frac{1}{j}), \ v_j(z) = \max(v(z),-\frac{1}{j}) $

    Then we have $ u_j, v_j\in \mathcal{E}_{0, m}(\mathbb{B}) $ for every $ j $, and $ e_{p, m}(u_j) = c(n, m), e_{0, m}(v_j) = 1 $. These show that the sequence $ \{u_j\} $ and $ \{v_j\} $ do not converge to $ 0 $ in $ \delta \mathcal{E}_{p, m}(\mathbb{B}) $ and $ \delta \mathcal{F}_m(\mathbb{B}) $ respectively as $ j\to \infty $. Moreover, for fixed $ \epsilon > 0 $ and $ K\Subset \mathbb{B} $ there exists $ j_0 $ such that for all $ j\geq j_0 $ we have

    $ u_j = v_j = -\frac{1}{j} \ \text{on } K. $

    This infers that both sets $ K\cap \{u_j < -\epsilon\} $ and $ K\cap \{v_j < -\epsilon\} $ are empty. Hence $ u_j $ and $ v_j $ tend to $ 0 $ in $ m $-capacity.

    The authors would like to thank Rafał Czyż for many valuable comments and suggestions for this manuscript. We are grateful to the referee whose remarks and comments helped to improve the paper.

    The authors declare no conflict of interest.

    [1] Oxford University, Oxford, 1991.
    [2] Springer-Verlag, Berlin, 1984.
    [3] Circ. Syst. Signal. Pr., 12 (1993), 375-390.
    [4] IEEE Press, New York, 2001.
    [5] SIAM Review, 50 (2008), 629-701.
    [6] Int. J. Bifurcat. Chaos, 11 (2001), 1121-1140.
    [7] Cochrane Db. Syst. Rev., 11 (2010) , 1-47.
    [8] Springer-Verlag, London, 1999.
    [9] Adv. Comput. Math., 5 (1996), 329-359.
    [10] J. Dynam. Diff. Eqns., 20 (2008), 31-53.
    [11] Rocky Mt. J. Math., 38 (2008), 1323-1334.
    [12] Walter De Gruyter, Berlin, 1992.
    [13] Kluwer Academic, Dordrecht, The Netherlands, 1988.
    [14] J. R. Soc. Interface, 7 (2010), 1247-1256
    [15] PNAS, 106 (2009), 6872-6877.
    [16] Proc. R. Soc. B., 274 (2007), 833-837.
    [17] European J. Control., 9 (2003), 311-319.
    [18] PLoS ONE, 4 (2009), e8032.
    [19] Comput. Math. Methods Med., 8 (2007), 153-164.
    [20] Int. J. Biomath., 1 (2008), 65-74.
    [21] J. Math. Anal. Appl., 303 (2004), 81-89.
    [22] Health Technol. Assess., 14 (2010), 183-266.
    [23] Math. Biosci., 230 (2011), 87-95.
    [24] BMC Public Health, 11 (2011), S5.
    [25] ISRN Biomathematics, (2012).
    [26] Math. Biosci., 201 (2006), 58-71.
    [27] IMA. J. Math. appl. Med. Biol., 28 (2010), 227-244.
    [28] Nonl. Anal. RWA., 11 (2010), 4154-4163.
    [29] Bull. Math. Biol., 74 (2012), 2403-2422.
  • This article has been cited by:

    1. Shengqiu Liu, Wei Wang, On Pluripotential Theory Associated to Quaternionic m-Subharmonic Functions, 2023, 33, 1050-6926, 10.1007/s12220-023-01197-x
  • Reader Comments
  • © 2013 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(4998) PDF downloads(825) Cited by(114)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog