Dynamic Inversion Design II

Similar documents
Gain Scheduling and Dynamic Inversion

Lecture 5 Classical Control Overview III. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore

Classical Numerical Methods to Solve Optimal Control Problems

Lecture 6 Classical Control Overview IV. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore

Using Lyapunov Theory I

Neuro-Adaptive Design II:

Linear Quadratic Regulator (LQR) Design I

Control Systems I. Lecture 4: Diagonalization, Modal Analysis, Intro to Feedback. Readings: Emilio Frazzoli

Linear Quadratic Regulator (LQR) Design II

Chapter 2 Review of Linear and Nonlinear Controller Designs

Lecture 4 Classical Control Overview II. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore

Adaptive Dynamic Inversion Control of a Linear Scalar Plant with Constrained Control Inputs

Constrained Optimal Control. Constrained Optimal Control II

Linear Quadratic Regulator (LQR) I

Constrained Optimal Control I

Multi-Robotic Systems

Nonlinear Control Lecture 9: Feedback Linearization

AUTOMATIC CONTROL. Andrea M. Zanchettin, PhD Spring Semester, Introduction to Automatic Control & Linear systems (time domain)

EE C128 / ME C134 Feedback Control Systems

Operator based robust right coprime factorization and control of nonlinear systems

Control of Mobile Robots

Neuro-Adaptive Design - I:

Gradient-Based Learning. Sargur N. Srihari

Ascent Phase Trajectory Optimization for a Hypersonic Vehicle Using Nonlinear Programming

1 Loop Control. 1.1 Open-loop. ISS0065 Control Instrumentation

Linear Quadratic Regulator (LQR) II

Chapter 13 Digital Control

Lecture 1: Feedback Control Loop

Exam. 135 minutes, 15 minutes reading time

Advanced Aerospace Control. Marco Lovera Dipartimento di Scienze e Tecnologie Aerospaziali, Politecnico di Milano

CDS 101/110a: Lecture 10-1 Robust Performance

Decentralized Control of Nonlinear Multi-Agent Systems Using Single Network Adaptive Critics

Slovak University of Technology in Bratislava Institute of Information Engineering, Automation, and Mathematics PROCEEDINGS

Topic # Feedback Control Systems

Control System Design

EECE 460. Decentralized Control of MIMO Systems. Guy A. Dumont. Department of Electrical and Computer Engineering University of British Columbia

Nonlinear Optimization for Optimal Control

Lecture 15 Review of Matrix Theory III. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore

Optimal Control Design

MCE693/793: Analysis and Control of Nonlinear Systems

Control Systems I. Lecture 7: Feedback and the Root Locus method. Readings: Jacopo Tani. Institute for Dynamic Systems and Control D-MAVT ETH Zürich

Ordinary Differential Equations Prof. A. K. Nandakumaran Department of Mathematics Indian Institute of Science Bangalore

Autonomous Helicopter Landing A Nonlinear Output Regulation Perspective

Autonomous navigation of unicycle robots using MPC

EECE Adaptive Control

EN Nonlinear Control and Planning in Robotics Lecture 10: Lyapunov Redesign and Robust Backstepping April 6, 2015

Switched Mode Power Conversion Prof. L. Umanand Department of Electronics Systems Engineering Indian Institute of Science, Bangalore

Control Systems I. Lecture 2: Modeling. Suggested Readings: Åström & Murray Ch. 2-3, Guzzella Ch Emilio Frazzoli

Chapter 1 Mathematical Preliminaries and Error Analysis

NEW DEVELOPMENTS IN PREDICTIVE CONTROL FOR NONLINEAR SYSTEMS

CS 450 Numerical Analysis. Chapter 5: Nonlinear Equations

Dependence. Practitioner Course: Portfolio Optimization. John Dodson. September 10, Dependence. John Dodson. Outline.

EE Control Systems LECTURE 9

A sub-optimal second order sliding mode controller for systems with saturating actuators

Riccati difference equations to non linear extended Kalman filter constraints

A Discrete Robust Adaptive Iterative Learning Control for a Class of Nonlinear Systems with Unknown Control Direction

Pattern Recognition Prof. P. S. Sastry Department of Electronics and Communication Engineering Indian Institute of Science, Bangalore

Uncertainty and Robustness for SISO Systems

Operation and Supply Chain Management Prof. G. Srinivasan Department of Management Studies Indian Institute of Technology, Madras

Control Systems I. Lecture 6: Poles and Zeros. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich

CS261: A Second Course in Algorithms Lecture #11: Online Learning and the Multiplicative Weights Algorithm

Queens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 2018 Instructor: Dr. Sateesh Mane.

Nonlinear Single-Particle Dynamics in High Energy Accelerators

Robust Adaptive Attitude Control of a Spacecraft

Lecture 25: Tue Nov 27, 2018

A Cross-Associative Neural Network for SVD of Nonsquared Data Matrix in Signal Processing

Introduction. Performance and Robustness (Chapter 1) Advanced Control Systems Spring / 31

6.4 Input-Output Linearization of SISO Systems

Adaptive Control with a Nested Saturation Reference Model

STATISTICAL ORBIT DETERMINATION

Scenario Optimization for Robust Design

CDS 101/110a: Lecture 10-2 Control Systems Implementation

An Iteration-Domain Filter for Controlling Transient Growth in Iterative Learning Control

Stochastic dynamical modeling:

Jim Lambers MAT 610 Summer Session Lecture 1 Notes

Linear State Feedback Controller Design

A Novel Integral-Based Event Triggering Control for Linear Time-Invariant Systems

1.10 Continuity Brian E. Veitch

INTEGRATED ARCHITECTURE OF ACTUATOR FAULT DIAGNOSIS AND ACCOMMODATION

Image restoration: numerical optimisation

UDE-based Dynamic Surface Control for Strict-feedback Systems with Mismatched Disturbances

Control Systems I. Lecture 2: Modeling and Linearization. Suggested Readings: Åström & Murray Ch Jacopo Tani

Dependence. MFM Practitioner Module: Risk & Asset Allocation. John Dodson. September 11, Dependence. John Dodson. Outline.

DISTURBANCE OBSERVER BASED CONTROL: CONCEPTS, METHODS AND CHALLENGES

NDI-BASED STRUCTURED LPV CONTROL A PROMISING APPROACH FOR AERIAL ROBOTICS

2tdt 1 y = t2 + C y = which implies C = 1 and the solution is y = 1

Robust Internal Model Control for Impulse Elimination of Singular Systems

Design of a Missile Autopilot using Adaptive Nonlinear Dynamic Inversion

Stability of Dynamical systems

Rigid-Body Attitude Control USING ROTATION MATRICES FOR CONTINUOUS, SINGULARITY-FREE CONTROL LAWS

Nonlinear Adaptive Flight Control for the X-38 Vehicle

Lecture 6 Asymptotic Structure for Four-Step Premixed Stoichiometric Methane Flames

MIT Weakly Nonlinear Things: Oscillators.

Lecture 12. Upcoming labs: Final Exam on 12/21/2015 (Monday)10:30-12:30

Math 273a: Optimization The Simplex method

Indirect Model Reference Adaptive Control System Based on Dynamic Certainty Equivalence Principle and Recursive Identifier Scheme

Short Questions (Do two out of three) 15 points each

Lecture 6 Optimization for Deep Neural Networks

Real-Time Feasibility of Nonlinear Predictive Control for Semi-batch Reactors

Computer Intensive Methods in Mathematical Statistics

Transcription:

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 and Possible Remedies 2

Summary of Dynamic Inversion Design Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore

Philosophy of Dynamic Inversion Carryout a co-ordinate transformation such that the problem appears to be linear in the transformed co-ordinates. Design the control for the linear-looking system using the linear control design techniques. Obtain the controller for the original system using an inverse transformation. Intuitively, the control design is carried out by enforcing a stable linear error dynamics. 4

Problem System Dynamics: Goal (Tracking): Special Class: (control affine & square) (, ) ( ) X = f X U Y = h X X R, U R, Y R n m p () * Y Y t, t () t * Assumption: Y is smooth ( ) ( ) ( ) X = f X + g X U p= m, gy X non-singular t 5

Dynamic Inversion Design Derive the output dynamics: Y h = X X h = X ( ) + ( ) { f X g X U} ( ) ( ) = fy X + gy X U Define Error of Tracking: Known: ( ) ( ) ( ) X = f X + g X U Y = h X h1 h1 x1 x n h X hp hp x1 x n * () ( ) ( ) E t Y t Y t 6

Dynamic inversion Select a fixed gain K > 0 such that: Kt E + K E = 0 E = e E 0, t Carry out the algebra: + = 0 ( *) ( * Y Y K Y Y ) * ( ) ( ) ( * + ) Y = fy X g X U Y K Y Y Solve for the controller: { } Y 1 * ( ) ( * ) ( ) U = gy X Y K Y Y f X 0 ( ) Usually K = diag 1/ τ, τ > 0 i i 7

Is a first-order error dynamics always enforced? Answer: No! The order of the error dynamics is dictated by the relative degree of the problem, which is defined as the number of times the output needs to be differentiated so that the control variable appears explicitly. If a second-order error dynamics needs to be enforced, then the corresponding equation is E + K E + K E = V P 0 ( ) ( 2 ξω ω ) ( ξ ω ) Usually K = diag 2, K = diag,, > 0 V i n P n i n i i i 8

When Does the Dynamic Inversion Fail? Fundamental Principle of Dynamic Inversion : 1. Differentiate y repeatedly until the input u appears. 2. Design u to cancel the nonlinearity. Q : Is it always possible to design u this way? Ans: Not necessarily! It is possible only if the relative degree is "well-defined". 9

Undefined Relative degree Undefined relative degree : It may so happen that upon sucessive differentiation of y, u appears. However, the coefficient of u may vanish at X 0, whereas it is non-zero at points arbitrarly close to X In such cases, the relative degree is undefined at X. Ex : x x 1 2 0 = + u x ρ 1 ( x, x ) 2 1 2 y = x 2 1 Then y = 2x x = 2x x, ρ : Some nonlinear function 1 1 1 2 1 2 1 2 1 1 ( ) [ ] 2 y = 2x x + 2x x = 2 x ρ x, x + u 2 + 2x 2 0 0. 10

Undefined Relative degree ( ) 2 1 ρ 1 2 2 1 y = 2 x x, x + 2 x + 2x u As x1 = 0, gy ( x) = 0. 1 ( ) g ( X ) fy X y Hence, at x = 0 the relative degree is. Note: If one choses y = x then y = x = x 1 2 (, ) y = x = ρ x x + u 1 2 2 and the coefficient of u = 1 0 globally. 1, not defined globally In this case the relative degree is well defined. 11

Advantages of DI design Simple design: No need of tedious gain scheduling (hence, sometimes dynamic inversion is known as a universal gain scheduling design). Easy online implementation: It leads to a closed form solution for the controller. Asymptotic (rather exponential) stability is guaranteed for the Error dynamics (subjected to the control availability, this is true globally ). No problem if parameters are updated (the updated values can simply be used in the formula). 12

Issues and Remedies in Dynamic Inversion Design Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore

Issues of DI design Does the inverse exist for all time? What if the problem is non-square? (i.e. the number of controllers are not equal to the number of performance outputs)? Is the internal dynamics stable? Is it sensitive to modeling inaccuracies (e.g. parameter inaccuracies)? 14

Does the inverse exist for all time? Answer: Not Necessarily! Strategies: Do not update the controller unless ( ( )) abs gy X > ε (it may lead to performance degradation locally) Reformulate the problem Some advanced techniques do address this issue partly in the control design process itself 15

What if the output dynamics is non-square? Case - 1: m < p No. of controllers < No. of objectives Theorem: Perfect tracking is not possible in general. Case - 2: m > p No. of controllers > No. of objectives Additional objectives can be demanded in the design One such approach: Optimal Dynamic Inversion 16

Optimal Dynamic Inversion: Meeting the tracking objective with minimum Control Select a (fixed) gain K > 0 such that: Kt E + K E = 0 E = e E 0, t Carry out the algebra: Hence: + = 0 ( *) ( * Y Y K Y Y ) * ( ) ( ) ( * + ) Y = fy X g X U Y K Y Y { } Y * ( ) ( * = ) ( ) g X U Y K Y Y f X Y A b 0 ( ) Usually K = diag 1/ τ, τ > 0 i i 17

Optimal Dynamic Inversion: Meeting the tracking objective with minimum Control Minimize: Subjected to: Solution: J = AU 1 2 = ( T U RU) b A g ( ) Y X ( ) ( ) Y * * b Y K Y Y f X 1 Augmented Performance Index: J = U RU + AU b 2 J J Necessary Conditions: = 0, = 0 U λ ( T ) T λ ( ) 18

Optimal Dynamic Inversion: Meeting the tracking objective with minimum Control T Solution: RU + A λ = 0 U AU = b = 0 1 T R A λ ( 1 T λ ) A R A = b λ = ( 1 T AR A ) 1 b 1 T 1 1 U = R A AR A b ( T) 19

Useful Features of ODI The enforced error dynamics is satisfied exactly in ODI: Asymptotic tracking is guaranteed ODI gives a platform for optimal control allocation as well (both time-wise as well as location-wise). Pseudo inverse 1 T ( 1 T) 1 T ( T) 1 U R A AR A + = b= A AA b= A b Special case: if R = I, then + { } Y U= gy X Y K Y Y f X ( ) * ( *) ( ) 20

References on ODI Radhakant Padhi and Mangal Kothari, An Optimal Dynamic Inversion Based Neuro-Adaptive Approach for Treatment of Chronic Myelogenous Leukemia, Computer Methods and Programs in Biomedicine, Vol.87, 2007, pp.208-228. Radhakant Padhi and S. N. Balakrishnan, Optimal Dynamic Inversion Control Design for a Class of Nonlinear Distributed Parameter Systems with Continuous and Discrete Actuators, IET Control Theory and Applications, Vol. 1, Issue 6, Nov. 2007, pp.1662-1671. 21

Is the internal dynamics always stable? Answer: Not Necessarily! The control solution is meaningless, unless this issue is addressed explicitly. Standard Results: If the relative degree of the problem is equal to the number of states, then the internal dynamics is stable. Asymptotic stability of zero dynamics is sufficient for local input-to-state stability of internal dynamics. If analytical justification is not possible, extensive simulation studies must be carried out. A nonlinear system is said to be minimum phase if its zero dynamics is stable. 22

Is there a procedure for tracking control design for non-minimum phase systems? Answer: Fortunately YES! An Output Redefinition Technique is available, where the internal dynamics is first made locally stable before achieving the tracking objective. References: E. M. Wallner and K. H. Well, Attitude Control of a Re-entry Vehicle with Internal Dynamic, Journal of Guidance, Control, and Dynamics, Vol. 26, No. 6, Nov-Dec 2003. (Application oriented: Good for practicing engineers. It also contains information about zero dynamics and the relationship of its stability with the stability of internal dynamics) S. Gopalswamy and J. K. Hedrick, Tracking Nonlinear Non-Minimum Phase Systems Using Sliding Control, International Journal of Control, Vol. 57, No. 5, 1993, pp.1141-1158. (Theoretical foundations for output redefinition technique: Good for academicians and further research) 23

Philosophy of Output Redefinition Technique Key Idea: Do not aim for perfect tracking! Approximate tracking is OK. ( ) Let the original output y = h X have unstable zero dynamics. Then redefine output y1 = h1( X) such that the resulting zero-dynamics for y tracking is stable. 1 Next, design the controller such that y ( t ) y ( t ) exactly. Then this implies good tracking of the original output provided certain conditions are met. 1 d 24

Output Redefinition Technique: An Illustrative Simple Example ( 1 s/ b) B0 ( s) As ( ) Let y = u where b> 0 and zeros of B s are in the left half plane. 0 ( ) This system has a right-half plane zero at s= b. To avoid the problem of unstable zero dynamics, let us consider the control of a nominal output y, which is defined as y 1 B ( s) ( ) u 0 = As 1 25

Output Redefinition Technique: An Illustrative Simple Example Design ut ( ) such that y( t) y( t) asymptotically. Let us then analyze the case when y ( t) = y ( t). y() t = 1 s/ b y () t = 1 s/ b y () t ( ) ( ) 1 1 [ ] ( ) Error e() t y() t y () t = y ()/ t b d Hence et ( ) is bounded as long as y d ( t) is bounded. Further, if b is large (as compared to y ( t)) et () remains small. d d d 1 d d 26

Output Redefinition Technique: An Illustrative Simple Example B ( ) 0 s An alternate choice: y2 = u As ( )( 1 + sb / ) Let us then analyze the case when y ( t) = y ( t). y() t = 1 s/ b 1 + s/ b y () t = 1 s / b y () t ( )( ) ( 2 2) 2 [ ] ( 2 = ) Error et () yt () y() t y()/ t b Hence et ( ) is bounded as long as y( t) d d 2 Furthermore, if b is large (as compared to yd ( t)) et ( ) remains small. d 2 d d is bounded. 27

Output Redefinition Technique: An Illustrative Simple Example Out of the two choices, y 2 choice leads to a better tracking performance, provided and vice-versa. y d < b y d 28

Another Approach to Deal with Non-minimum phase Systems Philosophy: Neglect the terms containing the input u, while doing successive iterations of the output until the number of successive differentiations becomes n. So there is "approximately" no internal dynamics. Design controller based on this "approximate" input-outpt dynamics. Note: this approach works as long as the coefficients of u at the intermediate steps are small i.e. the system is "weakly non-minimum phase". 29

Other Ideas to Deal with Non-minimum phase Systems (1) Modify the desired trajectories (this may not be feasible and/or it may limit the operating zone. (2) Modify the plant itself. This may be possible by relocation/addition of actuators/sensors. May also be possible by physical modification of the plant (e.g. placing the control surfaces at a different location in an aircraft). 30

Is DI design sensitive to modeling and parameter inaccuracies? Answer: Very much YES! Solution: Augmenting the Dynamic Inversion design with other Robust and/or Adaptive control design techniques. One Popular Approach: Neuro-Adaptive Control Design 31

Motivations for Neuro-Adaptive Design Perfect system modeling is difficult Sources of imperfection can arise from: Unmodelled dynamics (missing algebraic terms in the model) Inaccurate knowledge of system parameters Change of system parameters/dynamics during operation The adaptive control design should be able to learn the unknown function through neural network(s) and then compensate for this unknown error. 32

A Demonstrative Example ( ) 0.1sin( ) x = 2sin x + x Known part of actual system (nominal system) Unknown part of actual system ( x) csin( x) = 2sin +Δ Weight Basis Function 33

Summary: Dynamic Inversion Design DI design offers several advantages: Its is nonlinear design: No linearization of system dynamics is necessary It is a promising substitute to gain scheduling philosophy. Quite often a constant gain is found to be satisfactory, even though gain scheduling can still be an option for performance improvement It assures perfect tracking (i.e. asymptotic stability of error dynamics) under ideal assumption of perfect knowledge of the system dynamics DI offers a closed form solution for the control and hence it can easily be implemented online 34

Summary: Dynamic Inversion Design In DI design, one has to be careful about the following important issues: Non-existence of the matrix inversion If this is a local problem, do not update the control Stability of internal dynamics If the analysis fails, then either reformulate the problem or opt for output redefinition technique Robustness with respect to modeling inaccuracies Augment the DI design with either robust or adaptive control design ideas). Neuro-adaptive technique has recently evolved as one such promising idea. 35

36