New Trends on Modeling, Design, and Control of Chaotic SystemsView this Special Issue
Research Article | Open Access
Edgar-Cristian Díaz-González, Baltazar Aguirre-Hernández, Jorge Antonio López-Rentería, Eric Campos-Cantón, Carlos Arturo Loredo-Villalobos, "Stability and Multiscroll Attractors of Control Systems via the Abscissa", Mathematical Problems in Engineering, vol. 2017, Article ID 6743734, 9 pages, 2017. https://doi.org/10.1155/2017/6743734
Stability and Multiscroll Attractors of Control Systems via the Abscissa
We present an approach to generate multiscroll attractors via destabilization of piecewise linear systems based on Hurwitz matrix in this paper. First we present some results about the abscissa of stability of characteristic polynomials from linear differential equations systems; that is, we consider Hurwitz polynomials. The starting point is the Gauss–Lucas theorem, we provide lower bounds for Hurwitz polynomials, and by successively decreasing the order of the derivative of the Hurwitz polynomial one obtains a sequence of lower bounds. The results are extended in a straightforward way to interval polynomials; then we apply the abscissa as a measure to destabilize Hurwitz polynomial for the generation of a family of multiscroll attractors based on a class of unstable dissipative systems (UDS) of affine linear type.
Consider the parametric dynamical system where is the state vector, is a parameter vector, and is an enough smooth vector field. Several techniques have been proposed in the analysis of the solutions behavior of a dynamical system. The Hartman–Grobman theorem establishes that its internal evolution is determined by its Jacobian matrix. That is, the behavior of its solutions is described by the spectrum of its linearization. If all of the solutions of a dynamical system converge to an equilibrium point then it is said to be a locally asymptotically stable system. The importance of studying Hurwitz polynomials is due to its usefulness in the stability analysis of linear systems: if the characteristic polynomial of a linearized system is Hurwitz (roots with negative real part) then it is asymptotically stable. This has motivated researchers working on applications seeking such polynomials. Maxwell  posed the problem in the following way: How can one find the necessary and sufficient conditions to decide whether a polynomial has all its roots with negative real part? A solution was given by Hurwitz  and it is known as the Routh–Hurwitz criterion. Related information about Hurwitz polynomials can be found in [3–6].
The study of stability with a polynomial approach had an important impulse when Kharitonov’s theorem was published in 1978. This theorem gives conditions for the stability of an interval family of polynomials (see ). Since then, a lot of works related to this theorem have been published (see, e.g., [8–12]). The importance of studying the stability of families of polynomials can be appreciated in applications where the presence of uncertainties in the polynomial coefficients has to be taken into account. Other families of polynomials that have been investigated are the segments of polynomials (see [13–16]). Good references on families of stable polynomials are [3, 17–19]. The importance of knowing the abscissa of stability has been pointed out in [20–22]. Lower bounds were reported in [23, 24]; these are the first works about the abscissa of stability; and upper bounds were obtained in Bialas , Henrici , and Olifirov .
However, stability is not always required. For example, there is a class of chaotic dynamical systems based on unstable equilibria. Several times a structural change is given by one bifurcation parameter of that generates bifurcation in the solutions of the system. Generating chaotic behavior is the subject of interest in several areas in mathematics and engineering insomuch that researchers have taken the task of design systems with diverse techniques undergoing chaotic behavior with and without equilibria. One of the different chaotic behaviors is the presence of multiscroll attractor. Good references where the generation of multiscrolls has been studied are the works [28–35]. In this paper we use the abscissa of stability of Hurwitz polynomials to study the stability of systems in order to generate multiscroll attractors. To achieve the design of a chaotic system, a technique involving lower bounds for stabilizing and breaking down the stability to make multiscroll attractors arise is described. The rest of the paper is organized as follows: In Section 2, basic definitions and results needed for the development of our technique are given. In Section 3, the relation between the abscissas of stability and , of a Hurwitz polynomial and its derivative polynomial , respectively, is studied. Therein the relationship is the following inequality which is used to obtain a lower bound for the abscissa of stability of a polynomial or an interval family of Hurwitz polynomials. We use the Gauss–Lucas Theorem 2 to analyze the Hurwitz stability of a polynomial and its derivative. Finally, in Section 4 an application of the lower bound to generate chaos is given.
Consider an asymptotically stable linear system given bywhere is the state vector of the system and is a linear operator. Let be the characteristic polynomial of . The abscissa of stability of polynomial is given by the following definition.
Definition 1. If is a Hurwitz polynomial and are its zeros then the abscissa of stability of is defined byIf and are numbers such that , then they are named lower and upper bound, respectively.
In Section 4 we consider a polynomial , so that, varying the parameter , then we get destabilization of the polynomial and we get the generation of multiscroll. In Section 4 we give the details. Now we present a useful theorem in our results.
Theorem 2 (Gauss–Lucas ). Let be any convex polygon enclosing all the zeros of the polynomial . Then the zeros of lie in .
Remark 3. Let us recall that a set of points is convex if it contains, with any two points , in the set, the line segment joining and .
The abscissa of stability of the characteristic polynomial of system (2) gives certain minimum rate of decay. Zakian and Al-Naib indicated that in computer-aided design of dynamical and control systems the numerical computation of the abscissa of stability is required (see [21, 37–39]) to warrant stability under perturbations.
3. Main Results
3.1. Abscissa of Hurwitz Polynomials: An Inequality between and
Consider the polynomial with for all .
Theorem 4. If is a Hurwitz polynomial and and are the abscissas of stability of and , respectively, then .
Proof. Let be a Hurwitz polynomial. If are the roots of then its abscissa of stability is given by , where . Then , where . Now, by the Gauss–Lucas Theorem 2 if is Hurwitz then is Hurwitz. Consequently, if is a Hurwitz polynomial then is a Hurwitz polynomial. This implies that , as we claim.
Example 5. Consider the polynomial . The abscissa of stability of is and the abscissa of stability of is . We see that .
Example 6. Consider . The abscissa of stability of is and the abscissa of stability of is . Therefore .
Example 7. Let . The abscissa of stability of is and the abscissa of stability of is . In this case we have that .
Remark 8. Theorem 4 leads to glimpsing the following open problem: if is a Hurwitz polynomial, find necessary and sufficient conditions to make the equality hold.
3.2. A Lower Bound of the Abscissa of Stability of a Polynomial
Theorem 9. Let be a Hurwitz polynomial with positive coefficients and denote . The following inequalities hold: (a)If , then(b)If , then .
Proof. If is a Hurwitz polynomial then is a Hurwitz polynomial. By Theorem 4 we have that . But if and only if . If , thenand (a) is established. The proof of (b) follows in the same way.
Example 10. For the polynomial we have that , , , , and . By part of Theorem 9 we have thatis a lower bound of ; that is, .
Example 11. Consider . Here , , , , and . By part (b) of Theorem 9 we have that is a lower bound of ; that is, .
In fact, another way to obtain is by mean of the abscissa of stability of the th derivative . Note that Theorem 9(a) is a new lower bound for the abscissa of stability and since it depends on three coefficients of while the lower bounds in (7) only depend on two coefficients of , the bound in Theorem 9(a) is in some cases better than the bound in Theorem 9(b) as is illustrated by the following example.
Example 13. Consider the following polynomial . Here , , , and . By item (a) from Theorem 9 we have thatis a lower bound of and .
Example 14. Let . Here , , , and . By item (a) from Theorem 9 we have thatis a lower bound of and .
3.3. Lower Bounds for the Abscissa of Stability of an Interval Family of Hurwitz Polynomials
For a family of Hurwitz polynomials of degree of the formthe abscissa of stability is defined by .
Theorem 15. Consider the family of Hurwitz polynomials with , ; we have that (a) is a lower bound for the abscissa of stability of the family of polynomials; (b)if , then and are lower bounds for the abscissa of stability of the family of polynomials.
Proof. From item (b) of Theorem 9, is a lower bound for the abscissa of stability of . On the other hand, since and , we have that .
For item (b) of Theorem 9 suppose that . From the next inequalities are obtained: (1),(2),(3),(4),(5),(6).ThusThis proves Theorem 15.
Remark 16. Note that for every interval family of Hurwitz polynomials we give the lower bound . If additionally the family satisfies then we can give a second lower bound given by
Remark 17. In Theorem 15 we have two lower bounds, but there could be more lower bounds. The abscissa is the maximum of all of them. That is, another way of obtaining the abscissa of stability is to take the maximum of the lower bounds.
Example 18. Consider the family of Hurwitz polynomialswhere , , , , and . Here , , , and .
Since by part (a) of Theorem 15 we have that is a lower bound of the abscissa of stability of the family Hurwitz polynomials.
Example 19. Consider the family of Hurwitz polynomialswhere , , , and . Here , , , and . Since , by item (b) from Theorem 15 we have thatwhich is a lower bound of the abscissa of stability of the family of Hurwitz polynomials.
4. The Abscissa to Generate Instability and Multiscrolls Attractors
In the study of multiscroll attractors different aspects are interesting and one of them is when the multiscroll attractor exists for a particular set of system’s parameters; then the interest is about robustness against parametric perturbation. For instance, we would like to know the variation of the values of parameters of a given system in order to preserve the multiscroll attractor. In this direction a polynomial approach has been used to find the maximal robust dynamics  and for studying the maximum range for a set of parameters to preserve the useful instability for the generation of multiscroll attractors . Now, let us apply the abscissa approach for finding the lower bound of the abscissa of hyperbolicity and instability needed in UDS to generate multiscroll attractors. The linear system (2) under a control action is given as follows:with Hurwitz characteristic polynomial of , . Define , with . Note that is a set of polynomials such that is a Hurwitz polynomial and the abscissa of stability can be calculated byNow, by Taylor’s theorem can be rewritten asIf then has roots in the imaginary axis. Thence, the system is unstable in the interval . Let us describe the class of instabilities by considering the following system in .
Definition 20. We have the following system:where is the state vector, is a linear operator with eigenvalues , and is said to be dissipative if . The system is said to be unstable and dissipative of type I (UDS-I) if one of its eigenvalues is a negative real number and the other two are complex conjugate numbers with positive real part; and it is said to be of type II (UDS-II) if one of its eigenvalues is a positive real number and the other two are complex conjugate numbers with negative real part.
This work is based on UDS-I, so a generalization of the above definition for UDS-I with dimension greater than three can be given as follows.
Definition 21. The system given by (21) where , , and eigenvalues , , is said to be dissipative if . The system is said to be unstable and dissipative of type I (UDS-I) if of its eigenvalues are negative real numbers and the other two are complex conjugate numbers with positive real part.
Due to the relation between the linear system like (21) and its characteristic polynomial, we shall say that an -degree polynomial is dissipative if the sum of its roots is negative. In a similar way, will be a UDS-I polynomial if its roots satisfy Definition 21 for systems of type I. Notice that Definition 21 is only one possibility to define UDS considering negative real numbers.
Lemma 22. Let be a real -degree Hurwitz polynomial with roots . If is unstable and dissipative, then the following conditions are satisfied: (i).(ii).
Proof. The proof of (i) is obvious. We will focus on the proof of (ii). Firstly, it is not too hard to see that if the root of has nonzero imaginary part, then its translation and its conjugate are roots of , with . Namely, by writing , then Thence, and since , then and Therefore, if is unstable and dissipative, then , as we claim.
Remark 23. The previous lemma provides an upper bound for dissipativity. However, it may happen that in the case when , for all .
Given the fact that a Hurwitz polynomial can be perturbed to be unstable for and that is an upper bound for the dissipativity, it is possible to carry the system from stability to instability in the sense of UDS if at least one of its roots has different real part than the others. The following result is immediate from the aforementioned discussion.
Corollary 24. Consider the Hurwitz polynomial , with real roots and a pair of conjugate complex roots, say, , , for some . Then (i) is Hurwitz if and only if .(ii)If , , then is UDS if and only if .
In order to generate multiscroll attractors, let us consider the control systemwhere is the state vector, stands for a real affine vector, and with denotes a nonsingular linear matrix.
Let be the characteristic polynomial of the system, , and is the following step function:where the values ’s must be chosen in a suitable way that will be explained below. Define the linear control , where . Then the controlled system isThus, the closed-loop characteristic polynomial is given by When , is a stable matrix and but when we can obtain dissipative systems with unstable dynamics and the possibility of generating multiscroll attractors. As described in Definition 21, a system with stability index will be addressed as a system of the UDS type I. Besides, the following considerations have to be made in order to call (25) a UDS of type I that in addition generates an attractor .(a)The linear part of the system must satisfy the dissipative condition , where , , are eigenvalues of . Consider also that eigenvalues are negative real numbers, and two values are complex conjugate eigenvalues with positive real part , resulting in an unstable focus-saddle equilibrium . This type of equilibria presents a stable manifold with a fast eigendirection and an unstable manifold with a slow spiral eigendirection, where corresponds to the eigenvector of regarding the eigenvalue .(b)The affine vector must be considered as a discrete function that changes depending on which domain the trajectory is located at. Accordingly . Then a switching system based on (25) is given byThe equilibria of system (29) are , with , and each entry of the switching system is considered in order to preserve bounded trajectories of system (29). Thence, the choice of ’s in the definition of the step function will determine the commutation regions ’s that enclose each equilibrium .
The commuting system given by (29) induces in phase space the flow , , such that each forward trajectory of the initial point is the set . Furthermore, these systems have a dissipative bounded region named basin of attraction, such that the flow for every . The attractor is the largest attracting invariant subset of .
Definition 25. Consider a system given by (29) in and equilibrium points , with and . We say that system (29) can generate multiscroll attractors with the minimum of equilibrium points, if for any initial condition in the basin of attraction the orbit generates an attractor with oscillations around each .
We exemplify the theory by presenting a case in where the following theorem holds.
Theorem 26. Consider system (25) for the particular case where the dimension is three. That is, consider a 3D-control system with characteristic Hurwitz polynomial , where . If , then the closed-loop system with the control is UDS for all .
A system satisfying the previous theorem is candidate to generate multiscroll attractors emerging from its equilibria with a suitable step function . The number of scrolls in the attractor is due to the step function . Next, let us illustrate the generation of multiscroll attractors. Consider the system with step function, , , where is Hurwitz.
The controlled system isDenote for . is a Hurwitz polynomial and there is no multiscroll. Figure 1 shows the projection of the stable solution onto the planes: (a) ; (b) ; and (c) .
The abscissa of is . Then other behavior could appear when . For example, for , ; hence ; consequently system (32) is dissipative when and in Figure 2 the generation of multiscroll attractor is illustrated. Another reference where multiscroll attractors have been studied is .
The equilibria of the system for are given by , , , and . Between equilibria, the commutation surfaces at the planes are as follows: , with , dividing the space into four domains given by , , , and . Notice two important facts about the system; first the scrolls are increasing their size due to the unstable manifold; this can be better appreciated at the projection of the attractor onto the plane from Figure 2(b). Second, the trajectory of the system oscillating around the equilibrium point in escapes from the domain located in the left side of the commutation surface. This occurs near the unstable manifold where it crosses the commutation surface and it is attracted by the stable manifold to the equilibrium point in the domain located at the right side of the commutation surface . The process is repeated in the inverse way forming scrolls around each equilibrium point.
In this paper we use the Gauss–Lucas theorem for obtaining an inequality between the abscissas of stability of a Hurwitz polynomial and its derivative. Then we use such inequality for getting a lower bound for the abscissa of a polynomial and for an interval family of polynomials. We have compared the lower bounds obtained with other works and we can say that the obtained bounds in this paper are easy to calculate and sometimes are better that others. Based on the aforementioned results, an approach to generate multiscroll attractors was presented. We consider that this result is important to help in understanding the emergence of chaos in stable systems. Using the abscissa of stability we can generate multiscroll attractors from a Hurwitz polynomial. One interesting aspect is that we can generate multiscroll attractor with the change of only one parameter.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
The first author wishes to thank CONACYT for its Ph.D. scholarship support. Jorge Antonio López-Rentería also wishes to thank CONACYT for the postdoctoral grant (Grant no. 290941-UIA) and the Iberoamerican University for the support in the realization of this paper. C. A. Loredo-Villalobos also wishes to acknowledge the support of CONACYT through the postdoctoral fellowship. Eric Campos-Cantón acknowledges the CONACYT financial support for sabbatical. He would also like to thank the University of Houston for his sabbatical support and Professor Matthew Nicol for allowing him to work with him and his valuable discussions on dynamic systems.
- J. C. Maxwell, “On Governors,” Proceedings of the Royal Society of London, vol. 16, pp. 270–283, 1868.
- A. Hurwitz, “Ueber die Bedingungen, unter welchen eine Gleichung nur Wurzeln mit negativen reellen Theilen besitzt,” Mathematische Annalen, vol. 46, no. 2, pp. 273–284, 1895.
- S. P. Bhattacharyya, H. Chapellat, and L. H. Keel, Robust Control: The parametric approach, Prentice-Hall, Upper Saddle River, NJ, USA, 1995.
- F. R. Gantmacher, The Theory of Matrices, vol. 1, Chelsea Publishing, New York, NY, USA, 1959.
- P. Lancaster and M. Tismenetsky, The Theory of Matrices, Academic Press, New York, NY, USA, 2nd edition, 1985.
- M. Marden, Geometry of polynomials, Mathematical Surveys, No. 3, American Mathematical Society, Providence, R.I., USA, 2nd edition, 1966.
- V. L. Kharitonov, “Asymptotic stability of an equilibrium position of a family of systems of linear differential equations,” Differentsial'nye Uravneniya, vol. 14, pp. 2086–2088, 1978.
- F. Blanchini, R. Tempo, and F. Dabbene, “Computation of the minimum destabilizing volume for interval and affine families of polynomials,” Institute of Electrical and Electronics Engineers. Transactions on Automatic Control, vol. 43, no. 8, pp. 1159–1163, 1998.
- F. Dabbene, B. T. Polyak, and R. Tempo, “On the complete instability of interval polynomials,” Systems & Control Letters, vol. 56, no. 6, pp. 431–438, 2007.
- A. A. Kale and A. L. Tits, “On Kharitonov's theorem without invariant degree assumption,” Automatica. A Journal of IFAC, the International Federation of Automatic Control, vol. 36, no. 7, pp. 1075-1076, 2000.
- R. Tempo, “A dual result to Kharitonov's theorem,” Institute of Electrical and Electronics Engineers. Transactions on Automatic Control, vol. 35, no. 2, pp. 195–198, 1990.
- J. C. Willems and R. Tempo, “The Kharitonov theorem with degree drop,” Institute of Electrical and Electronics Engineers. Transactions on Automatic Control, vol. 44, no. 11, pp. 2218–2220, 1999.
- B. Aguirre, C. Ibarra, and R. Suárez, “Sufficient algebraic conditions for stability of cones of polynomials,” Systems & Control Letters, vol. 46, no. 4, pp. 255–263, 2002.
- B. Aguirre and R. Suárez, “Algebraic test for the Hurwitz stability of a given segment of polynomials,” Boletín de la Sociedad Matemática Mexicana 1 No. 2, vol. 12, no. 2, pp. 261–275, 2006.
- S. Bialas, “A necessary and sufficient condition for the stability of convex combinations of stable polynomials or matrices,” Bulletin of the Polish Academy of Sciences. Technical Sciences, vol. 33, no. 9-10, pp. 473–480, 1985.
- J.-A. López-Renteria, B. Aguirre-Hernández, and F. Verduzco, “The boundary crossing theorem and the maximal stability interval,” Mathematical Problems in Engineering, Article ID 123403, 13 pages, 2011.
- J. Ackerman, Robust Control. The Parameter Space Approach, Springer-Verlag, New York, 2002.
- R. B. Barmish, New Tools for Robustness of Linear Systems, MacMillan Publishing Co, New York, 1994.
- P. Dorato, R. Tempo, and G. Muscato, “Bibliography on robust control,” Automatica. A Journal of IFAC, the International Federation of Automatic Control, vol. 29, no. 1, pp. 201–213, 1993.
- O. Taussky-Todd, “On stable matrices,” in Colloques Internationaux Du Centre National De La Recherche Scientifique, vol. 165, pp. 75–88, Paris, France, 1968.
- V. Zakian and U. Al-Naib, “Design of dynamical and control systems by the method of inequalities,” Proceedings of the Institution of Electrical Engineers, vol. 120, no. 11, pp. 1421–1427, 1973.
- V. Zakian, “New formulation for the method of inequalities,” Proceedings of the Institution of Electrical Engineers, vol. 126, no. 6, pp. 579–584, 1979.
- G. Schrack F, Lower Bounds for the Abscissa of Stability of Stable Polynomials. Dissertation 4065 [M.sc. thesis], Eidgenossische Technische Hochschule, Zurich, 1967.
- G. F. Schrack, “Lower bounds to the abscissa of stability of a stable polynomial from symmetric functions,” SIAM Journal on Applied Mathematics, vol. 21, pp. 373–379, 1971.
- S. Bialas, “Upper bounds for the abscissa of stability of a stable interval polynomial,” Bulletin of the Polish Academy of Sciences. Mathematics, vol. 32, no. 1-2, pp. 1–9, 1984.
- P. Henrici, “Upper bounds for the abscissa of stability of a stable polynomial,” SIAM Journal on Numerical Analysis, vol. 7, pp. 538–544, 1970.
- K. L. Olifirov, “Determination of a neighborhood of the imaginary axis which is disjoint from the spectrum of a real polynomial,” Mathematical Notes of the Academy of Sciences of the USSR, vol. 22, no. 2, pp. 581–584, 1977.
- E. Campos-Cantón, J. G. Barajas-Ramírez, G. Solís-Perales, and R. Femat, “Multiscroll attractors by switching systems.,” Chaos (Woodbury, N.Y.), vol. 20, no. 1, p. 013116, 2010.
- E. Campos-Cantón, R. Femat, and G. Chen, “Attractors generated from switching unstable dissipative systems,” Chaos, vol. 22, no. 3, Article ID 033121, 2012.
- J. Lü, F. Han, X. Yu, and G. Chen, “Generation 3-D multi-scroll chaotic attractors: a hysteresis series switching methods,” Automatica, vol. 40, no. 10, pp. 1677–1687, 2014.
- L. J. Ontanon-García, E. Jimenez-López, E. Campos-Cantón, and M. Basin, “A family of hyperchaotic multi-scroll attractors in Rn,” Applied Mathematics and Computation, vol. 233, pp. 522–533, 2014.
- C. Sánchez-López, R. Trejo-Guerra, J. M. Muñoz-Pacheco, and E. Tlelo-Cuautle, “N-scroll chaotic attractors from saturated function series employing CCII+s,” Nonlinear Dynamics, vol. 61, no. 1-2, pp. 331–341, 2010.
- J. A. K. Suykens and J. Vandewalle, “Generation of n-Double Scrolls (n = 1, 2, 3, 4,…),” IEEE Transactions on Circuits and Systems I: Fundamental Theory and Applications, vol. 40, no. 11, pp. 861–867, 1993.
- J. A. K. Suykens, A. Huang, and L. O. Chua, “A family of n-scrolls attractors from a generalized Chuas circuit,” International Journal of Electronics and Communications, vol. 51, no. 3, pp. 131–138, 1997.
- M. s. Yalcin, J. A. Suykens, J. Vandewalle, and S. Özoguz, “Families of scroll grid attractors,” International Journal of Bifurcation and Chaos in Applied Sciences and Engineering, vol. 12, no. 1, pp. 23–41, 2002.
- H. S. Wilf, Mathematics for the physical sciences, Dover Publications, Inc., 180 Varick Street, New York, NY, USA, 1978.
- V. Zakian, “Computation of the abscissa of stability by repeated use of the Routh test,” Institute of Electrical and Electronics Engineers. Transactions on Automatic Control, vol. 24, no. 4, pp. 604–607, 1979.
- V. Zakian, “A criterion of approximation for the method of inequalities,” International Journal of Control, vol. 37, no. 5, pp. 1103–1111, 1983.
- V. Zakian, Control Systems Design: A New Framework, Springer–Verlag, London, UK, 2005.
- B. Aguirre-Hernández, E. Campos-Cantón, J. A. López-Renteria, and E. C. Díaz González, “A polynomial approach for generating a monoparametric family of chaotic attractors via switched linear systems,” Chaos, Solitons & Fractals, vol. 71, pp. 100–106, 2015.
- E. C. Díaz-González, J.-A. López-Rentería, E. Campos-Cantón, and B. Aguirre-Hernández, “Maximal unstable dissipative interval to preserve multi-scroll attractors via multi-saturated functions,” Journal of Nonlinear Science, vol. 26, no. 6, pp. 1833–1850, 2016.
Copyright © 2017 Edgar-Cristian Díaz-González 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.