MEM 355 Performance Enhancement of Dynamical Systems MIMO Introduction
|
|
- Emma Walton
- 6 years ago
- Views:
Transcription
1 MEM 355 Performance Enhancement of Dynamical Systems MIMO Introduction Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 11/2/214
2 Outline Solving State Equations Variation of parameters formula Matrix exponential Controllability/Observability/Kalman Decomposition System Zeros SISO/MIMO Closed Loops
3 Solving Linear State Equations n x = Ax + Bu t, x R, u R find : x t given : x t = x, u t for t t for t t x = Ax + b t, b t : = Bu t, forced or nonhomogeneous x = Ax homogeneous m
4 Approach Basic Properties Variation of Parameters Formula State Transition Matrix Matrix Exponential Laplace Transform
5 Basic properties p x t, x t sol'ns of homog., c, c constants = + x t cx t cx t = 1 2 any sol'n of forced, is a sol'n of homog. x t, x t sol'ns of forced x t x t x t x t x t is a sol'n of homog. any sol'n of homog. x t = x t + x t is a sol'n of forced h p h
6 Variation of Parameters Formula Recall, any sol'n of (forced) satisfies where = + x t x t x t h p h = p ( t) is any (particular) sol'n of (forced) xp( t) xp( t) = X ( t) c( t) = + and 1 = x c ( t) = X ( t) B( t) u( t) x t X t c for constant vector c, satisfies (homog.) x We seek. Assume the form. x Ax Bu X A p p = c( t ) c t + t t 1 X B u d τ τ τ τ
7 Variation of Parameters Formula, Cont d h = x t Ax t assume ( ) ( k k = ) h 2 = a + 2at+ 3at + = A a + at+ at a = Aa, a = Aa, a = Aa,, a = Aa k k k 1 a = Aa, a = A a, a = A a,, a = A a h x t a at at k k k! x t a Aat Aat Aat 2 k! = ( + ) k k xh t I At 2 A t k! A t ( ) e I At A t A t At k k 2 k! initial condition = = At a c xh t e c a
8 Variation of Parameters Formula, Cont d We seek x t. Assume the form x t = e c t. p At d At x = Ax + Bu and e = Ae p p dt + = + At At At Ae c t e c t Ae c t B t u t c t = e B t u t c t At = c + t At Aτ e B u d t At ( τ ) x p t = e B u d p τ τ τ τ τ τ
9 Variations of Parameters Formula Cont d Thus, we have x t = x t + x t h t At At ( τ ) x t = e c + e Bu ( τ) d Now impose: x = x c= x p t At At ( τ ) x t = e x + e Bu τ d τ τ
10 Transition Matrix via Laplace Transform d e At At At = Ae, e = I dt t= L = L At At s e I A e ( ) At 1 1 e si A e si A L = = [ ] At 1, L [ ]
11 Caley-Hamilton Theorem Cayley-Hamilton Theorem: Every square matrix satisfies its own characteristic equation, i.e. φλ = λ = λ + λ + + n n 1 I A an 1 a n φ( A) = A + a A + + ai= n 1 n 1 From this we obtain: At e = α () ti+ α () ta+ + α ( t) A 1 n 1 n 1
12 System x = Ax + Bu y = Cx x R, u R, y R n m p t At A(t-s) x(t; x,u) = e x + e Bu(s)ds
13 Summary: Controllability/Observability Controllability rank = C CA Observability rank = n 1 CA Kalman Decomposition, x z such that d dt n 1 B AB A B n z1 A11 A12 A13 A14 z1 B1 z1 z 2 A22 A 24 z 2 B 2 z = 2 + u, y = [ C 2 C ] 4 z 3 A33 A 34 z 3 z 3 z4 A44 z4 z 4 z, z controllable z, z observable Notice that the substate z 2 is both controllable and observable
14 Example x =, [ 1 ] 4 x+ u y = x C = [ B AB] = not controllable 1 4 s s G s = C si A B = [ 1 ] [ 1 ] ( s 2 )( s 4 ) ( s + 2) 1 = = ( s+ 2)( s+ 4) s+ 4
15 Transfer Functions x = Ax + Bu y = Cx = [ ] = [ ] x R, u R, y R n m p 1 1 G s C si A B C si A B where A, B, C are those of the Kalman decomposition -, i.e., parameters of a minimal realization. So, only the controllable and observable part of the system is characterized by its transfer function. Definition: G s is called a complete characterization of the system if the system is completely observable and controllabl e.
16 Some MATLAB Functions Function canon ctrb ctrbf gram obsv obsvf ss2ss ssbal minreal canonical state-space realizations controllability matrix controllability staircase form controllability and observability gramians observability matrix observability staircase form state coordinate transformation diagonal balancing of state-space realizations returns a minimal realization
17 Example, Continued >> A=[-2-2;,-4]; >> B=[1;1]; >> C=[1 ]; >> sys=ss(a,b,c,); >> ctrb(sys) ans = >> gram(sys,'c') ans = >> tf(sys) Transfer function: s s^2 + 6 s + 8 >> tf(minreal(sys)) 1 state removed. Transfer function: s + 4
18 System Poles & Zeros Two descriptions of linear time-invariant systems state space and transfer function. x = Ax + Bu y = Cx + Du [ ] 1 G s = C si A B + D Assumption: Gis a complete characterization of ABCD,,, or, equivalently, ABCD,,, is a minimal realization of G. Defining poles via state space is very easy: the eigenvalues of A. the poles of Defining zeros is more complicated. We do it via state space in the followin g. G are
19 SISO System Zeros recall: x = Ax + Bu n x R, u Ry, R y = Cx + Du { } [ ] [ ] 1 1 Y () s = C si A x + C si A B + D U () s suppose: { [ ] } λt 1 ns () u t = e, G( s): = C si A B + D = k, d ( s) = si A ds () Y s [ ] Adj[ ] CAdj si A x ns () 1 C si A x ns () G( λ) = + k = + + ds () ds () s λ ds () ds () s λ (), Partial fraction expansion N( s) A( s) B = + ( λ) ( λ) ( s λ ) N( s) ( s λ ) A( s) lim λ s λ d( s) ( s λ ) d( s) N ( λ ) = d ( λ ) d s s d s s lim s B = + lim s λ ( s λ ) ( s λ ) B
20 SISO System Zeros, Cont d x can always be chosen so that [ ] CAdj si A x ns () + ds () ds () in which case ( λ ) G Y( s) =, if λ is a zero of G( λ ), then Y( s) =, y( t) = s λ In summary: if λ is a system zero, there exists x such that x t = x and u( t) = e λ t y t
21 MIMO System Zeros x = Ax + Bu y = Cx + Du x R, u R, y R n m p m n Does there exist g R and x R such that λt λt u() t = ge x() t = x e and y() t? The assumed solution must satisfy λt λt λt λxe = Axe + Bge λi A B x = C D g λt λt = Cxe + Dge The solution is nontrivial only if the system matrix is singular
22 MIMO System Zeros, Cont d This represents n+ p equations in n+ m unknowns. Suppose λi A B r = rank C D r = n+ min( mp, ) max r< n+ min( mp, ) nontrivial sol'ns p< m always nontrivial sol'ns r= r = n+ p m p max max independent sol'ns p> m and r= r = n+ m there are no nontrivial sol'ns
23 Square MIMO Systems (p=m) Nondegenerate case: If for typical λ, λi A B r = rank n m C D = + Those specific values of λ for which r < n+ m invariant zeros. Invariant zeros consist of are called input decoupling zeros (uncontrollable modes), λ satisfies rank [ λ ] I A B < n output decoupling zeros (unobservable modes), λ satisfies λi A rank n C < transmission zeros, all other invariant zeros.
24 Square MIMO Systems, Cont d degenerate case: ( λ ) For typical λ, λi A B r = rank < n+ m C D G insufficient independent controls,rank B< m insufficient independent outputs, rank C < p
25 Example - GTM θ = q x = V cosγ z = Vsinγ V = ( T cos α ρv SC 2 D( αδ, e, q) mg sin γ) m γ= ( T sin α+ ρv SC 2 L( αδ, e, q) mg cosγ) mv M q =, M = V ScC,, q + V ScC,, q x x mgx + lt I y α = θ γ 2 2 ( ) 1 ρ 1 2 m αδe ρ 2 Z αδe cgref cg cg t
26 Choose Equilibrium & Linearize Equilibrium Requires: = θ = q = x = V cosγ = z = Vsinγ = V = ( T cos α ρv SC 2 D( αδ, e, q) mg sin γ), α= θ γ m = γ= ( T sin α+ ρv SC 2 L( αδ, e, q) mg cosγ) mv M = q =, M = V ScC,, q + V ScC,, q x x mgx + lt I y 2 2 ( ) 1 ρ 1 2 m αδe ρ 2 Z αδe cgref cg cg t Also require: horizontal flight: γ = specify speed: V = V Solve 6 equations for remaining 4 states xz,, θ, q and 2 controls δ, T e
27 GTM Equilibrium Surface Straight and Level Flight Analysis 27
28 Controllability Around Bifurcation Point GTM Example: Straight Wings-Level Climb Note: Even though the system loses linear controllability at the bifurcation point, it is locally (nonlinearly) controllable in the sense that the controllability distribution has full generic rank around the bifurcation 28 point. Thus, any stabilizing feedback controller would have to be nonlinear and probably non-smooth (and possibly non-intuitive)
29 GTM Cont d Coordinated Turn Analysis Coordinated turn of 85 fps Coordinated turn of 9 fps Coordinated turn of 87 fps 29
30 GTM Cont d Coordinated Turn Analysis Techno-Sciences, Inc. Proprietary 3
31 Example: GTM Phugoid Dynamics T V δ e γ d V V T dt γ = γ δ e Normal trim d V V T dt γ = γ δ e High AOA trim Stall trim d V V T dt γ = γ δ e
32 Multivariable Interconnections & Feedback Loops 11/2/214
33 Well-Posed Loops: Example 1 G( s) - 1 s G( s) = 1 s 2 s+ 1 s+ 1 1 s+ 1 s 1 Gcl ( s) = G( s) I + G( s) = improper s 1 Theorem: Let GH, be proper rational transfer matrices. Then cl [ ] G = G I + HG 1 is proper and rational iff I + H G is nonsingular.
34 Poles of Closed Loops - G( s) H( s) cl [ ] G = G I + HG 1 It might be anticipated that the poles of the roots of det [ I + HG]. Not True!! G cl are s s s 1 s+ 1 G( s) =, H( s) = I2 2 1 s + 1 2s s s+ 1 s+ 1 Gcl ( s) =, det [ I + HG] 1 s 2 1 s 1 obviously, has poles at s = ± 1.
35 Poles, Cont d Theore m: If GH, are proper, reational matrices and det I + H G, then the poles of Gcl are the roots of the polynomial G s ( s) ( s) ( s) det I H( s) G( s) = G H + Example: = s ( s.5)( s 1.5) ( s. 5) ( s.5)( s+ 1.5 ) ( s.5), H = I 2
36 Poles: Example Cont d G cl ( s) ( s 1.5)( s.5 ), ( s) 1,det I G( s) G = + H = + = ( s) ( s) ( s) ( s) det I G( s) ( s 1.5)( s.5) = G H + = ( s.5) = 1 1 ( s 1.5) ( s ) s = ( s+ 1.5)( s+.5) s+.5 s s +.5 s.5
37 System Interconnections G 1 G 2 parallel G 1 series - G 1 G 2 feedback G 2 Systems G, G are complete characterizations, with Gi = Dli Nli = NriDri, i = 1, 2 coprime fractions. parallel connection controllable D, D left coprime r1 r2 observable Dl1, Dl2 right coprime series connection controllable D, N, or DD, N or D, N N are left coprime r2 r1 l1 r2 l1 l2 l2 r1 observable D, N, or DD, N or D, N N right coprime det I + G2 G1. Then controllable GG 1 2 controllable observable G 2 G 1 observable l1 r2 l1 r2 r2 r1 l2 r1
38 Stability x = Ax + Bu y = Cx + Du [ ] 1 x R, u R, y R G s = C si A B + D n m p The basic idea is that stable system responds to a perturbation by remaining within a neighborhood of its equilibrium point. a state perturbation with zero input (Lyapunov/Asymptotic) an input perturbation with zero state (BIBO) simultaneous state and input perturbation (Total) Lyapunov: Reλ A, eigenvalues with Reλ = have full set of eigenvectors. Asymptotic: Reλ A < BIBO: poles of G s < Total: Lyapunov + poles of G s <
39 Poles & Zeros from Transfer Functions 11/2/214
40 Numbers: Prime & Coprime A prime number (or integer) is a positive integer p > 1 that has no positive integer divisors other than 1 and itself. Two integers are relatively prime or coprime if they share no positive integer factors (divisors) other than 1 - i.e., their greatest common divisor is 1. Bezout's identity: exist integers x and y such that (, ) GCD a b = ax + by If a and b are integers not both zero, then there If a and b are coprime then there exist integers x and y such that 1 = ax + by
41 Polynomials These ideas have been extended to polynomials, matrices with polynomial elements matrices with rational elements Two polynomials ns = as + a s + + a, a m m 1 m m 1 m d s = bs + b s + + b, b n 1 n n 1 n are coprime if their greatest common divisor is a nonzero constant, i.e., they have no common factors.
42 Example F 1 = s 4 + 2s 3 + s + 2 = (s + 1) (s + 2) (s 2 - s + 1) F 2 = s 5 + s 4 + 2s 3 + 3s 2 + 3s + 2 = (s + 1) (s 2 - s + 2) (s 2 + s + 1) gcd(f 1,F 2 ) = (s + 1) Bezout relation: (5/24s 3 + 1/12s 2 + 1/4s + 5/24) F 1 + (-5/24s 2-7/24s + 7/24) F 2 = s + 1
43 Polynomial Matrices Two matrix polynomials m m 1 = n 1 i = n + n right common divisor R( s) R =, = D( s) right coprime R( s) N s As A s A have a m m 1 pa, q pb, q A R, Bi R D s Bs B s B N s N s R s D s D s R s N s, are if the only right common qq, divisors are unimodular, i.e., det = c. Similary, left coprimeness can be defined for polynomial matrices with the same number of rows. if
44 Poles Suppose the q m transfer matrix G s is a complete characterization of x = Ax + Bu, y = Cx + Du can always be factored into 1 1 G( s) = D ( s) N ( s) = N ( s) D ( s) G s l l r r where D, N and D, N are coprime pairs of polynomial matrices. N l l r r and D, D are called numerator, denominator matrices, respectively. Theorem: [ si A] = D s = D ( s) det α det α det α, α are constants. Definition: Poles are the roots of: det r r 1 l 2 r 1 2 [ si A], or det D ( s), or det D ( s) = = = l r r, N l
45 Poles & Zeros from Transfer Functions Assume G(s) is a complete characterization. Theorem: The pole polynomial (s) is the least common denominator of all non-identically-zero minors of all orders of G(s). Theorem: The zero polynomial is the greatest common divisor of all numerators of all order-r minors of G(s), where r is the generic rank of G(s), provided that these minors have been adjusted to have (s) as there denominator.
46 Example Recall a minor of a matrix is the determinant of a matrix obtained by deleting rows and columns. Consider the transfer function: 1 s 1 4 G( s) = ( s 1 s + ) To determine poles we need all minors of all orders. The 4 minors of order 1 are s s 1,,,2 s+ 2 s+ 2 s+ 2 s+ 2 rank G s =2. The single minor of order 2 is det G s = 2 ( s) pole polynomial: φ = s+ 2 zero polynomial: z s = s 4 s 4 ( s + 2)
Robust Control 2 Controllability, Observability & Transfer Functions
Robust Control 2 Controllability, Observability & Transfer Functions Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University /26/24 Outline Reachable Controllability Distinguishable
More informationJoão P. Hespanha. January 16, 2009
LINEAR SYSTEMS THEORY João P. Hespanha January 16, 2009 Disclaimer: This is a draft and probably contains a few typos. Comments and information about typos are welcome. Please contact the author at hespanha@ece.ucsb.edu.
More informationMULTIVARIABLE ZEROS OF STATE-SPACE SYSTEMS
Copyright F.L. Lewis All rights reserved Updated: Monday, September 9, 8 MULIVARIABLE ZEROS OF SAE-SPACE SYSEMS If a system has more than one input or output, it is called multi-input/multi-output (MIMO)
More informationControllability, Observability, Full State Feedback, Observer Based Control
Multivariable Control Lecture 4 Controllability, Observability, Full State Feedback, Observer Based Control John T. Wen September 13, 24 Ref: 3.2-3.4 of Text Controllability ẋ = Ax + Bu; x() = x. At time
More informationLetting be a field, e.g., of the real numbers, the complex numbers, the rational numbers, the rational functions W(s) of a complex variable s, etc.
1 Polynomial Matrices 1.1 Polynomials Letting be a field, e.g., of the real numbers, the complex numbers, the rational numbers, the rational functions W(s) of a complex variable s, etc., n ws ( ) as a
More informationContents. 1 State-Space Linear Systems 5. 2 Linearization Causality, Time Invariance, and Linearity 31
Contents Preamble xiii Linear Systems I Basic Concepts 1 I System Representation 3 1 State-Space Linear Systems 5 1.1 State-Space Linear Systems 5 1.2 Block Diagrams 7 1.3 Exercises 11 2 Linearization
More informationControl Systems Design, SC4026. SC4026 Fall 2010, dr. A. Abate, DCSC, TU Delft
Control Systems Design, SC4026 SC4026 Fall 2010, dr. A. Abate, DCSC, TU Delft Lecture 4 Controllability (a.k.a. Reachability) and Observability Algebraic Tests (Kalman rank condition & Hautus test) A few
More informationEEE582 Homework Problems
EEE582 Homework Problems HW. Write a state-space realization of the linearized model for the cruise control system around speeds v = 4 (Section.3, http://tsakalis.faculty.asu.edu/notes/models.pdf). Use
More informationSMITH MCMILLAN FORMS
Appendix B SMITH MCMILLAN FORMS B. Introduction Smith McMillan forms correspond to the underlying structures of natural MIMO transfer-function matrices. The key ideas are summarized below. B.2 Polynomial
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 informationMultivariable Control. Lecture 05. Multivariable Poles and Zeros. John T. Wen. September 14, 2006
Multivariable Control Lecture 05 Multivariable Poles and Zeros John T. Wen September 4, 2006 SISO poles/zeros SISO transfer function: G(s) = n(s) d(s) (no common factors between n(s) and d(s)). Poles:
More informationAdvanced Control Theory
State Space Solution and Realization chibum@seoultech.ac.kr Outline State space solution 2 Solution of state-space equations x t = Ax t + Bu t First, recall results for scalar equation: x t = a x t + b
More informationFull State Feedback for State Space Approach
Full State Feedback for State Space Approach State Space Equations Using Cramer s rule it can be shown that the characteristic equation of the system is : det[ si A] 0 Roots (for s) of the resulting polynomial
More informationControl Systems Design, SC4026. SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft
Control Systems Design, SC4026 SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft Lecture 4 Controllability (a.k.a. Reachability) vs Observability Algebraic Tests (Kalman rank condition & Hautus test) A few
More informationModule 03 Linear Systems Theory: Necessary Background
Module 03 Linear Systems Theory: Necessary Background Ahmad F. Taha EE 5243: Introduction to Cyber-Physical Systems Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ taha/index.html September
More informationChap 3. Linear Algebra
Chap 3. Linear Algebra Outlines 1. Introduction 2. Basis, Representation, and Orthonormalization 3. Linear Algebraic Equations 4. Similarity Transformation 5. Diagonal Form and Jordan Form 6. Functions
More informationLinear System Fundamentals
Linear Sytem Fundamental MEM 355 Performance Enhancement of Dynamical Sytem Harry G. Kwatny Department of Mechanical Engineering & Mechanic Drexel Univerity Content Sytem Repreentation Stability Concept
More informationIntro. Computer Control Systems: F8
Intro. Computer Control Systems: F8 Properties of state-space descriptions and feedback Dave Zachariah Dept. Information Technology, Div. Systems and Control 1 / 22 dave.zachariah@it.uu.se F7: Quiz! 2
More informationControl Systems Design
ELEC4410 Control Systems Design Lecture 14: Controllability Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 14: Controllability p.1/23 Outline
More informationMultivariable Control. Lecture 03. Description of Linear Time Invariant Systems. John T. Wen. September 7, 2006
Multivariable Control Lecture 3 Description of Linear Time Invariant Systems John T. Wen September 7, 26 Outline Mathematical description of LTI Systems Ref: 3.1-3.4 of text September 7, 26Copyrighted
More informationCONTROL DESIGN FOR SET POINT TRACKING
Chapter 5 CONTROL DESIGN FOR SET POINT TRACKING In this chapter, we extend the pole placement, observer-based output feedback design to solve tracking problems. By tracking we mean that the output is commanded
More informationRepresent this system in terms of a block diagram consisting only of. g From Newton s law: 2 : θ sin θ 9 θ ` T
Exercise (Block diagram decomposition). Consider a system P that maps each input to the solutions of 9 4 ` 3 9 Represent this system in terms of a block diagram consisting only of integrator systems, represented
More information16.31 Fall 2005 Lecture Presentation Mon 31-Oct-05 ver 1.1
16.31 Fall 2005 Lecture Presentation Mon 31-Oct-05 ver 1.1 Charles P. Coleman October 31, 2005 1 / 40 : Controllability Tests Observability Tests LEARNING OUTCOMES: Perform controllability tests Perform
More informationRaktim Bhattacharya. . AERO 422: Active Controls for Aerospace Vehicles. Dynamic Response
.. AERO 422: Active Controls for Aerospace Vehicles Dynamic Response Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. . Previous Class...........
More informationẋ 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)
EEE582 Topical Outline A.A. Rodriguez Fall 2007 GWC 352, 965-3712 The following represents a detailed topical outline of the course. It attempts to highlight most of the key concepts to be covered and
More informationPOLE PLACEMENT. Sadegh Bolouki. Lecture slides for ECE 515. University of Illinois, Urbana-Champaign. Fall S. Bolouki (UIUC) 1 / 19
POLE PLACEMENT Sadegh Bolouki Lecture slides for ECE 515 University of Illinois, Urbana-Champaign Fall 2016 S. Bolouki (UIUC) 1 / 19 Outline 1 State Feedback 2 Observer 3 Observer Feedback 4 Reduced Order
More informationMEM 355 Performance Enhancement of Dynamical Systems
MEM 355 Performance Enhancement of Dynamical Systems State Space Design Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 11/8/2016 Outline State space techniques emerged
More informationEE221A Linear System Theory Final Exam
EE221A Linear System Theory Final Exam Professor C. Tomlin Department of Electrical Engineering and Computer Sciences, UC Berkeley Fall 2016 12/16/16, 8-11am Your answers must be supported by analysis,
More informationUniversity of Toronto Department of Electrical and Computer Engineering ECE410F Control Systems Problem Set #3 Solutions = Q o = CA.
University of Toronto Department of Electrical and Computer Engineering ECE41F Control Systems Problem Set #3 Solutions 1. The observability matrix is Q o C CA 5 6 3 34. Since det(q o ), the matrix is
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 informationModule 09 From s-domain to time-domain From ODEs, TFs to State-Space Modern Control
Module 09 From s-domain to time-domain From ODEs, TFs to State-Space Modern Control Ahmad F. Taha EE 3413: Analysis and Desgin of Control Systems Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/
More informationModule 07 Controllability and Controller Design of Dynamical LTI Systems
Module 07 Controllability and Controller Design of Dynamical LTI Systems Ahmad F. Taha EE 5143: Linear Systems and Control Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ataha October
More informationDynamical Systems & Lyapunov Stability
Dynamical Systems & Lyapunov Stability Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University Outline Ordinary Differential Equations Existence & uniqueness Continuous dependence
More informationStability of Parameter Adaptation Algorithms. Big picture
ME5895, UConn, Fall 215 Prof. Xu Chen Big picture For ˆθ (k + 1) = ˆθ (k) + [correction term] we haven t talked about whether ˆθ(k) will converge to the true value θ if k. We haven t even talked about
More informationEL2520 Control Theory and Practice
So far EL2520 Control Theory and Practice r Fr wu u G w z n Lecture 5: Multivariable systems -Fy Mikael Johansson School of Electrical Engineering KTH, Stockholm, Sweden SISO control revisited: Signal
More informationState Space Design. MEM 355 Performance Enhancement of Dynamical Systems
State Space Design MEM 355 Performance Enhancement of Dynamical Systems Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University Outline State space techniques emerged around
More informationDifferential equations
Differential equations Math 7 Spring Practice problems for April Exam Problem Use the method of elimination to find the x-component of the general solution of x y = 6x 9x + y = x 6y 9y Soln: The system
More informationLinear Algebra III Lecture 11
Linear Algebra III Lecture 11 Xi Chen 1 1 University of Alberta February 13, 2015 Outline Minimal Polynomial 1 Minimal Polynomial Minimal Polynomial The minimal polynomial f (x) of a square matrix A is
More informationA Necessary and Sufficient Condition for High-Frequency Robustness of Non-Strictly-Proper Feedback Systems
A Necessary and Sufficient Condition for High-Frequency Robustness of Non-Strictly-Proper Feedback Systems Daniel Cobb Department of Electrical Engineering University of Wisconsin Madison WI 53706-1691
More informationECE 388 Automatic Control
Controllability and State Feedback Control Associate Prof. Dr. of Mechatronics Engineeering Çankaya University Compulsory Course in Electronic and Communication Engineering Credits (2/2/3) Course Webpage:
More informationCDS Solutions to the Midterm Exam
CDS 22 - Solutions to the Midterm Exam Instructor: Danielle C. Tarraf November 6, 27 Problem (a) Recall that the H norm of a transfer function is time-delay invariant. Hence: ( ) Ĝ(s) = s + a = sup /2
More informationFEL3210 Multivariable Feedback Control
FEL3210 Multivariable Feedback Control Lecture 8: Youla parametrization, LMIs, Model Reduction and Summary [Ch. 11-12] Elling W. Jacobsen, Automatic Control Lab, KTH Lecture 8: Youla, LMIs, Model Reduction
More informationControl Systems I. Lecture 4: Diagonalization, Modal Analysis, Intro to Feedback. Readings: Emilio Frazzoli
Control Systems I Lecture 4: Diagonalization, Modal Analysis, Intro to Feedback Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 13, 2017 E. Frazzoli (ETH)
More informationState will have dimension 5. One possible choice is given by y and its derivatives up to y (4)
A Exercise State will have dimension 5. One possible choice is given by y and its derivatives up to y (4 x T (t [ y(t y ( (t y (2 (t y (3 (t y (4 (t ] T With this choice we obtain A B C [ ] D 2 3 4 To
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 informationLec 6: State Feedback, Controllability, Integral Action
Lec 6: State Feedback, Controllability, Integral Action November 22, 2017 Lund University, Department of Automatic Control Controllability and Observability Example of Kalman decomposition 1 s 1 x 10 x
More informationSolution of Linear State-space Systems
Solution of Linear State-space Systems Homogeneous (u=0) LTV systems first Theorem (Peano-Baker series) The unique solution to x(t) = (t, )x 0 where The matrix function is given by is called the state
More informationLinear State Feedback Controller Design
Assignment For EE5101 - Linear Systems Sem I AY2010/2011 Linear State Feedback Controller Design Phang Swee King A0033585A Email: king@nus.edu.sg NGS/ECE Dept. Faculty of Engineering National University
More informationLecture 3. Chapter 4: Elements of Linear System Theory. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.
Lecture 3 Chapter 4: Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 3 p. 1/77 3.1 System Descriptions [4.1] Let f(u) be a liner operator, u 1 and u
More informationChapter 5. Linear Algebra. A linear (algebraic) equation in. unknowns, x 1, x 2,..., x n, is. an equation of the form
Chapter 5. Linear Algebra A linear (algebraic) equation in n unknowns, x 1, x 2,..., x n, is an equation of the form a 1 x 1 + a 2 x 2 + + a n x n = b where a 1, a 2,..., a n and b are real numbers. 1
More informationStability, Pole Placement, Observers and Stabilization
Stability, Pole Placement, Observers and Stabilization 1 1, The Netherlands DISC Course Mathematical Models of Systems Outline 1 Stability of autonomous systems 2 The pole placement problem 3 Stabilization
More informationLecture 7 (Weeks 13-14)
Lecture 7 (Weeks 13-14) Introduction to Multivariable Control (SP - Chapters 3 & 4) Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 7 (Weeks 13-14) p.
More informationMA 527 first midterm review problems Hopefully final version as of October 2nd
MA 57 first midterm review problems Hopefully final version as of October nd The first midterm will be on Wednesday, October 4th, from 8 to 9 pm, in MTHW 0. It will cover all the material from the classes
More informationModule 08 Observability and State Estimator Design of Dynamical LTI Systems
Module 08 Observability and State Estimator Design of Dynamical LTI Systems Ahmad F. Taha EE 5143: Linear Systems and Control Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ataha November
More informationEE Control Systems LECTURE 9
Updated: Sunday, February, 999 EE - Control Systems LECTURE 9 Copyright FL Lewis 998 All rights reserved STABILITY OF LINEAR SYSTEMS We discuss the stability of input/output systems and of state-space
More informationECEN 605 LINEAR SYSTEMS. Lecture 7 Solution of State Equations 1/77
1/77 ECEN 605 LINEAR SYSTEMS Lecture 7 Solution of State Equations Solution of State Space Equations Recall from the previous Lecture note, for a system: ẋ(t) = A x(t) + B u(t) y(t) = C x(t) + D u(t),
More informationCayley-Hamilton Theorem
Cayley-Hamilton Theorem Massoud Malek In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n Let A be an n n matrix Although det (λ I n A
More information1. Find the solution of the following uncontrolled linear system. 2 α 1 1
Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +
More informationLinear System Theory
Linear System Theory Wonhee Kim Chapter 6: Controllability & Observability Chapter 7: Minimal Realizations May 2, 217 1 / 31 Recap State space equation Linear Algebra Solutions of LTI and LTV system Stability
More information7 Minimal realization and coprime fraction
7 Minimal realization and coprime fraction 7 Introduction If a tranfer function i realizable, what i the mallet poible dimenion? Realization with the mallet poible dimenion are called minimal-dimenional
More informationSolution via Laplace transform and matrix exponential
EE263 Autumn 2015 S. Boyd and S. Lall Solution via Laplace transform and matrix exponential Laplace transform solving ẋ = Ax via Laplace transform state transition matrix matrix exponential qualitative
More informationTheorem 1. ẋ = Ax is globally exponentially stable (GES) iff A is Hurwitz (i.e., max(re(σ(a))) < 0).
Linear Systems Notes Lecture Proposition. A M n (R) is positive definite iff all nested minors are greater than or equal to zero. n Proof. ( ): Positive definite iff λ i >. Let det(a) = λj and H = {x D
More informationA matrix is a rectangular array of. objects arranged in rows and columns. The objects are called the entries. is called the size of the matrix, and
Section 5.5. Matrices and Vectors A matrix is a rectangular array of objects arranged in rows and columns. The objects are called the entries. A matrix with m rows and n columns is called an m n matrix.
More information0.1 Rational Canonical Forms
We have already seen that it is useful and simpler to study linear systems using matrices. But matrices are themselves cumbersome, as they are stuffed with many entries, and it turns out that it s best
More informationLecture plan: Control Systems II, IDSC, 2017
Control Systems II MAVT, IDSC, Lecture 8 28/04/2017 G. Ducard Lecture plan: Control Systems II, IDSC, 2017 SISO Control Design 24.02 Lecture 1 Recalls, Introductory case study 03.03 Lecture 2 Cascaded
More informationA matrix is a rectangular array of. objects arranged in rows and columns. The objects are called the entries. is called the size of the matrix, and
Section 5.5. Matrices and Vectors A matrix is a rectangular array of objects arranged in rows and columns. The objects are called the entries. A matrix with m rows and n columns is called an m n matrix.
More information6 Linear Equation. 6.1 Equation with constant coefficients
6 Linear Equation 6.1 Equation with constant coefficients Consider the equation ẋ = Ax, x R n. This equating has n independent solutions. If the eigenvalues are distinct then the solutions are c k e λ
More informationMath Exam 3 Solutions
Math 6 - Exam 3 Solutions Thursday, July 3rd, 0 Recast the following higher-order differential equations into first order systems If the equation is linear, be sure to give the coefficient matrix At and
More informationRobust Control 3 The Closed Loop
Robust Control 3 The Closed Loop Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University /2/2002 Outline Closed Loop Transfer Functions Traditional Performance Measures Time
More informationEquality: Two matrices A and B are equal, i.e., A = B if A and B have the same order and the entries of A and B are the same.
Introduction Matrix Operations Matrix: An m n matrix A is an m-by-n array of scalars from a field (for example real numbers) of the form a a a n a a a n A a m a m a mn The order (or size) of A is m n (read
More informationModeling and Analysis of Dynamic Systems
Modeling and Analysis of Dynamic Systems Dr. Guillaume Ducard Fall 2017 Institute for Dynamic Systems and Control ETH Zurich, Switzerland G. Ducard c 1 / 57 Outline 1 Lecture 13: Linear System - Stability
More informationBalanced Truncation 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 2004: MODEL REDUCTION Balanced Truncation This lecture introduces balanced truncation for LTI
More informationThe Important State Coordinates of a Nonlinear System
The Important State Coordinates of a Nonlinear System Arthur J. Krener 1 University of California, Davis, CA and Naval Postgraduate School, Monterey, CA ajkrener@ucdavis.edu Summary. We offer an alternative
More informationMATH 3321 Sample Questions for Exam 3. 3y y, C = Perform the indicated operations, if possible: (a) AC (b) AB (c) B + AC (d) CBA
MATH 33 Sample Questions for Exam 3. Find x and y so that x 4 3 5x 3y + y = 5 5. x = 3/7, y = 49/7. Let A = 3 4, B = 3 5, C = 3 Perform the indicated operations, if possible: a AC b AB c B + AC d CBA AB
More informationECE557 Systems Control
ECE557 Systems Control Bruce Francis Course notes, Version.0, September 008 Preface This is the second Engineering Science course on control. It assumes ECE56 as a prerequisite. If you didn t take ECE56,
More informationEL 625 Lecture 10. Pole Placement and Observer Design. ẋ = Ax (1)
EL 625 Lecture 0 EL 625 Lecture 0 Pole Placement and Observer Design Pole Placement Consider the system ẋ Ax () The solution to this system is x(t) e At x(0) (2) If the eigenvalues of A all lie in the
More informationZeros and zero dynamics
CHAPTER 4 Zeros and zero dynamics 41 Zero dynamics for SISO systems Consider a linear system defined by a strictly proper scalar transfer function that does not have any common zero and pole: g(s) =α p(s)
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 informationObjective Type Questions
DISTANCE EDUCATION, UNIVERSITY OF CALICUT NUMBER THEORY AND LINEARALGEBRA Objective Type Questions Shyama M.P. Assistant Professor Department of Mathematics Malabar Christian College, Calicut 7/3/2014
More informationCDS Solutions to Final Exam
CDS 22 - Solutions to Final Exam Instructor: Danielle C Tarraf Fall 27 Problem (a) We will compute the H 2 norm of G using state-space methods (see Section 26 in DFT) We begin by finding a minimal state-space
More informationMath 489AB Exercises for Chapter 2 Fall Section 2.3
Math 489AB Exercises for Chapter 2 Fall 2008 Section 2.3 2.3.3. Let A M n (R). Then the eigenvalues of A are the roots of the characteristic polynomial p A (t). Since A is real, p A (t) is a polynomial
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science : Dynamic Systems Spring 2011
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science 6.4: Dynamic Systems Spring Homework Solutions Exercise 3. a) We are given the single input LTI system: [
More informationACM/CMS 107 Linear Analysis & Applications Fall 2016 Assignment 4: Linear ODEs and Control Theory Due: 5th December 2016
ACM/CMS 17 Linear Analysis & Applications Fall 216 Assignment 4: Linear ODEs and Control Theory Due: 5th December 216 Introduction Systems of ordinary differential equations (ODEs) can be used to describe
More informationNonlinear Control Systems
Nonlinear Control Systems António Pedro Aguiar pedro@isr.ist.utl.pt 7. Feedback Linearization IST-DEEC PhD Course http://users.isr.ist.utl.pt/%7epedro/ncs1/ 1 1 Feedback Linearization Given a nonlinear
More informationPerspective. ECE 3640 Lecture 11 State-Space Analysis. To learn about state-space analysis for continuous and discrete-time. Objective: systems
ECE 3640 Lecture State-Space Analysis Objective: systems To learn about state-space analysis for continuous and discrete-time Perspective Transfer functions provide only an input/output perspective of
More informationLinear Algebra: Matrix Eigenvalue Problems
CHAPTER8 Linear Algebra: Matrix Eigenvalue Problems Chapter 8 p1 A matrix eigenvalue problem considers the vector equation (1) Ax = λx. 8.0 Linear Algebra: Matrix Eigenvalue Problems Here A is a given
More informationChap. 3. Controlled Systems, Controllability
Chap. 3. Controlled Systems, Controllability 1. Controllability of Linear Systems 1.1. Kalman s Criterion Consider the linear system ẋ = Ax + Bu where x R n : state vector and u R m : input vector. A :
More informationc c c c c c c c c c a 3x3 matrix C= has a determinant determined by
Linear Algebra Determinants and Eigenvalues Introduction: Many important geometric and algebraic properties of square matrices are associated with a single real number revealed by what s known as the determinant.
More information3.1 Overview 3.2 Process and control-loop interactions
3. Multivariable 3.1 Overview 3.2 and control-loop interactions 3.2.1 Interaction analysis 3.2.2 Closed-loop stability 3.3 Decoupling control 3.3.1 Basic design principle 3.3.2 Complete decoupling 3.3.3
More informationTopic # /31 Feedback Control Systems
Topic #7 16.30/31 Feedback Control Systems State-Space Systems What are the basic properties of a state-space model, and how do we analyze these? Time Domain Interpretations System Modes Fall 2010 16.30/31
More informationOPTIMAL CONTROL SYSTEMS
SYSTEMS MIN-MAX Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University OUTLINE MIN-MAX CONTROL Problem Definition HJB Equation Example GAME THEORY Differential Games Isaacs
More informationFall 線性系統 Linear Systems. Chapter 08 State Feedback & State Estimators (SISO) Feng-Li Lian. NTU-EE Sep07 Jan08
Fall 2007 線性系統 Linear Systems Chapter 08 State Feedback & State Estimators (SISO) Feng-Li Lian NTU-EE Sep07 Jan08 Materials used in these lecture notes are adopted from Linear System Theory & Design, 3rd.
More informationControl Systems. System response. L. Lanari
Control Systems m i l e r p r a in r e v y n is o System response L. Lanari Outline What we are going to see: how to compute in the s-domain the forced response (zero-state response) using the transfer
More informationECEEN 5448 Fall 2011 Homework #4 Solutions
ECEEN 5448 Fall 2 Homework #4 Solutions Professor David G. Meyer Novemeber 29, 2. The state-space realization is A = [ [ ; b = ; c = [ which describes, of course, a free mass (in normalized units) with
More information1. What is the determinant of the following matrix? a 1 a 2 4a 3 2a 2 b 1 b 2 4b 3 2b c 1. = 4, then det
What is the determinant of the following matrix? 3 4 3 4 3 4 4 3 A 0 B 8 C 55 D 0 E 60 If det a a a 3 b b b 3 c c c 3 = 4, then det a a 4a 3 a b b 4b 3 b c c c 3 c = A 8 B 6 C 4 D E 3 Let A be an n n matrix
More informationChapter 30 Minimality and Stability of Interconnected Systems 30.1 Introduction: Relating I/O and State-Space Properties We have already seen in Chapt
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A. Dahleh George Verghese Department of Electrical Engineering and Computer Science Massachuasetts Institute of Technology 1 1 c Chapter
More informationDecay rates for partially dissipative hyperbolic systems
Outline Decay rates for partially dissipative hyperbolic systems Basque Center for Applied Mathematics Bilbao, Basque Country, Spain zuazua@bcamath.org http://www.bcamath.org/zuazua/ Numerical Methods
More informationLecture 4 Stabilization
Lecture 4 Stabilization This lecture follows Chapter 5 of Doyle-Francis-Tannenbaum, with proofs and Section 5.3 omitted 17013 IOC-UPC, Lecture 4, November 2nd 2005 p. 1/23 Stable plants (I) We assume that
More informationRobust Control 5 Nominal Controller Design Continued
Robust Control 5 Nominal Controller Design Continued Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 4/14/2003 Outline he LQR Problem A Generalization to LQR Min-Max
More informationProblem Set 5 Solutions 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Problem Set 5 Solutions The problem set deals with Hankel
More information