# Modeling and Stability Analysis of Worm Propagation in Wireless Sensor Network

Authors: Liping Feng; Lipeng Song; Qingshan Zhao; Hongbin Wang
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

An improved SIRS model considering communication radius and distributed density of nodes is proposed. The proposed model captures both the spatial and temporal dynamics of worms spread process. Using differential dynamical theories, we investigate dynamics of worm propagation to time in wireless sensor networks (WSNs). Reproductive number which determines global dynamics of worm propagation in WSNs is obtained. Equilibriums and their stabilities are also found. If reproductive number is less than one, the infected fraction of the sensor nodes disappears and if the reproduction number is greater than one, the infected fraction asymptotically stabilizes at the endemic equilibrium. Based on the reproduction number, we discuss the threshold of worm propagation about communication radius and distributed density of nodes in WSNs. Finally, numerical simulations verify the correctness of theoretical analysis.


---

## Body

## 1. Introduction

A sensor network is composed of hundreds or even thousands of sensor nodes that are allowed random deployment in inaccessible terrains or disaster relief operations [1]. Wireless sensor networks (WSNs), as a kind of new information and communication network, have gained worldwide attention owing to their potential in civil and military applications, for instance, intrusion detection, perimeter monitoring, information gathering, and smart logistics support in an unknown deployed area [2–4].With widespread applications of WSNs, research on WSNs has been a hot topic. Some methods have been proposed for prolonging the lifetime of WSNs focusing on energy consumption [5–7], device placement [8], and topology management [9]. Because sensor nodes are constrained sources, they have weak defenses and are attacking targets for worms. Injecting malware into some nodes has become a serious threat [10]. Recently, malicious codes targeting wireless devices have emerged, which can spread directly from device to device using wireless communication technology, such as Wi-Fi and Bluetooth [11–14]. For instance, computer worm like Cabir uses the Bluetooth interface to spread among cellphones, which means that worms have committed the wireless domain and WSNs are also extremely vulnerable to malware.Actions of malicious objects on the Internet have been studied by using epidemical models and have provided insights for controlling worm prevalence in networks [15–20]. In [15], the authors presented an E-mail virus model that accounts for behaviors of E-mail users and analyzed propagation features of E-mail viruses in different network topologies. In [16–20], the authors proposed epidemic models with time delay and analyzed dynamical features of worm prevalence. To effectively defend against worm intrusions, it is necessary to deeply understand dynamical features of worm propagation in WSNs. Existing research prove that epidemical models are valuable for portraying characteristics of worm propagation. Since there is a basic similarity between worm propagation through wireless devices and traditional worm spread on the Internet, the epidemical models extensively are applied to study worm spread in WSNs by some researchers in recent years [11, 21–23]. In [21], the authors proposed a SIRS malware propagation model with feedback controller and analyzed Hopf bifurcation dynamics of malware prevalence in mobile wireless networks. The authors in [22] presented an epidemic model with vaccination compartment which captures both the spatial and temporal dynamics of worm spread process, and some mathematical analyses and numerical simulations were performed based on this model. The common problem of the above models is that the characteristics of wireless sensor networks like energy consumption, communication radius, and distributed density of nodes have not been considered in models. In [11], the authors developed a mathematical model for the propagation that incorporates important parameters derived from the communication patterns of the protocol under test. Based on this model, the authors analyzed the propagation rate and the extent of spread of a malware over typical broadcast protocols proposed in the paper. Wang and Li derived an iSIR model describing the process of worm propagation with energy consumption of nodes in WSNs [23]. Numerical simulations are performed to observe the effects of the network topology and energy consumption of nodes on worm spread in WSNs. However, the authors have not performed mathematical analyses based on this model. In fact, key parameters of affecting worm spread can be found by explicit mathematical analyses.To better portrait the features of worm propagation in WSNs, in this paper, we study the attacking behavior of possible worms in WSNs by constructing an improved SIRS epidemic model. In this model, the following three factors are considered: (i) energy consumptions of nodes; (ii) communication radius of nodes; and (iii) distributed density of nodes in WSNs. Based on this model, we analyze the stability of worm prevalence through finding the equilibriums of model.The rest of this paper is organized as follows: in Section2, we analyze topology of WSNs and present the model formulation. Section 3 derives the equilibriums of the model and discusses the stability of worm propagation at the equilibriums. In Sections 4 and 5, numerical simulations are performed to verify the correctness of theoretical analyses and some conclusions are given, respectively.


## 2. The Proposed Model

### 2.1. System Description

We model a wireless sensor network composed ofN nodes. The nodes are uniformly distributed in L
×
L area (nodes average density is ρ
=
N
/
L
2) and the wireless communication range of every node is r. The topological structure of a WSN is shown in Figure 1.Figure 1
The topology structure of a WSN.Based on the existingS
I
R epidemic model [24, 25], the nodes in WSNs are classified into three states:(i)
Susceptible state (S): nodes in S have not been infected by worms, and these nodes are vulnerable to worms in WSNs.
(ii)
Infected state (I): nodes have been infected by worms and have the ability to infect other nodes in WSNs.
(iii)
Recover state (R): nodes have installed a detection tool that can identify and remove worms, or nodes have installed a software patch to eliminate the node vulnerability exploited by worms.We consider the following state transitions among these three states:(i)
Users may immunize their nodes with countermeasures in statesS and I with probabilities ω and γ, respectively.
(ii)
As the energy of nodes is exhausted, some nodes become dead nodes with probabilityμ.
(iii)
Infected nodesI infect susceptible S with effective infection rate β.
(iv)
Some recovered nodes may become susceptible nodes with probabilityε.Transition relationships among node states are described in Figure 2.Figure 2
Transition relationship of states of nodes.


### 2.2. Model Derivation

The communication area of a node is denoted byS
r, and the density of susceptible nodes in a unit area in WSNs is denoted by p
t. Then the following equations hold:(1)
p
t
=
S
t
L
2
,
S
r
=
π
r
2
,
S
′
t
=
S
r
p
t
.From (1), we can get(2)
S
′
t
=
π
r
2
S
t
L
2
.According to state transition relationships in Figure2, the mathematical model of worm propagation in WSNs can be derived as follows:(3)
d
S
t
d
t
=
μ
N
-
π
r
2
L
2
β
S
t
I
t
-
μ
+
ω
S
+
ε
R
,
d
I
t
d
t
=
π
r
2
L
2
β
S
t
I
t
-
μ
+
γ
I
,
d
R
t
d
t
=
ω
S
+
γ
I
-
μ
+
ε
R
.For convenience, let(4)
ξ
=
π
r
2
L
2
β
.Then, system (3) can be written as(5)
d
S
d
t
=
μ
N
-
ξ
S
t
I
t
-
μ
+
ω
S
+
ε
R
,
d
I
d
t
=
ξ
S
t
I
t
-
μ
+
γ
I
,
d
R
d
t
=
ω
S
+
γ
I
-
μ
+
ε
R
.


## 2.1. System Description

We model a wireless sensor network composed ofN nodes. The nodes are uniformly distributed in L
×
L area (nodes average density is ρ
=
N
/
L
2) and the wireless communication range of every node is r. The topological structure of a WSN is shown in Figure 1.Figure 1
The topology structure of a WSN.Based on the existingS
I
R epidemic model [24, 25], the nodes in WSNs are classified into three states:(i)
Susceptible state (S): nodes in S have not been infected by worms, and these nodes are vulnerable to worms in WSNs.
(ii)
Infected state (I): nodes have been infected by worms and have the ability to infect other nodes in WSNs.
(iii)
Recover state (R): nodes have installed a detection tool that can identify and remove worms, or nodes have installed a software patch to eliminate the node vulnerability exploited by worms.We consider the following state transitions among these three states:(i)
Users may immunize their nodes with countermeasures in statesS and I with probabilities ω and γ, respectively.
(ii)
As the energy of nodes is exhausted, some nodes become dead nodes with probabilityμ.
(iii)
Infected nodesI infect susceptible S with effective infection rate β.
(iv)
Some recovered nodes may become susceptible nodes with probabilityε.Transition relationships among node states are described in Figure 2.Figure 2
Transition relationship of states of nodes.


## 2.2. Model Derivation

The communication area of a node is denoted byS
r, and the density of susceptible nodes in a unit area in WSNs is denoted by p
t. Then the following equations hold:(1)
p
t
=
S
t
L
2
,
S
r
=
π
r
2
,
S
′
t
=
S
r
p
t
.From (1), we can get(2)
S
′
t
=
π
r
2
S
t
L
2
.According to state transition relationships in Figure2, the mathematical model of worm propagation in WSNs can be derived as follows:(3)
d
S
t
d
t
=
μ
N
-
π
r
2
L
2
β
S
t
I
t
-
μ
+
ω
S
+
ε
R
,
d
I
t
d
t
=
π
r
2
L
2
β
S
t
I
t
-
μ
+
γ
I
,
d
R
t
d
t
=
ω
S
+
γ
I
-
μ
+
ε
R
.For convenience, let(4)
ξ
=
π
r
2
L
2
β
.Then, system (3) can be written as(5)
d
S
d
t
=
μ
N
-
ξ
S
t
I
t
-
μ
+
ω
S
+
ε
R
,
d
I
d
t
=
ξ
S
t
I
t
-
μ
+
γ
I
,
d
R
d
t
=
ω
S
+
γ
I
-
μ
+
ε
R
.


## 3. Stability Analysis of Equilibriums

In this section, we will find the equilibriums of system (5) and investigate their stability. The equilibriums of system (5) are given by solutions of(6)
μ
N
-
ξ
S
t
I
t
-
μ
+
ω
S
+
ε
R
=
0
,
ξ
S
t
I
t
-
μ
+
γ
I
=
0
,
ω
S
+
γ
I
-
μ
+
ε
R
=
0
.Let ξ
S
t
I
t
-
μ
+
γ
I
=
0 (the second equation of (6)); we have I
∗
=
0 or I
∗
>
0 and S
∗
=
(
μ
+
γ
)
/
ξ.  For the case of I
∗
=
0, we have worm-free equilibrium (7)
Q
0
=
S
0
∗
,
I
0
∗
,
R
0
∗
=
μ
+
ε
μ
+
ε
+
ω
N
,
0
,
ω
S
μ
+
ε
.For the case of I
∗
>
0, we have endemic equilibrium(8)
Q
∗
=
S
1
∗
,
I
1
∗
,
R
1
∗
=
μ
+
γ
ξ
,
ξ
μ
+
ε
N
-
μ
+
γ
μ
+
ε
+
ω
ξ
μ
+
ε
+
γ
,
ω
S
1
∗
+
γ
I
1
∗
μ
+
ε
.Let(9)
R
0
=
ξ
μ
+
ε
μ
+
γ
μ
+
ε
+
ω
N
.Notably, the endemic equilibrium is meaningful only if R
0
>
1.


### 3.1. Worm-Free Equilibrium and Its Stability

Lemma 1.
The worm-free equilibrium is locally asymptotically stable ifR
0
<
1 and unstable if R
0
>
1.Proof.
According toQ
0
=
(
(
μ
+
ε
/
μ
+
ε
+
ω
)
N
,
0
,
ω
S
/
(
μ
+
ε
)
), the characteristic equation of system (5) at worm-free equilibrium Q
0 is (10)
det
⁡
-
μ
+
ω
-
λ
-
ξ
S
0
∗
ε
0
ξ
S
-
μ
+
γ
-
λ
0
ω
γ
μ
+
ε
-
λ
=
0
,which is equivalent to(11)
ξ
S
0
∗
-
μ
+
γ
-
λ
λ
2
-
2
μ
+
ε
+
ω
λ
+
μ
2
+
ε
μ
+
μ
ω
=
0
.Equation (11) has a characteristic root λ
1
=
ξ
S
0
∗
-
μ
+
γ
=
μ
+
γ
R
0
-
1 and the roots of equation(12)
λ
2
+
2
μ
+
ε
+
ω
λ
+
μ
2
+
ε
μ
+
μ
ω
=
0
.
Obviously, in accordance with the relationship between roots and coefficients of quadratic equation, there is no positive real part characteristic root of (12). Hence, when R
0
<
1, (11) has no positive real root, and worm-free equilibrium Q
0 is locally asymptotically stable. When R
0
>
1, (11) has a positive root; thus, worm-free equilibrium Q
0 is an unstable saddle-point.Furthermore, the following theorem holds.Theorem 2.
The worm-free equilibrium is globally asymptotically stable ifR
0
≤
1
.Proof.
From the first equation of system (5)(13)
S
˙
t
≤
μ
+
ε
N
-
μ
+
ε
+
ω
S
.Thus, S
t
≤
μ
+
ε
N
/
μ
+
ε
+
ω
+
S
0
-
(
μ
+
ε
)
/
μ
+
ε
+
ω
e
x
p
-
μ
+
ε
+
ω
t. When t
→
∞, we obtain(14)
S
t
≤
μ
+
ε
μ
+
ε
+
ω
N
.Consider a Lyapunov function(15)
L
t
=
E
t
,
L
˙
t
=
ξ
S
t
I
t
-
μ
+
γ
I
t
≤
ξ
μ
+
ω
μ
+
ω
+
ε
N
-
μ
+
γ
I
t
=
1
μ
+
γ
R
0
-
1
I
t
≤
0
.So, we prove that worm-free equilibrium Q
0 is globally asymptotically stable.


### 3.2. Endemic Equilibrium and Its Stability

Now, we investigate the local stability of endemic equilibriumQ
∗. The characteristic equation of system (5) at endemic equilibrium Q
∗ is (16)
det
⁡
-
ξ
I
1
∗
-
μ
+
ω
-
λ
-
ξ
S
1
∗
ε
ξ
I
1
∗
ξ
S
1
∗
-
μ
+
γ
-
λ
0
ω
γ
μ
+
ε
-
λ
=
0
,which is equivalent to(17)
λ
3
+
p
0
λ
2
+
p
1
λ
+
p
2
=
0
,where p
0
=
3
μ
+
ε
+
ω
+
(
μ
+
ε
/
μ
+
ε
+
γ
)
N, p
1
=
μ
+
ε
2
μ
+
γ
+
(
μ
+
ε
/
μ
+
ε
+
γ
)
N
+
μ
ω
+
ξ
2
S
1
∗
I
1
∗, and p
2
=
μ
ξ
μ
+
ε
I
1
∗
-
S
1
∗
+
μ
μ
+
γ
μ
+
ε
+
ω
+
μ
ξ
γ
I
1
∗
-
ω
S
1
∗.Obviously,p
0
>
0, p
1
>
0, p
2
>
0, and p
0
p
1
-
p
2
>
0. According to the theorem of Routh-Hurwitz [26, 27], it follows that the roots of (17) have negative real parts. Therefore, the endemic equilibrium is locally asymptotically stable.From the above discussion, we can summarize the following conclusion.Lemma 3.
IfR
0
>
1, then endemic equilibrium is locally asymptotically stable.Note that the number of nodes in WSNs is relatively stable; that is, at timet, the number of nodes S
t, I
t, and R
t in states S
,
I, and R, respectively, satisfies(18)
S
t
+
I
t
+
R
t
=
N
.Hence, the dynamics of system (5) is equivalent to the following system:(19)
d
S
t
d
t
=
μ
+
ω
N
-
ξ
S
t
I
t
-
μ
+
ω
+
ε
S
t
-
ε
I
t
,
d
I
t
d
t
=
ξ
S
t
I
t
-
μ
+
γ
I
t
.Clearly, when R
0
>
1, system (19) has a unique positive equilibrium Q
∗
S
1
∗
,
I
1
∗, where (20)
S
1
∗
,
I
1
∗
=
μ
+
γ
ξ
,
ξ
μ
+
ε
N
-
μ
+
γ
μ
+
ε
+
ω
ξ
μ
+
ε
+
γ
.Now, we state and prove a result on the global stability of the endemic equilibrium S
1
∗
,
I
1
∗ of system (19).Theorem 4.
WhenR
0
>
1, the endemic equilibrium S
1
∗
,
I
1
∗ of system (19) is globally asymptotically stable.Proof.
Consider the following Lyapunov function [28]:(21)
L
t
=
∫
S
1
∗
S
x
-
S
1
∗
x
d
x
+
∫
I
1
∗
I
x
-
I
1
∗
x
d
x
.The time derivative of L
t along the solution of system (19) is given by(22)
L
˙
t
=
S
-
S
1
∗
S
S
′
+
I
-
I
1
∗
I
I
′
=
1
-
S
1
∗
S
μ
+
ε
N
-
ξ
S
t
I
t
-
μ
+
ε
+
ω
S
t
-
ε
I
t
+
1
-
I
1
∗
I
ξ
S
t
I
t
-
μ
+
γ
I
t
≤
1
-
S
1
∗
S
μ
+
ε
N
-
ξ
S
t
I
t
-
μ
+
ε
+
ω
S
t
+
1
-
I
1
∗
I
ξ
S
t
I
t
-
μ
+
γ
I
t
=
-
μ
+
ε
N
S
S
1
∗
S
1
∗
S
-
1
2
≤
0
.The proof is completed.Remark 5.
By Theorem4, we obtain that R
0
>
1; then endemic equilibrium Q
∗ of system (5) is globally asymptotically stable.


## 3.1. Worm-Free Equilibrium and Its Stability

Lemma 1.
The worm-free equilibrium is locally asymptotically stable ifR
0
<
1 and unstable if R
0
>
1.Proof.
According toQ
0
=
(
(
μ
+
ε
/
μ
+
ε
+
ω
)
N
,
0
,
ω
S
/
(
μ
+
ε
)
), the characteristic equation of system (5) at worm-free equilibrium Q
0 is (10)
det
⁡
-
μ
+
ω
-
λ
-
ξ
S
0
∗
ε
0
ξ
S
-
μ
+
γ
-
λ
0
ω
γ
μ
+
ε
-
λ
=
0
,which is equivalent to(11)
ξ
S
0
∗
-
μ
+
γ
-
λ
λ
2
-
2
μ
+
ε
+
ω
λ
+
μ
2
+
ε
μ
+
μ
ω
=
0
.Equation (11) has a characteristic root λ
1
=
ξ
S
0
∗
-
μ
+
γ
=
μ
+
γ
R
0
-
1 and the roots of equation(12)
λ
2
+
2
μ
+
ε
+
ω
λ
+
μ
2
+
ε
μ
+
μ
ω
=
0
.
Obviously, in accordance with the relationship between roots and coefficients of quadratic equation, there is no positive real part characteristic root of (12). Hence, when R
0
<
1, (11) has no positive real root, and worm-free equilibrium Q
0 is locally asymptotically stable. When R
0
>
1, (11) has a positive root; thus, worm-free equilibrium Q
0 is an unstable saddle-point.Furthermore, the following theorem holds.Theorem 2.
The worm-free equilibrium is globally asymptotically stable ifR
0
≤
1
.Proof.
From the first equation of system (5)(13)
S
˙
t
≤
μ
+
ε
N
-
μ
+
ε
+
ω
S
.Thus, S
t
≤
μ
+
ε
N
/
μ
+
ε
+
ω
+
S
0
-
(
μ
+
ε
)
/
μ
+
ε
+
ω
e
x
p
-
μ
+
ε
+
ω
t. When t
→
∞, we obtain(14)
S
t
≤
μ
+
ε
μ
+
ε
+
ω
N
.Consider a Lyapunov function(15)
L
t
=
E
t
,
L
˙
t
=
ξ
S
t
I
t
-
μ
+
γ
I
t
≤
ξ
μ
+
ω
μ
+
ω
+
ε
N
-
μ
+
γ
I
t
=
1
μ
+
γ
R
0
-
1
I
t
≤
0
.So, we prove that worm-free equilibrium Q
0 is globally asymptotically stable.


## 3.2. Endemic Equilibrium and Its Stability

Now, we investigate the local stability of endemic equilibriumQ
∗. The characteristic equation of system (5) at endemic equilibrium Q
∗ is (16)
det
⁡
-
ξ
I
1
∗
-
μ
+
ω
-
λ
-
ξ
S
1
∗
ε
ξ
I
1
∗
ξ
S
1
∗
-
μ
+
γ
-
λ
0
ω
γ
μ
+
ε
-
λ
=
0
,which is equivalent to(17)
λ
3
+
p
0
λ
2
+
p
1
λ
+
p
2
=
0
,where p
0
=
3
μ
+
ε
+
ω
+
(
μ
+
ε
/
μ
+
ε
+
γ
)
N, p
1
=
μ
+
ε
2
μ
+
γ
+
(
μ
+
ε
/
μ
+
ε
+
γ
)
N
+
μ
ω
+
ξ
2
S
1
∗
I
1
∗, and p
2
=
μ
ξ
μ
+
ε
I
1
∗
-
S
1
∗
+
μ
μ
+
γ
μ
+
ε
+
ω
+
μ
ξ
γ
I
1
∗
-
ω
S
1
∗.Obviously,p
0
>
0, p
1
>
0, p
2
>
0, and p
0
p
1
-
p
2
>
0. According to the theorem of Routh-Hurwitz [26, 27], it follows that the roots of (17) have negative real parts. Therefore, the endemic equilibrium is locally asymptotically stable.From the above discussion, we can summarize the following conclusion.Lemma 3.
IfR
0
>
1, then endemic equilibrium is locally asymptotically stable.Note that the number of nodes in WSNs is relatively stable; that is, at timet, the number of nodes S
t, I
t, and R
t in states S
,
I, and R, respectively, satisfies(18)
S
t
+
I
t
+
R
t
=
N
.Hence, the dynamics of system (5) is equivalent to the following system:(19)
d
S
t
d
t
=
μ
+
ω
N
-
ξ
S
t
I
t
-
μ
+
ω
+
ε
S
t
-
ε
I
t
,
d
I
t
d
t
=
ξ
S
t
I
t
-
μ
+
γ
I
t
.Clearly, when R
0
>
1, system (19) has a unique positive equilibrium Q
∗
S
1
∗
,
I
1
∗, where (20)
S
1
∗
,
I
1
∗
=
μ
+
γ
ξ
,
ξ
μ
+
ε
N
-
μ
+
γ
μ
+
ε
+
ω
ξ
μ
+
ε
+
γ
.Now, we state and prove a result on the global stability of the endemic equilibrium S
1
∗
,
I
1
∗ of system (19).Theorem 4.
WhenR
0
>
1, the endemic equilibrium S
1
∗
,
I
1
∗ of system (19) is globally asymptotically stable.Proof.
Consider the following Lyapunov function [28]:(21)
L
t
=
∫
S
1
∗
S
x
-
S
1
∗
x
d
x
+
∫
I
1
∗
I
x
-
I
1
∗
x
d
x
.The time derivative of L
t along the solution of system (19) is given by(22)
L
˙
t
=
S
-
S
1
∗
S
S
′
+
I
-
I
1
∗
I
I
′
=
1
-
S
1
∗
S
μ
+
ε
N
-
ξ
S
t
I
t
-
μ
+
ε
+
ω
S
t
-
ε
I
t
+
1
-
I
1
∗
I
ξ
S
t
I
t
-
μ
+
γ
I
t
≤
1
-
S
1
∗
S
μ
+
ε
N
-
ξ
S
t
I
t
-
μ
+
ε
+
ω
S
t
+
1
-
I
1
∗
I
ξ
S
t
I
t
-
μ
+
γ
I
t
=
-
μ
+
ε
N
S
S
1
∗
S
1
∗
S
-
1
2
≤
0
.The proof is completed.Remark 5.
By Theorem4, we obtain that R
0
>
1; then endemic equilibrium Q
∗ of system (5) is globally asymptotically stable.


## 4. Worm Propagation Threshold Analysis and Numerical Simulations

We have proved that the basic reproductive numberR
0 equaling zero is the threshold whether worms are eliminated. When R
0
≤
1, worms in WSNs can be eliminated, and system (5) will stabilize at worm-free equilibrium. When R
0
>
1, worms in WSNs will exist consistently, and system (5) will stabilize at the endemic equilibrium. For verifying the correctness of theoretical analysis, we perform worm propagation threshold analyses and numerical simulations from the following two sides.


### 4.1. Communication Radius of Nodesr

LetR
0
=
ξ
μ
+
ε
/
μ
+
γ
μ
+
ε
+
ω
N
=
1; we can get the threshold of worm propagation about communication radius of nodes r
c
=
L
μ
+
γ
μ
+
ε
+
ω
/
π
β
μ
+
ε
N; that is, when r
≤
r
c, R
0
≤
1, according to Theorem 2, worms in WSNs can be eliminated, and system (5) will stabilize at the worm-free equilibrium; when r
>
r
c, R
0
>
1, according to Remark 5, worms in WSNs will exist consistently, and system (5) will stabilize at the endemic equilibrium.We choose a set of simulation parameters as follows:N
=
1000, L
=
10, μ
=
0.001, β
=
0.0003, ω
=
0.001, ε
=
0.0003, and γ
=
0.002. By calculation, we have r
c
=
0.7506. Initial values of susceptible, infected, and recovered nodes in WSNs are S
0
=
990, I
0
=
10, and R
0
=
0. When r takes different values, simulation results are depicted in Figures 3(a)–3(d).Figure 3
Simulation results of different communication radius of nodes: (a)r
=
0.1; (b) r
=
0.3; (c) r
=
0.8; and (d) r
=
1.
(a)
(b)
(c)
(d)Whenr
=
0.1
<
r
c and r
=
0.3
<
r
c, Figures 3(a) and 3(b) show that system (5) stabilizes at worm-free equilibrium. The simulation results are consistent with Theorem 2.Whenr
=
0.8
>
r
c and r
=
1
>
r
c, Figures 3(c) and 3(d) show that system (5) stabilizes at the endemic equilibrium. Simulation results are consistent with Remark 5.


### 4.2. Nodes Distributed Densityρ

LetR
0
=
ξ
μ
+
ε
/
μ
+
γ
μ
+
ε
+
ω
N
=
π
r
2
β
μ
+
ε
/
μ
+
γ
μ
+
ε
+
ω
ρ
=
1; we can get the threshold of worm propagation about node distributed density ρ
c
=
μ
+
γ
μ
+
ε
+
ω
/
π
r
2
β
μ
+
ε. When ρ
<
ρ
c, R
0
≤
1, system (5) has only a worm-free equilibrium and is globally asymptotically stable; when ρ
>
ρ
c, R
0
>
1, system (5) has an endemic equilibrium besides the worm-free equilibrium, and endemic equilibrium is globally asymptotically stable.We choose a set of simulation parameters as follows:N
=
1000, r
=
1, μ
=
0.001, β
=
0.0003, ω
=
0.001, ε
=
0.0003, and γ
=
0.002. By calculation, we get ρ
c
=
5.6345. Initial values of system (5) are set as S
0
=
990, I
0
=
10, and R
0
=
0. When L
=
20, 25, 12, and 10, we can get ρ
=
N
/
L
2
=
2.5, 1.6, 6.944, and 10. Simulation results are depicted in Figures 4(a)–4(d).Figure 4
Simulation results of different nodes density: (a)ρ
=
2.5; (b) ρ
=
1.6; (c) ρ
=
6.944; and (d) ρ
=
10.
(a)
(b)
(c)
(d)Whenρ
=
2.5
<
ρ
c and ρ
=
1.6
<
ρ
c, Figures 4(a) and 4(b) indicate that system (5) stabilizes at the worm-free equilibrium, and worm propagation is controlled finally. Simulation results are consistent with theoretical analysis.Whenρ
=
6.944
>
ρ
c and ρ
=
10
>
ρ
c, Figures 4(c) and 4(d) show that the trajectories converge to the endemic equilibrium. The conclusions agree with theoretical analysis.


## 4.1. Communication Radius of Nodesr

LetR
0
=
ξ
μ
+
ε
/
μ
+
γ
μ
+
ε
+
ω
N
=
1; we can get the threshold of worm propagation about communication radius of nodes r
c
=
L
μ
+
γ
μ
+
ε
+
ω
/
π
β
μ
+
ε
N; that is, when r
≤
r
c, R
0
≤
1, according to Theorem 2, worms in WSNs can be eliminated, and system (5) will stabilize at the worm-free equilibrium; when r
>
r
c, R
0
>
1, according to Remark 5, worms in WSNs will exist consistently, and system (5) will stabilize at the endemic equilibrium.We choose a set of simulation parameters as follows:N
=
1000, L
=
10, μ
=
0.001, β
=
0.0003, ω
=
0.001, ε
=
0.0003, and γ
=
0.002. By calculation, we have r
c
=
0.7506. Initial values of susceptible, infected, and recovered nodes in WSNs are S
0
=
990, I
0
=
10, and R
0
=
0. When r takes different values, simulation results are depicted in Figures 3(a)–3(d).Figure 3
Simulation results of different communication radius of nodes: (a)r
=
0.1; (b) r
=
0.3; (c) r
=
0.8; and (d) r
=
1.
(a)
(b)
(c)
(d)Whenr
=
0.1
<
r
c and r
=
0.3
<
r
c, Figures 3(a) and 3(b) show that system (5) stabilizes at worm-free equilibrium. The simulation results are consistent with Theorem 2.Whenr
=
0.8
>
r
c and r
=
1
>
r
c, Figures 3(c) and 3(d) show that system (5) stabilizes at the endemic equilibrium. Simulation results are consistent with Remark 5.


## 4.2. Nodes Distributed Densityρ

LetR
0
=
ξ
μ
+
ε
/
μ
+
γ
μ
+
ε
+
ω
N
=
π
r
2
β
μ
+
ε
/
μ
+
γ
μ
+
ε
+
ω
ρ
=
1; we can get the threshold of worm propagation about node distributed density ρ
c
=
μ
+
γ
μ
+
ε
+
ω
/
π
r
2
β
μ
+
ε. When ρ
<
ρ
c, R
0
≤
1, system (5) has only a worm-free equilibrium and is globally asymptotically stable; when ρ
>
ρ
c, R
0
>
1, system (5) has an endemic equilibrium besides the worm-free equilibrium, and endemic equilibrium is globally asymptotically stable.We choose a set of simulation parameters as follows:N
=
1000, r
=
1, μ
=
0.001, β
=
0.0003, ω
=
0.001, ε
=
0.0003, and γ
=
0.002. By calculation, we get ρ
c
=
5.6345. Initial values of system (5) are set as S
0
=
990, I
0
=
10, and R
0
=
0. When L
=
20, 25, 12, and 10, we can get ρ
=
N
/
L
2
=
2.5, 1.6, 6.944, and 10. Simulation results are depicted in Figures 4(a)–4(d).Figure 4
Simulation results of different nodes density: (a)ρ
=
2.5; (b) ρ
=
1.6; (c) ρ
=
6.944; and (d) ρ
=
10.
(a)
(b)
(c)
(d)Whenρ
=
2.5
<
ρ
c and ρ
=
1.6
<
ρ
c, Figures 4(a) and 4(b) indicate that system (5) stabilizes at the worm-free equilibrium, and worm propagation is controlled finally. Simulation results are consistent with theoretical analysis.Whenρ
=
6.944
>
ρ
c and ρ
=
10
>
ρ
c, Figures 4(c) and 4(d) show that the trajectories converge to the endemic equilibrium. The conclusions agree with theoretical analysis.


## 5. Conclusions

In this paper, we have proposed an improved SIRS model for analyzing dynamics of worm propagation in WSNs. This model can describe the process of worm propagation with the energy consumption and different distributed density of nodes. Based on this model, a control parameterR
0 that completely determines the global dynamics of worm propagation has been obtained by the explicit mathematical analyses. From Theorems 4 and 2, we learn out that worm will be controlled in WSNs when R
0
<
1, and they will be prevalent otherwise. Finally, based on R
0, we discuss the threshold of worm propagation about communication radius and distributed densities of nodes in WSNs. Numerical simulations verify the correctness of theoretical analysis. Research results show that decreasing the value of communication radius or reducing distributed density of nodes is an effective method to prevent worms spread in WSNs. Research of this paper provides the theoretical basis for predicting and controlling worm propagation in WSNs. It is worth pointing out that we do not consider physical effects like “collisions” and heterogeneous distribution of nodes on infection rate when modeling, which is our focus in the future.



---

*Source: 129598-2015-09-02.xml*