NEW DEVELOPMENTS IN PREDICTIVE CONTROL FOR NONLINEAR SYSTEMS
|
|
- Hilary Owen
- 5 years ago
- Views:
Transcription
1 NEW DEVELOPMENTS IN PREDICTIVE CONTROL FOR NONLINEAR SYSTEMS M. J. Grimble, A. Ordys, A. Dutka, P. Majecki University of Strathclyde Glasgow Scotland, U.K
2 Introduction Model Predictive Control (MPC) is one of the most popular advanced control techniques The MPC algorithms are well established for linear systems Recent developments extended this methodology to the Non-linear systems control Techniques developed at the University of Strathclyde are presented
3 Linear Quadratic Gaussian Predictive Control LQGPC
4 Preview Control The algorithm was first presented by Tomizuka M. and D.E. Whitney. The algorithm uses the LQG approach to optimisation and a stochastic model for the reference signal beyond the preview horizon. State space model: x( t ) = Ax( t ) Bu( t ) Gw( t ) y( t) = Cx( t) v( t) Performance index: R S M L T J( t) = Ε y( t j ) r( t j ) y( t j ) r( t j ) T j= NM b g b g T λ ut ( j) ut ( j) } Reference generator: R ( t ) = Θ R ( t) η ξ( t N ) N N N N
5 Preview Control L NM x x( t ) ( t ) O QP = =L N M A Ο B G Ο w( t) χ( t ) χ( t ) ut ( ) Ο Θ Ο Ο µ ξ( t N ) RN, N N L NM R y( t) N ( t) O L QP = N M C Ο C Ο Ι OL QP NM x x( t) RN, A ( t) O QP O L QP N M Ι Ο O QP L N M M v( t) O QP B L N M M G OL QP NM O QP A standard (infinite horizon) LQG approach can now be used.
6 Linear Quadratic Gaussian Predictive Control LQGPC The standard predictive control performance index: J = Y( t) R( t) Λ Y( t) R( t) U( t) Λ U( t) t b g b g b g b g Starting with the state-space model: T e T u re-formulate it to be compatible with the above index: state equation: x( t ) = A x( t ) β U( t ) Γ W ( t ) output equation: ~ Y( t) = Φ Ax( t) S U( t) S W( t) V( t) N N N x( t ) = Ax( t ) Bu( t ) Gw( t ) y( t ) = Dx( t ) v( t ) The control vector U(t) contains all control actions within the horizon N into the future
7 Linear Quadratic Gaussian Predictive Control LQGPC Define the LQG-type performance index (finite or infinite) as a sum of predictive performance indices: J DPC J t = Ε R S T R S T T t T l= t T l= t J GPC U V W ( l) J DPC = Ε R S T t Τ lim JGPC ( l) T Τ l= t Substituting the criterion from the previous slide obtains: L O NM b g b g b g b gqp T T = Ε Y( l) R( l) Λ e Y( l) R( l) U( l) Λ u U( l) T The solution can be obtained through Dynamic Programming with two Riccati equations involved. The reference generator is used in a similar way as in Preview Control U V W U V W
8 Solution of Quadratic Gaussian problem for nonlinear systems and Non-Linear Quadratic Gaussian Predictive Control
9 NLQG and NLQGPC NLQG - an extension of SDRE method SDRE method State Dependent Riccati Equations The NLQGPC algorithm: predictive extension of SDRE
10 System representation The model: xt = f ( xt) g( xt) ut Gξt y = h( x ) t t Linear State Dependent form of the model: xt = A( xt) xt B( xt) ut Gξt y = C( x ) x t t t Assumption on the parameterisation of the model: xu, ( ( ), ( )) A x B x is controllable t t
11 System representation Simplified notation: A = Ax ( ), B = Bx ( ), C = Cx ( ) t t t t t t The Prediction of the trajectory: u, u,..., u x, x,..., x t t t N t t 2 t N Assumption on the trajectory after the prediction horizon.
12 The NLQG algorithm Estimate (or measure) the state x(t) Use previous feedback gain K(t-) to calculate prediction of current control u(t) ut () = Kt ( ) xt () Use current control prediction u(t) and the model re-calculated at time instant t (with the state x(t)) to obtain future state prediction x(t).
13 The NLQG algorithm The state prediction x(t) together with the state feedback gain K(t-) from previous iteration of the algorithm is used for a calculation of the future control prediction u(t). ut () = Kt ( ) xt () The model once again is re-calculated using future state prediction, stored and sequence is repeated n times.
14 The NLQG algorithm Use the model prediction for time instant tn and solve Algebraic Riccati Equation. The solution at time instant tn is obtained: Pt ( n, )
15 The NLQG algorithm Use as a boundary condition Pt ( n) = Pt ( n, ) for iterations of the Riccati Difference Equation and use appropriate prediction of the model throughout iterations of Riccati Equation. Use Pk ( ) to calculate the feedback control gain and calculate the current control.
16 The NLQG algorithm Use Pk ( ) to calculate the feedback control gain and calculate the current control. ut () = Ktxt () () K() t = function( A(), t B(), t P( t )) Calculated current control is used for the plant input signal manipulation.
17 The NLQG cost function The following cost function is minimised: J t T T t Th x ( k) Qc( k) x( k) u ( k) Rcu( k) = E lim Th 2 T k t T T T h = u ( k) Mc x( k) x ( k) Mcu( k) The cost function may be split in two parts: J = E lim J J Th 2 Th ( finite infinite ) t t t
18 The NLQG cost function First part is an infinite cost function: J infinite t T T t Th x ( k) Qc( t np) x( k) u ( k) Rcu( k) = k= t n T T T u ( k) Mc x( k) x ( k) Mcu( k) and this part is minimised by Algebraic Riccati Equation. Second part is a finite cost function J finite t T T t n x ( k) Qc( k) x( k) u ( k) Rcu( k) = k= t T T T u ( k) Mc x( k) x ( k) Mcu( k) And is minimised by Difference Riccati Equation with border condition given by the solution of ARE.
19 Remarks To obtain accurate results it should be assumed that system will remain time invariant after tn time instant: Therefore if real behaviour of the system is closer to the assumption results are more accurate
20 Example
21 Example Plant model Reference model ( x ) p,2 t atan ( ).7. xp( t ) = x ( ) ( ) ( ) p,2() t xp t u t ξ p t.3 [ ] y () t = y() t = x () t h x ( t ) = [] x ( t) [] ξ ( t) r r r r () t = [] x () t h r p
22 Results output SDRE Proposed algorithm Constant gain feedback
23 Results 2.5 control SDRE Proposed algorithm Constant gain feedback
24 The NLQGPC control law derivation Re-written state equation: where xt = Atxt βtut, N Gξt β = [ B,,,..., ] t t 2 N State-space model with prediction-output equation x = A x β U Gξ t t t t t, N t Y =Φ A x S U G Ξ t, N tn, t t tn, tn, tn, tn, Note: that state equation is identical to the state equation of the controlled system.
25 The NLQGPC control law derivation Reference signal model: R R R R R t = t ξt X A X G R R R t, N = C Xt Augmented system : A t = Θ t t Ω tut, N Γ t χ χ ξ with Ψ = ϒ χ S U G Ξ t, N t t tn, tn, tn, tn, xt At A ξt βt G χt =, t, ξt, t,, R Θ = R = R Ω = Γ = R X t A ξ t G Ψ t, N = Yt, N Rt, N, ϒ t = Φt, NAt C R
26 The NLQGPC control law derivation The cost function: t Th N N J = E lim ( y r ) Λ ( y r ) u Λ u or: Introduce notation: Final form: { T i } { T i } t k i k i E k i k i k i U k i Th 2Th k= t i= i= T T T T χk ϒk ΛEϒk χk Uk, N Sk, N ΛEϒk χk t T h T T T T Jt= E lim χk k ESkN, UkN, UkN, SkN, ESkN, U ϒ Λ Λ kn, Th 2 J T h k= t T UkN, ΛU U kn, T T T k =ϒk ΛEϒ k, k =ϒk Λ E k, N, k = k, N Λ E k, N ΛU Q M S R S S Q M J E U J t Th T T k k χk t = lim k k, N Th 2 χ T T h U k= t Mk Rk kn,
27 NLQGPC control law derivation - Algorithm The control law minimising the cost function: ( T ) ( T T = Ω ) Ω Ω Θ P U P R P M χ where t, : solution of Algebraic Riccati Equation t, N t t, t t t t, t t t ( )( ) ( T T ) T T T Pt, = Qt Θt Pt, Θt Mt Θt Pt, Ω t Rt Ωt Pt, Ω t Mt Ωt Pt, Θ t This Algebraic Riccati Equation contains state dependent matrices Qt, Mt, Rt calculated at time t, which contain the prediction of future system behaviour
28 NLQGPC control law derivation - Algorithm 2 A more accurate solution of the minimisation problem: t N T T Qk Mk χk χk U k, N T U k t Mk R k kn, = Jt = E lim J Th 2 Th t N T h Q T T k Mk χk χk U k, N T U k= t N Mk Rk kn, Cost function split in two parts finite horizon infinite horizon.
29 NLQGPC control law derivation - Algorithm 2 Difference Riccati Equation : T T T k = k Θk k Θk k Θk k Ω k k Ωk k Ω k k Ωk k Θk ( )( ) ( T T ) P Q P M P R P M P boundary condition iterated backwards P = P t N t N, for k = t N, t N,..., t The control vector minimising cost function: ( T ) ( T T ) U = Ω P Ω R Ω P Θ M χ t, N t t t t t t t t t
30 Example The model is given by the following non-linear state space equations y ( ) ( ) 2 t = t.3 sin t t g t ζ ζ ζ ζ ξ t = t.3 t t g t ζ ζ ζ υ ξ t = ζ t Next the model is re-arranged into Linear State Dependent form: [ ] ζ [ ] ( ζ t ).3 sin g ζ ζ υ ξ y t = ζ t t t t g 2 = υ, t t t.3 2 ( ζ t )
31 Example The system is controllable since rank ( ( 2) ) 2 ζ =.3 ζ t 2
32 Example The step response for two NLQGPC algorithms is compared with SDRE with g=. (noise level) 4 output Alg. 2, Setpoint=: J=2.578Setpoint=3: J= Alg., Setpoint=: J=2.524 Setpoint=3: J= SDRE, Setpoint=: J= Setpoint=3: J=
33 Example Now compare noise rejection (for two levels of process noise) 4 output 3 2 Alg. 2, g=.: J= Alg., g=.: J= SDRE, g=.: J= output 3 2 Alg. 2, g=.: J= Alg., g=.: J= SDRE, g=.: J=
34 Advantages & Disadvantages of NLQGPC Advantages:. Controls based on solutions to the NLQGPC have been shown to offer high performance. 2. Less computational burden than other non-linear predictive control techniques. Disadvantages:. Since NLQGPC utilizes the Riccati equation, it is an unconstrained predictive control technique. 2. Like SDRE, NLQGPC doesn t guarantee closedloop global stability.
35 Dealing with constraints in NLQGPC The input constraints can be approximated by means of smooth limiting functions, and then included into the dynamics of the plant in a state-dependent state-space form. u α α u α
36 Improving stability via Satisficing Satisficing is based on a point-wise cost/benefit comparison of an action. The benefits are given by the Selectability function P s (u,x), while the costs are given by the Rejectability function P r (u,x). The satisficing set is those options for which selectability exceeds rejectability: i.e., S( x,b ) = { u : P ( u,x ) bp ( u,x )} s r
37 CLF-Based Satisficing Technique The selectability criteria is defined to be: Ps ( u,x ) = V x ( f gu ) The rejectability criteria is defined to be: P ( ux, ) = lx ( ) xrx r b=, therefore, the satisficing set S: { T } x S( x, b) = u: V ( f gu) T T
38 Augmenting NLQGPC with Satisficing By projecting the NLQGPC controller point-wise onto the satisficing set, the good properties of the NLQGPC approach are combined with the analytical properties of satisficing.
39 Example: Control of F-8 aircraft The non-linear dynamical model of the F-8 fighter aircraft: rad. u u. u x u x. u. x. x. x. x. x, x x, u. u x. u x. u. x. x x x. x. x x. x x. x = = =
40 Angle of Attack (rad) Blue lines Unconstrained NLQGPC. Black lines Constrained NLQGPC. Magenta lines Constrained NLQGPC with guarantee of global asymptotic stability Time
41 .4 Elevator Deflection.2 u (rad) u rad Time
42 NLQGPC High Performance Deals with input Constraints Low Computational Burden Guarantee of Robustness & Asymptotic Stability
43 Non-Linear Generalized Minimum Variance Control
44 Contents Introduction Nonlinear GMV control problem and solution Relationship to the Smith Predictor Incorporating future information: Feedforward and Tracking Simulation example
45 Introduction 969: Åström introduces Minimum Variance (MV) controller assuming linear minimum phase plant. Successful applications in pulp and paper industry. 97s: Clarke and Hastings-James modify the MV control law by adding a control costing term. This is termed a Generalized Minimum Variance (GMV) control law and is the basis for their later self-tuning controller. The GMV control law has similar characteristics to LQG design in some cases and is much simpler to implement However, when the control weighting tends to zero the control law reverts to the initial algorithm of Åström, which is unstable for non-minimum phase processes.
46 Introduction Aim: introduce a GMV controller for nonlinear, multivariable, possibly time-varying processes The structure of the system is defined so that a simple solution is obtained. When the system is linear the results revert to those for the linear GMV controller. There is some loss of generality in assuming the reference and disturbance models are represented by linear subsystems. However, plant model can be in a very general nonlinear operator form, which might involve state-space, transfer operators, neural networks or even nonlinear function look-up tables.
47 Nonlinear system description φ = PeF u c c ξ Error weighting P c F c Control weighting W d Disturbance model ω Reference W r Controller r e C u - Nonlinear plant W m d y Nonlinear plant model: Linear disturbance model: Linear reference model: -k ( Wu)( t) = z ( W u)( t) d = f W A C r f d W = A E r k
48 Plant model Nonlinear plant model can be given in a very general form, e.g.: state-space formulation neural network / neuro-fuzzy model look-up table Fortran/C code It can include both linear and nonlinear components, e.g. Hammerstein model: ζ u y f( u, y ) = Plant subsystems W d d u m y W W k k Control Output Nonlinear Linear Delay z k Just need to obtain the output to given input signal
49 NGMV problem formulation To minimize: variance of the generalized output φ (t): J = E[ φ ( t)] NGMV 2 with ( ) ( )( ) φ t = Pe t F u t () c c P c = PcnP cd k ( u)() t = z ( u)( t) c - linear error weighting F F - control weighting (possibly nonlinear) ck Control weighting assumed invertible and potentially nonlinear to compensate for plant nonlinearities in appropriate cases The weighting selection is restricted by closed-loop stability
50 NGMV problem solution The approach also similar: ( W )( ) ε ( F )( ) k c k f c φ () t = P ( z u t Y ()) t u t k () = z ( F PW ) u t PY ε( t) ck c k c f ( ) φ () t = Fε() t ( Fck Pc Wk ) u t k Rε( t k) k PY c f = F z R Diophantine equation YY = WW WW * * * f f d d r r Spectral factorization statistically independent ε(t) white noise (sequence of independent random variables) Optimal control: NGMV u () t = ( F PW ) Rε () t ck c k stable causal nonlinear operator inverse
51 Controller implementation NGMV ck f Wk f u () t = [( F FY ) RY e]() t Disturbance d Controller Reference r - e - RY f Fck u Plant W Output y linear blocks FY f W k
52 Existence of a Stable Operator Inverse Necessary condition for optimality: operator ( PW F ) c k ck must have a stable inverse For linear systems: the operator must be strictly minimum-phase. To show this is satisfied for a very wide class of systems consider the case where F ck is linear and negative so that F ck = -F k. Then obtain: ( ) ( ) cwk k = k k cwk P F u F F P I u return-difference operator for a feedback system with A delay-free plant and controller c = k c. K F P
53 PID-based initial design Consider the delay-free plant W k and assume a PID controller K PID exists to stabilize the closed-loop system. Then a starting point for the weighting choice that will ensure the operator ( PW F ) c k k is stably invertible is P = K, F = c PID k To demonstrate this selection reasonable consider scalar case and let controller ( ) 2 ( k k k ) ( k 2k ) z k z k Kc = k k z = z z Assume the PID gains are positive numbers, with small derivative gain. Then simple to confirm if F k = the P cn term is minimum phase and has real zeros.
54 Relationship to the Smith Predictor The optimal controller can be expressed in a similar form to that of a Smith Predictor. This provides a new nonlinear version of the Smith Predictor. Compensator d Reference r _ ( AP ) p cd GY y f F ck u Plant W W k F Y f ( AP ) GY D p cd f k m k D k -
55 Smith Predictor form of NGMV controller The system may be redrawn and the compensator rearranged as shown below. This structure is essential if P c includes an integrator. Compensator Disturbance Reference r - p f A GY - ψ F ck P cd u Plant W Output y W k P cn D k p _
56 Comments This last structure is intuitively reasonable. With no plant-model mismatch, the control is not due to feedback but involves an open-loop stable compensator. The nonlinear inner-loop has weightings Fck Pc acting like an inner-loop controller. If weightings are chosen to be of usual form this will represent a filtered PID controller. Such a choice of weightings is only a starting point, since stability is easier to achieve. However, control weighting can have additional lead term and high frequency characteristics of optimal controller will then have more realistic roll off. Stability: Under the given assumptions the resulting Smith system is stable. This follows because the plant is stable, the inner-loop is stable and there are only stable terms in the input block.
57 Feedback Feedforward Tracking Control future reference information ( ) ( )( ) φ t = Pe t F u t () c c η W r Setpoint Error weighting Reference P c r r y f C 2 C Controller F c Control weighting Nonlinear plant Measurable ζ W d d Un-measurable ξ W d d Disturbance models ω W w w H f - e C u W m y Scaling H f Feedback gain/dynamics
58 Feedback, Tracking and Feedforward Control Signal Generation Control Modules Measured disturbance y f i ( ) d f d d di d3 GP D z GP D W W u Reference r Nonlinear ( t p) GP 2 rd E r r plant u - Setpoint Hf w GP D e d f Fck - W m u f Total disturbance d Output y f F Y W k Feedback, Feedforward and Tracking Controller H f
59 Simulated example 2-by-2 model given in the non-linear state space form: x () t x ( t ) = u ( t) 2 2 x ( t) x2 () t 2 = 2 2 x ( t ).9 x () t e u () t yt () = xt () Both outputs are followed by a transport delay of k = 6 samples so the timedelay matrix: Models: D k 6 z = 6 z W r z = z reference W d z 6 = z z measurable disturbance W d..5z =..5z unmeasurable disturbance
60 Transient responses Feedforward action 2.5 Output setpoint NGMV NGMVFF NGMVFFTR.5 Control Output Control 2 NGMV NGMVFF NGMVFFTR Future reference information incorporated
61 Stochastic performance 3 2 Output setpoint NGMV NGMVFF NGMVFFTR 2 Control Output Control 2 NGMV NGMVFF NGMVFFTR Controller Var[e] Var[u] Var[φ ] NGMV FB NGMV FBFF NGMV FBFFTR
62 Concluding Remarks State dependent models gives a useful structure for NL predictive controllers NGMV has much potential for development with multi - step predictive control an obvious development. Key to success in NL control is to show works and practical on real processes - why NGMV seems great potential.
Nonlinear Model Predictive Control Tools (NMPC Tools)
Nonlinear Model Predictive Control Tools (NMPC Tools) Rishi Amrit, James B. Rawlings April 5, 2008 1 Formulation We consider a control system composed of three parts([2]). Estimator Target calculator Regulator
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 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 informationProcess Modelling, Identification, and Control
Jan Mikles Miroslav Fikar 2008 AGI-Information Management Consultants May be used for personal purporses only or by libraries associated to dandelon.com network. Process Modelling, Identification, and
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 informationIterative Learning Control Analysis and Design I
Iterative Learning Control Analysis and Design I Electronics and Computer Science University of Southampton Southampton, SO17 1BJ, UK etar@ecs.soton.ac.uk http://www.ecs.soton.ac.uk/ Contents Basics Representations
More informationQuis custodiet ipsos custodes?
Quis custodiet ipsos custodes? James B. Rawlings, Megan Zagrobelny, Luo Ji Dept. of Chemical and Biological Engineering, Univ. of Wisconsin-Madison, WI, USA IFAC Conference on Nonlinear Model Predictive
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 informationModeling and Control Overview
Modeling and Control Overview D R. T A R E K A. T U T U N J I A D V A N C E D C O N T R O L S Y S T E M S M E C H A T R O N I C S E N G I N E E R I N G D E P A R T M E N T P H I L A D E L P H I A U N I
More informationCourse on Model Predictive Control Part II Linear MPC design
Course on Model Predictive Control Part II Linear MPC design Gabriele Pannocchia Department of Chemical Engineering, University of Pisa, Italy Email: g.pannocchia@diccism.unipi.it Facoltà di Ingegneria,
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More 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 informationOptimal Polynomial Control for Discrete-Time Systems
1 Optimal Polynomial Control for Discrete-Time Systems Prof Guy Beale Electrical and Computer Engineering Department George Mason University Fairfax, Virginia Correspondence concerning this paper should
More informationSteady State Kalman Filter
Steady State Kalman Filter Infinite Horizon LQ Control: ẋ = Ax + Bu R positive definite, Q = Q T 2Q 1 2. (A, B) stabilizable, (A, Q 1 2) detectable. Solve for the positive (semi-) definite P in the ARE:
More informationDesign Methods for Control Systems
Design Methods for Control Systems Maarten Steinbuch TU/e Gjerrit Meinsma UT Dutch Institute of Systems and Control Winter term 2002-2003 Schedule November 25 MSt December 2 MSt Homework # 1 December 9
More informationControl System Design
ELEC ENG 4CL4: Control System Design Notes for Lecture #36 Dr. Ian C. Bruce Room: CRL-229 Phone ext.: 26984 Email: ibruce@mail.ece.mcmaster.ca Friday, April 4, 2003 3. Cascade Control Next we turn to an
More informationEML5311 Lyapunov Stability & Robust Control Design
EML5311 Lyapunov Stability & Robust Control Design 1 Lyapunov Stability criterion In Robust control design of nonlinear uncertain systems, stability theory plays an important role in engineering systems.
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 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 information(RPG) (2017) IET,
Savvidis, Petros and Grimble, Michael and Majecki, Pawel and Pang, Yan (2017) Nonlinear predictive generalized minimum variance LPV control of wind turbines. In: 5th IET International Conference on Renewable
More informationMultivariable MRAC with State Feedback for Output Tracking
29 American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June 1-12, 29 WeA18.5 Multivariable MRAC with State Feedback for Output Tracking Jiaxing Guo, Yu Liu and Gang Tao Department
More informationPart II: Model Predictive Control
Part II: Model Predictive Control Manfred Morari Alberto Bemporad Francesco Borrelli Unconstrained Optimal Control Unconstrained infinite time optimal control V (x 0 )= inf U k=0 [ x k Qx k + u k Ru k]
More informationEC Control Engineering Quiz II IIT Madras
EC34 - Control Engineering Quiz II IIT Madras Linear algebra Find the eigenvalues and eigenvectors of A, A, A and A + 4I Find the eigenvalues and eigenvectors of the following matrices: (a) cos θ sin θ
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 informationPerformance assessment of MIMO systems under partial information
Performance assessment of MIMO systems under partial information H Xia P Majecki A Ordys M Grimble Abstract Minimum variance (MV) can characterize the most fundamental performance limitation of a system,
More informationCHAPTER 3 TUNING METHODS OF CONTROLLER
57 CHAPTER 3 TUNING METHODS OF CONTROLLER 3.1 INTRODUCTION This chapter deals with a simple method of designing PI and PID controllers for first order plus time delay with integrator systems (FOPTDI).
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 informationA FAST, EASILY TUNED, SISO, MODEL PREDICTIVE CONTROLLER. Gabriele Pannocchia,1 Nabil Laachi James B. Rawlings
A FAST, EASILY TUNED, SISO, MODEL PREDICTIVE CONTROLLER Gabriele Pannocchia, Nabil Laachi James B. Rawlings Department of Chemical Engineering Univ. of Pisa Via Diotisalvi 2, 5626 Pisa (Italy) Department
More informationFINITE HORIZON ROBUST MODEL PREDICTIVE CONTROL USING LINEAR MATRIX INEQUALITIES. Danlei Chu, Tongwen Chen, Horacio J. Marquez
FINITE HORIZON ROBUST MODEL PREDICTIVE CONTROL USING LINEAR MATRIX INEQUALITIES Danlei Chu Tongwen Chen Horacio J Marquez Department of Electrical and Computer Engineering University of Alberta Edmonton
More informationINVERSE MODEL APPROACH TO DISTURBANCE REJECTION AND DECOUPLING CONTROLLER DESIGN. Leonid Lyubchyk
CINVESTAV Department of Automatic Control November 3, 20 INVERSE MODEL APPROACH TO DISTURBANCE REJECTION AND DECOUPLING CONTROLLER DESIGN Leonid Lyubchyk National Technical University of Ukraine Kharkov
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 informationUCLA Chemical Engineering. Process & Control Systems Engineering Laboratory
Constrained Innite-time Optimal Control Donald J. Chmielewski Chemical Engineering Department University of California Los Angeles February 23, 2000 Stochastic Formulation - Min Max Formulation - UCLA
More informationModel predictive control of industrial processes. Vitali Vansovitš
Model predictive control of industrial processes Vitali Vansovitš Contents Industrial process (Iru Power Plant) Neural networ identification Process identification linear model Model predictive controller
More informationControl Systems I. Lecture 2: Modeling. Suggested Readings: Åström & Murray Ch. 2-3, Guzzella Ch Emilio Frazzoli
Control Systems I Lecture 2: Modeling Suggested Readings: Åström & Murray Ch. 2-3, Guzzella Ch. 2-3 Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich September 29, 2017 E. Frazzoli
More informationOptimal control and estimation
Automatic Control 2 Optimal control and estimation Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011
More informationRECURSIVE ESTIMATION AND KALMAN FILTERING
Chapter 3 RECURSIVE ESTIMATION AND KALMAN FILTERING 3. The Discrete Time Kalman Filter Consider the following estimation problem. Given the stochastic system with x k+ = Ax k + Gw k (3.) y k = Cx k + Hv
More informationChapter 2 Optimal Control Problem
Chapter 2 Optimal Control Problem Optimal control of any process can be achieved either in open or closed loop. In the following two chapters we concentrate mainly on the first class. The first chapter
More informationOptimal Control. Lecture 18. Hamilton-Jacobi-Bellman Equation, Cont. John T. Wen. March 29, Ref: Bryson & Ho Chapter 4.
Optimal Control Lecture 18 Hamilton-Jacobi-Bellman Equation, Cont. John T. Wen Ref: Bryson & Ho Chapter 4. March 29, 2004 Outline Hamilton-Jacobi-Bellman (HJB) Equation Iterative solution of HJB Equation
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 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 informationControl Systems Lab - SC4070 Control techniques
Control Systems Lab - SC4070 Control techniques Dr. Manuel Mazo Jr. Delft Center for Systems and Control (TU Delft) m.mazo@tudelft.nl Tel.:015-2788131 TU Delft, February 16, 2015 (slides modified from
More informationBasic Concepts in Data Reconciliation. Chapter 6: Steady-State Data Reconciliation with Model Uncertainties
Chapter 6: Steady-State Data with Model Uncertainties CHAPTER 6 Steady-State Data with Model Uncertainties 6.1 Models with Uncertainties In the previous chapters, the models employed in the DR were considered
More informationCALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems. CDS 110b
CALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems CDS 110b R. M. Murray Kalman Filters 14 January 2007 Reading: This set of lectures provides a brief introduction to Kalman filtering, following
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 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 informationAdaptive Nonlinear Control A Tutorial. Miroslav Krstić
Adaptive Nonlinear Control A Tutorial Miroslav Krstić University of California, San Diego Backstepping Tuning Functions Design Modular Design Output Feedback Extensions A Stochastic Example Applications
More informationFUZZY CONTROL OF NONLINEAR SYSTEMS WITH INPUT SATURATION USING MULTIPLE MODEL STRUCTURE. Min Zhang and Shousong Hu
ICIC Express Letters ICIC International c 2008 ISSN 1881-803X Volume 2, Number 2, June 2008 pp. 131 136 FUZZY CONTROL OF NONLINEAR SYSTEMS WITH INPUT SATURATION USING MULTIPLE MODEL STRUCTURE Min Zhang
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 informationAC&ST AUTOMATIC CONTROL AND SYSTEM THEORY SYSTEMS AND MODELS. Claudio Melchiorri
C. Melchiorri (DEI) Automatic Control & System Theory 1 AUTOMATIC CONTROL AND SYSTEM THEORY SYSTEMS AND MODELS Claudio Melchiorri Dipartimento di Ingegneria dell Energia Elettrica e dell Informazione (DEI)
More informationOptimization-Based Control
Optimization-Based Control Richard M. Murray Control and Dynamical Systems California Institute of Technology DRAFT v1.7a, 19 February 2008 c California Institute of Technology All rights reserved. This
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 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 informationLinear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard June 15, 2013
Linear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard June 15, 2013 Abstract As in optimal control theory, linear quadratic (LQ) differential games (DG) can be solved, even in high dimension,
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 informationLifted approach to ILC/Repetitive Control
Lifted approach to ILC/Repetitive Control Okko H. Bosgra Maarten Steinbuch TUD Delft Centre for Systems and Control TU/e Control System Technology Dutch Institute of Systems and Control DISC winter semester
More informationFall 線性系統 Linear Systems. Chapter 08 State Feedback & State Estimators (SISO) Feng-Li Lian. NTU-EE Sep07 Jan08
Fall 2007 線性系統 Linear Systems Chapter 08 State Feedback & State Estimators (SISO) Feng-Li Lian NTU-EE Sep07 Jan08 Materials used in these lecture notes are adopted from Linear System Theory & Design, 3rd.
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 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 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 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 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 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 informationIntegrator Backstepping using Barrier Functions for Systems with Multiple State Constraints
Integrator Backstepping using Barrier Functions for Systems with Multiple State Constraints Khoi Ngo Dep. Engineering, Australian National University, Australia Robert Mahony Dep. Engineering, Australian
More informationAnalysis and Synthesis of Single-Input Single-Output Control Systems
Lino Guzzella Analysis and Synthesis of Single-Input Single-Output Control Systems l+kja» \Uja>)W2(ja»\ um Contents 1 Definitions and Problem Formulations 1 1.1 Introduction 1 1.2 Definitions 1 1.2.1 Systems
More informationGoodwin, Graebe, Salgado, Prentice Hall Chapter 11. Chapter 11. Dealing with Constraints
Chapter 11 Dealing with Constraints Topics to be covered An ubiquitous problem in control is that all real actuators have limited authority. This implies that they are constrained in amplitude and/or rate
More informationRobust Internal Model Control for Impulse Elimination of Singular Systems
International Journal of Control Science and Engineering ; (): -7 DOI:.59/j.control.. Robust Internal Model Control for Impulse Elimination of Singular Systems M. M. Share Pasandand *, H. D. Taghirad Department
More informationOutline. Linear regulation and state estimation (LQR and LQE) Linear differential equations. Discrete time linear difference equations
Outline Linear regulation and state estimation (LQR and LQE) James B. Rawlings Department of Chemical and Biological Engineering 1 Linear Quadratic Regulator Constraints The Infinite Horizon LQ Problem
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 information10/8/2015. Control Design. Pole-placement by state-space methods. Process to be controlled. State controller
Pole-placement by state-space methods Control Design To be considered in controller design * Compensate the effect of load disturbances * Reduce the effect of measurement noise * Setpoint following (target
More informationWannabe-MPC for Large Systems Based on Multiple Iterative PI Controllers
Wannabe-MPC for Large Systems Based on Multiple Iterative PI Controllers Pasi Airikka, Mats Friman Metso Corp., Finland 17th Nordic Process Control Workshop Jan 26-27 2012 DTU Denmark Content Motivation
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 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 informationChapter 2. Classical Control System Design. Dutch Institute of Systems and Control
Chapter 2 Classical Control System Design Overview Ch. 2. 2. Classical control system design Introduction Introduction Steady-state Steady-state errors errors Type Type k k systems systems Integral Integral
More informationCALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems. CDS 110b
CALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems CDS 110b R. M. Murray Kalman Filters 25 January 2006 Reading: This set of lectures provides a brief introduction to Kalman filtering, following
More informationH-infinity Model Reference Controller Design for Magnetic Levitation System
H.I. Ali Control and Systems Engineering Department, University of Technology Baghdad, Iraq 6043@uotechnology.edu.iq H-infinity Model Reference Controller Design for Magnetic Levitation System Abstract-
More informationOn-off Control: Audio Applications
On-off Control: Audio Applications Graham C. Goodwin Day 4: Lecture 3 16th September 2004 International Summer School Grenoble, France 1 Background In this lecture we address the issue of control when
More informationControl of Electromechanical Systems
Control of Electromechanical Systems November 3, 27 Exercise Consider the feedback control scheme of the motor speed ω in Fig., where the torque actuation includes a time constant τ A =. s and a disturbance
More information1 (30 pts) Dominant Pole
EECS C8/ME C34 Fall Problem Set 9 Solutions (3 pts) Dominant Pole For the following transfer function: Y (s) U(s) = (s + )(s + ) a) Give state space description of the system in parallel form (ẋ = Ax +
More 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 informationANALYSIS AND IMPROVEMENT OF THE INFLUENCE OF MEASUREMENT NOISE ON MVC BASED CONTROLLER PERFORMANCE ASSESSMENT
International Journal of Innovative Computing, Information and Control ICIC International c 2018 ISSN 1349-4198 Volume 14, Number 2, April 2018 pp. 697 716 ANALYSIS AND IMPROVEMENT OF THE INFLUENCE OF
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 informationOutline. 1 Full information estimation. 2 Moving horizon estimation - zero prior weighting. 3 Moving horizon estimation - nonzero prior weighting
Outline Moving Horizon Estimation MHE James B. Rawlings Department of Chemical and Biological Engineering University of Wisconsin Madison SADCO Summer School and Workshop on Optimal and Model Predictive
More informationNonlinear Control Systems
Nonlinear Control Systems António Pedro Aguiar pedro@isr.ist.utl.pt 7. Feedback Linearization IST-DEEC PhD Course http://users.isr.ist.utl.pt/%7epedro/ncs1/ 1 1 Feedback Linearization Given a nonlinear
More informationCDS 101/110a: Lecture 8-1 Frequency Domain Design
CDS 11/11a: Lecture 8-1 Frequency Domain Design Richard M. Murray 17 November 28 Goals: Describe canonical control design problem and standard performance measures Show how to use loop shaping to achieve
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 informationA unified double-loop multi-scale control strategy for NMP integrating-unstable systems
Home Search Collections Journals About Contact us My IOPscience A unified double-loop multi-scale control strategy for NMP integrating-unstable systems This content has been downloaded from IOPscience.
More informationLinear Systems. Manfred Morari Melanie Zeilinger. Institut für Automatik, ETH Zürich Institute for Dynamic Systems and Control, ETH Zürich
Linear Systems Manfred Morari Melanie Zeilinger Institut für Automatik, ETH Zürich Institute for Dynamic Systems and Control, ETH Zürich Spring Semester 2016 Linear Systems M. Morari, M. Zeilinger - Spring
More informationAERT 2013 [CA'NTI 19] ALGORITHMES DE COMMANDE NUMÉRIQUE OPTIMALE DES TURBINES ÉOLIENNES
AER 2013 [CA'NI 19] ALGORIHMES DE COMMANDE NUMÉRIQUE OPIMALE DES URBINES ÉOLIENNES Eng. Raluca MAEESCU Dr.Eng Andreea PINEA Prof.Dr.Eng. Nikolai CHRISOV Prof.Dr.Eng. Dan SEFANOIU Eng. Raluca MAEESCU CONEN
More informationEECE Adaptive Control
EECE 574 - Adaptive Control Recursive Identification in Closed-Loop and Adaptive Control Guy Dumont Department of Electrical and Computer Engineering University of British Columbia January 2010 Guy Dumont
More informationCALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems
CDS 101 1. Åström and Murray, Exercise 1.3 2. Åström and Murray, Exercise 1.4 3. Åström and Murray, Exercise 2.6, parts (a) and (b) CDS 110a 1. Åström and Murray, Exercise 1.4 2. Åström and Murray, Exercise
More informationAutomatic Control II Computer exercise 3. LQG Design
Uppsala University Information Technology Systems and Control HN,FS,KN 2000-10 Last revised by HR August 16, 2017 Automatic Control II Computer exercise 3 LQG Design Preparations: Read Chapters 5 and 9
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 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 informationHigh-Gain Observers in Nonlinear Feedback Control. Lecture # 3 Regulation
High-Gain Observers in Nonlinear Feedback Control Lecture # 3 Regulation High-Gain ObserversinNonlinear Feedback ControlLecture # 3Regulation p. 1/5 Internal Model Principle d r Servo- Stabilizing u y
More informationFeedback Control of Turbulent Wall Flows
Feedback Control of Turbulent Wall Flows Dipartimento di Ingegneria Aerospaziale Politecnico di Milano Outline Introduction Standard approach Wiener-Hopf approach Conclusions Drag reduction A model problem:
More informationOn Stochastic Adaptive Control & its Applications. Bozenna Pasik-Duncan University of Kansas, USA
On Stochastic Adaptive Control & its Applications Bozenna Pasik-Duncan University of Kansas, USA ASEAS Workshop, AFOSR, 23-24 March, 2009 1. Motivation: Work in the 1970's 2. Adaptive Control of Continuous
More informationJUSTIFICATION OF INPUT AND OUTPUT CONSTRAINTS INCORPORATION INTO PREDICTIVE CONTROL DESIGN
JUSTIFICATION OF INPUT AND OUTPUT CONSTRAINTS INCORPORATION INTO PREDICTIVE CONTROL DESIGN J. Škultéty, E. Miklovičová, M. Mrosko Slovak University of Technology, Faculty of Electrical Engineering and
More informationSystem Identification: From Data to Model
1 : From Data to Model With Applications to Aircraft Modeling Division of Automatic Control Linköping University Sweden Dynamic Systems 2 A Dynamic system has an output response y that depends on (all)
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 informationNMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582
NMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING NMT EE 589 & UNM ME 482/582 Simplified drive train model of a robot joint Inertia seen by the motor Link k 1 I I D ( q) k mk 2 kk Gk Torque amplification G
More information