# Finite-Time Projective Lag Synchronization and Identification between Multiple Weights Markovian Jumping Complex Networks with Stochastic Perturbations

Authors: Qian Xie; Changhui Mu; Tong Wang; Gang Wu; Rong Jia
Journal: Complexity
 (2020)

---

## Abstract

Two nonidentical dimension Markovian jumping complex networks with stochastic perturbations are taken as objects. The network models under two conditions including single weight and double weights are established, respectively, to study the problem of synchronization and identification. A finite-time projection lag synchronization method is proposed and the unknown parameters of the network are identified. First of all, based on Itô’s formula and the stability theory of finite-time, a credible finite-time adaptive controller is presented to guarantee the synchronization of two nonidentical dimension Markovian jumping complex networks with stochastic perturbations under both conditions. Meanwhile, in order to identify the uncertain parameters of the network with stochastic perturbations accurately, some corresponding sufficient conditions are given. Finally, numerical simulations under two working conditions are given to demonstrate the effectiveness and feasibility of the main theory result.


---

## Body

## 1. Introduction

Complex network, a form of system structure and function, is an essential abstraction of the interaction of the complex system and has attracted substantial number of interest from researchers in varied realm [1, 2]. The World Wide Web, power network, and neural network as well as social network all belong to the complex network of our life. Synchronization is a momentous nonlinear phenomenon of nature. As a most valuable topic in many dynamic behaviors of the complex network, it has attracted a growing number of concerns in the investigation [3, 4]. If the dynamic behavior of the coupling node in the network evolves over time and eventually reaches the same state, then it will be called synchronization. At present, a series of related research results have been established, including complete synchronization, exponential synchronization, asymptotic synchronization, and generalized synchronization. In [5], a method of complete synchronization of complex networks is proposed by constructing adaptive control technology. A sufficient condition for exponential synchronization of a complex network with time-varying delays is given in [6]. Among these research studies, the dimensions of networks are identical. However, in some practical situations this hypothesis is inappropriate and there are still a large number of networks formed with nonidentical dimensions. Projective synchronization, an important synchronization phenomenon, has received adequate attention which can be used as an appropriate synchronization scheme for the nonidentical dimension complex networks.Moreover, the abovementioned network synchronization assumes that the synchronization time tends to be infinite [7, 8]. However, in the engineering field, it is often needed to achieve synchronization as quickly as possible, which means finite-time synchronization. In order to accomplish faster synchronization in complex dynamical networks, an effective approach is adopting finite-time control methods. In addition, as an important factor, stochastic perturbations will ineluctably affect the dynamical of the complex network and even destroy the stability of the system in some kinds of practical situations [9–11]. Therefore, in complex dynamical networks, both the stochastic perturbations and how to achieve finite-time synchronization cannot be ignored. Withal, at reality operation conditions, the topology of the network would be randomly switched over time [12, 13]. To this end, a Markovian process is introduced to describe. Because of the powerful mathematical modeling ability of the Markovian process in many situations, this kind of network has been developed and studied promptly.At the meantime, in many complex network synchronization research studies, the system parameters are known. In [14], the asymptotic synchronization of complex networks with nonlinear nodes of different dimensions is studied, but the parameters of the node system are known. However, this situation is not sufficient in the process of practical. The problem of Markovian jumping complex network synchronization with uncertain system parameters needs to be considered.The research studies on complex networks at home and abroad mainly focus on the synchronization analysis of networks with single weight [15–18]. However, there are few studies on complex networks with multiple weights and seldom involve the synchronization and parameter identification of Markovian jumping complex networks with multiple weights. It is widespread in our life that there are many complex networks with multiple weights which are composed of many subnetworks of different nature, such as communication networks, traffic networks, and interpersonal networks as well as power networks [19, 20]. Taking the power network as an example, considering each city as a node and the power transmission line between two cities as an edge, there may be a variety of voltage levels of the power transmission network between the two cities, such as high-voltage network and low-voltage network (as shown in Figure 1). Such a power network constitutes a complex network with multiple weights.Figure 1
Schematic diagram of the power network.Being different from the single weight, its network topology and node dynamics will be more complex. When dealing with multiweighted networks, the method of merging the edges of different properties into one edge for processing is the simplest and most common. However, some edges of the network in real life cannot be simply merged and processed, which cannot truly reflect the characteristics of the network [21].Inspired by relevant research, this paper discusses the identification and synchronization problem of single-weight [22] and double-weights Markovian jumping complex networks with stochastic perturbations in finite time. This paper will divide double-weights network into two single-weight subnetworks. Then, it is superimposed to establish the network model by the thought of network decomposition. As far as the author knows, there are rarely studies on the identification and synchronization of nonidentical dimension Markovian jumping complex networks with uncertain parameters in finite time; hence, it is worthwhile to conduct research from the practice and theory.


## 2. Preliminaries

Firstly, some significantly mathematical notations are introduced as follows, which will be applied throughout this paper. Let matrixAT (or xT) means the transpose of the A (or x). Denote x as the 2-norm. Indicate In ∈ Rn×n as the n-dimensional identity matrix. Let ⊗ represent the Kronecker product of the matrix. λmax⋅ is the largest eigenvalue of a matrix.Definition 1.
(see [4]). Define f:Rn×R+×S⟶Rn and g:Rn×R+×S⟶Rn×m are differentiable functions satisfying the following n-dimensional functional stochastic differential equation:(1)dxt=fxt,t,rtdt+gxt,t,rtdωt.
According to Itô’s formula, ifVxt,t,r∈C2,1Rn×R+×S;R+, the operator ℒV is defined as follows:(2)ℒVxt,t,p=Vtxt,t,p+Vxxt,t,pfxt,t,p+12trancegTxtt,pVxxgxt,t,p+∑q=1mπpqVxt,t,q,where Vtxt,t,p=∂Vxt,t,p/∂t, Vxxt,t,p=∂Vxt,t,p/∂x1,…,∂Vxt,t,p/∂xn, and Vxxxt,t,p=∂2Vxt,t,p/∂xi∂xjn×n.Assumption 1.
The uniform Lipschitz condition is satisfied by the noise intensity functionσi(t, ei(t)) and the constant ρi ≥ 0 is existed to make the following formula hold:(3)tranceσiTt,eit,rtσit,eit,rt≤ρirteiTteit.Lemma 1.
(see [23]). As for any vectors x, y ∈ Rn, the following matrix inequality will be hold, where the P ∈ Rn×n is a positive definite matrix:(4)2xTy≤xTPx+yTP−1y.Lemma 2.
To suppose thatN(t), a continuous, positive definite function, satisfies the following differential inequality N(t) ≤ −αNη(t), ∀t ≥ t0, and N(t0) ≥ 0, where constants α > 0 and 0 < η < 1. As for arbitrary t0 and N(t), the following differential inequality will be satisfied:(5)N1−ηt≤N1−ηt0−α1−ηt−t0,t0≤t≤t1,N(t) ≡ 0, ∀t ≥ t1. With t1, it is given as follows:(6)t1=t0+N1−ηt0α1−η.Lemma 3.
(7)
a
1
γ
+
a
2
γ
+
⋅
⋅
⋅
+
a
n
γ
≥
a
1
2
+
a
2
2
+
⋅
⋅
⋅
+
a
n
2
γ
/
2
.
As well as(8)∑i=1nanγ≥∑i=1nan2r/2,where ∀a1, a2, …, an ∈ Rn are any vectors and γ is a real number satisfying 0 < γ < 2.


## 3. Finite-Time-Generalized Matrix Projective Synchronization with Single-Weight Networks

The single-weight Markovian jumping complex network with random perturbations is taken as the research object and its parameter identification and finite-time synchronization are studied. The corresponding synchronization criteria and updating rules are obtained as well as the reliability and validity of the method are illustrated by numerical simulation. This is a preparation for further research on synchronization of complex networks with double weights.


### 3.1. Network Models

The drive-coupled complex network is described as follows, and its number of dynamic nodes isN:(9)x˙t=fi1xitαi+fi2xit+∑j=1NaijrtΓ1xjt,where xit=xi1t,xi2t,…,xintT∈Rn are the state vectors of the ith node; fi1: Rn ⟶ Rn×l is a continuous matrix function and fi2: Rn ⟶ Rn is a continuous vector function, respectively; αi ∈ Rl are the unidentified node dynamic parameters; Γ1 ∈ Rn×n represents the inner coupling matrix; the coupling configuration matrix is expressed by Art=aijrtN×N∈Rn×n, which represents the topological structure and the coupling strength of the network at time t in mode r(t).{r(t, t ≥ 0)} is defined as a right-continuous Markovian process in probability space, describing the switching state between different parameters of time t. It takes values from the finite space S=1,2,…,m with generator Π=πpqm×mp,q∈S. The transition probability from pth mode at time t to the qth mode at time t + Δt will be defined as follows:(10)Prt+Δt=q|rt=p=πpqΔt+oΔt,ifq≠p,1+πpqΔt+oΔt,ifq=p,where Δt > 0:(11)limΔt⟶0oΔtΔt=0.As the transition rateπpq ≥ 0 which can be represented from mode p at time t to mode q at time t + Δt,(12)πpp=−∑q=1,q≠pmπpq.Therefore, the elements of matrixA(r(t)) will be defined as follows. If it exists as a link from node j to node i (j ≠ i), then aij ≠ 0; otherwise, aij = 0. The diagonal elements of matrix A(r(t)) can be determined by aii=−∑j=1,j≠iNaij.Describe the response-coupled network as follows:(13)y˙it=gi1yitβi+gi2yit+∑j=1NbijrtΓ2yit+uit+σit,eit,rtω˙t,where yit=yi1t,yi2t,…,yimtT∈Rm are the state vectors of node i; gi1: Rm ⟶ Rm×o is a continuous matrix function and gi2: Rm ⟶ Rm is a continuous vector function; the node dynamic parameter vectors which need to be identified are βi ∈ Ro; Γ2 ∈ Rm×m represents the inner coupling matrix; the coupling configuration matrix is defined by Brt=bijrtM×M∈Rm×m as well as in network (9); ui(t) is the finite-time nonlinear feedback controller which needs to be devised; σi(t, ei(t), and r(t)) are the noisy intensity function which mainly describe the impact of environmental fluctuations or imprecise design of coupling strength on network synchronization; and ei(t) is the error of system synchronization.The error of system synchronization can be defined as follows:(14)eit=yit−Mxit−τ,where M ∈ Rm×n is the scaling matrix (i.e., generalized matrix) consisting of constant; eit=ei1t,ei2t,…,eimtT; and τ is the time delay.Definition 2.
The driving network and response network can be realized synchronization through the designed finite-time nonlinear controllers. As for scaling matrixMij = (mij) ∈ Rm×n, if the condition can be satisfied that each row of the element cannot be 0 at a time, it will hold that(15)limt⟶∞Eeit2=limt⟶∞Eyit−Mxit−τ2=0.


### 3.2. Main Result

In this section, the synchronization and identification of unknown parameters for Markovian jumping complex networks with nonidentical dimensions and stochastic perturbations are investigated by utilizing the control theory of finite-time. The following theorem and remarks are the main results.First, from (9) and (13), the synchronization error dynamic system can be obtained as below:(16)e˙it=y˙it−Mx˙it−τ=gi1yitβi+gi2yit+uit+σit,eit,rtω˙t+∑j=1NbijrtΓ2eit+∑j=1NbijrtΓ2Mxit−τ−Mfi1xit−ταi+fi2xit−τ−M∑j=1NaijrtΓ1xjt−τ.Theorem 1.
With Assumption1 holding, the single-weight Markovian jumping complex network (9) and (13) can be realized synchronization in a finite-time. Meanwhile, with the following controllers and undated laws, the value of uncertain parameters can be simultaneously identified:(17)uit=−gi1yitβ^i−gi2yit+Mfi1xit−τα^i+Mfi2xit−τ−kireit+M∑j=1NaijrtΓ1−bijrtΓ2xjt−τ−ψdiagei1tγ,…,eimtγsigneit,α^˙it=−θifi1Txit−τMTteit+ψθiα^i−αiγsignα^i−αi,β^˙it=λigi1Txiteit−ψλiβ^i−βiγsignβ^i−βi,where ki(r) > 0, 0 < γ < 2, ψ > 0, θi and λi are any positive constant, and sign(·) represents symbolic function.Remark 1.
When Markovian jumping complex network (9) and (13) achieve projection lag synchronization in a finite time, the uncertain parameter vectors α^i and β^i of the system are identified.
If the following inequalities hold,(18)∑q=1mπpqDq−Cp≤0,B^r−Ξr+Θr⊗Im≤0,where D(q) is a positive definite matrix; C(p) is a symmetrical matrix; B^r=B˜rB˜Tr/2; B˜r=Br⊗Γ2; Ξr=diagk1r,k2r,…,kNr; Θr=diagρ1r,ρ2r,…,ρNr; and r ∈ S.
Therefore, the synchronization between single-weight Markovian jumping complex network (9) and (13) with stochastic perturbations can be achieved global asymptotic stability in a finite-time under this condition:(19)t1≤t0+2Vt01−γ/22ψ1−γ,where Vt0=1/2∑i=1NeiTt0eit0+1/2∑i=1N1/θiα˜iTα˜i+1/2∑i=1N1/λiβ˜iTβ˜i.Proof.
The Lyapunov function is defined as below:(20)Vt,et,p=12∑i=1NeiTteit+12∑i=1N1θiα˜iTα˜i+12∑i=1N1λiβ˜iTβ˜i.
By Definition1, we can obtain the differential operator ℒV as follows:(21)ℒVxt,t,p=∑i=1NeiTte˙it+∑i=1N1θiα^˙iTα˜i+∑i=1N1λiβ^˙iTβ˜i+∑q=1mπpq∑i=1NeiTtDqet+∑i=1NtranceσiTt,eit,pσit,eit,p.
Substitute (16) and feedback controllers (17) into (21), resulting in(22)ℒVxt,t,p=∑i=1NeiTtgi1yitβi+σit,eit,rtω˙t+∑j=1NbijrtΓ2eit+∑j=1NbijrtΓ2Mxit−τ+uit+gi2yit−Mfi1xit−ταi+fi2xit−τ−M∑j=1NaijrtΓ1xjt−τ−∑i=1NeiTtMfi1xit−τα˜i+∑i=1NeiTtgi1xitβ˜i−ψ∑i=1N1θiα˜iα^i−αiγsignα^i−αi−ψ∑i=1N1λiβ˜iβ^i−βiγsignβ^i−βi+∑i=1NtranceσiTt,eit,pσit,eit,p+∑q=1mπpq∑i=1NeiTtDqet,uit=−gi1yitβ^i−gi2yit+Mfi1xit−τα^i+Mfi2xit−τ−kireit+M∑j=1NaijrtΓ1−bijrtΓ2xjt−τ−ψdiagei1tγ,…,eimtγsigneit,(23)ℒVxt,t,p=∑i=1NeiTtgi1yitβi+σit,eit,rtω˙t+∑j=1NbijrtΓ2eit+∑j=1NbijrtΓ2Mxit−τ+gi2yit−Mfi1xit−ταi+fi2xit−τ−M∑j=1NaijrtΓ1xjt−τ−gi1yitβ^i−gi2yit+Mfi1xit−τα^i+Mfi2xit−τ−kireit+M∑j=1NaijrtΓ1−bijrtΓ2xjt−τ−ψdiagei1tγ,…,eimtγsigneit,−∑i=1NeiTtMfi1xit−τα˜i+∑i=1NeiTtgi1xitβ˜i−ψ∑i=1N1θiα˜iα^i−αiγsignα^i−αi−ψ∑i=1N1λiβ˜iβ^i−βiγsignβ^i−βi+∑i=1NtranceσiTt,eit,pσit,eit,p+∑q=1mπpq∑i=1NeiTtDqet,(24)ℒVxt,t,p=∑i=1NeiTtMfi1xit−τα^it−αi−gi1yitβ^it−βi+∑j=1NbijrtΓ2eit−kireit−ψdiagei1tγ,…,eimtγsigneit−∑i=1NeiTtMfi1xit−τα˜i+∑i=1NeiTtgi1yitβ˜i−ψ∑i=1N1θiα˜iα^i−αiγsignα^i−αi−ψ∑i=1N1λiβ˜iβ^i−βiγsignβ^i−βi+∑i=1NtranceσiTt,eit,pσit,eit,p+∑q=1mπpq∑i=1NeiTtDqet,(25)ℒVxt,t,p=∑i=1NeiTt∑j=1NbijrtΓ2eit−kireit−ψdiagei1tγ,…,eimtγsigneit+∑i=1NtranceσiTt,eit,pσit,eit,p−ψ∑i=1N1θiα˜iα^i−αiγsignα^i−αi−ψ∑i=1N1λiβ˜iβ^i−βiγsignβ^i−βi+∑q=1mπpq∑i=1NeiTtDqet.
According to Lemma3, we obtain that(26)∑i=1NeiTteitφsigneit=∑i=1NeiTtdiagei1tγ,…,eimtγsigneit=∑i=1N∑j=1Neijteijtγsigneijt=∑i=1N∑j=1Neijt1+γ≥∑i=1N∑j=1Neijt21+γ/2=∑i=1NeiTteit1+γ/2⟹∑i=1NeiTteitφsigneit≥∑i=1NeiTteit1+γ/2.
From Assumption1, the conclusion is drawn:(27)ℒVxt,t,p≤∑i=1NeiTtB^r−Ξr+Θr⊗Imeit+∑q=1mπpq∑i=1NeiTtDq−Cpet−ψ∑i=1NeiTteit1+γ/2−ψ∑i=1N1θiα˜iTα˜i1+γ/2−ψ∑i=1N1λiβ˜iTβ˜i1+γ/2≤−ψ∑i=1NeiTteit+∑i=1N1θiα˜iTα˜i+∑i=1N1λiβ˜iTβ˜i1+γ/2≤−2ψ12∑i=1NeiTteit+12∑i=1N1θiα˜iTα˜i+12∑i=1N1λiβ˜iTβ˜i1+γ/2=−2ψVt1+γ/2⟹ℒVxt,t,p≤−2ψVt1+γ/2.
From Lemma2, the system error of synchronization ei(t) will converge to a steady state with a finite-time, which can be estimated by(28)t1≤t0+Vt01−1+γ/22ψ1−1+γ/2=t0+2Vt01−γ/22ψ1−γ⟹t1≤t0+2Vt01−γ/22ψ1−γ,where Vt0=1/2∑i=1NeiTt0eit0+1/2∑i=1N1/θiα˜iTα˜i+1/2∑i=1N1/λiβ˜iTβ˜i.
Hence, networks (9) and (13) can be realized synchronization with finite time t1. Also, the uncertain system parameter vectors α^i and β^i are adapted to the true value of parameters. The proof is completed.Remark 2.
The proposed approach is applicable to the finite-time-generalized matrix projection lag synchronization of any two single-weight Markovian jumping complex networks with uncertain parameters, stochastic perturbations, and different initial conditions. In addition, the Markovian jumping complex network can be identical or nonidentical dimensions.Remark 3.
The synchronization and identification speed can be adjusted through choosing the constantski(r), θi(r), and λi(r), adequately.Remark 4.
Inequality (18) in Theorem 1 is only just sufficient but not a necessary condition for single-weight Markovian jumping complex networks (9) and (13) to realize projection lag synchronization and parameter identification.


### 3.3. Simulation Results

Some numerical simulations are performed to prove the feasibility and effectiveness of the abovementioned result. The Chen system is taken as the node dynamic of the drive complex network:(29)x˙i1tx˙i2tx˙i3t=0−xi1txi3txi1txi2t+xi2t−xi1t00−xi1t0xi1t+xi2t0−xi3t0αi1αi2αi3.The identifying value of the parameter vectorαi is αi=αi1,αi2,αi3T=35,3,28T.The four-dimensional hyperchaotic Lüchaotic system will be the node dynamic of the response complex network:(30)x˙i1tx˙i2tx˙i3tx˙i4t=0−xi1txi3txi1txi2txi1txi3t+xi2t−xi1t00xi4t00xi2t00−xi3t00000xi4tβi1βi2βi3βi4.And the identifying value of the parameter vectorβi is βi=β1,β2,β3,β4T=36,3,20,1T.The size of network is taken asN = 10. The coupling configuration matrices of the network are separately given, as shown in Figures 2 and 3. Figure 2 is the topological structure of the driving system with mode 1 and mode 2. Figure 3 is the response system. Among them, the structure of network topology changing with two modes is shown with the dotted line.Figure 2
(a) Topological structure of the driving system in mode 1 and (b) topological structure of the driving system in mode 2.
(a)
(b)Figure 3
(a) Topological structure of the response system in mode 1 and (b) topological structure of the response system in mode 2.
(a)
(b)In order to judge the quality of the network synchronization, we define the system error value as follows:(31)Et=∑i=1N∑j=1Neij2tN.Whenlimt⟶∞Et=0 and ∀t > t1, systems (9) and (13) will implement global synchronization within the finite time t1.The other relevant parameters are as below:(32)Γ1=0.10000.10000.1,Γ2=0.200000.200000.200000.2,M=0.20−0.10.1−0.55.5−2−13−0.51.2−2,σit,eit,1=diag22ei1t,22ei1t,22ei1t,σit,eit,2=diag2ei1t,2ei1t,2ei1t.The switching of two different modes in the complex system can be shown at Figure4. α^i0∈R3×1 and β^i0∈R4×1 are selected arbitrary integer in the interval 0,1 and τ = 0.01. The parameters of controller (17) are chosen as follows: ki(1) = ki(2) = 15, γ = 0.6, and λi(1) = λi(2) = 12. Through Lemma 2 and (28), we can obtain t1 = 7.4729 by simple calculation.Figure 4
Switching of the system mode.The time-varying curves of the system synchronization errorsei(t)(1 ≤ i ≤ 4) are demonstrated in Figure 5. From that, one can see that the error converge to zero which means that the two networks are in the synchronization state. Figure 6 is the system error values E(t). The conclusion of Figure 6 matches that of Figure 5. Figures 7 and 8 give the identification of parameters α^ii=1,2,3 and β^ii=1,2,3,4. It is clear that the estimated parameters are approached to the value (35, 3, 28)T and (36, 3, 20, 1)T. They all attest that the control scheme can realize the identification and synchronization of two single-weight Markonian jumping complex networks, effectively.Figure 5
Synchronization errorsei(t) with time evolution.
(a)
(b)
(c)
(d)Figure 6
Synchronization error valueE(t) with varying time.Figure 7
Identification of the uncertain parameterα^it.
(a)
(b)
(c)Figure 8
Identification of the uncertain parameterβ^it.
(a)
(b)
(c)
(d)


## 3.1. Network Models

The drive-coupled complex network is described as follows, and its number of dynamic nodes isN:(9)x˙t=fi1xitαi+fi2xit+∑j=1NaijrtΓ1xjt,where xit=xi1t,xi2t,…,xintT∈Rn are the state vectors of the ith node; fi1: Rn ⟶ Rn×l is a continuous matrix function and fi2: Rn ⟶ Rn is a continuous vector function, respectively; αi ∈ Rl are the unidentified node dynamic parameters; Γ1 ∈ Rn×n represents the inner coupling matrix; the coupling configuration matrix is expressed by Art=aijrtN×N∈Rn×n, which represents the topological structure and the coupling strength of the network at time t in mode r(t).{r(t, t ≥ 0)} is defined as a right-continuous Markovian process in probability space, describing the switching state between different parameters of time t. It takes values from the finite space S=1,2,…,m with generator Π=πpqm×mp,q∈S. The transition probability from pth mode at time t to the qth mode at time t + Δt will be defined as follows:(10)Prt+Δt=q|rt=p=πpqΔt+oΔt,ifq≠p,1+πpqΔt+oΔt,ifq=p,where Δt > 0:(11)limΔt⟶0oΔtΔt=0.As the transition rateπpq ≥ 0 which can be represented from mode p at time t to mode q at time t + Δt,(12)πpp=−∑q=1,q≠pmπpq.Therefore, the elements of matrixA(r(t)) will be defined as follows. If it exists as a link from node j to node i (j ≠ i), then aij ≠ 0; otherwise, aij = 0. The diagonal elements of matrix A(r(t)) can be determined by aii=−∑j=1,j≠iNaij.Describe the response-coupled network as follows:(13)y˙it=gi1yitβi+gi2yit+∑j=1NbijrtΓ2yit+uit+σit,eit,rtω˙t,where yit=yi1t,yi2t,…,yimtT∈Rm are the state vectors of node i; gi1: Rm ⟶ Rm×o is a continuous matrix function and gi2: Rm ⟶ Rm is a continuous vector function; the node dynamic parameter vectors which need to be identified are βi ∈ Ro; Γ2 ∈ Rm×m represents the inner coupling matrix; the coupling configuration matrix is defined by Brt=bijrtM×M∈Rm×m as well as in network (9); ui(t) is the finite-time nonlinear feedback controller which needs to be devised; σi(t, ei(t), and r(t)) are the noisy intensity function which mainly describe the impact of environmental fluctuations or imprecise design of coupling strength on network synchronization; and ei(t) is the error of system synchronization.The error of system synchronization can be defined as follows:(14)eit=yit−Mxit−τ,where M ∈ Rm×n is the scaling matrix (i.e., generalized matrix) consisting of constant; eit=ei1t,ei2t,…,eimtT; and τ is the time delay.Definition 2.
The driving network and response network can be realized synchronization through the designed finite-time nonlinear controllers. As for scaling matrixMij = (mij) ∈ Rm×n, if the condition can be satisfied that each row of the element cannot be 0 at a time, it will hold that(15)limt⟶∞Eeit2=limt⟶∞Eyit−Mxit−τ2=0.


## 3.2. Main Result

In this section, the synchronization and identification of unknown parameters for Markovian jumping complex networks with nonidentical dimensions and stochastic perturbations are investigated by utilizing the control theory of finite-time. The following theorem and remarks are the main results.First, from (9) and (13), the synchronization error dynamic system can be obtained as below:(16)e˙it=y˙it−Mx˙it−τ=gi1yitβi+gi2yit+uit+σit,eit,rtω˙t+∑j=1NbijrtΓ2eit+∑j=1NbijrtΓ2Mxit−τ−Mfi1xit−ταi+fi2xit−τ−M∑j=1NaijrtΓ1xjt−τ.Theorem 1.
With Assumption1 holding, the single-weight Markovian jumping complex network (9) and (13) can be realized synchronization in a finite-time. Meanwhile, with the following controllers and undated laws, the value of uncertain parameters can be simultaneously identified:(17)uit=−gi1yitβ^i−gi2yit+Mfi1xit−τα^i+Mfi2xit−τ−kireit+M∑j=1NaijrtΓ1−bijrtΓ2xjt−τ−ψdiagei1tγ,…,eimtγsigneit,α^˙it=−θifi1Txit−τMTteit+ψθiα^i−αiγsignα^i−αi,β^˙it=λigi1Txiteit−ψλiβ^i−βiγsignβ^i−βi,where ki(r) > 0, 0 < γ < 2, ψ > 0, θi and λi are any positive constant, and sign(·) represents symbolic function.Remark 1.
When Markovian jumping complex network (9) and (13) achieve projection lag synchronization in a finite time, the uncertain parameter vectors α^i and β^i of the system are identified.
If the following inequalities hold,(18)∑q=1mπpqDq−Cp≤0,B^r−Ξr+Θr⊗Im≤0,where D(q) is a positive definite matrix; C(p) is a symmetrical matrix; B^r=B˜rB˜Tr/2; B˜r=Br⊗Γ2; Ξr=diagk1r,k2r,…,kNr; Θr=diagρ1r,ρ2r,…,ρNr; and r ∈ S.
Therefore, the synchronization between single-weight Markovian jumping complex network (9) and (13) with stochastic perturbations can be achieved global asymptotic stability in a finite-time under this condition:(19)t1≤t0+2Vt01−γ/22ψ1−γ,where Vt0=1/2∑i=1NeiTt0eit0+1/2∑i=1N1/θiα˜iTα˜i+1/2∑i=1N1/λiβ˜iTβ˜i.Proof.
The Lyapunov function is defined as below:(20)Vt,et,p=12∑i=1NeiTteit+12∑i=1N1θiα˜iTα˜i+12∑i=1N1λiβ˜iTβ˜i.
By Definition1, we can obtain the differential operator ℒV as follows:(21)ℒVxt,t,p=∑i=1NeiTte˙it+∑i=1N1θiα^˙iTα˜i+∑i=1N1λiβ^˙iTβ˜i+∑q=1mπpq∑i=1NeiTtDqet+∑i=1NtranceσiTt,eit,pσit,eit,p.
Substitute (16) and feedback controllers (17) into (21), resulting in(22)ℒVxt,t,p=∑i=1NeiTtgi1yitβi+σit,eit,rtω˙t+∑j=1NbijrtΓ2eit+∑j=1NbijrtΓ2Mxit−τ+uit+gi2yit−Mfi1xit−ταi+fi2xit−τ−M∑j=1NaijrtΓ1xjt−τ−∑i=1NeiTtMfi1xit−τα˜i+∑i=1NeiTtgi1xitβ˜i−ψ∑i=1N1θiα˜iα^i−αiγsignα^i−αi−ψ∑i=1N1λiβ˜iβ^i−βiγsignβ^i−βi+∑i=1NtranceσiTt,eit,pσit,eit,p+∑q=1mπpq∑i=1NeiTtDqet,uit=−gi1yitβ^i−gi2yit+Mfi1xit−τα^i+Mfi2xit−τ−kireit+M∑j=1NaijrtΓ1−bijrtΓ2xjt−τ−ψdiagei1tγ,…,eimtγsigneit,(23)ℒVxt,t,p=∑i=1NeiTtgi1yitβi+σit,eit,rtω˙t+∑j=1NbijrtΓ2eit+∑j=1NbijrtΓ2Mxit−τ+gi2yit−Mfi1xit−ταi+fi2xit−τ−M∑j=1NaijrtΓ1xjt−τ−gi1yitβ^i−gi2yit+Mfi1xit−τα^i+Mfi2xit−τ−kireit+M∑j=1NaijrtΓ1−bijrtΓ2xjt−τ−ψdiagei1tγ,…,eimtγsigneit,−∑i=1NeiTtMfi1xit−τα˜i+∑i=1NeiTtgi1xitβ˜i−ψ∑i=1N1θiα˜iα^i−αiγsignα^i−αi−ψ∑i=1N1λiβ˜iβ^i−βiγsignβ^i−βi+∑i=1NtranceσiTt,eit,pσit,eit,p+∑q=1mπpq∑i=1NeiTtDqet,(24)ℒVxt,t,p=∑i=1NeiTtMfi1xit−τα^it−αi−gi1yitβ^it−βi+∑j=1NbijrtΓ2eit−kireit−ψdiagei1tγ,…,eimtγsigneit−∑i=1NeiTtMfi1xit−τα˜i+∑i=1NeiTtgi1yitβ˜i−ψ∑i=1N1θiα˜iα^i−αiγsignα^i−αi−ψ∑i=1N1λiβ˜iβ^i−βiγsignβ^i−βi+∑i=1NtranceσiTt,eit,pσit,eit,p+∑q=1mπpq∑i=1NeiTtDqet,(25)ℒVxt,t,p=∑i=1NeiTt∑j=1NbijrtΓ2eit−kireit−ψdiagei1tγ,…,eimtγsigneit+∑i=1NtranceσiTt,eit,pσit,eit,p−ψ∑i=1N1θiα˜iα^i−αiγsignα^i−αi−ψ∑i=1N1λiβ˜iβ^i−βiγsignβ^i−βi+∑q=1mπpq∑i=1NeiTtDqet.
According to Lemma3, we obtain that(26)∑i=1NeiTteitφsigneit=∑i=1NeiTtdiagei1tγ,…,eimtγsigneit=∑i=1N∑j=1Neijteijtγsigneijt=∑i=1N∑j=1Neijt1+γ≥∑i=1N∑j=1Neijt21+γ/2=∑i=1NeiTteit1+γ/2⟹∑i=1NeiTteitφsigneit≥∑i=1NeiTteit1+γ/2.
From Assumption1, the conclusion is drawn:(27)ℒVxt,t,p≤∑i=1NeiTtB^r−Ξr+Θr⊗Imeit+∑q=1mπpq∑i=1NeiTtDq−Cpet−ψ∑i=1NeiTteit1+γ/2−ψ∑i=1N1θiα˜iTα˜i1+γ/2−ψ∑i=1N1λiβ˜iTβ˜i1+γ/2≤−ψ∑i=1NeiTteit+∑i=1N1θiα˜iTα˜i+∑i=1N1λiβ˜iTβ˜i1+γ/2≤−2ψ12∑i=1NeiTteit+12∑i=1N1θiα˜iTα˜i+12∑i=1N1λiβ˜iTβ˜i1+γ/2=−2ψVt1+γ/2⟹ℒVxt,t,p≤−2ψVt1+γ/2.
From Lemma2, the system error of synchronization ei(t) will converge to a steady state with a finite-time, which can be estimated by(28)t1≤t0+Vt01−1+γ/22ψ1−1+γ/2=t0+2Vt01−γ/22ψ1−γ⟹t1≤t0+2Vt01−γ/22ψ1−γ,where Vt0=1/2∑i=1NeiTt0eit0+1/2∑i=1N1/θiα˜iTα˜i+1/2∑i=1N1/λiβ˜iTβ˜i.
Hence, networks (9) and (13) can be realized synchronization with finite time t1. Also, the uncertain system parameter vectors α^i and β^i are adapted to the true value of parameters. The proof is completed.Remark 2.
The proposed approach is applicable to the finite-time-generalized matrix projection lag synchronization of any two single-weight Markovian jumping complex networks with uncertain parameters, stochastic perturbations, and different initial conditions. In addition, the Markovian jumping complex network can be identical or nonidentical dimensions.Remark 3.
The synchronization and identification speed can be adjusted through choosing the constantski(r), θi(r), and λi(r), adequately.Remark 4.
Inequality (18) in Theorem 1 is only just sufficient but not a necessary condition for single-weight Markovian jumping complex networks (9) and (13) to realize projection lag synchronization and parameter identification.


## 3.3. Simulation Results

Some numerical simulations are performed to prove the feasibility and effectiveness of the abovementioned result. The Chen system is taken as the node dynamic of the drive complex network:(29)x˙i1tx˙i2tx˙i3t=0−xi1txi3txi1txi2t+xi2t−xi1t00−xi1t0xi1t+xi2t0−xi3t0αi1αi2αi3.The identifying value of the parameter vectorαi is αi=αi1,αi2,αi3T=35,3,28T.The four-dimensional hyperchaotic Lüchaotic system will be the node dynamic of the response complex network:(30)x˙i1tx˙i2tx˙i3tx˙i4t=0−xi1txi3txi1txi2txi1txi3t+xi2t−xi1t00xi4t00xi2t00−xi3t00000xi4tβi1βi2βi3βi4.And the identifying value of the parameter vectorβi is βi=β1,β2,β3,β4T=36,3,20,1T.The size of network is taken asN = 10. The coupling configuration matrices of the network are separately given, as shown in Figures 2 and 3. Figure 2 is the topological structure of the driving system with mode 1 and mode 2. Figure 3 is the response system. Among them, the structure of network topology changing with two modes is shown with the dotted line.Figure 2
(a) Topological structure of the driving system in mode 1 and (b) topological structure of the driving system in mode 2.
(a)
(b)Figure 3
(a) Topological structure of the response system in mode 1 and (b) topological structure of the response system in mode 2.
(a)
(b)In order to judge the quality of the network synchronization, we define the system error value as follows:(31)Et=∑i=1N∑j=1Neij2tN.Whenlimt⟶∞Et=0 and ∀t > t1, systems (9) and (13) will implement global synchronization within the finite time t1.The other relevant parameters are as below:(32)Γ1=0.10000.10000.1,Γ2=0.200000.200000.200000.2,M=0.20−0.10.1−0.55.5−2−13−0.51.2−2,σit,eit,1=diag22ei1t,22ei1t,22ei1t,σit,eit,2=diag2ei1t,2ei1t,2ei1t.The switching of two different modes in the complex system can be shown at Figure4. α^i0∈R3×1 and β^i0∈R4×1 are selected arbitrary integer in the interval 0,1 and τ = 0.01. The parameters of controller (17) are chosen as follows: ki(1) = ki(2) = 15, γ = 0.6, and λi(1) = λi(2) = 12. Through Lemma 2 and (28), we can obtain t1 = 7.4729 by simple calculation.Figure 4
Switching of the system mode.The time-varying curves of the system synchronization errorsei(t)(1 ≤ i ≤ 4) are demonstrated in Figure 5. From that, one can see that the error converge to zero which means that the two networks are in the synchronization state. Figure 6 is the system error values E(t). The conclusion of Figure 6 matches that of Figure 5. Figures 7 and 8 give the identification of parameters α^ii=1,2,3 and β^ii=1,2,3,4. It is clear that the estimated parameters are approached to the value (35, 3, 28)T and (36, 3, 20, 1)T. They all attest that the control scheme can realize the identification and synchronization of two single-weight Markonian jumping complex networks, effectively.Figure 5
Synchronization errorsei(t) with time evolution.
(a)
(b)
(c)
(d)Figure 6
Synchronization error valueE(t) with varying time.Figure 7
Identification of the uncertain parameterα^it.
(a)
(b)
(c)Figure 8
Identification of the uncertain parameterβ^it.
(a)
(b)
(c)
(d)


## 4. Finite-Time-Generalized Matrix Projective Synchronization with Double-Weight Networks

In order to build a practical Markovian jumping complex network model with double-weights, this paper introduces a method of network splitting, which divides the different coupling configuration into different subnetworks. Each network after splitting has its own nature and structure. In this paper, different side voltage levels in the power network are defined as different properties, which are divided into two subnetworks (as shown in Figure9). Then, the finite-time generalized function lag projection synchronization is studied. The synchronization criteria are given and the values of the related unknown parameters are identified further.Figure 9
Network splitting diagram.


### 4.1. Network Models

In this section, the driving system is defined as follows and its dynamic node number isN:(33)xt=fi1xitαi+fi2xit+∑j=1NaijrtΓ1xjt+∑j=1NbijrtΓ2xjt,where Γi(i = 1, 2) are the inner coupling matrices; Art=aijrtN×N and Brt=bijrtN×N represent the coupling configuration matrices which are represented as different coupling relationships. Its elements satisfy the following. If there exists a link from node j to node i (j ≠ i), then aij ≠ 0(bij ≠ 0); otherwise, aij = 0(bij = 0). The diagonal elements of matrix A(r(t))(B(r(t))) are defined as(34)aii=−∑j=1,j≠iNaijbii=−∑j=1,j≠iNbij.The response coupled network can be described as follows:(35)y˙it=gi1yitβi+gi2yit+∑j=1NmijrtΓ3yit+∑j=1NnijrtΓ4yit+uit+σit,eit,rtω˙t,where Γi(i = 3, 4) represent the inner coupling matrices; Mrt=mijrtN×N and Nrt=nijrtN×N are the coupling configuration matrices, same as A(r(t))(B(r(t))) in (33); ui(t) is the finite-time nonlinear feedback controller which needs to be designed; and σi(t, ei(t), r(t)) is the noisy intensity function and ei(t) is the synchronization error.The system synchronization error is defined as follows:(36)eit=yit−Mtxit−τ,where M(t) ∈ Rm×n is the scaling matrix consisting of functions (i.e., generalized function matrix).Definition 3.
The driving network (33) and response network (35) can be realized synchronization with the designed finite-time nonlinear controllers. As for scaling matrix M(t) = (mij(t)) ∈ Rm×n, if the condition can be satisfied that the each row of element cannot be 0 at a time, it will hold that(37)limt⟶∞Eeit2=limt⟶∞Eyit−Mtxit−τ2=0.


### 4.2. Main Result

The synchronization error of the dynamic system is obtained as the following from (33) and (35):(38)e˙it=y˙it−Mx˙it−τ−M˙xit−τ=gi1yitβi+gi2yit+uit+σit,eit,rtω˙t+∑j=1NmijrtΓ3eit+∑j=1NmijrtΓ3Mtxit−τ+∑j=1NnijrtΓ4eit−M˙txit−τ+∑j=1NnijrtΓ4Mtxit−τ−Mtfi1xit−ταi+fi2xit−τ−Mt∑j=1NaijrtΓ1xjt−τ−Mt∑j=1NbijrtΓ2xjt−τ.Theorem 2.
With Assumption1 holding, the double-weights Markovian jumping complex networks (33) and (35) can be realized synchronization in a finite-time. Meanwhile, with following controllers and undated laws, the value of uncertain network parameters can be simultaneously identified:(39)uit=M˙txit−τ−gi1yitβ^i−gi2yit+Mfi1xit−τα^i+Mfi2xit−τ−kireit+∑j=1NaijrtMtΓ1−mijrtMtΓ3xjt−τ+∑j=1NbijrtMtΓ2−nijrtMtΓ4xjt−τ−ψdiagei1tγ,…,eimtγsigneit,α^˙it=−hirfi1Txit−τMTteit+ψhirα^i−αiγsignα^i−αi,β^˙it=sirgi1Txiteit−ψsirβ^i−βiγsignβ^i−βi.Remark 5.
When Markovian jumping complex networks (33) and (35) achieve projection lag synchronization in a finite time, the uncertain system parameter vectors α^i and β^i are identified.
If the following inequalities hold,(40)∑q=1mπpqDq−Cp≤0,Q˜r+P˜r+Θr−Ξr⊗Im≤0,where D(q) is a positive definite matrix; C(p) is a symmetrical matrix; Q˜r=Qr⊗Γ3; P˜r=Pr⊗Γ4; Θr=diagρ1r,ρ2r,…,ρNr; Ξr=diagk1r,k2r,…,kNr; and r ∈ S.
Therefore, the synchronization between double-weights Markovian jumping complex networks (33) and (35) with stochastic perturbations can be achieved global asymptotic stability in a finite-time under this condition:(41)t1≤t0+2Vt01−γ/22ψ1−γ,where Vt0=1/2∑i=1NeiTt0eit0+1/2∑i=1N1/hirα˜iTα˜i+1/2∑i=1N1/sirβ˜iTβ˜i.Proof.
The Lyapunov function is defined as follows:(42)Vt,et,p=12∑i=1NeiTteit+12∑i=1N1hipα˜iTα˜i+12∑i=1N1sipβ˜iTβ˜i.
By Definition1, the differential operator ℒV is given as(43)ℒVxt,t,p=∑i=1NeiTte˙it+∑i=1N1hipα^˙iTα˜i+∑i=1N1sipβ^˙iTβ˜i+∑q=1mπpq∑i=1NeiTtDqet+∑i=1NtranceσiTt,eit,pσit,eit,p.
Substitute (38) and feedback controllers (39) into (43), resulting in(44)ℒVxt,t,p=∑i=1NeiTtgi1yitβi+gi2yit+uit+σit,eit,rtω˙t+∑j=1NmijrtΓ3eit+∑j=1NmijrtΓ3Mtxit−τ+∑j=1NnijrtΓ4eit+∑j=1NnijrtΓ4Mtxit−τ−Mtfi1xit−ταi+fi2xit−τ−Mt∑j=1NaijrtΓ1xjt−τ−Mt∑j=1NbijrtΓ2xjt−τ−∑i=1NeiTtMfi1xit−τα˜i+∑i=1NeiTtgi1xitβ˜i−ψ∑i=1N1hirα˜iα^i−αiγsignα^i−αi−ψ∑i=1N1sirβ˜iβ^i−βiγsignβ^i−βi+∑i=1NtranceσiTt,eit,pσit,eit,p+∑q=1mπpq∑i=1NeiTtDqetuit=M˙txit−τ−gi1yitβ^i−gi2yit+Mfi1xit−τα^i+Mfi2xit−τ−kireit+∑j=1NaijrtMtΓ1−mijrtMtΓ3xjt−τ+∑j=1NbijrtMtΓ2−nijrtMtΓ4xjt−τ−ψdiagei1tγ,…,eimtγsigneit,(45)ℒVxt,t,p=∑i=1NeiTtgi1yitβi+gi2yit+σit,eit,rtω˙t+∑j=1NmijrtΓ3eit∑j=1NmijrtΓ3Mtxit−τ+∑j=1NnijrtΓ4eit+∑j=1NnijrtΓ4Mtxit−τ−Mtfi1xit−ταi+fi2xit−τ−Mt∑j=1NaijrtΓ1xjt−τ−Mt∑j=1NbijrtΓ2xjt−τM˙txit−τ−gi1yitβ^i−gi2yit+Mfi1xit−τα^i+Mfi2xit−τ−kireit+∑j=1NaijrtMtΓ1−mijrtMtΓ3xjt−τ+∑j=1NbijrtMtΓ2−nijrtMtΓ4−xjt−τ−ψdiagei1tγ,…,eimtγsigneit−∑i=1NeiTtMfi1xit−τα˜i+∑i=1NeiTtgi1xitβ˜i−ψ∑i=1N1hirα˜iα^i−αiγsignα^i−αi−ψ∑i=1N1sirβ˜iβ^i−βiγsignβ^i−βi+∑i=1NtranceσiTt,eit,pσit,eit,p+∑q=1mπpq∑i=1NeiTtDqet,(46)ℒVxt,t,p=∑i=1NeiTtMtfi1xit−τα^it−αi−gi1yitβ^it−βi−kireit+∑j=1NmijrtΓ3+nijrtΓ4eit−ψdiagei1tγ,…,eimtγsigneit−∑i=1NeiTtMtfi1xit−τα˜i+∑i=1NeiTtgi1yitβ˜i−ψ∑i=1N1hirα˜iα^i−αiγsignα^i−αi−ψ∑i=1N1sirβ˜iβ^i−βiγsignβ^i−βi+∑i=1NtranceσiTt,eit,pσit,eit,p+∑q=1mπpq∑i=1NeiTtDqet,(47)ℒVxt,t,p=∑i=1NeiTtMtfi1xit−τα^it−αi−∑i=1NeiTtgi1yitβ^it−βi+∑i=1NeiTt−kireit+∑j=1NmijrtΓ3+nijrtΓ4eit−ψdiagei1tγ,…,eimtγsigneit−∑i=1NeiTtMtfi1xit−τα˜i+∑i=1NeiTtgi1yitβ˜i−ψ∑i=1N1hirα˜iα^i−αiγsignα^i−αi−ψ∑i=1N1sirβ˜iβ^i−βiγsignβ^i−βi+∑i=1NtranceσiTt,eit,pσit,eit,p+∑q=1mπpq∑i=1NeiTtDqet,(48)ℒVxt,t,p=∑i=1NeiTt∑j=1NmijrtΓ3eit+∑j=1NnijrtΓ4eit−kireit−ψdiagei1tγ,…,eimtγsigneit+∑i=1NtranceσiTt,eit,pσit,eit,p−ψ∑i=1N1hirα˜iα^i−αiγsignα^i−αi−ψ∑i=1N1sirβ˜iβ^i−βiγsignβ^i−βi+∑q=1mπpq∑i=1NeiTtDqet.
According to Lemma3,(49)∑i=1NeiTteitφsigneit=∑i=1NeiTtdiagei1tγ,…,eimtγsigneit=∑i=1N∑j=1Neijteijtγsigneijt=∑i=1N∑j=1Neijt1+γ≥∑i=1N∑j=1Neijt21+γ/2=∑i=1NeiTteit1+γ/2⟹∑i=1NeiTteitφsigneit≥∑i=1NeiTteit1+γ/2.
From Assumption1, the conclusion is drawn:(50)ℒVxt,t,p≤eTtQ˜r+P˜r+Θr−Ξr⊗Imet+∑q=1mπpq∑i=1NeiTtDq−Cpet−ψ∑i=1NeiTteit1+γ/2−ψ∑i=1N1hirα˜iTα˜i1+γ/2−ψ∑i=1N1sirβ˜iTβ˜i1+γ/2≤−ψ∑i=1NeiTteit+∑i=1N1θiα˜iTα˜i+∑i=1N1λiβ˜iTβ˜i1+γ/2≤−2ψ12∑i=1NeiTteit+12∑i=1N1θiα˜iTα˜i+12∑i=1N1λiβ˜iTβ˜i1+γ/2=−2ψVt1+γ/2⟹ℒVxt,t,p≤−2ψVt1+γ/2.
From Lemma2, the system error of synchronization ei(t) will converge to a steady state with a finite-time which can be estimated by(51)t1≤t0+Vt01−1+γ/22ψ1−1+γ/2=t0+2Vt01−γ/22ψ1−γ⟹t1≤t0+2Vt01−γ/22ψ1−γ,where Vt0=1/2∑i=1NeiTt0eit0+1/2∑i=1N1/hirα˜iTα˜i+1/2∑i=1N1/sirβ˜iTβ˜i.
Hence, networks (33) and (35) could be realized synchronization with finite time t1. As well as, the unknown parameters vector α^i and β^i of the system are adapted to the true value with themselves. The proof is completed.Remark 6.
The proposed approach is applicable to the finite-time-generalized matrix projection lag synchronization of any two double-weights Markovian jumping complex networks with uncertain parameters, stochastic perturbations, and different initial conditions. In addition, the dimensions of the Markovian jumping complex network can be identical or nonidentical.Remark 7.
The synchronization and identification speed can be adjusted by choosing constantski(r), hi(r), and si(r), adequately.Remark 8.
Inequality (40) in Theorem 2 is only just the sufficient but not necessary condition for double-weights Markovian jumping complex networks (33) and (35) to realize projection lag synchronization and parameter identification.


### 4.3. Simulation Results

#### 4.3.1. Simulation Analysis I

The Chen system is taken as the node dynamic of drive complex network:(52)x˙i1tx˙i2tx˙i3t=0−xi1txi3txi1txi2t+xi2t−xi1t00−xi1t0xi1t+xi2t0−xi3t0αi1αi2αi3.The identifying value of the parameter vectorαi is αi=αi1,αi2,αi3T=35,3,28T.The four-dimensional Lüchaotic system will be the node dynamic of the response complex network:(53)x˙i1tx˙i2tx˙i3tx˙i4t=0−xi1txi3txi1txi2txi1txi3t+xi2t−xi1t00xi4t00xi2t00−xi3t00000xi4tβi1βi2βi3βi4.The identifying value of the parameter vectorβi is βi=β1,β2,β3,β4T=36,3,20,1T.The size of network is taken asN = 10. And the coupling configuration matrices of the network are separately given as Figures 10 and 11. Figure 10 is the topological structure of the driving system with mode 1 and mode 2. And Figure 11 is the response systems. The complex coupling network is divided into two simple coupling structures according to the theory of network splitting. Among the two subgraphs, the structure of network topology change with two modes is shown with the dotted line. The switching of the system mode is shown in Figure 12.Figure 10
(a) Topological structure of the driving system in mode 1 and (b) topological structure of the driving system in mode 2.
(a)
(b)Figure 11
(a) Topological structure of the response system in mode 1 and (b) topological structure of the response system in mode 2.
(a)
(b)Figure 12
Switching of the system mode.The other relevant parameters are as below:(54)Γ1=100010001,Γ2=0.50000.50000.5,Γ3=0.200000.200000.200000.2,Γ4=0.500000.500000.500000.5,M=1+sint2−cos2t−13+sin3t0.53cos2t−3cos3t2sint1−23sin3t,σit,eit,1=diag22ei1t,22ei1t,22ei1t,σit,eit,2=diag2ei1t,2ei1t,2ei1t,τ = 0.01, t0 = 0 s, ψ = 15, ki(1) = ki(2) = 15, γ = 0.6, hi(1) = hi(2) = 10, and si(1) = si(2) = 12. Through Lemma 2 and (51), we can obtain t1 = 14.427 by simple calculation.The time-varying curves of the system synchronization errorsei(t)(1 ≤ i ≤ 4) are demonstrated in Figure 13 from which one can see that the error converge to zero. Figures 14 and 15 give the identification of unknown parameters α^ii=1,2,3 and β^i=i=1,2,3,4. It is clear that the estimated parameters are approached to values (35, 3, 28)T and (36, 3, 20, 1)T. Figure 16 is the system error values E(t). The conclusion of Figure 16 matches that of Figure 13. The all attest that the control scheme can realize the identification and synchronization of two double-weights Markonian jumping complex networks, effectively.Figure 13
Synchronization errorsei(t) with time evolution.
(a)
(b)
(c)
(d)Figure 14
Identification of the uncertain parametersα^it.
(a)
(b)
(c)Figure 15
Identification of the uncertain parametersβ^it.
(a)
(b)
(c)
(d)Figure 16
Synchronization error valueE(t) with varying time.


#### 4.3.2. Simulation Analysis II

In order to further verify that the synchronization criterion can be applied to the synchronization research of double-weights Markovian jumping complex networks, different driving systems are selected for simulation analysis.The Liu system is taken as the node dynamic of the drive complex network. Among them,k = 1 and h = 4:(55)x˙i1tx˙i2tx˙i3t=0−kxi1txi3thxi1t2+xi2t−xi1t000xi1t000−xi3tαi1αi2αi3.The identifying value of the parameter vectorαi is αi=αi1,αi2,αi3T=10,40,2.5T.The four-dimensional hyperchaotic Lüchaotic system will be the node dynamic of the response complex network.Other relevant parameters are the same as those in Section4.3.1. t1 = 17.806 is calculated.The time-varying curves of the system synchronization errorsei(t)(1 ≤ i ≤ 4) are demonstrated in Figure 17. From that, one can see that the error converge to zero. Figure 18 is the identification of unknown parameters α^ii=1,2,3. It is clear that the estimated parameters are approached to the value (10, 40, 2.5)T.Figure 17
Synchronization errorsei(t) with time evolution.
(a)
(b)
(c)
(d)Figure 18
Identification of the uncertain parametersα^it.
(a)
(b)
(c)Under the same research background and system parameters, the control scheme proposed in this paper is compared with the control scheme without finite time to demonstrate the better performance of the proposed control scheme. Figures19 and 20 are the simulation results under the control scheme without finite time. Figure 19 reveals the time evolution curves of synchronization error ei(t). Figure 20 indicates that estimated parameter α^i converge to their truth values (10, 40, 2.5)T. By comparing with the proposed control scheme in Figures 17 and 18, it is obvious that the proposed control scheme can achieve synchronization and identification more effectively and cost less time. Figure 21 is the synchronization error value E(t) with time-varying under the control scheme with finite time and without finite time.Figure 19
Synchronization errorsei(t) with time evolution under the control scheme without finite time.
(a)
(b)
(c)
(d)Figure 20
Identification of the uncertain parameterα^it under the control scheme without finite time.
(a)
(b)
(c)Figure 21
Synchronization error valueE(t) with time-varying under the control scheme with finite time and without finite time.These illustrate that the synchronization and identification of the drive and response networks have been effectively achieved in finite time.


## 4.1. Network Models

In this section, the driving system is defined as follows and its dynamic node number isN:(33)xt=fi1xitαi+fi2xit+∑j=1NaijrtΓ1xjt+∑j=1NbijrtΓ2xjt,where Γi(i = 1, 2) are the inner coupling matrices; Art=aijrtN×N and Brt=bijrtN×N represent the coupling configuration matrices which are represented as different coupling relationships. Its elements satisfy the following. If there exists a link from node j to node i (j ≠ i), then aij ≠ 0(bij ≠ 0); otherwise, aij = 0(bij = 0). The diagonal elements of matrix A(r(t))(B(r(t))) are defined as(34)aii=−∑j=1,j≠iNaijbii=−∑j=1,j≠iNbij.The response coupled network can be described as follows:(35)y˙it=gi1yitβi+gi2yit+∑j=1NmijrtΓ3yit+∑j=1NnijrtΓ4yit+uit+σit,eit,rtω˙t,where Γi(i = 3, 4) represent the inner coupling matrices; Mrt=mijrtN×N and Nrt=nijrtN×N are the coupling configuration matrices, same as A(r(t))(B(r(t))) in (33); ui(t) is the finite-time nonlinear feedback controller which needs to be designed; and σi(t, ei(t), r(t)) is the noisy intensity function and ei(t) is the synchronization error.The system synchronization error is defined as follows:(36)eit=yit−Mtxit−τ,where M(t) ∈ Rm×n is the scaling matrix consisting of functions (i.e., generalized function matrix).Definition 3.
The driving network (33) and response network (35) can be realized synchronization with the designed finite-time nonlinear controllers. As for scaling matrix M(t) = (mij(t)) ∈ Rm×n, if the condition can be satisfied that the each row of element cannot be 0 at a time, it will hold that(37)limt⟶∞Eeit2=limt⟶∞Eyit−Mtxit−τ2=0.


## 4.2. Main Result

The synchronization error of the dynamic system is obtained as the following from (33) and (35):(38)e˙it=y˙it−Mx˙it−τ−M˙xit−τ=gi1yitβi+gi2yit+uit+σit,eit,rtω˙t+∑j=1NmijrtΓ3eit+∑j=1NmijrtΓ3Mtxit−τ+∑j=1NnijrtΓ4eit−M˙txit−τ+∑j=1NnijrtΓ4Mtxit−τ−Mtfi1xit−ταi+fi2xit−τ−Mt∑j=1NaijrtΓ1xjt−τ−Mt∑j=1NbijrtΓ2xjt−τ.Theorem 2.
With Assumption1 holding, the double-weights Markovian jumping complex networks (33) and (35) can be realized synchronization in a finite-time. Meanwhile, with following controllers and undated laws, the value of uncertain network parameters can be simultaneously identified:(39)uit=M˙txit−τ−gi1yitβ^i−gi2yit+Mfi1xit−τα^i+Mfi2xit−τ−kireit+∑j=1NaijrtMtΓ1−mijrtMtΓ3xjt−τ+∑j=1NbijrtMtΓ2−nijrtMtΓ4xjt−τ−ψdiagei1tγ,…,eimtγsigneit,α^˙it=−hirfi1Txit−τMTteit+ψhirα^i−αiγsignα^i−αi,β^˙it=sirgi1Txiteit−ψsirβ^i−βiγsignβ^i−βi.Remark 5.
When Markovian jumping complex networks (33) and (35) achieve projection lag synchronization in a finite time, the uncertain system parameter vectors α^i and β^i are identified.
If the following inequalities hold,(40)∑q=1mπpqDq−Cp≤0,Q˜r+P˜r+Θr−Ξr⊗Im≤0,where D(q) is a positive definite matrix; C(p) is a symmetrical matrix; Q˜r=Qr⊗Γ3; P˜r=Pr⊗Γ4; Θr=diagρ1r,ρ2r,…,ρNr; Ξr=diagk1r,k2r,…,kNr; and r ∈ S.
Therefore, the synchronization between double-weights Markovian jumping complex networks (33) and (35) with stochastic perturbations can be achieved global asymptotic stability in a finite-time under this condition:(41)t1≤t0+2Vt01−γ/22ψ1−γ,where Vt0=1/2∑i=1NeiTt0eit0+1/2∑i=1N1/hirα˜iTα˜i+1/2∑i=1N1/sirβ˜iTβ˜i.Proof.
The Lyapunov function is defined as follows:(42)Vt,et,p=12∑i=1NeiTteit+12∑i=1N1hipα˜iTα˜i+12∑i=1N1sipβ˜iTβ˜i.
By Definition1, the differential operator ℒV is given as(43)ℒVxt,t,p=∑i=1NeiTte˙it+∑i=1N1hipα^˙iTα˜i+∑i=1N1sipβ^˙iTβ˜i+∑q=1mπpq∑i=1NeiTtDqet+∑i=1NtranceσiTt,eit,pσit,eit,p.
Substitute (38) and feedback controllers (39) into (43), resulting in(44)ℒVxt,t,p=∑i=1NeiTtgi1yitβi+gi2yit+uit+σit,eit,rtω˙t+∑j=1NmijrtΓ3eit+∑j=1NmijrtΓ3Mtxit−τ+∑j=1NnijrtΓ4eit+∑j=1NnijrtΓ4Mtxit−τ−Mtfi1xit−ταi+fi2xit−τ−Mt∑j=1NaijrtΓ1xjt−τ−Mt∑j=1NbijrtΓ2xjt−τ−∑i=1NeiTtMfi1xit−τα˜i+∑i=1NeiTtgi1xitβ˜i−ψ∑i=1N1hirα˜iα^i−αiγsignα^i−αi−ψ∑i=1N1sirβ˜iβ^i−βiγsignβ^i−βi+∑i=1NtranceσiTt,eit,pσit,eit,p+∑q=1mπpq∑i=1NeiTtDqetuit=M˙txit−τ−gi1yitβ^i−gi2yit+Mfi1xit−τα^i+Mfi2xit−τ−kireit+∑j=1NaijrtMtΓ1−mijrtMtΓ3xjt−τ+∑j=1NbijrtMtΓ2−nijrtMtΓ4xjt−τ−ψdiagei1tγ,…,eimtγsigneit,(45)ℒVxt,t,p=∑i=1NeiTtgi1yitβi+gi2yit+σit,eit,rtω˙t+∑j=1NmijrtΓ3eit∑j=1NmijrtΓ3Mtxit−τ+∑j=1NnijrtΓ4eit+∑j=1NnijrtΓ4Mtxit−τ−Mtfi1xit−ταi+fi2xit−τ−Mt∑j=1NaijrtΓ1xjt−τ−Mt∑j=1NbijrtΓ2xjt−τM˙txit−τ−gi1yitβ^i−gi2yit+Mfi1xit−τα^i+Mfi2xit−τ−kireit+∑j=1NaijrtMtΓ1−mijrtMtΓ3xjt−τ+∑j=1NbijrtMtΓ2−nijrtMtΓ4−xjt−τ−ψdiagei1tγ,…,eimtγsigneit−∑i=1NeiTtMfi1xit−τα˜i+∑i=1NeiTtgi1xitβ˜i−ψ∑i=1N1hirα˜iα^i−αiγsignα^i−αi−ψ∑i=1N1sirβ˜iβ^i−βiγsignβ^i−βi+∑i=1NtranceσiTt,eit,pσit,eit,p+∑q=1mπpq∑i=1NeiTtDqet,(46)ℒVxt,t,p=∑i=1NeiTtMtfi1xit−τα^it−αi−gi1yitβ^it−βi−kireit+∑j=1NmijrtΓ3+nijrtΓ4eit−ψdiagei1tγ,…,eimtγsigneit−∑i=1NeiTtMtfi1xit−τα˜i+∑i=1NeiTtgi1yitβ˜i−ψ∑i=1N1hirα˜iα^i−αiγsignα^i−αi−ψ∑i=1N1sirβ˜iβ^i−βiγsignβ^i−βi+∑i=1NtranceσiTt,eit,pσit,eit,p+∑q=1mπpq∑i=1NeiTtDqet,(47)ℒVxt,t,p=∑i=1NeiTtMtfi1xit−τα^it−αi−∑i=1NeiTtgi1yitβ^it−βi+∑i=1NeiTt−kireit+∑j=1NmijrtΓ3+nijrtΓ4eit−ψdiagei1tγ,…,eimtγsigneit−∑i=1NeiTtMtfi1xit−τα˜i+∑i=1NeiTtgi1yitβ˜i−ψ∑i=1N1hirα˜iα^i−αiγsignα^i−αi−ψ∑i=1N1sirβ˜iβ^i−βiγsignβ^i−βi+∑i=1NtranceσiTt,eit,pσit,eit,p+∑q=1mπpq∑i=1NeiTtDqet,(48)ℒVxt,t,p=∑i=1NeiTt∑j=1NmijrtΓ3eit+∑j=1NnijrtΓ4eit−kireit−ψdiagei1tγ,…,eimtγsigneit+∑i=1NtranceσiTt,eit,pσit,eit,p−ψ∑i=1N1hirα˜iα^i−αiγsignα^i−αi−ψ∑i=1N1sirβ˜iβ^i−βiγsignβ^i−βi+∑q=1mπpq∑i=1NeiTtDqet.
According to Lemma3,(49)∑i=1NeiTteitφsigneit=∑i=1NeiTtdiagei1tγ,…,eimtγsigneit=∑i=1N∑j=1Neijteijtγsigneijt=∑i=1N∑j=1Neijt1+γ≥∑i=1N∑j=1Neijt21+γ/2=∑i=1NeiTteit1+γ/2⟹∑i=1NeiTteitφsigneit≥∑i=1NeiTteit1+γ/2.
From Assumption1, the conclusion is drawn:(50)ℒVxt,t,p≤eTtQ˜r+P˜r+Θr−Ξr⊗Imet+∑q=1mπpq∑i=1NeiTtDq−Cpet−ψ∑i=1NeiTteit1+γ/2−ψ∑i=1N1hirα˜iTα˜i1+γ/2−ψ∑i=1N1sirβ˜iTβ˜i1+γ/2≤−ψ∑i=1NeiTteit+∑i=1N1θiα˜iTα˜i+∑i=1N1λiβ˜iTβ˜i1+γ/2≤−2ψ12∑i=1NeiTteit+12∑i=1N1θiα˜iTα˜i+12∑i=1N1λiβ˜iTβ˜i1+γ/2=−2ψVt1+γ/2⟹ℒVxt,t,p≤−2ψVt1+γ/2.
From Lemma2, the system error of synchronization ei(t) will converge to a steady state with a finite-time which can be estimated by(51)t1≤t0+Vt01−1+γ/22ψ1−1+γ/2=t0+2Vt01−γ/22ψ1−γ⟹t1≤t0+2Vt01−γ/22ψ1−γ,where Vt0=1/2∑i=1NeiTt0eit0+1/2∑i=1N1/hirα˜iTα˜i+1/2∑i=1N1/sirβ˜iTβ˜i.
Hence, networks (33) and (35) could be realized synchronization with finite time t1. As well as, the unknown parameters vector α^i and β^i of the system are adapted to the true value with themselves. The proof is completed.Remark 6.
The proposed approach is applicable to the finite-time-generalized matrix projection lag synchronization of any two double-weights Markovian jumping complex networks with uncertain parameters, stochastic perturbations, and different initial conditions. In addition, the dimensions of the Markovian jumping complex network can be identical or nonidentical.Remark 7.
The synchronization and identification speed can be adjusted by choosing constantski(r), hi(r), and si(r), adequately.Remark 8.
Inequality (40) in Theorem 2 is only just the sufficient but not necessary condition for double-weights Markovian jumping complex networks (33) and (35) to realize projection lag synchronization and parameter identification.


## 4.3. Simulation Results

### 4.3.1. Simulation Analysis I

The Chen system is taken as the node dynamic of drive complex network:(52)x˙i1tx˙i2tx˙i3t=0−xi1txi3txi1txi2t+xi2t−xi1t00−xi1t0xi1t+xi2t0−xi3t0αi1αi2αi3.The identifying value of the parameter vectorαi is αi=αi1,αi2,αi3T=35,3,28T.The four-dimensional Lüchaotic system will be the node dynamic of the response complex network:(53)x˙i1tx˙i2tx˙i3tx˙i4t=0−xi1txi3txi1txi2txi1txi3t+xi2t−xi1t00xi4t00xi2t00−xi3t00000xi4tβi1βi2βi3βi4.The identifying value of the parameter vectorβi is βi=β1,β2,β3,β4T=36,3,20,1T.The size of network is taken asN = 10. And the coupling configuration matrices of the network are separately given as Figures 10 and 11. Figure 10 is the topological structure of the driving system with mode 1 and mode 2. And Figure 11 is the response systems. The complex coupling network is divided into two simple coupling structures according to the theory of network splitting. Among the two subgraphs, the structure of network topology change with two modes is shown with the dotted line. The switching of the system mode is shown in Figure 12.Figure 10
(a) Topological structure of the driving system in mode 1 and (b) topological structure of the driving system in mode 2.
(a)
(b)Figure 11
(a) Topological structure of the response system in mode 1 and (b) topological structure of the response system in mode 2.
(a)
(b)Figure 12
Switching of the system mode.The other relevant parameters are as below:(54)Γ1=100010001,Γ2=0.50000.50000.5,Γ3=0.200000.200000.200000.2,Γ4=0.500000.500000.500000.5,M=1+sint2−cos2t−13+sin3t0.53cos2t−3cos3t2sint1−23sin3t,σit,eit,1=diag22ei1t,22ei1t,22ei1t,σit,eit,2=diag2ei1t,2ei1t,2ei1t,τ = 0.01, t0 = 0 s, ψ = 15, ki(1) = ki(2) = 15, γ = 0.6, hi(1) = hi(2) = 10, and si(1) = si(2) = 12. Through Lemma 2 and (51), we can obtain t1 = 14.427 by simple calculation.The time-varying curves of the system synchronization errorsei(t)(1 ≤ i ≤ 4) are demonstrated in Figure 13 from which one can see that the error converge to zero. Figures 14 and 15 give the identification of unknown parameters α^ii=1,2,3 and β^i=i=1,2,3,4. It is clear that the estimated parameters are approached to values (35, 3, 28)T and (36, 3, 20, 1)T. Figure 16 is the system error values E(t). The conclusion of Figure 16 matches that of Figure 13. The all attest that the control scheme can realize the identification and synchronization of two double-weights Markonian jumping complex networks, effectively.Figure 13
Synchronization errorsei(t) with time evolution.
(a)
(b)
(c)
(d)Figure 14
Identification of the uncertain parametersα^it.
(a)
(b)
(c)Figure 15
Identification of the uncertain parametersβ^it.
(a)
(b)
(c)
(d)Figure 16
Synchronization error valueE(t) with varying time.


### 4.3.2. Simulation Analysis II

In order to further verify that the synchronization criterion can be applied to the synchronization research of double-weights Markovian jumping complex networks, different driving systems are selected for simulation analysis.The Liu system is taken as the node dynamic of the drive complex network. Among them,k = 1 and h = 4:(55)x˙i1tx˙i2tx˙i3t=0−kxi1txi3thxi1t2+xi2t−xi1t000xi1t000−xi3tαi1αi2αi3.The identifying value of the parameter vectorαi is αi=αi1,αi2,αi3T=10,40,2.5T.The four-dimensional hyperchaotic Lüchaotic system will be the node dynamic of the response complex network.Other relevant parameters are the same as those in Section4.3.1. t1 = 17.806 is calculated.The time-varying curves of the system synchronization errorsei(t)(1 ≤ i ≤ 4) are demonstrated in Figure 17. From that, one can see that the error converge to zero. Figure 18 is the identification of unknown parameters α^ii=1,2,3. It is clear that the estimated parameters are approached to the value (10, 40, 2.5)T.Figure 17
Synchronization errorsei(t) with time evolution.
(a)
(b)
(c)
(d)Figure 18
Identification of the uncertain parametersα^it.
(a)
(b)
(c)Under the same research background and system parameters, the control scheme proposed in this paper is compared with the control scheme without finite time to demonstrate the better performance of the proposed control scheme. Figures19 and 20 are the simulation results under the control scheme without finite time. Figure 19 reveals the time evolution curves of synchronization error ei(t). Figure 20 indicates that estimated parameter α^i converge to their truth values (10, 40, 2.5)T. By comparing with the proposed control scheme in Figures 17 and 18, it is obvious that the proposed control scheme can achieve synchronization and identification more effectively and cost less time. Figure 21 is the synchronization error value E(t) with time-varying under the control scheme with finite time and without finite time.Figure 19
Synchronization errorsei(t) with time evolution under the control scheme without finite time.
(a)
(b)
(c)
(d)Figure 20
Identification of the uncertain parameterα^it under the control scheme without finite time.
(a)
(b)
(c)Figure 21
Synchronization error valueE(t) with time-varying under the control scheme with finite time and without finite time.These illustrate that the synchronization and identification of the drive and response networks have been effectively achieved in finite time.


## 4.3.1. Simulation Analysis I

The Chen system is taken as the node dynamic of drive complex network:(52)x˙i1tx˙i2tx˙i3t=0−xi1txi3txi1txi2t+xi2t−xi1t00−xi1t0xi1t+xi2t0−xi3t0αi1αi2αi3.The identifying value of the parameter vectorαi is αi=αi1,αi2,αi3T=35,3,28T.The four-dimensional Lüchaotic system will be the node dynamic of the response complex network:(53)x˙i1tx˙i2tx˙i3tx˙i4t=0−xi1txi3txi1txi2txi1txi3t+xi2t−xi1t00xi4t00xi2t00−xi3t00000xi4tβi1βi2βi3βi4.The identifying value of the parameter vectorβi is βi=β1,β2,β3,β4T=36,3,20,1T.The size of network is taken asN = 10. And the coupling configuration matrices of the network are separately given as Figures 10 and 11. Figure 10 is the topological structure of the driving system with mode 1 and mode 2. And Figure 11 is the response systems. The complex coupling network is divided into two simple coupling structures according to the theory of network splitting. Among the two subgraphs, the structure of network topology change with two modes is shown with the dotted line. The switching of the system mode is shown in Figure 12.Figure 10
(a) Topological structure of the driving system in mode 1 and (b) topological structure of the driving system in mode 2.
(a)
(b)Figure 11
(a) Topological structure of the response system in mode 1 and (b) topological structure of the response system in mode 2.
(a)
(b)Figure 12
Switching of the system mode.The other relevant parameters are as below:(54)Γ1=100010001,Γ2=0.50000.50000.5,Γ3=0.200000.200000.200000.2,Γ4=0.500000.500000.500000.5,M=1+sint2−cos2t−13+sin3t0.53cos2t−3cos3t2sint1−23sin3t,σit,eit,1=diag22ei1t,22ei1t,22ei1t,σit,eit,2=diag2ei1t,2ei1t,2ei1t,τ = 0.01, t0 = 0 s, ψ = 15, ki(1) = ki(2) = 15, γ = 0.6, hi(1) = hi(2) = 10, and si(1) = si(2) = 12. Through Lemma 2 and (51), we can obtain t1 = 14.427 by simple calculation.The time-varying curves of the system synchronization errorsei(t)(1 ≤ i ≤ 4) are demonstrated in Figure 13 from which one can see that the error converge to zero. Figures 14 and 15 give the identification of unknown parameters α^ii=1,2,3 and β^i=i=1,2,3,4. It is clear that the estimated parameters are approached to values (35, 3, 28)T and (36, 3, 20, 1)T. Figure 16 is the system error values E(t). The conclusion of Figure 16 matches that of Figure 13. The all attest that the control scheme can realize the identification and synchronization of two double-weights Markonian jumping complex networks, effectively.Figure 13
Synchronization errorsei(t) with time evolution.
(a)
(b)
(c)
(d)Figure 14
Identification of the uncertain parametersα^it.
(a)
(b)
(c)Figure 15
Identification of the uncertain parametersβ^it.
(a)
(b)
(c)
(d)Figure 16
Synchronization error valueE(t) with varying time.


## 4.3.2. Simulation Analysis II

In order to further verify that the synchronization criterion can be applied to the synchronization research of double-weights Markovian jumping complex networks, different driving systems are selected for simulation analysis.The Liu system is taken as the node dynamic of the drive complex network. Among them,k = 1 and h = 4:(55)x˙i1tx˙i2tx˙i3t=0−kxi1txi3thxi1t2+xi2t−xi1t000xi1t000−xi3tαi1αi2αi3.The identifying value of the parameter vectorαi is αi=αi1,αi2,αi3T=10,40,2.5T.The four-dimensional hyperchaotic Lüchaotic system will be the node dynamic of the response complex network.Other relevant parameters are the same as those in Section4.3.1. t1 = 17.806 is calculated.The time-varying curves of the system synchronization errorsei(t)(1 ≤ i ≤ 4) are demonstrated in Figure 17. From that, one can see that the error converge to zero. Figure 18 is the identification of unknown parameters α^ii=1,2,3. It is clear that the estimated parameters are approached to the value (10, 40, 2.5)T.Figure 17
Synchronization errorsei(t) with time evolution.
(a)
(b)
(c)
(d)Figure 18
Identification of the uncertain parametersα^it.
(a)
(b)
(c)Under the same research background and system parameters, the control scheme proposed in this paper is compared with the control scheme without finite time to demonstrate the better performance of the proposed control scheme. Figures19 and 20 are the simulation results under the control scheme without finite time. Figure 19 reveals the time evolution curves of synchronization error ei(t). Figure 20 indicates that estimated parameter α^i converge to their truth values (10, 40, 2.5)T. By comparing with the proposed control scheme in Figures 17 and 18, it is obvious that the proposed control scheme can achieve synchronization and identification more effectively and cost less time. Figure 21 is the synchronization error value E(t) with time-varying under the control scheme with finite time and without finite time.Figure 19
Synchronization errorsei(t) with time evolution under the control scheme without finite time.
(a)
(b)
(c)
(d)Figure 20
Identification of the uncertain parameterα^it under the control scheme without finite time.
(a)
(b)
(c)Figure 21
Synchronization error valueE(t) with time-varying under the control scheme with finite time and without finite time.These illustrate that the synchronization and identification of the drive and response networks have been effectively achieved in finite time.


## 5. Conclusions

In this paper, the models of nonidentical dimension Markovian jumping complex networks with stochastic perturbations under the conditions of single weight and double weights are established, firstly. Then, by means of Itô’s formula and finite-time convergence theory, feedback controllers and some update laws are given under corresponding conditions to realize the synchronization, and the value of uncertain parameters in networks can be identified successfully. Finally, the simulation results are shown and the effectiveness and feasibility of the proposed method and the method can be applied to various complex dynamic networks.



---

*Source: 9713652-2020-04-08.xml*