POSITIVE PARTIAL REALIZATION PROBLEM FOR LINEAR DISCRETE TIME SYSTEMS
|
|
- Angelina Shields
- 5 years ago
- Views:
Transcription
1 Int J Appl Math Comput Sci 7 Vol 7 No 65 7 DOI: 478/v POSITIVE PARTIAL REALIZATION PROBLEM FOR LINEAR DISCRETE TIME SYSTEMS TADEUSZ KACZOREK Faculty of Electrical Engineering Białystok Technical University Wiejska 45D 5 35 Białystok kaczorek@iseppwedupl A partial realization problem for positive linear discrete-time systems is addressed Sufficient conditions for the existence of its solution are established A procedure for the computation of a positive partial realization for a given finite sequence of the values of the impulse response is proposed The procedure is illustrated by four numerical examples Keywords: positive systems partial realization existence computation linear systems discrete-time systems Introduction In positive systems the inputs state variables outputs take only non-negative values Examples of positive systems are industrial processes involving chemical reactors heat exchangers distillation columns storage systems compartmental systems or water atmospheric pollution models A variety of models having positive linear system behaviour can be found in engineering management science economics social sciences biology medicine etc Positive linear systems are defined on cones not on linear spaces Therefore the theory of positive systems is more complicated less advanced Overviews of the state of the art in positive systems theory are given in the monographs (Farina Rinaldi ; Kaczorek ) The realization problem for positive linear systems was considered in many papers books (Benvenuti Farina 4; Farina Rinaldi ; Kaczorek 99; Kaczorek ; Kaczorek 6a; Kaczorek 6b; Kaczorek 4; Kaczorek 5; Kaczorek 6c; Kaczorek 6d; Kaczorek 6e; Kaczorek Busłowicz 4) The problem of constructing linear state variable models from given impulse responses (or Markov parameters) was considered in (Ho Kalman 966; Kaczorek 99) The partial realization problem consists in constructing state variables models from given finite sequences of impulse responses In this paper the partial realization problem will be addressed for positive linear discrete-time systems Sufficient conditions for the existence of its solution will be established a procedure for the computation of its positive realization for a given finite sequence of the values of the impulse response will be proposed To the best of the author s knowledge the positive partial realization problem for discrete-time linear system has not been considered yet Preliminaries Problem Formulation Let R n m be the set of n m real matrices R n : R n Thesetofn m real matrices with non-negative entries will be denoted by R n m + the set of nonnegative integers by Z + Consider the linear discrete-time system x i+ Ax i + Bu i i Z + y i Cx i + Du i () where x i R n u i R m y i R p are the state input output vectors A R n n B R n m C R p n D R p m The system () is called (internally) positive if only if x i R n + y i R p + i Z + for every x R n + any input sequence u i R m + i Z + The system () is (internally) positive if only if (Kaczorek ; Fa-
2 66 T Kaczorek rina Rinaldi ): A R n n + B Rn m + C R p n + D Rp m + () The transfer matrix of () is given by T (z)ciz A B+D CI Ad db+d (3) where d z is the delay operator Let G(d) G i d i (4) be the impulse response matrix of the system () From the well-known equality { D for i G i CA i (5) B for i () it follows that for the positive system () we have that G i R p m + for i Z + The partial positive realization problem can be formulated as follows: Given a finite sequence of nonnegative matrices Ĝ i R p m + for i N (6) Subproblem Given a finite sequence ĝ i for i N () find n N the transfer function (8) such that ĝ i g i for i N () where g i is defined by (9) Subproblem Given the transfer function (8) find its positive realization () 3 Solution of Subproblem Taking into account that T (d) using (8) we can write g i d i () a(d)g(d) b(d) (3) find the matrices () of the positive system () such that Ĝ i G i for i N (7) a(d) ĝ(d)+ḡ(d) b(d) (4) with G i defined by (4) In this paper sufficient conditions for the existence of a solution to the above problem will be established a procedure for the computation of the matrices () for a given finite sequence G i R p m + N will be proposed 3 Problem Solution First the idea of the proposed method will be outlined for a single-input single-output (SISO m p )systemin this case the transfer function of () can be written down in the form T (d) b(d) a(d) its impulse response is g(d) b i d i n (8) a i d i g i d i g i R + i Z + (9) The partial positive realization problem can be decomposed into the following two subproblems: where the polynomial is known the sum ĝ(d) ḡ(d) N ĝ i d i (5) in+ ḡ i d i (6) is unknown Equation (4) can be rewritten in the form where a(d)ĝ(d) b(d)+ g(d) (7) g(d) a(d)ḡ(d) in+ g i d i (8) The subproblem has thus been reduced to the following one: Given the polynomial (5) find the polynomials a(d) b(d) of (8) such that (7) holds for some g(d)
3 Positive partial realization problem for linear discrete-time systems 67 Using (8) (5) we obtain a(d)ĝ(d) where ( j a i d i)( N ĝ j d j) j N i+n ĝ j d j a i ĝ k i d k N k + ( ĝ k N+n kn+ N q k d k + k ki q k ĝ k ( n a i ĝ k i ) d k ik N N+n kn+ a i ĝ k i ) d k ( n ik N a i ĝ k i )d k (9) ĝ k i a i () Let q k for k n +N Then from () we have Ĝ nn a ĝ nn () where ĝ n ĝ n ĝ ĝ n ĝ Ĝ nn ĝ N ĝ N ĝ N n a a n We assume that ĝ nn ĝ n+ ĝ N () rank ĜnN rankĝnn ĝ nn (3) Eqn () has a non-negative solution ie a R n +In (Kaczorek 6e) a method was proposed for the computation of a positive solution to () Note that (9) has the form (7) if only if q k for k n +N or equivalently the coefficients a n of the polynomial a(d) are a solution to () Knowing the coefficients a i i n we can compute the coefficients b j j nof the polynomial b(d) using the formula b k ĝ k ĝ k i a i for k n (4) It follows from the comparison of the right-h sides of (7) (9) for q k k n +N In summary we have the following result: Theorem Let n N Given a finite sequence () there exists a transfer function of the form T (z) b z n + + b n z + b n z n z n a n z a n D + b z n + + b n z + b n z n z n a n z a n (5) with non-negative coefficients if Eqn () has a nonnegative solution ie a R n + k bk ĝ k ĝ k i a i k n (6) 3 Solution of Subproblem It is well known (Kaczorek ; Benvenuti Farina 4; Kaczorek 4; Kaczorek 6a; Kaczorek 6c) that if the coefficients a n b b b n are nonnegative then there always exists a positive realization () of the transfer function (5) For example we may choose the positive realization as follows: A B a n a n a n C bn bn b D b (7a) or a n bn a n bn Ā a n B b C D b (7b) From (5) we have D b b b + b b b + b b n b n + a n b (8) using (4) we obtain (6) Therefore if the coefficients a n b b b n are non-negative then so are the coefficients b b b n the realizations (7) are then positive From (8) it follows that b k a k do not imply that b k kn
4 68 T Kaczorek If () has a non-negative solution a R n +thena positive realization () for a given finite sequence (9) can be computed using the following procedure: Procedure (SISO systems) Step Given a sequence () find a non-negative solution a R n + to () Step Using (4) compute the coefficients b b b n find the transfer function (5) Step 3 Using one of the well-known methods (Benvenuti Farina 4; Farina Rinaldi ; Kaczorek ; Kaczorek 6a; Kaczorek 6b; Kaczorek 4; Kaczorek 5; Kaczorek 6c; Kaczorek 6d; Kaczorek 6e; Kaczorek Busłowicz 4) compute a desired positive realization () of the transfer function (5) eg the positive realization (7) Remark If n is not known apriori then it is recommended to start the procedure with its small value to increase it in the next step Find a positive realization () for the se- Example quence ĝ ĝ 3 ĝ 8 ĝ 3 9 ĝ 4 46 (9) Using the procedure for N 4 n we obtain the following: Step In this case Eqn () for (9) has the form its solution is Step Using (4) we compute the coefficients b ĝ b ĝ ĝ b ĝ ĝ ĝ The desired transfer function has the form T (z) z + z + 3z + z + z z (3) z 3z + z z 3z +8z +9z 3 +46z 4 + Step 3 Using (7a) for (3) we obtain the desired positive realization A B C b b 3 D (3) Now let us assume that only the first four numbers of the sequence (9) are given ie N 3 n Inthis case Eqn () takes the form (3) One of the coefficients can be chosen arbitrarily If we choose then we obtain the same transfer function (3) its positive realization (3) If we choose then /3 b b b 4/3 The corresponding transfer function has the form T (z) 3z +6z +4 3z +5 3z + 3z z z 3 its positive realization is given by In this case 3z +5 z z 3 A 3 B (33) C 5 3 D (34) 3z +8z +9z 3 +48z 4 + Note that g 4 48is different from ĝ 4 46in the previous case Theorem Let a finite sequence ĝ i > for i n (35) be given Then a positive realization of the form (7) exists if ĝ n > ĝn+ ĝ n+ ĝ n ĝ n+ > ĝn+ ĝ n+3 ĝ n ĝ n > > > ĝ(n ) ĝ n > > ĝ(n ) ĝ n ĝ > > ĝn ĝ (n ) > ĝn ĝ n > ĝn > ĝn ĝ n (36) k bk ĝ k ĝ k i a i k n (37) with a i i n constituting the solution to () for N n
5 Positive partial realization problem for linear discrete-time systems 69 Proof Given a sequence (35) Eqn () takes the form ĝ n ĝ n ĝ ĝ ĝ n ĝ 3 ĝ ĝ (n ) ĝ n 3 ĝ n ĝ n a n ĝ n ĝ (n ) ĝ n a n ĝ n+ ĝ n ĝ n Therefore ĝ n ĝ n ĝ ĝ ĝ n ĝ 3 ĝ ĝ n+ ĝ n+ ĝ n+ ĝ n+ ĝ (n ) ĝ n 3 ĝ n ĝ n ĝ n ĝ n ĝ n ĝ n ĝ n ĝ (n ) ĝ n ĝ n ĝ n ĝ n ĝ n (38) a n a n If the conditions (36) are met then the coefficient matrix of (38) has a cyclic structure by Theorem A (see the Appendix) Eqn (38) has a positive solution a i > for i n In this case using (37) we may find the coefficients b k for k n the positive realization (7) Example (Continuation of Example ) For the sequence (9) we have that N n 4 (38) takes the form (39) It is easy to see that the coefficient matrix of (39) satisfies the conditions (36) since 8 9 > > 3 9 The transfer function has the form (3) the desired positive realization is given by (3) Theorem 3 A positive asymptotically stable realization (7) of the finite sequence (35) exists if ĝ i+ < ĝ i for i n (4) the conditions (36) (37) satisfied Proof If the conditions (36) are met then by Theorem A Eqn (38) has a positive solution a i > for i n We shall show that if (4) holds then a i < (4) Let S i j g n i+j g n+j be the sum of the entries of the i-th column of the coefficient matrix in (38) i n Adding the n equations of (38) we obtain S i a i n (4) Note that if (4) holds then each entry g n i+j / g n+j of the coefficient matrix is greater than Write S mins i Then from (4) we have i S a i n a i n S < since S >n From the well-known final value theorem we obtain g lim g i lim (z ) T (z) i z lim (z ) C I n z A B+D z b z n ++b z+b lim (z ) z z n z n a n z a since (4) holds Therefore the positive realization (7) is asymptotically stable Example 3 We wish to find a positive realization () for the sequence g g g 9 g 3 6 g 4 5 (43)
6 7 T Kaczorek This sequence satisfies the conditions (4) (36) (37) for n by Theorem 3 there exists a positive asymptotically stable realization Using Procedure we obtain the following: Step We have the system of equations g g 9 g 3 g g 3 g 4 g g 4 its solution is 4/ 9/ Step Using (37) we compute the coefficients b ĝ b ĝ ĝ 49 The desired transfer function has the form z + 49 z +49 T (z) z 4 z 9 z 4z 9 (44) Step 3 Using (7a) for (44) we obtain the asymptotically stable realization A 9 4 B 49 C b b D 4 Extension for Multi Input Multi Output Systems The method presented in Section 3 can be easily extended for linear discrete-time systems with m inputs p outputs (MIMO systems) In this case the modified procedure has the following form: Procedure (MIMO systems) Step Given a sequence (6) find a non-negative solution to () with Ĝ nn a I m Ĝ nn I m a n I m T Ĝ n Ĝ n Ĝ Ĝ n+ Ĝ n Ĝ Ĝ N Ĝ N Ĝ N n Ĝ n+ Ĝ n+ (45) Ĝ N Step Using B k Ĝk Ĝ k i a i for k n (46) compute the matrices B B B n the transfer matrix T (z) z n z n a n z a n B z n +B z n ++B n z+b n (47) Step 3 Using one of the well-known methods (Benvenuti Farina 4; Farina Rinaldi ; Kaczorek ; Kaczorek 6a; Kaczorek 6b; Kaczorek 4; Kaczorek 5; Kaczorek 6c; Kaczorek 6d; Kaczorek 6e; Kaczorek Busłowicz 4) compute the desired positive realization () of the transfer matrix (47) Example 4 We wish to find a positive realization () for the sequence Ĝ Ĝ 3 7 Ĝ Ĝ Ĝ 5 7 (48) We have m p Using Procedure for N 4 n we obtain the following: Step In this case Eqn () with (45) has the form its solution is Step Using (46) we compute the matrices B Ĝ ˆB Ĝ Ĝ ˆB Ĝ Ĝ Ĝ The desired transfer matrix (47) has the form T (z) z z z z z + z + z + z + z z +4 (49) (5)
7 Positive partial realization problem for linear discrete-time systems 7 Step 3 Using (7a) we obtain the desired positive realization () of (5) in the form A B C Concluding Remarks D (5) The partial realization problem for positive discrete-time linear systems has been formulated solved It was shown that it has a solution if Eqn () has a non-negative solution (6) holds (Theorem ) A procedure for SISO MIMO systems was proposed for the computation of positive partial realizations for a given finite sequence () (or (6)) The procedure was illustrated with four numerical examples It was shown that if the conditions (4) (36) (38) are satisfied then there exists a positive asymptotically stable realization of the form (7) The deliberations can be extended for linear continuous-time systems for D linear systems (Kaczorek ) Acknowledgment This work was supported by the Ministry of Science Higher Education in Pol under Grant No 3 TA 6 7 References Benvenuti L Farina L (4): A tutorial on the positive realization problem IEEE Trans Automat Contr Vol 49 No 5 pp Farina L Rinaldi S (): Positive Linear Systems: Theory Applications Wiley New York Ho BL Kalman RE (966): Effective construction of linear state-variable models from input/output functions Regelungstechnik Bd pp Kaczorek T (99): Linear Control Systems Vol Taunton: Research Studies Press Kaczorek T (): Positive D D Systems London: Springer Kaczorek T (4): Realization problem for positive discretetime systems with delay System Sci Vol 3 No 4 pp 7 3 Kaczorek T (5): Positive minimal realizations for singular discrete-time systems with delays in state control Bull Pol Acad Sci Techn Sci Vol 53 No 3 pp Kaczorek T (6a): Realization problem for positive multivariable discrete-time linear systems with delays in the state vector inputs Int J Appl Math Comput Sci Vol 6 No pp Kaczorek T (6b) A realization problem for positive continuous-time linear systems with reduced number of delay Int J Appl Math Comput Sci Vol 6 No 3 pp Kaczorek T (6c): Computation of realizations of discretetime cone-systems Bull Polish Acad Sci Techn Sci Vol 54 No 3 pp Kaczorek T (6d): A realization problem for positive continuous-time systems with reduced numbers of delays Int J Appl Math Comput Sci Vol 6 No 3 pp 7 Kaczorek T (6e): Minimal positive realization for discretetime systems with state time-delays Int J Comput Math Electr Eng COMPEL Vol 5 No 4 pp 8 86 Kaczorek T Busłowicz M (4): Minimal realization problem for positive multivariable linear systems with delay Int J Appl Math Comput Sci Vol 4 No pp 8 87 Merzyakov JI (963): On the existence of positive solutions of a system of linear equations Uspekhi Matematicheskikh Nauk Vol 8 No 3 pp (in Russian) Appendix Definition A A positive matrix n A a ij > a n a nn i n; j n (A) is called a matrix with a (strictly) cyclic structure if only if a ii >a i+i > >a ni >i > >a i i (A) for i n For example the matrix has a cyclic structure the matrix does not have a cyclic structure since <a 3 3 Theorem A (Merzyakov 963) The system of linear algebraic equations a i x + a i x + + a in x i n (A3) with a ij > ijnhas a positive solution x i > nif only if its coefficient matrix (A) has a cyclic structure Received: March 7
8
A new approach to the realization problem for fractional discrete-time linear systems
BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES Vol 64 No 26 DOI: 55/bpasts-26-2 A new approach to the realization problem for fractional discrete-time linear systems T KACZOREK Faculty of
More informationAPPROXIMATION OF FRACTIONAL POSITIVE STABLE CONTINUOUS TIME LINEAR SYSTEMS BY FRACTIONAL POSITIVE STABLE DISCRETE TIME SYSTEMS
Int. J. Appl. Math. Comput. Sci., 213, Vol. 23, No. 3, 51 56 DOI: 1.2478/amcs-213-38 APPROXIMATION OF FRACTIONAL POSITIVE STABLE CONTINUOUS TIME LINEAR SYSTEMS BY FRACTIONAL POSITIVE STABLE DISCRETE TIME
More informationStability of interval positive continuous-time linear systems
BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES, Vol. 66, No. 1, 2018 DOI: 10.24425/119056 Stability of interval positive continuous-time linear systems T. KACZOREK Białystok University of
More informationStabilization of fractional positive continuous-time linear systems with delays in sectors of left half complex plane by state-feedbacks
Control and Cybernetics vol. 39 (2010) No. 3 Stabilization of fractional positive continuous-time linear systems with delays in sectors of left half complex plane by state-feedbacks by Tadeusz Kaczorek
More informationAnalysis of the regularity, pointwise completeness and pointwise generacy of descriptor linear electrical circuits
Computer Applications in Electrical Engineering Vol. 4 Analysis o the regularity pointwise completeness pointwise generacy o descriptor linear electrical circuits Tadeusz Kaczorek Białystok University
More informationSIMPLE CONDITIONS FOR PRACTICAL STABILITY OF POSITIVE FRACTIONAL DISCRETE TIME LINEAR SYSTEMS
Int. J. Appl. Math. Comput. Sci., 2009, Vol. 19, No. 2, 263 269 DOI: 10.2478/v10006-009-0022-6 SIMPLE CONDITIONS FOR PRACTICAL STABILITY OF POSITIVE FRACTIONAL DISCRETE TIME LINEAR SYSTEMS MIKOŁAJ BUSŁOWICZ,
More informationLOCALLY POSITIVE NONLINEAR SYSTEMS
Int. J. Appl. Math. Comput. Sci. 3 Vol. 3 No. 4 55 59 LOCALLY POSITIVE NONLINEAR SYSTEMS TADEUSZ KACZOREK Institute of Control Industrial Electronics Warsaw University of Technology ul. Koszykowa 75 66
More informationEXTERNALLY AND INTERNALLY POSITIVE TIME-VARYING LINEAR SYSTEMS
Int. J. Appl. Math. Comput. Sci., 1, Vol.11, No.4, 957 964 EXTERNALLY AND INTERNALLY POSITIVE TIME-VARYING LINEAR SYSTEMS Tadeusz KACZOREK The notions of externally and internally positive time-varying
More informationOn Positive Stable Realization for Continuous Linear Singular Systems
Int. Journal of Math. Analysis, Vol. 8, 2014, no. 8, 395-400 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2014.4246 On Positive Stable Realization for Continuous Linear Singular Systems
More informationDESCRIPTOR FRACTIONAL LINEAR SYSTEMS WITH REGULAR PENCILS
Int. J. Appl. Math. Comput. Sci., 3, Vol. 3, No., 39 35 DOI:.478/amcs-3-3 DESCRIPTOR FRACTIONAL LINEAR SYSTEMS WITH REGULAR PENCILS TADEUSZ KACZOREK Faculty of Electrical Engineering Białystok Technical
More informationarxiv: v1 [math.ds] 27 Sep 2014
AUTOMATYZACJA PROCESÓW DYSKRETNYCH 2014 Konrad Andrzej MARKOWSKI Warsaw University of Technology, Electrical Department Institute of Control and Industrial Electronics Koszykowa 75, 00-662 Warsaw e-mail:
More informationReachability and controllability of time-variant discrete-time positive linear systems
Control and Cybernetics vol. 33 (2004) No. 1 Reachability and controllability of time-variant discrete-time positive linear systems by Ventsi G. Rumchev 1 and Jaime Adeane 2 1 Department of Mathematics
More informationA basic canonical form of discrete-time compartmental systems
Journal s Title, Vol x, 200x, no xx, xxx - xxx A basic canonical form of discrete-time compartmental systems Rafael Bru, Rafael Cantó, Beatriz Ricarte Institut de Matemàtica Multidisciplinar Universitat
More informationMinimal positive realizations of transfer functions with nonnegative multiple poles
1 Minimal positive realizations of transfer functions with nonnegative multiple poles Béla Nagy Máté Matolcsi Béla Nagy is Professor at the Mathematics Department of Technical University, Budapest, e-mail:
More informationCONTROLLABILITY OF NONLINEAR DISCRETE SYSTEMS
Int. J. Appl. Math. Comput. Sci., 2002, Vol.2, No.2, 73 80 CONTROLLABILITY OF NONLINEAR DISCRETE SYSTEMS JERZY KLAMKA Institute of Automatic Control, Silesian University of Technology ul. Akademicka 6,
More informationREACHABILITY OF CONE FRACTIONAL CONTINUOUS TIME LINEAR SYSTEMS
Int. J. Appl. Math. Comput. Sci., 29, Vol. 19, No. 1, 89 93 DOI: 1.2478/v16-9-8-4 REACHABILITY OF CONE FRACTIONAL CONTINUOUS TIME LINEAR SYSTEMS TADEUSZ KACZOREK Faculty of Electrical Engineering Białystok
More informationMinimum Energy Control of Fractional Positive Electrical Circuits with Bounded Inputs
Circuits Syst Signal Process (26) 35:85 829 DOI.7/s34-5-8-7 Minimum Energy Control o Fractional Positive Electrical Circuits with Bounded Inputs Tadeusz Kaczorek Received: 24 April 25 / Revised: 5 October
More informationA dynamic model of mobile source air pollution and its properties
Environmental Health Risk IV 4 A dynamic model of mobile source air pollution and its properties V. Rumchev & L. Caccetta Department of Mathematics and Statistics, Curtin University of Technology, Australia
More informationCANONICAL LOSSLESS STATE-SPACE SYSTEMS: STAIRCASE FORMS AND THE SCHUR ALGORITHM
CANONICAL LOSSLESS STATE-SPACE SYSTEMS: STAIRCASE FORMS AND THE SCHUR ALGORITHM Ralf L.M. Peeters Bernard Hanzon Martine Olivi Dept. Mathematics, Universiteit Maastricht, P.O. Box 616, 6200 MD Maastricht,
More informationPositive systems in the behavioral approach: main issues and recent results
Positive systems in the behavioral approach: main issues and recent results Maria Elena Valcher Dipartimento di Elettronica ed Informatica Università dipadova via Gradenigo 6A 35131 Padova Italy Abstract
More informationCONTROLLABILITY AND OBSERVABILITY OF 2-D SYSTEMS. Klamka J. Institute of Automatic Control, Technical University, Gliwice, Poland
CONTROLLABILITY AND OBSERVABILITY OF 2D SYSTEMS Institute of Automatic Control, Technical University, Gliwice, Poland Keywords: Controllability, observability, linear systems, discrete systems. Contents.
More informationANALYTIC SOLUTION OF TRANSCENDENTAL EQUATIONS
Int J Appl Math Comput Sci 00 ol 0 No 4 67 677 DOI: 0478/v0006-00-0050- ANALYTIC SOLUTION OF TRANSCENDENTAL EQUATIONS HENRYK GÓRECKI Faculty of Informatics Higher School of Informatics ul Rzgowska 7a 9
More informationASYMPTOTIC BEHAVIOR OF SOLUTIONS OF DISCRETE VOLTERRA EQUATIONS. Janusz Migda and Małgorzata Migda
Opuscula Math. 36, no. 2 (2016), 265 278 http://dx.doi.org/10.7494/opmath.2016.36.2.265 Opuscula Mathematica ASYMPTOTIC BEHAVIOR OF SOLUTIONS OF DISCRETE VOLTERRA EQUATIONS Janusz Migda and Małgorzata
More informationStability of fractional positive nonlinear systems
Archives of Cotrol Scieces Volume 5(LXI), 15 No. 4, pages 491 496 Stability of fractioal positive oliear systems TADEUSZ KACZOREK The coditios for positivity ad stability of a class of fractioal oliear
More informationPositive observers for positive interval linear discrete-time delay systems. Title. Li, P; Lam, J; Shu, Z
Title Positive observers for positive interval linear discrete-time delay systems Author(s) Li, P; Lam, J; Shu, Z Citation The 48th IEEE Conference on Decision and Control held jointly with the 28th Chinese
More informationIMPULSIVE CONTROL OF DISCRETE-TIME NETWORKED SYSTEMS WITH COMMUNICATION DELAYS. Shumei Mu, Tianguang Chu, and Long Wang
IMPULSIVE CONTROL OF DISCRETE-TIME NETWORKED SYSTEMS WITH COMMUNICATION DELAYS Shumei Mu Tianguang Chu and Long Wang Intelligent Control Laboratory Center for Systems and Control Department of Mechanics
More informationOn the Characterization Feedback of Positive LTI Continuous Singular Systems of Index 1
Advanced Studies in Theoretical Physics Vol. 8, 2014, no. 27, 1185-1190 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/astp.2014.48245 On the Characterization Feedback of Positive LTI Continuous
More informationCONTINUITY OF SOLUTIONS OF RICCATI EQUATIONS FOR THE DISCRETE TIME JLQP
Int. J. Appl. Math. Comput. Sci., 2002, Vol.12, No.4, 539 543 CONTINUITY O SOLUTIONS O RICCATI EQUATIONS OR THE DISCRETE TIME JLQP ADAM CZORNIK, ANDRZEJ ŚWIERNIAK Department of Automatic Control Silesian
More informationCONTROL SYSTEMS, ROBOTICS AND AUTOMATION - Vol. VII - System Characteristics: Stability, Controllability, Observability - Jerzy Klamka
SYSTEM CHARACTERISTICS: STABILITY, CONTROLLABILITY, OBSERVABILITY Jerzy Klamka Institute of Automatic Control, Technical University, Gliwice, Poland Keywords: stability, controllability, observability,
More informationBUMPLESS SWITCHING CONTROLLERS. William A. Wolovich and Alan B. Arehart 1. December 27, Abstract
BUMPLESS SWITCHING CONTROLLERS William A. Wolovich and Alan B. Arehart 1 December 7, 1995 Abstract This paper outlines the design of bumpless switching controllers that can be used to stabilize MIMO plants
More informationNewton, Fermat, and Exactly Realizable Sequences
1 2 3 47 6 23 11 Journal of Integer Sequences, Vol. 8 (2005), Article 05.1.2 Newton, Fermat, and Exactly Realizable Sequences Bau-Sen Du Institute of Mathematics Academia Sinica Taipei 115 TAIWAN mabsdu@sinica.edu.tw
More informationSpectral Properties of Matrix Polynomials in the Max Algebra
Spectral Properties of Matrix Polynomials in the Max Algebra Buket Benek Gursoy 1,1, Oliver Mason a,1, a Hamilton Institute, National University of Ireland, Maynooth Maynooth, Co Kildare, Ireland Abstract
More informationCanonical lossless state-space systems: staircase forms and the Schur algorithm
Canonical lossless state-space systems: staircase forms and the Schur algorithm Ralf L.M. Peeters Bernard Hanzon Martine Olivi Dept. Mathematics School of Mathematical Sciences Projet APICS Universiteit
More informationPOSITIVE REALNESS OF A TRANSFER FUNCTION NEITHER IMPLIES NOR IS IMPLIED BY THE EXTERNAL POSITIVITY OF THEIR ASSOCIATE REALIZATIONS
POSITIVE REALNESS OF A TRANSFER FUNCTION NEITHER IMPLIES NOR IS IMPLIED BY THE EXTERNAL POSITIVITY OF THEIR ASSOCIATE REALIZATIONS Abstract This letter discusses the differences in-between positive realness
More informationStability Theory for Nonnegative and Compartmental Dynamical Systems with Time Delay
1 Stability Theory for Nonnegative and Compartmental Dynamical Systems with Time Delay Wassim M. Haddad and VijaySekhar Chellaboina School of Aerospace Engineering, Georgia Institute of Technology, Atlanta,
More informationPositive and Monotone Systems
History of Control, 2012, Lund May 29, 2012 Outline. 1 Positive Systems Definition Fathers Occurrence Example Publications. 2 Monotone Systems Definition Early days Publications Positive System A continuous
More informationOn the connection between discrete linear repetitive processes and 2-D discrete linear systems
Multidim Syst Sign Process (217) 28:341 351 DOI 1.17/s1145-16-454-8 On the connection between discrete linear repetitive processes and 2-D discrete linear systems M. S. Boudellioua 1 K. Galkowski 2 E.
More informationStudies in Systems, Decision and Control. Series editor Janusz Kacprzyk, Polish Academy of Sciences, Warsaw, Poland
Studies in Systems, Decision and Control Volume 13 Series editor Janusz Kacprzyk, Polish Academy of Sciences, Warsaw, Poland e-mail: kacprzyk@ibspan.waw.pl About this Series The series "Studies in Systems,
More informationA q x k+q + A q 1 x k+q A 0 x k = 0 (1.1) where k = 0, 1, 2,..., N q, or equivalently. A(σ)x k = 0, k = 0, 1, 2,..., N q (1.
A SPECTRAL CHARACTERIZATION OF THE BEHAVIOR OF DISCRETE TIME AR-REPRESENTATIONS OVER A FINITE TIME INTERVAL E.N.Antoniou, A.I.G.Vardulakis, N.P.Karampetakis Aristotle University of Thessaloniki Faculty
More informationIN many practical systems, there is such a kind of systems
L 1 -induced Performance Analysis and Sparse Controller Synthesis for Interval Positive Systems Xiaoming Chen, James Lam, Ping Li, and Zhan Shu Abstract This paper is concerned with the design of L 1 -
More informationApplication of LMI for design of digital control systems
BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES Vol. 56, No. 4, 28 Application of LMI for design of digital control systems A. KRÓLIKOWSKI and D. HORLA Institute of Control and Information
More informationAn alternative proof of the Barker, Berman, Plemmons (BBP) result on diagonal stability and extensions - Corrected Version
1 An alternative proof of the Barker, Berman, Plemmons (BBP) result on diagonal stability and extensions - Corrected Version Robert Shorten 1, Oliver Mason 1 and Christopher King 2 Abstract The original
More informationControl for stability and Positivity of 2-D linear discrete-time systems
Manuscript received Nov. 2, 27; revised Dec. 2, 27 Control for stability and Positivity of 2-D linear discrete-time systems MOHAMMED ALFIDI and ABDELAZIZ HMAMED LESSI, Département de Physique Faculté des
More informationZero controllability in discrete-time structured systems
1 Zero controllability in discrete-time structured systems Jacob van der Woude arxiv:173.8394v1 [math.oc] 24 Mar 217 Abstract In this paper we consider complex dynamical networks modeled by means of state
More informationStructural identifiability of linear mamillary compartmental systems
entrum voor Wiskunde en Informatica REPORTRAPPORT Structural identifiability of linear mamillary compartmental systems JM van den Hof Department of Operations Reasearch, Statistics, and System Theory S-R952
More informationA Simple Derivation of Right Interactor for Tall Transfer Function Matrices and its Application to Inner-Outer Factorization Continuous-Time Case
A Simple Derivation of Right Interactor for Tall Transfer Function Matrices and its Application to Inner-Outer Factorization Continuous-Time Case ATARU KASE Osaka Institute of Technology Department of
More informationReachability of a class of discrete-time positive switched systems
Reachability of a class of discrete-time positive switched systems Ettore Fornasini and Maria Elena Valcher Abstract In this paper, monomial reachability and reachability of a special class of discrete-time
More informationThe Tuning of Robust Controllers for Stable Systems in the CD-Algebra: The Case of Sinusoidal and Polynomial Signals
The Tuning of Robust Controllers for Stable Systems in the CD-Algebra: The Case of Sinusoidal and Polynomial Signals Timo Hämäläinen Seppo Pohjolainen Tampere University of Technology Department of Mathematics
More informationPositive Markov Jump Linear Systems (PMJLS) with applications
Positive Markov Jump Linear Systems (PMJLS) with applications P. Bolzern, P. Colaneri DEIB, Politecnico di Milano - Italy December 12, 2015 Summary Positive Markov Jump Linear Systems Mean stability Input-output
More informationProcess Modelling, Identification, and Control
Jan Mikles Miroslav Fikar Process Modelling, Identification, and Control With 187 Figures and 13 Tables 4u Springer Contents 1 Introduction 1 1.1 Topics in Process Control 1 1.2 An Example of Process Control
More informationarxiv:math/ v1 [math.co] 13 Jul 2005
A NEW PROOF OF THE REFINED ALTERNATING SIGN MATRIX THEOREM arxiv:math/0507270v1 [math.co] 13 Jul 2005 Ilse Fischer Fakultät für Mathematik, Universität Wien Nordbergstrasse 15, A-1090 Wien, Austria E-mail:
More informationMonomial subdigraphs of reachable and controllable positive discrete-time systems
Monomial subdigraphs of reachable and controllable positive discrete-time systems Rafael Bru a, Louis Caccetta b and Ventsi G. Rumchev b a Dept. de Matemàtica Aplicada, Univ. Politècnica de València, Camí
More informationLinear Algebra and its Applications
Linear Algebra and its Applications 433 (200) 867 875 Contents lists available at ScienceDirect Linear Algebra and its Applications journal homepage: www.elsevier.com/locate/laa On the exponential exponents
More informationApplications of Controlled Invariance to the l 1 Optimal Control Problem
Applications of Controlled Invariance to the l 1 Optimal Control Problem Carlos E.T. Dórea and Jean-Claude Hennet LAAS-CNRS 7, Ave. du Colonel Roche, 31077 Toulouse Cédex 4, FRANCE Phone : (+33) 61 33
More informationG. NADIBAIDZE. a n. n=1. n = 1 k=1
GEORGIAN MATHEMATICAL JOURNAL: Vol. 6, No., 999, 83-9 ON THE ABSOLUTE SUMMABILITY OF SERIES WITH RESPECT TO BLOCK-ORTHONORMAL SYSTEMS G. NADIBAIDZE Abstract. Theorems determining Weyl s multipliers for
More informationH 2 Suboptimal Estimation and Control for Nonnegative
Proceedings of the 2007 American Control Conference Marriott Marquis Hotel at Times Square New York City, USA, July 11-13, 2007 FrC20.3 H 2 Suboptimal Estimation and Control for Nonnegative Dynamical Systems
More informationAssignment 10. Arfken Show that Stirling s formula is an asymptotic expansion. The remainder term is. B 2n 2n(2n 1) x1 2n.
Assignment Arfken 5.. Show that Stirling s formula is an asymptotic expansion. The remainder term is R N (x nn+ for some N. The condition for an asymptotic series, lim x xn R N lim x nn+ B n n(n x n B
More informationTuning regulators for tracking MIMO positive LTI systems
Tuning regulators for tracking MIMO positive LTI systems Bartek Roszak Edward J. Davison Edward S. Rogers Sr. Department of Electrical and Computer Engineering University of Toronto Systems Control Group
More information1 Multiply Eq. E i by λ 0: (λe i ) (E i ) 2 Multiply Eq. E j by λ and add to Eq. E i : (E i + λe j ) (E i )
Direct Methods for Linear Systems Chapter Direct Methods for Solving Linear Systems Per-Olof Persson persson@berkeleyedu Department of Mathematics University of California, Berkeley Math 18A Numerical
More informationMinimal-span bases, linear system theory, and the invariant factor theorem
Minimal-span bases, linear system theory, and the invariant factor theorem G. David Forney, Jr. MIT Cambridge MA 02139 USA DIMACS Workshop on Algebraic Coding Theory and Information Theory DIMACS Center,
More informationMultivariable ARMA Systems Making a Polynomial Matrix Proper
Technical Report TR2009/240 Multivariable ARMA Systems Making a Polynomial Matrix Proper Andrew P. Papliński Clayton School of Information Technology Monash University, Clayton 3800, Australia Andrew.Paplinski@infotech.monash.edu.au
More informationControl Systems. Laplace domain analysis
Control Systems Laplace domain analysis L. Lanari outline introduce the Laplace unilateral transform define its properties show its advantages in turning ODEs to algebraic equations define an Input/Output
More informationController synthesis for positive systems under l 1-induced performance
Title Controller synthesis for positive systems under l 1-induced performance Author(s) Chen, X; Lam, J; Li, P; Shu, Z Citation The 24th Chinese Control and Decision Conference (CCDC 212), Taiyuan, China,
More informationA Tropical Extremal Problem with Nonlinear Objective Function and Linear Inequality Constraints
A Tropical Extremal Problem with Nonlinear Objective Function and Linear Inequality Constraints NIKOLAI KRIVULIN Faculty of Mathematics and Mechanics St. Petersburg State University 28 Universitetsky Ave.,
More informationMATH 3511 Lecture 1. Solving Linear Systems 1
MATH 3511 Lecture 1 Solving Linear Systems 1 Dmitriy Leykekhman Spring 2012 Goals Review of basic linear algebra Solution of simple linear systems Gaussian elimination D Leykekhman - MATH 3511 Introduction
More informationALGEBRAIC CONDITION FOR DECOMPOSITION OF LARGE SCALE LINEAR DYNAMIC SYSTEMS
Int. J. Appl. Math. Comput. Sci., 2009, Vol. 19, No. 1, 107 111 DOI: 10.2478/v10006-009-0010-x ALGEBRAIC CONDITION FOR DECOMPOSITION OF LARGE SCALE LINEAR DYNAMIC SYSTEMS HENRYK GÓRECKI Faculty of Informatics
More informationON THE COMPUTATION OF THE MINIMAL POLYNOMIAL OF A POLYNOMIAL MATRIX
Int J Appl Math Comput Sci, 25, Vol 5, No 3, 339 349 ON THE COMPUTATION OF THE MINIMAL POLYNOMIAL OF A POLYNOMIAL MATRIX NICHOLAS P KARAMPETAKIS, PANAGIOTIS TZEKIS Department of Mathematics, Aristotle
More informationProblems of Eigenvalues/Eigenvectors
67 Problems of Eigenvalues/Eigenvectors Reveiw of Eigenvalues and Eigenvectors Gerschgorin s Disk Theorem Power and Inverse Power Methods Jacobi Transform for Symmetric Matrices Spectrum Decomposition
More informationThe Miracles of Tropical Spectral Theory
The Miracles of Tropical Spectral Theory Emmanuel Tsukerman University of California, Berkeley e.tsukerman@berkeley.edu July 31, 2015 1 / 110 Motivation One of the most important results in max-algebra
More informationObservations on the Stability Properties of Cooperative Systems
1 Observations on the Stability Properties of Cooperative Systems Oliver Mason and Mark Verwoerd Abstract We extend two fundamental properties of positive linear time-invariant (LTI) systems to homogeneous
More informationPOLYNOMIAL SOLUTIONS OF PELL S EQUATION AND FUNDAMENTAL UNITS IN REAL QUADRATIC FIELDS
J. London Math. Soc. 67 (2003) 16 28 C 2003 London Mathematical Society DOI: 10.1112/S002461070200371X POLYNOMIAL SOLUTIONS OF PELL S EQUATION AND FUNDAMENTAL UNITS IN REAL QUADRATIC FIELDS J. MCLAUGHLIN
More informationModel reduction of interconnected systems
Model reduction of interconnected systems A Vandendorpe and P Van Dooren 1 Introduction Large scale linear systems are often composed of subsystems that interconnect to each other Instead of reducing the
More informationFALSENESS OF THE FINITENESS PROPERTY OF THE SPECTRAL SUBRADIUS
Int J Appl Math Comput Sci, 2007, Vol 17, No 2, 173 178 DOI: 102478/v10006-007-0016-1 FALSENESS OF THE FINITENESS PROPERTY OF THE SPECTRAL SUBRADIUS ADAM CZORNIK, PIOTR JURGAŚ Department of Automatic Control,
More informationTHE MEASURE-THEORETIC ENTROPY OF LINEAR CELLULAR AUTOMATA WITH RESPECT TO A MARKOV MEASURE. 1. Introduction
THE MEASURE-THEORETIC ENTROPY OF LINEAR CELLULAR AUTOMATA WITH RESPECT TO A MARKOV MEASURE HASAN AKIN Abstract. The purpose of this short paper is to compute the measure-theoretic entropy of the onedimensional
More information9.5 The Transfer Function
Lecture Notes on Control Systems/D. Ghose/2012 0 9.5 The Transfer Function Consider the n-th order linear, time-invariant dynamical system. dy a 0 y + a 1 dt + a d 2 y 2 dt + + a d n y 2 n dt b du 0u +
More information1 Continuous-time Systems
Observability Completely controllable systems can be restructured by means of state feedback to have many desirable properties. But what if the state is not available for feedback? What if only the output
More informationCONTROL SYSTEMS, ROBOTICS AND AUTOMATION CONTENTS VOLUME VII
CONTENTS VOLUME VII Control of Linear Multivariable Systems 1 Katsuhisa Furuta,Tokyo Denki University, School of Science and Engineering, Ishizaka, Hatoyama, Saitama, Japan 1. Linear Multivariable Systems
More informationSums of diagonalizable matrices
Linear Algebra and its Applications 315 (2000) 1 23 www.elsevier.com/locate/laa Sums of diagonalizable matrices J.D. Botha Department of Mathematics University of South Africa P.O. Box 392 Pretoria 0003
More informationTrace Inequalities for a Block Hadamard Product
Filomat 32:1 2018), 285 292 https://doiorg/102298/fil1801285p Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://wwwpmfniacrs/filomat Trace Inequalities for
More informationOPTIMAL INPUT SIGNAL DESIGN FOR IDENTIFICATION OF MAX PLUS LINEAR SYSTEMS
OPTIMAL INPUT SIGNAL DESIGN FOR IDENTIFICATION OF MAX PLUS LINEAR SYSTEMS Gernot Schullerus, Volker Krebs, Bart De Schutter, Ton van den Boom Institut für Regelungs- und Steuerungssysteme, Universität
More informationReview of Vectors and Matrices
A P P E N D I X D Review of Vectors and Matrices D. VECTORS D.. Definition of a Vector Let p, p, Á, p n be any n real numbers and P an ordered set of these real numbers that is, P = p, p, Á, p n Then P
More informationControl Systems. Frequency domain analysis. L. Lanari
Control Systems m i l e r p r a in r e v y n is o Frequency domain analysis L. Lanari outline introduce the Laplace unilateral transform define its properties show its advantages in turning ODEs to algebraic
More informationConstrained controllability of semilinear systems with delayed controls
BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES Vol. 56, No. 4, 28 Constrained controllability of semilinear systems with delayed controls J. KLAMKA Institute of Control Engineering, Silesian
More informationThe PRIMES 2014 problem set
Dear PRIMES applicant! The PRIMES 24 problem set This is the PRIMES 24 problem set. Please send us your solutions as part of your PRIMES application by December, 24. For complete rules, see http://web.mit.edu/primes/apply.shtml
More informationMULTI-RESTRAINED STIRLING NUMBERS
MULTI-RESTRAINED STIRLING NUMBERS JI YOUNG CHOI DEPARTMENT OF MATHEMATICS SHIPPENSBURG UNIVERSITY SHIPPENSBURG, PA 17257, U.S.A. Abstract. Given positive integers n, k, and m, the (n, k)-th m- restrained
More informationThe characteristic equation and minimal state space realization of SISO systems in the max algebra
KULeuven Department of Electrical Engineering (ESAT) SISTA Technical report 93-57a The characteristic equation and minimal state space realization of SISO systems in the max algebra B De Schutter and B
More informationQuadratic and Copositive Lyapunov Functions and the Stability of Positive Switched Linear Systems
Proceedings of the 2007 American Control Conference Marriott Marquis Hotel at Times Square New York City, USA, July 11-13, 2007 WeA20.1 Quadratic and Copositive Lyapunov Functions and the Stability of
More informationComplements to Full Order Observers Design for Linear Systems with Unknown Inputs
Author manuscript, published in "Applied Mathematics Letters 22, 7 (29) 117-1111" DOI : 1.116/j.aml.28.11.4 omplements to Full Order Observers Design for Linear Systems with Unknown Inputs M. Darouach
More informationMinimal Positive Realizations of Transfer Functions with Positive Real Poles
1370 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS I: FUNDAMENTAL THEORY AND APPLICATIONS, VOL 47, NO 9, SEPTEMBER 2000 Minimal Positive Realizations of Transfer Functions with Positive Real Poles Luca Benvenuti,
More informationLifted approach to ILC/Repetitive Control
Lifted approach to ILC/Repetitive Control Okko H. Bosgra Maarten Steinbuch TUD Delft Centre for Systems and Control TU/e Control System Technology Dutch Institute of Systems and Control DISC winter semester
More informationUsing Matrix Analysis to Approach the Simulated Annealing Algorithm
Using Matrix Analysis to Approach the Simulated Annealing Algorithm Dennis I. Merino, Edgar N. Reyes Southeastern Louisiana University Hammond, LA 7040 and Carl Steidley Texas A&M University- Corpus Christi
More informationOUTPUT REGULATION OF RÖSSLER PROTOTYPE-4 CHAOTIC SYSTEM BY STATE FEEDBACK CONTROL
International Journal in Foundations of Computer Science & Technology (IJFCST),Vol., No., March 01 OUTPUT REGULATION OF RÖSSLER PROTOTYPE-4 CHAOTIC SYSTEM BY STATE FEEDBACK CONTROL Sundarapandian Vaidyanathan
More informationKalman Decomposition B 2. z = T 1 x, where C = ( C. z + u (7) T 1, and. where B = T, and
Kalman Decomposition Controllable / uncontrollable decomposition Suppose that the controllability matrix C R n n of a system has rank n 1
More informationInt. J. Contemp. Math. Sciences, Vol. 3, 2008, no. 35, Fermat-GCD Matrices
Int. J. Contemp. Math. Sciences, Vol. 3, 2008, no. 35, 745-753 Fermat-GCD Matrices Şerife Büyükköse The University of AhiEvran, Faculty of Science and Arts Department of Mathematics, 4000 Kirşehir, Turkey
More informationCommentationes Mathematicae Universitatis Carolinae
Commentationes Mathematicae Universitatis Carolinae Bogdan Szal Trigonometric approximation by Nörlund type means in L p -norm Commentationes Mathematicae Universitatis Carolinae, Vol. 5 (29), No. 4, 575--589
More informationON THE CONSTRUCTION OF GENERAL SOLUTION OF THE GENERALIZED SYLVESTER EQUATION
TWMS J App Eng Math V7, N1, 2017, pp 1-6 ON THE CONSTRUCTION OF GENERAL SOLUTION OF THE GENERALIZED SYLVESTER EQUATION FA ALIEV 1, VB LARIN 2, Abstract The problem of construction the general solution
More informationSubquadratic Computational Complexity Schemes for Extended Binary Field Multiplication Using Optimal Normal Bases
1 Subquadratic Computational Complexity Schemes for Extended Binary Field Multiplication Using Optimal Normal Bases H. Fan and M. A. Hasan March 31, 2007 Abstract Based on a recently proposed Toeplitz
More informationResearch Article Stabilization Analysis and Synthesis of Discrete-Time Descriptor Markov Jump Systems with Partially Unknown Transition Probabilities
Research Journal of Applied Sciences, Engineering and Technology 7(4): 728-734, 214 DOI:1.1926/rjaset.7.39 ISSN: 24-7459; e-issn: 24-7467 214 Maxwell Scientific Publication Corp. Submitted: February 25,
More informationPersistence and global stability in discrete models of Lotka Volterra type
J. Math. Anal. Appl. 330 2007 24 33 www.elsevier.com/locate/jmaa Persistence global stability in discrete models of Lotka Volterra type Yoshiaki Muroya 1 Department of Mathematical Sciences, Waseda University,
More informationOn the Regularized Trace of a Fourth Order. Regular Differential Equation
Int. J. Contemp. Math. Sci., Vol., 26, no. 6, 245-254 On the Regularized Trace of a Fourth Order Regular Differential Equation Azad BAYRAMOV, Zerrin OER, Serpil ÖZTÜRK USLU and Seda KIZILBUDAK C. ALIṢKAN
More information