1 Steady State Error (30 pts)
|
|
- Douglas Porter
- 6 years ago
- Views:
Transcription
1 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 ] u(t), y = ] x () a) Given error e(t) = r(t) y(t) where r(t) is a scalar, find lim t e(t) for input r(t) a unit step, and control law u(t) = r(t) y(t). Recall that ( ) e ss = I + CĀ B r ss Ā = A BKC B = BK when K is the controller gain. In this case, K =. Therefore: ] Ā = 2 7 = 2 7 ] = 22 7 ] Ā = B = = I + CĀ B = + ] 22 = 5 e ss = 6 r ss = b) Evaluate the steady state error lim t e(t) for input r(t) a unit step, with state feedback, that is, u = K x + r, where K is chosen so that the closed loop poles are at s i = 5, 5. First, find K. We want the eigenvalues of A B K to be ( 5, 5). Find the characteristic
2 equation, then make it equal (s + 5)(s + 5): ] A B K = k 2 7 k 2 ] = 2 k 7 k 2 ] (s) = s 2 k 7 k 2 s k + 2 = 25 = s(s k 2 ) + k + 2 = s 2 + (7 + k 2 )s + (k + 2) = (s + 5)(s + 5) = s 2 + s + 25 k = 3 k = k 2 = 3 Now find the steady state error. Note that the method we used in (a) had u = K(r y). This time, u = r K x. The difference is that r does not get multiplied by K. We need to re-derive the steady state error relationship. ẋ = (A B K )x + B r Note there s no C, because we re doing full state feedback. Assume that ẋ = in steady state. = (A B K )x ss + B r ss x ss = (A B K ) B r ss e ss = r ss C x ss ) = (I + C (A B K ) B r ss Now we can evaluate e ss. ] A B K = = ] = 25 2
3 (actually, we could have jumped right to this, since A, B are in controllable canonical form, so A B K will have the coefficients of the characteristic equation in the bottom row.) (A B K ) = 25 I + C (A B K ) B = + ] 25 = 2 5 e ss = 3 5 r ss = 3 5 ] For parts c) and d) below, we want to add an integral controller to the system ( ) to get zero steady state error, but still have state feedback to allow arbitrary pole placement. c) Let u = K 2 x+k e (r Cx)dτ. Write the augmented state and output equations for this system with integral control, for example, ż = A 2 z + B 2 r, where A 2 is 3 3. Find gains K 2, k e such that the system has closed-loop poles at s = 4, 5,, and evaluate the steady state error for a step input r(t). We are adding a new state variable to represent the integrator. To avoid confusing it with the original x state variables, we ll call it w. w z = Manipulate the definition of u to get dynamics for w. u = K 2 x + k e (r C x) dτ x x 2 u = K 2 x + k e w w = (r C x) dτ ẇ = (r C x) Notice that the new form, ż = A 2 z + B 2 r, doesn t have u. We need to make A 2 capture the inside loop of the full state feedback. Manipulate the original dynamics. ẋ = A x + B u = A x + B ( K 2 x + k e w) = (A B K 2 )x + B k e w 3
4 Now we can write the state and output equations. ż = A 2 z + B 2 r y = C 2 z ] ] C ż = w w + r y = C B k e A B K 2 x x w ] w ż = x + r y = x k e 2 k 7 k 2 x 2 x 2 To find the gains K 2, k e, find the characteristic equation for A 2, and make it equal (s + 4)(s + 5)(s + ): (s) = s s k e 2 k s 7 k 2 ] ] ] = s s 2 k s 7 k 2 s 7 k 2 k e + s k e 2 k = s(s(s k 2 ) k ) (k e ) = s 3 (7 + k 2 )s 2 (2 + k )s k e We can flip the sign of (s) for convenience, since we only care about its roots = s 3 + (7 + k 2 )s 2 + (2 + k )s + k e = (s + 4)(s + 5)(s + ) = s 3 + 9s 2 + s k 2 = 9 k 2 = k = k = 98 k e = 2 k e = 2 Now we can evaluate the steady state error. ż = A 2 z + B 2 r = A 2 z ss + B 2 r ss z ss = A 2 B 2r ss y ss = C 2 A 2 B 2r ss e ss = r ss y ss ( ) = I + C 2 A 2 B 2 r ss ] = = 4
5 d) Augment the system matrix A to include an integrator (keeping the phase variable form), and then choose u = K 3 (r ] T x). Write the augmented state and output equations for this system for example, ż = A 3 z + B 3 r, where A 3 is 3 3. Find gains K 3 such that the system has closed-loop poles at s = 4, 5,, and evaluate the steady state error for a step input r(t). Here is the signal flow diagram of the original system. Notice that we start in phase variable form. u x 2 x y 2 7 There is only one way to add an integrator to this system that keeps the phase variable form and that does not modify the dynamics of the original system: after x. u x 2 x x INT y 2 7 Notice that this new integrator doesn t actually affect anything (both gains leading away from x INT are zero). That s OK; when we do full state feedback, this new state will start affecting the rest of the system. We can write the pre-state-feedback system equations based on this diagram: z = x INT x x 2 ż = A 3z + B 3u = z + u 2 7 ] y = C 3 z = z Next, we have to bring the inside loop into A 3, like we did in part (c). To do this, we write 5
6 out the dynamics, then get rid of u by substitution. ż = A 3z + B 3u = A 3z + B 3K 3 (r ] T x) = (A 3 B 3K ] 3 )z + k k 2 k 3 r = z + r k 2 k 2 7 k 3 k ] y = z Now we have A 3, B 3, and C 3. Next, we need to find the K 3 values. This time it s easier, because A 3, B 3 are in phase variable form. We want the characteristic equation to equal (s + 4)(s + 5)(s + ) = s 3 + 9s 2 + s + 2, so the bottom row of A 3 must be ( 2,, 9). k = 2 k = 2 2 k 2 = k 2 = 98 7 k 3 = 9 k 3 = 2 Finally, we evaluate the steady state error. Because the system in (c) and (d) have the same form: ż = A i z + B i r: we can re-use the derivation of the steady state error from (c). ( ) e ss = I + C 3 A 3 B 3 r ss ] = = What s going on here? The open loop system (defined by A 3, B 3 ) clearly has a pole at the origin (one of the elements of the bottom row of A 3 is zero). Doesn t that mean it s a Type system, and therefore zero steady state error when driven by a unit step? There are a couple of ways of looking at this. As mentioned earlier, that the new integrator isn t actually doing anything in the open loop form. It s integrating x, but the result of the integration doesn t affect the other states or the output. Another way to look at it is to convert the open loop state space system into a transfer function. (See section 3.6 of Nise). When you do this, you find that the TF representation of the open loop system doesn t have a pole at the origin. 2 Linear Quadratic Regulator (35 pts) STOP STOP x3 safety zone x safety zone u2 car 2 u 2 car t = sec car 2 car t = sec 6
7 u Consider a pair of cars in a platoon regulated to a stopping position, as illustrated above. The dynamics of car are (x 2 = ẍ = u and car 2 has a plant model x 4 = ẍ 3 =.5u 2 where u and u 2 are the car s thrust due to engine and braking. (The offset of x from car s rear bumper prevents the cars from colliding when x 3 = x.) The outputs of the system are y = x and y 2 = x 3 x. Note that if y 2 > then car 2 has intruded on the safety zone of car. Initial conditions are car at - m, 2 m/s, and car 2 at -2 m, 25 m/s. a) Write the system equations in state space form. ẋ ẋ 2 ẋ 3 ẋ 4 y y 2 x = x 2 x 3 + x 4.5 ] ] x x = 2 2 x() = 2 25 b) Use the LQR method (from Matlab function lqr(sys,q,r), with Q=diag(2.5,,2.5,]) and R=diag(,]) to find an optimal K for the state feedback control u = K b x. Plot x(t) and u(t) for the given initial condition (Matlab initial) and state feedback with gain K b. Response to Initial Conditions x 3 x 4 ] u u 2 2 To: Out() Control effort 8 Amplitude To: Out(2) u u Time (seconds) K b = t ] Notice that the safety zone is violated by about 5cm at t = 7s. c) In part b) both cars get within m of the stop sign at a speed less than m/s at seconds. However, car 2 lags more than 4 meters behind car during the stopping trajectory. Find a new 7
8 u cost function Q = diag(q q 2 q 3 q 4 ]) which keeps car 2 within 2 m of car for the whole trajectory, while maintaining y 2 < to prevent a collision, and distance from stop sign less than m at sec with velocity less than m/s. (Suggestion: let q = 2.5 and q 2 = and search for new q 3 and q 4.) Plot x(t) and u(t) for the given initial condition (Matlab initial) and state feedback with new gain K c. There are many solutions that satisfy the constraints. The gray region in the following graph shows q 3, q 4 values that work Safety 45 q Final Distance q 3 Taking Q=diag(2.5 4]) gives the following response: Response to Initial Conditions 4 Control effort 2 2 Amplitude To: Out() To: Out(2) u u Time (seconds) K c = t ] d) Find the solution to the Riccati equation P using Matlab function care(a,b,q,r) and estimate the cost J = (x T P x)() for each of b) and c). 8
9 P b = x T ()P b x() = P c = x T ()P c x() = 2683 Answers for P c will vary depending on the choice of Q. Notice that the top left corner of P c is the same, since q and q 2 were not changed. e) Briefly compare the tradeoffs between control effort and time response between the two cases. Notice that the time response, control effort, and K entries associated with the first car are identical in (b) and (c). This will be true as long as the suggestion to keep q and q 2 the same is followed. Car 2 s response, control effort, and K entries are different between the two cases. The control effort used in (c) is bigger in magnitude than in (b). Correspondingly, the time response of y 2 (the distance between the two cars) is much different as well: in (c), it seems to be decaying to zero much faster than in (b). y 2 in (b) shows evidence of several modes with similar time constants at work (the back-and-forth motion). In (c), y 2 seems to be dominated by a single, faster mode. In general, we get better performance (faster convergence to zero) in (c) compared to (b), at the cost of greater control effort. 3 Discrete Time Control (35 pts) Given the following continuous time (CT) system ẋ = Ax + Bu = x u(t), y = 2 ] x the corresponding discrete time (DT) system is xn + ] = Gxn] + Hun], yn] = Cxn] which can be found using the Matlab function c2d(sys,t, zoh ). a) With initial condition x = ], plot the ZIR using Matlab function initial() for the CT system and the DT system (with T =.4 sec). >> A = ; ; ]; >> B = ;;]; >> C = 2 ]; >> D = ; >> ctsys = ss(a,b,c,d); 9
10 >> dtsys = c2d(ctsys,.4, zoh ); >> initial(ctsys,;;],:.:4]); >> initial(dtsys,;;],:.4:4]); Continuous time response: Zero-Order Hold Discrete time response: b) For output feedback control u = k(r y), sketch the root locus for the equivalent transfer function for the continuous time (CT) system. >> NUM,DEN] = ss2tf(a,b,c,d) NUM =
11 DEN = Root locus for TF: s + 2 s 3 + 5s 2 + 2s + 8 = s + 2 (s + )(s j)(s + 2 2j) c) Determine the closed loop pole locations for the CT system for k = 5 and plot the closed-loop step response using Matlab. >> z,p,k] = zpkdata(feedback(5*ss(a,b,c,d),)); >> p{} ans = i i >> step(feedback(5*ss(a,b,c,d),),:.:4]); Continuous step response for system with feedback:
12 d) The closed loop DT system has state equation xn + ] = (G khc)xn] + khrn], yn] = Cxn] (which can be found using the Matlab feedback function). Using Matlab, determine the closed loop pole locations for the DT system for k = 5 and sampling period T =.4 sec and plot the step response. >> z,p,k] = zpkdata(feedback(5*dtsys,)); >> p{} ans = i i >> step(feedback(5*dtsys,),:.4:4]); Discrete step response for system with feedback: 2
13 e) Use Matlab (iteratively if necessary) to find a sampling period T which gives a closed-loop step response that is reasonably close to the CT closed-loop step response. Determine closed-loop pole locations, and plot the DT step response. >> step(feedback(5*c2d(ctsys,.4, zoh ),)) >> step(feedback(5*c2d(ctsys,.3, zoh ),)) >> step(feedback(5*c2d(ctsys,.2, zoh ),)) >> step(feedback(5*c2d(ctsys,., zoh ),),:.:4]) >> z,p,k] = zpkdata(feedback(5*c2d(ctsys,., zoh ),)); >> p{} ans = i i With sampling period of.s, there is at least 8 samples for the first hump and the tail settles quickly to the equilibrium value without jittering around. Discrete step response for system with feedback (Sampling Period changed to T =.s) 3
14 f) Briefly explain why the CT and DT ZIR responses from a) above are reasonably close, but the closed loop responses from c) and d) (with T =.4 sec) do not agree at all. (Hint, consider e AT.) We know there will be some error in the open-loop responses due to discritization error. Intuitively, controlling the system based on an erroneous output would cause even more error in the responses. Lets look at a quantitative analysis of the differences: For the discrete system: G = e A.4 = and, in feedback, G 5HC = Suppose we discretized the closed loop continous system instead so that the sampler is outside the controller: G cl = e (A 5BC).4 = Note how these two matrices are similar, but there are a few entries that are off by a few tenths. In the first case, the controller is acting on observations containing discretization error, wheras in the second case, the controller is acting continously, and then the system is discretized. The difference in the matrices show that a controller can, in effect, amplify discretization error when modelling a continuous system. 4
1 (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 informationProfessor Fearing EE C128 / ME C134 Problem Set 10 Solution Fall 2010 Jansen Sheng and Wenjie Chen, UC Berkeley
Professor Fearing EE C28 / ME C34 Problem Set Solution Fall 2 Jansen Sheng and Wenjie Chen, UC Berkeley. (5 pts) Final Value Given the following continuous time (CT) system ẋ = Ax+Bu = 5 9 7 x+ u(t), y
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 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 informationMODERN CONTROL DESIGN
CHAPTER 8 MODERN CONTROL DESIGN The classical design techniques of Chapters 6 and 7 are based on the root-locus and frequency response that utilize only the plant output for feedback with a dynamic controller
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 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 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 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 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 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 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 informationProfessor Fearing EE C128 / ME C134 Problem Set 7 Solution Fall 2010 Jansen Sheng and Wenjie Chen, UC Berkeley
Professor Fearing EE C8 / ME C34 Problem Set 7 Solution Fall Jansen Sheng and Wenjie Chen, UC Berkeley. 35 pts Lag compensation. For open loop plant Gs ss+5s+8 a Find compensator gain Ds k such that the
More informationECEEN 5448 Fall 2011 Homework #4 Solutions
ECEEN 5448 Fall 2 Homework #4 Solutions Professor David G. Meyer Novemeber 29, 2. The state-space realization is A = [ [ ; b = ; c = [ which describes, of course, a free mass (in normalized units) with
More 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 informationEE C128 / ME C134 Final Exam Fall 2014
EE C128 / ME C134 Final Exam Fall 2014 December 19, 2014 Your PRINTED FULL NAME Your STUDENT ID NUMBER Number of additional sheets 1. No computers, no tablets, no connected device (phone etc.) 2. Pocket
More informationME 132, Fall 2017, UC Berkeley, A. Packard 317. G 1 (s) = 3 s + 6, G 2(s) = s + 2
ME 132, Fall 2017, UC Berkeley, A. Packard 317 Be sure to check that all of your matrix manipulations have the correct dimensions, and that the concatenations have compatible dimensions (horizontal concatenations
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 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 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 informationX 2 3. Derive state transition matrix and its properties [10M] 4. (a) Derive a state space representation of the following system [5M] 1
QUESTION BANK 6 SIDDHARTH GROUP OF INSTITUTIONS :: PUTTUR Siddharth Nagar, Narayanavanam Road 5758 QUESTION BANK (DESCRIPTIVE) Subject with Code :SYSTEM THEORY(6EE75) Year &Sem: I-M.Tech& I-Sem UNIT-I
More informationEE221A Linear System Theory Final Exam
EE221A Linear System Theory Final Exam Professor C. Tomlin Department of Electrical Engineering and Computer Sciences, UC Berkeley Fall 2016 12/16/16, 8-11am Your answers must be supported by analysis,
More informationUniversity of Toronto Department of Electrical and Computer Engineering ECE410F Control Systems Problem Set #3 Solutions = Q o = CA.
University of Toronto Department of Electrical and Computer Engineering ECE41F Control Systems Problem Set #3 Solutions 1. The observability matrix is Q o C CA 5 6 3 34. Since det(q o ), the matrix is
More 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 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 informationSuppose that we have a specific single stage dynamic system governed by the following equation:
Dynamic Optimisation Discrete Dynamic Systems A single stage example Suppose that we have a specific single stage dynamic system governed by the following equation: x 1 = ax 0 + bu 0, x 0 = x i (1) where
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 informationDepartment of Electronics and Instrumentation Engineering M. E- CONTROL AND INSTRUMENTATION ENGINEERING CL7101 CONTROL SYSTEM DESIGN Unit I- BASICS AND ROOT-LOCUS DESIGN PART-A (2 marks) 1. What are the
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 informationLast week: analysis of pinion-rack w velocity feedback
Last week: analysis of pinion-rack w velocity feedback Calculation of the steady state error Transfer function: V (s) V ref (s) = 0.362K s +2+0.362K Step input: V ref (s) = s Output: V (s) = s 0.362K s
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 informationExtensions and applications of LQ
Extensions and applications of LQ 1 Discrete time systems 2 Assigning closed loop pole location 3 Frequency shaping LQ Regulator for Discrete Time Systems Consider the discrete time system: x(k + 1) =
More informationEE C128 / ME C134 Fall 2014 HW 6.2 Solutions. HW 6.2 Solutions
EE C28 / ME C34 Fall 24 HW 6.2 Solutions. PI Controller For the system G = K (s+)(s+3)(s+8) HW 6.2 Solutions in negative feedback operating at a damping ratio of., we are going to design a PI controller
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 informationLINEAR QUADRATIC GAUSSIAN
ECE553: Multivariable Control Systems II. LINEAR QUADRATIC GAUSSIAN.: Deriving LQG via separation principle We will now start to look at the design of controllers for systems Px.t/ D A.t/x.t/ C B u.t/u.t/
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 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 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 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 informationFull State Feedback for State Space Approach
Full State Feedback for State Space Approach State Space Equations Using Cramer s rule it can be shown that the characteristic equation of the system is : det[ si A] 0 Roots (for s) of the resulting polynomial
More informationProfessor Fearing EECS120/Problem Set 2 v 1.01 Fall 2016 Due at 4 pm, Fri. Sep. 9 in HW box under stairs (1st floor Cory) Reading: O&W Ch 1, Ch2.
Professor Fearing EECS120/Problem Set 2 v 1.01 Fall 20 Due at 4 pm, Fri. Sep. 9 in HW box under stairs (1st floor Cory) Reading: O&W Ch 1, Ch2. Note: Π(t) = u(t + 1) u(t 1 ), and r(t) = tu(t) where u(t)
More informationMS-E2133 Systems Analysis Laboratory II Assignment 2 Control of thermal power plant
MS-E2133 Systems Analysis Laboratory II Assignment 2 Control of thermal power plant How to control the thermal power plant in order to ensure the stable operation of the plant? In the assignment Production
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 informationPitch Rate CAS Design Project
Pitch Rate CAS Design Project Washington University in St. Louis MAE 433 Control Systems Bob Rowe 4.4.7 Design Project Part 2 This is the second part of an ongoing project to design a control and stability
More informationSolution to HW State-feedback control of the motor with load (from text problem 1.3) (t) I a (t) V a (t) J L.
EE/ME 7: Advanced Linear Systems Solution Solution to HW Name On each item, grading is, or (full points, /2 points, points) EE/ME 7: Advanced Linear Systems Solution State-feedback control of the motor
More informationMassachusetts Institute of Technology Department of Mechanical Engineering Dynamics and Control II Design Project
Massachusetts Institute of Technology Department of Mechanical Engineering.4 Dynamics and Control II Design Project ACTIVE DAMPING OF TALL BUILDING VIBRATIONS: CONTINUED Franz Hover, 5 November 7 Review
More informationControl Systems. Design of State Feedback Control.
Control Systems Design of State Feedback Control chibum@seoultech.ac.kr Outline Design of State feedback control Dominant pole design Symmetric root locus (linear quadratic regulation) 2 Selection of closed-loop
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 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 informationInverted Pendulum: State-Space Methods for Controller Design
1 de 12 18/10/2015 22:45 Tips Effects TIPS ABOUT BASICS HARDWARE INDEX NEXT INTRODUCTION CRUISE CONTROL MOTOR SPEED MOTOR POSITION SYSTEM MODELING ANALYSIS Inverted Pendulum: State-Space Methods for Controller
More informationSUCCESSIVE POLE SHIFTING USING SAMPLED-DATA LQ REGULATORS. Sigeru Omatu
SUCCESSIVE POLE SHIFING USING SAMPLED-DAA LQ REGULAORS oru Fujinaka Sigeru Omatu Graduate School of Engineering, Osaka Prefecture University, 1-1 Gakuen-cho, Sakai, 599-8531 Japan Abstract: Design of sampled-data
More informationNote. Design via State Space
Note Design via State Space Reference: Norman S. Nise, Sections 3.5, 3.6, 7.8, 12.1, 12.2, and 12.8 of Control Systems Engineering, 7 th Edition, John Wiley & Sons, INC., 2014 Department of Mechanical
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 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 informationLinear-Quadratic Optimal Control: Full-State Feedback
Chapter 4 Linear-Quadratic Optimal Control: Full-State Feedback 1 Linear quadratic optimization is a basic method for designing controllers for linear (and often nonlinear) dynamical systems and is actually
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 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 informationCDS 110b: Lecture 2-1 Linear Quadratic Regulators
CDS 110b: Lecture 2-1 Linear Quadratic Regulators Richard M. Murray 11 January 2006 Goals: Derive the linear quadratic regulator and demonstrate its use Reading: Friedland, Chapter 9 (different derivation,
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 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 informationDr Ian R. Manchester Dr Ian R. Manchester AMME 3500 : Review
Week Date Content Notes 1 6 Mar Introduction 2 13 Mar Frequency Domain Modelling 3 20 Mar Transient Performance and the s-plane 4 27 Mar Block Diagrams Assign 1 Due 5 3 Apr Feedback System Characteristics
More informationDigital Control Engineering Analysis and Design
Digital Control Engineering Analysis and Design M. Sami Fadali Antonio Visioli AMSTERDAM BOSTON HEIDELBERG LONDON NEW YORK OXFORD PARIS SAN DIEGO SAN FRANCISCO SINGAPORE SYDNEY TOKYO Academic Press is
More informationDO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED.
EE 537 Homewors Friedland Text Updated: Wednesday November 8 Some homewor assignments refer to Friedland s text For full credit show all wor. Some problems require hand calculations. In those cases do
More informationSubject: Optimal Control Assignment-1 (Related to Lecture notes 1-10)
Subject: Optimal Control Assignment- (Related to Lecture notes -). Design a oil mug, shown in fig., to hold as much oil possible. The height and radius of the mug should not be more than 6cm. The mug must
More informationDigital Control: Part 2. ENGI 7825: Control Systems II Andrew Vardy
Digital Control: Part 2 ENGI 7825: Control Systems II Andrew Vardy Mapping the s-plane onto the z-plane We re almost ready to design a controller for a DT system, however we will have to consider where
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 informationME 132, Fall 2017, UC Berkeley, A. Packard 334 # 6 # 7 # 13 # 15 # 14
ME 132, Fall 2017, UC Berkeley, A. Packard 334 30.3 Fall 2017 Final # 1 # 2 # 3 # 4 # 5 # 6 # 7 # 8 NAME 20 15 20 15 15 18 15 20 # 9 # 10 # 11 # 12 # 13 # 14 # 15 # 16 18 12 12 15 12 20 18 15 Facts: 1.
More informationControl Systems I Lecture 10: System Specifications
Control Systems I Lecture 10: System Specifications Readings: Guzzella, Chapter 10 Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich November 24, 2017 E. Frazzoli (ETH) Lecture
More informationRobotics. Control Theory. Marc Toussaint U Stuttgart
Robotics Control Theory Topics in control theory, optimal control, HJB equation, infinite horizon case, Linear-Quadratic optimal control, Riccati equations (differential, algebraic, discrete-time), controllability,
More informationState space control for the Two degrees of freedom Helicopter
State space control for the Two degrees of freedom Helicopter AAE364L In this Lab we will use state space methods to design a controller to fly the two degrees of freedom helicopter. 1 The state space
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 information= 0 otherwise. Eu(n) = 0 and Eu(n)u(m) = δ n m
A-AE 567 Final Homework Spring 212 You will need Matlab and Simulink. You work must be neat and easy to read. Clearly, identify your answers in a box. You will loose points for poorly written work. You
More informationx(n + 1) = Ax(n) and y(n) = Cx(n) + 2v(n) and C = x(0) = ξ 1 ξ 2 Ex(0)x(0) = I
A-AE 567 Final Homework Spring 213 You will need Matlab and Simulink. You work must be neat and easy to read. Clearly, identify your answers in a box. You will loose points for poorly written work. You
More informationChapter 3. State Feedback - Pole Placement. Motivation
Chapter 3 State Feedback - Pole Placement Motivation Whereas classical control theory is based on output feedback, this course mainly deals with control system design by state feedback. This model-based
More informationLecture 10: Powers of Matrices, Difference Equations
Lecture 10: Powers of Matrices, Difference Equations Difference Equations A difference equation, also sometimes called a recurrence equation is an equation that defines a sequence recursively, i.e. each
More informationSAMPLE SOLUTION TO EXAM in MAS501 Control Systems 2 Autumn 2015
FACULTY OF ENGINEERING AND SCIENCE SAMPLE SOLUTION TO EXAM in MAS501 Control Systems 2 Autumn 2015 Lecturer: Michael Ruderman Problem 1: Frequency-domain analysis and control design (15 pt) Given is a
More informationOften, in this class, we will analyze a closed-loop feedback control system, and end up with an equation of the form
ME 32, Spring 25, UC Berkeley, A. Packard 55 7 Review of SLODEs Throughout this section, if y denotes a function (of time, say), then y [k or y (k) denotes the k th derivative of the function y, y [k =
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 information1 x(k +1)=(Φ LH) x(k) = T 1 x 2 (k) x1 (0) 1 T x 2(0) T x 1 (0) x 2 (0) x(1) = x(2) = x(3) =
567 This is often referred to as Þnite settling time or deadbeat design because the dynamics will settle in a Þnite number of sample periods. This estimator always drives the error to zero in time 2T or
More informationD(s) G(s) A control system design definition
R E Compensation D(s) U Plant G(s) Y Figure 7. A control system design definition x x x 2 x 2 U 2 s s 7 2 Y Figure 7.2 A block diagram representing Eq. (7.) in control form z U 2 s z Y 4 z 2 s z 2 3 Figure
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 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 informationExam. 135 minutes, 15 minutes reading time
Exam August 6, 208 Control Systems II (5-0590-00) Dr. Jacopo Tani Exam Exam Duration: 35 minutes, 5 minutes reading time Number of Problems: 35 Number of Points: 47 Permitted aids: 0 pages (5 sheets) A4.
More informationUsing Theorem Provers to Guarantee Closed-Loop Properties
Using Theorem Provers to Guarantee Closed-Loop Properties Nikos Aréchiga Sarah Loos André Platzer Bruce Krogh Carnegie Mellon University April 27, 2012 Aréchiga, Loos, Platzer, Krogh (CMU) Theorem Provers
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering Dynamics and Control II Fall 2007
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering.4 Dynamics and Control II Fall 7 Problem Set #9 Solution Posted: Sunday, Dec., 7. The.4 Tower system. The system parameters are
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 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 informationESC794: Special Topics: Model Predictive Control
ESC794: Special Topics: Model Predictive Control Nonlinear MPC Analysis : Part 1 Reference: Nonlinear Model Predictive Control (Ch.3), Grüne and Pannek Hanz Richter, Professor Mechanical Engineering Department
More informationGEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL & COMPUTER ENGINEERING FINAL EXAM. COURSE: ECE 3084A (Prof. Michaels)
GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL & COMPUTER ENGINEERING FINAL EXAM DATE: 30-Apr-14 COURSE: ECE 3084A (Prof. Michaels) NAME: STUDENT #: LAST, FIRST Write your name on the front page
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 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 informationLQR, Kalman Filter, and LQG. Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin
LQR, Kalman Filter, and LQG Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin May 2015 Linear Quadratic Regulator (LQR) Consider a linear system
More informationChap. 3. Controlled Systems, Controllability
Chap. 3. Controlled Systems, Controllability 1. Controllability of Linear Systems 1.1. Kalman s Criterion Consider the linear system ẋ = Ax + Bu where x R n : state vector and u R m : input vector. A :
More informationEE451/551: Digital Control. Chapter 8: Properties of State Space Models
EE451/551: Digital Control Chapter 8: Properties of State Space Models Equilibrium State Definition 8.1: An equilibrium point or state is an initial state from which the system nevers departs unless perturbed
More informationSeparation Principle & Full-Order Observer Design
Separation Principle & Full-Order Observer Design Suppose you want to design a feedback controller. Using full-state feedback you can place the poles of the closed-loop system at will. U Plant Kx If the
More informationLectures 25 & 26: Consensus and vehicular formation problems
EE 8235: Lectures 25 & 26 Lectures 25 & 26: Consensus and vehicular formation problems Consensus Make subsystems (agents, nodes) reach agreement Distributed decision making Vehicular formations How does
More informationExtreme Values and Positive/ Negative Definite Matrix Conditions
Extreme Values and Positive/ Negative Definite Matrix Conditions James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 8, 016 Outline 1
More informationAutonomous Mobile Robot Design
Autonomous Mobile Robot Design Topic: Guidance and Control Introduction and PID Loops Dr. Kostas Alexis (CSE) Autonomous Robot Challenges How do I control where to go? Autonomous Mobile Robot Design Topic:
More informationLinear-Quadratic-Gaussian (LQG) Controllers and Kalman Filters
Linear-Quadratic-Gaussian (LQG) Controllers and Kalman Filters Emo Todorov Applied Mathematics and Computer Science & Engineering University of Washington Winter 204 Emo Todorov (UW) AMATH/CSE 579, Winter
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 information