# On Stability of Parametrized Families of Polynomials and Matrices

Authors: Handan Akyar; Taner Büyükköroğlu; Vakıf Dzhafarov
Journal: Abstract and Applied Analysis
 (2010)

---

## Abstract

The Schur and Hurwitz stability problems for a parametric polynomial
family as well as the Schur stability problem for a compact set of real
matrix family are considered. It is established that the Schur stability
of a family of real matrices𝒜 is equivalent to the nonsingularity
of the family {A2−2tA+I:A∈𝒜,t∈[−1,1]} if 𝒜 has at
least one stable member. Based on the Bernstein expansion of a
multivariable polynomial and extremal properties of a multilinear
function, fast algorithms are suggested.


---

## Body

## 1. Introduction

Letℝn (ℝ) be the set of real n vectors (numbers), ℂ the set of complex numbers. Let a polynomial family be defined by(1.1)p(z,q)=a0(q)+a1(q)z+⋯+an(q)zn,
where the uncertainty vector q belongs to a box Q(1.2)Q={(q1,q2,…,ql)∈ℝl:qi-≤qi≤qi+,i=1,2,…,l}.
Denote the set of all polynomials p(z,q) by 𝒫, that is,(1.3)𝒫={p(z,q):q∈Q}.
The family 𝒫 is said to be Schur (Hurwitz) stable if every polynomial in this family is Schur (Hurwitz) stable, that is, all roots lie in the open unit disc (open left half plane). A similar definition is valid for a matrix family where the word “roots” is replaced by “eigenvalues”.Ifai(q)=qi(i=0,1,…,n), then the family (1.3) is called an interval polynomial family. The Hurwitz stability problem of interval polynomials is solved by Kharitonov Theorem [1]. The Schur stability problem of interval polynomials has been studied in many works (see [2–5] and references therein). In [3, 5] using techniques from complex analysis, necessary conditions for the Schur stability of interval polynomials are obtained.A functiona(·):Q→ℝk is said to be multilinear if it is affine-linear with respect to each component of q∈Q. The polynomial (1.1) is called multilinear if all coefficient functions ai(q)(i=0,1,…,n) are multilinear. The family (1.1) is called polynomially parameter dependent if all coefficient functions are depending polynomially on parameters qi, (i=1,2,…,l).In [6] an algorithm for the robust Schur stability verification of polynomially parameter dependent families is given. This algorithm relies on the Bernstein expansion of a multivariable polynomial and is based on the decomposition of a polynomial into its symmetric and antisymmetric parts, and on the Chebyshev polynomials of the first and second kinds. In this paper we investigate the robust Schur stability of polynomially dependent families without employing Chebyshev polynomials (cf. [6]) (see Sections 2 and 4).The following theorems express the well-known properties of a multilinear function defined on a box.Theorem 1.1 (see [2]).
Suppose thatQ⊂ℝl is a box with extreme points {qi}, (i=1,2,…,k) and f(·):Q→ℝ is multilinear. Then both the maximum and the minimum of f(·) are attained at extreme points of Q. That is,
(1.4)maxq∈Qf(q)=maxif(qi),minq∈Qf(q)=minif(qi).Theorem1.1 leads to the following sufficient condition for stability of the multilinear family (1.3).Theorem 1.2 (see [2]).
Consider a family of polynomials𝒫 (1.3) with invariant degree, multilinear uncertainty structure, and at least one Schur (Hurwitz) stable member p(z,q0). In addition, Q is the box with extreme points {qi}, (i=1,2,…,k) given by (1.2). Then the family of multilinear polynomials 𝒫 is Schur stable (Hurwitz stable) if the Zero Exclusion Condition
(1.5)0∉conv{p(z,qi):i=1,2,…,k}
is satisfied for all z∈∂𝒟, where ∂𝒟 denotes the boundary of the unit disc (imaginary axis).It is well known that a multilinear polynomial family appears quite frequently in practical applications. The characteristic polynomial of an interval matrix is a multilinear polynomial. The mapping between the coefficient vector of a polynomial and its reflection vector is also multilinear. In [7], using this multilinearity, as well as Theorem 1.1 and a splitting procedure of the box of reflection coefficients, new conditions for the Schur stability are given.The application of Theorem1.2 for determining the stability of polynomials with multilinear coefficients yield conservative results. In [8, 9] sufficient conditions are given for ensuring that the image of a multilinear function over the box Q is a convex polygon whose edges are images of the edges of the box Q. In this case the stability can be tested by the Edge Theorem [2]. In [10], the notion of the principal point of Q is introduced, and it is shown that for a multilinear mapping f:Q→ℂ, the boundary of f(Q) is a subset of f(P), where P is the set of all principal points. Application of the results obtained to stability of one important subclass of multilinear systems is given in [10]. In [11], using the notion of generalized principal point the characterization of the smallest set of regions in the complex plane within which the roots of (1.1) lie, is given.As far as the recent works on the stability of multilinear families are concerned, we can refer the reader to [12, 13]. In [13], a multilinear family which can be expressed as the product of independent linear structures is considered. The paper suggests an elimination approach which eliminates the vertices qi∈Q that are not useful for the construction of the boundary of the value set f(Q). In [12] a sufficient condition for the zero inclusion of the value set f(Q) is given, where f:Q→ℂ is multilinear. On the basis of this condition a numerical procedure for testing the whether or not f(Q) includes the origin is given. The procedure uses the iterative subdivision of the box Q.In this paper we suggest a new simple algorithm for testing Schur stability of a multilinear family. This algorithm is based on Theorem1.1 and is sufficiently fast.The Schur stability (rather than Hurwitz stability) has the following advantage. In the obtained results, the segment[-2,2] arises naturally (see Theorem 1.1, Algorithm 3.1, Theorem 5.1) whereas the cutoff frequencies should be calculated in the Hurwitz stability problems. On the other hand the Hurwitz stability can also be tested by this algorithm, since by using the well-known transformation s=(z+1)/(z-1), the Hurwitz stability problem can be transformed into Schur stability problem (see Example 3.3 taken from [14]).In the second part of the paper, we consider the application of our approach to matrix Schur stability problem. Stability problem of matrix families has been studied in many works (see [2, 4, 15–18] and references therein). Naturally, a great deal of research has been devoted to interval matrices. Interval matrix structures are ubiquitous in nature and engineering. In [15, 18] extreme point results for Hurwitz stability are obtained which expresses the stability conditions in terms of extreme matrices. In [17], using the notion of a block P-matrix a characterization of the Schur stability of all convex combinations of Schur stable matrices is derived.We consider the Schur stability problem for a family𝒜, which is a pathwise connected real matrix family. We show that Schur stability of 𝒜 is equivalent to the nonsingularity of an extended family. A similar problem for the Hurwitz stability is considered in [16] where 𝒜 is a polytope.As pointed out above this paper addresses the following points:(1)
Robust Schur stability of polynomially dependent polynomials without involving Chebyshev polynomials (see (2.3)).(2)
A new algorithm (multilinearization) for a multilinear family (Algorithm3.1).(3)
Robust Schur stability criteria for a real matrix family via nonsingularity of an extended family.In the computational procedures, we use Theorem1.1 and the Bernstein expansion of a multivariable polynomial developed in [19–21].The examples were run on a 2.53 GHz Intel Core2 Duo with 4 GB of RAM.


## 2. Stability of a Polynomially Dependent Family

Consider the family (1.1), where ai(q) are polynomials. Here we give two polynomial equations defined on a box and show that the Schur stability is equivalent to the nonexistence of common solutions to these equations (Theorem 2.1).Suppose that the pointsz=±1 are not the roots of 𝒫 and 𝒫 has at least one stable member. Suppose that the family (1.3) is not Schur stable. Then, by continuity of roots (cf. page 52 in [2]), there exists θ∈(0,π) such that z=ejθ is a root of 𝒫, where j2=-1. Then z=e-jθ is also a root, and there exist b0,b1,…,bn-2 such that(2.1)a0(q)+a1(q)z+⋯+an(q)zn=(z-ejθ)(z-e-jθ)·(b0+b1z+⋯+bn-2zn-2)=(z2-2cosθz+1)·(b0+b1z+⋯+bn-2zn-2)
is valid. Taking t=2cosθ in (2.1), it follows that the equalities (2.2)b0=a0(q)b1-tb0=a1(q)b2-tb1+b0=a2(q)b3-tb2+b1=a3(q)⋮⋮⋮bk-tbk-1+bk-2=ak(q)⋮⋮⋮bn-2-tbn-3+bn-4=an-2(q)bn-3-tbn-2=an-1(q)bn-2=an(q)
are satisfied.Eliminatingb0,b1,…,bn-2 in (2.2) reduces the system of equations (2.2) to the system(2.3)f1(t,q)=0,f2(t,q)=0.
Thus we obtain the following result.Theorem 2.1.
Assume that the pointsz=±1 are not the roots of 𝒫 and 𝒫 has at least one stable member. Then 𝒫 is Schur stable if and only if the system (2.3) has no solution on [-2,2]×Q⊂ℝl+1.


## 3. Special Case: Stability of a Multilinear Family, Multilinearization

Let the family (1.1) be given, where the functions ai(q) are multilinear (i=0,1,2,…,n). In the case of a multilinear family, using Theorem 1.1, we can easily check whether the family (1.1) has roots z=±1.Consider the system (2.3). In (2.3) both f1 and f2 are multilinear on q and polynomially depend on t. More precisely, the greatest powers of t in (f1,f2) are(3.1)(t,t)forn=3,(t2,t)forn=4,(t2,t2)forn=5,(t3,t2)forn=6,(t3,t3)forn=7,
and so on. The advantage of system (2.3) is that it is “almost” multilinear and the variables (t,q) vary on the box [-2,2]×Q. This system can be transformed into a multilinear system by introducing new variables. Indeed, if system (2.3) contains tk as the greatest power of t, we can replace tk by the product t1t2⋯tk and add new equations t2-t1=0,t3-t1=0,…,tk-t1=0 to (2.3) (we set t=t1) (cf. [22]). This extended system will then be a multilinear system defined on a box, and Theorem 1.1 will be applicable. For example, assume that n=4. Then system (2.3) becomes (3.2)a0(q)-a2(q)+a4(q)-ta1(q)-t2a0(q)=0,a1(q)-a3(q)+ta0(q)-ta4(q)=0,
and the extended system will be(3.3)a0(q)-a2(q)+a4(q)-t1a1(q)-t1t2a0(q)=0,a1(q)-a3(q)+t1a0(q)-t1a4(q)=0,t2-t1=0,
where (t1,t2,q)∈[-2,2]×[-2,2]×Q. The family 𝒫 is Schur stable if and only if the system of multilinear equations (3.3) has no solution on the box Q̃=[-2,2]×[-2,2]×Q. By Theorem 1.1 the range of each function from (3.3) over the box Q̃ can be easily and exactly calculated. If one of these ranges does not contain zero, then system (3.3) has no solution on this box. Otherwise, the initial box Q̃ should be divided into small subboxes, and the new ranges over these small boxes should be calculated. A subbox on which at least one range does not contain zero will be eliminated since the system (3.3) has no solution on this subbox.If the family𝒫 is Schur stable, then all subboxes will be eliminated after a finite number of steps.Using the above procedure, we get the following algorithm for checking the Schur stability of a multilinear family.Algorithm 3.1.
Let the multilinear family𝒫  be given, where  Q⊂ℝl.(1)
Using Theorem1.1, check the nonexistence of the roots z=±1 (according to this theorem, only extreme points are sufficient for this checking). Otherwise, 𝒫 is unstable.(2)
Obtain the equations(3.4)f1(t,q)=0,f2(t,q)=0.(3)
Multilinearize this system by replacingt=t1 and introducing new variables t2,t3,…,tk and new equations
(3.5)f1(t1,…,tk,q)=0,f2(t1,…,tk,q)=0,t2-t1=0,t3-t1=0,…,tk-t1=0,
where (t1,t2,…,tk,q)∈Q̃=[-2,2]×⋯×[-2,2]︸k×Q⊂ℝk+l.(4)
Check for stability of (1.3) in all extreme points of the box Q̃. If there is an extreme point qi such that the polynomial p(z,qi) is not Schur stable, then stop. The family 𝒫 is unstable. Otherwise, apply the next step.(5)
Using Theorem1.1, calculate the ranges of all functions in (3.5). If at least one range does not contain zero, then stop. The family of 𝒫 is stable. Otherwise, apply the next step.(6)
Divide the boxQ̃ into two boxes in the chosen coordinate direction. For each subbox repeat the steps (2)–(5). Eliminate a subbox on which one range does not contain zero.The algorithm is finished if  (a)
all subboxes are eliminatedor(b)
for a givenɛ>0 the size of all remaining boxes are less then ɛ. In case (a), the family is stable. In case (b), it is expected that the family is unstable, and check for unstability of the family for q=q*, where q* is the center of a remaining box.We have solved a number of examples using Algorithm3.1. These examples show that this algorithm is sufficiently fast.Example 3.2.
Consider the family(3.6)p(z,q)=(30q1+40q1q2+65)z5+(5.1q1+0.01q2+q1q2+26.1)z4+(-22+0.2q1+4q1q2)z3+(6q1q2-0.02q2-10q1-10.2)z2+(-18+4q1q2-0.2q1)z+0.01q2+4.9q1+q1q2+24.1,q1∈[1,2],q2∈[1,2].
The system (3.5), in this case, becomes
(3.7)t1a5(q1,q2)+a4(q1,q2)-t1t2a0(q1,q2)-t1a1(q1,q2)+a0(q1,q2)-a2(q1,q2)=0,t1a0(q1,q2)+a1(q1,q2)-t1t2a5(q1,q2)-t1a4(q1,q2)-a3(q1,q2)+a5(q1,q2)=0,t2-t1=0,
and (t1,t2,q1,q2)∈Q̃=[-2,2]×[-2,2]×[1,2]×[1,2]⊂ℝ4. Algorithm 3.1 reports after 0.078 sec in 11 steps that this family is Schur stable.The solution processes show that if a family is unstable, the division of the boxQ̃ in all directions is more effective. In this case the division process will lead to a such value of q∈Q, for which the polynomial p(z,q) is unstable.The following example is taken from [14].Example 3.3 (see [14]).
Consider, Hurwitz stability problem for the interval matrix family(3.8)A(q)=(q1-12.06-0.060-0.25-0.0310.50.25-4-1.03000.50q2),q1∈[-1.5,-0.5],q2∈[-4,-1].
The characteristic polynomial of A(q) is
(3.9)p̃(s,q)=s4+(-q1-q2+1.06)s3+(-1.06q2+0.7809-1.06q1+q1q2)s2+(-3.7809q1-1.0309q2-0.2875+1.06q1q2)s+0.2575q1-0.00375+0.03q2+4.0309q1q2.
We apply the transformation s=(z+1)/(z-1) to the polynomial (3.9) and obtain
(3.10)p(z,q)=(z-1)4p̃(z+1z-1,q)=(-5.5834q1-3.0609q2+6.0909q1q2+2.54965)z4+(4.5318q1-0.0582q2-18.2436q1q2+6.71)z3+(22.1854q1q2+4.4157+3.665q1+2.3q2)z2+(1.32-0.1818q2-14.0036q1q2-6.5918q1)z+3.9784q1+1.0009q2+3.9709q1q2+1.00465.
Therefore, Hurwitz stability of the family (3.9) is equivalent to the Schur stability of the family (3.10) (see [2], page 221). For this example, the box Q̃ is [-2,2]×[-2,2]×[-1.5,-0.5]×[-4,-1]⊂ℝ4. By dividing the box Q̃ in all directions, Algorithm 3.1 reports after 0.7 sec that this family is not Schur stable. The polynomial p(z,q) (3.10) becomes Schur unstable for an extreme point in which q1=-1, q2=-2.5. Therefore the interval matrix family A(q) is not Hurwitz stable. Note that, in [14], this unstability has been established through the solutions of at least 1200 linear programming problems.


## 4. Bernstein Expansion

One of the methods for checking the positivity of a multivariable polynomial on a box is the method of Bernstein expansion developed in [19–21, 23]. Let us briefly describe this algorithm.LetL=(i1,i2,…,im) be an m-tuple of nonnegative integers, and for x=(x1,x2,…,xm)(4.1)xL=x1i1x2i2⋯xmim.
For N=(n1,…,nm),(4.2)L≤N⇔0≤ik≤nk(k=1,2,…,m).
An m-variate polynomial p(x) is defined as(4.3)p(x)=∑L≤NaLxL,(x∈ℝm).
Here N is called the degree of the polynomial p(x).Theith Bernstein polynomial of degree d is defined as(4.4)bd,i(x)=(di)xi(1-x)d-i,0≤i≤d.In the multivariate case, theLth Bernstein polynomial of degree N is defined by(4.5)BN,L(x)=bn1,i1(x1)⋯bnm,im(xm)(x∈ℝm).
The transformation of a polynomial from its power form (4.3) into its Bernstein form results in(4.6)p(x)=∑L≤NpL(U)BN,L(x),
where the Bernstein coefficients pL(U) of p over the m-dimensional unit box U=[0,1]×⋯×[0,1] are given by(4.7)pL(U)=∑J≤L(LJ)(NJ)aJ,(L≤N).
Here (NL) is defined as the product (n1i1)⋯(nmim). In [20], a matrix method for computing the Bernstein coefficients efficiently is described.Denote(4.8)m̲=min{p(x):x∈U},m¯=max{p(x):x∈U},α=min{pL(U):L≤N},β=max{pL(U):L≤N}.Theorem 4.1 (see [19]).
The inequalities(4.9)α≤m̲≤m¯≤β
are satisfied.Theorem4.1 gives the bounds for the range of the multivariate polynomial (4.3) over the unit box U. In order to obtain the Bernstein coefficients and bounds over an arbitrary box D, the box D should be affinely mapped onto U. As a result a new polynomial is obtained, and its Bernstein coefficients are the Bernstein coefficients of the initial polynomial p(x) (4.3) over D.In order to obtain convergent bounds for the range of polynomial (4.3) over a box D, the box D should be divided into two boxes. If the division is continued and one calculates the minimal and maximal Bernstein coefficients in each subdivision step, the calculated bounds converge to the exact bounds (provided that the diameter of subboxes tends to zero). Note that by the sweep procedure the explicit transformation of the subboxes generated by sweeps back to D is avoided.It is known that the number of Bernstein coefficients is often very large. In [23], a new technique which avoids the computation of all the Bernstein coefficients is represented for evaluating m̲ and m¯ from Theorem 4.1. The computational complexity of this technique is only nearly linear in the number of the terms of the polynomial.Bernstein expansion can also be applied in order to check whether two polynomial equations (for example, system (2.3)) have a common zero in a box D [6]. Indeed, consider system (2.3) on the box Q̃=[-2,2]×Q. By expanding f1 and f2 simultaneously into their Bernstein forms, we obtain a set of points (bI(f1,Q̃),bI(f2,Q̃)) in the plane, denoted by bI(Q̃). Then(4.10)convp(Q̃)⊆convB(Q̃),
where p(Q̃)={(f1(x),f2(x)):x∈Q̃}, B(Q̃)={bI(Q̃):I∈S}. By any standard convex hull algorithm, it is easy to check whether the origin belongs to convB(Q̃). If it is outside, then the family (1.1) is stable. Otherwise the sweep procedure is applied splitting the domain D into two subdomains on which we proceed as before.In all examples below the Bernstein coefficients have been evaluated on the basis of the results in [20].Example 4.2.
Consider the Schur stability problem for the following polynomially parameter dependent family:(4.11)p(z,q)=(q22-0.8q1)z4+(-0.07q1q23-0.05q15-0.25q2)z3+(-0.8q1+0.57q2-0.548)z2+(0.1q12+0.18q1q2-0.1)z+0.425q1∈[-0.7,0.5],q2∈[1,2].
This family is stable for q1=0,q2=1. System (2.3) now reduces to
(4.12)f1(t,q1,q2)=-0.425t2-t(0.1q12+0.18q1q2-0.1)+q22-0.57q2+0.973=0,f2(t,q1,q2)=t(-q22+0.8q1+0.425)+0.05q15+0.1q12+0.18q1q2+0.07q1q23+0.25q2-0.1=0
and (t,q1,q2)∈Q̃=[-2,2]×[-0.7,0.5]×[1,2]⊂ℝ3. We have to test for the nonexistence of a solution for this system.
After calculating the Bernstein coefficient onQ̃ for f1 and f2, we see that zero is contained in the convex hull of the Bernstein coefficients (see (4.10)). Therefore the bisection procedure must be applied to this problem. After 8 bisections in 0.421 sec, we conclude that the family is Schur stable.
On the other hand, the algorithm from [6] (by using Chebyshev polynomials) gives result after 0.51 sec in 10 steps.Example 4.3.
Consider the following multilinear family:(4.13)p(z,q)=(4q2q3+26+5q1q2-5q1)z8+(q1q2+3q1)z7+(-4q1+3+q2)z6+(-2q3+q1-6q2q3+4)z5+(2q3-5-8q1-3q1q3)z4+(-4q1q2+q1q3)z3+(3q1q2-1-4q2-q3)z2+(-2q1q2+2+12q2q3)z-2q1q2q3+q1q2-4q2q3+2,q1∈[-1,0],q2∈[-0.5,0],q3∈[-1,0]
which is Schur stable for q1=q2=q3=0. For this family, system (2.3) becomes
(4.14)f1(t,q1,q2,q3)=(2q1q2q3-q1q2+4q2q3-2)t4+(2q1q2-12q2q3-2)t3+(-6q1q2q3+5q1q2-8q2q3-5q1+4q2+q3+33)t2+(q1q2+24q2q3-q1q3+3q1+4)t+2q1q2q3-3q1q2+3q1q3+9q1-3q2-3q3-21,f2(t,q1,q2,q3)=(-2q1q2q3-4q1q2-8q2q3+5q1-24)t3+(12q2q3-3q1q2-3q1+2)t2+(4q1q2q3+11q1q2+16q2q3-6q1-5q2-q3+44)t-q1q2-6q2q3+q1q3+2q1+2q3-6,
and (t,q1,q2,q3)∈Q̃=[-2,2]×[-1,0]×[-0.5,0]×[-1,0]⊂ℝ4. We have to test for the nonexistence of a solution for this system.
After calculating the Bernstein coefficient onQ̃ for f1 and f2, we see that zero is contained in the convex hull of the Bernstein coefficients (see (4.10)). After 255 bisections in 8.892 sec, we conclude that the family is Schur stable.
On the other hand, the multilinearization algorithm gives result after 92.758 sec in 2143 steps.Example 4.4.
The family(4.15)p(z,q)=z6+(-0.2-q2-q3)z5+(0.2q3-0.1q1+0.2q2+q2q3)z4+(0.1q1q2+0.1q1q3-0.01q3-0.1q1q4-0.2q2q3+0.001)z3+(-0.001q2+0.01q1q4-0.1q1q2q3+0.01q2q3+0.001q1+0.1q1q2q4)z2+(-0.01q1q2q4-0.001q1q2-0.001q1q3)z+0.001q1q2q3,qi∈[0.1,0.2](i=1,2,3,4)
is Schur stable.
(4.16)Results.Bymultilinearizationin19steps,after1.15sec,Bernsteinexpansionin31steps,after3.744sec,usingChebyshevpolynomialsin527steps,after69.717sec.Example 4.5.
The family(4.17)p(z,q)=(200q1q2+q3-q4q6q7)z6+(30q1+40q1q2-q7+65)z5+(5.1q1q4+0.01q2+q1q2-2q6+26.1)z4+(0.2q1+4q1q2q3-22)z3+(6q1q2-0.02q2-10q1-q7-10.2)z2+(4q1q2-0.2q1q5q6-18)z+0.01q2q3+4.9q1+q4+q5+q1q2+24.1,q1∈[1.8,2],q2∈[1.5,2],q3∈[-0.5,0],q4∈[-0.5,0],q5∈[0.5,1],q6∈[0.5,1],q7∈[0,0.5]
is Schur stable.
(4.18)Results:Bymultilinearizationin13steps,after6.188sec,Bernsteinexpansionin255steps,after289.797sec.Increase in the number of parametersqi essentially increases the computational time for the Bernstein approach in comparison with multilinearization.


## 5. Stability of a Matrix Family

In this section, we apply our approach (i.e., the approach which has been applied in Sections2, 3, 4 to polynomial stability) to the matrix stability problem. Let 𝒜 be a pathwise connected set of real n×n matrices and contain at least one Schur stable member. In this section, we obtain a criterion for the Schur stability of 𝒜. A similar Hurwitz stability criterion of 𝒜 is obtained in [16], where 𝒜 is a polytope.Theorem 5.1.
Let𝒜 be a pathwise connected set of real n×n matrices and contain at least one Schur stable member A0. Then 𝒜 is Schur stable if and only if
(5.1)det(A2-2tA+I)≠0
for all A∈𝒜, t∈[-1,1].Proof.
⇒) Assume that 𝒜 is Schur stable and A∈𝒜. Then
(5.2)det(A-ejθI)(A-e-jθI)≠0
for all θ∈[0,π], where j2=-1, gives
(5.3)det(A2-2tA+I)≠0,
where t=cosθ.
⇐) Assume that (5.1) is satisfied. By contradiction if 𝒜 contains a matrix which is not Schur stable, then by continuity there exists A*∈𝒜 such that A* has an eigenvalue λ* with |λ*|=1.
Ifλ*=1 (λ*=-1), then  det(A*-I)=0 (det(A*+I)=0) which contradicts (5.1).
Ifλ*≠±1, then λ*=ejθ*, θ*∈(0,π), and we have
(5.4)0=det((A*-ejθ*I)(A*-e-jθ*I))=det(A*2-2cosθ*A*+I)
which contradicts (5.1).Let a family𝒜 be defined as(5.5)𝒜={A(q):q∈Q},
where Q is a box and A(q) polynomially depends on q. In this case by Theorem 5.1 Schur stability of 𝒜 can be tested via positivity (or negativity) of the scalar multivariate polynomial(5.6)f(t,q)=det(A2(q)-2tA(q)+I)
on Q̃=[-1,1]×Q. For this problem, Theorem 4.1 and the splitting procedure described in the previous section can be applied.The following example is taken from [24].Example 5.2 (see [24]).
Consider the interval matrix family(5.7)𝒜={(0.6q1q2q3):q1∈[0,0.2],q2∈[-0.78,0],q3∈[-0.6,0.6]}.
The matrix A0=(0.6000) belongs to 𝒜 and is Schur stable. The determinant function (5.1) is
(5.8)f(t,q1,q2,q3)=det(A2-2tA+I)=2q1q2+1.36q32-2.72tq3+1.36+q12q22+2tq1q2q3+1.2tq1q2-1.2tq32+2.4t2q3-1.2t-1.2q1q2q3-4t2q1q2.
After 22 bisections and eliminations in 1.1 sec, we decide that f(t,q1,q2,q3)>0 on [-1,1]×[0,0.2]×[-0.78,0]×[-0.6,0.6]⊂ℝ4, and by Theorem 5.1 this family is Schur stable.Example 5.3.
Consider the following quadratic family:(5.9)𝒜={A0+λA1+λ2A2:λ∈[0,1]},
where
(5.10)A0=(-0.100.20.10.50.5100.3),A1=(0-21-10-2-110),A2=(00.20.10.10.2-0.50.30.2-0.4).
The matrix A0 is Schur stable. The determinant function (5.1)
(5.11)f(t,λ)=det(A2(λ)-2tA(λ)+I)
is a two-variable polynomial having 40 terms. After 17 bisections and eliminations in 1.75 sec we conclude that f(t,λ) is positive on [-1,1]×[0,1]⊂ℝ2 and this family is Schur stable.The investigation of Schur stability of a matrix family by the guardian maps [25] requires calculating of eigenvalues of matrices having high dimensions. In Example 5.3, if the determinant of Kronecker sum A→det(A⊕B) is used as a guardian map, then a matrix whose eigenvalues should be calculated would have dimension 36×36.Theorem5.1 shows that the map(5.12)A→mint∈[-1,1]|det(A2-2tA+I)|
is a semiguardian map for the family of n×n dimensional Schur stable matrices.Note that an alternative method for checking the positivity of a multivariable polynomial on a box is given in [26].


## 6. Conclusions

We consider stability problems for a multilinear and polynomially dependent polynomial families. Two algorithms such as multilinearization and the Bernstein expansion algorithm are suggested. If the number of parameters is increased then multilinearization gives a better result. In the case of unstability, the multilinearization algorithm leads to an unstable point.A new result on Schur stability of an compact matrix family is obtained. Based on this result and the Bernstein expansion, a fast algorithm for Schur stability is given.



---

*Source: 687951-2010-10-11.xml*