Discrete-time linear systems
|
|
- Suzanna Gregory
- 6 years ago
- Views:
Transcription
1 Automatic Control Discrete-time linear systems Prof. Alberto Bemporad University of Trento Academic year 2-2 Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year 2-2 / 34
2 Introduction Introduction 3.5 y(t), y(kt s ) 4 u(kt s ) time t Sampling of a continuous signal time t Discrete-time signal Discrete-time models describe relationships between sampled variables x(kt s ), u(kt s ), y(kt s ), k =,,... The value x(kt s ) is kept constant during the sampling interval [kt s, (k + )T s ) A discrete-time signal can either represent the sampling of a continuous-time signal, or be an intrinsically discrete signal Discrete-time signals are at the basis of digital controllers (as well as of digital filters in signal processing) Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
3 Difference equations Difference equation Consider the first order difference equation (autonomous system) x(k + ) = ax(k) x() = x The solution is x(k) = a k x a> 2 x(k).5 a=.5 <a< k Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
4 Difference equations Difference equation First-order difference equation with input (non-autonomous system) The solution has the form x(k + ) = ax(k) + bu(k) x() = x x(k) = a k x }{{} natural response a i bu(k i) k + i= } {{ } forced response The natural response depends on the initial condition x(), the forced response on the input signal u(k) Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
5 Difference equations Example - Wealth of a bank account k: year counter ρ: interest rate x(k): wealth at the beginning of year k u(k): money saved at the end of year k x : initial wealth in bank account Discrete-time model: x k u(k) 5 k ρ % x(k + ) = ( + ρ)x(k) + u(k) x() = x Stored amount of money (keur) x(k) 25 2 x(k) = (.) k + (.)k. 5 = 6(.)k k (years) Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
6 Difference equations Linear discrete-time system Consider the set of n first-order linear difference equations forced by the input u(k) x (k + ) = a x (k) a n x n (k) +b u(k) x 2 (k + ) = a 2 x (k) a 2n x n (k) +b 2 u(k)... x n (k + ) = a n x (k) a nn x n (k) +b n u(k) x () = x,... x n () = x n In compact matrix form: x(k + ) = Ax(k) + Bu(k) where x = x.. n. x n x() = x Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
7 Difference equations Linear discrete-time system The solution is x(k) = A k x }{{} natural response If matrix A is diagonalizable, A = TΛT Λ =. λ... λ A i Bu(k i) k + i= A k = T.... λ n } {{ } forced response λ k... λ k λ k n where T = [v... v n ] collects n independent eigenvectors. T Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
8 Modal response Modal response Assume input u(k) =, k Assume A is diagonalizable, A = TΛT The state trajectory (natural response) is where x(k) = A k x = TΛ k T x = n α i λ k i v i λ i = eigenvalues of A v i = eigenvectors of A α i = coefficients that depend on the initial condition x() α = i= α.. = T x(), T = [v... v n ] α n The system modes depend on the eigenvalues of A, as in continuous-time Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
9 Modal response Example Consider the linear discrete-time system x (k + ) = 2 x (k) + 2 x 2(k) x 2 (k + ) = x 2 (k) + u(k) x () = x 2 () = x(k + ) = x() = 2 2 x(k) + u(k) Eigenvalues of A: λ = 2, λ 2 = Solution: x(k) = = = 2 2 k k + 2 i= 2 k 2 k k + i= k k 2 + i= }{{} natural response 2 i u(k i) 2 i u(k i) 2 i u(k i) } {{ } forced response x (k),x 2 (k) step k simulation for u(k) Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
10 Linear difference equations n th -order difference equation Consider the n th -order difference equation forced by u a n y(k n) + a n y(k n + ) + + a y(k ) + y(k) = b n u(k n) + + b u(k ) + b u(k) Equivalent linear discrete-time system in canonical state matrix form x(k + ) =... x(k) +. u(k)... a n a n a n 2... a y(k) = (b n b a n )... (b b a ) x(k) + b u(k) There are infinitely many state-space realizations MATLAB tf2ss n th -order difference equations are very useful for digital filters, digital controllers, and to reconstruct models from data (system identification) Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year 2-2 / 34
11 Discrete-time linear systems Discrete-time linear system x(k + ) = Ax(k) + Bu(k) y(k) = Cx(k) + Du(k) x() = x Given the initial condition x() and the input sequence u(k), k, it is possible to predict the entire sequence of states x(k) and outputs y(k), k The state x() summarizes all the past history of the system The dimension n of the state x(k) n is called the order of the system The system is called proper (or strictly causal) if D = General multivariable case: x(k) n u(k) m y(k) p A n n B n m C p n D p m Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year 2-2 / 34
12 Discrete-time linear systems Example - Student dynamics Problem Statement: 3-years undergraduate course percentages of students promoted, repeaters, and dropouts are roughly constant direct enrollment in 2nd and 3rd academic year is not allowed students cannot enrol for more than 3 years Notation: k Year x i (k) Number of students enrolled in year i at year k, i =, 2, 3 u(k) Number of freshmen at year k y(k) Number of graduates at year k α i promotion rate during year i, α i β i failure rate during year i, β i γ i dropout rate during year i, γ i = α i β i 3 rd -order linear discrete-time system: x (k + ) = β x (k) + u(k) x 2 (k + ) = α x (k) + β 2 x 2 (k) x 3 (k + ) = α 2 x 2 (k) + β 3 x 3 (k) y(k) = α 3 x 3 (k) Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
13 Discrete-time linear systems Example - Student dynamics In matrix form β x(k + ) = α β 2 x(k) + u(k) α 2 β 3 y(k) = α 3 x(k) Simulation 4 y(k) α =.6 β =.2 α 2 =.8 β 2 =.5 α 3 =.9 β 3 =.8 u(k) 5, k = 22, step k u(k) step k Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
14 Lecture: Discrete-time linear systems Discrete-time linear systems Example - Supply chain!x(k)!2x2(k) "2x2(k) "x(k) u(k) #3x3(k) P S x(k) y(k) $3x3(k) R x2(k) x3(k) Problem Statement: S purchases the quantity u(k) of raw material at each month k a fraction δ of raw material is discarded, a fraction α is shipped to producer P a fraction α2 of product is sold by P to retailer R, a fraction δ2 is discarded retailer R returns a fraction β3 of defective products every month, and sells a fraction γ3 to customers Mathematical model: x (k + ) x2 (k + ) x3 (k + ) y(k) = = = = ( α δ )x (k) + u(k) α x (k) + ( α2 δ2 )x2 (k) + β3 x3 (k) α2 x2 (k) + ( β3 γ3 )x3 (k) γ3 x3 (k) Prof. Alberto Bemporad (University of Trento) Automatic Control k x (k) x2 (k) x3 (k) y(k) month counter raw material in S products in P products in R products sold to customers Academic year / 34
15 Stability of discrete-time linear systems Equilibrium Consider the discrete-time nonlinear system x(k + ) = f(x(k), u(k)) y(k) = g(x(k), u(k)) Definition A state x r n and an input u r m are an equilibrium pair if for initial condition x() = x r and constant input u(k) u r, k, the state remains constant: x(k) x r, k Equivalent definition: (x r, u r ) is an equilibrium pair if f(x r, u r ) = x r x r is called equilibrium state, u r equilibrium input The definition generalizes to time-varying discrete-time nonlinear systems Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
16 Stability of discrete-time linear systems Stability Consider the nonlinear system x(k + ) = f(x(k), ur ) y(k) = g(x(k), u r ) and let x r an equilibrium state, f(x r, u r ) = x r Definition The equilibrium state x r is stable if for each initial conditions x() close enough to x r, the corresponding trajectory x(k) remains near x r for all k a a Analytic definition: ε > δ > : x() x r < δ x(k) x r < ε, k The equilibrium point x r is called asymptotically stable if it is stable and x(k) x r for k Otherwise, the equilibrium point x r is called unstable Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
17 Stability of discrete-time linear systems Stability of first-order linear systems Consider the first-order linear system x(k + ) = ax(k) + bu(k) x r =, u r = is an equilibrium pair For u(k), k =,,..., the solution is The origin x r = is x(k) = a k x unstable if a > stable if a asymptotically stable if a < x(k) x a> a= <a<.5 x r a= k Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
18 Stability of discrete-time linear systems Stability of discrete-time linear systems Since the natural response of x(k + ) = Ax(k) + Bu(k) is x(k) = A k x, the stability properties depend only on A. We can therefore talk about system stability of a discrete-time linear system (A, B, C, D) Theorem: Let λ,..., λ m, m n be the eigenvalues of A n n. The system x(k + ) = Ax(k) + Bu(k) is asymptotically stable iff λ i <, i =,..., m (marginally) stable if λ i, i =,..., m, and the eigenvalues with unit modulus have equal algebraic and geometric multiplicity a unstable if i such that λ i > a Algebraic multiplicity of λ i = number of coincident roots λ i of det(λi A). Geometric multiplicity of λ i = number of linearly independent eigenvectors v i, Av i = λ i v i The stability properties of a discrete-time linear system only depend on the modulus of the eigenvalues of matrix A Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
19 Stability of discrete-time linear systems Stability of discrete-time linear systems Proof: The natural response is x(k) = A k x If matrix A is diagonalizable, A = TΛT, λ... λ 2... Take any eigenvalue λ = ρe jθ : Λ = A k = T λ n λ k... λ k λ k = ρ k e jkθ = ρ k λ k n T A is always diagonalizable if algebraic multiplicity - geometric multiplicity If A is not diagonalizable, it can be transformed to Jordan form. In this case the natural response x(t) contains modes k j λ k, j =,,..., alg. multiplicity = geom. multiplicity Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
20 Stability of discrete-time linear systems Example x(k + ) = x(k) 2 x() = x x 2 eigenvalues of A:, 2 solution: x (k) =, k =, 2,... x 2 (k) = k 2 x + k 2 x2, k =, 2, x (k) x 2 (k).5 asymptotically stable k k Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
21 Stability of discrete-time linear systems Example 2 x(k + ) = x() = x x 2 x(k) eigenvalues of A: {+j, j} solution: x (k) = x cos kπ 2 + x 2 sin kπ 2, k =,,... x 2 (k) = x sin kπ 2 + x 2 cos kπ 2, k =,,..., x (k).2 x 2 (k).2 marginally stable k k Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
22 Stability of discrete-time linear systems Example 3 x(k + ) = x(k) x() = x x 2 x (k) = x + x 2 k, k =,,... x 2 (k) = x 2, k =,,... eigenvalues of A: {, } Note: A is not diagonalizable! x (k) 2 x 2 (k) unstable k k Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
23 Stability of discrete-time linear systems Example 4 2 x(k + ) = x(k) x() = x x 2 eigenvalues of A: {, 2} solution: x (k) = 2 k x, k =,,... x 2 (k) = 2 k x, k =, 2, x (k) 2 x 2 (k) 2 unstable k k Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
24 Stability of discrete-time linear systems Zero eigenvalues Modes λ i = determine finite-time convergence to zero. This has no continuous-time counterpart, where instead all converging modes tend to zero in infinite time (e λ it ) Example: dynamics of a buffer u(k) x 3(k) x 2(k) x (k) y(k) x (k + ) = x 2 (k) x 2 (k + ) = x 3 (k) x 3 (k + ) = u(k) y(k) = x (k) x(k + ) = x(k) + u(k) y(k) = x(k) Natural response: A 3 x() = for all x() 3 For u(k), the buffer deploys after at most 3 steps! Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
25 Stability of discrete-time linear systems Summary of stability conditions for linear systems system continuous-time discrete-time asympt. stable i =,..., n R(λ i ) < λ i < unstable i such that R(λ i ) > λ i > stable i,..., n R(λ i ) λ i and λ i such that R(λ i ) = λ i = algebraic = geometric mult. Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
26 Sampling continuous-time systems Exact sampling Consider the continuos-time system ẋ(t) = Ax(t) + Bu(t) y(t) = Cx(t) + Du(t) x() = x We want to characterize the value of x(t), y(t) at the time instants t =, T s, 2T s,..., kt s,..., under the assumption that the input u(t) is constant during each sampling interval (zero-order hold, ZOH).5 y(t), y(kt s ) 2 u(t), u(kt s ).5 u(t) = ū(k), kt s t < (k + )T s.5.5 x(k) x(kt s ) and ȳ(k) y(kt s ) are the state and the output samples at the k th sampling instant, respectively time t time t Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
27 Sampling continuous-time systems Exact sampling Let us evaluate the response of the continuous-time system between time t = kt s and t = (k + )T s from the initial condition x(t ) = x(kt s ) using Lagrange formula: x(t) = e A(t t ) x(t )+ t t e A(t σ) Bu(σ)dσ = e A((k+)Ts kts) x(kt s )+ (k+)ts kt s e A((k+)Ts σ) Bu(σ)dσ Since the input u(t) is piecewise constant, u(σ) ū(k), kt s σ < (k + )T s. By setting τ = σ kt s we get Ts x((k + )T s ) = e AT s x(kt s ) + e A(T s τ) dτ Bu(kT s ) and hence Ts x(k + ) = e ATs x(k) + e A(T s τ) dτ Bū(k) which is a linear difference relation between x(k) and ū(k)! Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
28 Sampling continuous-time systems Exact sampling The discrete-time system x(k + ) = Ā x(k) + Bū(k) ȳ(k) = C x(k) + Dū(k) depends on the original continuous-time system through the relations Ā e AT s, B Ts e A(T s τ) dτ B, C C, D D If u(t) is piecewise constant, (Ā, B, C, D) provides the exact evolution of state and output samples at discrete times kt s.4 y(t), y(kt s ).2 MATLAB sys=ss(a,b,c,d); sysd=c2d(sys,ts); [Ab,Bb,Cb,Db]=ssdata(sysd); time t.5 u(t), u(kt s ) time t Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
29 Sampling continuous-time systems Exact sampling Rule of thumb: T s of the rise time = time to move from % to 9% of the steady-state value, for input u(t), x() =.9 9 % rise time.8 y(t), y(kt s ) output time t u(t), u(kt s ) % rise time time t time t More on the choice of sampling time in the second part of the course... Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
30 Sampling continuous-time systems Approximate sampling - Euler s method ẋ(kt s ) x((k + )T s) x(kt s ) T s x(t) x((k x((k+)t + )Ts) ) x((k + x((k+)t )Ts) ) x(kt x(kts) ) T Ts x(kts) ) For nonlinear systems ẋ(t) = f(x(t), u(t)): kt kts x(k + ) = x(k) + T s f( x(k), ū(k)) Ts T ẋ(kts) _x(t) (k (k+)t + )Ts t Leonhard Paul Euler (77-783) For linear systems ẋ(t) = Ax(t) + Bu(t): x((k + )T s ) = (I + T s A)x(kT s ) + T s Bu(kT s ) Ā I + AT s, B T s B, C C, D D Note that e T sa = I + T s A Tn s An n! +... Therefore when T s is small Euler s method and exact sampling are similar Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
31 Sampling continuous-time systems Example - Hydraulic system Continuous time: d dt h(t) = a 2g A h(t) + A u(t) q(t) = a 2g h(t) Discrete time: T h(k + ) = h(k) s a 2g T A h(k) + s A ū(k) q(k) = a 2g h(k) u 7 6 level h(t) (m) 5 4 h q 3 2 continuous time Euler approximation time (s) Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
32 Sampling continuous-time systems Tustin s discretization method Assume u(k) constant within the sampling interval. Given the linear system ẋ = Ax + Bu, apply the trapezoidal rule to approximate the integral and therefore x(k + ) x(k) = (k+)ts kt s ẋ(t)dt = (k+)ts kt s (Ax(t) + Bu(t))dt T s (Ax(k) + Bu(k) + Ax(k + ) + Bu(k)) (trapezoidal rule) 2 (I T s 2 A)x(k + ) = (I + T s 2 )x(k) + T sbu(k) x(k + ) = I T s 2 A I + T s 2 A x(k) + I T s 2 A T s Bu(k) Advantage: simpler to compute than exponential matrix, without too much loss of approximation quality Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
33 Sampling continuous-time systems N-steps Euler method We can obtain the matrices A, B of the discrete-time linearized model while integrating the nonlinear continuous-time dynamic equations ẋ = f(x, u) N-steps explicit forward Euler method: given x(k), u(k), execute the following steps x = x(k), A = I, B = 2 for n=:n do 3 end A (I + T s N B (I + T s N x x + T s N f x (x, u(k))a f x (x, u(k))b + T s N f(x, u(k)) f u (x, u(k))a 4 return x(k + ) x and matrices A, B such that x(k + ) Ax(k) + Bu(k). Property: the difference between the state x(k + ) and its approximation x computed by the above iterations satisfies x(k + ) x) = O T s N Explicit forward Runge-Kutta 4 method also available Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
34 Sampling continuous-time systems English-Italian Vocabulary discrete-time linear systems sampling interval difference equation zero-order hold piecewise constant rise time sistemi lineari a tempo discreto tempo (o intervallo) di campionamento equazione alle differenze mantenitore di ordine zero costante a tratti tempo di salita Translation is obvious otherwise. Prof. Alberto Bemporad (University of Trento) Automatic Control Academic year / 34
Automatic Control 2. Model reduction. Prof. Alberto Bemporad. University of Trento. Academic year
Lecture: Automatic Control 2 Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011 1 / 17 Lecture:
More information4F3 - Predictive Control
4F3 Predictive Control - Discrete-time systems p. 1/30 4F3 - Predictive Control Discrete-time State Space Control Theory For reference only Jan Maciejowski jmm@eng.cam.ac.uk 4F3 Predictive Control - Discrete-time
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 informationIntegral action in state feedback control
Automatic Control 1 in state feedback control Prof. Alberto Bemporad University of Trento Academic year 21-211 Prof. Alberto Bemporad (University of Trento) Automatic Control 1 Academic year 21-211 1 /
More informationAutomatic Control 2. Nonlinear systems. Prof. Alberto Bemporad. University of Trento. Academic year
Automatic Control 2 Nonlinear systems Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011 1 / 18
More informationIdentification Methods for Structural Systems
Prof. Dr. Eleni Chatzi System Stability Fundamentals Overview System Stability Assume given a dynamic system with input u(t) and output x(t). The stability property of a dynamic system can be defined from
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 informationHomogeneous and particular LTI solutions
Homogeneous and particular LTI solutions Daniele Carnevale Dipartimento di Ing. Civile ed Ing. Informatica (DICII), University of Rome Tor Vergata Fondamenti di Automatica e Controlli Automatici A.A. 2014-2015
More informationDiscrete and continuous dynamic systems
Discrete and continuous dynamic systems Bounded input bounded output (BIBO) and asymptotic stability Continuous and discrete time linear time-invariant systems Katalin Hangos University of Pannonia Faculty
More informationOptimal control and estimation
Automatic Control 2 Optimal control and estimation Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011
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 informationIntroduction to Modern Control MT 2016
CDT Autonomous and Intelligent Machines & Systems Introduction to Modern Control MT 2016 Alessandro Abate Lecture 2 First-order ordinary differential equations (ODE) Solution of a linear ODE Hints to nonlinear
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationControl Systems Design
ELEC4410 Control Systems Design Lecture 13: Stability Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 13: Stability p.1/20 Outline Input-Output
More information16.30 Estimation and Control of Aerospace Systems
16.30 Estimation and Control of Aerospace Systems Topic 5 addendum: Signals and Systems Aeronautics and Astronautics Massachusetts Institute of Technology Fall 2010 (MIT) Topic 5 addendum: Signals, Systems
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 informationContinuous Dynamics Solving LTI state-space equations גרא וייס המחלקה למדעי המחשב אוניברסיטת בן-גוריון
Continuous Dynamics Solving LTI state-space equations גרא וייס המחלקה למדעי המחשב אוניברסיטת בן-גוריון 2 State Space Models For a causal system with m inputs u t R m and p outputs y t R p, an nth-order
More informationLecture: Sampling. Automatic Control 2. Sampling. Prof. Alberto Bemporad. University of Trento. Academic year
Automatic Control 2 Sampling Prof. Alberto Bemporad University of rento Academic year 2010-2011 Prof. Alberto Bemporad (University of rento) Automatic Control 2 Academic year 2010-2011 1 / 31 ime-discretization
More information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 7: State-space Models Readings: DDV, Chapters 7,8 Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology February 25, 2011 E. Frazzoli
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 informationECE504: Lecture 8. D. Richard Brown III. Worcester Polytechnic Institute. 28-Oct-2008
ECE504: Lecture 8 D. Richard Brown III Worcester Polytechnic Institute 28-Oct-2008 Worcester Polytechnic Institute D. Richard Brown III 28-Oct-2008 1 / 30 Lecture 8 Major Topics ECE504: Lecture 8 We are
More informationModelling and Control of Dynamic Systems. Stability of Linear Systems. Sven Laur University of Tartu
Modelling and Control of Dynamic Systems Stability of Linear Systems Sven Laur University of Tartu Motivating Example Naive open-loop control r[k] Controller Ĉ[z] u[k] ε 1 [k] System Ĝ[z] y[k] ε 2 [k]
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationControl System Design
ELEC4410 Control System Design Lecture 19: Feedback from Estimated States and Discrete-Time Control Design Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science
More informationAutomatique. A. Hably 1. Commande d un robot mobile. Automatique. A.Hably. Digital implementation
A. Hably 1 1 Gipsa-lab, Grenoble-INP ahmad.hably@grenoble-inp.fr Commande d un robot mobile (Gipsa-lab (DA)) ASI 1 / 25 Outline 1 2 (Gipsa-lab (DA)) ASI 2 / 25 of controllers Signals must be sampled and
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 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 informationControl Systems. Internal Stability - LTI systems. L. Lanari
Control Systems Internal Stability - LTI systems L. Lanari outline LTI systems: definitions conditions South stability criterion equilibrium points Nonlinear systems: equilibrium points examples stable
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016
Prof. Dr. Eleni Chatzi Lecture 4-09. March, 2016 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural
More informationELEC 3035, Lecture 3: Autonomous systems Ivan Markovsky
ELEC 3035, Lecture 3: Autonomous systems Ivan Markovsky Equilibrium points and linearization Eigenvalue decomposition and modal form State transition matrix and matrix exponential Stability ELEC 3035 (Part
More informationLinear dynamical systems with inputs & outputs
EE263 Autumn 215 S. Boyd and S. Lall Linear dynamical systems with inputs & outputs inputs & outputs: interpretations transfer function impulse and step responses examples 1 Inputs & outputs recall continuous-time
More informationEE263: Introduction to Linear Dynamical Systems Review Session 6
EE263: Introduction to Linear Dynamical Systems Review Session 6 Outline diagonalizability eigen decomposition theorem applications (modal forms, asymptotic growth rate) EE263 RS6 1 Diagonalizability consider
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 informationLinearization problem. The simplest example
Linear Systems Lecture 3 1 problem Consider a non-linear time-invariant system of the form ( ẋ(t f x(t u(t y(t g ( x(t u(t (1 such that x R n u R m y R p and Slide 1 A: f(xu f(xu g(xu and g(xu exist and
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 informationStabilisation of Linear Time-Invariant Systems subject to Output Saturation
Stabilisation of Linear Time-Invariant Systems subject to Output Saturation by Gijs Hilhorst, BSc A Thesis in Partial Fulfillment of the Requirement for the Degree of Master of Science at the Department
More informationIngegneria dell Automazione - Sistemi in Tempo Reale p.1/28
Ingegneria dell Automazione - Sistemi in Tempo Reale Selected topics on discrete-time and sampled-data systems Luigi Palopoli palopoli@sssup.it - Tel. 050/883444 Ingegneria dell Automazione - Sistemi in
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 informationEE363 homework 7 solutions
EE363 Prof. S. Boyd EE363 homework 7 solutions 1. Gain margin for a linear quadratic regulator. Let K be the optimal state feedback gain for the LQR problem with system ẋ = Ax + Bu, state cost matrix Q,
More informationRural/Urban Migration: The Dynamics of Eigenvectors
* Analysis of the Dynamic Structure of a System * Rural/Urban Migration: The Dynamics of Eigenvectors EGR 326 April 11, 2019 1. Develop the system model and create the Matlab/Simulink model 2. Plot and
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 informationFrequency domain analysis
Automatic Control 2 Frequency domain analysis Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011
More informationSampling of Linear Systems
Sampling of Linear Systems Real-Time Systems, Lecture 6 Karl-Erik Årzén January 26, 217 Lund University, Department of Automatic Control Lecture 6: Sampling of Linear Systems [IFAC PB Ch. 1, Ch. 2, and
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 informationMATH4406 (Control Theory) Unit 6: The Linear Quadratic Regulator (LQR) and Model Predictive Control (MPC) Prepared by Yoni Nazarathy, Artem
MATH4406 (Control Theory) Unit 6: The Linear Quadratic Regulator (LQR) and Model Predictive Control (MPC) Prepared by Yoni Nazarathy, Artem Pulemotov, September 12, 2012 Unit Outline Goal 1: Outline linear
More informationAC&ST AUTOMATIC CONTROL AND SYSTEM THEORY SYSTEMS AND MODELS. Claudio Melchiorri
C. Melchiorri (DEI) Automatic Control & System Theory 1 AUTOMATIC CONTROL AND SYSTEM THEORY SYSTEMS AND MODELS Claudio Melchiorri Dipartimento di Ingegneria dell Energia Elettrica e dell Informazione (DEI)
More informationDepartment of Electrical and Computer Engineering ECED4601 Digital Control System Lab3 Digital State Space Model
Department of Electrical and Computer Engineering ECED46 Digital Control System Lab3 Digital State Space Model Objectives. To learn some MATLAB commands that deals with the discrete time systems.. To give
More informationComplex Dynamic Systems: Qualitative vs Quantitative analysis
Complex Dynamic Systems: Qualitative vs Quantitative analysis Complex Dynamic Systems Chiara Mocenni Department of Information Engineering and Mathematics University of Siena (mocenni@diism.unisi.it) Dynamic
More informationTheory of Linear Systems Exercises. Luigi Palopoli and Daniele Fontanelli
Theory of Linear Systems Exercises Luigi Palopoli and Daniele Fontanelli Dipartimento di Ingegneria e Scienza dell Informazione Università di Trento Contents Chapter. Exercises on the Laplace Transform
More informationInput-to-state stability and interconnected Systems
10th Elgersburg School Day 1 Input-to-state stability and interconnected Systems Sergey Dashkovskiy Universität Würzburg Elgersburg, March 5, 2018 1/20 Introduction Consider Solution: ẋ := dx dt = ax,
More informationLinear System Theory. Wonhee Kim Lecture 1. March 7, 2018
Linear System Theory Wonhee Kim Lecture 1 March 7, 2018 1 / 22 Overview Course Information Prerequisites Course Outline What is Control Engineering? Examples of Control Systems Structure of Control Systems
More informationControl Systems. Time response
Control Systems Time response L. Lanari outline zero-state solution matrix exponential total response (sum of zero-state and zero-input responses) Dirac impulse impulse response change of coordinates (state)
More informationEE451/551: Digital Control. Chapter 7: State Space Representations
EE45/55: Digital Control Chapter 7: State Spae Representations State Variables Definition 7.: The system state is a minimal set of { variables x ( t ), i =, 2,, n needed together with the i input ut, and
More informationCONTROL OF DIGITAL SYSTEMS
AUTOMATIC CONTROL AND SYSTEM THEORY CONTROL OF DIGITAL SYSTEMS Gianluca Palli Dipartimento di Ingegneria dell Energia Elettrica e dell Informazione (DEI) Università di Bologna Email: gianluca.palli@unibo.it
More informationAutonomous system = system without inputs
Autonomous system = system without inputs State space representation B(A,C) = {y there is x, such that σx = Ax, y = Cx } x is the state, n := dim(x) is the state dimension, y is the output Polynomial representation
More informationEC744 Lecture Notes: Economic Dynamics. Prof. Jianjun Miao
EC744 Lecture Notes: Economic Dynamics Prof. Jianjun Miao 1 Deterministic Dynamic System State vector x t 2 R n State transition function x t = g x 0 ; t; ; x 0 = x 0 ; parameter 2 R p A parametrized dynamic
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 informationChapter III. Stability of Linear Systems
1 Chapter III Stability of Linear Systems 1. Stability and state transition matrix 2. Time-varying (non-autonomous) systems 3. Time-invariant systems 1 STABILITY AND STATE TRANSITION MATRIX 2 In this chapter,
More informationME Fall 2001, Fall 2002, Spring I/O Stability. Preliminaries: Vector and function norms
I/O Stability Preliminaries: Vector and function norms 1. Sup norms are used for vectors for simplicity: x = max i x i. Other norms are also okay 2. Induced matrix norms: let A R n n, (i stands for induced)
More informationAutomatic Control Systems theory overview (discrete time systems)
Automatic Control Systems theory overview (discrete time systems) Prof. Luca Bascetta (luca.bascetta@polimi.it) Politecnico di Milano Dipartimento di Elettronica, Informazione e Bioingegneria Motivations
More information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 8: Solutions of State-space Models Readings: DDV, Chapters 10, 11, 12 (skip the parts on transform methods) Emilio Frazzoli Aeronautics and Astronautics Massachusetts
More informationSTATE-SPACE DYNAMIC SYSTEMS (DISCRETE-TIME)
ECE452/552: Multivariable Control Systems I. 3 1 STATE-SPACE DYNAMIC SYSTEMS (DISCRETE-TIME) 3.1: The transform Computer control requires analog-to-digital (A2D) and digital-toanalog (D2A) conversion.
More informationLecture 2 and 3: Controllability of DT-LTI systems
1 Lecture 2 and 3: Controllability of DT-LTI systems Spring 2013 - EE 194, Advanced Control (Prof Khan) January 23 (Wed) and 28 (Mon), 2013 I LTI SYSTEMS Recall that continuous-time LTI systems can be
More informationLecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.
Lecture 4 Chapter 4: Lyapunov Stability Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 4 p. 1/86 Autonomous Systems Consider the autonomous system ẋ
More informationControl Systems. Time response. L. Lanari
Control Systems Time response L. Lanari outline zero-state solution matrix exponential total response (sum of zero-state and zero-input responses) Dirac impulse impulse response change of coordinates (state)
More informationESC794: Special Topics: Model Predictive Control
ESC794: Special Topics: Model Predictive Control Discrete-Time Systems Hanz Richter, Professor Mechanical Engineering Department Cleveland State University Discrete-Time vs. Sampled-Data Systems A continuous-time
More informationLinear ODEs. Existence of solutions to linear IVPs. Resolvent matrix. Autonomous linear systems
Linear ODEs p. 1 Linear ODEs Existence of solutions to linear IVPs Resolvent matrix Autonomous linear systems Linear ODEs Definition (Linear ODE) A linear ODE is a differential equation taking the form
More informationDynamic interpretation of eigenvectors
EE263 Autumn 2015 S. Boyd and S. Lall Dynamic interpretation of eigenvectors invariant sets complex eigenvectors & invariant planes left eigenvectors modal form discrete-time stability 1 Dynamic interpretation
More informationSymbolic Dynamics of Digital Signal Processing Systems
Symbolic Dynamics of Digital Signal Processing Systems Dr. Bingo Wing-Kuen Ling School of Engineering, University of Lincoln. Brayford Pool, Lincoln, Lincolnshire, LN6 7TS, United Kingdom. Email: wling@lincoln.ac.uk
More informationME 234, Lyapunov and Riccati Problems. 1. This problem is to recall some facts and formulae you already know. e Aτ BB e A τ dτ
ME 234, Lyapunov and Riccati Problems. This problem is to recall some facts and formulae you already know. (a) Let A and B be matrices of appropriate dimension. Show that (A, B) is controllable if and
More informationChapter 7. Canonical Forms. 7.1 Eigenvalues and Eigenvectors
Chapter 7 Canonical Forms 7.1 Eigenvalues and Eigenvectors Definition 7.1.1. Let V be a vector space over the field F and let T be a linear operator on V. An eigenvalue of T is a scalar λ F such that there
More informationWhy do we need a model?
Why do we need a model? For simulation (study the system output for a given input) Ex. Thermal study of a space shuttle when it enters the atmosphere (étude thermique concernant la rentrée dans l atmosphère
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 informationCDS 101 Precourse Phase Plane Analysis and Stability
CDS 101 Precourse Phase Plane Analysis and Stability Melvin Leok Control and Dynamical Systems California Institute of Technology Pasadena, CA, 26 September, 2002. mleok@cds.caltech.edu http://www.cds.caltech.edu/
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 informationDesigning Information Devices and Systems II Fall 2015 Note 22
EE 16B Designing Information Devices and Systems II Fall 2015 Note 22 Notes taken by John Noonan (11/12) Graphing of the State Solutions Open loop x(k + 1) = Ax(k) + Bu(k) y(k) = Cx(k) Closed loop x(k
More information9 Controller Discretization
9 Controller Discretization In most applications, a control system is implemented in a digital fashion on a computer. This implies that the measurements that are supplied to the control system must be
More informationLecture 4 Continuous time linear quadratic regulator
EE363 Winter 2008-09 Lecture 4 Continuous time linear quadratic regulator continuous-time LQR problem dynamic programming solution Hamiltonian system and two point boundary value problem infinite horizon
More informationLinear Systems. Manfred Morari Melanie Zeilinger. Institut für Automatik, ETH Zürich Institute for Dynamic Systems and Control, ETH Zürich
Linear Systems Manfred Morari Melanie Zeilinger Institut für Automatik, ETH Zürich Institute for Dynamic Systems and Control, ETH Zürich Spring Semester 2016 Linear Systems M. Morari, M. Zeilinger - Spring
More informationEntrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4.
Entrance Exam, Differential Equations April, 7 (Solve exactly 6 out of the 8 problems). Consider the following initial value problem: { y + y + y cos(x y) =, y() = y. Find all the values y such that the
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 informationControl Systems Lab - SC4070 Control techniques
Control Systems Lab - SC4070 Control techniques Dr. Manuel Mazo Jr. Delft Center for Systems and Control (TU Delft) m.mazo@tudelft.nl Tel.:015-2788131 TU Delft, February 16, 2015 (slides modified from
More informationData Driven Discrete Time Modeling of Continuous Time Nonlinear Systems. Problems, Challenges, Success Stories. Johan Schoukens
1/51 Data Driven Discrete Time Modeling of Continuous Time Nonlinear Systems Problems, Challenges, Success Stories Johan Schoukens fyuq (,, ) 4/51 System Identification Data Distance Model 5/51 System
More informationLogistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations
Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations S. Y. Ha and J. Park Department of Mathematical Sciences Seoul National University Sep 23, 2013 Contents 1 Logistic Map 2 Euler and
More information21 Linear State-Space Representations
ME 132, Spring 25, UC Berkeley, A Packard 187 21 Linear State-Space Representations First, let s describe the most general type of dynamic system that we will consider/encounter in this class Systems may
More informationECE 602 Exam 2 Solutions, 3/23/2011.
NAME: ECE 62 Exam 2 Solutions, 3/23/211. You can use any books or paper notes you wish to bring. No electronic devices of any kind are allowed. You can only use materials that you bring yourself. You are
More informationNonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems. p. 1/1
Nonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems p. 1/1 p. 2/1 Converse Lyapunov Theorem Exponential Stability Let x = 0 be an exponentially stable equilibrium
More informationDifferential Equations
Differential Equations Overview of differential equation! Initial value problem! Explicit numeric methods! Implicit numeric methods! Modular implementation Physics-based simulation An algorithm that
More informationAutomatic Control 2. Loop shaping. Prof. Alberto Bemporad. University of Trento. Academic year
Automatic Control 2 Loop shaping Prof. Alberto Bemporad University of Trento Academic year 21-211 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 21-211 1 / 39 Feedback
More information2 Lyapunov Stability. x(0) x 0 < δ x(t) x 0 < ɛ
1 2 Lyapunov Stability Whereas I/O stability is concerned with the effect of inputs on outputs, Lyapunov stability deals with unforced systems: ẋ = f(x, t) (1) where x R n, t R +, and f : R n R + R n.
More informationProperties of Open-Loop Controllers
Properties of Open-Loop Controllers Sven Laur University of Tarty 1 Basics of Open-Loop Controller Design Two most common tasks in controller design is regulation and signal tracking. Regulating controllers
More informationPoincaré Map, Floquet Theory, and Stability of Periodic Orbits
Poincaré Map, Floquet Theory, and Stability of Periodic Orbits CDS140A Lecturer: W.S. Koon Fall, 2006 1 Poincaré Maps Definition (Poincaré Map): Consider ẋ = f(x) with periodic solution x(t). Construct
More informationSolving Dynamic Equations: The State Transition Matrix
Overview Solving Dynamic Equations: The State Transition Matrix EGR 326 February 24, 2017 Solutions to coupled dynamic equations Solutions to dynamic circuits from EGR 220 The state transition matrix Discrete
More informationLyapunov Stability Theory
Lyapunov Stability Theory Peter Al Hokayem and Eduardo Gallestey March 16, 2015 1 Introduction In this lecture we consider the stability of equilibrium points of autonomous nonlinear systems, both in continuous
More informationIntroduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems
p. 1/5 Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems p. 2/5 Time-varying Systems ẋ = f(t, x) f(t, x) is piecewise continuous in t and locally Lipschitz in x for all t
More informationDepartment of Mathematics IIT Guwahati
Stability of Linear Systems in R 2 Department of Mathematics IIT Guwahati A system of first order differential equations is called autonomous if the system can be written in the form dx 1 dt = g 1(x 1,
More information4F3 - Predictive Control
4F3 Predictive Control - Lecture 2 p 1/23 4F3 - Predictive Control Lecture 2 - Unconstrained Predictive Control Jan Maciejowski jmm@engcamacuk 4F3 Predictive Control - Lecture 2 p 2/23 References Predictive
More informationAutonomous Systems and Stability
LECTURE 8 Autonomous Systems and Stability An autonomous system is a system of ordinary differential equations of the form 1 1 ( 1 ) 2 2 ( 1 ). ( 1 ) or, in vector notation, x 0 F (x) That is to say, an
More informationNonlinear Control Lecture # 14 Input-Output Stability. Nonlinear Control
Nonlinear Control Lecture # 14 Input-Output Stability L Stability Input-Output Models: y = Hu u(t) is a piecewise continuous function of t and belongs to a linear space of signals The space of bounded
More informationReview: control, feedback, etc. Today s topic: state-space models of systems; linearization
Plan of the Lecture Review: control, feedback, etc Today s topic: state-space models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered
More information