1166
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS—I: REGULAR PAPERS, VOL. 53, NO. 5, MAY 2006
On Schur Stable Multivariate Polynomials J. A. Torres-Muñoz, E. Rodríguez-Angeles, and V. L. Kharitonov
Abstract—The class of stable multivariate polynomials, recently introduced by Kaczorek (1985), is the largest class of polynomials preserving the stability property under small coefficient variations. The principal goal of the contribution is to show that the class of Schur stable multivariate polynomials is the Moebius transformation of the latter one. This fundamental relation provides a vehicle to translate results known for one class to the other one. Index Terms—Multivariate polynomials, stability.
I. INTRODUCTION
T
HE CONTRIBUTION deals with the stability of multivariate polynomials in the continuous, as well as in the discrete domains. In both domains there are different stability concepts which generate different classes of stable polynomials. A detailed account of the classes may be found in [3]–[5], [8]–[10]. In the continuous domain our main attention is oriented to the class of stable multivariate polynomials introduced in [15]. Polynomials of the class possess several properties which allow an extension of results, known for univariate Hurwitz stable polynomials, see [18], [20]. Among various stability concepts used in the discrete domain for the stability analysis of multidimensional digital filters, see [1], [11], [13], one can distinguishes the class of strict sense Schur stable polynomials, see [3], [4], [10]. But, for some technical reasons this class has been defined in such a way that, in the univariate case, it does not coincide with the standard class of Schur stable (SS) polynomials, but rather with the class of anti-Schur polynomials. This fact does not permit an extension of habitual properties of SS univariate polynomials to the multivariate ones. On the other hand, the class of SS multivariate polynomials, see [13], [14], may be considered as a natural multivariate extension of the class of SS univariate polynomials. A fundamental relation of univariate Schur and Hurwitz polynomials is that they can be transformed one into the other by means of the Moebius transformation. The transformation allows to translate any stability criterion known for Hurwitz polynomials to a similar one for Schur polynomials. Such translation becomes much more complicated when one addresses the multivariate polynomials because of the diversity of the stability
notions in this case, see [5], [12]. As a consequence, several results obtained for discrete stable multivariate polynomials can not be easily transformed to the continuous counterparts. The main goal of the contribution is to show that the class of stable multivariate polynomials introduced in [15] is the Moebius transformation of the class of SS multivariate polynomials. Once the relation is established, it provides a vehicle to connect continuous and discrete stable polynomials. In the next section, some notations and preliminary statements concerning multivariate polynomials are provided. Our main result is presented in Section III. Section IV is dedicated to reformulation of some basic properties of stable multivariate polynomials in terms of SS polynomials. Some new sets of convex directions for SS multivariate polynomials are given there as well. An illustrative example, and concluding remarks end the paper. II. MULTIVARIATE POLYNOMIALS In this section, some basic notions and concepts for multivariate polynomials are recalled, see [5]. A multivariate polynomial is a finite sum of the form (1) where , and are real (or complex) coefficients. One may define, following the lexicographic order of the indices, the coefficient vector
The maximal degree of with a non zero coefficient in (1) . The vector of partial is known as its partial degree, degrees is denoted by . Polynomial (1) may be also written, with respect to the variable , as (2)
Manuscript received March 10, 2005; revised August 30, 2005. This paper was supported by CONACyT-Mexico.This paper was recommended by Associate Editor C. Hadjicostis. J. A. Torres-Muñoz and V. L. Kharitonov are with the Automatic Control Department, Center of Research and Advanced Studies (CINVESTAV-IPN) Mexico, Mexico city, Mexico (e-mail:
[email protected];
[email protected]). E. Rodríguez-Angeles is with the Engineering Faculty, Mexico State Autonomous University, CP 50130 Toluca, Mexico (e-mail:
[email protected]). Digital Object Identifier 10.1109/TCSI.2006.870219
, polynomials. Here and are named the main and free polynomial coefficients with respect to the variable , respectively. where
1057-7122/$20.00 © 2006 IEEE
coefficients
are
-variate
TORRES-MUÑOZ et al.: ON SS MULTIVARIATE POLYNOMIALS
1167
Given any vector with nonnegative integer components, , let us define the set
of constant degree -variate polynomials. With polynomial one can associate the reciprocal polynomial
Here means that all coefficients of complex conjugate numbers.
B. Discrete Stable Polynomials Hereafter, a discrete multivariate polynomial of degree will be denoted as (5) where
and the coefficient vector is . Let us introduce the polydomain
are changed by the
along with its essential boundary
A. Continuous Stable Polynomials Let us define the polydomain (3) and its essential boundary (4) Definition 1: A polynomial strict sense stable (SSS) if
is called
The class of SSS polynomials plays an important role in the stability and robust stability analysis of multidimensional systems, see [2], [6], [7]. Nevertheless, polynomials of this class may lose the stability property under arbitrary small coefficient variations. To remedy this situation, the following stability concept has been recently introduced. : Definition 2: [15] Given a polynomial 1) for : polynomial of degree is called stable if it is Hurwitz stable; : polynomial is called stable if it is SSS 2) for and additionally satisfies the following conditions: • main coefficients
•
are and
Definition 3: A polynomial
is called SS
if
In the next section, we will need the following two auxiliary statements. be a SS polynomial, and Lemma 3: Let . Then
is a SS polynomial from . Proof: To simplify our notations, let us introduce vectors and . Without loss of generality we may assume that all partial degrees of are positive. We show first that is SS. If it is not the case, has a root . This implies that has a then root . Last conclusion contradicts the Schur . stability of with respect to the Assume now that partial degree of is less than . The main coefficient in the decomvariable position
-variate stable polynomials;
The class of stable polynomials is the largest class of polynomials preserving the stability property under small coefficient variations. It turns out that several useful properties of Hurwitz stable univariate polynomials have been extended to stable multivariate polynomials, see [15]. Among them the following two will be used later on. be a stable polyTheorem 1: [15] Let such that every polynomial with nomial. There exists the coefficient vector from the -neighborhood of the coefficient vector of is also stable. be a stable polyTheorem 2: [15] Let such that polynomial has no nomial. There exists . roots in the -neighborhood of the essential boundary
is not equal to zero, , otherwise the correis less than . On the other sponding partial degree of hand, according to our degree assumption the main coefficient in the decomposition
is identically zero, . Certainly, at least one of the other coefficients in the last decomposition, say , is not equal to zero, otherwise will have a . Let us select such root in the polydomain that . As is not trivial, for any in a small neighborhood of there exists a such that point and
1168
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS—I: REGULAR PAPERS, VOL. 53, NO. 5, MAY 2006
As a result, for sufficiently small values of polynomial
the univariate
has a root of magnitude greater than 1. And we arrive at has a root . This the conclusion that contradicts the Schur stability of . So, the partial degree with respect to the variable is equal to . Same of arguments can be applied to any other of the partial degrees, too. Corollary 4: If we fix in a SS polynomial some of the variables in such a way that the magnitude of the fixed variables is greater or equal to one, then the resulting polynomial will be SS and its partial degrees will coincide with the corresponding . partial degrees of the polynomial be a SS polynomial. Lemma 5: Let in the decompoThen, the main coefficient sition
is a SS polynomial from . Proof: Ones again let us set vectors and . Observe first that polynomial is not identically zero, otherwise the partial degree of with respect to the variable is less than . Assume by contradiction that admits a root . Define the univariate polynomial
If for all , then fixing such that we obtain a root of . This contradiction in the case when for proves the Schur stability of is not trivial, then at least one of the coefficients all . If , , say is not equal to zero, . For any in a small neighborhood say there exists a point , such that of the point and It means that for sufficiently small at least one of the zeros of lies outside the unit disc of the complex has a root . plane. As a result in this case This contradicts the Schur stability of , and we arrive at the . Schur stability of the main coefficient Let us address now the degree condition. Assume by contradiction that at least one of the partial degrees of the main , is less than the corresponding partial decoefficient, , say . On the gree of the polynomial , other hand, we know that at least one of the coefficients , say , should have the partial degree such that the folequal to . Let us select a point lowing conditions hold.
•
•
the main coefficient
in the de-
is not equal to zero. the main coefficient At the point composition
in the de-
At the point composition
is not equal to zero. Now one can factorize polynomial
as follows:
For sufficiently large all coefficients of the polynomial in the square brackets are bounded, and additionally
and
As
, then for any
there exists
such that
and It means that for a sufficiently small the univariate polyhas a root which lies outside the unit nomial disc of the complex plane. And we arrive at the contradiction because SS polynomial has a root . The contradiction proves the fact that . Same , arguments can be applied to any other partial degree of is a SS polynomial of degree so the main coefficient . Remark 1: Lemma 5 also holds when one considers decomwith respect to any other of position of the SS polynomial . the variables , for III. MAIN RESULT The main goal of the contribution is to show that, in the multivariate case, the class of stable polynomials is the Moebius transformation of the class of SS polynomials. To motivate the further analysis let us recall, firstly, that it is well known that such intimate relationship holds when dealing with univariate case. be a Hurwitz stable polynomial of deTheorem 6: Let gree , then
is a SS polynomial of degree , and vice versa.
TORRES-MUÑOZ et al.: ON SS MULTIVARIATE POLYNOMIALS
1169
On the other hand, it is also true that discrete SS and continuous SSS polynomials are not associated, via the Moebius transformation, as it can be seen from the following example. Example 1: Polynomial is SSS of degree (1,1), but it is not stable in the sense of Definition 2. The transformed polynomial
is not SS because it allows a root . Even more, the class of continuous SSS polynomials is not associated in this sense with any other class of discrete stable polynomials, see [5], [12]. Before introducing our main result we need the following two auxiliary statements. be a stable polynomial. Lemma 7: Let Then
(6) . is a SS polynomial from Proof: For conciseness, define . Let in the decomposition
and
be the coefficients of
, be the coefficients of
in the decomposition
Let us consider the case when some of the components of the are equal to one, say the first of them. Polynomial root can be decomposed as follows:
where all terms constituting have no factor . It -variate polynomial is the main cois clear that with respect to the variefficient in the decomposition of ables , and by Definition it is a stable polynomial of degree . Now, we can easily check that
Last equality implies that stable polynomial
has a root
This contradiction proves the Schur stability of . be a SS polynomial. Lemma 8: Let Then
(7) . is a stable polynomial from , Proof: Let , and . We prove the statement by induction the with respect to the number of the variables. For statement follows from Theorem 6. Assume now that the variables, statement is true for polynomials with less than and consider the case of -variate polynomials, i.e., , . It follows from (7) that
It follows from (6) that
The fact that the inequality
This implies that
Polynomial is stable, so . As a result, the is equal to . degree of . Assume by contraNow we prove the Schur stability of diction that has a root . Observe first that not all of the components of the root are equal to one. This directly follows from the formula
If all the components of the root are different from one then has a root
This contradicts the stability of
.
directly follows from
Let us check now the strict sense stability of . Assume by contradiction that has a root . If all components has a root of the root are different of one, then
and we arrive at a contradiction with the Schur stability of . Consider now the case when some components of the root are equal to one, say the first of them. As components have positive real parts, one may assume, without any loss of generality, that , . Define the polynomial
1170
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS—I: REGULAR PAPERS, VOL. 53, NO. 5, MAY 2006
By definition . If for all , , then selecting in the right half complex plane values , we obtain a new root of which all components are different of one. Then, applying the previous arguments we arrive at the contradiction with the Schur . Let now polynomial be not trivial. Then stability of there exists a vector all components of which are different of zero, such that for sufficiently small
Now, polynomial
has a root
in a small neighborhood of the point , so is a root of . When is small enough we may conclude that this root belongs to the , and all components of the root are now difpolydomain ferent of one. Then, this again contradicts the Schur stability of , so polynomial is SSS. The main coefficient, , in the decomposition
can be written as
Polynomial
is -variate SS of degree , see Lemma 3. Therefore, by the inducis a -variate stable polynomial tion hypothesis of degree . The same is true for the main with respect to any coefficient in the decomposition of other of the variables. This ends the proof of the Lemma. Now we are ready to formulate the main result of the contribution. Theorem 9: Polynomial is stable if and only if the polynomial
is SS polynomial from . Proof: The statement directly follows from Lemmas 7 8.
be a SS polynomial. Theorem 10: Let such that every polynomial with a coefficient There exists vector from the -neighborhood of the coefficient vector of is SS. Proof: The property holds for continuous stable polynomials, see Theorem 1. This and Theorem 9 imply the result. It is worth to be mentioned that this result can be also derived from the Rudin’s Theorem, see [22]. The following result is known for the univariate Schur polynomials. be a SS polynomial Theorem 11: Let in the representation (5), then . Proof: The statement follows from Lemma 7 and Theorem 9. This theorem is also a consequence of the Rudin’s Theorem [22]. be a SS polynomial. Theorem 12: Let such that has no roots in the -neighThere exists . borhood of the essential boundary Proof: Continuous stable polynomials have no zeros in a small neighborhood of the essential boundary, see Theorem 2. Joining this with the basic Theorem 9 allow to arrive at the statement. Once again, the result can be also deduced from The Rudin’s Theorem [22]. Following statement shows the invariance of the Schur stability under differentiation as in the univariate case. be a SS polynomial. Theorem 13: Let , then If
is a SS polynomial from . Proof: Assume by contradiction that polynomial has a root . On the other hand, it is is a Schur univariate polynoknown that mial, with , see Corollary 4. Then, polynomial
is a SS univariate polynomial, so has no roots in . This contradiction proves the statement. Corollary 14: Let be a SS polynomial. , , then If
IV. SOME BASIC RESULTS FOR SCHUR STABLE POLYNOMIALS The statement of Theorem 9 provides a bridge to translate properties and stability results between the above mentioned polynomials classes. In this Section it is shown how the translation works. A. General Statements In this Subsection we apply Theorem 9 in order to derive several basic properties of SS -variate polynomials.
is a SS polynomial from . Observe that Theorem 13 does not imply that other coeffishould also cients in decomposition (2) of a SS polynomial be SS. is SS of Example 2: Polynomial degree (1,1). Its free coefficient with respect to the variable , , is not a SS univariate polynomial.
TORRES-MUÑOZ et al.: ON SS MULTIVARIATE POLYNOMIALS
1171
B. Discrete Convex Directions In this subsection, we use Theorem 9 in order to derive a concept of convex directions for SS -variate polynomials. Definition 4: An -variate real (or complex) polynomial with , , is said to be a real (or complex) convex direction for SS polynomials from if the following condition is satisfied: Given any SS real (or complex) polynomial such that is also SS, then is SS for all . Following statement is an extension of a known result for univariate polynomials, see [21]. , for . A polynomial Lemma 15: Let is a convex direction for SS polynomials if and only if the polynomial from
is a convex direction for stable polynomials from . Proof: The statement follows from Theorem 9 and Theorem 27 in [17]. Theorem 16: Let , for . A real (or is a real (or complex) complex) polynomial if the convex direction for SS polynomials from inequality
is a convex direction for . SS polynomials from , for . Then any Theorem 20: Let -variate polynomial of the form , are such that for at least one , is a where , . real convex direction for SS polynomials from Proof: The statement follows from Corollary 17 in [19], and Theorem 9. is called A polynomial anti-Schur if its reciprocal polynomial is SS. Theorem 21: Let , . Then any antiis a convex direction for Schur polynomial . SS polynomials from Proof: From [7, Th. 2] and Theorem 9 one may derive a necessary phase condition for SS polynomials. The statement directly follows from this condition and the definition of antiSchur polynomials. . Theorem 22: Given a SS polynomial Polynomial is a convex . direction for -variate SS polynomials from . Proof: To simplify let us denote such Assume by contradiction there is a SS polynomial that polynomial is SS, but polynomial is not SS for some , i.e., it has a root . Point , so both polynomials and are SS univariate polynomials, see Corollary 4. By construction polynomial is antiSchur univariate polynomial. So, is a convex direction for SS univariate polynomials. This implies that polynomial may have no zeros in . The contradiction proves the statement.
(8) V. EXAMPLE , where . holds for all Proof: The statement is a consequence of the corresponding result for continuous polynomials, see Theorems 36 37 in [16], and Theorem 9. A real polynomial is said to be auto-reciprocal if , where is a complex number such than . From Theorem one may derive the following class of convex directions for SS polynomials. , for . Any autoLemma 17: Let is a complex convex reciprocal polynomial direction for SS polynomials from . More classes of convex directions are described in the following statements. . Then every -variate Theorem 18: Let polynomial with for , is a convex direction for SS -variate polynomials from . Proof: The statement follows from [17, Th. 27], and Theorem 9. Corollary 19: Let , for . If , then any -variate polynomial
The aim is to illustrate the use of convex directions. For this, let us consider the following plant:
The plant is not SS because for and the denominator, ,, has a root . The goal is to find a proportional controller such that the closed-loop system is stable. In this case, the closed-loop characteristic polynomial is given by
This polynomial has a polytopic structure, so the edge theorem may be applied. Theorem 23: [23] Let a family of invariant degree, , of bivariate polynomials with set of uncertainty, , convex polytope. Then the family is SS if and only if all its edge polynomials are SS.
1172
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS—I: REGULAR PAPERS, VOL. 53, NO. 5, MAY 2006
Recall that an edge is the convex combination of two vertex and we get first the vertex polynomials polynomials
introduced in [15]. This statement allows to transform stability and robust stability results obtained for one class to the polynomials of the other one.
ACKNOWLEDGMENT The authors would like to thank an anonymous reviewer for pointing out that several properties of SS polynomials can be easily deduced from the Rudin’s Theorem. so, the edge polynomials are REFERENCES
In all of these edges, the term multiplied by is a constant or a univariate polynomial, so it is a convex direction for SS bivariate polynomials, see Theorem 18. So one only has to prove Schur stability of the vertex polynomials. , For the vertex polynomials at a value , such that we obtain the roots
then the following inequalities are obtained:
Inequalities , , 2, 3, 4, are necessary and suffi, therecient in order to guarantee that there are no roots in fore the proportional controller capable to stabilize the proposal . plant must have a gain
VI. CONCLUDING REMARKS It is shown that the class of SS multivariate polynomials is the Moebius transformation of the class of stable polynomials
[1] H. G. Ansell, “On certain two-variable generalizations of circuit theory, with applications to networks of transmission lines and lumped reactances,” IEEE Trans. Circuit Theory, vol. CT-11, no. 2, pp. 214–223, Feb. 1964. [2] S. Basu, “On the Multidimensional Generalization of robustness of scattering Hurwitz property of complex polynomials,” IEEE Trans. Circuits Syst., vol. 36, no. 11, pp. 1159–1167, Nov. 1989. [3] S. Basu and A. Fettweis, “Tests for discrete scattering Hurwitz polynomials,” Circuits Syst. Signal Process., vol. 3, pp. 225–242, 1984. [4] , “New results on stable multidimensional polynomials Part II: Discrete case,” IEEE Trans. Circuits Syst., vol. CAS-34, no. 12, pp. 1264–1274, Dec. 1987. [5] N. K. Bose, Applied Multidimensional Systems Theory. New York: Van Nostrand Reinhold, 1982. [6] , “Robust multivariate scattering Hurwitz interval polynomial,” Linear Algebr. Appl., vol. 98, pp. 123–136, 1988. [7] , “Edge property from end-points scattering Hurwitz polynomials,” Automatica, vol. 32, pp. 655–657, 1996. [8] A. Fettweis, “A new approach to Hurwitz polynomials in several variables,” Circuits Syst. Signal Process., vol. 5, pp. 405–417, 1986. [9] A. Fettweis and S. Basu, “New results on stable multidimensional polynomials Part I: Continous case,” IEEE Trans. Circuits Syst., vol. CAS-34, no. 12, pp. 1221–1232, Dec. 1987. [10] , “On discrete scattering Hurwitz polynomials,” Int. J. Circuit Theory Appl., vol. 34, pp. 47–59, 1985. [11] E. Fornasini and G. Marchesini, “Doubly-indexed dynamical systems: state-space models and structural properties,” Math. Syst. Theory, vol. 12, pp. 59–72, 1978. [12] D. Goodman, “Some difficulties with double bilinear transformations in 2-D filters design,” Proc. IEEE, vol. 66, no. 4, pp. 796–797, Apr. 1978. [13] T. S. Huang, “Stability of two-dimensional recursive digital filters,” IEEE Trans. Audio Electroacoust., vol. 20, pp. 158–163, 1972. [14] T. Kaczorek, “Two-dimensional linear systems,” in Lectures Notes in Control and Information Sciences. New York: Springer-Verlag, 1985. [15] V. L. Kharitonov and J. A. Torres Muñoz, “Robust stability of multivariate polynomials—Part 1: Small coefficient perturbations,” Multidim. Syst. Signal Process., vol. 10, pp. 7–20, 1999. [16] , “Recent results on the robust stability of multivariate polynomials,” IEEE Trans. Circuits Syst. I, Fundam. Theory Appl., vol. 49, no. 6, pp. 715–724, Jun. 2002. [17] , “Robust stability of multivariate polynomials—Part 4: Conic sets,” Multidim. Syst. Signal Process., vol. 14, pp. 343–363, 2003. [18] V. L. Kharitonov, J. A. Torres Muñoz, and M. I. Ramirez-Sosa, “Robust stability of multivariate polynomials—Part 2: Polytopic coefficient variations,” Multidim. Syst. Signal Process., vol. 10, pp. 21–32, 1999. [19] V. L. Kharitonov, J. A. Torres-Muñoz, and E. Rodríguez-Angeles, “On robust stability of multivariate interval plants,” Int. J. Robust Nonlinear Contr., vol. 13, pp. 939–950, 2003. [20] V. L. Kharitonov, M. I. Ramirez-Sosa, and J. A. Torres Muñz, “Robust stability of multivariate polynomials—Part 3: Frequency domain approach,” Multidim. Syst. Signal Process., vol. 11, pp. 213–231, 2000. [21] A. Rantzer, “Stability conditions for polytopes of polynomials,” IEEE Trans. Autom. Contr., vol. 37, no. 1, pp. 79–89, Jan. 1992. [22] P. Delsarte, Y. V. Genin, and Y. G. Kamp, “A simple proof of Rudin’s multivariable stability theorem,” IEEE Trans. Acoust. Speech Signal Process., vol. ASSP-28, pp. 701–705, [AU: MONTH/ISSUE NUMBER?] 1980. [23] Y. Xiao and R. Unbehauen, “Schur stability of polytopes of bivariate polynomials,” IEEE Trans. Circuits Syst. I, Fundam. Theroy Appl., vol. 49, no. 9, pp. 1020–1023, Sep. 2002.
TORRES-MUÑOZ et al.: ON SS MULTIVARIATE POLYNOMIALS
J. A. Torres-Muñoz received the M.S. degree in electrical engineering from the Center of Research and Advanced Studies (CINVESTAV-IPN), Mexico city, Mexico, and the Ph. D. degree in automatic control from the Polytechnic National Institute Grenoble (INPG), Grenoble, France, in 1985 and 1990, respectively. He has been a Professor in the Department of Automatic Control, CINVESTAV-IPN since 1990. He served as the Head of the Department of Automatic Control from 1999 to 2003. He was Visiting Researcher at IRCCYN, Ecole Centrale de Nantes, Nantes, France (1997–1998). His main research interests are in the structural approach of linear systems and in the stability and robust stability of multidimensional polynomials.
E. Rodríguez-Angeles was born in Hidalgo, Mexico, in 1976. He received the B.S. degree in control and automation engineering from the Escuela Superior de Ingeniería Mecánica y Eléctrica (ESIME-IPN), Mexico City, Mexico, in 1999, and the M.S. and Ph.D. degrees in in automatic control from de Center of Research and Advanced Studies (CINVESTAV-IPN), Mexico City, Mexico, in 2001 and 2004, respectively. Since early 2005, he is working in the Engineering Faculty, Autonomous University of Mexico State (UAEM), Toluca, Mexico. His research interests are in the area of robust stability of multivariate polynomials.
1173
V. L. Kharitonov received the M.A. degree in 1973, the Candidate of Science degree in automatic control in 1977, the Doctor of Science degree in 1990, all from the Leningrad State University, Leningrad, Russia. He is a Professor in the Department of Applied mathematics and Control Theory St. Petersburg State University, St Petersburg, Russia, and is currently also with the Center of Research and Advanced Studies (CINVESTAV-IPN), Mexico City, Mexico. His scientific interest includes control, stability, and robust stability.