Abstract
The state observer for dynamic links in complex dynamical networks (CDNs) is investigated by using the adaptive method whether the networks are undirected or directed. In this paper, a complete network model is proposed, which is composed of two coupled subsystems called nodes subsystem and links subsystem, respectively. Especially, for the links subsystem, associated with some assumptions, the state observer with parameter adaptive law is designed. Compared to the existing results about the state observer design of CDNs, the advantage of this method is that a estimation problem of dynamic links is solved in directed networks for the first time. Finally, the results obtained in this paper are demonstrated by performing a numerical example.
1. Introduction
In recent past decades, the research on CDNs has become a hot topic in many fields [1–4]. From the perspective of large system, a complete CDN contains many nodes and links (weights of connections between nodes), which implies that a complete CDN is composed of the nodes subsystem and links subsystem, and the two subsystems are usually coupled with each other [5–7]. It is worth noting that the existing researches mainly focus on the nodes subsystem because some behaviors are reflected by nodes such as synchronization [8, 9], stabilization [10, 11], and consensus [12, 13].
From the above results about the synchronization, stabilization, consensus, or other problems of CDNs, it is easy to see all states in CDNs, including the states of nodes and links, are required to be measured accurately. However, this assumption is too hard to be satisfied in practice because of the influence of external environment, measurement costs, and technical constraints [14]. Thus, constructing state observers for the CDNs to estimate the unknown states is very necessary and important. Fortunately, some scholars have discussed the state estimation problems of CDNs and obtained some research results, including cases with the coupling time delays [15, 16], packet loss [17, 18], stochastic noisy disturbance [19], and uncertain coupling strength [20].
However, the above results only consider the estimation problems of the states in nodes subsystem, and assume that the links between nodes are known. It implies that the measurement and state estimation problems of links in the CDNs are ignored. In fact, due to the limitation of measurement methods, the state values of links in CDNs are more difficult to be measured accurately in practical situation, compared to the states of nodes. Hence, only a few papers have studied and discussed the effective measurement problem of the links between individuals (nodes), and the measurement method mainly depends on the physical interaction between individuals [21] or the adaptive weights of links [22]. Similar to the state values of nodes, not all state values of links’ weights can be measured and obtained. Therefore, it is necessary to design observers to estimate the unmeasured state values of links. As we know, there is only one paper to have discussed the state estimation problem of dynamic links in CDNs [23]. Unfortunately, the method proposed in [23] is only effective for undirected networks and cannot solve the estimation problem of dynamic links in directed networks.
Inspired by the above discussions, this paper mainly focuses on the state observer design for dynamic links in directed networks. Specifically, a mathematical model for a class of directed CDNs is proposed, which is described by both the nodes subsystem and links subsystem with coupling between the two subsystems, and we have designed a state observer for the links subsystem by using the adaptive method. This means that a state estimation problem of dynamic links in directed networks is solved for the first time, which is also regarded as the biggest contribution of this paper.
The rest of this paper is organized as follows: in Section 2, a complete CDN model is proposed, which is composed of the nodes subsystem and links subsystem with outputs; Section 3 introduces the design process of state observer for the links subsystem; in Section 4, the simulation example is presented and used to demonstrate the effectiveness of this method; finally, the conclusions are given in Section 5.
1.1. Notations
The -dimensional Euclidean space is denoted as , the set of real matrices is denoted as , the Euclidean norm of a vector or a matrix is denoted as , and the transpose of matrix and -dimensional identity matrix is denoted as and , respectively.
2. Preliminaries and Model Description
If the states of nodes and links in CDNs evolve over time, then the mathematical model of CDNs, including directed and undirected networks, can be described by both the nodes subsystem and links subsystem, where the two subsystems are coupled with each other. In this paper, we only consider the case that each node is -dimensional continuous system in CDNs with nodes, then the nodes subsystem and links subsystem can be described by vector differential equations and matrix differential equation as follows, respectively:where is the state vector of node ; the constant matrices and ; the vector functions and ; is a known constant, which denotes the common connection strength of node in the network; the constant matrices and ; the coupling matrix , and , where is a bounded and closed set in ; the output matrix is given; and the links matrix , where the state variable denotes the weight of link from node to node . Especially, for undirected networks, and at least, one pair such that for directed networks. In addition, if , then denotes the link strength of node itself.
For the CDNs composed of subsystems (1) and (2), the following instructions are given:(1)The background of dynamic links is given as follows. For example, the biological neural networks consist of neurons (nodes) and synapse (links), and Gamma oscillations in neurons may cause the synaptic facilitation, which is regarded as a dynamic behavior of the links [5, 6, 24]. Similarly, the web winding systems can be regarded to be composed of motors (nodes) and the web (links), and the regulation values of web tensions vary with the speed of the motors, which is also regarded as a dynamic behavior of links [25]. In the above examples, the state values of links need to be measured by some sensors.(2)The existing research results show that the nodes in networks can emerge synchronization or stabilization phenomenon with the help of the links, which mean that the nodes are the main body of synchronization and stabilization [8–11]. In contrast, the links as another part of networks can also emerge some characteristic phenomena in many real networks, such as the structural balance in social networks [5, 6, 26]. It is worth noting that the paper [26] has researched on structural balance by using the Riccati matrix differential equation, and the reason is that this type of equation is more easily to emerge the phenomenon of structural balance. In view of this, we choose linear Riccati matrix differential equation to describe the links subsystem. Clearly, the model of CDNs, composed of both nodes subsystem (1) and links subsystem (2), can help us to understand and explain the dynamic behaviors of networks in a better way.(3)The subsystem (2) is used to describe dynamic change of links’ weights in the CDNs, and in general, the CDNs are directed. However, if and , then we can obtain , while the initial value of the state in subsystem (2) satisfies . Hence, the subsystems (1) and (2) can be used to describe both undirected and directed networks (the undirected networks can be regarded as a special case of directed networks). To the best of my knowledge, there is only one paper to have solved the state estimation problem of links subsystem [23]. However, this method is only effective for undirected networks, but not for directed networks. This drives us to study estimation problems of dynamic links in directed networks.(4)It is difficult to accurately measure all states of the links between individuals (nodes) in practical applications, which imply that only partial states in (2) can be measured accurately and made available (). It is worth noting that the precise measurement of the partial states is bidirectional; that is, if is measurable, then must also be measurable. That is why the two outputs and appear in (2).
Now, some useful definitions and operators involved in this paper will be introduced as follows.
Definition 1. (see [27]). The application is defined bywhere the matrix is called the vectorization operator.
Definition 2. (see [27]). If there are two matrices and , then the Kronecker product of and is denoted as and defined as follows:By using Definitions 1 and 2, the following basic properties about Kronecker product and operator can be obtained and shown as follows [27]:(1)(2)(3)(4)(5) and are the matrices with compatible dimensions, and represents the identity matrix with compatible dimensions. Especially, it is assumed that both and are invertible in property (3).
According to Definitions 1 and 2 and their corresponding properties, the Riccati differential equation (2) can be rewritten aswhere and .
Assumption 1. For the links subsystem (2), the double matrices and are completely stable.
If Assumption 1 is true, then we can obtain matrices and , which can make and to be Hurwitz stable, respectively. Thus, as long as any matrices and are given, there must be positive definite matrices and that satisfy the following two Lyapunov equations, respectively:
Lemma 1. If Assumption 1 is true, then the following Lyapunov equationshold, where , , , , , and . Clearly, , , and .
Proof. If Assumption 1 holds, then the following equations can be obtain from (6) and (7):Using the properties of Kronecker product, (9) and (10) can be rewritten asThus, we can get If we multiply both sides of the equalities (12) and (13) by and from right, respectively, then we get thatIt is noticed that . Therefore, the equalities (14) and (15) can be rewritten as follows:Thus, Lemma 1 is completely proved.
Assumption 2. For subsystem (2), in which the coupling matrix satisfies that , where .
If Assumption 2 holds, then we can get that ≤ . Meanwhile, we note that is a bounded and closed set in , and . Thus, there exists a positive constant to satisfy the inequality .
General speaking, is unknown. However, we can use the adaptive method to estimate it. In this paper, we use to denote the estimated value of . Hence, the estimation error is denoted as .
3. Main Results
Definition 3. Designing a matrix differential system , if the state satisfies , then the matrix differential system can be regarded as a state observer of the links subsystem (2).
If Assumptions 1 and 2 hold, the state observer of the links subsystem (2) can be designed and presented as follows:with the following adaptive lawwhere denotes the estimated value of the state in (2); the robust term , where , is a given positive constant, and the matrices , , , and can be obtained by solving the Lyapunov equations (6) and (7), respectively.
According to (3) and (4), we can deduce from (17) thatClearly, ; thus, we get .
In this paper, the estimation error is denoted by . By using (3), (4), and properties about Kronecker product and operator, we can get the following error system:
Theorem 1. If Assumptions 1 and 2 are true, then the matrix differential system (17) with the parameter adaptive law (18) is the state observer of the links subsystem (2).
Proof. Consider the following Lyapunov function:Calculating the orbit derivative of along (20) gives thatFrom inequality (22), we can obtain that the estimation error matrix is bounded and . Thus, Theorem 1 is completely proved.
4. Simulation Example
In this paper, we consider a continuous analog Hopfield network with 10 neurons () [23, 28], which is composed of nodes subsystem and links subsystem, where the nodes subsystem is described as follows:where , , , and .
Meanwhile, we assume that the changes in the links’ weights satisfy the Riccati differential equation (2). If we choose and and randomly select matrices , and satisfying Assumption 1, then the matrices and can be obtained by solving the Lyapunov equations (6) and (7), respectively. Thus, we can get the coupling matrix in (2) satisfying Assumption 2.
Finally, randomly select the initial values of states , and , in the range , and the numerical results are shown in Figures 1–5:(i)From Figures 2–4, we can see that the estimation error converges asymptotically to zero. According to Definition 3, we know that the Riccati dynamical equation (17) with the adaptive law (18) is a state observer of the subsystem (2), and the state observer is effective.(ii)Compared to the results in [23], our advantage is that the result about the state observer of the subsystem (2) is true whatever the network is directed or undirected. Meanwhile, it is worth noting that, due to the effect of the parameter adaptive law (18), the state observer (17) does not contain the states of the nodes. This shows that the state observer is less affected by the dynamic changes in the nodes and thus improves the robustness of the state observer.
5. Conclusions
In this paper, a complete model of CDNs is proposed, which is composed of two coupled subsystems, called nodes subsystem and links subsystem, respectively. Contrary to the existing results on the state estimation problem of nodes subsystem, we mainly focus on the state estimation of the links subsystem with outputs and have designed a state observer with the parameter adaptive law to estimate the state of the links subsystem in this paper. In particular, this method solves the estimation problem of dynamic links in directed networks for the first time. Meanwhile, it implies that we can use the state estimation information of the links to directly design a controller for the links subsystem; thus, some control problems may be solved effectively. Therefore, the design method of state observer for dynamic links proposed in this paper can enrich the achievements about the state estimation of CDNs.
Data Availability
In this paper, we submitted data mainly related to theoretical proof and numerical simulation, in which the part of numerical simulation is realized by Matlab software; if necessary, we can provide simulation source program and relevant data at any time.
Conflicts of Interest
The authors declare that there are no conflicts of interest.
Acknowledgments
This work was supported by the Key Laboratory of Chongqing Municipal Institutions of Higher Education ([2017]3), the Program of Chongqing Development and Reform Commission (2017[1007]), the Scientific and Technological Research Program of Chongqing Municipal Education Commission (KJ1710244, KJQN201801215, KJQN201801209, KJQN201901236, and KJ1710241), and the Scientific Research Project of Chongqing Three Gorges University (19QN07).