H 2 Optimal State Feedback Control Synthesis. Raktim Bhattacharya Aerospace Engineering, Texas A&M University
|
|
- Jodie Short
- 5 years ago
- Views:
Transcription
1 H 2 Optimal State Feedback Control Synthesis Raktim Bhattacharya Aerospace Engineering, Texas A&M University
2 Motivation
3 Motivation w(t) u(t) G K y(t) z(t) w(t) are exogenous signals reference, process noise, sensor noise z(t) are signals we want to keep small e(t), u(t) Find controller K that minimizes G w z 2 Appropriate when spectral density function of w(t) is known For example w(t) can be stationary noise H2 optimal is then linear quadratic Gaussian (LQG) AERO 632, Instructor: Raktim Bhattacharya 3 / 38
4 Motivation Stochastic Input For stationary stochastic process w(t), autocorrelation matrix is R w (τ) := E [w(t + τ)w (t)] The Fourier transform of R w (τ) is the spectral density Ŝw(jω) Signal power is related to it by E [ w(t) ] 2 := tr [R w ()] = 1 ] tr [Ŝw (jω) dω 2π If z and w are related by z = P w, for a stable LTI system P, then Ŝ z (jω) = ˆP (jω)ŝw(jω) ˆP (jω) AERO 632, Instructor: Raktim Bhattacharya 4 / 38
5 Motivation Stochastic Input (contd.) Therefore, E [ z(t) ] 2 = 1 2π [ tr ˆP (jω) Ŝ w (jω) ˆP ] (jω) dω. Looks like weighted H 2 norm of P. If w(t) is white noise = Ŝw(jω) = I. H 2 norm is output variance with white noise input. For any other Ŝw(jω), Therefore, E [ z(t) ] 2 = 1 2π Ŝ w (jω) = Ŵ (jω)ŵ (jω). Think of ˆP (jω)ŵ (jω) as weighted system [( ) tr ˆP (jω) Ŵ (jω) I (Ŵ (jω) ˆP )] (jω) dω. AERO 632, Instructor: Raktim Bhattacharya 5 / 38
6 Motivation K Stochastic Input (contd.) P z(t) w(t) z(t) P w(t) z(t) P w(t) W white noise E [ z(t) ] 2 = 1 P z(t) 2π w(t) [( ) tr ˆP (jω) Ŵ (jω) W white noise Think of ˆP (jω)ŵ (jω) as weighted system I (Ŵ (jω) ˆP )] (jω) dω. AERO 632, Instructor: Raktim Bhattacharya 6 / 38
7 Motivation Impulse Response Input signal is known in advance Tracking a fixed signal e.g. step Consider a special case: Scalar w(t) = δ(t). Implies z 2 2 = = 1 2π = 1 2π = ˆP 2 2 z (t)z(t)dt ẑ (jω)ẑ(jω)dω ˆP (jω) ˆP (jω)dω AERO 632, Instructor: Raktim Bhattacharya 7 / 38
8 P W z(t) w(t) white noise Motivation Impulse Response (contd.) What about any other referencez(t) r(t)? P r(t) z(t) P r(t) z(t) P r(t) h(t) (t) r(t) can bepreplaced byh(t) the impulse response h(t) of a known z(t) r(t) (t) filter W (jω) z(t) becomes impulse response of a weighted plant 1 AERO 632, Instructor: Raktim Bhattacharya 8 / 38
9 Computation of H 2 Norm Best computed in state-space realization of system State Space Model: General MIMO LTI system modeled as ẋ = Ax + Bu, y = Cx + Du, where x R n, u R m, y R p. Transfer Function Ĝ(s) = D + C(sI A) 1 B strictly proper when D = Impulse Response G(t) = L 1 { C(sI A) 1 B } = Ce ta B. AERO 632, Instructor: Raktim Bhattacharya 9 / 38
10 H 2 Norm MIMO Systems Ĝ(jω) 2 2 = 1 [Ĝ ] tr 2π (jω)ĝ(jω) for matrix transfer function = G(t) 2 2 Parseval [ ] = tr Ce ta BB T e tat C T dt ( = tr C ) e ta BB T e tat dt } {{ } L c C T Lc = controllability Gramian = tr [ CL c C T ] AERO 632, Instructor: Raktim Bhattacharya 1 / 38
11 H 2 Norm (contd.) MIMO Systems For any matrix M tr [M M] = tr [MM ] = Ĝ(jω) 2 2 = tr BT ( ) e tat C T Ce ta dt B }{{} L o = tr [ B T L o B ] L o = observability Gramian H 2 Norm of Ĝ(jω) Ĝ(jω) 2 2 = tr [ CL c C T ] = tr [ B T L o B ]. AERO 632, Instructor: Raktim Bhattacharya 11 / 38
12 H 2 Norm How to determine L c and L o? They are solutions of the following equation AL c + L c A T + BB T =, A T L o + L o A + C T C =. AERO 632, Instructor: Raktim Bhattacharya 12 / 38
13 Proposition 1 Suppose P is a state-space system with realization (A, B, C). Then iff X = X T > such that A is Hurwitz and ˆP 2 1, tr [CXC ] < 1 and AX + XA + BB <. Proof Only If Recall ˆP 2 = tr [CL c C ]. Therefore for X = L c and A is Hurwitz ˆP 2 is finite tr [CL c C ] < 1 = ˆP 2 < 1, AERO 632, Instructor: Raktim Bhattacharya 13 / 38
14 contd. Consider X = e ta (BB + ϵi n )e ta dt, is continuous in ϵ equals to L c when ϵ =. It can be shown that this X satisfies Lyapunov equation follow the derivation AX + XA + BB + ϵi n =, or AX + XA + BB <. AERO 632, Instructor: Raktim Bhattacharya 14 / 38
15 contd. Proof If X satisfies and Implies, A is Hurwitz. AX + XA + BB < tr [CXC ] < 1. Proposition: It can be shown that if L c satisfies and X satisfies A L c + L c A + Q =, then X L c. A X + XA + Q, AERO 632, Instructor: Raktim Bhattacharya 15 / 38
16 contd. Inequality X L c implies ˆP 2 = tr [CL c C ] tr [CXC ] < 1. AERO 632, Instructor: Raktim Bhattacharya 16 / 38
17 System Dynamics w(t) u(t) G y(t) z(t) K Dynamics ẋ = Ax + B w w + B u u z = C z x + D u u y = x With u = Kx, F l (Ĝ, K) := Ĝw z(s) = = Ĝ(s) = [ A + Bu K B w C z + D u K A B w B u C z D u I ]. AERO 632, Instructor: Raktim Bhattacharya 17 / 38
18 Controller Synthesis Proposition There exists feedback gain K that internally stabilizes G and satisfies F l (Ĝ, K) 2 < 1 iff Z R m n such that K = ZX 1, where X > satisfies inequalities [ ] [ ] X A Bu + [ [ X Z Z ] A ] B u + B w B w <, tr [ (C z X + D u Z)X 1 (C z X + D u Z) ] < 1. Proof: Use closed-loop state-space data and earlier proposition. Not an LMI. AERO 632, Instructor: Raktim Bhattacharya 18 / 38
19 Controller Synthesis (contd.) Apply Schur complement to get following convex problem. There exists feedback gain K that internally stabilizes G satisfying F l (Ĝ, K) 2 < 1 iff X R n n, W R q q, and Z R m n, such that and K = ZX 1, min tr [W ] W [ ] [ ] X A Bu + [ [ X Z Z ] A ] Bu + B w Bw <, [ ] W (C z X + D u Z) (C z X + D u Z) > X AERO 632, Instructor: Raktim Bhattacharya 19 / 38
20 Very simple example System Dynamics ẋ = x + 2 u + w, [ ] [ ] z = x + u Compare with LQR AERO 632, Instructor: Raktim Bhattacharya 2 / 38
21 Very simple example 6 4 H 2 LQR Output z Time(s) Output z Time(s) tr [W ] = e 4 Good disturbance attenuation. AERO 632, Instructor: Raktim Bhattacharya 21 / 38
22 MATLAB Code clear; clc; A = [-3-2 1; 1 2 1; ;]; Bu = [2 ; 2; 1]; Bw = [3;;1]; Cz = [1 1; 1 1]; Du = [1 1; 1]; nx = 3; nu = 2; nz = 2; nw = 1; cvx_begin sdp variable X(nx,nx) symmetric variable W(nz,nz) symmetric variable Z(nu,nx) [A Bu]*[X;Z] + [X Z ]*[A ;Bu ] + Bw*Bw < [W (Cz*X + Du*Z); (Cz*X + Du*Z) X] > minimize trace(w) cvx_end h2k = Z*inv(X); [lqrk,s,e] = lqr(a,bu,cz *Cz,Du *Du); h2g = ss(a+bu*h2k, Bw, Cz + Du*h2K, zeros(nw,1)) lqrg = ss(a-bu*lqrk, Bw, Cz + Du*lqrK,... zeros(nw,1)); T = [:.1:2]; w = 5*randn(length(T),1); [y1,t1,x1] = lsim(h2g,w,t); [y2,t2,x2] = lsim(lqrg,w,t); f1=figure(1); clf; set(f1, defaulttextinterpreter, latex ); for i=1:2 subplot(2,1,i); plot(t1,y1(:,i), r,t2,y2(:,i), b:,... linewidth,2); xlabel( Time(s) ); ylabel(sprintf( Output $z_%d$,i)); end subplot(2,1,1); legend( H_2, LQR ); print -depsc h2ex1.eps AERO 632, Instructor: Raktim Bhattacharya 22 / 38
23 Regulator with disturbance Active Suspension Equations of motion k1 m1 m2 u(t) x1(t) c1 x2(t) m 1 ẍ 1 = k 1 (x 1 x 2 ) c 1 (ẋ 1 ẋ 2 ) + u, m 2 ẍ 2 = k 1 (x 1 x 2 ) + c 1 (ẋ 1 ẋ 2 ) k 2 (x 2 d) u. State Variables k2 q 1 := x 1, q 2 := x 2, q 3 := ẋ 1, q 4 := ẋ 2. d(t) AERO 632, Instructor: Raktim Bhattacharya 23 / 38
24 Regulator with disturbance Active Suspension Linear System q 1 1 q 1 q 2 q 3 = 1 q 2 k 1 /m 1 k 1 /m 1 c 1 /m 1 c 1 /m 1 q 3 + q 4 k 1 /m 2 (k 1 + k 2 )/m 2 c 1 /m 2 c 1 /m 2 q 4 1/m 1 u + d. 1/m 2 k 2 /m 2 Output z(t) z = q 1 + u. AERO 632, Instructor: Raktim Bhattacharya 24 / 38
25 Regulator with disturbance Active Suspension.8 d(t) = white noise [ 1, 1] m H2 LQR.8 d(t) =.5 sin (5 T) m H2 LQR output z(t) output z(t) Time Time tr [W ] = e 1 Very good disturbance rejection. AERO 632, Instructor: Raktim Bhattacharya 25 / 38
26 MATLAB Code clear; clc; % System Parameters m1 = 29; % kg -- Body mass m2 = 6; % kg -- suspension mass k1 = 162; % N/m k2 = 191; % N/m c1 = 1; % Ns/m A = [ 1 ; 1; -k1/m1 k1/m1 -c1/m1 c1/m1; k1/m1 -(k1+k2)/m2 c1/m2 -c1/m2]; Bu = [;;1/m1;-1/m2]; Bw = [;;;k2/m2]; nx = 4; nu = 1; nz = 1; nw = 1; Cz = [1,,,]; Du = 1*ones(nz,nu); cvx_begin sdp variable X(nx,nx) symmetric variable W(nz,nz) symmetric variables Z(nu,nx) gam [A Bu]*[X;Z] + [X Z ]*[A ;Bu ] + Bw*Bw < [W (Cz*X + Du*Z); (Cz*X + Du*Z) X] > minimize trace(w) cvx_end h2k = Z*inv(X); [lqrk,s,e] = lqr(a,bu,cz *Cz,Du *Du); % Simulation h2g = ss(a+bu*h2k, Bw, Cz + Du*h2K, zeros(nz,nw)); lqrg = ss(a-bu*lqrk, Bw, Cz + Du*lqrK, zeros(nz,nw)); T = [:.1:2]/1; w = 2*rand(length(T),1)-1; w =.5*sin(1*T); [y1,t1,x1] = lsim(h2g,w,t); [y2,t2,x2] = lsim(lqrg,w,t); f1 = figure(1); clf; set(f1, defaulttextinterpreter, latex ); plot(t1,y1, r,t2,y2, b:, Linewidth,2); xlabel( Time ); ylabel( output $z(t)$ ); title( $d(t) =.5*\sin\;\; (5*T)$ m ); title( $d(t)$ = white noise $\in [-1,\;1]$ m ); legend( H2, LQR ); print -depsc h2qcar2.eps AERO 632, Instructor: Raktim Bhattacharya 26 / 38
27 Tracking with disturbance Let dynamical system be ẋ = Ax + B u u + B d d, z = [ ] x, y = C u y x. Design a controller of the form u = Kx + Gr K is designed in H 2 optimal sense G is regular tracking gain AERO 632, Instructor: Raktim Bhattacharya 27 / 38
28 Tracking with disturbance Closed-loop system with K is therefore ẋ = (A + BK)x + B u Gr, y = C y x. Ignore disturbance when determining G Steady-state response to constant r is = (A + BK)x ss + BGr, y ss = C y x ss = r. Or x ss = (A + BK) 1 BGr, implies C y (A + BK) 1 BGr = r, or C y (A + BK) 1 BG = I. Solve for G AERO 632, Instructor: Raktim Bhattacharya 28 / 38
29 Tracking with disturbance Existence of solution of C y (A + BK) 1 BG = I is necessary and sufficient condition for existence of a tracking controller. In general, when y = C y x + D y u, the equation becomes [ Dy (C y + D y K)(A + BK) 1 B ] G = I. Can be rewritten in terms of Π Π = (A + BK) 1 BG = (C y + D y K)Π + D y G = I, AERO 632, Instructor: Raktim Bhattacharya 29 / 38
30 Tracking with disturbance Or Or rearranged to (A + BK)Π + BG =, (C y + D y K)Π + D y G = I. AΠ + B(KΠ + G) =, C y Π + D y (KΠ + G) = I. Therefore, with Γ := KΠ + G, we get [ ] [ ] [ ] A B Π nx n + y C y D y Γ I ny ny =. This is the so called regulator equation. Get G = Γ KΠ. AERO 632, Instructor: Raktim Bhattacharya 3 / 38
31 Tracking with disturbance Solve regulator equation [ A B C y D y ] [ Π Γ ] + [ nx n y I ny n y ] =. Get G = Γ KΠ. Control law is u = Kx + Gr. AERO 632, Instructor: Raktim Bhattacharya 31 / 38
32 Tracking with disturbance Longitudinal F16 Control Law Longitudinal Motion States [V (ft/s) α(rad) θ(rad) q(rad/s)] T Controls [T (lb) δ e (deg)] Constraints on u [ ] 1 u 25 Trimmed at steady-level flight at 932 ft/s x trim = Disturbance in α equation [ ], u trim = [ ] [ ] AERO 632, Instructor: Raktim Bhattacharya 32 / 38
33 Tracking with disturbance Longitudinal F16 Control Law d(t) System Interconnection r G d x K + F 16 Cy y Control Objective Design K to minimize G d x 2 [ ] V Design G to track r := reference 1γ AERO 632, Instructor: Raktim Bhattacharya 33 / 38
34 Tracking with disturbance Longitudinal F16 Control Law Disturbance Rejection d(t) U [ 1,1] rad, tr [W ] = V ft/s Time AERO 632, Instructor: Raktim Bhattacharya 34 / 38
35 Tracking with disturbance Longitudinal F16 Control Law Tracking Performance V ref = 5 ft/s step, γ ref = V (ft/s) Gamma (deg) Thrust (lb) Elevator (deg) Time Time AERO 632, Instructor: Raktim Bhattacharya 35 / 38
36 Tracking with disturbance Longitudinal F16 Control Law Tracking Performance V ref =, γ ref = 45 deg step V (ft/s) Gamma (deg) Thrust (lb) Elevator (deg) Time Time AERO 632, Instructor: Raktim Bhattacharya 36 / 38
37 Tracking with disturbance Longitudinal F16 Control Law Tracking Performance V ref = 5 ft/s step, γ ref = 45 deg step V (ft/s) Gamma (deg) AERO 632, Instructor: Raktim Bhattacharya 37 / 38
38 Tracking with disturbance Longitudinal F16 Control Law Simulink Cy* u F16 Output V, gamma ysim To Workspace Vref Br* u Gain1 Gref x' = Ax+Bu y = Cx+Du F16 Longi h2k* u Gain Uniform Random Number AERO 632, Instructor: Raktim Bhattacharya 38 / 38
Raktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Norms for Signals and Systems
. AERO 632: Design of Advance Flight Control System Norms for Signals and. Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. Norms for Signals ...
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 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 informationRaktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Preliminaries
. AERO 632: of Advance Flight Control System. Preliminaries Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. Preliminaries Signals & Systems Laplace
More informationRaktim Bhattacharya. . AERO 422: Active Controls for Aerospace Vehicles. Dynamic Response
.. AERO 422: Active Controls for Aerospace Vehicles Dynamic Response Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. . Previous Class...........
More informationLecture 15: H Control Synthesis
c A. Shiriaev/L. Freidovich. March 12, 2010. Optimal Control for Linear Systems: Lecture 15 p. 1/14 Lecture 15: H Control Synthesis Example c A. Shiriaev/L. Freidovich. March 12, 2010. Optimal Control
More informationModern Optimal Control
Modern Optimal Control Matthew M. Peet Arizona State University Lecture 22: H 2, LQG and LGR Conclusion To solve the H -optimal state-feedback problem, we solve min γ such that γ,x 1,Y 1,A n,b n,c n,d
More informationPrinciples of Optimal Control Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 16.323 Principles of Optimal Control Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 16.323 Lecture
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 informationRaktim Bhattacharya. . AERO 422: Active Controls for Aerospace Vehicles. Frequency Response-Design Method
.. AERO 422: Active Controls for Aerospace Vehicles Frequency Response- Method Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. ... Response to
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 informationLecture 10: Linear Matrix Inequalities Dr.-Ing. Sudchai Boonto
Dr-Ing Sudchai Boonto Department of Control System and Instrumentation Engineering King Mongkuts Unniversity of Technology Thonburi Thailand Linear Matrix Inequalities A linear matrix inequality (LMI)
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 informationLecture 7 : Generalized Plant and LFT form Dr.-Ing. Sudchai Boonto Assistant Professor
Dr.-Ing. Sudchai Boonto Assistant Professor Department of Control System and Instrumentation Engineering King Mongkuts Unniversity of Technology Thonburi Thailand Linear Quadratic Gaussian The state space
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 Matrix Inequality (LMI)
Linear Matrix Inequality (LMI) A linear matrix inequality is an expression of the form where F (x) F 0 + x 1 F 1 + + x m F m > 0 (1) x = (x 1,, x m ) R m, F 0,, F m are real symmetric matrices, and the
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 informationRaktim Bhattacharya. . AERO 422: Active Controls for Aerospace Vehicles. Basic Feedback Analysis & Design
AERO 422: Active Controls for Aerospace Vehicles Basic Feedback Analysis & Design Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University Routh s Stability
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 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 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 informationRobust and Optimal Control, Spring 2015
Robust and Optimal Control, Spring 2015 Instructor: Prof. Masayuki Fujita (S5-303B) D. Linear Matrix Inequality D.1 Convex Optimization D.2 Linear Matrix Inequality(LMI) D.3 Control Design and LMI Formulation
More informationRiccati Equations and Inequalities in Robust Control
Riccati Equations and Inequalities in Robust Control Lianhao Yin Gabriel Ingesson Martin Karlsson Optimal Control LP4 2014 June 10, 2014 Lianhao Yin Gabriel Ingesson Martin Karlsson (LTH) H control problem
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 informationNonlinear Control Systems
Nonlinear Control Systems António Pedro Aguiar pedro@isr.ist.utl.pt 5. Input-Output Stability DEEC PhD Course http://users.isr.ist.utl.pt/%7epedro/ncs2012/ 2012 1 Input-Output Stability y = Hu H denotes
More informationProblem Set 5 Solutions 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Problem Set 5 Solutions The problem set deals with Hankel
More 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 informationChapter 3. LQ, LQG and Control System Design. Dutch Institute of Systems and Control
Chapter 3 LQ, LQG and Control System H 2 Design Overview LQ optimization state feedback LQG optimization output feedback H 2 optimization non-stochastic version of LQG Application to feedback system design
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 informationTopic # Feedback Control Systems
Topic #20 16.31 Feedback Control Systems Closed-loop system analysis Bounded Gain Theorem Robust Stability Fall 2007 16.31 20 1 SISO Performance Objectives Basic setup: d i d o r u y G c (s) G(s) n control
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 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design.
ISS0031 Modeling and Identification Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design. Aleksei Tepljakov, Ph.D. September 30, 2015 Linear Dynamic Systems Definition
More informationAn Introduction to Linear Matrix Inequalities. Raktim Bhattacharya Aerospace Engineering, Texas A&M University
An Introduction to Linear Matrix Inequalities Raktim Bhattacharya Aerospace Engineering, Texas A&M University Linear Matrix Inequalities What are they? Inequalities involving matrix variables Matrix variables
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 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 8. Chapter 5: Input-Output Stability Chapter 6: Passivity Chapter 14: Passivity-Based Control. Eugenio Schuster.
Lecture 8 Chapter 5: Input-Output Stability Chapter 6: Passivity Chapter 14: Passivity-Based Control Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture
More informationProblem Set 2 Solutions 1
Massachusetts Institute of echnology Department of Electrical Engineering and Computer Science 6.45: MULIVARIABLE CONROL SYSEMS by A. Megretski Problem Set Solutions Problem. For each of the statements
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 informationObservability. It was the property in Lyapunov stability which allowed us to resolve that
Observability We have seen observability twice already It was the property which permitted us to retrieve the initial state from the initial data {u(0),y(0),u(1),y(1),...,u(n 1),y(n 1)} It was the property
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 informationThe norms can also be characterized in terms of Riccati inequalities.
9 Analysis of stability and H norms Consider the causal, linear, time-invariant system ẋ(t = Ax(t + Bu(t y(t = Cx(t Denote the transfer function G(s := C (si A 1 B. Theorem 85 The following statements
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 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 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 informationGrammians. Matthew M. Peet. Lecture 20: Grammians. Illinois Institute of Technology
Grammians Matthew M. Peet Illinois Institute of Technology Lecture 2: Grammians Lyapunov Equations Proposition 1. Suppose A is Hurwitz and Q is a square matrix. Then X = e AT s Qe As ds is the unique solution
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 informationStationary trajectories, singular Hamiltonian systems and ill-posed Interconnection
Stationary trajectories, singular Hamiltonian systems and ill-posed Interconnection S.C. Jugade, Debasattam Pal, Rachel K. Kalaimani and Madhu N. Belur Department of Electrical Engineering Indian Institute
More informationTopic # Feedback Control
Topic #7 16.31 Feedback Control State-Space Systems What are state-space models? Why should we use them? How are they related to the transfer functions used in classical control design and how do we develop
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 information14 - Gaussian Stochastic Processes
14-1 Gaussian Stochastic Processes S. Lall, Stanford 211.2.24.1 14 - Gaussian Stochastic Processes Linear systems driven by IID noise Evolution of mean and covariance Example: mass-spring system Steady-state
More informationẋ n = f n (x 1,...,x n,u 1,...,u m ) (5) y 1 = g 1 (x 1,...,x n,u 1,...,u m ) (6) y p = g p (x 1,...,x n,u 1,...,u m ) (7)
EEE582 Topical Outline A.A. Rodriguez Fall 2007 GWC 352, 965-3712 The following represents a detailed topical outline of the course. It attempts to highlight most of the key concepts to be covered and
More informationIndustrial Model Predictive Control
Industrial Model Predictive Control Emil Schultz Christensen Kongens Lyngby 2013 DTU Compute-M.Sc.-2013-49 Technical University of Denmark DTU Compute Matematiktovet, Building 303B, DK-2800 Kongens Lyngby,
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 informationCore Concepts Review. Orthogonality of Complex Sinusoids Consider two (possibly non-harmonic) complex sinusoids
Overview of Continuous-Time Fourier Transform Topics Definition Compare & contrast with Laplace transform Conditions for existence Relationship to LTI systems Examples Ideal lowpass filters Relationship
More information9 The LQR Problem Revisited
9 he LQR Problem Revisited Problem: Compute a state feedback controller u(t) = Kx(t) that stabilizes the closed loop system and minimizes J := lim E [ z(t) z(t) ] t for the LI system Assumptions: a) D
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 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 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 informationL2 gains and system approximation quality 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 24: MODEL REDUCTION L2 gains and system approximation quality 1 This lecture discusses the utility
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 informationJoão P. Hespanha. January 16, 2009
LINEAR SYSTEMS THEORY João P. Hespanha January 16, 2009 Disclaimer: This is a draft and probably contains a few typos. Comments and information about typos are welcome. Please contact the author at hespanha@ece.ucsb.edu.
More informationTopics in control Tracking and regulation A. Astolfi
Topics in control Tracking and regulation A. Astolfi Contents 1 Introduction 1 2 The full information regulator problem 3 3 The FBI equations 5 4 The error feedback regulator problem 5 5 The internal model
More informationto have roots with negative real parts, the necessary and sufficient conditions are that:
THE UNIVERSITY OF TEXAS AT SAN ANTONIO EE 543 LINEAR SYSTEMS AND CONTROL H O M E W O R K # 7 Sebastian A. Nugroho November 6, 7 Due date of the homework is: Sunday, November 6th @ :59pm.. The following
More informationChapter 9 Observers, Model-based Controllers 9. Introduction In here we deal with the general case where only a subset of the states, or linear combin
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A. Dahleh George Verghese Department of Electrical Engineering and Computer Science Massachuasetts Institute of Technology c Chapter 9 Observers,
More informationHankel Optimal Model Reduction 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 2004: MODEL REDUCTION Hankel Optimal Model Reduction 1 This lecture covers both the theory and
More information2006 Fall. G(s) y = Cx + Du
1 Class Handout: Chapter 7 Frequency Domain Analysis of Feedback Systems 2006 Fall Frequency domain analysis of a dynamic system is very useful because it provides much physical insight, has graphical
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 informationContents lecture 5. Automatic Control III. Summary of lecture 4 (II/II) Summary of lecture 4 (I/II) u y F r. Lecture 5 H 2 and H loop shaping
Contents lecture 5 Automatic Control III Lecture 5 H 2 and H loop shaping Thomas Schön Division of Systems and Control Department of Information Technology Uppsala University. Email: thomas.schon@it.uu.se,
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 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 informationLecture 3. Chapter 4: Elements of Linear System Theory. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.
Lecture 3 Chapter 4: Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 3 p. 1/77 3.1 System Descriptions [4.1] Let f(u) be a liner operator, u 1 and u
More 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 informationIntroduction. Performance and Robustness (Chapter 1) Advanced Control Systems Spring / 31
Introduction Classical Control Robust Control u(t) y(t) G u(t) G + y(t) G : nominal model G = G + : plant uncertainty Uncertainty sources : Structured : parametric uncertainty, multimodel uncertainty Unstructured
More informationAnalysis of systems containing nonlinear and uncertain components by using Integral Quadratic Constraints
Analysis of systems containing nonlinear and uncertain components by using Integral Quadratic Constraints Tamas Peni MTA-SZTAKI October 2, 2012 Tamas Peni (MTA-SZTAKI) Seminar October 2, 2012 1 / 42 Contents
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 informationPole placement control: state space and polynomial approaches Lecture 2
: state space and polynomial approaches Lecture 2 : a state O. Sename 1 1 Gipsa-lab, CNRS-INPG, FRANCE Olivier.Sename@gipsa-lab.fr www.gipsa-lab.fr/ o.sename -based November 21, 2017 Outline : a state
More informationTime Response Analysis (Part II)
Time Response Analysis (Part II). A critically damped, continuous-time, second order system, when sampled, will have (in Z domain) (a) A simple pole (b) Double pole on real axis (c) Double pole on imaginary
More information1 Controllability and Observability
1 Controllability and Observability 1.1 Linear Time-Invariant (LTI) Systems State-space: Dimensions: Notation Transfer function: ẋ = Ax+Bu, x() = x, y = Cx+Du. x R n, u R m, y R p. Note that H(s) is always
More informationRobust Control 9 Design Summary & Examples
Robust Control 9 Design Summary & Examples Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 5/6/003 Outline he H Problem Solution via γ-iteration Robust stability via
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 informationEvery real system has uncertainties, which include system parametric uncertainties, unmodeled dynamics
Sensitivity Analysis of Disturbance Accommodating Control with Kalman Filter Estimation Jemin George and John L. Crassidis University at Buffalo, State University of New York, Amherst, NY, 14-44 The design
More informationH 2 Suboptimal Estimation and Control for Nonnegative
Proceedings of the 2007 American Control Conference Marriott Marquis Hotel at Times Square New York City, USA, July 11-13, 2007 FrC20.3 H 2 Suboptimal Estimation and Control for Nonnegative Dynamical Systems
More informationProblem Set 4 Solutions 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6245: MULTIVARIABLE CONTROL SYSTEMS by A Megretski Problem Set 4 Solutions Problem 4T For the following transfer
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 informationModern Optimal Control
Modern Optimal Control Matthew M. Peet Arizona State University Lecture 19: Stabilization via LMIs Optimization Optimization can be posed in functional form: min x F objective function : inequality constraints
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 informationNonlinear Control. Nonlinear Control Lecture # 18 Stability of Feedback Systems
Nonlinear Control Lecture # 18 Stability of Feedback Systems Absolute Stability + r = 0 u y G(s) ψ( ) Definition 7.1 The system is absolutely stable if the origin is globally uniformly asymptotically stable
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 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 informationNonlinear Control. Nonlinear Control Lecture # 6 Passivity and Input-Output Stability
Nonlinear Control Lecture # 6 Passivity and Input-Output Stability Passivity: Memoryless Functions y y y u u u (a) (b) (c) Passive Passive Not passive y = h(t,u), h [0, ] Vector case: y = h(t,u), h T =
More informationÜbersetzungshilfe / Translation aid (English) To be returned at the end of the exam!
Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 3.. 24 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid -
More informationLMIs for Observability and Observer Design
LMIs for Observability and Observer Design Matthew M. Peet Arizona State University Lecture 06: LMIs for Observability and Observer Design Observability Consider a system with no input: ẋ(t) = Ax(t), x(0)
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 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 informationNonlinear Observers. Jaime A. Moreno. Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México
Nonlinear Observers Jaime A. Moreno JMorenoP@ii.unam.mx Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México XVI Congreso Latinoamericano de Control Automático October
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 informationChapter 7 Interconnected Systems and Feedback: Well-Posedness, Stability, and Performance 7. Introduction Feedback control is a powerful approach to o
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A. Dahleh George Verghese Department of Electrical Engineering and Computer Science Massachuasetts Institute of Technology c Chapter 7 Interconnected
More informationLecture 8. Applications
Lecture 8. Applications Ivan Papusha CDS270 2: Mathematical Methods in Control and System Engineering May 8, 205 / 3 Logistics hw7 due this Wed, May 20 do an easy problem or CYOA hw8 (design problem) will
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 information