Linear Quadratic Regulator (LQR) I
|
|
- Avis Brianna Austin
- 6 years ago
- Views:
Transcription
1 Optimal Control, Guidance and Estimation Lecture Linear Quadratic Regulator (LQR) I Pro. Radhakant Padhi Dept. o Aerospace Engineering Indian Institute o Science - Bangalore Generic Optimal Control Problem Perormance Index (to minimize / maximize): Path Constraint: Xɺ = t X U Boundary Conditions: X = ϕ (, ) + (,, ) J t X L t X U dt (,, ) t t t = X : Speciied ( ) : Fixed, X t : Free Pro. Radhakant Padhi, AE Dept., IISc-Bangalore
2 LQR Design: Problem Objective o drive the state X o a linear (rather linearized) system Xɺ = A X + BU to the origin by minimizing the ollowing quadratic perormance index (cost unction) J = ( X S X ) + ( X Q X + U RU ) dt where t S, Q (psd), R > (pd) t Pro. Radhakant Padhi, AE Dept., IISc-Bangalore 3 LQR Design: Guideline or Selection o Weighting Matrices S (psd), Q (psd), R > (pd) hese are usually chosen as diagonal matrices, with s q r i i i = maximum expected/acceptable value o / ( xi ) = maximum expected/acceptable value o ( / xi ) = maximum expected/acceptable value o ( / u i ) Pro. Radhakant Padhi, AE Dept., IISc-Bangalore 4
3 LQR Design: Some Facts to Remember { } he pair A, B needs to be controllable and the pair A, Q needs to be detectable { } S (psd), Q (psd), R > (pd) (these are usually chosen as diagonal matrices) By deault, it is assumed that t Constrained problems (i.e. problems with state and control inequality constraints) are not considered here. hose will be considered later. Pro. Radhakant Padhi, AE Dept., IISc-Bangalore 5 LQR Design: Problem Statement Perormance Index (to minimize): Path Constraint: t J = ( X S X ) + X Q X + U RU dt t ϕ( X ) L( X, U ) Xɺ = A X + BU Boundary Conditions: X t = X : Speciied ( ) : Fixed, X t : Free Pro. Radhakant Padhi, AE Dept., IISc-Bangalore 6
4 LQR Design: Necessary Conditions o Optimality erminal penalty: ϕ ( X ) = ( X S X ) Hamiltonian: H = X Q X + U RU + AX + BU State Equation: Xɺ = AX + BU Costate Equation: Optimal Control Eq.: Boundary Condition: ( H / X ) ( QX A ) ɺ = = + ( / ) = = H U U R B ( ϕ / X ) = = S X Pro. Radhakant Padhi, AE Dept., IISc-Bangalore 7 LQR Design: Derivation o Riccati Equation Guess: Justiication: ( t) P ( t) X ( t) = From unctional analysis theory o normed linear space, lies in the "dual space" o X t o all continuous linear unctionals o X t. ( t), which is the space consisting Reerence: Optimization by Vector Space Methods D. G. Luenberger, John Wiley & Sons, 969. Pro. Radhakant Padhi, AE Dept., IISc-Bangalore 8
5 LQR Design: Derivation o Riccati Equation Guess ( t) = P ( t) X ( t) ɺ = PX ɺ + PXɺ = PX ɺ + P AX + BU = PX ɺ + P( AX BR B ) = PX ɺ + P( AX BR B PX ) ( QX + A PX ) = ( Pɺ + PA PBR B P) X ( Pɺ + PA + A P PBR B P + Q) X = Pro. Radhakant Padhi, AE Dept., IISc-Bangalore 9 LQR Design: Derivation o Riccati Equation Riccati equation ɺ = P PA A P PBR B P Q Boundary condition ( ) = ( is ree) P t X S X X ( ) P t = S Pro. Radhakant Padhi, AE Dept., IISc-Bangalore
6 LQR Design: Solution Procedure Use the boundary condition P t = S and integrate the Riccati Equation backwards rom to t t Store the solution history or the Riccati matrix Compute the optimal control online U = R B P X = K X Pro. Radhakant Padhi, AE Dept., IISc-Bangalore LQR Design: Ininite ime Regulator Problem heorem (By Kalman) As t, or constant Q and R matrices, Pɺ t Algebraic Riccati Equation (ARE) Note: PA A P PBR B P Q + + = ARE is still a nonlinear equation or the Riccati matrix. It is not straightorward to solve. However, eicient numerical methods are now available. A positive deinite solution or the Riccati matrix is needed to obtain a stabilizing controller. Pro. Radhakant Padhi, AE Dept., IISc-Bangalore
7 Example : Stabilization o Inverted Pendulum Pro. Radhakant Padhi Dept. o Aerospace Engineering Indian Institute o Science - Bangalore A Motivating Example: Stabilization o Inverted Pendulum System dynamics: ɺɺ θ = ω θ u, ω = g / L n n ( Linearized about vertical equilibrium point) u θ L mg System dynamics (state space orm): Deine: x = θ, x = ɺ θ xɺ x u x = ωn + x ɺ Xɺ A X B Perormance Index (to minimize): J = θ u dt + c Q =, R = c Pro. Radhakant Padhi, AE Dept., IISc-Bangalore 4
8 A Motivating Example: Stabilization o Inverted Pendulum ARE: PA A P PBR B P Q Let + + = P p p = p p 3 (a symmetric matrix) pω n p pωn p3ω n c p c p p p = 3ωn p p p c p p3 c p3 Equations: pωn c p + = p = 4 ω n ± ωn + c c p + p ω c p p = (repeated) 3 n 3 p c p = p = ± p c 3 3 u θ L mg Pro. Radhakant Padhi, AE Dept., IISc-Bangalore 5 A Motivating Example: Stabilization o Inverted Pendulum However, Hence, p p 3 4 n n 3 3 n n is a diagonal term, which needs to be real and positive. needs to be positive. hereore p = ω + ω + c, p = p c c Moreover, p + p ω c p p = p = c p p p ω (not needed in this problem) u θ L mg Gain Matrix: K = R B P = c p c p 3 Control: = = + ɺ ( θ 3θ ) u K X c p p Pro. Radhakant Padhi, AE Dept., IISc-Bangalore 6
9 A Motivating Example: Stabilization o Inverted Pendulum Analysis Open-Loop System: I A = = ωn = ω = ± ω n n ( right hal pole: unstable system) u θ L mg Closed-Loop System: Deine: ACL = A BK = ωn c p c p 3 Closed-Loop Poles: I A = CL ω = ω + c 4 n p = + c ( ωn ω ) p3 = p = + c c ( ωn ω ) / Pro. Radhakant Padhi, AE Dept., IISc-Bangalore 7 A Motivating Example: Stabilization o Inverted Pendulum Analysis Closed-Loop Poles: ( n ) / + ω + ω + ω =, = ω + ω ± ω ω 4 ( Note: ω = ωn + c > ωn ) ( n ) j ( n ) / / u θ L mg Both o the closed-loop poles are strictly in the let-hal plane. Hence, the closed-loop is guaranteed to be asymptotically stable. Pro. Radhakant Padhi, AE Dept., IISc-Bangalore 8
10 Example : Finite-time time emperature Control Pro. Radhakant Padhi Dept. o Aerospace Engineering Indian Institute o Science - Bangalore Example: Finite ime emperature Control Problem System dynamics : ɺ θ = a θ θ + bu where a, b : Constants θ θ u a a : emperature : Ambient temperature (Constant = C) : Heat input Pro. Radhakant Padhi, AE Dept., IISc-Bangalore
11 Problem ormulations Case : Case : Cost Function: Cost Function: J θ = ( t ) t t u dt J = s ( θ 3) + u dt = θ = 3 C (Hard constraint) s > : Weightage ( t ) i.e. θ 3 C (Sot Constraint) Pro. Radhakant Padhi, AE Dept., IISc-Bangalore Solution: Solution: x ( θ θa ), θ ( ) = θa ( θ θ ) xɺ = ax + bu, x = = Η = + + ɺ Η = = a x Η = u = b u u ax bu a Pro. Radhakant Padhi, AE Dept., IISc-Bangalore a Necessary conditions xɺ = ax + bu ɺ = a u = b
12 Solution: Case (Hard constraint) a( t t ) a( t t) = e = e u = be a t ( t) a t xɺ = ax b e ( t) aking laplace transorm: sx ( s) x( ) at ax s b e = s a Pro. Radhakant Padhi, AE Dept., IISc-Bangalore 3 Solution: Case (Hard constraint) s a at = b e a s a s + a Hence = a at X s = b e at at at x t b e e e Unknown ( θ θ ) However, x( t ) = = C a Pro. Radhakant Padhi, AE Dept., IISc-Bangalore 4
13 Solution: Case (Hard constraint) at ( e ) at ( e ) at at at x( t ) = = b e e e b = a a = b x( t) = b a at at a e e at at at e ( e e ) = b ( e ) a e e at at at Pro. Radhakant Padhi, AE Dept., IISc-Bangalore 5 Solution: Case (Hard constraint) Note : x( t ) = at ( e at e ) at ( e at e ) = (i.e. he boundary condition is "exactly met"!) Controller : u( t) = be at ( ) ae = b ( at ) ( at at e b e e ) a t t a Pro. Radhakant Padhi, AE Dept., IISc-Bangalore 6
14 Solution: Case (Sot constraint) θ 3 C x C. Hence the cost unction is t J = s ( x ) + u dt = s ( x ) x = + s However, we have b x t e e e a at at at = ( ) Pro. Radhakant Padhi, AE Dept., IISc-Bangalore 7 Solution: Case (Sot constraint) At t = t, x t = e = + a s b at b at i. e. + ( e ) = s a s a i. e. = at a + s b ( e ) Hence a( t ) t a t t s a = e = e a + s b e at Pro. Radhakant Padhi, AE Dept., IISc-Bangalore 8
15 Solution: Case (Sot constraint) u t = b = be = be s a a + s b e a t t at at s abe s b ae + e e a t t at at at Pro. Radhakant Padhi, AE Dept., IISc-Bangalore 9 Correlation between hard and sot constraint results As s, at abe lim u( t) = lim b ae + e e s s ( e ) S. C. s at at at at ae = = u t b e at at H. C. i. e. he "sot constraint" problem behaves like the "hard constraint" problem when s. Pro. Radhakant Padhi, AE Dept., IISc-Bangalore 3
16 hanks or the Attention.!! Pro. Radhakant Padhi, AE Dept., IISc-Bangalore 3
Linear Quadratic Regulator (LQR) Design I
Lecture 7 Linear Quadratic Regulator LQR) Design I Dr. Radhakant Padhi Asst. Proessor Dept. o Aerospace Engineering Indian Institute o Science - Bangalore LQR Design: Problem Objective o drive the state
More informationLinear Quadratic Regulator (LQR) II
Optimal Control, Guidance and Estimation Lecture 11 Linear Quadratic Regulator (LQR) II Pro. Radhakant Padhi Dept. o Aerospace Engineering Indian Institute o Science - Bangalore Outline Summary o LQR design
More informationLinear Quadratic Regulator (LQR) Design II
Lecture 8 Linear Quadratic Regulator LQR) Design II Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Outline Stability and Robustness properties
More informationConstrained Optimal Control I
Optimal Control, Guidance and Estimation Lecture 34 Constrained Optimal Control I Pro. Radhakant Padhi Dept. o Aerospace Engineering Indian Institute o Science - Bangalore opics Motivation Brie Summary
More informationConstrained Optimal Control. Constrained Optimal Control II
Optimal Control, Guidance and Estimation Lecture 35 Constrained Optimal Control II Prof. Radhakant Padhi Dept. of Aerospace Engineering Indian Institute of Science - Bangalore opics: Constrained Optimal
More informationUsing Lyapunov Theory I
Lecture 33 Stability Analysis of Nonlinear Systems Using Lyapunov heory I Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Outline Motivation Definitions
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 informationFeedback Optimal Control for Inverted Pendulum Problem by Using the Generating Function Technique
(IJACSA) International Journal o Advanced Computer Science Applications Vol. 5 No. 11 14 Feedback Optimal Control or Inverted Pendulum Problem b Using the Generating Function echnique Han R. Dwidar Astronom
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 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 informationGain Scheduling and Dynamic Inversion
Lecture 31 Gain Scheduling and Dynamic Inversion Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Topics A Brief Overview of Gain Scheduling Philosophy
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 informationChapter 2 SDOF Vibration Control 2.1 Transfer Function
Chapter SDOF Vibration Control.1 Transfer Function mx ɺɺ( t) + cxɺ ( t) + kx( t) = F( t) Defines the transfer function as output over input X ( s) 1 = G( s) = (1.39) F( s) ms + cs + k s is a complex number:
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 # /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 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 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 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 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 information4F3 - Predictive Control
4F3 Predictive Control - Lecture 2 p 1/23 4F3 - Predictive Control Lecture 2 - Unconstrained Predictive Control Jan Maciejowski jmm@engcamacuk 4F3 Predictive Control - Lecture 2 p 2/23 References Predictive
More informationLecture 6 Classical Control Overview IV. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore
Lecture 6 Classical Control Overview IV Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Lead Lag Compensator Design Dr. Radhakant Padhi Asst.
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 information4F3 - Predictive Control
4F3 Predictive Control - Lecture 3 p 1/21 4F3 - Predictive Control Lecture 3 - Predictive Control with Constraints Jan Maciejowski jmm@engcamacuk 4F3 Predictive Control - Lecture 3 p 2/21 Constraints on
More informationFeedback Linearization
Feedback Linearization Peter Al Hokayem and Eduardo Gallestey May 14, 2015 1 Introduction Consider a class o single-input-single-output (SISO) nonlinear systems o the orm ẋ = (x) + g(x)u (1) y = h(x) (2)
More informationContents. 1 State-Space Linear Systems 5. 2 Linearization Causality, Time Invariance, and Linearity 31
Contents Preamble xiii Linear Systems I Basic Concepts 1 I System Representation 3 1 State-Space Linear Systems 5 1.1 State-Space Linear Systems 5 1.2 Block Diagrams 7 1.3 Exercises 11 2 Linearization
More 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 informationLecture 5 Classical Control Overview III. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore
Lecture 5 Classical Control Overview III Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore A Fundamental Problem in Control Systems Poles of open
More informationLecture 9: Discrete-Time Linear Quadratic Regulator Finite-Horizon Case
Lecture 9: Discrete-Time Linear Quadratic Regulator Finite-Horizon Case Dr. Burak Demirel Faculty of Electrical Engineering and Information Technology, University of Paderborn December 15, 2015 2 Previous
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 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 informationThe Inverted Pendulum
Lab 1 The Inverted Pendulum Lab Objective: We will set up the LQR optimal control problem for the inverted pendulum and compute the solution numerically. Think back to your childhood days when, for entertainment
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 informationOptimal Control of Nonlinear Systems using RBF Neural Network and Adaptive Extended Kalman Filter
9 American Control Conerence Hyatt Regency Riverront, St. Louis, MO, USA June -, 9 WeA. Optimal Control o Nonlinear Systems using RBF Neural Network and Adaptive Extended Kalman Filter Peda V. Medagam,
More informationLinear control of inverted pendulum
Linear control of inverted pendulum Deep Ray, Ritesh Kumar, Praveen. C, Mythily Ramaswamy, J.-P. Raymond IFCAM Summer School on Numerics and Control of PDE 22 July - 2 August 213 IISc, Bangalore http://praveen.cfdlab.net/teaching/control213
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 information6. Linear Quadratic Regulator Control
EE635 - Control System Theory 6. Linear Quadratic Regulator Control Jitkomut Songsiri algebraic Riccati Equation (ARE) infinite-time LQR (continuous) Hamiltonian matrix gain margin of LQR 6-1 Algebraic
More informationECE7850 Lecture 7. Discrete Time Optimal Control and Dynamic Programming
ECE7850 Lecture 7 Discrete Time Optimal Control and Dynamic Programming Discrete Time Optimal control Problems Short Introduction to Dynamic Programming Connection to Stabilization Problems 1 DT nonlinear
More informationLecture 2: Discrete-time Linear Quadratic Optimal Control
ME 33, U Berkeley, Spring 04 Xu hen Lecture : Discrete-time Linear Quadratic Optimal ontrol Big picture Example onvergence of finite-time LQ solutions Big picture previously: dynamic programming and finite-horizon
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 informationLecture : Feedback Linearization
ecture : Feedbac inearization Niola Misovic, dipl ing and Pro Zoran Vuic June 29 Summary: This document ollows the lectures on eedbac linearization tought at the University o Zagreb, Faculty o Electrical
More informationClassical Numerical Methods to Solve Optimal Control Problems
Lecture 26 Classical Numerical Methods to Solve Optimal Control Problems Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Necessary Conditions
More information1 Relative degree and local normal forms
THE ZERO DYNAMICS OF A NONLINEAR SYSTEM 1 Relative degree and local normal orms The purpose o this Section is to show how single-input single-output nonlinear systems can be locally given, by means o a
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 informationAustralian Journal of Basic and Applied Sciences, 3(4): , 2009 ISSN Modern Control Design of Power System
Australian Journal of Basic and Applied Sciences, 3(4): 4267-4273, 29 ISSN 99-878 Modern Control Design of Power System Atef Saleh Othman Al-Mashakbeh Tafila Technical University, Electrical Engineering
More informationSolutions for Homework #8. Landing gear
Solutions or Homewor #8 PROBEM. (P. 9 on page 78 in the note) An airplane is modeled as a beam with masses as shown below: m m m m π [rad/sec] anding gear m m.5 Find the stiness and mass matrices. Find
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 informationHere represents the impulse (or delta) function. is an diagonal matrix of intensities, and is an diagonal matrix of intensities.
19 KALMAN FILTER 19.1 Introduction In the previous section, we derived the linear quadratic regulator as an optimal solution for the fullstate feedback control problem. The inherent assumption was that
More informationDynamic Inversion Design II
Lecture 32 Dynamic Inversion Design II Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Topics Summary of Dynamic Inversion Design Advantages Issues
More informationTime-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 2015
Time-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 15 Asymptotic approach from time-varying to constant gains Elimination of cross weighting
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 informationNumerics and Control of PDEs Lecture 1. IFCAM IISc Bangalore
1/1 Numerics and Control of PDEs Lecture 1 IFCAM IISc Bangalore July 22 August 2, 2013 Introduction to feedback stabilization Stabilizability of F.D.S. Mythily R., Praveen C., Jean-Pierre R. 2/1 Q1. Controllability.
More informationTime-Invariant Linear Quadratic Regulators!
Time-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 17 Asymptotic approach from time-varying to constant gains Elimination of cross weighting
More informationOptimization Tutorial 1. Basic Gradient Descent
E0 270 Machine Learning Jan 16, 2015 Optimization Tutorial 1 Basic Gradient Descent Lecture by Harikrishna Narasimhan Note: This tutorial shall assume background in elementary calculus and linear algebra.
More informationMCE/EEC 647/747: Robot Dynamics and Control. Lecture 8: Basic Lyapunov Stability Theory
MCE/EEC 647/747: Robot Dynamics and Control Lecture 8: Basic Lyapunov Stability Theory Reading: SHV Appendix Mechanical Engineering Hanz Richter, PhD MCE503 p.1/17 Stability in the sense of Lyapunov A
More informationStatic and Dynamic Optimization (42111)
Static and Dynamic Optimization (421) Niels Kjølstad Poulsen Build. 0b, room 01 Section for Dynamical Systems Dept. of Applied Mathematics and Computer Science The Technical University of Denmark Email:
More informationOptimal control of nonlinear systems with input constraints using linear time varying approximations
ISSN 1392-5113 Nonlinear Analysis: Modelling and Control, 216, Vol. 21, No. 3, 4 412 http://dx.doi.org/1.15388/na.216.3.7 Optimal control of nonlinear systems with input constraints using linear time varying
More informationOptimal Control. Quadratic Functions. Single variable quadratic function: Multi-variable quadratic function:
Optimal Control Control design based on pole-placement has non unique solutions Best locations for eigenvalues are sometimes difficult to determine Linear Quadratic LQ) Optimal control minimizes a quadratic
More informationReverse Order Swing-up Control of Serial Double Inverted Pendulums
Reverse Order Swing-up Control of Serial Double Inverted Pendulums T.Henmi, M.Deng, A.Inoue, N.Ueki and Y.Hirashima Okayama University, 3-1-1, Tsushima-Naka, Okayama, Japan inoue@suri.sys.okayama-u.ac.jp
More informationDesign and Comparison of Different Controllers to Stabilize a Rotary Inverted Pendulum
ISSN (Online): 347-3878, Impact Factor (5): 3.79 Design and Comparison of Different Controllers to Stabilize a Rotary Inverted Pendulum Kambhampati Tejaswi, Alluri Amarendra, Ganta Ramesh 3 M.Tech, Department
More informationDissipativity. Outline. Motivation. Dissipative Systems. M. Sami Fadali EBME Dept., UNR
Dissipativity M. Sami Fadali EBME Dept., UNR 1 Outline Differential storage functions. QSR Dissipativity. Algebraic conditions for dissipativity. Stability of dissipative systems. Feedback Interconnections
More informationLecture 6. Foundations of LMIs in System and Control Theory
Lecture 6. Foundations of LMIs in System and Control Theory Ivan Papusha CDS270 2: Mathematical Methods in Control and System Engineering May 4, 2015 1 / 22 Logistics hw5 due this Wed, May 6 do an easy
More informationTheory in Model Predictive Control :" Constraint Satisfaction and Stability!
Theory in Model Predictive Control :" Constraint Satisfaction and Stability Colin Jones, Melanie Zeilinger Automatic Control Laboratory, EPFL Example: Cessna Citation Aircraft Linearized continuous-time
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 informationRobust feedback linearization
Robust eedback linearization Hervé Guillard Henri Bourlès Laboratoire d Automatique des Arts et Métiers CNAM/ENSAM 21 rue Pinel 75013 Paris France {herveguillardhenribourles}@parisensamr Keywords: Nonlinear
More informationProducts and Convolutions of Gaussian Probability Density Functions
Tina Memo No. 003-003 Internal Report Products and Convolutions o Gaussian Probability Density Functions P.A. Bromiley Last updated / 9 / 03 Imaging Science and Biomedical Engineering Division, Medical
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 informationOPTIMAL CONTROL AND ESTIMATION
OPTIMAL CONTROL AND ESTIMATION Robert F. Stengel Department of Mechanical and Aerospace Engineering Princeton University, Princeton, New Jersey DOVER PUBLICATIONS, INC. New York CONTENTS 1. INTRODUCTION
More informationLyapunov Stability Analysis: Open Loop
Copyright F.L. Lewis 008 All rights reserved Updated: hursday, August 8, 008 Lyapunov Stability Analysis: Open Loop We know that the stability of linear time-invariant (LI) dynamical systems can be determined
More informationLecture 10 Linear Quadratic Stochastic Control with Partial State Observation
EE363 Winter 2008-09 Lecture 10 Linear Quadratic Stochastic Control with Partial State Observation partially observed linear-quadratic stochastic control problem estimation-control separation principle
More information522. Optimal robust control of aeroelastic system vibrations
522. Optimal robust control of aeroelastic system vibrations K. Bousson University of Beira Interior, Avionics and Control Laboratory, Department of Aerospace Sciences 621-1 Covilhã Portugal E-mail: k1bousson@yahoo.com,
More information18-660: Numerical Methods for Engineering Design and Optimization
8-66: Numerical Methods or Engineering Design and Optimization Xin Li Department o ECE Carnegie Mellon University Pittsburgh, PA 53 Slide Overview Linear Regression Ordinary least-squares regression Minima
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 informationSemidefinite Programming Basics and Applications
Semidefinite Programming Basics and Applications Ray Pörn, principal lecturer Åbo Akademi University Novia University of Applied Sciences Content What is semidefinite programming (SDP)? How to represent
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 informationHamilton-Jacobi-Bellman Equation Feb 25, 2008
Hamilton-Jacobi-Bellman Equation Feb 25, 2008 What is it? The Hamilton-Jacobi-Bellman (HJB) equation is the continuous-time analog to the discrete deterministic dynamic programming algorithm Discrete VS
More informationEG4321/EG7040. Nonlinear Control. Dr. Matt Turner
EG4321/EG7040 Nonlinear Control Dr. Matt Turner EG4321/EG7040 [An introduction to] Nonlinear Control Dr. Matt Turner EG4321/EG7040 [An introduction to] Nonlinear [System Analysis] and Control Dr. Matt
More informationExample: When describing where a function is increasing, decreasing or constant we use the x- axis values.
Business Calculus Lecture Notes (also Calculus With Applications or Business Math II) Chapter 3 Applications o Derivatives 31 Increasing and Decreasing Functions Inormal Deinition: A unction is increasing
More informationEE C128 / ME C134 Feedback Control Systems
EE C128 / ME C134 Feedback Control Systems Lecture Additional Material Introduction to Model Predictive Control Maximilian Balandat Department of Electrical Engineering & Computer Science University of
More informationLinearization problem. The simplest example
Linear Systems Lecture 3 1 problem Consider a non-linear time-invariant system of the form ( ẋ(t f x(t u(t y(t g ( x(t u(t (1 such that x R n u R m y R p and Slide 1 A: f(xu f(xu g(xu and g(xu exist and
More informationOutline. 1 Linear Quadratic Problem. 2 Constraints. 3 Dynamic Programming Solution. 4 The Infinite Horizon LQ Problem.
Model Predictive Control Short Course Regulation James B. Rawlings Michael J. Risbeck Nishith R. Patel Department of Chemical and Biological Engineering Copyright c 217 by James B. Rawlings Outline 1 Linear
More informationIntroduction to centralized control
ROBOTICS 01PEEQW Basilio Bona DAUIN Politecnico di Torino Control Part 2 Introduction to centralized control Independent joint decentralized control may prove inadequate when the user requires high task
More informationMulti-Model Adaptive Regulation for a Family of Systems Containing Different Zero Structures
Preprints of the 19th World Congress The International Federation of Automatic Control Multi-Model Adaptive Regulation for a Family of Systems Containing Different Zero Structures Eric Peterson Harry G.
More informationComparison of LQR and PD controller for stabilizing Double Inverted Pendulum System
International Journal of Engineering Research and Development ISSN: 78-67X, Volume 1, Issue 1 (July 1), PP. 69-74 www.ijerd.com Comparison of LQR and PD controller for stabilizing Double Inverted Pendulum
More informationStabilization of a 3D Rigid Pendulum
25 American Control Conference June 8-, 25. Portland, OR, USA ThC5.6 Stabilization of a 3D Rigid Pendulum Nalin A. Chaturvedi, Fabio Bacconi, Amit K. Sanyal, Dennis Bernstein, N. Harris McClamroch Department
More informationSemidefinite Programming Duality and Linear Time-invariant Systems
Semidefinite Programming Duality and Linear Time-invariant Systems Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University 2 July 2004 Workshop on Linear Matrix Inequalities in Control LAAS-CNRS,
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 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 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 informationy,z the subscript y, z indicating that the variables y and z are kept constant. The second partial differential with respect to x is written x 2 y,z
8 Partial dierentials I a unction depends on more than one variable, its rate o change with respect to one o the variables can be determined keeping the others ied The dierential is then a partial dierential
More informationLecture 8 Optimization
4/9/015 Lecture 8 Optimization EE 4386/5301 Computational Methods in EE Spring 015 Optimization 1 Outline Introduction 1D Optimization Parabolic interpolation Golden section search Newton s method Multidimensional
More informationDenis ARZELIER arzelier
COURSE ON LMI OPTIMIZATION WITH APPLICATIONS IN CONTROL PART II.2 LMIs IN SYSTEMS CONTROL STATE-SPACE METHODS PERFORMANCE ANALYSIS and SYNTHESIS Denis ARZELIER www.laas.fr/ arzelier arzelier@laas.fr 15
More informationSIO 211B, Rudnick. We start with a definition of the Fourier transform! ĝ f of a time series! ( )
SIO B, Rudnick! XVIII.Wavelets The goal o a wavelet transorm is a description o a time series that is both requency and time selective. The wavelet transorm can be contrasted with the well-known and very
More informationSwitched systems: stability
Switched systems: stability OUTLINE Switched Systems Stability of Switched Systems OUTLINE Switched Systems Stability of Switched Systems a family of systems SWITCHED SYSTEMS SWITCHED SYSTEMS a family
More informationOptimization of Linear Systems of Constrained Configuration
Optimization of Linear Systems of Constrained Configuration Antony Jameson 1 October 1968 1 Abstract For the sake of simplicity it is often desirable to restrict the number of feedbacks in a controller.
More informationChapter 8 Stabilization: State Feedback 8. Introduction: Stabilization One reason feedback control systems are designed is to stabilize systems that m
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A. Dahleh George Verghese Department of Electrical Engineering and Computer Science Massachuasetts Institute of echnology c Chapter 8 Stabilization:
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 informationEL2520 Control Theory and Practice
EL2520 Control Theory and Practice Lecture 8: Linear quadratic control Mikael Johansson School of Electrical Engineering KTH, Stockholm, Sweden Linear quadratic control Allows to compute the controller
More informationSWEEP METHOD IN ANALYSIS OPTIMAL CONTROL FOR RENDEZ-VOUS PROBLEMS
J. Appl. Math. & Computing Vol. 23(2007), No. 1-2, pp. 243-256 Website: http://jamc.net SWEEP METHOD IN ANALYSIS OPTIMAL CONTROL FOR RENDEZ-VOUS PROBLEMS MIHAI POPESCU Abstract. This paper deals with determining
More information10. Joint Moments and Joint Characteristic Functions
10. Joint Moments and Joint Characteristic Functions Following section 6, in this section we shall introduce various parameters to compactly represent the inormation contained in the joint p.d. o two r.vs.
More informationMath Review and Lessons in Calculus
Math Review and Lessons in Calculus Agenda Rules o Eponents Functions Inverses Limits Calculus Rules o Eponents 0 Zero Eponent Rule a * b ab Product Rule * 3 5 a / b a-b Quotient Rule 5 / 3 -a / a Negative
More information