Some solutions of the written exam of January 27th, 2014
|
|
- Christian Daniels
- 5 years ago
- Views:
Transcription
1 TEORIA DEI SISTEMI Systems Theory) Prof. C. Manes, Prof. A. Germani Some solutions of the written exam of January 7th, 0 Problem. Consider a feedback control system with unit feedback gain, with the following transfer function in open-loop Ws) = K s+) s +).. Draw the amplitude and phase Bode diagrams, and the polar diagram for K = ;. Compute the denominator of the closed-loop transfer function; 3. Compute the number of poles with negative real part of the closed loop transfer function as a function of the gain K, ), using both the Nyquist criterion and the Routh criterion. Solution of problem. Let Ws) denote Ws) for K = : Ws) = s+) s +) We consider the Bode plots and the polar plot of Ws). In Bode form we have Ws) = K Ws) = K ) ) +s) + s The low frequency gain is K 0 = lim s 0 Ws) = K 0 db = 0log 0 K 0 = 0log 0 = 0log 0 = 6 db. Thus Wjω) = +jω +jω/) = +jω ω / Wjω) = +jω ω /. Note that the term ω / is real, and therefore its phase is 0 when ω < because ω / > 0) and is π when ω > because ω / < 0). the Bode plots and the Nyquist plot of the open loop transfer function are in the enclosed file). The closed-loop transfer function is s+) s +) W CL s) = K Ws) +K Ws) = K +K s+) s +) = K s+) s +) K and the denominator of W CL s) is: d CL s) = s+) s +) K = s +s+)s +) K = s +s 3 +5s +8s+ K NYQUIST ANALYSIS Let N be the number of times that the Nyquist plot of Wjω) encircles the point in the counterclockwise direction. From the plot it is clear that For K 0,) we have N = 0: the Nyquist plot does not encircle the point ; For K > we have N = : the Nyquist plot encircles one time the point in the clockwise negative) direction; For K < 0 we have N = : the Nyquist plot encircles two times the point in the clockwise negative) direction;
2 Recall the Nyquist formula in the form p CL = p OL N where, p CL is the number of poles with positive real part of the Closed Loop CL) system, and p OL is the number of poles with positive real part of the Open Loop OL) system. Since for the given Ws) we have p OL = 0 there isn t any unstable pole in the open-loop system) we have p CL = N, and therefore For K 0,) we have N = 0, and thus p CL = 0 stable closed loop system); For K > we have N =, and thus p CL = unstable closed loop system); For K < 0 we have N =, and thus p CL = unstable closed loop system). Note that for K = the denominator of W CL s) is 0 for s = 0, and thus s = 0 is a pole of W CL s). ROUTH ANALYSIS The case K = 0 will be not analyzed because it corresponds to the trivial case where the open-loop transfer function is zero. The characteristic polynomial of the closed-loop system is the denominator of W CL s): d CL s) = s +s 3 +5s +8s+ K The first two rows rows and 3) of the Routh table are: 5 K 3 8 The computation of the elements in the third row row number ) gives a, = 5 8 = 8 0 =, a, = -K 0 = K) = K Thus we have 5 K 3 8 K The computation of the element in the fourth row row number ) gives a, = 8 K = K) 8 = K, and the last element, a 0, is K 5 K 3 8 K K 0 K Analyzing the signs of the first column we have: For K 0,) we have no sign variation, so that p CL = 0; For K > we have one sign variation 0), so that p CL = ; For K < we have two sign variation and 0) so that p CL =. For the particular case of K =, the characteristic polynomial of the closed-loop system is d CL s) == s +s 3 +5s +8s = ss 3 +s +5s+8), thus, we have a pole in the origin. The sign of the remaining poles can be studied again by means of the Routh criterion:
3 There is no sign variation in the first column, and therefore the remaining poles have all negative real part simple stability of the origin of the state space). These results coincide with those obtained using the Nyquist analysis.
4 Problem. Given the following discrete-time system [ ] xt+) = Axt)+But) dove A = yt) = Cxt), discuss the properties of the natural modes;. compute the state-transition matrix A t ; 3. compute the input-output impulse response and transfer function B = [ ] 0 C = [ 0 ] Solution of problem. Let s start computing the eigenvalues of A, which are the roots of the characteristic polynomial of A. The computation of the characteristic polynomial gives A = [ ] λi A = λ ) + ) The eigenvalues of A are computed solving λi A = 0. Instead of using the classical formula for second-degree equations, we can proceed as follows λ ) + ) = 0 From this we have the two complex-conjugates eigenvalues: λ ) = λ ) = ±. λ = +j, λ = j We denote r C and l C the right and left eigenvalues associated to λ. Their computation is made by solving the two systems of equation λ I A)r = 0, l λ I A) = 0, and taking, among the infinite solutions, one solution such that l r =. The right and left eigenvalues associated to λ are the conjugates of r and l so that λ = λ, r = r and l = l ). We choose the following set of eigenvalues λ = +j r = [ ], l j = [ ] j, λ = j r = [ ], l j = [ ] j. check that l r = l r = and l r = l r = 0). The properties of the natural modes to be discussed are the following: - Stability, Eccitability, Observability. Both modes are asymptotically stable, because they satisfy the stability condition λ k < note that the system is a discrete-time system). In particular ) λ = λ = + ) = 8 = < Moreover, both natural modes are eccitable because l B = j/ 0 and l B = j/ 0, and are observable because Cr = Cr = 0. The transition matrix can be computed using the spectral decomposition A t = λ t r l +λ t r l = R λ t r l ) The expression A t = R +j ) t [ ] [ ]) j = R j +j ) t [ ]) j j is not useful if we can not expand the power +j ) t. Thus, it is important to write λ using the magnitude and phase form Euler formula) λ = λ e j λ = λ cos λ )+j sin λ ))
5 so that λ t = λ t e j λ t = λ t cos λ t)+j sin λ t)) The magnitude and phase of λ = +j are λ = / 8 and λ = arctan) = π. Thus ) t ) t λ t = e j π t = 8 cos π ) 8 t)+j sinπ t) Thus the transition matrix A t becomes ) t A t = R cos π ) [ ] ) j 8 t)+j sinπ t) j ) t = R cos π )[ ]) j 8 t)+j sinπ t). j Taking the real part we have ) [ t cos π A t = t) sin π t) ] 8 sin π t) cos π t) The impulse response is w0) = D and wt) = CA t B for t > 0. Then we have w0) = 0 and fort > 0 wt) = [ 0 ] ) [ t cos π 8 t )) sin π t )) ] [0 ] ) t π ) sin π t )) cos π t )) = sin 8 t ). The transfer function Wz) can be computed in two ways: ) by taking the Z-transform of the impulse response wt); ) by using the formula Wz) = CzI A) B +D, with D = 0. Using the second method we easily get Wz) = z ) + = z ) + ) [ 0 ] [ z z ) ] [0 ] As an exercise you can compute Wz) by transforming wt). Of course, you should obtain the same result! I suggest to develop the Z-transform of a generic function ft) = ca t sin bt ) and then to apply the delay theorem: Consider a function ft), t 0, and its Z-transform Fz) = Zft)). Let gt) = ft ). Then, the Z-transform of gt) is such that Gz) = z Fz). When ft) = ca t sin bt ), then gt) = ca t sin bt ) ). Notice that the impulse response wt) has the same structure of gt) when c =, a = / 8 and b = π/. Thus, the computation of Wz) from Gz) and from Fz) is straightforward. The computation of the Z-transform Fz) is easy if we rewrite the sinus in the exponential form: ft) = ca t e jbt e jbt )/j, so that ft) = c ae jb ) t ae jb ) t) /j. The only transformation we need to remember is Zp t ) = z/z p). From this, Z ft) ) = c ) z j z ae z jb z ae jb Replacing c =, a = / 8 and b = π/, and noting that ae jb = / 8)e jπ/ = / 8)+j)/ = +j)/, after few computation we get the expression for Wz).
6 Problem 3. Consider a continuous-time system characterized by the following impulse response: wt) = 3e t. Compute the forced response and the harmonic response to the input ut) = 5 sint). Solution of problem 3. Both responses are easily computed in the Laplace domain. First of all we need to compute the transfer function of the system by transforming the impulse response: Ws) = L wt) ). Recalling that L e at) = /s a) we get Ws) = 3 s+. Forced response. The forced response in the Laplace domain is easily computed as Ys) = Ws)Us), where Us) is the Laplace transform of the input function. Recalling that L c sinωt) ) = cω/s +ω ), we have and then Ys) = Ws)Us) = Us) = 0 s + 30 s+)s +) = 30 s+)s j)s+j) Now must compute the partial fraction decomposition of Ys), and we can choose one of the two forms: or Ys) = R s+) + R s j) + R 3 s+j) Ys) = R s+) + as+b s +) recall that R 3 = R )...to be continued... the back-transformation of Ys) gives the forced response yt) =......to be computed... Harmonic response. The harmonic response exists if an only if the system is asymptotically stable or at least the controllable and observable subsystem should be asymptotically stable). This means that all poles of the transfer function that are also the eigenvalues of the controllable and observable subsystem), must have negative part. Since this condition is verified in this problem the pole is -), then the harmonic response exists and we can proceed in its computation. The general formula for the harmonic response of continuous-time systems to an input of the form ut) = c sinωt) is y h t) = c Wjω) sinωt+ Wjω) ). In our problem Ws) = 3/s+), and ω =. Thus we need to compute Wj) and Wj) ): Wjω) = 3 Wjω) = jω + = 3 ω + 3 ω = 3 jω + = jω + = arctan. jω + ) From these, when ω = we have 3 Wj) = = 3 ) +, Wj) = arctan = π. Thus, the harmonic response to the input ut) = 5 sint) is y h t) = 5 Wj) sint+ Wj) ) = 5 sint π/)
7 Problem. Given the system xt+) = Axt)+But), yt) = Cxt) with matrices 0 0 A = B = 0 C = [ 0 0 ]. Find a basis for the space of reachable states, and a basis for the space of unobservable states;. Find the subspaces X, X, X 3 and X of the Kalman structural decomposition; 3. Find an input sequence that brings the state from x0) = [0 0 0] T to x a = [ 0] T.. Find an input sequence that brings the state from x0) = [0 0 0] T to x b = [ 0 ] T. Solution of problem. The computation of the reachability matrix gives P 3 = [ B AB A B ] 0 0 = 0 0 P = RP 3 ) the range of P 3 ) is the space of reachable states. Note that only the first two columns of P 3 are independent, so the rank of P 3 is, and is the dimension of P. The first two columns of P 3 are two admissible basis vectors: 0 P = R[B AB A B]) = R[B AB]) = span 0, 0, v = 0, v = The computation of the observability matrix gives C 0 0 Q 3 = CA = CA 0 0 Q = NQ 3 ) the null-space of Q 3 ) is the space of unobservable states. Note that only the first two rows of Q 3 are independent, so the rank of Q 3 is, and is the dimension of Q is the rank drop of Q 3 ). Thus, we need to find one basis vector for Q, that is a solution of the homogeneous system Q 3 b = 0 3. Observing that the first two columns of Q 3 are equal, we have the following solution b = [ 0] T, which is a basis for Q. Thus C Q = N CA = N CA [ ]) C CA = span, b = 0 Now we are ready to find the four subspaces X, X, X 3 and X. Let s start with X = Q P. Being Q a one-dimensional space and P a two-dimensional space, the intersections is either 0 or Q itself. Thus, we must only check if b spanv,v ) or not. This can be made by checking whether rank[v v ]) = rank[v v b]) or not. It is easy to check that 0 0 rank 0 = rank 0 0 In other words, b is linearly dependent on v,v it is clear that b = v +v ). Thus, b is a basis for X. The space X is any space such that X X = P. We can choose either X = spanv ) or X = spanv ), because we have P = span[b v )], and also P = span[b v )]. Let s choose X = spanv ), so that a basis for X is v = [0 0 ] T. The space X 3 is any space such that X X 3 = Q. However, since we have X = Q, then X 3 =. The space X is any space such that X X X 3 X = C 3. Then, we only need to find a vector linearly independent of b,v. Let s choose X = spand), with d = [ 0 0] T, such that X =spanb), rank [ b v d ]) 0 X = rank 0 0 =spanv ), = 3, and then 0 0 X 3 =, X =spand).. 0.
8 Note that other choices for X and X are possible. Question. and.: - Find an input sequence that brings the state from x0) = [0 0 0] T to x a = [ 0] T ; - Find an input sequence that brings the state from x0) = [0 0 0] T to x b = [ 0 ] T. Note that x a P, and therefore we know that there exists an input sequence that brings x0) = 0 3 to x a. From the general explicit form of a discrete-time system: ut ) t xt) = A t x0)+ A t τ Buτ) = A t x0)+ [ B AB A t B ]. τ=0 u). u0) We know that the matrix P 3 transforms an input sequence [u0) u) u)] T to x3), starting with x0): x3) = [ B AB A B ] u) u). u0) Thus, if a vector x P 3, then there exists a triple u0),u),u) ) that transfers x0) = 0 to x. However, since only the first two columns of P 3 are independent, when x P, then we can also look for a pair u0),u) ) that transfers x0) = 0 to x, by solving x = [ B AB ][ ] u). u0) We see that x a P x a is reachable) while x b P x b is not reachable). Thus, for x a = [ 0] T we can find a pair u0),u) ) that transfers x0) = 0 to x a. We only have to solve: x a = [ B AB ][ ] [ ] 0 u) u) = 0. u0) 0 u0) Solving this vector equation it is easy to see that u0) = and u) = bring the state x0) = 0 to the state x) = x a.
9 Problem 5. Given the system {ẋ t) = α)x t)+αx t)x t) ẋ t) = α x t) x t). study the stability of the equilibrium point x e = 0,0) for all the values of the parameter α, ), using the method of linear approximation at the equilibrium point, and the Lyapunov function method;. study for what values of α the system admits other equilibrium points, and study their stability properties. Suggestion for the Lyapunov function: Vx) = x x e, ) +βx x e, ), with suitable β > 0.) Solution of problem 5. The system considered is of the form ẋt) = f xt);α ), where xt) is the state and α is a constant parameter. The vector function fx;α) is as follows { f x;α) = α)x +αx x The Jacobian is The Jacobian computed at x e = 0,0) is f x;α) = α x x [ ] [ ] x f x f α+αx αx Jx) = = x f x f x α Jx e ) = [ α ] 0 0 α Since Jx e ) is diagonal, the two eigenvalues are the diagonal terms: Thus: λ = α, λ = α. 0,0) is asymptotically stable when both eigenvalues are strictly negative, that is when α < 0) and α < 0). That means that 0,0) is asymptotically stable for α > ; 0,0) is unstable when at least one eigenvalue is strictly positive i.e. when α > 0) or α > 0). That means that 0,0) is unstable for α <. The only case where we can not conclude anything about the stability of 0,0) is when α = in this case the origin is a simply stable equilibrium point of the linear approximation of the nonlinear system, but nothing can be concluded about the nonlinear system). In this case we must study the stability of the point x e = 0,0) by using a suitable Lyapunov function. For α = the function fx;α) is { f x;) = x x f x;) = x x. Following the suggestion we consider the Lyapunov function Vx) = x +βx, which is positive definite for any β > 0. According to the Lyapunov theorem, if the derivative Vx) = dv/dx)fx;) is semidefinite negative, then the equilibrium x e is simply) stable, while if Vx) is definite negative, then the equilibrium is asymptotically stable. The computation of Vx) gives Vx) = dv dx fx;) = x V)f x;)+ x V)f x;) = x f x;)+βx f x;) = x x βx βx x.
10 By choosing β = we can eliminate the term x x which is sign-indefinite, so that we get Vx) = x. which is semidefinite negative it is never positive, but it is zero for any pair x,0), and not only at the equilibrium point x e = 0,0)). As a consequence, when α = the equilibrium 0,0) is simply stable stable, but not asymptotically stable). All other equilibrium points can be found by solving the system of equations fx;α) = 0, i.e. { α)x +αx x = 0 From the first equation we have α x x = 0. x α)+αx ) = 0 α+αx = 0 x = α α, α 0 the other solution is x = 0, and leads to the equilibrium point 0,0) already studied in the previous question). From the second equation we have α x x = 0 α α α x = 0 x = α α) Thus, there exists other equilibrium points in R only if α 0 and α α) > 0, that is for α 0,). In this case we have two more equilibrium points: ) α α), α x e, =, x e, = α α), α ). α α Note that when α = all the equilibrium points collapse to 0,0). The Jacobian computed in the two equilibrium points is [ 0 ±α ] α α) Jx e,, ) = x e,, = ± α α), α ). α α) α α The characteristic polynomial is λi Jx e,, ) [ λ α ] α α) = ± α α) λ+α = λ +α λ+α α). Note that the two Jacobians have the same characteristic polynomial. Recall that the roots of a second order polynomial all have negative real part if and only if the three coefficients have the same sign check this property using the Routh table). Thus, both equilibrium points are asymptotically stable, because they exists if and only if α 0 and α 0,), which imply that α > 0 and α α) > 0.
Theory 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 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 information1 (s + 3)(s + 2)(s + a) G(s) = C(s) = K P + K I
MAE 43B Linear Control Prof. M. Krstic FINAL June 9, Problem. ( points) Consider a plant in feedback with the PI controller G(s) = (s + 3)(s + )(s + a) C(s) = K P + K I s. (a) (4 points) For a given constant
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 I. Lecture 6: Poles and Zeros. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich
Control Systems I Lecture 6: Poles and Zeros Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 27, 2017 E. Frazzoli (ETH) Lecture 6: Control Systems I 27/10/2017
More informationFull-State Feedback Design for a Multi-Input System
Full-State Feedback Design for a Multi-Input System A. Introduction The open-loop system is described by the following state space model. x(t) = Ax(t)+Bu(t), y(t) =Cx(t)+Du(t) () 4 8.5 A =, B =.5.5, C
More informationControl Systems I. Lecture 5: Transfer Functions. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich
Control Systems I Lecture 5: Transfer Functions Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 20, 2017 E. Frazzoli (ETH) Lecture 5: Control Systems I 20/10/2017
More informationÜbersetzungshilfe / Translation aid (English) To be returned at the end of the exam!
Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 5. 2. 2 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid -
More informationTime Response Analysis (Part II)
Time Response Analysis (Part II). A critically damped, continuous-time, second order system, when sampled, will have (in Z domain) (a) A simple pole (b) Double pole on real axis (c) Double pole on imaginary
More informationAndrea Zanchettin Automatic Control AUTOMATIC CONTROL. Andrea M. Zanchettin, PhD Spring Semester, Linear systems (frequency domain)
1 AUTOMATIC CONTROL Andrea M. Zanchettin, PhD Spring Semester, 2018 Linear systems (frequency domain) 2 Motivations Consider an LTI system Thanks to the Lagrange s formula we can compute the motion of
More informationControl Systems I. Lecture 7: Feedback and the Root Locus method. Readings: Jacopo Tani. Institute for Dynamic Systems and Control D-MAVT ETH Zürich
Control Systems I Lecture 7: Feedback and the Root Locus method Readings: Jacopo Tani Institute for Dynamic Systems and Control D-MAVT ETH Zürich November 2, 2018 J. Tani, E. Frazzoli (ETH) Lecture 7:
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 informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Illinois Institute of Technology Lecture 22: The Nyquist Criterion Overview In this Lecture, you will learn: Complex Analysis The Argument Principle The Contour
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 informationThe stability of linear time-invariant feedback systems
The stability of linear time-invariant feedbac systems A. Theory The system is atrictly stable if The degree of the numerator of H(s) (H(z)) the degree of the denominator of H(s) (H(z)) and/or The poles
More informationLecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design.
ISS0031 Modeling and Identification Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design. Aleksei Tepljakov, Ph.D. September 30, 2015 Linear Dynamic Systems Definition
More informationAnalysis of Discrete-Time Systems
TU Berlin Discrete-Time Control Systems 1 Analysis of Discrete-Time Systems Overview Stability Sensitivity and Robustness Controllability, Reachability, Observability, and Detectabiliy TU Berlin Discrete-Time
More informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Illinois Institute of Technology Lecture 23: Drawing The Nyquist Plot Overview In this Lecture, you will learn: Review of Nyquist Drawing the Nyquist Plot Using
More informationSTABILITY OF CLOSED-LOOP CONTOL SYSTEMS
CHBE320 LECTURE X STABILITY OF CLOSED-LOOP CONTOL SYSTEMS Professor Dae Ryook Yang Spring 2018 Dept. of Chemical and Biological Engineering 10-1 Road Map of the Lecture X Stability of closed-loop control
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 informationClass 13 Frequency domain analysis
Class 13 Frequency domain analysis The frequency response is the output of the system in steady state when the input of the system is sinusoidal Methods of system analysis by the frequency response, as
More informationSolutions of Spring 2008 Final Exam
Solutions of Spring 008 Final Exam 1. (a) The isocline for slope 0 is the pair of straight lines y = ±x. The direction field along these lines is flat. The isocline for slope is the hyperbola on the left
More informationFrequency methods for the analysis of feedback systems. Lecture 6. Loop analysis of feedback systems. Nyquist approach to study stability
Lecture 6. Loop analysis of feedback systems 1. Motivation 2. Graphical representation of frequency response: Bode and Nyquist curves 3. Nyquist stability theorem 4. Stability margins Frequency methods
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 informationAPPM 2360: Midterm exam 3 April 19, 2017
APPM 36: Midterm exam 3 April 19, 17 On the front of your Bluebook write: (1) your name, () your instructor s name, (3) your lecture section number and (4) a grading table. Text books, class notes, cell
More informationStabilization and Passivity-Based Control
DISC Systems and Control Theory of Nonlinear Systems, 2010 1 Stabilization and Passivity-Based Control Lecture 8 Nonlinear Dynamical Control Systems, Chapter 10, plus handout from R. Sepulchre, Constructive
More informationAnalysis of Discrete-Time Systems
TU Berlin Discrete-Time Control Systems TU Berlin Discrete-Time Control Systems 2 Stability Definitions We define stability first with respect to changes in the initial conditions Analysis of Discrete-Time
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 informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Arizona State University Lecture 23: Drawing The Nyquist Plot Overview In this Lecture, you will learn: Review of Nyquist Drawing the Nyquist Plot Using the
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. 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 informationr + - FINAL June 12, 2012 MAE 143B Linear Control Prof. M. Krstic
MAE 43B Linear Control Prof. M. Krstic FINAL June, One sheet of hand-written notes (two pages). Present your reasoning and calculations clearly. Inconsistent etchings will not be graded. Write answers
More informationLecture Note #6 (Chap.10)
System Modeling and Identification Lecture Note #6 (Chap.) CBE 7 Korea University Prof. Dae Ryoo Yang Chap. Model Approximation Model approximation Simplification, approximation and order reduction of
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 informationControl Systems I. Lecture 9: The Nyquist condition
Control Systems I Lecture 9: The Nyquist condition Readings: Åstrom and Murray, Chapter 9.1 4 www.cds.caltech.edu/~murray/amwiki/index.php/first_edition Jacopo Tani Institute for Dynamic Systems and Control
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 informationMath 215/255 Final Exam (Dec 2005)
Exam (Dec 2005) Last Student #: First name: Signature: Circle your section #: Burggraf=0, Peterson=02, Khadra=03, Burghelea=04, Li=05 I have read and understood the instructions below: Please sign: Instructions:.
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 informationExamples include: (a) the Lorenz system for climate and weather modeling (b) the Hodgkin-Huxley system for neuron modeling
1 Introduction Many natural processes can be viewed as dynamical systems, where the system is represented by a set of state variables and its evolution governed by a set of differential equations. Examples
More informationCopyright (c) 2006 Warren Weckesser
2.2. PLANAR LINEAR SYSTEMS 3 2.2. Planar Linear Systems We consider the linear system of two first order differential equations or equivalently, = ax + by (2.7) dy = cx + dy [ d x x = A x, where x =, and
More informationECEN 605 LINEAR SYSTEMS. Lecture 20 Characteristics of Feedback Control Systems II Feedback and Stability 1/27
1/27 ECEN 605 LINEAR SYSTEMS Lecture 20 Characteristics of Feedback Control Systems II Feedback and Stability Feedback System Consider the feedback system u + G ol (s) y Figure 1: A unity feedback system
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 informationECE317 : Feedback and Control
ECE317 : Feedback and Control Lecture : Routh-Hurwitz stability criterion Examples Dr. Richard Tymerski Dept. of Electrical and Computer Engineering Portland State University 1 Course roadmap Modeling
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 informationLinear Differential Equations. Problems
Chapter 1 Linear Differential Equations. Problems 1.1 Introduction 1.1.1 Show that the function ϕ : R R, given by the expression ϕ(t) = 2e 3t for all t R, is a solution of the Initial Value Problem x =
More informationLecture Note #6 (Chap.10)
System Modeling and Identification Lecture Note #6 (Chap.) CBE 7 Korea University Prof. Dae Ryook Yang Chap. Model Approximation Model approximation Simplification, approximation and order reduction of
More informationSolutions to Math 53 Math 53 Practice Final
Solutions to Math 5 Math 5 Practice Final 20 points Consider the initial value problem y t 4yt = te t with y 0 = and y0 = 0 a 8 points Find the Laplace transform of the solution of this IVP b 8 points
More informationProblem 2 (Gaussian Elimination, Fundamental Spaces, Least Squares, Minimum Norm) Consider the following linear algebraic system of equations:
EEE58 Exam, Fall 6 AA Rodriguez Rules: Closed notes/books, No calculators permitted, open minds GWC 35, 965-37 Problem (Dynamic Augmentation: State Space Representation) Consider a dynamical system consisting
More informationSTABILITY ANALYSIS. Asystemmaybe stable, neutrallyormarginallystable, or unstable. This can be illustrated using cones: Stable Neutral Unstable
ECE4510/5510: Feedback Control Systems. 5 1 STABILITY ANALYSIS 5.1: Bounded-input bounded-output (BIBO) stability Asystemmaybe stable, neutrallyormarginallystable, or unstable. This can be illustrated
More information1 The Observability Canonical Form
NONLINEAR OBSERVERS AND SEPARATION PRINCIPLE 1 The Observability Canonical Form In this Chapter we discuss the design of observers for nonlinear systems modelled by equations of the form ẋ = f(x, u) (1)
More information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 24: H2 Synthesis Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology May 4, 2011 E. Frazzoli (MIT) Lecture 24: H 2 Synthesis May
More informationReturn Difference Function and Closed-Loop Roots Single-Input/Single-Output Control Systems
Spectral Properties of Linear- Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 2018! Stability margins of single-input/singleoutput (SISO) systems! Characterizations
More informationSYSTEMTEORI - ÖVNING Stability of linear systems Exercise 3.1 (LTI system). Consider the following matrix:
SYSTEMTEORI - ÖVNING 3 1. Stability of linear systems Exercise 3.1 (LTI system. Consider the following matrix: ( A = 2 1 Use the Lyapunov method to determine if A is a stability matrix: a: in continuous
More informationDesign Methods for Control Systems
Design Methods for Control Systems Maarten Steinbuch TU/e Gjerrit Meinsma UT Dutch Institute of Systems and Control Winter term 2002-2003 Schedule November 25 MSt December 2 MSt Homework # 1 December 9
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 informationExam. 135 minutes + 15 minutes reading time
Exam January 23, 27 Control Systems I (5-59-L) Prof. Emilio Frazzoli Exam Exam Duration: 35 minutes + 5 minutes reading time Number of Problems: 45 Number of Points: 53 Permitted aids: Important: 4 pages
More information1 (20 pts) Nyquist Exercise
EE C128 / ME134 Problem Set 6 Solution Fall 2011 1 (20 pts) Nyquist Exercise Consider a close loop system with unity feedback. For each G(s), hand sketch the Nyquist diagram, determine Z = P N, algebraically
More informationÜbersetzungshilfe / Translation aid (English) To be returned at the end of the exam!
Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 3.. 24 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid -
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 informationCME 345: MODEL REDUCTION
CME 345: MODEL REDUCTION Balanced Truncation Charbel Farhat & David Amsallem Stanford University cfarhat@stanford.edu These slides are based on the recommended textbook: A.C. Antoulas, Approximation of
More informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Arizona State University Lecture 21: Stability Margins and Closing the Loop Overview In this Lecture, you will learn: Closing the Loop Effect on Bode Plot Effect
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 informationSolution for Homework 5
Solution for Homework 5 ME243A/ECE23A Fall 27 Exercise 1 The computation of the reachable subspace in continuous time can be handled easily introducing the concepts of inner product, orthogonal complement
More informationsc Control Systems Design Q.1, Sem.1, Ac. Yr. 2010/11
sc46 - Control Systems Design Q Sem Ac Yr / Mock Exam originally given November 5 9 Notes: Please be reminded that only an A4 paper with formulas may be used during the exam no other material is to be
More informationTime Response of Systems
Chapter 0 Time Response of Systems 0. Some Standard Time Responses Let us try to get some impulse time responses just by inspection: Poles F (s) f(t) s-plane Time response p =0 s p =0,p 2 =0 s 2 t p =
More informationThe goal of this chapter is to study linear systems of ordinary differential equations: dt,..., dx ) T
1 1 Linear Systems The goal of this chapter is to study linear systems of ordinary differential equations: ẋ = Ax, x(0) = x 0, (1) where x R n, A is an n n matrix and ẋ = dx ( dt = dx1 dt,..., dx ) T n.
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 informationExercise 8: Frequency Response of MIMO Systems
Exercise 8: Frequency Response of MIMO Systems 8 Singular Value Decomposition (SVD The Singular Value Decomposition plays a central role in MIMO frequency response analysis Let s recall some concepts from
More informationAbstract Vector Spaces
CHAPTER 1 Abstract Vector Spaces 1.1 Vector Spaces Let K be a field, i.e. a number system where you can add, subtract, multiply and divide. In this course we will take K to be R, C or Q. Definition 1.1.
More informationThird In-Class Exam Solutions Math 246, Professor David Levermore Thursday, 3 December 2009 (1) [6] Given that 2 is an eigenvalue of the matrix
Third In-Class Exam Solutions Math 26, Professor David Levermore Thursday, December 2009 ) [6] Given that 2 is an eigenvalue of the matrix A 2, 0 find all the eigenvectors of A associated with 2. Solution.
More informationDigital Control Systems
Digital Control Systems Lecture Summary #4 This summary discussed some graphical methods their use to determine the stability the stability margins of closed loop systems. A. Nyquist criterion Nyquist
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 informationMAS107 Control Theory Exam Solutions 2008
MAS07 CONTROL THEORY. HOVLAND: EXAM SOLUTION 2008 MAS07 Control Theory Exam Solutions 2008 Geir Hovland, Mechatronics Group, Grimstad, Norway June 30, 2008 C. Repeat question B, but plot the phase curve
More informationPrüfung Regelungstechnik I (Control Systems I) Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam!
Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 29. 8. 2 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid
More informationMath 312 Lecture Notes Linear Two-dimensional Systems of Differential Equations
Math 2 Lecture Notes Linear Two-dimensional Systems of Differential Equations Warren Weckesser Department of Mathematics Colgate University February 2005 In these notes, we consider the linear system of
More informationDifferential Equations 2280 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 2015 at 12:50pm
Differential Equations 228 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 25 at 2:5pm Instructions: This in-class exam is 5 minutes. No calculators, notes, tables or books. No answer check is
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 informationROOT LOCUS. Consider the system. Root locus presents the poles of the closed-loop system when the gain K changes from 0 to. H(s) H ( s) = ( s)
C1 ROOT LOCUS Consider the system R(s) E(s) C(s) + K G(s) - H(s) C(s) R(s) = K G(s) 1 + K G(s) H(s) Root locus presents the poles of the closed-loop system when the gain K changes from 0 to 1+ K G ( s)
More informationUnderstand the existence and uniqueness theorems and what they tell you about solutions to initial value problems.
Review Outline To review for the final, look over the following outline and look at problems from the book and on the old exam s and exam reviews to find problems about each of the following topics.. Basics
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 informationMassachusetts Institute of Technology
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.011: Introduction to Communication, Control and Signal Processing QUIZ 1, March 16, 2010 ANSWER BOOKLET
More informationIntroduction to Nonlinear Control Lecture # 4 Passivity
p. 1/6 Introduction to Nonlinear Control Lecture # 4 Passivity È p. 2/6 Memoryless Functions ¹ y È Ý Ù È È È È u (b) µ power inflow = uy Resistor is passive if uy 0 p. 3/6 y y y u u u (a) (b) (c) Passive
More informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Systems of Differential Equations Phase Plane Analysis Hanz Richter Mechanical Engineering Department Cleveland State University Systems of Nonlinear
More informationCDS 101: Lecture 4.1 Linear Systems
CDS : Lecture 4. Linear Systems Richard M. Murray 8 October 4 Goals: Describe linear system models: properties, eamples, and tools Characterize stability and performance of linear systems in terms of eigenvalues
More information1. The Transition Matrix (Hint: Recall that the solution to the linear equation ẋ = Ax + Bu is
ECE 55, Fall 2007 Problem Set #4 Solution The Transition Matrix (Hint: Recall that the solution to the linear equation ẋ Ax + Bu is x(t) e A(t ) x( ) + e A(t τ) Bu(τ)dτ () This formula is extremely important
More informationControls Problems for Qualifying Exam - Spring 2014
Controls Problems for Qualifying Exam - Spring 2014 Problem 1 Consider the system block diagram given in Figure 1. Find the overall transfer function T(s) = C(s)/R(s). Note that this transfer function
More informationLaplace Transforms Chapter 3
Laplace Transforms Important analytical method for solving linear ordinary differential equations. - Application to nonlinear ODEs? Must linearize first. Laplace transforms play a key role in important
More informationH(s) = s. a 2. H eq (z) = z z. G(s) a 2. G(s) A B. s 2 s(s + a) 2 s(s a) G(s) 1 a 1 a. } = (z s 1)( z. e ) ) (z. (z 1)(z e at )(z e at )
.7 Quiz Solutions Problem : a H(s) = s a a) Calculate the zero order hold equivalent H eq (z). H eq (z) = z z G(s) Z{ } s G(s) a Z{ } = Z{ s s(s a ) } G(s) A B Z{ } = Z{ + } s s(s + a) s(s a) G(s) a a
More informationEECS C128/ ME C134 Final Thu. May 14, pm. Closed book. One page, 2 sides of formula sheets. No calculators.
Name: SID: EECS C28/ ME C34 Final Thu. May 4, 25 5-8 pm Closed book. One page, 2 sides of formula sheets. No calculators. There are 8 problems worth points total. Problem Points Score 4 2 4 3 6 4 8 5 3
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 informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi System Stability - 26 March, 2014
Prof. Dr. Eleni Chatzi System Stability - 26 March, 24 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
More informationTable of Laplacetransform
Appendix Table of Laplacetransform pairs 1(t) f(s) oct), unit impulse at t = 0 a, a constant or step of magnitude a at t = 0 a s t, a ramp function e- at, an exponential function s + a sin wt, a sine fun
More informationLearn2Control Laboratory
Learn2Control Laboratory Version 3.2 Summer Term 2014 1 This Script is for use in the scope of the Process Control lab. It is in no way claimed to be in any scientific way complete or unique. Errors should
More informationExam 2 Study Guide: MATH 2080: Summer I 2016
Exam Study Guide: MATH 080: Summer I 016 Dr. Peterson June 7 016 First Order Problems Solve the following IVP s by inspection (i.e. guessing). Sketch a careful graph of each solution. (a) u u; u(0) 0.
More informationGATE EE Topic wise Questions SIGNALS & SYSTEMS
www.gatehelp.com GATE EE Topic wise Questions YEAR 010 ONE MARK Question. 1 For the system /( s + 1), the approximate time taken for a step response to reach 98% of the final value is (A) 1 s (B) s (C)
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 informationEigenspaces and Diagonalizable Transformations
Chapter 2 Eigenspaces and Diagonalizable Transformations As we explored how heat states evolve under the action of a diffusion transformation E, we found that some heat states will only change in amplitude.
More informatione st f (t) dt = e st tf(t) dt = L {t f(t)} s
Additional operational properties How to find the Laplace transform of a function f (t) that is multiplied by a monomial t n, the transform of a special type of integral, and the transform of a periodic
More information= A(λ, t)x. In this chapter we will focus on the case that the matrix A does not depend on time (so that the ODE is autonomous):
Chapter 2 Linear autonomous ODEs 2 Linearity Linear ODEs form an important class of ODEs They are characterized by the fact that the vector field f : R m R p R R m is linear at constant value of the parameters
More informationBIBO STABILITY AND ASYMPTOTIC STABILITY
BIBO STABILITY AND ASYMPTOTIC STABILITY FRANCESCO NORI Abstract. In this report with discuss the concepts of bounded-input boundedoutput stability (BIBO) and of Lyapunov stability. Examples are given to
More information