Abstract

This paper is devoted to the study of synchronization problems in uncertain dynamical networks with time-delayed nodes and coupling delays. First, a complex dynamical network model with time-delayed nodes and coupling delays is given. Second, for a complex dynamical network with known or unknown but bounded nonlinear couplings, an adaptive controller is designed, which can ensure that the state of a dynamical network asymptotically synchronizes at the individual node state locally or globally in an arbitrary specified network. Then, the Lyapunov-Krasovskii stability theory is employed to estimate the network coupling parameters. The main results provide sufficient conditions for synchronization under local or global circumstances, respectively. Finally, two typical examples are given, using the M-G system as the nodes of the ring dynamical network and second-order nodes in the dynamical network with time-varying communication delays and switching communication topologies, which illustrate the effectiveness of the proposed controller design methods.

1. Introduction

A complex network refers to a large-scale network with a complex topology and dynamic behavior. It consists of a large set of interconnected nodes, in which the node is a fundamental unit with specific contents [1]. Complex networks exist in many fields of science, engineering, and society and have attracted much attention in recent years [25]. As one of the most important collective behaviors, synchronization phenomena have been a topic of research. There are a lot of synchronization phenomena in physical sciences and in mathematics. Several books and reviews [69] have also dealt with this topic. Such applications are pervasive and include clock synchronization in complex networks [1012], coordination of unmanned aerial vehicles [13], and fair allocation of network resources [14].

In the past decade, special attention has been paid to the synchronization of dynamical systems, particularly large-scale and complex networks of chaotic oscillators [15, 16]. With the discovery of the small-world effect and scale-free feature of most complex networks, a lot of works have been concentrated on synchronization in small-world networks and scale-free networks [1720]. In these investigations, an essential requirement is that the structure of the networks and the coupling functions are known beforehand. In order to overcome the aforementioned constraints, synchronization in complex networks by controller methods has also been investigated [2124]. Zhou et al. [21] proposed some synchronization criteria and designed simple controllers for several specific complex networks. In [22], for a class of complex networks with uncertain inner couplings, the event-triggered synchronization control problem was put forward. The uncertain inner coupling under consideration was characterized in terms of the interval matrix. In order to save the communication and computation resources, the event-based mechanism was adopted and an event-triggered synchronization control scheme was proposed for the complex networks. Based on the leader–follower model, Zhou et al. [23] proposed an improved network structure model for realizing the cluster synchronization on multiple subnetworks of complex networks, and some suitable pinning controllers on the chosen nodes of each follower’s subnetwork are designed. Yu et al. [24] employed a decomposition approach to incorporate the nodes’ inertial effects into the distributed control design for second-order nodes in a dynamical network with communication delays and switching communication topologies.

However, the following problems still exist in the control of dynamical networks and have not been solved in a good way. It is often difficult to get an exact estimation of coupling coefficients in reality. Therefore, it is preferable if the coupling functions in a dynamical network are not restricted to be completely known. If the network structure is only partially known or completely unknown, it is very difficult or even impossible to design a controller to achieve the intended network synchronization. Besides the problem above, the nodes in the dynamical network usually have the characteristics of time delay. Because of the finite speeds of transmission as well as traffic congestion, a signal or influence traveling through a complex network is often associated with time delays. Time delay is the main factor which gives rise to instability in dynamical networks. This issue has attracted the attention of many scholars recently. In [25, 26], the authors found that the solution space of time-delayed networks was infinite-dimensional, and its theoretical analysis was often very difficult. Dhamala et al. [27] assumed that the coupling was linear between nodes and explored the effect on the synchronization of time-delayed nodes in various networks. Then, Atay et al. [28] found that altering the time delay would be helpful to improve the capability of the network synchronization. Gao et al. [29] studied the uncertain time-delay systems and robust stability analysis. They solved the problems about a class of nonlinear time-delay systems with memory stabilization. In [30], a complex dynamical network model with time-varying inner coupling was introduced. And a sufficient condition was given to achieve the exponential synchronization.

In this paper, we will give positive answers to all the questions above. We combine the adaptive control method with the network with time-delayed nodes and coupling delays. The local and global features of adaptive synchronization of networks are studied based on the Lyapunov-Krasovskii stability theorem. The results we got are suitable for the general dynamical network; that is, they are independent of the network topology. The paper is organized as follows. The general network model and several mathematical preliminaries are described in Section 2. In Section 3, the analysis of adaptive synchronization is studied, and also two desirable controllers are designed. Section 4 gives two typical examples. The M-G system is used with the time-delayed nodes in the ring network and second-order nodes in the dynamical network with time-varying communication delays and switching communication topologies. Simulation results show the effectiveness and correctness of the proposed methods. Finally, conclusions are given in Section 5.

2. Network Model and Mathematical Preliminaries

2.1. A Complex Dynamical Network Model with Time-Delayed Nodes and Coupling Delays

Consider a nonlinear complex network consisting of identical delay nodes, which are coupled as follows [31]:where represents the state vector of the th node; represents a smooth nonlinear vector function; represents a smooth unknown nonlinear coupling function; , represents the control input.

When the network reaches synchronization, that is, , , the coupling control will disappear at this point; that is,This ensures that the solution of any solitary node is also the solution of the synchronous coupling network.

As can be seen from the Introduction, there is a delay in the actual network. After adding the time-delay term in the model of (1), the dynamical network with these time-delayed nodes and coupling delays can be written as where represents the state vector of the th node; , represents a continuous nonlinear function; , represents an unknown continuous nonlinear coupling function; is a time delay; and , represents a control input.

Obviously, when the complex network achieves synchronization, there is , and the coupling and control items will disappear at this time; that is, .

2.2. Mathematical Preliminaries

The following definitions and assumptions are necessary for discussing the network synchronization problem.

Definition 1. Let be a unique solution of a complex network (3) which satisfies the following equation: where can be an equilibrium point, a nontrivial periodic orbit, or even a chaotic orbit.

Definition 2. If there is a nonempty subset , for any , there is , which satisfies the following formula:Then, the solution of the complex network (3) is asymptotically synchronized with (5). is called the synchronous region.
Next, subtracting (2) from (4) yields the following error dynamical system: whereBy linearizing (6) which is evaluated at , we can get the following equation: where , , is a Jacobian matrix of on , is a Jacobian matrix of on , and there is .
In order to facilitate the follow-up study, we present the following assumptions.

Assumption 3. Suppose that and are bounded. That is, there exist nonnegative constants and satisfying and .

Assumption 4. Suppose that the nonlinear function is Lipschitz continuous. That is, , , the following formula holds:Therefore, (3) can be written as follows:where is a nonnegative Lipschitz constant, .

Lemma 5 (see [32]). For any , one has .

3. Synchronization Analysis and Design of Adaptive Controller

In this section, we apply the Lyapunov-Krasovskii stability theorem to the problems of local and global synchronization of coupled delay dynamical network with time-delay nodes. Then, we design an adaptive controller according to the different situations.

3.1. Local Synchronization of Coupled Delay Dynamical Network with Time-Delay Nodes

Consider all the nodes in the dynamical network without loss of energy and generality. Assume that the ith  () node is controlled by the following adaptive controller, which is described bywhere and ,  , are positive constants. Thus, the controlled network (2) can be rewritten as follows:When the parameters are known or unknown but bounded nonlinear couplings, there exists the following theorem.

Theorem 6. If Assumptions 3 and 4 are established, then the complex network (2) will achieve local asymptotical stability under the action of the adaptive controller (11).

Proof. Construct a Lyapunov candidate as follows:where , and , are positive constants, and is the positive definite matrix that needs to be determined. Thus, the time derivative of along (6) isSubstituting the adaptive controller (11) into (14) yields the following: According to Assumptions 3 and 4 and Lemma 5, selecting , where the constant , it can be obtained thatwhere , .
Therefore, we can choose the appropriate constants , and , such that the diagonal matrix is a negative matrix. This equation shows that , is uniformly bounded, and the error system (8) is asymptotically stable by the action of the adaptive controller (11). By using the Lyapunov method, it can be clearly seen that the error system (6) is also asymptotically stable under the action of the adaptive controller (11). It is concluded that is an asymptotically stable equilibrium point of the error system (8) and it is also deduced that is an asymptotically stable equilibrium point of the system.

The above stability criteria are delay-independent. In a coupled complex network with different delays, (17) can also be synchronized using an adaptive controller (11): where , is the coupling delay; the proof is similar to that of Theorem 6. We can derive from Theorem 6 that synchronization of such complex networks mainly depends on three basic parameters: the dynamic characteristics of complex network nodes , the dynamic parameters of the adaptive controller , and the network structure parameters .

In summary, we can see that the dynamical network with time-delay nodes and coupling delay structure can achieve local synchronization; that is, complex network synchronization can be realized under the action of adaptive controller (11).

3.2. Global Synchronization of Coupled Delay Dynamical Network with Time-Delay Nodes

In this section, we discuss the global synchronization problem with time-delay nodes and coupled delay networks. The network structure is as follows:where , is an unknown and smooth function. According to the derivation process of Section 2, we can compute the error dynamic equation of (11):where , .

The definitions of other sections are the same as those defined in Section 2.

Assumption 7 (see [33]). Suppose that there are unknown but nonnegative constants , satisfyingThe following theorem gives a method of global synchronization criterion for the dynamical network (18).

Theorem 8. If Assumptions 3, 4, and 7 are established, then the complex dynamical network (18) achieves global asymptotic stability under the action of the adaptive controller of (11).

Proof. Select the following Lyapunov function:where , and , are positive constants, and is the positive definite matrix that needs to be determined. Thus, the time derivative of along (6) is The substitution of the adaptive controller (3) into (13) yields the following:According to Assumptions 4 and 7 and Lemma 5, selecting , where the constant , the following can be obtained:where , .
Therefore, we can choose the appropriate constants , and , such that the diagonal matrix is a negative matrix. This equation shows that , is uniformly bounded. By using the Lyapunov method, the error system (18) is also asymptotically stable under the action of the adaptive controller (11) (e.g., ). It is concluded that is an asymptotically stable equilibrium point of the error system (19).

4. Numerical Simulation

In this section, two typical examples of the complex network are used to show the correctness and effectiveness of the adaptive controller (11). One example is a ring network with fifty nodes, in which the dynamic performance for each node is complex. The other example is a dynamical network with switching topologies and twelve second-order nodes, in which the dynamic performance for each node is more complex. And the network state equation satisfies (3).

4.1. Simulation I

Chaos often comes from the bifurcation control to the equation of the formwhere is a nonlinear function and , . Generally speaking, maximum production occurs at some intermediate value of , and thus the control displays both positive and negative feedback characteristics. The M-G system [34, 35]comes from the equation in case of controlling the density of circulating blood cells. We choose this model for our studies since it has been thoroughly investigated in the literature [34] and is easily implementable electronically [35]. The number of positive Lyapunov exponents as well as the dimension of the strange attractor of this system can easily be controlled by varying the delay time . For fixed values of the parameters, both of these quantities increase linearly with the increase of [34]. In order to analyze the system dynamics and the proposed controller design, it is adopted in (25) that a polynomial approximation of the M-G system shows the same qualitative behavior. The system equation thus becomesWhen , and , system (27) becomes chaotic. This system has three unstable fixed points:We choose a ring network with fifty nodes as our complex network, because the network is one regular network of all typical networks and is an ordinary mathematical model [36], which can represent many physical models. Assume that it satisfies the following conditions:

() The complex network coupling function form is

() It satisfies the previous conditions in Assumptions 3 and 4.

() The dynamic equations of each node are M-G systems, where , and the coupling strengths are all 0.1.

The network structure is shown in Figure 1.

Then, the state equation of this complex network can be written as From the knowledge of Section 2, it is well known that such networks are very hard to synchronize. Figure 2 shows the state responses of this ring complex network with and without control. It is obvious that the network is desynchronized.

The synchronization target node is selected as . The initial states of 50 nodes are selected randomly from 0 to 25. Let , , and then according to Theorem 8, it can be seen that this complex network is locally asymptotically synchronous under the action of adaptive controller (11). The adaptive synchronization analysis of complex networks is carried out by using MATLAB to prove that the adaptive controller can synchronize the complex network on the trajectory of the specified chaotic node. If the simulation result is correct, this will prove that the controller is effective. Figures 3 and 4 show the synchronous error , , for and with the same controller, respectively.

Obviously, it should be pointed out that the controlled network is asymptotically stable at zero under the effective adaptive controller (11). The simulation results indicate that the coupling delays create some obstacle to the synchronization of the complex networks. The response performance is deteriorated as the time delays become large. It is clear from the figures that the errors are very small.

4.2. Simulation II

The proposed theorems have been used for synchronizing the second-order nodes in the dynamical complex network. Consider twelve second-order nodes in the dynamical network with switching topologies. The state equation of each second-order node is described bythat is, where are the position vectors; stand for the velocity vectors; are symmetric positive definite matrices; and are the control inputs. Let , , , and . According to (32), for the arbitrary nth node, there is

Suppose that there exists a communication time delay between nodes and . The interaction topology of a dynamical network of nodes with second-order nodes is represented using a directed graph with the set of nodes , the set of directed edges , and the adjacency matrix . is the set of node subscripts, and denotes the directed edge from node to node . The elements in the adjacency matrix satisfy , (if and only if ). The set of neighbors of a node is . Let ; the degree matrix of a directed graph can be expressed as . The Laplacian matrix is defined by . It can be shown that, using the Gershgorin disc theorem [14], all of the eigenvalues of have a nonnegative real part. Furthermore, if is undirected, then the Laplacian matrix of is symmetric and satisfies , which means that is positive semidefinite. The in-degree and out-degree of node can be defined as and , respectively. If , then node is an equilibrium point. If there exists a direct graph between nodes in graph G, then is a strongly connected and balanced graph.

Figure 5 gives three strongly connected and balanced graphs with 0–2 weights. The initial values of and are selected randomly in the region and and , respectively.

The inertias of the twelve second-order nodes are , . We set ; Figures 6 and 7 show the numerical simulation results of the twelve second-order nodes in the dynamical network with an arbitrary switching signal. Figure 6 gives the inertial nodes’ position error curves. Figure 7 shows the inertial nodes’ velocity error curves. It is obvious that, for the dynamical network with communication time delay and switching topology, the control strategy proposed in this paper also achieves the exponential stability.

From the above, the synchronization process of the dynamical network can be described as follows: detect the designated node and the other nodes; then, the state variables of nodes are subtracted from the designed node, respectively, by using the subtraction device. This way, we can get the error signals , and then the error signals are sent to the adaptive controller. Furthermore, the adaptive controller chooses a one-dimensional variable as the control variable. According to , we can adjust the control outputs and make the variables of the designated node and other nodes synchronized. At the moment, the error is close to zero. A block diagram is expressed in Figure 8.

5. Conclusions

This paper studied the adaptive synchronization of uncertain dynamical networks with time-delayed nodes and coupling delays, with known or unknown but bounded nonlinear coupling functions. The parameters of the adaptive controllers can be chosen according to the presented theorems which can stabilize the network to the desired fixed points. The adaptive synchronization criteria by using the Lyapunov-Krasovskii stability theory also provide some new insights into the network synchronization and the possible applications in engineering systems. At last, an example of a dynamical network with fifty M-G system nodes has been simulated, demonstrating the effectiveness of the proposed controllers. Note that the adaptive controllers proposed are both delay-dependent, so there are no compensations to the delay couplings. However, how to solve the delay-dependent problem to reach synchronization seems to be new and useful. Furthermore, as another issue, in some real complex systems, interacting subjects are not identical and have different individual properties. Some synchronization criteria are expected to be more powerful and deserve further exploration.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The work of Miao Yu was supported by the National Natural Science Foundation of China (51407201, 61402032) and Doctoral Foundation of Beijing University of Civil Engineering and Architecture under Grant 00331614032.