ECEEN 5448 Fall 2011 Homework #4 Solutions
|
|
- Roxanne Lyons
- 6 years ago
- Views:
Transcription
1 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 input being a force on the mass and output being the position. The states are the mass position and velocity. (a) Demonstrate the transfer function is /s 2. This shows the realization is controllable and observable (why)? si A = [ s s (si A) b = s 2 [ s s c(si A) b = s 2 [ [ s (si A) = [ s s 2 s [ (si A) b = [ s 2 s = s 2 This shows controllability and observability because no cancellations (b) Find the state-feedback gain k that makes det (si A + bk)) = s 2 +5s +6 ([ ) s det (si A + bk) =det = s 2 + k k s + k2 2 s + k So answer is k = [6 5. (c) Find the closed-loop transfer function H yr (s) with the control u = kx + r operating and k being the gain vector you found in b. Zeros don t move, new denominator is det (si A + bk) (d) With the initial condition H yr (s) =c (si A + bk) b = x() = s 2 +5s +6 = (s + 2)(s +3) Find the response y(t) if r and sketch it. (No computer plots, please; just sketch the thing. You can use MATLAB to plot it st if you want and then sketch from the plot). [ ẋ =(A bk)x; y = cx Y (s) =c(si A + bk) x = And so My sketch shown in Figure y(t) =e 2t e 3t (s + 2)(s +3) = s +2 s +3
2 2 Figure : Sketch of y(t) when u = kx and x() = [ T (e) Now find the observer gain l that makes det (si A + lc)) = s 2 +25s + 5. Dual problem to finding k. Answer is [ 25 l = 5 (f) With the control u = kˆx + r operating where ˆx is, of course, the observer state, find the response y(t) ifr and [ x() = ˆx() = No observer mis-match to start here, so x(t) for all time and response will be same as in d. y(t) =e 2t e 3t (g) Now with u = kˆx + r operating and r find the response y(t) if [ [ x() = and ˆx() = Here x() = x() ˆx() = [ Closed loop state space equations, for r are [ [ [ d x(t) A bk bk x(t) = ; y = [ c [ x(t) dt x(t) A lc x(t) x(t) and so Y (s) = [ c ([ si A + bk bk si A + lc We recall [ A B C and so [ A A = BC C ) [ x() x() Y (s) =c (si A + bk) x() + c (si A + bk) bk (si A + lc) x() }{{} dynamics caused by x() and so we see, explicitly the extra dyanamics in the response when estimated state-feedback is used and the observer error is not initially zero. Y (s) = = (s + 2)(s +3) + 5s + 3 (s + 2)(s +3) (s + )(s + 5) = s 2 +3s + 28 (s + 2)(s + 3)(s + )(s + 5) 225 4(s +2) 5 2(s +3) (s + ) 4 95(s + 5)
3 and so y(t) = e 2t 5 2 e 3t e t 4 95 e 5t The two responses, for u = kx and u = kˆx are plotted with MATLAB in Figure 2. The response that peaks higher is the one with u = kˆx Figure 2: Response to x() = [ T when u = kx and when u = kˆx and ˆx() =
4 4 2. Given the -dimensional system ẋ = x+u with x() = x, you will find the control u(t) which minimizes u2 (t) dt subject to the constraint that lim t x(t) = using two different approaches. First you will do it via LQR theory and and then you will do it by the theory of pseudo-inverses. (a) Setup the LQR problem of minimizing J ɛ (x,u)= ɛx 2 (t) +u 2 (t) dt where ɛ>. Formulate the associated ARE and find its unique, positive-definite solution, P. The ARE is 2p p 2 + ɛ = p = 2 ± 4+4ɛ 2 and picking the positive definite solution shows p =+ +ɛ (b) From (2a) above, compute u ɛ (t) the optimal state feedback minimizing J ɛ. Verify that A bk has LHP eigenvalues (so, for this -dimensional example, verify that A bk has negative real part.) The optimal feedback gain is b T P/ρ =+ +ɛ so here the optimal control is u ɛ (t) = ( + +ɛ)x(t) This makes A bk = ( + +ɛ) = +ɛ which is clearly in the LHP. (c) Now let ɛ in order to find the control, u (t), which minimizes u 2 (t) dt subject to lim t x(t) = Taking ɛ to zero gives u (t) = 2x(t) This control flips the unstable pole across the jω-axis which is a solution we expect. (d) Show that a control u L 2 [,T yields x(t )=iffl[u = e T x where L[u = e A(T τ) bu(τ) dτ is the usual input to final-state map studied in class and the final time T>is given and fixed. The solution to ẋ = Ax + bu, x() = x is so x(t )=means x(t) =e At x + t e A(t τ) bu(τ) dτ =e AT x + e A(T τ) bu(τ) dτ which means L[u = e AT x Here A =, of course, so the condition is L[u = e T x. (e) Let L be the adjoint of L. Find an expression for L [x. x, L[u = L [x,u x T e A(T τ) bu(τ) dτ = x T e A(T τ) b = (L [x(τ)) T L [x =b T e AT (T τ) x For our system, A =, and b =so L [x(t) =e T t x (L [x(τ)) T u(τ) dτ
5 5 (f) Let L denote the pseudo-inverse of L. Find L [x. Plugging in and so LL [x = e A(T τ) bb T e AT (T τ) xdτ = { [ L [x(t)= b T e AT (T t) e Aτ bb T e AT τ xdτ } e Aτ bb T e AT τ dτ x In our case, for A =and b =this is { [ } { } L [x(t) = e (T t) e 2τ (T t) 2e dτ x = e 2T x () (g) Using your psuedo-inverse, explicitly find the the control, u T (t) which minimizes u2 (t) dt subject to x(t )= A control takes use to x(t )=if, and only if L[u = e T x. The minimum norm solution to this system of linear equations is u T (t) =L [ e T x. Using our expression for L this becomes { } 2e (T t) u T (t) = e 2T e T x (2) (h) Verify that as T u T (t) u (t) from part (2c). It is NOT enough to simply take T in (2) and show you get 2e t x!! If you did only this, you get zero for this problem and a wrist slap for being sloppy. { } Instead apply u T (t) = 2e (T t) e T x e 2T to the system ẋ = x + u, x() = x and we get and so while 2x(t) is x(t) = e t x 2e2T e 2T x [ ( 2e = e t 2T + e 2T t e (t τ) (e τ ) dτ ) e t ( e 2t 2 u T (t) = 2e2T ( e t ) e 2T x ) x [ 2e t e 2T (e 2T ) + ( 2e 2T ) e t ( e 2t ) x ( e t ) ( e x [e 2t 2T = 2e2T e 2T e 2T = 2e2T ( e t ) [ e 2T e 2t e 2T x e 2T [ e 2T e 2t = u T (t) e 2T ) +( e 2t ) and NOW we can see that as T the bracketed term in (3) goes to for any finite t and thus lim T u T (t) = lim T 2x(t) where x(t) is solution to ẋ(t) =x(t)+u T (t), x() = x (3)
6 6 [ 2 3. Suppose P = 2 (a) If x T x how big can x T Px be? Why? Eigenvalues of P are, 3. Minimax theory says x T Px λ max (P )x T x and so x T Px 3 when x T x (b) If x T x = how small can x T Px be? Why? Similar. Now λ min (P ) comes into play. Minimax theory says x T Px λ min x T x and so x T Px when x T x =because here λ min (P )=. Thus smallest it can be is.
7 7 4. Proof or counterexample. You may assume all PD matrices are also symmetric. (a) If P is n n and PD and P 2 is n n and PD, then P P 2 + P 2 P is PD. This seems fishy since x T (P P 2 + P 2 P )x = x T P P 2 x + x T P 2 P x =2x T P P 2 x and we saw last HW that product doesn t preserve PD. In fact P = [ and P 2 = [ 2 2 is one counter-example. (b) If U and V are orthogonal matrices of size n n, then so are UV and VU. This is true. Proofs are easy. (UV) T UV = V T U T UV = V T V = I and UV(UV) T = UVV T U T = UU T = I. This shows UV is orthogonal. The proof about VU is similarly easy. (c) For any two n n matrices A and B, (si A) (si B) =(si A) (A B)(sI B). (If true, what does this say about parallel and cascade connections of transfer functions?) This is true. Hence (si A) (si B) = (si A) [(si B) (si A) (si B) = (si A) (A B)(sI B) c(si A) b c(si B) b = c(si A) (A B)(sI B) b and so the parallel connection of two transfer functions with the same c and b vectors can be represented as a cascade connection of a multi-output transfer function with the B and b dynamics and a multi-input transfer function with the c and A.
8 5. Consider the mass-spring damper system shown in Figure 3 of a mass M sliding on a fritionless surface, connected by a spring and a damper to an immovable wall. We can apply an external force, u to the mass as shown. With output y = q and using the natural, physical state K B M q u Figure 3: Two Masses Connected by a Linear Spring x = [ q q a state-space realization for the system is [ A = K/M B/M [ ; b = /M ; c = [ Let us suppose M = 6, K = 2 and B = 3 for the rest of this problem. (a) Find the eigenvalues of A and show that this system is asymptotically stable. For the parameters given [ A = ; det(si A) =s 2 +3/6s +3/4 LHP roots by inspection 3/4 3/6 (b) Suppose V (x) =x T x = x 2. If u = (input is turned off) and x(t) is the (undriven) system behavior starting from x() = x, show that dt d (V (x(t))) t= can be positive for some x. Clearly x T x = x T Px for P = I. ForV (x) =x T Px we know V (x) =x T (PA+A T P )x by a calculation done about eleven-bajillion times in class now. So here Now, A+A T = [ /4 /4 3/8 d dt (V (x(t))) t= = x T (A + A T )x ; det ( si (A + A T ) ) = s 2 +3/8 /6 λ i ( A + A T ) =/8, /2 This means (A+A T ) is not PD so we can find a non-zero vector x such that x (A + A T )x >. A perfectly good one is the eigenvector, v, corresponding to λ =/8 which is [ 8 v = So, answer is d dt (V (x(t))) t= > for x = 4 [ 8 4
9 9 (c) Find an initial condition, x, and a time T> so that x(t ) > x() (with the input still set to zero, of course). How can this be if the system is asymptotically stable? Explain, both mathematically and (for this physical system) physically. Let s start from the x found above in part 5b. Since the derivative of x T x is > at that point, if we make T small enough, we will have it. So you could numerically just try smaller and smaller T until you get one that works. But I ll go the long way to show you a more complete approach. OK. The solution, as we know, is x(t) =e At x and e At = L {( si A) )} = and so with x = [ 8 4 [ /4 /4 3/8 we get { x(t) = L 6s 2 +3s +2 = e 3t/32 6 and hence, starting from x = { [ = L 6s s 2 +3s s [ 28s s 96 } [ sin [Θ(t) + 28 cos [Θ(t) 64 cos [Θ(t)) sin [Θ(t) [ 8 4 we find where Θ(t) = 759t 32 } x T (t)x(t) = x(t) 2 = e 3t/6 ( 6472 sin 2 [Θ(t) + 8 ) 759 cos [Θ(t) sin [Θ(t) cos 2 [Θ(t) t where Θ(t) = 32 which I have plotted below in Figure 4. You can see that any T <. works. PLEASE NOTE! x(t) =e At x and x T (t)x(t) =x T eat t e At x but this is NOT x T e(at +A)t x! Why? Because A and A T do not generally commute. (d) Show that if P = [ then for V (x) =x T Px it is true that dt d (V (x(t))) x 2. Here PA+ A T P = [ and so V (x) = x 2 x 2 as was to be shown. (e) (MATLAB useful). Find a change of coordinate matrix T so that if x = T x, Ā = T TAT, c = ct, and b = T b as usual, then x(t ) < x() for all T> and for all x()
10 85 Norm of state starting from x = [8,4 T 8 75 x T (t)x(t) Time (sec) Figure 4: x(t) 2 for Mass-Spring Damper System when M = 6, K = 2, B =3,x = [8 4 T and ẋ = Ax Can you interpret the new states, x physically? The change of coordinate is [ T = P / = I cannot interpret these new states physically. They are non-intuitive combinations of the position and velocity!
Module 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 informationECEEN 5448 Fall 2011 Homework #5 Solutions
ECEEN 5448 Fall 211 Homework #5 Solutions Professor David G. Meyer December 8, 211 1. Consider the 1-dimensional time-varying linear system ẋ t (u x) (a) Find the state-transition matrix, Φ(t, τ). Here
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 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 information5. Observer-based Controller Design
EE635 - Control System Theory 5. Observer-based Controller Design Jitkomut Songsiri state feedback pole-placement design regulation and tracking state observer feedback observer design LQR and LQG 5-1
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 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 informationTopic # Feedback Control. State-Space Systems Closed-loop control using estimators and regulators. Dynamics output feedback
Topic #17 16.31 Feedback Control State-Space Systems Closed-loop control using estimators and regulators. Dynamics output feedback Back to reality Copyright 21 by Jonathan How. All Rights reserved 1 Fall
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 informationCourse Outline. Higher Order Poles: Example. Higher Order Poles. Amme 3500 : System Dynamics & Control. State Space Design. 1 G(s) = s(s + 2)(s +10)
Amme 35 : System Dynamics Control State Space Design Course Outline Week Date Content Assignment Notes 1 1 Mar Introduction 2 8 Mar Frequency Domain Modelling 3 15 Mar Transient Performance and the s-plane
More informationControl Systems Design
ELEC4410 Control Systems Design Lecture 18: State Feedback Tracking and State Estimation Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 18:
More informationIntro. Computer Control Systems: F9
Intro. Computer Control Systems: F9 State-feedback control and observers Dave Zachariah Dept. Information Technology, Div. Systems and Control 1 / 21 dave.zachariah@it.uu.se F8: Quiz! 2 / 21 dave.zachariah@it.uu.se
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 information6 OUTPUT FEEDBACK DESIGN
6 OUTPUT FEEDBACK DESIGN When the whole sate vector is not available for feedback, i.e, we can measure only y = Cx. 6.1 Review of observer design Recall from the first class in linear systems that a simple
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 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 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 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 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 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 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 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 information1 Steady State Error (30 pts)
Professor Fearing EECS C28/ME C34 Problem Set Fall 2 Steady State Error (3 pts) Given the following continuous time (CT) system ] ẋ = A x + B u = x + 2 7 ] u(t), y = ] x () a) Given error e(t) = r(t) y(t)
More informationEE C128 / ME C134 Fall 2014 HW 9 Solutions. HW 9 Solutions. 10(s + 3) s(s + 2)(s + 5) G(s) =
1. Pole Placement Given the following open-loop plant, HW 9 Solutions G(s) = 1(s + 3) s(s + 2)(s + 5) design the state-variable feedback controller u = Kx + r, where K = [k 1 k 2 k 3 ] is the feedback
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 informationCDS 101: Lecture 5-1 Reachability and State Space Feedback
CDS 11: Lecture 5-1 Reachability and State Space Feedback Richard M. Murray 23 October 26 Goals: Define reachability of a control system Give tests for reachability of linear systems and apply to examples
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 informationEECS C128/ ME C134 Final Wed. Dec. 14, am. Closed book. One page, 2 sides of formula sheets. No calculators.
Name: SID: EECS C128/ ME C134 Final Wed. Dec. 14, 211 81-11 am Closed book. One page, 2 sides of formula sheets. No calculators. There are 8 problems worth 1 points total. Problem Points Score 1 16 2 12
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 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 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 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 informationEECS C128/ ME C134 Final Wed. Dec. 15, am. Closed book. Two pages of formula sheets. No calculators.
Name: SID: EECS C28/ ME C34 Final Wed. Dec. 5, 2 8- am Closed book. Two pages of formula sheets. No calculators. There are 8 problems worth points total. Problem Points Score 2 2 6 3 4 4 5 6 6 7 8 2 Total
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 informationState Feedback and State Estimators Linear System Theory and Design, Chapter 8.
1 Linear System Theory and Design, http://zitompul.wordpress.com 2 0 1 4 State Estimator In previous section, we have discussed the state feedback, based on the assumption that all state variables are
More informationModel reduction for linear systems by balancing
Model reduction for linear systems by balancing Bart Besselink Jan C. Willems Center for Systems and Control Johann Bernoulli Institute for Mathematics and Computer Science University of Groningen, Groningen,
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 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 informationQuadratic Stability of Dynamical Systems. Raktim Bhattacharya Aerospace Engineering, Texas A&M University
.. Quadratic Stability of Dynamical Systems Raktim Bhattacharya Aerospace Engineering, Texas A&M University Quadratic Lyapunov Functions Quadratic Stability Dynamical system is quadratically stable if
More informationObservability and state estimation
EE263 Autumn 2015 S Boyd and S Lall Observability and state estimation state estimation discrete-time observability observability controllability duality observers for noiseless case continuous-time observability
More information06/12/ rws/jMc- modif SuFY10 (MPF) - Textbook Section IX 1
IV. Continuous-Time Signals & LTI Systems [p. 3] Analog signal definition [p. 4] Periodic signal [p. 5] One-sided signal [p. 6] Finite length signal [p. 7] Impulse function [p. 9] Sampling property [p.11]
More informationState Feedback and State Estimators Linear System Theory and Design, Chapter 8.
1 Linear System Theory and Design, http://zitompul.wordpress.com 2 0 1 4 2 Homework 7: State Estimators (a) For the same system as discussed in previous slides, design another closed-loop state estimator,
More informationLecture 19 Observability and state estimation
EE263 Autumn 2007-08 Stephen Boyd Lecture 19 Observability and state estimation state estimation discrete-time observability observability controllability duality observers for noiseless case continuous-time
More informationState Regulator. Advanced Control. design of controllers using pole placement and LQ design rules
Advanced Control State Regulator Scope design of controllers using pole placement and LQ design rules Keywords pole placement, optimal control, LQ regulator, weighting matrixes Prerequisites Contact state
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 informationRobust Multivariable Control
Lecture 2 Anders Helmersson anders.helmersson@liu.se ISY/Reglerteknik Linköpings universitet Today s topics Today s topics Norms Today s topics Norms Representation of dynamic systems Today s topics Norms
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 informationHomework Solution # 3
ECSE 644 Optimal Control Feb, 4 Due: Feb 17, 4 (Tuesday) Homework Solution # 3 1 (5%) Consider the discrete nonlinear control system in Homework # For the optimal control and trajectory that you have found
More informationTopic # Feedback Control Systems
Topic #17 16.31 Feedback Control Systems Deterministic LQR Optimal control and the Riccati equation Weight Selection Fall 2007 16.31 17 1 Linear Quadratic Regulator (LQR) Have seen the solutions to the
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 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 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 informationTRACKING AND DISTURBANCE REJECTION
TRACKING AND DISTURBANCE REJECTION Sadegh Bolouki Lecture slides for ECE 515 University of Illinois, Urbana-Champaign Fall 2016 S. Bolouki (UIUC) 1 / 13 General objective: The output to track a reference
More information16.30/31, Fall 2010 Recitation # 13
16.30/31, Fall 2010 Recitation # 13 Brandon Luders December 6, 2010 In this recitation, we tie the ideas of Lyapunov stability analysis (LSA) back to previous ways we have demonstrated stability - but
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 informationRobust 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 informationTopic # Feedback Control
Topic #5 6.3 Feedback Control State-Space Systems Full-state Feedback Control How do we change the poles of the state-space system? Or,evenifwecanchangethepolelocations. Where do we put the poles? Linear
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 informationLinear-quadratic control problem with a linear term on semiinfinite interval: theory and applications
Linear-quadratic control problem with a linear term on semiinfinite interval: theory and applications L. Faybusovich T. Mouktonglang Department of Mathematics, University of Notre Dame, Notre Dame, IN
More informationSYSTEMTEORI - ÖVNING 5: FEEDBACK, POLE ASSIGNMENT AND OBSERVER
SYSTEMTEORI - ÖVNING 5: FEEDBACK, POLE ASSIGNMENT AND OBSERVER Exercise 54 Consider the system: ẍ aẋ bx u where u is the input and x the output signal (a): Determine a state space realization (b): Is the
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 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 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 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 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 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 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 information= m(0) + 4e 2 ( 3e 2 ) 2e 2, 1 (2k + k 2 ) dt. m(0) = u + R 1 B T P x 2 R dt. u + R 1 B T P y 2 R dt +
ECE 553, Spring 8 Posted: May nd, 8 Problem Set #7 Solution Solutions: 1. The optimal controller is still the one given in the solution to the Problem 6 in Homework #5: u (x, t) = p(t)x k(t), t. The minimum
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 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 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 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 informationOPTIMAL CONTROL. Sadegh Bolouki. Lecture slides for ECE 515. University of Illinois, Urbana-Champaign. Fall S. Bolouki (UIUC) 1 / 28
OPTIMAL CONTROL Sadegh Bolouki Lecture slides for ECE 515 University of Illinois, Urbana-Champaign Fall 2016 S. Bolouki (UIUC) 1 / 28 (Example from Optimal Control Theory, Kirk) Objective: To get from
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 informationEE363 homework 2 solutions
EE363 Prof. S. Boyd EE363 homework 2 solutions. Derivative of matrix inverse. Suppose that X : R R n n, and that X(t is invertible. Show that ( d d dt X(t = X(t dt X(t X(t. Hint: differentiate X(tX(t =
More informationMIT Final Exam Solutions, Spring 2017
MIT 8.6 Final Exam Solutions, Spring 7 Problem : For some real matrix A, the following vectors form a basis for its column space and null space: C(A) = span,, N(A) = span,,. (a) What is the size m n of
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 information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 12: I/O Stability Readings: DDV, Chapters 15, 16 Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology March 14, 2011 E. Frazzoli
More informationControl Systems Design, SC4026. SC4026 Fall 2010, dr. A. Abate, DCSC, TU Delft
Control Systems Design, SC426 SC426 Fall 2, dr A Abate, DCSC, TU Delft Lecture 5 Controllable Canonical and Observable Canonical Forms Stabilization by State Feedback State Estimation, Observer Design
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 informationTopic # Feedback Control Systems
Topic #16 16.31 Feedback Control Systems State-Space Systems Closed-loop control using estimators and regulators. Dynamics output feedback Back to reality Fall 2007 16.31 16 1 Combined Estimators and Regulators
More informationSYSTEMTEORI - KALMAN FILTER VS LQ CONTROL
SYSTEMTEORI - KALMAN FILTER VS LQ CONTROL 1. Optimal regulator with noisy measurement Consider the following system: ẋ = Ax + Bu + w, x(0) = x 0 where w(t) is white noise with Ew(t) = 0, and x 0 is a stochastic
More informationProblem Set 4 Solution 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Problem Set 4 Solution Problem 4. For the SISO feedback
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 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 informationProblem 1 Cost of an Infinite Horizon LQR
THE UNIVERSITY OF TEXAS AT SAN ANTONIO EE 5243 INTRODUCTION TO CYBER-PHYSICAL SYSTEMS H O M E W O R K # 5 Ahmad F. Taha October 12, 215 Homework Instructions: 1. Type your solutions in the LATEX homework
More informationNPTEL Online Course: Control Engineering
NPTEL Online Course: Control Engineering Ramkrishna Pasumarthy Assignment-11 : s 1. Consider a system described by state space model [ ] [ 0 1 1 x + u 5 1 2] y = [ 1 2 ] x What is the transfer function
More informationEE363 homework 8 solutions
EE363 Prof. S. Boyd EE363 homework 8 solutions 1. Lyapunov condition for passivity. The system described by ẋ = f(x, u), y = g(x), x() =, with u(t), y(t) R m, is said to be passive if t u(τ) T y(τ) dτ
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 informationFinal Exam Solutions
EE55: Linear Systems Final Exam SIST, ShanghaiTech Final Exam Solutions Course: Linear Systems Teacher: Prof. Boris Houska Duration: 85min YOUR NAME: (type in English letters) I Introduction This exam
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 information1 (30 pts) Dominant Pole
EECS C8/ME C34 Fall Problem Set 9 Solutions (3 pts) Dominant Pole For the following transfer function: Y (s) U(s) = (s + )(s + ) a) Give state space description of the system in parallel form (ẋ = Ax +
More informationMassachusetts Institute of Technology. Department of Electrical Engineering and Computer Science : MULTIVARIABLE CONTROL SYSTEMS by A.
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Q-Parameterization 1 This lecture introduces the so-called
More informationMAE 143B - Homework 8 Solutions
MAE 43B - Homework 8 Solutions P6.4 b) With this system, the root locus simply starts at the pole and ends at the zero. Sketches by hand and matlab are in Figure. In matlab, use zpk to build the system
More informationTopic # /31 Feedback Control Systems
Topic #16 16.30/31 Feedback Control Systems Add reference inputs for the DOFB case Reading: FPE 7.8, 7.9 Fall 2010 16.30/31 16 2 Reference Input - II On page 15-6, compensator implemented with reference
More informationTopic # /31 Feedback Control Systems. Analysis of Nonlinear Systems Lyapunov Stability Analysis
Topic # 16.30/31 Feedback Control Systems Analysis of Nonlinear Systems Lyapunov Stability Analysis Fall 010 16.30/31 Lyapunov Stability Analysis Very general method to prove (or disprove) stability of
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 informationLecture 9. Introduction to Kalman Filtering. Linear Quadratic Gaussian Control (LQG) G. Hovland 2004
MER42 Advanced Control Lecture 9 Introduction to Kalman Filtering Linear Quadratic Gaussian Control (LQG) G. Hovland 24 Announcement No tutorials on hursday mornings 8-9am I will be present in all practical
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 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 information