# The Positive Role of Multiplicative Noise in Complete Synchronization of Unidirectionally Coupled Ring with Three Nodes

Authors: Yuzhu Xiao; Sufang Tang; Zhongkui Sun
Journal: Journal of Applied Mathematics
 (2014)

---

## Abstract

The role of multiplicative noise in the synchronization of unidirectionally coupled ring with three nodes is studied. Based on the theory of stochastic differential equations, we demonstrate that noise plays a positive role in complete synchronization. In numerical simulations, the Lorenz system, Rössler like system, and Hindmarsh-Rose neuron model are employed to demonstrate the correctness of our theoretical result.


---

## Body

## 1. Introduction

Synchronization is an omnipresent phenomenon in nonlinear systems wherein two or more systems coupled by a suitable configuration or driven by external forcing adjust one property of their motion to a common behavior [1, 2]. With the development of chaos synchronization, different kinds of synchronization have been found, such as complete synchronization [2–5], phase synchronization [6], lag synchronization [7, 8], and generalized synchronization [9, 10]. The synchronization of nonlinear systems has many potential applications in different areas, such as secure communication [11, 12], biological systems [13, 14], chemical reactions [15], and physical systems [16]. Due to the potential applications of synchronization, some synchronization methods have been proposed, such as active control [17], adaptive control [18], sliding control [19], and integral observer [20].Since noise is ubiquitous in real-world systems, the synchronization of coupled systems is unavoidably affected by different kinds of noise. Therefore, the role of noise in synchronization is very important and has been extensively studied. Intuitively, noise is random disturbance, which should play a destroyed role in the synchronization of coupled systems. However, some recent studies have shown that synchronization can be induced by noise. Firstly, noise-induced synchronization was found in periodic systems [21, 22]. With the development of chaos synchronization, the phenomenon of noise-induced synchronization was extended to chaotic map [23]. The result in [23] has spurred a long-standing dispute since it was found. Firstly, Pikovsky pointed out that the result in [23] is not right, because the largest Lyapunov exponent is positive [24]. Some other authors argued that synchronization in [23] is caused by biased noise. One explanation is that biased noise can move the system into “convergence region,” so unbiased noise can not induce synchronization [25–27]. However, some recent work showed that some chaotic systems can be synchronized by unbiased additive noise [28, 29]. With the extensive study of noise’s effect on synchronization, the phenomenon of noise-induced synchronization has been extended to generalized synchronization [30, 31], phase synchronization [32, 33], and lag synchronization [34]. However, the above studies mostly concentrate on the case of additive noise, and most of results are based on numerical simulations and experiments.Recently, the study of noise’s effect on synchronization is extended to multiplicative noise. In [35], the authors studied two coupled chaotic systems perturbed by multiplicative noise. Based on the invariance principle of stochastic differential equations, the sufficient condition ensuring the complete synchronization was established. Sooner, above method was generalized to coupled chaotic systems with multiple time delays [36] and linear generalized synchronization [37]. In [38], the sufficient conditions of complete synchronization of linear stochastic coupled network under adaptive control were established. However, the role of noise in synchronization was not discussed. In [39], the positive role of multiplicative noise in synchronization was revealed in two unidirectionally coupled chaotic systems, and authors proposed a good analytical method to analyze this phenomenon. In [40], two bidirectionally coupled chaotic systems were studied by above method. In this paper, the method in [39] will be extended to three coupled chaotic systems with unidirectional ring structure. The coupled ring has been used widely in modeling physiological, biochemical, and biological phenomena [13, 14]. Though many real-world systems should be described by the network with many nodes and complex structure, the ring with three nodes is the network motifs of some complex networks. The investigation of the synchronization of ring with three nodes is helpful to understand the collective behavior in complex networks [41].


## 2. Theoretical Analysis

The unidirectionally coupled ring with three nodes can be described as follows:(1)x˙1=f(x1)+c(x2-x1),x˙2=f(x2)+c(x3-x2),x˙3=f(x3)+c(x1-x3).
Here xi=(xi1,xi2,…,xim)T∈Rm is the state vector of node i (i=1,2,3), and f=(f1,f2,…,fm)T is the vector function which describes the dynamic of uncoupled system. c is a positive constant, which is usually called coupling strength. Consider that the real-world systems are unavoidably disturbed by internal noise which may lead to the loss of information in transmission process [38]; then the coupled ring model (1) can be improved as follows:
(2)x˙1=f(x1)+(c+dξ(t))(x2-x1),x˙2=f(x2)+(c+dξ(t))(x3-x2),x˙3=f(x3)+(c+dξ(t))(x1-x3),
where d is a constant which describes the noise intensity, and ξ(t) is one-dimensional Gaussian white noise with statistical properties 〈ξ(t)〉=0 and 〈ξ(t),ξ(t′)〉=δ(t-t′). In fact, the assumption of Gaussian white noise is often impractical, but this assumption is a good approximation of many real-world situations and generates mathematically tractable models. The noise term in (2) is of multiplicative form, so (2) can be interpreted in the sense of Itô or Stratonovich. The Itô sense is employed here; though this interpretation seems to be mathematical limiting procedure, this is reasonable when the rapidity of fluctuation of environment is far less than the macroscopic time scale of the system, such as the external fluctuation in many biological systems [39, 42]. In most of real physical systems, the Stratonovich sense is usually employed. However, the recent study shows that the multiplicative noise in electric circuit, which usually obeys the Stratonovich convention, crosses over to obey the Itô convention when the certain parameters of the systems are changed [43]. In a word, the Itô interpretation of model is not totally idealizing but is realizable in many real-world systems.To begin the theoretical analysis, the coupled ring (2) and the vector function f=(f1,f2,…,fm)T need to satisfy the following assumption.Assumption.
Assume that the coupled ring (2) is bounded, and for any bounded set Ω∈Rm, there exists a positive constant l such that x, y∈Ω satisfies
(3)(x-y)T[f(x)-f(y)]≤l(x-y)T(x-y),
where l is a positive constant. Notice that this condition is very loose; for example, the condition (3) is satisfied as long as ∂fi/∂xi are bounded on Ω. Therefore, most of well-known chaotic systems satisfy condition (3).Actually, the coupled ring (2) is said to be synchronized if x1=x2=x3→s(t) as t→∞.s(t)∈Rm is called synchronized manifold which is usually a solution of uncoupled system; that is, s˙=f(s). Inspired by [44], the manifold s-(t)=(x1+x2+x3)/3 is introduced to study the stability of s(t). This manifold has some new properties which is helpful to theoretical analysis. It is easy to see that
(4)s-˙(t)=13∑i=13f(xi)≜g(s-).Defineei=xi-s-(t) to be the error state between the state of ith node and manifold s-(t). Obviously, the error states are dominated by the following equation:
(5)e˙1=f(x1)-g(s-(t))+(c+dξ(t))(-e1+e2),e˙2=f(x2)-g(s-(t))+(c+dξ(t))(-e2+e3),e˙3=f(x3)-g(s-(t))+(c+dξ(t))(-e3+e1).
Here, one should notice that error states ei (i=1,2,3) satisfy the following equality:
(6)e1+e2+e3=0.
This property is very useful in the following theoretical analysis.For convenience of expression, (5) is rewritten as the following matrix form:
(7)E˙(t)=F(x)+C(E)+H(E)W˙(t).
Here
(8)E(t)=(e1e2e3),F(x)=(f(x1)-g(s(t))f(x2)-g(s(t))f(x2)-g(s(t))),C(E)=c(-e1+e2-e2+e3-e3+e1),H(E)=d(-e1+e2-e2+e3-e3+e1),W˙(t)=ξ(t).According to the stability theory of stochastic differential equations [39, 45], the error equation (7) possesses a unique global solution with any initial state, without loss of generality, denoted by E(t;t0,E0). It is easy to see that E(t;0,0)≡0 is a trivial solution of the error system (7). The complete synchronization of dynamical network (2) is achieved, if the zero solution is asymptotically stable; that is, lim⁡t→∞⁡∥E(t)∥=0 with probability one.In the following, the role of noise in the synchronization of the coupled ring (2) will be analyzed by the theory of stochastic differential equations. To this end, the following function is defined:
(9)V(E)=12log⁡(ET(t)E(t)).
By [39], E(t;t0,E0) never reaches zero with probability one if E0≠0. This ensures the correctness of adopting the logarithm in (9). By applying Itô formula to (9) along with system (7) [39, 45], one has
(10)V(E(t))=V(E(t0))+∫t0t{12DV(E(s))[F(x)+C(E)]+12trace[HT(E(s))D2(V(E(s)))H(E(s))]}ds+M(t),
where
(11)DV(E)=ET(t)ET(t)E(t),D2V(E)=IET(t)E(t)-2E(t)ET(t)(ET(t)E(t))2,
and the continuous martingale is
(12)M(t)=∫t0tET(s)H(E(t))ET(s)E(s)dW(s).
The quadratic variation of (12) is
(13)[M(t),M(t)]=∫t0t[ET(s)H(E(t))]2[ET(s)E(s)]2ds=d2∫t0t[-∑i=13eiTei+e1Te2+e2Te3+e3Te1]2(∑i=13eiTei)2ds=d2∫t0t[-(3/2)∑i=13eiTei]2(∑i=13eiTei)2ds=94d2(t-t0).
Notice that the condition (6) and the equality
(14)eiTej=12(ei+ej)T(ei+ej)-12(eiTei+ejTej)
are used in above proof.By the strong law of large numbers [39, 45], one has that
(15)lim⁡t→∞M(t)t=0
holds almost sure.SubstitutingE(t), DV(E), D2V(E), F(x), C(E), and H(E) into (10) and applying equality (6), one has
(16)V(E(t))=V(E(t0))+∫t0t∑i=13eiT[f(xi)-f(s(t))]∑i=13eiTeids+c∫t0t-(3/2)∑i=13eiTei∑i=13eiTeids+12d2∫t0t{∑i=133eiTei∑i=13eiTei-2[-(3/2)∑i=13eiTei]2(∑i=13eiTei)2}ds+M(t)=V(E(t0))+∫t0t∑i=13eiT[f(xi)-f(s(t))]∑i=13eiTeids-32c(t-t0)-34d2(t-t0)+M(t).If the coupled ring (2) is bounded and vector function f satisfies the condition (3), it yields
(17)V(E(t))≤V(E(t0))+l(t-t0)-32c(t-t0)-34d2(t-t0)+M(t).
Combining (15), one obtains
(18)lim⁡t→∞sup⁡V(E(t))t≤l-32c-34d2.
Obviously, if l-(3/2)c-(3/4)d2<0, the logarithm function V(E(t)) approaches to -∞ as t→+∞; that is, E(t) approaches to zero as t→+∞. In other words, the trivial solution E(t;0,0)≡0 of error system (7) is asymptotically stable with probability one when l<(3/2)c+(3/4)d2. Moreover, the last term about the noise intensity in inequality (18) is negative, so the noise term has positive effect on the synchronization of the coupled ring (2). Particularly, when coupling strength c=0, the coupled ring (2) can synchronize with l<(3/4)d2. This means synchronization may be achieved when all nodes are just coupled by multiplicative noise.Remark 1.
One should notice that the above inequality (18) is acquired with bounded assumption of coupled ring (2) and the condition (3). For large noise intensity, the coupled ring (2) may become unbounded. Though it is difficult for us to acquire the bounded condition of coupled ring (2), coupled ring (2) is usually bounded when noise intensity d is small.


## 3. Numerical Examples

To demonstrate that noise term in coupled ring (2) plays a positive role in synchronization, the Lorenz systems [46], Hindmarsh-Rose neuron system [47], and Rössler like system [48] are employed to describe the dynamic of uncoupled system. In numerical simulations, the weak order 2 Runge-Kutta scheme with time step Δt=10-2 is adopted to solve the stochastic differential equations [49].To simplify the representation, a scale is defined to measure the synchronization of the coupled ring (2), which is
(19)Δ(d)=〈13∑k=13∑j=1m∫T1T2|xkj(t,ω,d)-S-j(t,ω,d)|dtT2-T1〉,
where ω∈Ω, Ω denotes the set of all elementary events, 〈·〉 represents sample average, and m is the dimension of node system. According the definition, Δ(d) can be used to exhibit effect of synchronization with respect to noise intensity.Example 1 (Lorenz system).
The famous Lorenz system can be described as follows:(20)x˙=-σ(x-y),y˙=-xz+rx-y,z˙=xy-bz.
Here x, y, and z make up the system state, and σ, r, b are positive system parameters. In this example, let σ=10, r=24.1, b=8/3. Under this group of parameters, system (20) has chaotic solution.In this example, the coupling strengthc is set to be zero. 103 sample paths are computed with initial state xi=(2.54+ui1,2.52+ui2,17.27+ui3), where uij are randomly chosen from [0,0.3]. In Figure 1, we plot the evolution of Δ with respect to d. Form this figure, one can see that the scale Δ decreases with d increasing, and three Lorenz systems just coupled by internal noise can achieve synchronization when the noise strength is large enough.Figure 1
The evolution ofΔ with respect to d in three coupled Lorenz systems, when c=0, T0=4×103, and T1=5×103.Example 2 (Rössler like system).
The Rössler like system is dominated by(21)x˙=Dx+g(x),
where x=(x1,x2,x3)T∈R3 is the state vector, and
(22)D=α(-Γ-β-λ1γ000-η),f(x)={0,x<θ,ζ(x-θ),x≥θ,g(x)=(00αηf(x1))T.
This system is a modified version of original Rössler system and has been used to study the synchronization of coupled chaotic systems in experiment. Under parameters α=1.0, β=1.5, γ=0.2, η=1.5, Γ=0.075, λ=0.75, ζ=18.43, and θ=0.5, the system (21) has a chaotic attractor.In this example, we set coupling strengthc=0.045, and 103 sample paths are computed with initial states xi=(-0.33+ui1,-0.03+ui2,0.01+ui3), where uij are randomly chosen from [0,0.3]. In Figure 2, we plot the evolution of Δ with respect to d. From this figure, one can see that Δ is order of 10-2 when noise intensity d=0. With increase of d, the order of Δ decrease form 10-2 to 10-11. This demonstrates that multiplicative noise really plays a positive role in the synchronization of three coupled Rössler like systems.Figure 2
The evolution ofΔ with respect to d in three coupled Rössler like systems, when c=0.045, T0=4×103, and T1=5×103.Example 3 (Hindmarsh-Rose (HR) neuron model).
The HR neuron model can be described as follows:(23)x˙=y-ax3+bx2-z+I,y˙=c-dx2-y,z˙=r(s(x-χ)-z),
where x is the membrane potential and y, z describe the transport of ions across the membrane through the ion channels. In this example, let the parameters a=1, b=3, c=1, d=5, s=4, r=0.015, χ=-1.6, and external forcing I=2.95. Under above parameters, the system (23) exhibits chaotic spike.In this example, we set coupling strengthc=0.002. 103 sample paths are computed with initial states xi=(1.05+ui1,-4.71+ui2,3.26+ui3) where uij are randomly chosen from [0,0.1]. The evolution of Δ with respect to d is plotted in Figure 3. Form this figure, one can see that Δ is order of 100 when noise intensity d=0. With increase of d, the order of Δ decreases to 10-11. This demonstrates that noise really plays a positive role in the synchronization of coupled HR neuron systems.Figure 3
The evolution ofΔ with respect to d in three coupled Lorenz systems, when c=0.002, T0=4×103, and T1=5×103.


## 4. Conclusions and Discussions

To summarize, we study the role of multiplicative noise in the synchronization of unidirectionally coupled ring with three nodes. Based on the stability theory of stochastic differential equations, we demonstrate that multiplicative noise really has positive effect on synchronization. To demonstrate the correctness of the theoretical result, the Lorenz system, HR neuron system, and Rössler like system are employed as numerical examples. In this paper, we just study the coupled ring with three nodes and common noise. To be more practical, different nodes should be subjected to different noise. Moreover, many real systems should be described by networks with many nodes and complex structures. The role of multidimensional noise in the synchronization of complex dynamical network is one of our further works. For unidirectionally coupled ring with more than three nodes, we can not demonstrate that noise plays a positive role in synchronization by present method. We performed some numerical simulations for ring with many nodes. We found that noise plays a positive role in the synchronization of some samples, but some samples become unbounded. Is this unboundedness raised by artifact of finite precision of numerical simulation or determined by the structure of node system and ring? The theoretical result of this paper is acquired under the assumption of boundedness of ring. Though the bounded condition of the ring is difficult to be acquired by theoretical method, the bounded condition of network (2) is a problem which is worthwhile to be investigated. Further, the circuit experiment [50] is also best way to demonstrate the theoretical results, but the design of circuit is beyond the extent of our knowledge.



---

*Source: 741961-2014-07-16.xml*