Design of Linear Time Varying Flatness-Based Control for Single-Input Single-Output Systems Marouen Sleimi Research Laboratory in Automatic Control - LARA University of Tunis EL Manar National Engineering School of Tunis, BP-36, le Bilvedere 002 Tunis Mohamed Ben Abdallah Higher Institute of Technology Study of Rades, BP 72, 2098, Radès Médina,Tunisia Mounir Ayadi Research Laboratory in Automatic Control - LARA University of Tunis EL Manar National Engineering School of Tunis, BP-36, le Bilvedere 002 Tunis Abstract In this paper, the control of linear discrete-time Varying Single-Input Single-Output systems is tacled By using flatness theory combined with a dead-beat observer, a two degree of freedom controller is designed with high performances in terms of trajectory tracing The aim of this wor is to avoid the choice of closed loop poles in linear discrete-time varying framewor which build a very serious problem in system control The effectiveness of this control law is highlighted by simulation results Keywords Flatness theory; discrete-time systems; linear time varying; single-input single-output; dead-beat observer; two degree of freedom controller I INTRODUCTION The theory of linear time-invariant systems gives a wide range of design methods and solutions to control problems including all inds of techniques such as state feedbac controllers and observers, Kalman filters, H 2 control and H control Since then, researchers have ept on maing efforts to extend the previous systems control approaches to timevarying systems and sampled-data systems LTV (Linear Time Varying) systems are of a great interest because of the fact that time invariant nonlinear systems can be approximated by LTV systems around desired trajectories after there linearization However, a fundamental part in the study of LTV systems is insured by the state transition matrix noted φ, which can be computed as the sum of the Peano- Baer series We can note here that not all arguments and assumptions used of time-invariant systems are useful in timevarying framewor To extend them, time-varying systems are examined carefully on their controllability and stability Moreover, many control design approaches use LTV systems, one of the most important way of systems control is SMC (Sliding mode control) which is very used in the case of dynamic uncertain systems Furthermore, we can evaluate Bac Stepping method which used to aim the objective of stabilizing nonlinear dynamic systems leading to an LTV system after linearization 2 Besides, we find fuzzy control operating nonlinear systems to mae nonlinear controllers via the use of heuristic information 3 Among these control design approaches, flatness-based control remains the most suitable method in trajectory tracing, then in the rest of the paper, we are interested to this ind of controller in the LTV case Previously, it is shown that flatness property considerably simplifies the 2DOF (Two Degree of Freedom) controllers design for continuous-time SISO (Single-Input Single-Output) systems for LTV framewor 4, 5 In these wors, the main feature of this flatness approach for 2DOF controllers design, using flatness-based control and dead-beat observer, is to avoid the choice of closed loop poles and no need to solve diophantine equation any more In this design the closed loop dynamics are related to the chosen tracing dynamics This approach was extended to discrete-time framewor for LTI SISO 6, 7 and LTI MIMO (Multiple-Input Multiple- Output) systems 8, 9, 0 This paper extend the previous approaches to deal with LTV flatness-based control in SISO discrete-time framewor This paper is organized as follows: In Section II, some preliminaries are presented Then in Section III, the definition of the canonical controllable form in discrete-time SISO LTV framewor is given Moreover, in Section IV, the new approach of control design is developed In Section V, effectiveness of this control law is illustrated using an academic discrete-time SISO system In the following, we will develop the paper in a discretetime formulation, using the shift forward operator q and the delay operator q II PRELIMINARIES Introducing a given problem concerning LTV systems, start generally by results given by LTI techniques and trying to adapt it in the new context Flatness-based control is introduced and developed by Fliess and co-researchers and used by many authors, firstly in the LTI framewor then in the LTV one This wor deals with flatness-based control in LTV discrete-time framewor So to aim this objective a canonical controllable form must be built and is exploited in the proposed control law design A Canonical controllable form for LTV discrete-time SISO systems Canonical forms are widely used in control theory In this section, the discrete-time LTV controllable canonical form is presented 2 598 P a g e
Considering the following discrete-time LTV system given by: x+ = A x + B u y = C x () where x is a n dimensional vector, A is a (n n) matrix, B and C are a m dimensional vectors If the given system is uniformly N-step controllable, then there exist a transformation applied to the state given by: Z = T x (2) where T is a Lyapunov transformation matrix, presented in Appendix A, obeying the conditions below T is defined each sample time T and T + are bounded each sample time There exist a constant m 0 < m < det( T ), > 0 This algorithm of controllable form leads to a new state vector equation written as following: Z+ = ĀZ + B ū y = C Z (3) where Ā, B and C are given by the new canonical controllable form for LTV systems, such as: Ā = T + A T (4) and u = f 2 ( z, z +,, z +m ) T (9) y = g ( x, u ) T (0) As nown, on the first hand, the state, the input vector and the output vector in continuous time framewor are defined as a successive derivatives of the flat output, on the other hand, in discrete-time framewor it is considered as a phase advance sequence of the flat output Marlait et al 3 results confirm that discrete-time LTV system is the equivalent of a controllable LTI system and noting that flatness of a given one is strongly related to the uniform controllability To build up flatness-based control it is necessary that the considered system is flat, so firstly its controllability 2, 4 must be checed out If the system is uniformly N-step controllable then we can apply the control law Let s consider the plant given by () and the transformation given by (2) Noting that z the first element of the state vector Z, then from (3) it s clear that Z can be written as : Z = z z + z +n T () where z is the Bronovsy output The expression of the input u is a function of z, and its forward terms: u = z +n + γ n ( + n )z +n + + γ ( + )z + (2) = 0 0 0 0 0 0 γ 0 () γ ( + ) γ n ( + n ) 0 B = T + B = 0 C = C T+ = ( α α (n ) ) B Flatness Properties Flatness can parametrize, in a very simple way, the dynamic aspect of a given system based on highlighting some endogenous fundamental variables called flat outputs In fact, the state, the input and the output using the flat output can be written as follows: (6) (7) x = f ( z, z +,, z +n ) T (8) (5) +γ 0 ()z The output y which depends on z is written as follows: y = C n Z = α i z +i (3) z can be considered as a flat output of the discrete LTV system From the canonical form given previously flatnessbased control for SISO LTV systems will be developed in the next section III FLATNESS-BASED CONTROL FOR SISO LTV SYSTEMS A Flatness-based Control Following () the flatness-based control law is given by: n n u = z+n d + κ i (z+i d z +i ) + γ i ( + i)z +i (4) κ i are derived from the following polynomial: n K(q) = q n + κ i q i (5) 599 P a g e
which must be shur 5, 6 To build this control it is necessary to estimate the state vector Z using an observer 7 In this paper the dead-beat observer is used in the design strategy leading to a 2DOF controller The structure of the flatness-based control is given by the following expression: u = K(q)z d + Λ Z (6) Λ = γ 0 κ 0 γ i(+i) κ i The control law schema is shown in Fig (7) Noting that: Y = ( y y +n ) T (23) U = ( u u +n 2 ) T (24) Equation (8) can be then written as: Y = L Z + M U (25) Supposing that the system is uniformly N-step observable then the ran of L is n for all Noting that L is the observability matrix, so: Z = L Y L M U (26) In a second h Z = q ( Ā Z + B u ) After a second iteration: Z = q ( A q (A Z + B u ) + B u ) (27) (28) Fig Control law design B Dead-beat Discrete-time Varying Observer Design Successive iterations of the output equation in (3) to the order (n ) 8, give: y u y + u + M = y +n L = = L Z + M C C (+) Ā C (+2) Ā + Ā u +n 2 C (+n ) Ā (+n 2) Ā M, 0 0 M 2, M,+ M 3, M 2,+ 0 M n, M n 2,+ M,+n (8) (9) (20) M, = C + B (2) M i, = C +i Ā +i B (22) The iteration of this expression to the order l gives: l Z = φ l q l Z + φ i B (i+) u (i+) + B u (29) i= φ l = l Ā j (30) j= φ 0 = I n (3) I is the identity matrix and φ l is (n n)-dimensional matrices By replacing Z in the right side of (29) by the expression in (26), Z can be written as follows: then: Z = φ l L Z = φ l q l (L Y L M U ) (32) φ i B (i+) u (i+) + B u l + l i= l + y l y φ l L l M l u l u +n 2 l φ i B (i+) u (i+) + B u i= As Z is written as: Z = φ l L l y l y φ l L l M l u l u +n 2 l (33) (34) 600 P a g e
+ ( φ l B l φ B 2 φ 0 B ) we obtain: u l u +n 2 l Z = φ l L l Y l φ l L l M l ΨU l (35) This system is an academic second order model, with single input and single output used to highlight the effectiveness of the discrete-time flatness-based control approach in LTV SISO case Noting that, is the iteration ran and T = 05s is the simple time The open loop step response is given by Fig 2 Supposing that l = n, the expression of Z is then given as follows: Z = φ (n ) L (n ) Y (n ) (36) φ l L (n ) M (n ) ΨU (n ) Ψ = Fig 2 φ (n ) B (n ) φ B 2 φ 0 B (37) Control law design C 2-DOF controller using flatness based control By replacing Z by its expression of (36) in the expression of the flatness-based control law given by (6) The flatness-based control can be then written in the following form: S(, q )u = K(q)z d R(, q)y (38) S(, q ) = ( + Λ φn L n M n Ψ )Π (39) R(, q) = Λ φ l L n Π (40) Π = ( q n+ ) T Π = ( q n+ q ) T (4) (42) The final form of the controller allows us to obtain a 2DOF controller in LTV framewor without need to define any observation dynamics and without resolution of the diophantine equation B Desired Trajectory To achieve the implementation of this control law, firstly the definition of the trajectory is needed to be traced by the considered system The aim is to force any system to trac the given trajectory, so it s necessary to choose the reference trajectory, then we can determine correctly the endogenous parameters satisfying this objective with a correct dynamics Let s consider the desired flat output defined as follows: z d (T ) = z in + ( 2 ( ) 5 ( ) 6 ( ) ) 7 T T T 35 + 5 (44) 50 50 50 (z fin z in ) z in and z fin are the initial and final values of the flat output that is determined endogenously from the initial and final values of the system outputs The trajectory presented in the discrete-time is used as a reference for the flatness-based control using dead-beat observer The discrete-time desired trajectory is represented in Fig 3 IV A Considered model ILLUSTRATIVE NUMERICAL EXAMPLE Let s consider the following discrete-time system defining by the matrices: B = 0 e T A = e T e (+)T, C = 0 (43) Fig 3 Desired trajectory for the discrete-time flat output 60 P a g e
C Control design From (4), the control law is given by: S(, q ) = ( s s 2 ) s 4 q (52) s 3 s 5 u = z d +2 + z d + + 0 z d + (γ 0 () 0 ) z (45) q R(, q) = r, r 2 (53) + (γ () ) z + After applying the controllable canonical form, the system matrices are given as follows: 0 Ā = (46) γ 0 () γ () B = e (+)T, C = 0 γ 0 () = e(+3) e (3) + e (3+) e (2+) + e (4+) e (47) γ () = f + f 2 f 3 (48) Λ = f4 f 5 f 6 f 7 f = e 3 e e e 2 e f 2 = e 2 e 4 + e 3 e 3 f 3 = e 2 e 2 e + e 4 e 4 f 4 = 00 ( e +3 e 3 + e 3+) + 7 ( e 2+ + e 4+ e ) f 5 = 00 ( e 2+ + e 4+ e ) f 6 = e+4 + e 3 + e 2+3 e +2 + e 3+4 e+ ( e 2 ) f 7 = e +2 + e 3+4 265 (49) s = e+ (e 2+ + ) (e +4 + e + + 20) e +2 + e 3+4 + 06 s 2 = e 2 e +3 + e 2 + e 2+ e 3+ e 2+3 + e 4+ e s 3 = e 2+ + e 4+ e s 4 = 00(e +3 + e 3+ ) + 3(e 2+ + e 4+ e ) 2 s 5 = 00(e 2+ + e 4+ e ) r = e+3 + 03e 2+ + e 3+ + 03e 4+ + 03e 20 e 2+ + e 4+ e r 2 = (r 3 r 4 r 5 ) r 6 r 3 = e 3+ (e 3+2 + e ) r 4 = e +5 e 6 + e 3+ r 5 = e+ (e 2 + ) e +4 + e 2+3 + 20 e +2 + e 3+4 + 06 r 6 = (e 5 e 2+5 + e 3+4 )(e 4 e +2 + e 2+2 ) A 2-DOF controller in LTV framewor is obtained without defining any observation dynamics, the effectiveness of this control method is shown in the next section D Simulation Results After applying the control law signal shown in Fig 4 to the considered system, the effectiveness of this method is proven and the system output follows the desired flat output with an error which tends to zero as shown in Fig 6 Both the system output and the desired flat output are represented in Fig 5 The transition matrix is: φ = Ā, φ 0 = Ψ = Ā B 2 B As noted previously, this control law can lead to a 2DOF controller with the following parameters: S (, q ) = ( + Λ φ L M Ψ ) Π (50) R (, q) = Λφ L Π (5) Fig 4 The control signal u 602 P a g e
Fig 5 System output and desired flat output: y and y d Fig 9 Tracing error for the second initial condition with error for the first one given by Fig 9 The linearization around a given trajectory applied to a nonlinear system leads to a linear time-varying system The real time simulation of such method poses a very important field of study because of the strong symbolic calculations leading to a large size of polynomial matrices Moreover, The robustness of such method touching the time-varying perturbations, the varying parameters and the practical implementation are prospects which will be carried out in the future wors Fig 6 Tracing error for the system output Fig 7 System output and desired flat output: y and y d with the first initial condition δy 0 = 02 V CONCLUSION In this paper, a discrete-time control law using the flatness property with trajectory tracing is proposed Many wors which have used the fatness theory to achieve a trajectory tracing by solving bezout equation and carrying out a pole placement Several difficulties gone through during the design of those methods Contrary, with a quite few calculations and avoiding the above-mentioned problems our developed approach seems more effective The design of this control law coupled with a dead-beat observer leads to a 2DOF controller without need to solve diophantine equation Such a controller based on flatness implies the use of a dead-beat observer, leading to a direct interpretation of the closed loop poles, which are naturally constituted by the poles of the the tracing polynomial K(q) This control law method shows good performances in terms of trajectory tracing in the case of LTV discrete-time systems APPENDIX The matrix T is called the transformation matrix, its the tool that transforms any system to its controllable form 2 In the LTV discrete-time case this matrix changes over time at each sample until convergence The plant () becomes: Z+ = T + A T Z + T + B U y = C T Z (54) Fig 8 System output and desired flat output: y and y d with the second initial condition δy 0 = 05 When applying two different initial conditions value to the desired trajectory, we obtain the results shown in Fig 7 and 8 with Ā = T + A T, B = T + B and C = C T In the following, we give the steps to obtain the controllable canonical form Consider the following vector sequence for i = n R 0 = B, i = 0 (55) R i = A R i ( + ), i = n 603 P a g e
the controllability matrix is then written as: R() = ( R 0 () R n () ) (56) if ran(r()) = n and to be uniformly N-step controllable means that, for all, the system is controllable on the interval N +, K + Then the effect of this variable change results in the following equation 2: then: R i () = T + R i () (57) R c () = T + R() (58) When the pair (Ā, B ) is deduced by the variables change of the pair (A, B ), and if R() is non-singular, then the following form is obtained 7: T + = R c ()R T () ( R()R T () ) otherwise we have the following form: (59) T + = R()R c () (60) The condition of singularity of the matrix R is a controllability criterion in the discrete-time case This property is preserved by the variables change 2 The algorithm for constructing the controllable form is then: Construction of γ() vector γ() = R ()R n () = ( γ 0 () γ n () ) T (6) Transformation T + then T which puts the pair (Ā, B ) as in (5), (6) and (7) REFERENCES Utin, VI, Variable structure systems wth sliding mode, IEEE Transaction on Automatic control, vol 22, n 2, pp22 222, 977 2 Carrol, JJ, Schneider, M, and Dawson, DM, Integrator bacstepping techniques for the traing control of permanent magnet bruch Dc motors, In conference Record of the 993 IEEE Industry Applications Society Annual Meeting, pp663 67, 993 3 Passino, KM and Yurovich, S, Fuyzzy control Addison-wesley Longman INC, 998 4 ROTELLA, F, CARRILLO, F J and AYADI, M, Polynomial controller design based on flatness, Kybernetia, vol 38, no 5, pp 57 584, 2002 5 BEN ABDALLAH, M, AYADI, M, ROTELLA, F et BENREJEB, M, Régulateurs polynomiaux par platitude pour la commande des systèmes non stationnaires, Conférence Internationale Francophone de l Automatique, CIFA 202, Grenoble 6 Ayadi, M, Contributions à la commande des systèmes linèaires plats de dimension finie, Thèse de Doctorat, Institut National Polytechnique de Toulouse, Tarbes 2002 7 Rotella, F, Carrillo, F and Ayadi, M, Digital flatness-based robust controller applied to a thermal process, IEEE International Conference on Control applications, Mexico, pp936 94, 200 8 BEN ABDALLAH, M, AYADI, M et BENREJEB, M, Flatness-based control of MIMO linear Systems, Transaction on Systems, Signals and Devices, vol 6, n, pp 24, 20 9 SLEIMI, M, BEN ABDALLAH, M and AYADI, M, Digital flatnessbased control design for LTI MIMO systems, IEEE 4th International Conference on control Engineering and Information Technology, CEIT, Hammamet, 206 0 BENABDALLAH, M, Sur la commande par platitude de systèmes dynamiques SISO et MIMO, Thèse de Doctorat en Génie Electrique, Ecole Nationale d Ingénieur de Tunis, 204 Fliess, M, Levine, J, Martin, P and Rouchon, P, Flatness and defect of non-linear systems: introductory and examples, International Journal of Control, vol 6, no 6, pp327 36, 995 2 Kamen, E W, Fundamentals of linear time-varying systems, The control handboo, Control System Advanced Methods, second edition, Taylor and Francis Group, New Yor 3 Malrait, F, Martin, Ph and Rouchon, P, Dynamic feedbac transformations of controllable linear time-varying systems, nonlinear Control in the year 2000 Edite par Isodori, A, Lamnabhi- Lagarrigue, F, Responde, W, Lecture Notes in Control and Information Sciences, vol 259, pp 55 62, Springer, London, 200 4 BINGULAC, S and MEADAWS, HE, Controllability ans observability in time-variable linear systems, SIMA Journal of control and optimization, no 5, pp64 73, 967 5 SUN, H, BUTT, SS and ASHERMAN, H, Discrete-time flatness-based control for a twin rotor helicopter with an extended alman ilter, IEEE International Conference on Advanced Intelligent Mechatronics, Banff, AB, Canada, 206 6 PAULO, S and ROUCHON, P, Flatnessbased control of a single qubit gate, IEEE transactions on Automatic control, vol 53, pp775 779, 2008 7 HUANG, R, Output feedbac tracing control of nonlinear time-varying systems by trajectory linearisation, Phd thesis, the Russ College of Engineering and Technology of Ohio University, 967 8 FIACCHINI, M and MILLERIOUX, G, Dead-Beat functional observers for discrete-time LPV systems with unnown inputs, IEEE transactions on Automatic control, vol 58, pp3230 3235, 203 604 P a g e