Research Article  Open Access
ShihYu Li, ShiAn Chen, ChinTeng Lin, LiWei Ko, ChengHsiung Yang, HengHui Chen, "Generalized Synchronization of Nonlinear Chaotic Systems through Natural Bioinspired Controlling Strategy", Abstract and Applied Analysis, vol. 2015, Article ID 725674, 14 pages, 2015. https://doi.org/10.1155/2015/725674
Generalized Synchronization of Nonlinear Chaotic Systems through Natural Bioinspired Controlling Strategy
Abstract
A novel bioinspired control strategy design is proposed for generalized synchronization of nonlinear chaotic systems, combining the bioinspired stability theory, fuzzy modeling, and a novel, simpleform Lyapunov control function design of derived high efficient, heuristic and bioinspired controllers. Three main contributions are concluded: (1) apply the bioinspired stability theory to further analyze the stability of fuzzy error systems; the high performance of controllers has been shown in previous study by Li and Ge 2009, (2) a new Lyapunov control function based on bioinspired stability theory is designed to achieve synchronization without using traditional LMI method, which is a simple linear homogeneous function of states and the process of designing controller to synchronize two fuzzy chaotic systems becomes much simpler, and (3) three different situations of synchronization are proposed; classical master and slave Lorenz systems, slave Chen’s system, and Rossler’s system as functional system are illustrated to further show the effectiveness and feasibility of our novel strategy. The simulation results show that our novel control strategy can be applied to different and complicated control situations with high effectiveness.
1. Introduction
Since Pecora and Carroll first obtained the synchronization of two identical chaotic systems under different initial conditions [1], synchronization of chaotic systems has received significant attention. Nowadays, more and more techniques of chaos synchronization were proposed, such as active control [2–4], backstepping control method [5–8], linear error feedback control [9–11], adaptive control [12–16], and sliding mode control [17–20].
In our real world, various kinds of nonlinear chaotic biosystems describing natural phenomenon are found to have some states always positive. It means these states are always in the first quadrant, such as the three species preypredator systems [15, 16, 21], double MackeyGlass systems [17, 22, 23], energy communication system in biological research [19, 24], and virusimmune system [20]. For the three species preypredator systems, which consist of two competing preys and one predator can be described by the following set of nonlinear differential equations:where , , , , and , are the model parameters assuming only positive values, and the functions , represent the densities of the two prey species and represents the density of the predator species. The predator consumes the preys according to the response functions [25]. The prescribed model characterized by nonlinear response since amount of food consumed by predator per unit time depends upon the available food sources from the two preys and . And for the virusimmune system, a mathematical model is considered as follows:where , , and represent the population concentrations of uninfected, infected target cells and virus, respectively. We denote by the constant supply of target cells from its precursor, where is the rate at which new T cells are created. These cells have finite life time and represent the average death rates of these cells, and represent the death rate of the infected cells before viral production commences, and is the viral clearance rate constant. These target cells are assumed to grow logistically with specific growth rate and carrying capacity . Above all, their states are all positive.
As a result, a natural bioinspired stability theory is proposed for generalized synchronization in this paper. It means that there exists a given functional relationship between the states of the master and that of the slave. Via using the natural bioinspired stability theory, the new Lyapunov function is a simple linear homogeneous function of states and introduces less simulation error.
In this paper, a new Lyapunov control function based on the bioinspired stability theory is proposed to synchronize fuzzy chaotic systems. This new Lyapunov control function is designed as a simple linear homogeneous function of states and the process of designing controller to synchronize two fuzzy chaotic systems becomes much simpler. Furthermore, we design three different synchronization situations to show the effectiveness of our strategy— CASE I, synchronization of identical systems: synchronization of master fuzzy Lorenz system and slave fuzzy Lorenz system, CASE II, synchronization of two different systems: synchronization of master fuzzy Lorenz system and slave fuzzy Chen system, and CASE III, synchronization of complicated different systems: synchronization of master fuzzy Lorenz system and complicated fuzzy chaotic system (slave fuzzy Chen system and functional fuzzy Rossler’s system).
The layout of the rest of the paper is as follows. In Section 2, the natural bioinspired stability theory is introduced. In Section 3, modeling of chaotic systems is proposed. In Section 4, three simulation examples are given. In Section 5, conclusions are given.
2. Natural Bioinspired Controlling Strategy
2.1. Definition of the Stability on Partial Region
Consider the differential equations of disturbed motion of a nonautonomous system in the normal formwhere the function is defined on the intersection of the partial region (shown in Figure 1) andand , where and are certain positive constants. , which vanishes when the variables are all zero, is a real valued function of , . It is assumed that is smooth enough to ensure the existence and uniqueness of the solution of the initial value problem. When does not contain explicitly, the system is autonomous.
Obviously, () is a solution of (3). We are interested in the asymptotical stability of this zero solution on partial region (including the boundary) of the neighborhood of the origin which in general may consist of several subregions (Figure 1).
Definition 1. For any given number , if there exists a , such that on the closed given partial region whenfor all , the inequalityis satisfied for the solutions of (3) on , then the disturbed motion () is stable on the partial region .
Definition 2. If the undisturbed motion is stable on the partial region , and there exists a , so that on the given partial region whenthe equalityis satisfied for the solutions of (3) on , then the undisturbed motion ( is asymptotically stable on the partial region .
The intersection of and region defined by (7) is called the region of attraction.
Definition of Functions . Let us consider the functions given on the intersection of the partial region and the regionfor , where and are positive constants. We suppose that the functions are singlevalued and have continuous partial derivatives and become zero when .
Definition 3. If there exists and a sufficiently small , so that on partial region and , (or ), then is a positive (or negative) semidefinite, in general semidefinite, function on the and .
Definition 4. If there exists a positive (negative) definitive function on , so that on the partial region and then is a positive definite function on the partial region and .
Definition 5. If is neither definite nor semidefinite on and , then is an indefinite function on partial region and . That is, for any small and any large , can take either positive or negative value on the partial region and .
Definition 6 (bounded function ). If there exists , , so that on the partial region , we havewhere is a positive constant, then is said to be bounded on .
Definition 7 (function with infinitesimal upper bound). If is bounded, and for any , there exists , so that on when , and , we havethen admits an infinitesimal upper bound on .
2.2. Bioinspired GYC Theorem of Stability and of Asymptotical Stability on Partial Region
Theorem 8. If there can be found a definite function on the partial region for (3) and the derivative with respect to time based on these equations isthen, it is a semidefinite function on the partial region whose sense is opposite to that of , or if it becomes zero identically, then the undisturbed motion is stable on the partial region.
Proof. Let us assume for the sake of definiteness that is a positive definite function. Consequently, there exists a sufficiently large number and a sufficiently small number , such that on the intersection of partial region andand , the following inequality is satisfiedwhere is a certain positive definite function which does not depend on . Besides that, (13) may assume only negative or zero value in this region.
Let be an arbitrarily small positive number. We will suppose that in any case . Let us consider the aggregation of all possible values of the quantities , which are on the intersection of andand let us designate by the precise lower limit of the function under this condition. By virtue of (10), we will haveWe will now consider the quantities as functions of time which satisfy the differential equations of disturbed motion. We will assume that the initial values of these functions for lie on the intersection of and the regionwhere is so small thatBy virtue of the fact that , such a selection of the number is obviously possible. We will suppose that in any case the number is smaller than . Then the inequalitybeing satisfied at the initial instant will be satisfied, in the very least, for a sufficiently small , since the functions very continuously with time. We will show that these inequalities will be satisfied for all values . Indeed, if these inequalities were not satisfied at some time, there would have to exist such an instant for which this inequality would become an equality. In other words, we would haveand consequently, on the basis of (19)On the other hand, since , the inequality (9) is satisfied in the entire interval of time , and consequently, in this entire time interval . This yieldswhich contradicts (20) on the basis of (19). Thus, the inequality (23) must be satisfied for all values of .
Finally, we must point out that from the viewpoint of mathematics, the stability on partial region in general does not be related logically to the stability on whole region. If an undisturbed solution is stable on a partial region, it may be either stable or unstable on the whole region and vice versa. In specific practical problems, we do not study the solution starting within and running out of .
Theorem 9. If in satisfying the conditions of Theorem 8 the derivative is a definite function on the partial region with opposite sign to that of and the function itself permits an infinitesimal upper limit, then the undisturbed motion is asymptotically stable on the partial region.
Proof. Let us suppose that is a positive definite function on the partial region and that, consequently, is negative definite. Thus on the intersection of and the region defined by (9) and there will be satisfied not only the inequality (10), but also the following inequality: where is a positive definite function on the partial region independent of .
Let us consider the quantities as functions of time which satisfy the differential equations of disturbed motion assuming that the initial values of these quantities satisfy the inequalities (18). Since the undisturbed motion is stable in any case, the magnitude may be selected so small that for all values of the quantities remain within . Then, on the basis of (24) the derivative of function will be negative at all times and, consequently, this function will approach a certain limit, as increases without limit, remaining larger than this limit at all times. We will show that this limit is equal to some positive quantity different from zero. Then for all values of the following inequality will be satisfied:where .
Since permits an infinitesimal upper limit, it follows from this inequality thatwhere is a certain sufficiently small positive number. Indeed, if such a number did not exist, that is, if the quantity was smaller than any preassigned number no matter how small, then the magnitude , as follows from the definition of an infinitesimal upper limit, would also be arbitrarily small, which contradicts (25).
If for all values of the inequality (26) is satisfied, then (24) shows that the following inequality will be satisfied at all times:where is positive number different from zero which constitutes the precise lower limit of the function under condition (26). Consequently, for all values of we will havewhich is, obviously, in contradiction with (25). The contradiction thus obtained shows that the function approached zero as increase without limit. Consequently, the same will be true for the function as well, from which it follows directly thatwhich proves the theorem.
Discussion and Statement. Theorems 8 and 9 in Bioinspired GYC Theorem of Stability of Asymptotical Stability on Partial Region provides a new way to investigate the stability problem—to solve the stability problem in the first quadrant, called Partial Region, which is inspired via the biological behavior in nature and is used to control the states of a system which exist in the first quadrant to achieve the goal system with a set of simple controllers through designing a simpler and more convenient Lyapunov function.
3. Modeling of Chaotic Systems
In this section, the wellfamous TakagiSugeno fuzzy model is applied to model the classical Lorenz system, Chen system, and Rossler’s system for further designing the natural bioinspired controller.
3.1. Fuzzy Modeling of Master Lorenz System
For master Lorenz systemwhere , , are the parameters. When , , , and initial states are , the dynamic behavior is chaotic. Assume that and ; then Lorenz system can be exactly represented by TS fuzzy model as follows:whereand . and are fuzzy set of Lorenz system. Here, we call (31) the first liner subsystem under the fuzzy rule and (32) the second liner subsystem under the fuzzy rule. The final output of the fuzzy Lorenz system is inferred as follows and the chaotic behavior is shown in Figure 2whereWe further summarize all the parameters in Section 3.1 into Table 1.

3.2. Fuzzy Modeling of Slave Lorenz System
For slave Lorenz systemwhere , , are the parameters. When , , , and initial states are , the dynamic behavior is chaotic. Assume that and ; then Lorenz system can be exactly represented by TS fuzzy model as follows:whereand . and are fuzzy set of Lorenz system. Here, we call (37) the first linear subsystem under the fuzzy rule and (38) the second linear subsystem under the fuzzy rule. As a result, and are controllers of the first and second subsystems. The final output of the fuzzy Lorenz system is inferred as follows and the chaotic behavior without using controllers is shown in Figure 3whereWe further summarize all the parameters in Section 3.2 into Table 2.

3.3. Fuzzy Modeling of Chen System as Slave System
For Chen system as slave systemwhere , , are the parameters. When , , , and initial states are , the dynamic behavior is chaotic. Assume that and ; then Chen system can be exactly represented by TS fuzzy model as follows:Hereand . and are fuzzy set of Lorenz system. Here, we call (43) the first liner subsystem under the fuzzy rule and (44) the second liner subsystem under the fuzzy rule. As a result, and are controllers of the first and second subsystems. The final output of the fuzzy Chen system is inferred as follows and the chaotic behavior without using controllers is shown in Figure 4whereWe further summarize all the parameters in Section 3.3 into Table 3.

3.4. Fuzzy Modeling of Rossler’s System as Functional System
For Rossler’s system as functional systemwhere , , are the parameters. When , , , and initial states are , the dynamic behavior is chaotic. Assume that and ; then Chen system can be exactly represented by TS fuzzy model as follows:Hereand . and are fuzzy set of Lorenz system. Here, we call (49) the first liner subsystem under the fuzzy rule and (50) the second liner subsystem under the fuzzy rule. The final output of the fuzzy Rossler’s system is inferred as follows and the chaotic behavior is shown in Figure 5whereWe further summarize all the parameters in Section 3.4 into Table 4.

4. Numerical Simulation Results of the Bioinspired Strategy
In this section, there are three main CASEs for illustrations, Case 1: the generalized synchronization error function is , (), Case 2: the generalized synchronization error function is , (), and Case 3: the generalized synchronization error function is , (). A set of simpleform and efficient Lyapunov function based on the natural Bioinspired stability theory is directly used to achieve generalized synchronization of two fuzzy chaotic systems.
Case 1. The generalized synchronization error function is , ().
The addition of 100 makes the error dynamics always happen in first quadrant. Our goal is ; that is, The error and error dynamics areBefore control action, the error dynamics always happens in first quadrant as shown in Figure 6. By the natural bioinspired stability theory, one can choose a Lyapunov function in the form of a positive definite function in first quadrant:Its time derivative through error dynamics (55) isChoose We obtainwhich is negative definite function in the first quadrant. The numerical results are shown in Figures 7 and 8. In Figures 7 and 8, the master and slave systems are two identical Lorenz systems with different initial conditions, where the parameters , , , , , initial states are for master system, and for slave system. The control procedure is arranged that, in the initial 30 sec, the system is in a nature situation without any control. The designed controllers operate after 30 sec to show the high efficiency of our bioinspired controllers. In Figure 7, it is clear that the error states , , almost immediately smoothly approach zero after the controllers put in at 30 sec. Also, in Figure 8, the states of the slave system are all synchronized to the states of the target system with high performance.
Case 2. The generalized synchronization error function is , ().
The addition of 100 makes the error dynamics always happen in first quadrant. Our goal is ; that is, The error and error dynamics areBefore control action, the error dynamics always happens in first quadrant as shown in Figure 9. By the natural bioinspired stability theory, one can choose a Lyapunov function in the form of a positive definite function in first quadrant:Its time derivative through error dynamics (61) isChoose We obtainwhich is negative definite function in the first quadrant. The numerical results are shown in Figures 10 and 11. In Figures 10 and 11, Lorenz system with initial state in (34) is considered as master system and Chen system with initial state in (46) is referring to the slave system, where the parameters , , , , , , , . The control procedure is arranged that, in the initial 30 sec, the system is in a nature situation without any control. The designed controllers operate after 30 sec to show the high efficiency of our bioinspired controllers. In Figure 10, it is clear that the error states , , almost immediately smoothly approach zero after the controllers put in at 30 sec. Also, in Figure 11, the states of the slave system are all synchronized to the states of the target system with high performance.
Case 3. The generalized synchronization error function is , (), where the states of the slave Chen system in (46) are going to achieve generalized synchronization of a complicated designed goal system—the states of Lorenz system in (32) with the states of the functional system in (52).
The addition of 100 makes the error dynamics always happen in first quadrant. Our goal is ; that is, The error and error dynamics areBefore control action, the error dynamics always happens in first quadrant as shown in Figure 12. By the natural bioinspired stability theory, one can choose a Lyapunov function in the form of a positive definite function in first quadrant:Its time derivative through error dynamics (67) isChoose We obtainwhich is negative definite function in the first quadrant. The numerical results are shown in Figures 13 and 14. In Figures 13 and 14, a complicated goal system—Lorenz system with initial state in (34) combined with the functional system with initial state in (52), is considered as the master system and Chen system with initial state in (46) is referring to the slave system, where the parameters , , , , , , , , , , , . The control procedure is arranged that, in the initial 30 sec, the system is in a nature situation without any control. The designed controllers operate after 30 sec to show the high efficiency of our bioinspired controllers. In Figure 12, it is clear that the error states , , almost immediately smoothly approach zero after the controllers put in at 30 sec. Also, in Figure 13, the states of the slave system are all synchronized to the states of the target system with high performance.
5. Conclusions
In our daily life, there are lots of natural systems whose certain states are always positive, that is, exist in the first quadrant. Consequently, a novel control strategy by using bioinspired stability theory is proposed to achieve generalized fuzzy chaos synchronization. In this paper, a new Lyapunov function based on bioinspired control strategy is proposed to directly achieve chaos synchronization of chaotic systems. Via using this strategy, the new Lyapunov function used is a simple linear homogeneous function of states and the process of synchronization of two fuzzy chaotic systems becomes more simple and the controllers operate after 30 sec; the error states , , almost immediately smoothly approach zero. Classical Lorenz, Chen, and Rossler’s systems are used in numerical examples and are given to demonstrate the effectiveness of the proposed new strategy. Also, this paper provides a new sight of view for researchers to reference and we are going to further apply the novel bioinspired controlling strategy to other realworld applications in the following stage, such as energy communication system in biological research, virusimmune system, and HIV virusimmune system.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors would like to thank Professor HsienKeng Chen for the great assistance and precious advice, especially for the explanations and applications of the GYC Partial Region Stability Theory. Professor Chen is the Professor of the Department of Mechanical Engineering, Hsiuping University of Science and Technology. Also, this work was supported in part by the USTUCSD International Center of Excellence in Advanced BioEngineering sponsored by the Ministry of Science and Technology IRiCE Program under Grant no. NSC1032911I009101. This work was also supported in part by the Aiming for the Top University Plan of National Chiao Tung University, the Ministry of Education, Taiwan, under Contract 104W963, and supported in part by the Army Research Laboratory and was accomplished under Cooperative Agreement no. W911NF1020022.
References
 L. M. Pecora and T. L. Carroll, “Synchronization in chaotic systems,” Physical Review Letters, vol. 64, no. 8, pp. 821–824, 1990. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. R. Paul, S. Joshi, A. Jindal, S. P. Maurya, and A. Jain, “Experimental studies of active and passive flow control techniques applied in a twin airintake,” The Scientific World Journal, vol. 2013, Article ID 523759, 8 pages, 2013. View at: Publisher Site  Google Scholar
 B. Zhu, Y. Chen, and J. Zhao, “Integrated chassis control of active front steering and yaw stability control based on improved inverse nyquist array method,” The Scientific World Journal, vol. 2014, Article ID 919847, 14 pages, 2014. View at: Publisher Site  Google Scholar
 Y. Lan and Q. Li, “Chaos synchronization of a new hyperchaotic system,” Applied Mathematics and Computation, vol. 217, no. 5, pp. 2125–2132, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 S.Y. Li, C.H. Yang, C.T. Lin, L.W. Ko, and T.T. Chiu, “Adaptive synchronization of chaotic systems with unknown parameters via new backstepping strategy,” Nonlinear Dynamics, vol. 70, no. 3, pp. 2129–2143, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 H. Yang, H. Wang, and H. R. Karimi, “Robust adaptive neural backstepping control for a class of nonlinear systems with dynamic uncertainties,” Abstract and Applied Analysis, vol. 2014, Article ID 658671, 12 pages, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 X. Wu, G. Chen, and J. Cai, “Chaos synchronization of the masterslave generalized Lorenz systems via linear state error feedback control,” Physica D: Nonlinear Phenomena, vol. 229, no. 1, pp. 52–80, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 X. Yang, K.X. Peng, and C.N. Tong, “Robust backstepping control for cold rolling main drive system with nonlinear uncertainties,” Abstract and Applied Analysis, vol. 2013, Article ID 387890, 7 pages, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 Z.M. Ge and S.Y. Li, “Yang and Yin parameters in the Lorenz system,” Nonlinear Dynamics, vol. 62, no. 12, pp. 105–117, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 S.Y. Li and Z.M. Ge, “Pragmatical adaptive synchronization of different orders chaotic systems with all uncertain parameters via nonlinear control,” Nonlinear Dynamics, vol. 64, no. 12, pp. 77–87, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 S.Y. Li and Z.M. Ge, “A novel study of parity and attractor in the time reversed Lorentz system,” Physics Letters A, vol. 373, no. 44, pp. 4053–4059, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 M. R. Soltanpour and M. H. Khooban, “A particle swarm optimization approach for fuzzy sliding mode control for tracking the robot manipulator,” Nonlinear Dynamics, vol. 74, no. 12, pp. 467–478, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 C.H. Yang, Z.M. Ge, C.M. Chang, and S.Y. Li, “Chaos synchronization and chaos control of quantumCNN chaotic system by variable structure control and impulse control,” Nonlinear Analysis: Real World Applications, vol. 11, no. 3, pp. 1977–1985, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Shang, Y. Yuan, and F. Yuan, “Adaptive output feedback stabilization of nonholonomic systems with nonlinear parameterization,” Abstract and Applied Analysis, vol. 2014, Article ID 623605, 9 pages, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 K. Chakraborty, S. Jana, and T. K. Kar, “Effort dynamics of a delayinduced preypredator system with reserve,” Nonlinear Dynamics, vol. 70, no. 3, pp. 1805–1829, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 S. Jana and T. K. Kar, “A mathematical study of a preypredator model in relevance to pest control,” Nonlinear Dynamics, vol. 74, no. 3, pp. 667–683, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 Z. M. Ge, Y. T. Wong, and S.Y. Li, “Temporary lag and anticipated synchronization and antisynchronization of uncoupled timedelayed chaotic systems,” Journal of Sound and Vibration, vol. 318, no. 12, pp. 267–278, 2008. View at: Publisher Site  Google Scholar
 S. S. D. Xu, “Supertwistingalgorithmbased terminal sliding mode control for a bioreactor system,” Abstract and Applied Analysis, vol. 2014, Article ID 495680, 9 pages, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 Y. G. Zu, Z. H. Zhao, J. H. Yu, and F. J. Yang, NonLinear Ecological Modelling, Science Press, Beijing, China, 2004.
 B. Mukhopadhyay and R. Bhattacharyya, “Analysis of a virusimmune system model with two distinct delays,” Nonlinear Studies, vol. 15, no. 1, pp. 37–50, 2008. View at: Google Scholar  MathSciNet
 A. ElGohary and A. S. AlRuzaiza, “Chaos and adaptive control in two prey, one predator system with nonlinear feedback,” Chaos, Solitons & Fractals, vol. 34, no. 2, pp. 443–453, 2007. View at: Publisher Site  Google Scholar
 F. A. Khasawneh, D. A. Barton, and B. P. Mann, “Periodic solutions of nonlinear delay differential equations using spectral element method,” Nonlinear Dynamics, vol. 67, no. 1, pp. 641–658, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 D. Ghosh, “Projectivedual synchronization in delay dynamical systems with timevarying coupling delay,” Nonlinear Dynamics, vol. 66, no. 4, pp. 717–730, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 A. C. Chamgoué, R. Yamapi, and P. Woafo, “Bifurcations in a birhythmic biological system with timedelayed noise,” Nonlinear Dynamics, vol. 73, no. 4, pp. 2157–2173, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 A. S. Perelson, “Modelling viral and immune system dynamics,” Nature Reviews Immunology, vol. 2, no. 1, pp. 28–36, 2002. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 ShihYu Li 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.