Next Article in Journal
Some Properties of a Function Originating from Geometric Probability for Pairs of Hyperplanes Intersecting with a Convex Body
Previous Article in Journal
A Note on Some Solutions of Copper-Water (Cu-Water) Nanofluids in a Channel with Slowly Expanding or Contracting Walls with Heat Transfer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fixed Order Controller for Schur Stability

by
Taner Büyükköroğlu
Department of Mathematics, Faculty of Science, Anadolu University, Eskisehir 26470, Turkey
Math. Comput. Appl. 2016, 21(2), 25; https://doi.org/10.3390/mca21020025
Submission received: 25 April 2016 / Revised: 30 May 2016 / Accepted: 13 June 2016 / Published: 20 June 2016

Abstract

:
If the characteristic polynomial of a discrete-time system has all its roots in the open unit disc of the complex plane, the system is called Schur stable. In this paper, the Schur stabilization problem of closed loop discrete-time system by affine compensator is considered. For this purpose, the distance function between the Schur stability region and the affine controller subset is investigated.

1. Introduction

A polynomial p ( s ) = s n + a n s n 1 + + a 2 s + a 1 with real coefficients is called Schur stable if all its roots lie in the open unit disc of the complex plane. Schur stability is very important in the investigation of discrete-time systems. One of the ways to investigate of Schur stability is the use of the well-known bilinear transformation between the open left half plane and the open unit disc. By using this transformation, the Schur stability problem of a given system can be transformed into Hurwitz stability problem of transformed polynomial (Recall that a system is called Hurwitz stable if all roots of its characteristic polynomial belong to the open left half plane) [1,2,3,4,5].
Each polynomial p ( s ) corresponds to an n-dimensional vector p = ( a 1 , a 2 , , a n ) T R n . The vector p is called stable if the corresponding monic polynomial p ( s ) is Schur stable. Denote the set of such stable vectors by D . In other words D = { ( a 1 , a 2 , , a n ) T R n : The polynomial p ( s ) = s n + a n s n 1 + + a 2 s + a 1 is Schur stableg}. It is well-known that the set D is non-convex for n > 2 and open. The closed convex hull of D is a polytope [6,7].
Many stabilization problems of discrete-time systems by feedback controller can be reduced to the following: Consider a family
p ( s , c ) = p 0 ( s ) + c 1 p 1 ( s ) + c 2 p 2 ( s ) + + c l p l ( s )
where degree ( p 0 ) = n , degree ( p i ) < n ( i = 1 , 2 , , l ) and l < n . Find l-dimensional vector c = ( c 1 , c 2 , , c l ) T R l such that the corresponding polynomial p ( s , c ) is Schur stable. In this case, the vector c is called a stabilizing vector. The polynomial equality (1) can be written in vector-matrix form as follows. Consider column vectors p 0 , p 1 , , p l where p 0 corresponds to p 0 ( s ) and p i ( i = 1 , 2 , , l ) correspond to p i ( s ) with the property that added zero component have dimension n. For instance, assume p 0 ( s ) = s 4 + 5 s 3 + 2 s 2 + s + 1 and p 1 ( s ) = 2 s 3 + s . Then p 0 = ( 1 , 1 , 2 , 5 ) T , p 1 = ( 0 , 1 , 0 , 2 ) T .
Define n × l dimensional matrix A = p 1 , p 2 , , p l . Then relation (1) can be written in vector-matrix form as follows:
p ( c ) = A c + p 0 .
The stabilization problem is the determination of a parameter c R l with the property that A c + p 0 D , namely, find conditions under which P D , where P = { A c + p 0 : c R l } is an affine subset of R n .
In [8,9] for determination of a stabilizing parameter c , the method of random generation of stable vectors and stable segments are suggested. In [10], for stabilization of continuous-time systems part of specially selected parameters are chosen randomly and another part deterministically. The combination of sum of square and linear matrix inequalities techniques for approximation of the set of stabilizing controllers for continuous time system has been considered in [11]. For such systems, a linear programming approach is considered in [12]. Stabilizing controls based on matrix inequalities have been considered in [13,14,15].
In this paper, we consider the problem of existence of an affine controller for discrete-time systems. In Section 2, we give a simple condition for the existence of a stabilizing parameter c. In Section 3, the Euclidean distance between the sets D and P is investigated, from which some existence and nonexistence conditions are obtained.

2. Rank Condition

In this section we give condition for the existence and evaluating of a stabilizing parameter c and corresponding examples. The condition is based on Rouché-Capelli theorem [16]. Without loss of generality we can assume that p 1 , , p l are linearly independent. If these vectors are linearly dependent then number of uncertain parameters c 1 , c 2 , , c l will be reduced.
Proposition 1. 
Assume that the vectors p 1 , p 2 , , p l are linearly independent. Then there exists a stabilizing parameter c R l if and only if there exists q D such that
rank p 1 , p 2 , , p l , q p 0 = l .
Proof. 
There exists a stabilizing vector if and only if there exist q D , c R l such that
A c + p 0 = q   or   A c = q p 0 .
By Rouché-Capelli theorem, for given q D , the linear system (4) on unknown vector c has unique solution c R l if and only if
rank p 1 , p 2 , , p l , q p 0 = rank p 1 , p 2 , , p l = l .
Corollary 2. 
Assume that p 1 , p 2 , , p l are linearly independent and
rank p 0 , p 1 , , p l = l .
Then there exists a stabilizing parameter c R l .
Proof. 
In Proposition 1, take q = 0 which corresponds to the Schur stable polynomial p ( s ) = s n . Then (4) is satisfied for q = 0 D .   □
Example 1. 
Consider the proper plant
G ( s ) = s 2 + 0 . 8 s 0 . 5 s 4 + 0 . 7 s 3 + 0 . 16 s 2 0 . 67 s + 0 . 2
with PID controller C ( s ) = c 3 s 2 + c 2 s + c 1 s .
The closed loop characteristic polynomial is
p ( s , c ) = ( s 2 + 0 . 8 s 0 . 5 ) ( c 3 s 2 + c 2 s + c 1 ) + s ( s 4 + 0 . 7 s 3 + 0 . 16 s 2 0 . 67 s + 0 . 2 ) = s 5 + 0 . 7 s 4 + 0 . 16 s 3 0 . 67 s 2 + 0 . 2 s + c 1 ( s 2 + 0 . 8 s 0 . 5 ) + c 2 ( s 3 + 0 . 8 s 2 0 . 5 s ) + c 3 ( s 4 + 0 . 8 s 3 0 . 5 s 2 ) .
The polynomial p 0 ( s ) = s 5 + 0 . 7 s 4 + 0 . 16 s 3 0 . 67 s 2 + 0 . 2 s is not Schur stable. Here
p 0 = ( 0 , 0 . 2 , 0 . 67 , 0 . 16 , 0 . 7 ) T , p 1 = ( 0 . 5 , 0 . 8 , 1 , 0 , 0 ) T , p 2 = ( 0 , 0 . 5 , 0 . 8 , 1 , 0 ) T , p 3 = ( 0 , 0 , 0 . 5 , 0 . 8 , 1 ) T .
The vectors p 1 , p 2 , p 3 are linearly independent, whereas p 0 , p 1 , p 2 and p 3 are linearly dependent:
p 0 + 0 . 4 p 2 0 . 7 p 3 = 0 .
Therefore, by Corollary 2 the vector c = ( 0 , 0 . 4 , 0 . 7 ) T is a stabilizing parameter.
Example 2. 
Consider the fourth-order polynomial family
p ( s , c ) = s 4 1 . 8 s 3 + 1 . 86 s 2 0 . 52 s + c 1 ( s 3 3 s 2 s 0 . 4 ) + c 2 ( s 2 s 0 . 5 ) + c 3 ( s 3 s 2 + 0 . 5 s 0 . 6 )
where the polynomial p 0 ( s ) = s 4 1 . 8 s 3 + 1 . 86 s 2 0 . 52 s is not Schur stable. For this family, corresponding vectors are
p 0 = ( 0 , 0 . 52 , 1 . 86 , 1 . 8 ) T , p 1 = ( 0 . 4 , 1 , 3 , 1 ) T , p 2 = ( 0 . 5 , 1 , 1 , 0 ) T , p 3 = ( 0 . 6 , 0 . 5 , 1 , 1 ) T
and the vectors p 1 , p 2 , p 3 are linearly independent. For the Schur stable polynomial q ( s ) = s 4 2 . 1 s 3 + 2 . 16 s 2 1 . 02 s + 0 . 1 , the vectors p 0 q p 1 , p 2 , p 3 are linearly dependent, where q = ( 0 . 1 , 1 . 02 , 2 . 16 , 2 . 1 ) T . We have 0 . 1 p 1 + 0 . 2 p 2 0 . 4 p 3 = q p 0 and by Proposition 1 the vector c = ( 0 . 1 , 0 . 2 , 0 . 4 ) T is a stabilizing parameter.

3. Schur-Szegö Parameters and Polynomial Optimization

Schur-Szegö parameters [17] or reflection coefficients are widely used in the stability problems of discrete time systems. Let us briefly recall these coefficients. In [18], for k i R ( i = 1 , 2 , , n ) the reflection map f = ( f 1 , f 2 , , f n ) T is defined recursively:
( f 1 , f 2 , , f n ) T = R n ( k n ) 0 T R n 1 ( k n 1 ) 0 T R 1 ( k 1 ) 0 1
where R j ( k j ) = I j + 1 + k j E j + 1 , I j is j × j unit matrix and, E j is j × j dimensional, defined by
E j = 0 1 1 0 .
In the above formula, the first zero is n-dimensional, whereas the last is one-dimensional. The above formula gives
for n = 3 ,
f 1 ( k 1 , k 2 , k 3 ) = k 3 , f 2 ( k 1 , k 2 , k 3 ) = k 1 k 2 k 3 + k 1 k 3 k 2 , f 3 ( k 1 , k 2 , k 3 ) = k 1 k 2 + k 2 k 3 k 1 ,
for n = 4 ,
f 1 ( k 1 , k 2 , k 3 , k 4 ) = k 4 , f 2 ( k 1 , k 2 , k 3 , k 4 ) = k 1 k 2 k 4 k 2 k 3 k 4 + k 1 k 4 k 3 , f 3 ( k 1 , k 2 , k 3 , k 4 ) = k 1 k 2 k 3 k 4 k 1 k 2 k 3 k 1 k 3 k 4 + k 1 k 3 + k 2 k 4 k 2 , f 4 ( k 1 , k 2 , k 3 , k 4 ) = k 1 k 2 + k 2 k 3 + k 3 k 4 k 1
and so on. The polynomial p ( s ) = s n + f n s n 1 + + f 2 s + f 1 is Schur stable if and only if | k i | < 1 for all i = 1 , 2 , , n . The map f ( k ) between the reflection vector k = ( k 1 , k 2 , , k n ) T and the coefficient vector p = ( f 1 , f 2 , , f n ) T is multilinear. By known property of multilinear maps [1] (p. 435) namely, the convex hull of the image of a multilinear map defined on a box Q is equal the convex hull of the images of vertices of Q. In [6], it has been proven that there is no need to take all vertices of Q, the convex hull of the image of f is a polytope, whose vertices correspond to the n + 1 polynomials for which zeros are equal to 1 or 1.
By using Schur-Szegö parameters, we can generate an arbitrary number of Schur stable polynomials: For this purpose it is sufficient to choose any vector k = ( k 1 , k 2 , , k n ) ( 1 , 1 ) × ( 1 , 1 ) × × ( 1 , 1 ) and find the image f ( k ) . The obtained vector f ( k ) is stable.
Consider the distance function between the sets P and D n . Recall that P is closed affine set whereas D n is an open set. We show that the distance function between them is n-variable polynomial with total degree 2 n . The variables of this polynomial are Schur-Szegö parameters k = ( k 1 , k 2 , , k n ) ( 1 , 1 ) n . We consider the minimization of this polynomial over the closed cube [ 1 , 1 ] n and use the Bernstein expansion for the outer approximation of the range.
Consider the number
α : = inf c R l , k [ 1 , 1 ] n A c + p 0 f ( k 1 , k 2 , , k n )
where c = ( c 1 , c 2 , , c l ) T , A is n × l matrix with column vectors p 1 , p 2 , , p l , i.e., A = p 1 , p 2 , , p l , x is the Euclidean norm of x, f is multilinear reflection map.
Theorem 3. 
The equality
α = min k [ 1 , 1 ] n A ( A T A ) 1 A T ( f ( k 1 , k 2 , , k n ) p 0 ) + p 0 f ( k 1 , k 2 , , k n )
is satisfied.
Proof. 
Write
α = inf k inf c A c y = inf k F ( k 1 , k 2 , , k n )
where y = f ( k 1 , k 2 , , k n ) p 0 , F ( k 1 , k 2 , , k n ) = inf c A c y .
Define l-dimensional subspace W = span { p 1 , p 2 , , p l } . Then F ( k 1 , k 2 , , k n ) is the distance from the point y to W. By the well-known theorem of functional analysis [19], the nearest point y ^ W exists and satisfies ( y y ^ ) W , where x W means x , w = 0 for all w W and the symbol · , · stands for the scalar product.
We have
( y y ^ ) p i , ( p i ) T ( y y ^ ) = 0 ( i = 1 , 2 , , l )
or
A T ( y y ^ ) = 0 , A T y = A T y ^ .
Since y ^ W , where y ^ = c 1 * p 1 + c 2 * p 2 + + c l * p l = A c * , where c * = ( c 1 * , c 2 * , , c l * ) T , we have
A T y = A T A c * .
The matrix A T A is nonsingular. By contradiction, assume that there exists nonzero x R l such that A T A x = 0 . Then
A x 2 = A x , A x = ( A x ) T A x = x T A T A x = x T 0 = 0
and A x = 0 . The last equality contradicts to the linear independence of p 1 , p 2 , , p l . Therefore c * = ( A T A ) 1 A T y and
F ( k 1 , k 2 , , k n ) = inf c A c y = A c * y = A ( A T A ) 1 A T ( f ( k 1 , k 2 , , k n ) p 0 ) ( f ( k 1 , k 2 , , k n ) p 0 )
and the equality (6) follows.   □
Proposition 4. 
If α > 0 then there is no stabilizing parameter c.
If α = 0 then P D ¯ and either there exists a stabilizing parameter c or there is no stabilizing parameter c, but there exists a parameter c such that p ( s , c ) is marginally stable, i.e., has all roots in the closure of the open unit disc. It should be noted that the last case is rather a pathological than a typical one.
Consider the number
α 2 = min k [ 1 , 1 ] n F 2 ( k 1 , k 2 , , k n ) .
The function ( k 1 , k 2 , , k n ) F 2 ( k 1 , k 2 , , k n ) is a multivariable polynomial defined on the box [ 1 , 1 ] n . The range F 2 ( k 1 , k 2 , , k n ) can be estimated by the Bernstein expansion. Let us briefly describe this expansion for n-variate polynomials [20].
An n-variate polynomial v ( x ) is defined as
v ( x ) = L N a L x L ( x R n ) ,
where L = ( i 1 , i 2 , , i n ) , N = ( s 1 , s 2 , , s n ) , x L = x 1 i 1 x 2 i 2 x n i n and
L n 0 i j s j ( j = 1 , 2 , , n ) .
The Lth Bernstein polynomial of degree d is defined by
B N , L ( x ) = b s 1 , i 1 ( x 1 ) b s n , i n ( x n ) ( x R n )
where b s , i ( t ) = n i t i ( 1 t ) n i . The transformation of a polynomial from its power from (7) into its Bernstein form result in
v ( x ) = L N v L ( U ) B N , L ( x ) ,
where the Bernstein coefficients v L ( U ) of v over the n-dimensional unit box U = [ 0 , 1 ] × × [ 0 , 1 ] are given by
v L ( U ) = j N L J N J a J ( L N ) .
Here N L is defined as s 1 i 1 s n i n .
Let
m ̲ = min { v ( x ) : x U } , m ¯ = max { v ( x ) : x U } , α = min { v L ( U ) : L N } , β = max { v L ( U ) : L N } .
Theorem 5. 
[20] The following inequalities
α m ̲ m ¯ β
are satisfied.
Theorem 5 gives outer approximation of the range of v ( x ) over the unit box U.
In order to obtain the Bernstein coefficients and bounds over an arbitrary box D rather the unit box U, the box D should be affinely mapped onto U. To obtain convergent bounds for the range of the polynomial (7) over the box U, the box U should be divided into small boxes.
If α > 0 ( β < 0 ) then the polynomial is positive (negative) on U. If α 0 , β 0 then by the bisection in the chosen coordinate direction, the box U is divided into two boxes. A new box on which the inequality α > 0 or β < 0 is satisfied should be eliminated, since our polynomial has constant sign on this box. Otherwise, the box should be divided into two new boxes.
If at some step of the Bernstein expansion the lower bound is positive then α 2 > 0 and α > 0 and there is no stabilizing c. If this is not the case, an additional investigation on existence is required.
Taking into account the definition of α and Proposition 4, we can suggest the following algorithm for a stabilizing vector.
Algorithm 6. 
(1)
Given family (1), explicitly calculate the multivariable polynomial F 2 ( k 1 , k 2 , , k n ) .
(2)
Calculate step-by-step the Bernstein coefficients for the function F 2 over the box [ 1 , 1 ] n . If at some step, the minimal Bernstein coefficient is positive then stop, there is no stabilizing parameter c.
(3)
If after a sufficiently large number of steps the lower Bernstein coefficient remains negative, then stop and carry out an additional test for the existence of a stabilizing vector c. A cycling of the calculation indicates the existence of a stabilizing parameter. In this case, we can proceed to a random search if number of remaining boxes is small. For example, choose the center point k * of a remaining box, calculate c * = ( A T A ) 1 ( f ( k * ) U 0 ) and test c * for a stabilizing vector, i.e., test Schur stability of p ( s , c * ) . Alternatively, the gradient minimization method of the smooth function F 2 could be applied.
Example 3. 
Consider stabilizing problem for the following polynomial family:
p ( s , c ) = s 5 + ( c 1 + c 2 ) s 4 + ( 2 + 3 c 1 c 2 ) s 3 + ( 2 c 1 + 4 c 2 ) s 2 + ( 1 c 1 5 c 2 ) s 1 + c 1 3 c 2 .
The matrix A (see (2)) and the vector p 0 are
A = 1 3 1 5 2 4 3 1 1 1 , p 0 = 1 1 0 2 0
and we have
F 2 ( k 1 , k 2 , k 3 , k 4 , k 5 ) = A ( A T A ) 1 A T ( f ( k 1 , k 2 , k 3 , k 4 , k 5 ) p 0 ) + p 0 f ( k 1 , k 2 , k 3 , k 4 , k 5 ) 2 = 83 192 k 1 2 k 2 2 k 4 2 k 5 2 + 9 16 k 1 2 k 3 2 k 4 2 k 5 2 + + 11 8 k 1 2 k 2 k 3 k 4 2 k 5 2 61 24 k 1 k 2 k 3 2 k 4 2 k 5 + + 9 16 k 2 2 k 3 2 k 4 2 11 16 k 1 2 k 3 2 k 5 11 16 k 1 2 k 3 k 4 2 7 16 k 1 2 k 3 k 5 2 + 1 2 k 1 k 2 k 5 9 4 k 1 k 2 k 4 29 24 k 1 k 3 k 5 9 4 k 2 k 3 k 4 29 24 k 2 k 4 k 5 31 16 k 2 k 3 k 5 1 2 k 5 + 5 3 .
When Algorithm 6 is applied, 477 subboxes remained at the end of 1000 steps after 48 seconds. Consider one of these subboxes:
[ 0 . 625 , 0 . 75 ] × [ 0 . 625 , 0 . 75 ] × [ 0 . 875 , 1 ] × [ 0 . 5 , 0 . 75 ] × [ 0 . 5 , 0 . 75 ] .
The center of this subbox is k * = ( 0 . 6875 , 0 . 6875 , 0 . 9375 , 0 . 625 , 0 . 625 ) T and the corresponding c * can be calculated as
c * = ( A T A ) 1 ( f ( k * ) U 0 ) = ( 0 . 992173 , 0 . 241873 ) T .
The corresponding polynomial p ( s , c * ) is Schur stable, i.e., c * is a stabilizing vector.
Example 4 
(There is no stabilizing c) Consider the family
p ( s , c ) = s 5 + ( c 1 + c 2 ) s 4 + ( 0 . 9 3 c 1 2 c 2 ) s 3 + ( 3 c 1 4 c 2 ) s 2 + ( 1 . 8 3 c 1 + 6 c 2 ) s 1 . 1 + c 1 c 2 .
Here p 0 = ( 1 . 1 , 1 . 8 , 0 , 0 . 9 , 0 ) T , p 1 = ( 1 , 3 , 3 , 3 , 1 ) T and p 2 = ( 1 , 6 , 4 , 2 , 1 ) T . The polynomial F 2 is 5-variate quadratic polynomial. We apply the Bernstein expansion, splitting-elimination procedure and after 698 steps (during 34 . 5 second) conclude that F 2 > 0 on [ 1 , 1 ] 5 . Therefore, there is no stabilizing parameter c.
Example 5. 
Let the control system shown in the Figure 1 be given. Assume that transfer function and controller are
G ( s ) = s 3 0 . 3 s 2 0 . 1 s + 1 . 6 s 5 + 0 . 4 s 4 1 . 8 s 3 0 . 7 s 2 + 0 . 2 s 0 . 6 and C ( s ) = c 3 s 2 + c 2 s + c 1 s .
The characteristic polynomial of the closed-loop system becomes
p ( s , c ) = s 6 + ( 0 . 4 + c 3 ) s 5 + ( 1 . 8 + c 2 0 . 3 c 3 ) s 4 + ( 0 . 7 + c 1 0 . 3 c 2 0 . 1 c 3 ) s 3 + ( 0 . 2 0 . 3 c 1 0 . 1 c 2 + 1 . 6 c 3 ) s 2 + ( 0 . 6 0 . 1 c 1 + 1 . 6 c 2 ) s + 1 . 6 c 1
Let c 3 = 1 . We apply Algorithm 6, part 2 and after 3726 steps conclude that there is no stabilizing parameter c.
Let c 3 = 0 . 5 . In this case, Algorithm 6 applied to the problem of stabilization does not give a negative result after 20000 steps with 6538 subboxes, and
[ 0 . 9375 , 1 ] × [ 0 . 1875 , 0 . 25 ] × [ 0 , 0 . 125 ] × [ 0 , 0 . 125 ] × [ 0 , 0 . 125 ] × [ 0 . 625 , 0 . 75 ]
is one of them. The center of this subbox is k * = ( 0 . 96875 , 0 . 21875 , 0 . 0625 , 0 . 0625 , 0 . 0625 , 0 . 6875 ) T and calculations give c * = ( 0 . 05756065 , 0 . 91839451 , 0 . 5 ) T . Since the polynomial p ( s , c * ) is Schur stable, c * is a stabilizing vector.
Remark 1. 
Here, we indicate advantages and disadvantages of our results given in this paper. Firstly, note that Proposition 1 has a simple form. On the other hand the results existing in the literature on the existence and evaluation of Schur stable element in an affine polynomial family are mainly random search methods (see [9] and references therein). Algorithm 6 gives an answer when there does not exist a stabilizing parameter; however, it gives a useful hint as to whether such a parameter exists. The main setback of Algorithm 6 is that high-dimensional systems require long calculations.

4. Conclusions

In this paper, the stabilization problem of closed-loop discrete-time systems by affine compensator is considered. A condition based on linear dependence of system vectors is obtained. The distance function between the Schur stability region and the affine controller subset is investigated. The optimal value of this function gives the conditions for stabilizability of the system.

Acknowledgments

The author thanks the referees for many valuable comments.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Bhattacharya, S.P.; Chapellat, H.; Keel, L. Robust Control: The Parametric Approach; Prentice-Hall: Englewood Cliffs, NJ, USA, 1995. [Google Scholar]
  2. Ackermann, J. Robust Control: The Parametric Space Approach; Springer-Verlag: London, UK, 2002. [Google Scholar]
  3. Hernández, B.A.; Frías-Armenta, M.E.; Verduzco, F. On differential structures of polynomial spaces in control theory. J. Syst. Sci. Syst. Eng. 2012, 21, 372–382. [Google Scholar] [CrossRef]
  4. López-Renteria, J.A.; Aguirre-Hernández, B.; Verduzco, F. Stable polynomial curves and some properties with application in control. Bol. Soc. Matemática Mex. 2016, 1–21. [Google Scholar] [CrossRef]
  5. García, F.R.; Aguirre, B.; Suárez, R. Stabilization of linear sampled-data systems by a time-delay feedback control. Math. Probl. Eng. 2008, 2008. [Google Scholar] [CrossRef]
  6. Fam, A.T.; Meditch, J.S. A canonical parameter space for linear systems design. IEEE Trans. Autom. Control 1978, 23, 454–458. [Google Scholar]
  7. Tesi, A.; Vicino, A.; Zappa, G. Convexity properties of polynomials with assigned root location. IEEE Trans. Autom. Control 1994, 39, 668–672. [Google Scholar] [CrossRef]
  8. Petrikevich, Y.I. Randomized methods of stabilization of the discrete linear systems. Autom. Remote Control 2008, 69, 1911–1921. [Google Scholar] [CrossRef]
  9. Nurges, Ü.; Avanessov, S. Fixed-order stabilising controller design by a mixed randomized/deterministic method. Int. J. Control 2015, 88, 335–346. [Google Scholar] [CrossRef]
  10. Fujisaki, Y.; Oishi, Y.; Tempo, R. Mixed deterministic/randomized methods for fixed order controller design. IEEE Trans. Autom. Control 2008, 53, 2033–2047. [Google Scholar] [CrossRef]
  11. Henrion, D.; Sebek, M.; Kucera, V. Positive polynomials and robust stabilization with fixed-order controllers. IEEE Trans. Autom. Control 2003, 48, 1178–1186. [Google Scholar] [CrossRef]
  12. Malik, W.A.; Darbha, S.; Bhattacharyya, S.P. A linear programming approach to the synthesis of fixed-structure controllers. IEEE Trans. Autom. Control 2008, 53, 1341–1352. [Google Scholar] [CrossRef]
  13. Aguirre, B.; Ibarra, C.; Suárez, R. Stability of cones of polynomials. An application to the design of high-gain controllers for saturated systems. In Proceedings of the 39th IEEE Conference on Decision and Control, Sydney, NSW, Australia, 12–15 December 2000; pp. 1087–1092.
  14. Aguirre, B.; Ibarra, C.; Suárez, R. Sufficient algebraic conditions for stability of cones of polynomials. Syst. Control Lett. 2002, 46, 255–263. [Google Scholar] [CrossRef]
  15. Aguirre, B.; Solis-Daun, J.; Suárez, R. Stabilization of Linear Systems: A polynomial approach. In Advances in Dynamics, Instrumentation and Control: Proceedings of the 2006 International Conference (Cdic ’06), Queretaro, Mexico, 13–16 August 2006; Volume II, pp. 21–28.
  16. Klafszky, E.; Terlaky, T. The role of pivoting in proving some fundamental theorems of linear algebra. Linear Algebra Its Appl. 1991, 151, 97–118. [Google Scholar] [CrossRef]
  17. Levinson, N. The Wiener RMS error criterion in filter design and prediction. J. Math. Phys. 1946, 25, 261–278. [Google Scholar] [CrossRef]
  18. Nurges, Ü. New stability conditions via reflection coefficients of polynomial. IEEE Trans. Autom. Control 2005, 50, 1354–1360. [Google Scholar] [CrossRef]
  19. Rynne, B.P.; Youngson, M.A. Linear Functional Analysis; Springer: London, UK, 2000. [Google Scholar]
  20. Garloff, J. The Bernstein algorithm. Interval Comput. 1993, 2, 154–168. [Google Scholar]
Figure 1. A closed loop system for Example 5.
Figure 1. A closed loop system for Example 5.
Mca 21 00025 g001

Share and Cite

MDPI and ACS Style

Büyükköroğlu, T. Fixed Order Controller for Schur Stability. Math. Comput. Appl. 2016, 21, 25. https://doi.org/10.3390/mca21020025

AMA Style

Büyükköroğlu T. Fixed Order Controller for Schur Stability. Mathematical and Computational Applications. 2016; 21(2):25. https://doi.org/10.3390/mca21020025

Chicago/Turabian Style

Büyükköroğlu, Taner. 2016. "Fixed Order Controller for Schur Stability" Mathematical and Computational Applications 21, no. 2: 25. https://doi.org/10.3390/mca21020025

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop