# Global Dynamics of SIRS Model with No Full Immunity on Semidirected Networks

Authors: Jiawei Huo; Yimin Li; Jing Hua
Journal: Mathematical Problems in Engineering
 (2019)

---

## Abstract

In this paper, an epidemic model with no full immunity is analyzed on semidirected networks. Directed networks led into previous scale-free networks, and we consider that some infectious diseases do not have full immunity. So we use strong self-protection instead of immunity and establish a semidirected network infectious disease model without full immunity. The basic reproduction numberR0 is calculated. If R0<1, the disease-free equilibrium E0 is locally and globally asymptotically stable. And the endemic equilibrium E∗ is globally asymptotically stable in some condition. A large number of simulation results in this paper verify the correctness of the above conclusions and provide a solution for controlling disease transmission in the future.


---

## Body

## 1. Introduction

Based on mathematical models, especially the research on infectious diseases based on dynamics models has a history over 100 years. In 1873–1894, P D EN’KO established a model of modern mathematical infectious diseases [1]. In 1927, Kermack and McKendrick studied the SIR compartmental model after studying the Black Death and plague [2] and proposed the SIS compartmental model in 1932. Based on the compartmental model, threshold theory to distinguish whether epidemics exist was proposed [3]. After that, especially in the past 30 years, biomathematicians established and studied various infectious disease dynamic models for infectious disease compartmental models [4–6].In recent years, with the rapid development of complex networks, many practical problems can be abstracted into complex network models for research. The infectious disease model is also applied to the scale-free network of complex networks. In 2001, Pastor-Satorras and Vespignani used the average field theory to study the SIS infectious disease model on the general network, applied it to the scale-free network, and proved that the scale-free network does not have a threshold under the appropriate parameters [7]. In the same year, May and Lloyd gave the basic regeneration number of the general network in the scale-free SIR model [8]. In 2002, Pastor-Satorras and Vespignani studied the transmission threshold of finite-scale network-free infectious diseases and proposed consistent immunity and optimized immunity for SIS infectious disease models [9]. In 2004, Liu et al. considered the dynamic network model of birth-death infectious disease with the static network [10]. In the same year, Hayashi et al. proposed the SIR virus propagation model under the linear growth scale-free network [11]. With the deep understanding of the network, this symmetry hypothesis is often not correct in the process of propagation research, such as mother-to-child transmission, virus transmission on computer networks, and information dissemination are dissymmetric. A large number of networks in nature involve directed networks [12–14]. Li et al. established a directed network propagation model [15]. In a real network, the contact between nodes is not all symmetric and asymmetric, but a situation of directed and undirected coexistence. Sharkey et al. established pair-level approximations to the spatiotemporal dynamics of epidemics on asymmetric contact networks [16]. Meyers et al. distinguished from previous work by using the probabilistic parent function method to study semidirected network propagation problems [17]. Zhang et al. studied the SIS model of semidirected networks and gave detailed dynamic analysis [18]. In many communication processes, not all communication processes gain immunity, such as the spread of mobile viruses [19], as do many diseases. However, although the node has not obtained immunity, it is less likely to be transmitted than before. To research this propagation, this paper establishes a new propagation model on a semidirected network. This article replaces immunity with strong self-protection, and the expected conclusion is closer to the actual situation. And Liu et al. [20] and Huang [21] has researched this propagation in the scale-free network but not on semidirected networks. This paper set up this propagation model on the semidirected network. This is a new propagation model. This study wants to use the model to find the methods and final range of the propagation.Based on the aforementioned research results, this paper uses the idea of semidirected network and self-protection to establish a semidirected network epidemic model without full immunity, which can be used for the transmission of some certain diseases or computer viruses or mobile viruses. It may help to prevent the propagation of diseases in our daily living networks. And this study wants to research the methods of the propagation and the range affected by diseases. This paper systematically analyzes the dynamic properties of the model. The basic regeneration number and equilibrium point expression are calculated, and the stability of the equilibrium point is studied. The correctness of the result is proved by a large number of simulations, which can provide suggestions for the control of propagation.


## 2. Model Description

In this section, an epidemic model with no full immunity on semidirected complex networks is described. In the semidirected network, each node sends a directed or undirected connection to other nodes. The connection status of the node can be expressed byi,j;n, where i represents the in-degree, j represents the out-degree, and n represents the undirected degree. In the semidirected network, it is assumed that the disease propagates only along the out-degree edge and the undirected edge. The undirected degree and in-degree of each node indicate the possibility that the node will be propagated. The undirected degree and out-degree of each node indicate the possibility that the node propagates the disease to others. We use Min, Mout, and Mu to represent the max in-degree, the max out-degree, and the max undirected degree, respectively. Ni,j;n represents the number of nodes with in-degree i, out-degree j, and undirected degree n. Pi,j;n represents the possibility that the node with in-degree i, out-degree j, and undirected degree n is selected randomly.The states of node in the propagation model are divided intoS (susceptible), I (infected), and R (strong self-protection) states. Not all diseases will equip immune after healing, so this article introduces a strong self-protection state instead of immune status. The nodes which in state S or R may be infected by disease. And the probability of infection is represented by αS and αR. In order to reflect the self-protection of state R, we guarantee αS>αR. Using τd and τu, respectively, to indicate the probability the infected node will transmit disease to adjacent noninfected nodes through the directed and undirected edges. The susceptible nodes are transformed into strong self-protection nodes by the probability γ due to the influence of the surrounding strong self-protection nodes. The strong self-protection nodes are also transformed into susceptible nodes by probability η because it is not infected by disease for a long time. The recovery rate of the infected nodes is β. The number of susceptible nodes, infected nodes, and strong self-protection nodes with degree i,j;n is recorded as Si,j;n, Ii,j;n, and Ri,j;n, respectively. So Si,j;n+Ii,j;n+Ri,j;n=Ni,j;n. Ni,j;n is constant. We can express it by the relative density method. Si,j;nt+Ii,j;nt+Ri,j;nt=1. In order to reduce the length of the equation, we use Si,j;n, Ii,j;n, and Ri,j;n to replace Si,j;nt, Ii,j;nt, and Ri,j;nt. Based on this, we can get the following dynamic model of the semidirected network:(1)dSi,j;ndt=−αSτdiSi,j;nΘd+τunSi,j;nΘu+ηRi,j;n−γSi,j;n,dIi,j;ndt=αSτdiSi,j;nΘd+τunSi,j;nΘu+αRτdiRi,j;nΘd+τunRi,j;nΘu−βIi,j;n,dRi,j;ndt=−αRτdiRi,j;nΘd+τunRi,j;nΘu−ηRi,j;n+γSi,j;n+βIi,j;n,where 0≤i≤Min, 0≤j≤Mout, 0≤n≤Mu, Θd indicates the probability that the node in the network is connected to infected nodes through in-degree, and Θu indicates the probability that the node in the network is connected to infected nodes through undirected degree. Their mathematical expressions are as follows:(2)Θd=∑jPi,j;nIi,j;ntdout,Θu=∑nPi,j;nIi,j;ntdu,where dout=∑jPi,j;n, du=∑nPi,j;n, and din=∑iPi,j;n. Because this network is a semidirected network, dout=din is correct. The in-degree of the node corresponds to the out-degree of the other node.


## 3. Positive Invariant Set

Lemma 1.
System (1) has the positive invariant set:(3)Ω=S1,1;1,I1,1;1,S1,1;2,I1,1;2,…,S1,1;n,I1,1;n,…,SMin,Mout,Mu,IMin,Mout,Mu∈R+:Si,j;n+Ii,j;n≤1.Proof.
Rewrite the above invariant set into the following form:(4)Ω∗=z1,z2,…,z2i+j+n∈R+:zk+l+m+z2i+j+n−k−l−m≤1.
The boundary ofΩ∗ consists of the following three kinds of hyperplanes:(5)Va=z∈Ω∗za=0,Hb=z∈Ω∗zb=0,Qa=z∈Ω∗za+zb=1,where a+b=2Min+Mout+Mu, a,b>0 which have(6)as their outer normal vectors, respectively.
Next, consider system (1), for 0≤i≤Min, 0≤j≤Mout, and 0≤n≤Mu, calculations yield(7)dzdtz∈Va,ϕa=−η1−zb<0,dzdtz∈Hb,φb=−αS−αRzbτdiΘd+τumΘu−αRτdiΘd+τumΘu<0,dzdtz∈Va,ϕa=−γza−βzb<0.
Consider a systemdx/dt=fx which is defined at least in a compact set C. Then, C is invariant if for every point y on ∂C (the boundary of C), the vector fy is tangent to or pointing into C [22, 23]. According to this theorem, we can obtain the conclusion that Ω is positively invariant.


## 4. Equilibria and Basic Reproduction Number

Obviously, system (1) has a disease-free equilibrium Ei,j;n0=η/η+r,0,r/η+r for 0≤i≤Min, 0≤j≤Mout, 0≤n≤Mu. So, the basic reproduction number R0 can be calculated. Basic reproduction number is the expected number of secondary cases produced by a typical infected individual during its entire period of infectiousness in a completely susceptible population [24]. It can be clearly seen that system (1) is a closed system. According to the calculation method given in [25, 26], R0=ρFV−1. After calculation, it can be concluded as follows:(8)V=−β0⋯00−β⋯0⋮⋮⋱⋮00⋯−β=−βE.So,V−1=−1/βE. The matrix F can be expressed in the following form:(9)F=F1,1,11,1,1⋯F1,1,11,1,du⋯F1,1,1din,dout,1⋯F1,1,1din,dout,du⋮⋱⋮⋱⋮⋱⋮F1,1,du1,1,1⋯F1,1,du1,1,du⋯F1,1,dudin,dout,1⋯F1,1,dudin,dout,du⋮⋱⋮⋱⋮⋱⋮Fdin,dout,11,1,1⋯Fdin,dout,11,1,du⋯Fdin,dout,1din,dout,1⋯Fdin,dout,1din,dout,du⋮⋱⋮⋱⋮⋱⋮Fdin,dout,du1,1,1⋯Fdin,dout,du1,1,du⋯Fdin,dout,dudin,dout,1⋯Fdin,dout,dudin,dout,du,where(10)Fk,l;mi,j,n=αSτdkηη+γjPi,j;ndout+αSτumηη+γnPi,j;ndu+αRτdkγη+γjPi,j;ndout+αRτumγη+γnPi,j;ndu=αSη+αRγη+γτdkjPi,j;ndout+αSη+αRγη+γτumnPi,j;ndu=αSη+αRγη+γτdkjPi,j;ndout+τumnPi,j;ndu.Bring the above results to the original matrix.FV−1=−1/βFk,l;mi,j,n. Next, we calculate its spectral radius. After calculation, it is found that the matrix has a similar matrix FS=αSη+αRγ/η+γF∗000. The F∗ is(11)F∗=τudu2duτudinduduτddoutdudoutτddoutdindout,where dindout=∑i,j,nijPi,j;n, dindu=∑i,j,ninPi,j;n, and dudout=∑i,j,njnPi,j;n. The characteristic polynomial of matrix F∗ is(12)λ2−τddindoutdout+τudu2duλ+τuτddoutdudindoutdu2−dindududout=0.The discriminant of the equation is(13)Δ=τddindoutdout+τudu2du2−4τuτddoutdudindoutdu2−dindududout=τddindoutdout−τudu2du2+4τuτddindududoutdoutdu>0.So, the matrixF∗ has two unequal real roots. And the basic reproduction number of the model is(14)R0=αSη+αRγ2βη+γτddindoutdout+τudu2du+Δ,where Δ=τddindout/dout−τudu2/du2+4τuτddindududout/doutdu.Lemma 2.
IfR0>1, then exists a endemic equilibrium E∗Si,j;n∗,Ii,j;n∗,Ri,j;n∗, where(15)Si,j;n∗=η1−Ii,j;n∗γ+η+αSτdiΘ∗d+τunΘ∗u,Ri,j;n∗=β−γIi,j;n∗+γγ+η+αRτdiΘ∗d+τunΘ∗u,Ii,j;n∗=1−ββ+αRτdiΘ∗d+τunΘ∗u+ηαS−αRτdiΘ∗d+τunΘ∗u/γ+η+αSτdiΘ∗d+τunΘ∗u,Θ∗d=1dout∑jPi,j;nIi,j;n∗,Θ∗u=1du∑nPi,j;nIi,j;n∗.


## 5. Dynamical Analysis of the Model

Theorem 1.
IfR0<1, the disease-free equilibrium E0 of system (1) is stable. If R0>1, E0 is unstable.Theorem 2.
IfR0<1, the disease-free equilibrium E0 of system (1) is globally asymptotically stable.Proof.
In the feasible region, we consider the first equation of system (1):(16)dSi,j;ntdt=−αSτdiSi,j;ntΘd+τunSi,j;ntΘu+ηRi,j;nt−γSi,j;nt.
Obviously,dSi,j;n/dt≤η−η+γSi,j;n. So, we can obtain limt⟶∞supSi,j;n≤η/η+γ=Si,j;n0 such that Si,j;n≤Si,j;n0+ξ1. Next, we analyze the second equation of system (1). When t>T1,(17)dIi,j;ndt≤αS−αRτdiΘd+τunΘuSi,j;n0+ξ1−βIi,j;n+αRτdiΘd+τunΘu1−Ii,j;n≤αS−αRτdiΘd+τunΘuSi,j;n0+ξ1−βIi,j;n+αRτdiΘd+τunΘu,dΘd+Θudt=dΘddt+dΘudt=1dout∑jPi,j;ndIi,j;ndt+1du∑nPi,j;ndIi,j;ndt=∑jdout+nduPi,j;ndIi,j;ndt=∑jdout+nduPi,j;nαS−αRτdiΘd+τunΘuSi,j;n0+ξ1+∑​jdout+nduPi,j;nαRτdiΘd+τunΘu−βIi,j;n=∑jdout+nduPi,j;nαS−αRτdiΘd+τunΘuSi,j;n0+ξ1+∑jdout+nduPi,j;nαRτdiΘd+τunΘu−∑jdout+nduPi,j;nβIi,j;n=∑jdout+nduPi,j;nαS−αRSi,j;n0+ξ1+αRτdiΘd+τunΘu−2βΘd+Θu.
In this,∑j/dout+n/duPi,j;nαS−αRSi,j;n0+ξ1+αRτdiΘd+τunΘu can be divided into the following four parts:
The first part:(18)Θdτd∑ijdoutPi,j;nαS−αRSi,j;n0+ξ1+αR=Θdτd∑​ijdoutPi,j;nαS−αRηη+γ+αR+αS−αRξ1=Θdτd∑ijdoutPi,j;nηη+γαS+γη+γαR+αS−αRξ1=ΘdηαS+γαRη+γτd∑​ijdoutPi,j;n+ε1=ΘdηαS+γαRη+γτddindoutdout+ε1.
The second part:(19)Θuτu∑n2duPi,j;nαS−αRSi,j;n0+ξ1+αR=Θuτu∑n2duPi,j;nηαS+γαRη+γ+αS−αRξ1=ΘuηαS+γαRη+γτudu2du+ε2.
The third part:(20)Θdτd∑induPi,j;nαS−αRSi,j;n0+ξ1+αR=Θdτd∑induPi,j;nηαS+γαRη+γ+αS−αRξ1=ΘdηαS+γαRη+γτddindudu+ε3.
The fourth part:(21)Θuτu∑jnduPi,j;nαS−αRSi,j;n0+ξ1+αR=Θuτu∑jnduPi,j;nηαS+γαRη+γ+αS−αRξ1=ΘuηαS+γαRη+γτudoutdudout+ε4.
So,(22)dΘd+Θudt=dΘddt+dΘudt=ηαS+γαRη+γΘdτddindoutdout+Θuτudu2du+Θdτddindudu+Θuτudoutdudout+ε−2βΘd+Θu=2βΘd+ΘuR0+ε∗−1,where ε∗=ε1+ε2+ε3+ε4. In order to prove that the last step of the above formula is established, it proves that Θd+ΘuR0−4parts=0 is established. Next, we prove it.(23)Θd+ΘuR0−4parts=ηαS+γαRη+γΘd+Θuτddindoutdout+τudu2du+Δ−Θdτddindoutdout−Θuτudu2du−Θdτddindudu−Θuτudoutdudout=Θuτddindoutdout+Θdτudu2du+Θd+ΘuΔ−Θdτddindudu−Θuτudoutdudout.
To prove that the formula is 0. That is, satisfy(24)Θuτddindoutdout+Θdτudu2du−Θdτddindudu−Θuτudoutdudout2=Θd+Θu2Δ.
After simplification, we can obtain(25)dindoutduτd−doutdu2τu2−dindudoutτd−dudoutduτu∗=0.
According to abovementioned equation,(26)dindoutduτd−doutdu2τu2−dindudoutτd−dudoutduτu=0.
SodΘd+Θu/dt≤2βΘd+ΘuR0+ε∗−1 is correct. It guarantees that if R0<1, dΘd+Θu/dt<0. Hence, limt⟶+∞Θd+Θu=0 and limt⟶+∞Ii,j;n=0 for all 1≤i≤Min, 1≤i≤Mout, and 1≤n≤Mu. For any ξ2, there exist T2>0 such that 0≤Ii,j;n≤ξ2. According to system (1), we can obtain(27)dSi,j;ndt≥η−η+γSi,j;n−η+αSτdiSi,j;nΘd+τunSi,j;nΘuξ2.
Settingξ2⟶0, then dSi,j;n/dt≥η−η+γSi,j;n. So, limt⟶∞supSi,j;n≥η/η+γ=Si,j;n0. And, we have obtained limt⟶∞supSi,j;n≤η/η+γ=Si,j;n0. According to limt⟶∞Si,j;n=η/η+γ=Si,j;n0 and limt⟶∞Ii,j;n=0=Ii,j;n0 for 0≤i≤Min, 0≤i≤Mout, and 0≤n≤Mu, the disease-free equilibrium E0 is globally asymptotically stable.Theorem 3.
IfR0>1, the endemic equilibrium E∗ of system (1) is globally asymptotically stable provided that 2γη+β>γ+η−β.Proof.
Let us consider the following four positive definite functions, which is defined along a solution of system (1):(28)VSi,j;n=12Si,j;n−Si,j;n∗2,VRi,j;n=12Ri,j;n−Ri,j;n∗2,VΘd=Θd−Θ∗d−Θ∗dlnΘdΘ∗d,VΘu=Θu−Θ∗u−Θ∗ulnΘuΘ∗u.
The following solves the function for the differential form:(29)dVSi,j;ndt=Si,j;n−Si,j;n∗dSi,j;ndt=Si,j;n−Si,j;n∗−αSτdiSi,j;nΘd+τunSi,j;nΘu+ηRi,j;n−γSi,j;n=Si,j;n−Si,j;n∗−αSτdiSi,j;nΘd−Si,j;n∗Θ∗d−αSτunSi,j;nΘu−Si,j;n∗Θ∗u+ηRi,j;n−Ri,j;n∗−γSi,j;n−Si,j;n∗=γSi,j;n−Si,j;n∗2+ηRi,j;n−Ri,j;n∗Si,j;n−Si,j;n∗−αSτdiΘdSi,j;n−Si,j;n∗2−αSτunΘuSi,j;n−Si,j;n∗2−αSτdiSi,j;n∗Si,j;n−Si,j;n∗Θd−Θ∗d−αSτunSi,j;n∗Si,j;n−Si,j;n∗Θu−Θ∗u,dVRi,j;ndt=Ri,j;n−Ri,j;n∗dRi,j;ndt=Ri,j;n−Ri,j;n∗−αRτdkRi,j;nΘd+τumRi,j;nΘu−ηRi,j;n+γSi,j;n+βIi,j;n=Ri,j;n−Ri,j;n∗−αRτdiRi,j;nΘd+τunRi,j;nΘu−η+βRi,j;n+γ−βSi,j;n+β=Ri,j;n−Ri,j;n∗−αRτdiRi,j;nΘd−Ri,j;n∗Θ∗d−αRτunRi,j;nΘu−Ri,j;n∗Θ∗u−η+βRi,j;n−Ri,j;n∗+γ−βSi,j;n−Si,j;n∗=−η+βRi,j;n−Ri,j;n∗2+γ−βSi,j;n−Si,j;n∗Ri,j;n−Ri,j;n∗−αRτdiΘdRi,j;n−Ri,j;n∗2−αRτunΘuRi,j;n−Ri,j;n∗2−αRτdiRi,j;n∗Ri,j;n−Ri,j;n∗Θd−Θ∗d−αRτunRi,j;n∗Ri,j;n−Ri,j;n∗Θu−Θ∗u,dVΘddt=Θd−Θ∗dΘddΘddt=Θd−Θ∗dΘd1dout∑jPi,j;nαSτdiSi,j;nΘd+τunSi,j;nΘu+Θd−Θ∗dΘd1dout∑jPi,j;nαRτdiRi,j;nΘd+τunRi,j;nΘu−βIi,j;n=Θd−Θ∗dΘd1dout∑jPi,j;nαSτdiSi,j;nΘd+τunSi,j;nΘu+Θd−Θ∗dΘd1dout∑jPi,j;nαRτdiRi,j;nΘd+τunRi,j;nΘu−βΘd=Θd−Θ∗d1dout∑jPi,j;nαSτdiSi,j;n+τunSi,j;nΘuΘd+Θd−Θ∗d1dout∑jPi,j;nαRτdiRi,j;n+τunRi,j;nΘuΘd−β=1dout∑jPi,j;nαSτdiSi,j;n−Si,j;n∗Θd−Θ∗d+αSτunSi,j;nΘuΘd−Si,j;n∗Θ∗uΘ∗dΘd−Θ∗d+αSτdiRi,j;n−Ri,j;n∗Θd−Θ∗d+αSτunRi,j;nΘuΘd−Ri,j;n∗Θ∗uΘ∗dΘd−Θ∗d.
For the same reason,(30)dVΘudt=1du∑nPi,j;nαSτdiSi,j;nΘdΘu−Si,j;n∗Θ∗dΘ∗uΘu−Θ∗u+αSτunSi,j;n−Si,j;n∗Θu−Θ∗u+αSτdiRi,j;nΘdΘu−Ri,j;n∗Θ∗dΘ∗u+αSτunRi,j;n−Ri,j;n∗Θu−Θ∗u.
Let us consider the following Lyapunov function:(31)V=∑jdout+nduVsi,j;n+∑jdout+nduVSi,j;n+VΘd+VΘu,dVdt=∑jdout+nduγSi,j;n−Si,j;n∗2+ηRi,j;n−Ri,j;n∗Si,j;n−Si,j;n∗−∑​jdout+nduαSτdiΘdSi,j;n−Si,j;n∗2+αSτunΘuSi,j;n−Si,j;n∗2+∑​jdout+ndu−η+βRi,j;n−Ri,j;n∗2+γ−βSi,j;n−Si,j;n∗Ri,j;n−Ri,j;n∗−∑​jdout+nduαRτdiΘdRi,j;n−Ri,j;n∗2+αRτunΘuRi,j;n−Ri,j;n∗2+A,where(32)A=1dout∑jPi,j;nαSτunSi,j;nΘuΘd−Si,j;n∗Θ∗uΘ∗dΘd−Θ∗d+1dout∑jPi,j;nαRτunRi,j;nΘuΘd−Ri,j;n∗Θ∗uΘ∗dΘd−Θ∗d+1du∑nPi,j;nαSτdiSi,j;nΘdΘu−Si,j;n∗Θ∗dΘ∗uΘu−Θ∗u+1du∑nPi,j;nαRτdiRi,j;nΘdΘu−Ri,j;n∗Θ∗dΘ∗uΘu−Θ∗u−1dout∑jPi,j;nαSτunSi,j;n−Si,j;n∗Θu−Θ∗u−1dout∑jPi,j;nαRτunRi,j;n−Ri,j;n∗Θu−Θ∗u−1du∑nPi,j;nαSτdiSi,j;n−Si,j;n∗Θd−Θ∗d−1du∑nPi,j;nαRτdiRi,j;n−Ri,j;n∗Θd−Θ∗d.
We simplify formulaA, and we combine the first and fifth items to get the following results:(33)A1=1dout∑jPi,j;nαSτunSi,j;nΘuΘd−Si,j;n∗Θ∗uΘ∗dΘd−Θ∗d−1dout∑jPi,j;nαSτunSi,j;n−Si,j;n∗Θu−Θ∗u=1dout∑jPi,j;nαSτunSi,j;n∗Θu+Si,j;nΘ∗u−Si,j;n∗Θ∗uΘ∗dΘd−Si,j;nΘuΘdΘ∗d=1dout∑jPi,j;nαSτun1ΘdΘ∗dΘdΘ∗u−ΘuΘ∗dSi,j;nΘ∗d−Si,j;n∗Θd.
According to the definition of the model, it is not difficult to findΘdΘ∗u−ΘuΘ∗d=0. So, A1=0. A=0 is correct for the same reason. So(34)dVdt≤∑jdout+ndu−γSi,j;n−Si,j;n∗2+ηRi,j;n−Ri,j;n∗Si,j;n−Si,j;n∗+∑jdout+ndu−η+βRi,j;n−Ri,j;n∗2+γ−βSi,j;n−Si,j;n∗Ri,j;n−Ri,j;n∗=∑jdout+ndu−γSi,j;n−Si,j;n∗2+η+γ−βRi,j;n−Ri,j;n∗Si,j;n−Si,j;n∗−η+βRi,j;n−Ri,j;n∗2.
When,2γη+β>γ+η−β, V′=0 if and only if Si,j;n=Si,j;n∗, Ii,j;n=Ii,j;n∗, and Ri,j;n=Ri,j;n∗. According to the Lyapunov theorem [27] and the LaSalle’s invariant principal [28], we can conclude that the endemic equilibrium E∗ of system (1) is globally asymptotically stable.


## 6. Numerical Simulation and Analysis

In this section, we give some numerical simulation and analysis to verify the theorems which is obtained in Sections4 and 5. All the simulations are based on the semidirected networks. The way of building the semidirected networks is as follows: First, we build a scale-free network. We randomly select the connected edges between nodes to become directed, and randomly specify the direct of connection to ensure the generality of the network. To fit the actual situation, the number of the undirected connection is much more than the directed connection. None of the nodes in this numerical simulation are isolated nodes. The number of all nodes on the semidirected network is 500.For the case that ifR0<1, Theorem 1 and Theorem 2 show that the disease-free E0 is locally and globally asymptotically stable. We will make use of the semidirected network to verify the correctness of it.Example 1.
The parameters in system (1) are taken as αS=0.2, αR=0.1, τd=0.02, τu=0.03, γ=0.1, β=0.1, and η=0.3. The semidirected network that is built randomly has the following property. din=1.536, dout=1.536, du=4.78, dindout=2.680, du2=43.492, dudout=10.282, and dindu=12.936. So the basic reproduction number R0=0.546<1. The disease-free equilibrium E0 is globally asymptotically stable. Therefore, the diseases will disappear in the semidirected network.
From Figure1 it is clear that nodes with degree (2, 10; 12) will recover from diseases. And all of the nodes will recover from the diseases. This is to say the entire semidirected network will not be affected by the disease.
Next, we observe an evolution of two different kinds of nodes. The degree of one kind node is (2, 6; 7). It is a kind of node with fewer connections. Another kind of nodes with degree (8, 22; 46) is the centre node in the semidirected network.
The first figure in Figure2 is the evolution of the centre node with degree (8, 22; 46), and the second is the node with degree (2, 6; 7). It is easy to find that the node with a different degree will reach the same situation. It fits the conclusion that all nodes will reach the disease-free equilibrium Ei,j;n0=η/η+r,0,r/η+r for 0≤i≤Min, 0≤j≤Mout, and 0≤n≤Mu. It is clear that the node with bigger degree has a greater magnitude of change than the node with fewer degree. The infectious rate of the node with degree (8, 22; 46) increases at the beginning, and the centre node is more susceptible to disease transmission. The disease propagates in the seminetwork mainly through nodes with bigger degree because they have enough edges which connect to the whole network. Finally, it becomes 0. And the disease disappears on the semidirected networks. When disease disappears, nodes in state S and state R will meet a balance according to the stability of the two-dimensional system. And in Figure 3, we will show the infectious evolution of the whole nodes. It meets our conclusion. The node with small degree in state I will disappear not increase like the nodes with big degree in the beginning. All the nodes will not be affected by the disease. So the disease disappears.
Example 2. The parameters in system (1) are taken as αS=0.2, αR=0.1, τd=0.15, τu=0.2, γ=0.4, β=0.2, and η=0.2. The semidirected network that is rebuilt randomly has the following property: din=1.536, dout=1.536, du=4.78, dindout=2.680, du2=43.492, dudout=10.282, and dindu=12.936. So the basic reproduction number R0=1.430>1. The epidemic equilibrium E∗ is globally asymptotically stable. And the disease will exist on the semidirected networks forever, and the number of infectious nodes is invariable.
From Figure4, we can find that the node with degree will be in a balanced state. The number of nodes in the infected state is invariable. All the nodes are in the dynamic equilibrium state. The number of nodes in the self-protection state increases in the beginning, and it will fall down after experiencing the highest peak, because the infectious diseases have not spread in the network and a large number of the nodes are in the state S. Almost all nodes with different degrees have such a trend.
Next, we compare three kinds of nodes with different degrees. One is (8, 22; 32) and another is (3, 8; 9). One node has a large degree, and it belongs to the central node. Other has the middle degree on the semidirected network. The second one is the node with degree (1, 2; 5). The results are shown in Figure5. It has the smaller degree, and it belongs to the edge node.
For three different degrees of nodes, the trend of change is almost the same. The difference is that the value of the final balance point is different. It is clear that the node with great degree will have the higher rate of the infected in the end. The result can be seen in expression ofIi,j;n∗ due to its strong connectivity. Higher degree may result in increase in the value of Ii,j;n∗ increases. The greater the degree of the node, the greater the probability of being infected. It means that if we control the rate of the infected nodes with bigger degree, the rate of the infected nodes with smaller degree will decrease. In this way, we can control the disease propagation in the semidirected network. As long as we control the nodes with high degrees, we can control the spread of the entire network. In Figure 6, we will give evolution of the whole nodes. In this figure, I=∑Pi,j;nIi,j;n, S=∑Pi,j;nSi,j;n, and R=∑Pi,j;nRi,j;n.
It is clear that the nodes on the semidirected networks reach equilibrium after a period of time. Its changing trace is similar to the figures given earlier. Disease will exist in the network forever. The state of the node reaches dynamic balance. Next, we rebuild the semidirected network with 1000 nodes, and the parameters are unchanged. The simulation is shown in Figure7.
There is no significant change between Figures6 and 7. It means that the increase of the node number affects the final result, unless changing the method of connection. Next, we present Figure 8 to show R0 under different τd and τu. It meets that the semidirected networks which we build has large undirected connection and few directed connection.Figure 1
Evolution of nodes with degree (2, 10; 12) on semidirected networks.Figure 2
Evolution of the nodes of different degrees.
(a)
(b)Figure 3
Evolution of the whole nodes.Figure 4
Evolution of nodes with degree (1, 2; 5) on semidirected networks.Figure 5
Evolution of the nodes of different degrees.
(a)
(b)Figure 6
Evolution of the whole nodes.Figure 7
Evolution of the whole nodes (1000 nodes).Figure 8
R
0 under different τd and τu.


## 7. Conclusion and Discussion

In this paper, we study the epidemic model with no full immunity on semidirected networks. In this model, we use self-protection to replace the immunity which is used in the past. Based on this model, the basic reproduction numberR0 can be used to control infectious diseases. We mainly analyze stability of disease-free equilibrium E0 and endemic equilibrium E∗. If R0<1, the disease-free equilibrium E0 is locally and globally asymptotically stable. It has been proved. And endemic equilibrium E∗ is globally asymptotically stable under some condition. It means that the infectious diseases will exist on the semidirected networks in a long time.According to the basic reproduction numberR0 that we calculated, the propagation of diseases is affected by many factors. And, it is also affected by the degree distribution of the networks. We can consider more about the recovery η. If the recovery is increasing, the diseases can be controlled effectively.When we prove the globally asymptotically stable, we give some condition on it. We cannot investigate the globally asymptotical stability only whenR0>1. We hope to improve in the future.



---

*Source: 8792497-2019-10-23.xml*