DESCRIPTOR FRACTIONAL LINEAR SYSTEMS WITH REGULAR PENCILS

Similar documents
APPROXIMATION OF FRACTIONAL POSITIVE STABLE CONTINUOUS TIME LINEAR SYSTEMS BY FRACTIONAL POSITIVE STABLE DISCRETE TIME SYSTEMS

SIMPLE CONDITIONS FOR PRACTICAL STABILITY OF POSITIVE FRACTIONAL DISCRETE TIME LINEAR SYSTEMS

Stabilization of fractional positive continuous-time linear systems with delays in sectors of left half complex plane by state-feedbacks

A new approach to the realization problem for fractional discrete-time linear systems

Analysis of the regularity, pointwise completeness and pointwise generacy of descriptor linear electrical circuits

REACHABILITY OF CONE FRACTIONAL CONTINUOUS TIME LINEAR SYSTEMS

Stability of interval positive continuous-time linear systems

POSITIVE PARTIAL REALIZATION PROBLEM FOR LINEAR DISCRETE TIME SYSTEMS

EXTERNALLY AND INTERNALLY POSITIVE TIME-VARYING LINEAR SYSTEMS

LOCALLY POSITIVE NONLINEAR SYSTEMS

CONTROL SYSTEMS, ROBOTICS AND AUTOMATION - Vol. VII - System Characteristics: Stability, Controllability, Observability - Jerzy Klamka

CONTROLLABILITY AND OBSERVABILITY OF 2-D SYSTEMS. Klamka J. Institute of Automatic Control, Technical University, Gliwice, Poland

Determining the order of minimal realization of descriptor systems without use of the Weierstrass canonical form

ALGEBRAIC CONDITION FOR DECOMPOSITION OF LARGE SCALE LINEAR DYNAMIC SYSTEMS

Minimum Energy Control of Fractional Positive Electrical Circuits with Bounded Inputs

On Positive Stable Realization for Continuous Linear Singular Systems

A 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.

Studies in Systems, Decision and Control. Series editor Janusz Kacprzyk, Polish Academy of Sciences, Warsaw, Poland

The Kalman-Yakubovich-Popov Lemma for Differential-Algebraic Equations with Applications

ANALYTIC SOLUTION OF TRANSCENDENTAL EQUATIONS

Structural Properties of LTI Singular Systems under Output Feedback

Full-State Feedback Design for a Multi-Input System

Reachability and controllability of time-variant discrete-time positive linear systems

Review of Controllability Results of Dynamical System

Control Systems. Frequency domain analysis. L. Lanari

Control Systems. Laplace domain analysis

1 Continuous-time Systems

PARAMETERIZATION OF STATE FEEDBACK GAINS FOR POLE PLACEMENT

Computers and Mathematics with Applications. The controllability of fractional control systems with control delay

CONTROL SYSTEMS, ROBOTICS AND AUTOMATION CONTENTS VOLUME VII

CONTROLLABILITY OF NONLINEAR DISCRETE SYSTEMS

On the Characterization Feedback of Positive LTI Continuous Singular Systems of Index 1

A Simple Derivation of Right Interactor for Tall Transfer Function Matrices and its Application to Inner-Outer Factorization Continuous-Time Case

Applied Mathematics Letters

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science : Dynamic Systems Spring 2011

Full State Feedback for State Space Approach

State Variable Analysis of Linear Dynamical Systems

POSITIVE REALNESS OF A TRANSFER FUNCTION NEITHER IMPLIES NOR IS IMPLIED BY THE EXTERNAL POSITIVITY OF THEIR ASSOCIATE REALIZATIONS

Peter C. Müller. Introduction. - and the x 2 - subsystems are called the slow and the fast subsystem, correspondingly, cf. (Dai, 1989).

On linear quadratic optimal control of linear time-varying singular systems

Intro. Computer Control Systems: F9

Solution for Homework 5

On the Interval Controllability of Fractional Systems with Control Delay

Lecture 2. Linear Systems

On the connection between discrete linear repetitive processes and 2-D discrete linear systems

Linear Systems. Linear systems?!? (Roughly) Systems which obey properties of superposition Input u(t) output

Model reduction via tangential interpolation

CONTROL SYSTEMS, ROBOTICS, AND AUTOMATION - Vol. V - Prediction Error Methods - Torsten Söderström

A dynamic model of mobile source air pollution and its properties

Analysis of undamped second order systems with dynamic feedback

1. Find the solution of the following uncontrolled linear system. 2 α 1 1

Finding Optimal Minors of Valuated Bimatroids

Some solutions of the written exam of January 27th, 2014

Module 07 Controllability and Controller Design of Dynamical LTI Systems

Research Article Stabilization Analysis and Synthesis of Discrete-Time Descriptor Markov Jump Systems with Partially Unknown Transition Probabilities

A BOUNDARY VALUE PROBLEM FOR LINEAR PDAEs

Linear Algebra and its Applications

Perspective. ECE 3640 Lecture 11 State-Space Analysis. To learn about state-space analysis for continuous and discrete-time. Objective: systems

EE C128 / ME C134 Final Exam Fall 2014

An Analytic-Numerical Method for Solving Difference Equations with Variable Coefficients by Discrete Multiplicative Integration

Stability of fractional positive nonlinear systems

Observability and Constructability

= AΦ, Φ(0) = I, = k! tk A k = I + ta t2 A t3 A t4 A 4 +, (3) e ta 1

ECEN 605 LINEAR SYSTEMS. Lecture 7 Solution of State Equations 1/77

EL2520 Control Theory and Practice

Constrained controllability of semilinear systems with delayed controls

On the computation of the Jordan canonical form of regular matrix polynomials

Design of Linear Functional Observer for MIMO LTI Systems

Chap 4. State-Space Solutions and

CONTROL SYSTEMS, ROBOTICS, AND AUTOMATION Vol. III Estimation and Compensation of Nonlinear Perturbations by Disturbance Observers - Peter C.

arxiv: v1 [math.ds] 27 Sep 2014

Research Article State-PID Feedback for Pole Placement of LTI Systems

EL2520 Control Theory and Practice

Control for Coordination of Linear Systems

Minimalsinglelinearfunctionalobserversforlinearsystems

L2 gains and system approximation quality 1

A Necessary and Sufficient Condition for High-Frequency Robustness of Non-Strictly-Proper Feedback Systems

LTI system response. Daniele Carnevale. Dipartimento di Ing. Civile ed Ing. Informatica (DICII), University of Rome Tor Vergata

Stability Theory for Nonnegative and Compartmental Dynamical Systems with Time Delay

ẋ n = f n (x 1,...,x n,u 1,...,u m ) (5) y 1 = g 1 (x 1,...,x n,u 1,...,u m ) (6) y p = g p (x 1,...,x n,u 1,...,u m ) (7)

Reachability and Controllability

Linear dynamical systems with inputs & outputs

Intro. Computer Control Systems: F8

Projection of state space realizations

Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design.

Linear System Theory

Time and Frequency domain design of Functional Filters

ACM/CMS 107 Linear Analysis & Applications Fall 2016 Assignment 4: Linear ODEs and Control Theory Due: 5th December 2016

Solutions of Generalized Linear Matrix Differential Equations which Satisfy Boundary Conditions at Two Points

Computation of a canonical form for linear differential-algebraic equations

On feedback stabilizability of time-delay systems in Banach spaces

Research Article Minor Prime Factorization for n-d Polynomial Matrices over Arbitrary Coefficient Field

First, we review some important facts on the location of eigenvalues of matrices.

9.2 The Input-Output Description of a System

Optimization Based Output Feedback Control Design in Descriptor Systems

Control Systems Design, SC4026. SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft

SAMPLE SOLUTION TO EXAM in MAS501 Control Systems 2 Autumn 2015

Tracking control for multi-agent consensus with an active leader and variable topology

ECE 388 Automatic Control

DETERMINATION OF AN UNKNOWN SOURCE TERM IN A SPACE-TIME FRACTIONAL DIFFUSION EQUATION

Transcription:

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 University, ul. Wiejska 45D, 5-35 Białystok, Pol e-mail: kaczorek@isep.pw.edu.pl Methods for finding solutions of the state equations of descriptor fractional discrete-time continuous-time linear systems with regular pencils are proposed. The derivation of the solution formulas is based on the application of the Z transform, the Laplace transform the convolution theorems. Procedures for computation of the transition matrices are proposed. The efficiency of the proposed methods is demonstrated on simple numerical examples. Keywords: descriptor system, fractional, system, regular pencil.. Introduction Descriptor (singular linear systems with regular pencils have been considered in many papers books (Dodig Stosic, 9; Wang,, Dai, 989; Fahmy O Reill, 989; Kaczorek, 4; 99; 7a; 7b; Kucera Zagalak, 988, Luenberger, 978; Van Dooren, 979. The eigenvalue invariant assignment by state output feedbacks was investigated by Dodig Stosic (9, Wang (, Dai (989, Fahmy O Reill (989 as well as Kaczorek (4; 99, while the realization problem for singular positive continuous-time systems with delays was discussed by Kaczorek (7b. The computation of Kronecker s canonical form of a singular pencil was analyzed by Van Dooren (979. A delay dependent criterion for a class of descriptor systems with delays varying in intervals was proposed by Wang (. Fractional positive continuous-time linear systems were addressed by Kaczorek (8, along with positive linear systems with different fractional order (Kaczorek, 7a. A new concept of practical stability of positive fractional D systems was proposed by Kaczorek (b, who also presented an analysis of fractional linear electrical circuits (Kaczorek, a some selected problems in the theory of fractional linear systems (Kaczorek, b. A new class of descriptor fractional linear systems electrical circuits was introduced, their solution of state equations was derived a method for decomposition of the descriptor fractional linear systems with regular pencils into dynamic static parts was proposed by Kaczorek (a, who also considered positive fractional continuous-time linear systems with singular pencils (Kaczorek, b. Fractional-order iterative learning control for fractional-order systems was addressed by Yan et al. (c. In this paper, methods of finding solutions of the state equations of descriptor fractional discrete-time continuous-time linear systems with regular pencils will be proposed. The paper is organized as follows. In Section the solution to the state equation of the descriptor system is derived using the method based on the Z transform the convolution theorem. A method for computation of the transition matrix is proposed illustrated on a simple numerical example in Section 3. In Section 4 the proposed method is extended to continuous-time linear systems. Concluding remarks are given in Section 5. The following notation will be used: R is the set of real numbers, R n m is the set of real n m matrices R n R n, Z is the set of nonnegative integers, I n is the n n identity matrix.. Discrete-time fractional linear systems Consider the descriptor fractional discrete-time linear system EΔ α x i Ax i Bu i, i Z {,,,...}, <α<, (

3 T. Kaczorek α is the fractional order, x i R n is the state vector u i R m is the input vector E,A R n n, B R n m. It is assumed that det E, but the pencil (E, A is regular, i.e., detez A for somez C. ( Without lost of generality we may assume that E E R n n, E R r r rank E ranke r<n. (3 Admissible boundary conditions for ( are given by x. The fractional difference of the order α, is defined by i Δ α x i c k x i k, (4a ( α k { c k ( k ( α k α(α...(α k k!, k,,... (4b for for k, k,,... (4c Substitution of (4a into ( yields i Ex i Fx i c k x i k Bu i, i Z, (5 k F A Ec A Eα. Applying to (5 the Z-transform taking into account that (Kaczorek, 99 Zx i p p z p X(zz p x j z j, p,,..., (6 j we obtain X(z Ex F {Ex z H(zBU(z}, X(z Zx i x i z i, U(z Zu i i u i z i, i i H(z Zh i, h i Ec k x i k. k (7a (7b Let Ex F j μ ψ j z (j, (8 μ is the positive integer defined by the pair (E,A (Kaczorek, 99; Van Dooren, 979. Comparison of the coefficients at the same powers of z of the equality Ex F j μ j μ ψ j z (j ψ j z (j Ex F I n (9a yields Eψ μ ψ μ E (9b Eψ k Eψ k ψ k E ψ k E { In for k, for k μ, μ,...,,,,... (9c From (9b (9c we have the matrix equation ψμ V G, (a ψ N G G R G G (Nμn (Nμn,... F... G...... RNn (μn,... E......... G............... R (μn (μn, G...... E.................. R Nn Nn,

Descriptor fractional linear systems with regular pencils ψ μ ψ μ ψ μ.. R(μn n, ψ N V ψ ψ ψ. ψ N. I n RNn n, R(μn n. Equation (a has the solution ψμ ψ N (b for given G V if only if { } V rank G, rankg. ( It is easy to show that the condition ( is satisfied if the condition ( is met. Substituting (8 into (7a we obtain X(z ψ j z (j Ex z H(zBU(z. j μ ( Applying the inverse transform Z the convolution theorem to (, we obtain iμ k x i ψ i Ex ψ i k c j x k j iμ j ψ i k Bu k. (3 To find the solution to Eqn. (, first we compute the transition matrices ψ j for j μ, μ,...,,,... next, using (3, we obtain the desired solution. 3. Computation of transition matrices To compute the transition matrices ψ k for k μ, μ,...,n,..., the following procedure is recommended. Procedure. Step. Find a solution ψ μ of the equation G ψ μ V, (4 G, ψ μ V are defined by (b. Note that, if the matrix E has the form (3, then the first r rows of the matrix ψ μ are zero its n r last rows are arbitrary. 3 Step. Choose n r arbitrary rows of the matrix ψ so that { E rank the equation E } In Fψ, E rank ψ ψ In Fψ (5 has a solution with arbitrary last n r rows of the matrix ψ. Step 3. Knowing ψ μ, choose the last n r rows of the matrix ψ so that { } E Fψ E rank, rank (6 the equation E ψ F ψ ψ (7 has a solution with arbitrary last n r rows of the matrix ψ. Repeating the last step for ψ ψ3,,... ψ 3 ψ 4 we may compute the desired matrices ψ k for k μ, μ,... The details of the procedure will be shown on the following example. Example. Find the solution to Eqn. ( for α.5 with the matrices E, A, B (8 the initial condition x u i, i Z. In this case the pencil ( of (8 is regular since detez A z z, (9 μ α A α.5. (

3 T. Kaczorek Using Procedure, we obtain the following. Step. In this case Eqn. (4 has the form ψ α ψ, ( its solution with the arbitrary second row ψ ψ of ψ is given by ψ ψ.5 ψ ψ. ( Step. We choose the row ψ ψ of ψ so that (5 holds, i.e., rank α the equation α has the solution ψ rank ψ ψ ψ α.5 α ψ ψ, (3 (4 (5 with the second arbitrary row ψ ψ of ψ. Step 3. We choose ψ ψ so that the equation α α ψ ψ α α ψ ψ (6 has the solution ψ ψ α.5α α ψ ψ (7 with arbitrary ψ ψ. Continuing the procedure, we obtain ψ, ψ.5 k ( k α k.5α k (8 for k,,... Using (3, (8 (, we obtain the desired solution of the form x i ψ i i i ψ i k k ψ i k j c j are defined by (4b. c j x k j u k, (9 4. Continuous-time fractional linear systems Consider the descriptor fractional continuous-time linear system described by the state equations ED α x(t Ax(tBu(t, n <α n {,,...} y(t Cx(tDu(t, (3a (3b D α is the Caputo differentiation operator, x(t R n,u(t R m, y(t R p are the state, input output vectors E,A R n n,b R n m,c R p n,d R p m. It is assumed that det E the pencil (E,A is regular, i.e., deteλ A for some λ C. (3 Admissible initial conditions for (3a are given by x (k ( x k, for k,,...,n. (3 Applying the Laplace transform (L to Eqn. (3a, we obtain (Kaczorek, b Es α AX(s BU(s s α k x k,, (33 X(s Lx(t U(s Lu(t. Ifthe condition (3 is satisfied, then from (33 we obtain X(s Es α A (BU(s s α k x k,. (34

Descriptor fractional linear systems with regular pencils In a particular case when det E, from (34 we have X(s (E A i E s (iα (35 since i ( BU(s s α k x k, Es α A Es α (I n (Es α A (E A i E s (iα. i Using the inverse Laplace transform (L the convolution theorem, we obtain x(t L X(s t (E i (t τ(iα A E Bu(τdτ Γ(i α i t iαk Γ(iα k (E A i x k,, (36 Γ(α is the gamma function (Kaczorek, b. Therefore, the following theorem for det E has been proved. Theorem. The solution of Eqn. (3a for det E the initial conditions (3 is given by (36. If E I n, then (36 takes the form (Kaczorek, b x(t n Φ l (tx (l ( l Φ l (t t A k t (kαl Γ(kα l, Φ(t τbu(τdτ, (37 (38a A k t (kα Φ(t Γ(k α. (38b If det E butthepencil is regular (the condition (3 is met, then Es α A T i satisfy the equality T i s (iα, (39 i μ ET i AT i T i E T i A { In for i, for i, (4 33 T i for i< μ, T μ E ET μ. The equality (4 follows from the comparison of the coefficients at the same powers of s in the equality Es α A T i s (iα i μ i μ Substitution of (39 into (34 yields X(s i μ i T i s (iα Es α A I n. T i s (iα BU(s s iα k x k, T i s (iα BU(s s iα k x k, i μ T i s (iα BU(s s iα k x k,. (4 Applying the inverse Laplace transform the convolution theorem to (4, we obtain x(t t (t τ(iα T i B Bu(τdτ Γ(i α i T i μ i i t iαk Γ(iα k Ex k, T i Bu(t (i α δ (iα Ex k, (4 or t x(t (T A i (t τ(iα T B u(τdτ i Γ(i α (T A i t iαk T Γ(iα k x k, μ T i Bu(t (i α δ (iα Ex k,, (43 since by (4 T i (T A i T for i,,... δ (k is the k-th derivative of the delta impulse function δ. Therefore, the following theorem has been proved. Theorem. If the condition (3 is satisfied, then the solution of Eqn. (3a with the admissible initial conditions (3 is given by (4 or (43.

34 T. Kaczorek In a particular case of <α, from (43 we have x(t t (T A i (t τ(iα T B u(τdτ i Γ(i α t iα Γ(iα Ex μ T i Bu(t (i α δ (iα Ex. i (44 To compute the matrices T i for i μ, μ,...,the procedure given in Section 3 is recommended. Example. Consider Eqn. (3a for α.5with the matrices E, A, B (45 the zero initial condition x. The pencil is regular since deteλ A λ λ, (λ sα (46 Eλ A T λ.5λ.5.5 T T λ,, T.5 (47a. (47b Using (4, (43 (47b we obtain t (t τ.5 x(t T B u(τdτ T Bu(t Γ(.5 t (t τ.5 u(τdτ u(t..5 Γ(.5 (48 5. Concluding remarks New methods of finding solutions of the state equations of descriptor fractional discrete-time continuous-time linear systems with regular pencils have been proposed. Derivation of the solution formulas has been based on the application of the Z-transform, the Laplace transform the convolution theorems. A procedure for computation of the transition matrices has been proposed its application has been demonstrated on simple numerical examples. An open problem is the extension of the method for D descriptor fractional discrete continuous-discrete linear systems. Acknowledgment This work was supported by the National Science Centre in Pol under the grant no. N N54 6389 4. References Dodig, M. Stosic, M. (9. Singular systems state feedbacks problems, Linear Algebra Its Applications 43(8: 67 9. Dai, L. (989. Singular Control Systems, Lectures Notes in Control Information Sciences, Vol. 8, Springer-Verlag, Berlin. Fahmy, M.H O Reill, J. (989. Matrix pencil of closed-loop descriptor systems: Infinite-eigenvalues assignment, International Journal of Control 49(4: 4 43. Gantmacher, F.R. (96. The Theory of Matrices, Chelsea Publishing Co., New York, NY. Kaczorek, T. (a. Descriptor fractional linear systems with regular pencils, Asian Journal of Control 5(4: 4. Kaczorek, T. (b. Positive fractional continuous-time linear systems with singular pencils, Bulletin of the Polish Academy of Sciences: Technical Sciences 6(: 9. Kaczorek, T. (a. Positive linear systems consisting of n subsystems with different fractional orders, IEEE Transactions on Circuits Systems 58(7: 3. Kaczorek, T. (b. Selected Problems of Fractional System Theory, Springer-Verlag, Berlin. Kaczorek, T. (a. Positive linear systems with different fractional orders, Bulletin of the Polish Academy of Sciences: Technical Sciences 58(3: 453 458. Kaczorek, T. (b. Practical stability asymptotic stability of positive fractional D linear systems, Asian Journal of Control (: 7. Kaczorek, T. (8. Fractional positive continuous-time linear systems their reachability, International Journal of Applied Mathematics Computer Science 8(: 3 8, DOI:.478/v6-8--. Kaczorek, T. (7a. Polynomial Rational Matrices. Applications in Dynamical Systems Theory, Springer-Verlag, London. Kaczorek, T. (7b. Realization problem for singular positive continuous-time systems with delays, Control Cybernetics 36(: 47 57. Kaczorek, T. (4. Infinite eigenvalue assignment by an output feedbacks for singular systems, International Journal of Applied Mathematics Computer Science 4(: 9 3.

Descriptor fractional linear systems with regular pencils Kaczorek, T. (99. Linear Control Systems, Vol., Research Studies Press J. Wiley, New York, NY. Kucera, V. Zagalak, P. (988. Fundamental theorem of state feedback for singular systems, Automatica 4(5: 653 658. Luenberger, D.G. (978. Time-invariant descriptor systems, Automatica 4: 473 48. Podlubny, I. (999. Fractional Differential Equations, Academic Press, New York, NY. Wang, C. (. New delay-dependent stability criteria for descriptor systems with interval time delay, Asian Journal of Control 4(: 97 6. Van Dooren, P. (979. The computation of Kronecker s canonical form of a singular pencil, Linear Algebra Its Applications 7: 3 4. Yan L., YangQuan C., Hyo-Sung A., (c. Fractional-order iterative learning control for fractional-order systems, Asian Journal of Control 3(: 54 63. 35 Tadeusz Kaczorek received the M.Sc., Ph.D. D.Sc. degrees in electrical engineering from the Warsaw University of Technology in 956, 96 964, respectively. In the years 968 69 he was the dean of the Electrical Engineering Faculty, in the period of 97 73 he was a deputy rector of the Warsaw University of Technology. In 97 he became a professor in 974 a full professor at the same university. Since 3 he has been a professor at Białystok Technical University. In 986 he was elected a corresponding member in 996 a full member of the Polish Academy of Sciences. In the years 988 99 he was the director of the Research Centre of the Polish Academy of Sciences in Rome. In 4 he was elected an honorary member of the Hungarian Academy of Sciences. He has been granted honorary doctorates by nine universities. His research interests cover systems theory, especially singular multidimensional systems, positive multidimensional systems, singular positive D D systems, as well as positive fractional D D systems. He initiated research in the field of singular D, positive D positive fractional linear systems. He has published 4 books (six in English over scientific papers. He has also supervised 69 Ph.D. theses. He is the editor-in-chief of the Bulletin of the Polish Academy of Sciences: Technical Sciences a member of editorial boards of ten international journals. Received: 8 May Revised: September