Research Article | Open Access
Lifu Wang, Yali Zhang, Jingxiao Han, Zhi Kong, "Quantitative Controllability Index of Complex Networks", Advances in Mathematical Physics, vol. 2018, Article ID 2586536, 9 pages, 2018. https://doi.org/10.1155/2018/2586536
Quantitative Controllability Index of Complex Networks
In this paper, the controllability issue of complex network is discussed. A new quantitative index using knowledge of control centrality and condition number is constructed to measure the controllability of given networks. For complex networks with different controllable subspace dimensions, their controllability is mainly determined by the control centrality factor. For the complex networks that have the equal controllable subspace dimension, their different controllability is mostly determined by the condition number of subnetworks’ controllability matrix. Then the effect of this index is analyzed based on simulations on various types of network topologies, such as ER random network, WS small-world network, and BA scale-free network. The results show that the presented index could reflect the holistic controllability of complex networks. Such an endeavour could help us better understand the relationship between controllability and network topology.
In recent years, the study of complex networks has drawn the attention of many scholars from both the science and the engineering communities [1–5]. Studies of theirs impact our understanding and control of a wide range of systems, from the Internet and the power-grid to cellular and ecological networks. While there are many challenges in the process of research due to the diversity of complex networks, one most challenging issue in modern science and engineering networks is the control of complex networks [6, 7]. The controllability of a dynamical system reflects the ability of external input information to influence the motion of the overall system. A complex network is controllable if, imposing appropriate external signals, the system can be driven from any initial state to any final state in finite time [8–10]. Although great effort has been devoted to understanding the interplay between complex networks and dynamical processes taking place on them in various natural and technological systems [11–17], the control of complex dynamical networks remains to be an outstanding problem. The Kalman controllability formed the foundation of the controllability theory by a set of algebraic criteria to check whether or not a given system is controllable. However, the Kalman controllability is qualitative; we cannot know how it is difficult or easy to control, that is, cannot know the size of network controllability for one network system. Therefore, it is necessary to find a method that can quantitatively measure controllability of a given system.
With deeper understanding of the controllability of complex networks [18–26], more scholars have begun to explore the ability to control complex networks through various indicators. For example, Liu et al.  addressed the structural controllability of arbitrary complex directed networks, identifying a minimal set of driver nodes that can guide the system to any desired state. They selected an index denoted by to quantitatively measure the extent of controllability of complex networks. Liu et al.  propose the concept of node control centrality to quantify the ability of a single node of complex network and calculate the distribution of control centrality for some networks. Jia et al.  proposed the concept of node control capability in the further study; that is, the possibility of node becoming a driving node was calculated. Many research works have been published under the theoretical framework of structural controllability [21–26]. However, these quantitative controllability methods only roughly measure the dimensionality of controllable subsystems; they could not better distinguish the controllability of systems with the same number of dimensions of controllable subsystems. Cai Ning  studied a way to quantitatively measure the extent of controllability of any given controllable network. This method is based on fully controllable networks; for uncontrollable networks, it could not make a judgement.
In the current paper, we will try to explore the possible ways to quantitatively measure the extent of controllability of any given network. An index will be proposed to assess the controllability of a dynamical network, that is, the network whether or not being easily or difficulty controlled via the input information. While the previous Kalman controllability index is qualitative, only defining whether or not a system is controllable, the one in our paper is quantitative. Particularly, for a not completely controllable system, we can quantitative measure controllability of the controllable subsystem by decomposing the network into a controllable subsystem and an uncontrollable subsystem, so that we can measure how far the network is from being uncontrollable. Such a route may be called quantitative index measure controllability of complex network. Then simulations will be performed to show the effect of controllability index between different network topologies on three distinct types of model networks, namely, the random networks, the small-world networks, and the BA scale-free networks. Through simulation results, it can be found that controllability of given network is related to their topologies, such as the number of nodes and edges and edge density, so it is possible to improve the controllability of the network by adjusting certain parameters, such as the connectivity probability and the number of nodes .
2. Controllability of Complex Networks
Consider a complex system described by a directed weighted network of nodes whose time evolution follows the linear time invariant dynamics .where is the state vector of each node at time . denotes the state matrix that depicts the linking strength between the state nodes. The matrix element gives the strength or weight that node can affect node . is the input matrix, indicating the nodes that are controlled by the time-dependent input vector with independent signals imposed by an outside controller. The element represents the control strength from the control node to the state node . A system is controllable if we can drive it from any initial state to any desired final state in finite time; otherwise, it is uncontrollable.
The rank of the controllability matrix C, denoted by rank(C), provides the dimension of the controllable subspace of system (). The system is not completely controllable if . Then the system can be decomposed into a controllable subsystem and an uncontrollable subsystem through a linear transformation , which iswhere is the controllable subsystem whose dimension is , and is the uncontrollable subsystem whose dimension is .
When we control node only ( is the ith node in the network), reduces to the vector with a single nonzero entry, and we denote with . We can therefore use rank as a natural measure of node ’s ability to control the system: if , then node alone can control the whole system; that is, it can drive the system to travel between any points in the N-dimensional state space in finite time. Any value of less than provides the dimension of the subspace that can control. In particular, if , then node can only control itself.
Kalman controllability theory provides great convenience for checking whether or not a given network is controllable. However, from the point of view of Kalman controllability, there are some problems that restrict the study of the controllability of complex networks to some extent. First, almost any arbitrary system is completely controllable in the sense of Kalman controllability . This fact reduces the significance of Kalman controllability theory. In addition, Kalman controllability of dynamical system is qualitative, which means that, in Kalman's theory, the system could be only defined as either controllable or uncontrollable. However, for many real-world network systems, people need to know more information about the controllability issue such as how much our control force should be to control the whole system or maybe part of the system. In the next section, one index is brought up with to measure the controllability of networks quantitatively.
3. Quantitative Index Measuring Controllability
Before we begin to analyze quantitative controllability, we first introduce the definition of conditional number.
Definition 1 (see ). The condition number of a square matrix iswhere denotes 2-norm of the square matrix .
The condition number can measure the nonsingularity of a matrix, with range of variation in . When is identity matrix, . On the contrary, when is nearly singular matrix, . That is, the greater the condition number is, the closer the matrix is to being singular.
3.1. Examples of Analysis Controllability
Let us start the analysis on quantitative index measuring controllability by investigating two typical examples.
Consider two simple networks of nodes with different structures. The controlled network has an input connecting to a state node . The network topology is shown in Figure 1. The controllability of node to the whole network is analyzed as in Figure 1.
For Figure 1(a), the state matrix of the network isWhen we control node only, the input matrix isThen the controllability matrix by controlling node isTherefore, the rank of controllability matrix , indicating that the whole system is not completely controllable by controlling node only. Yet, there exists a controllable subsystem with 6 nodes, and these nodes can be found by all stem-cycle disjoints . For example, constitute a controllable subsystem. In other words, by controlling node with a time-dependent signal , we can drive the subsystem from any initial state to any final state in finite time. According to linear system theory, the network system () can be decomposed by linear transformation [31, 32]. We can get the controllable subsystem () as follows:Then the controllability matrix of the controllable subsystem () isThen subsystem () can be completely controllable after the structural decomposing. We can compute that , indicating that the distance from the controllable subsystem to the uncontrollable ( singular) is .
Similarly, for Figure 1(b), we can get the state matrix of the network:And when the node is controlled only, the input matrix isThen the controllability matrix isSo we also can get when we control node only. Through structural decomposing, the controllable subsystem () isAnd the controllability matrix of the controllable subsystem isThe controllable subsystem is also . For the controllable subspace, we can get that , indicating that the distance from the controllable subsystem to the uncontrollable one ( singular) is . By comparing Figures 1(a) and 1(b), both controllable subsystems with the same size, we can see that they do not have the same distances from being uncontrollable. In these two examples, the two controllable subsystems both are equally distant from being uncontrollable. However, there are some networks that are controllable but have a closer distance to being uncontrollable; that is, the value of their is very large and is really close to being singular.
3.2. Introduction of New Index: Quantitative Measuring Controllability
In this subsection, a more precise index of evaluating the controllability is presented. More than just saying “controllable” or “uncontrollable” to a network system, it can quantify the controllability of the given network. Such a route may be called quantitative controllability index.
Definition 2 (quantitative controllability index). The quantitative index measuring controllability of single node can be computed bywhere is the rank of controllability matrix by controlling node only; ; thus provides the dimension of the subspace node can control; ; is the controllability matrix of controllable subsystem by controlling node only; then can measure the distance from the controllable subsystem to the uncontrollable one. Thus, captures the ability of node to control the whole network. In addition, can continuously characterize the controllability of the given network, while as a positive integer can only measure the dimension of the controllable subspace (the higher is, the stronger controllability given network achieves). For controllable subspaces with the same dimension, their controllability size is the same when is used to measure controllability (e.g., Figure 1), so this characterization is not perfect. Therefore, we measure the controllability with the index which adds to the base of and uses , which indicates how far the subsystem is away from being uncontrollable, to measure the controllability of subsytems with the same dimensional controllable subspace. When is the same, plays a key role in measuring the controllability. The larger is, the further the subsystem is away from the fully uncontrollable point, which means that the controllability of this subsystem is stronger. Therefore, the greater value of is, the stronger controllability the network achieves. In Figure 1(a), when the node is controlled separately, we can get . In Figure 1(b), similarly, when the node is controlled separately, we can get . Their controllable subspace is the same, but their distances from the fully uncontrollable point are different. Therefore, the network controllability of Figure 1(a) is stronger than the network controllability of Figure 1(b) by controlling node only.
By the above definition, we can calculate the controllability of each node for the given network. Then we can get the average quantitative controllability index of the whole network:Therefore, can be used to quantitatively measure the controllability of the whole network system. For Figures 1(a) and 1(b), the quantitative measuring the controllability of each node is shown in Table 1. Then, the average quantitative index measuring controllability of whole network is calculated: and as shown in Table 1. Thus, we can conclude that the average controllability of network in Figure 1(a) is slightly weaker than network in Figure 1(b); that is, the network shown in Figure 1(b) is easier to be controlled than the other one in Figure 1(a).
The change of number of nodes does influence the performance, while the increase of edges has no influence on the performance. We have also analyzed the index, especially calculating its time complexity. We estimate the time complexity in (16). For any complex network with nodes, when we use to calculate the controllability of the network, the time complexity of the algorithm is .
The controllability index consists of two parts: one is the control centrality, and the distribution of the control centrality is mainly determined by the degree distribution ; and the other part is taking condition number as a parameter with the range between 0 and 1, which is the indicator of controllability difference among networks with the same number of dimensions, and it acts more like a modifier in , having no decisive effect on distribution of . Therefore, we deduce that the distribution of is mostly determined by the network’s degree distribution. We use averaging of all nodes to assess the controllability of the whole network and to quantify the network’s controllability entirely, just like using average degree to analyze the whole network.
4. Empirical Analysis of Quantitative Controllability Index
The controllability of complex networks is related to their network topology. Therefore, it is significant to find the relationship between quantitative index measuring controllability and the structure parameters of complex network. In this section, we will conduct a series of simulation analyses for three different typical complex network models to observe the regularity with which the network controllability changes with the variations of parameters.
4.1. ER Network
Liu et al.  proved that the minimum number of nodes that the input signals are injected into, called driver nodes, could be determined by detecting the “maximum matching” in the network. A maximum matching is the maximum set of links that do not share start or end nodes. A node is said to be matched if a link of the maximum matching points at it; otherwise, it is unmatched. If the nodes in a network are all matched, then the network is perfectly matched and the number of driver nodes is 1; otherwise, , where is the number of matched nodes. The ratio of driving nodes to total number of nodes can be used to measure controllability of the network. Here, we compare the quantitative index Q with that of the number of driver nodes as shown in Figure 2.
Consider the ER network with twenty nodes, ; each time one of the nodes is controlled separately. Let the connectivity probability p vary across . Here we show the variation of controllability index Q with connection probability in Figure 2. Every data in the figure has been calculated 1000 times and we took the average value.
We can see from Figure 2 that, for small connecting probability , controllability index increases with and reaches the maximum at around . Meanwhile, for higher link density , Q begins to decrease with connection probability p. However, according to the red line in Figure 2, when , the number of driver nodes decreases and the connection possibility p increases. The number of drive nodes is larger for smaller connection probability , indicating the weak controllability of the network because more nodes need to be controlled. When , the number of driver nodes basically remains the same when the connection possibility increases. And when , the number of driver nodes increases with the increase of the connection possibility p, which indicates that the controllability is weaker in the meantime. Besides, the closer the network is to the complete graph, the more nodes need to be controlled and then the lower the controllability is. From the comparison of Q and , it is not hard to find out that, for the same ER networks, the controllability variation trend is basically the same. Index Q is positively correlated to the controllability, while is a negative correlation coefficient of controllability.
Then we fix the connection probability and let the number of nodes vary across . The variation of controllability index with the number of nodes is shown in Figure 3. We can see from Figure 3 that, for the ER network model, the controllability index Q of the network decreases and the number of nodes N increases. In addition, the index Q is very poor for the network with the large number of nodes, which indicates that when the number of the nodes of the network is huge, the controllability is low and the whole network is relatively hard to control. In this regard, it is consistent with our practical experience.
4.2. WS Small-World Network
For comparison, we apply our controllability index Q and driver nodes to the other class of networks. The WS small-world model is one of complex networks with short average path lengths and high clustering . The topology of the network is relatively homogeneous, meaning that all nodes are of similar degree. Roughly speaking, start with a nearest-neighbor coupled ring lattice with N nodes, in which each node is connected to its 2K neighboring nodes , K on each side. Randomly rewire each link of the network with probability p such that self-connections and duplicated links are excluded.
For WS small-world network model, it could been seen from Figure 4 that when the number of nodes N and the coefficient K are fixed, the connectivity possibility varies from 0 to 1, which indicates that there is no direct relation between controllability Q and p. In other words, under this circumstance, the change of parameter p does not influence the controllability.
Set and , and let K change from 2 to 11. We can clearly see from the simulation results in Figure 5 that K and controllability are negatively correlated. This means that the initial network with higher density edge and the controllability are relatively weak.
Then we set and . And when N increases from 20 to 60, we can see from Figure 6 that the controllability decreases as the number of nodes increases. Therefore, we can conclude that the number of nodes and edges of this network are negatively correlated to the controllability.
4.3. BA Scale-Free Network
The algorithm of the BA scale-free model  is generated as follows. Growth: starting with a small number () of nodes, at every time step, add a new node with () edges that link the new node to m different nodes already presented in the network. Preferential attachment: when choosing the nodes to which the new node connects, assume that the probability that a new node will be connected to node i depends on the degree of node i, in such way thatAfter time steps, we get a network having nodes and mt edges. This network evolves into a scale-invariant state with the probability that a node has edges following a power-law distribution. However, the degree distribution of BA model is heterogeneous.
For the controllability of BA scale-free network model, firstly, let the start nodes number ; the connections of nodes in the initial network are as follows: both are isolated points; constitute a complete graph; randomly connect some edges. The number of new edges added to each node varies across . When the number of added nodes , the total numbers of nodes N are 15,17,19, and 21, respectively. The relationship between the controllability index of the network and m is shown in Figure 7.
In Figure 7, the simulation shows that when the number of nodes N is relatively small (), the controllability Q is positively correlated to m (the number of new edges added when one node is connected to the previous network). However, as the number of added nodes () increases, the controllability becomes negatively correlated with as shown in Figure 7.
When setting and , let the number of added nodes change from 10 to 70. Thus, the number of nodes in the BA scale-free network varies across . For each case, the average controllability index of 1000 trials is recorded and shown in Figure 8. As can be seen from the simulation results, as the number of nodes increases (i.e., the size of the network increases), the controllability index decreases. This indicates that the controllability of large-scale networks is quite weak.
5. Conclusions and Further Work
In this paper, we mainly discussed the problem of quantitative measurement controllability for arbitrary given network. A quantitative index is presented based on the control centrality and conditional number of the controllability matrix of controllable network (or subnetwork) system. And the index can quantitatively measure the controllability of each node in a given network. The effect of this controllability index is observed and discussed mainly by a series of experiments on various types of networks, namely, E–R networks, WS small-world networks, and BA scale-free networks. For ER networks, the controllability of networks increases with the increase of connection probability p; however, when p reached a certain value, the controllability will decrease with the increase of connection probability. Besides, the controllability of the network is negatively correlated with the number of nodes when the connectivity probability is fixed, so the controllability of the large-scale network is poor. For WS small-world network model, the controllability of the network is negatively correlated to both K and the number of nodes. For the BA scale-free network, due to the different number of new edges added when one node is trying to connect to the previous network, the controllability of the network presents different trends, and specific reasons of this phenomenon need to be further studied. These findings can help us better understand the relationship between controllability and network topology.
Much work still remains to be done concerning the controllability issue of complex networks. There still exist abundant potential future extensions: How to apply this index to larger network? Our method has been estimated using MATLAB R2014a, which can execute the network with 1000 nodes for an hour. When the number of nodes in the network increases to a large complex network with millions of nodes, it is necessary to further optimize the algorithm and rewrite programs in C language or any other language to reduce the complexity of the algorithm and improve the computing speed. The measuring controllability index of this paper proposed is mainly based on given complex networks and linear time-invariant dynamic systems. It cannot be directly applied to the time-dynamic complex networks. Should this index be used on time-variant networks, we need to consider whether under this circumstance the “control centrality” could reflect the dimensional change of controllable subsystems and whether conditional number can still measure the singularity of matrix. These questions need further consideration. Further research can be made to compare the controllability between complex networks with different topology. More experiments can be conducted to discover more empirical regularities of real-world networks. Certain phenomenon observed could be explained analytically or even mathematically. To be brief, the current complex network research continues to develop, but it also faces many problems that needed to be solved. How to apply these theoretical concepts and equations into practice is also of great significance in future studies.
The simulation data used to support the findings of this study are included within the article. If other data or MATLAB programs used to support the findings of this study are needed, you can obtain them from the corresponding author.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
This work is supported by Natural Science Foundation of Hebei Province (Grants nos. F2016501023 and F2017501041), Fundamental Research Funds for the Central Universities (Grant no. N172304030), and National Natural Science Foundation of China (Grant no. 61402088).
- H. Jeong, B. Tombor, R. Albert, Z. N. Oltval, and A.-L. Barabásl, “The large-scale organization of metabolic networks,” Nature, vol. 407, no. 6804, pp. 651–654, 2000.
- E. Bullmore and O. Sporns, “Complex brain networks: graph theoretical analysis of structural and functional systems,” Nature Reviews Neuroscience, vol. 10, no. 3, pp. 186–198, 2009.
- D. Stauffer, A. Aharony, L. Da Fontoura Costa, and J. Adler, “Efficient Hopfield pattern recognition on a scale-free neural network,” The European Physical Journal B, vol. 32, no. 3, pp. 395–399, 2003.
- J. Shen and B. Zheng, “Cross-correlation in financial dynamics,” EPL (Europhysics Letters), vol. 86, no. 4, p. 48005, 2009.
- D. Liben-Nowell and J. Kleinberg, “The link-prediction problem for social networks,” Journal of the Association for Information Science and Technology, vol. 58, no. 7, pp. 1019–1031, 2007.
- A. L. Barabási, H. Jeong, Z. Néda, E. Ravasz, A. Schubert, and T. Vicsek, “Evolution of the social network of scientific collaborations,” Physica A: Statistical Mechanics and its Applications, vol. 311, no. 3-4, pp. 590–614, 2002.
- S. Fortunato, “Community detection in graphs,” Physics Reports, vol. 486, no. 3–5, pp. 75–174, 2010.
- R. E. Kalman, “Mathematical description of linear dynamical systems,” Journal of the Society for Industrial & Applied Mathematics, Series A: Control, vol. 1, no. 2, pp. 152–192, 1963.
- D. G. Luenberger, Introduction to Dynamic Systems: Theory, Models and Applications, Wiley, 1979.
- J.-J. Slotine and W. Li, Applied Nonlinear Control, Prentice-Hall, 1991.
- R. Albert and A. Barabási, “Statistical mechanics of complex networks,” Reviews of Modern Physics, vol. 74, no. 1, pp. 47–97, 2002.
- M. E. Newman, D. J. Watts, and S. H. Strogatz, “Random graph models of social networks,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 99, supplement 1, pp. 2566–2572, 2002.
- S. Boccaletti, V. Latora, Y. Moreno, M. Chavez, and D. W. Hwang, “Complex networks: Structure and dynamics,” Physics Reports, vol. 424, no. 4-5, pp. 175–308, 2006.
- G. Palla, I. Derényi, I. Farkas, and T. Vicsek, “Uncovering the overlapping community structure of complex networks in nature and society,” Nature, vol. 435, no. 7043, pp. 814–818, 2005.
- M. Nagy, Z. Ákos, D. Biro, and T. Vicsek, “Hierarchical group dynamics in pigeon flocks,” Nature, vol. 464, no. 7290, pp. 890–893, 2010.
- Y. Xiao, S. Lao, L. Hou, M. Small, L. Bai, and S. Hayasaka, “Effects of Edge Directions on the Structural Controllability of Complex Networks,” PLoS ONE, vol. 10, no. 8, p. e0135282, 2015.
- T. Jia and M. Pósfai, “Connecting Core Percolation and Controllability of Complex Networks,” Scientific Reports, vol. 4, no. 1, 2015.
- Y.-Y. Liu, J.-J. Slotine, and A.-L. Barabási, “Controllability of complex networks,” Nature, vol. 473, no. 7346, pp. 167–173, 2011.
- Y. Liu, J. Slotine, A. Barabási, and Y. Moreno, “Control Centrality and Hierarchical Structure in Complex Networks,” PLoS ONE, vol. 7, no. 9, p. e44459, 2012.
- T. Jia and A.-L. Barabási, “Control capacity and a random sampling method in exploring controllability of complex networks,” Scientific Reports, vol. 3, 2013.
- Z. Z. Yuan, C. Zhao, Z. R. Di, W.-X. Wang, and Y.-C. Lai, “Exact controllability of complex networks,” Nature Communications, vol. 4, article 2447, 2013.
- G. Yan, J. Ren, Y. Lai, C. Lai, and B. Li, “Controlling complex networks: how much energy is needed?” Physical Review Letters, vol. 108, no. 21, Article ID 218703, 2012.
- Y.-Y. Liu and A.-L. Barabási, “Control principles of complex systems,” Reviews of Modern Physics, vol. 88, no. 3, 2016.
- H. Yin and S. Zhang, “Minimum structural controllability problems of complex networks,” Physica A: Statistical Mechanics and its Applications, vol. 443, pp. 467–476, 2016.
- S. P. Chin, J. Cohen, A. Albin et al., “A Mathematical Analysis of Network Controllability Through Driver Nodes,” IEEE Transactions on Computational Social Systems, vol. 4, no. 2, pp. 40–51, 2017.
- C. Zhao, A. Zeng, R. Jiang, Z. Yuan, and W. Wang, “Controllability of flow-conservation networks,” Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 96, no. 1, 2017.
- N. Cai, “On quantitatively measuring controllability of complex networks,” Physica A: Statistical Mechanics and its Applications, vol. 474, pp. 282–292, 2017.
- W. J. Rugh, Linear System Theory, Prentice Hall, Upper Saddle River, NJ, USA, 1996.
- S. Hosoe, “Determination of generic dimensions of controllable subspaces and its application,” Institute of Electrical and Electronics Engineers Transactions on Automatic Control, vol. 25, no. 6, pp. 1192–1196, 1980.
- R. A. Horn and C. R. Johnson, Matrix Analysis, Cambridge University Press, Cambridge, UK, 1985.
- C. T. Chen, Linear System Theory and Design, Oxford University Press, Oxford, UK, 1995.
- B. M. Chen and Z. Lin, “Linear systems theory: a structural decomposition approach,” vol. 16, pp. 281–283, 2006.
- D. J. Watts and S. H. Strogatz, “Collective dynamics of “small-world” networks,” Nature, vol. 393, no. 6684, pp. 440–442, 1998.
- A. Barabasi and R. Albert, “Emergence of scaling in random networks,” Science, vol. 286, no. 5439, pp. 509–512, 1999.
Copyright © 2018 Lifu Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.