Contents lecture 6 2(17) Automatic Control III. Summary of lecture 5 (I/III) 3(17) Summary of lecture 5 (II/III) 4(17) H 2, H synthesis pros and cons:

Similar documents
Contents 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

Stability of Dynamical systems

Classification of Phase Portraits at Equilibria for u (t) = f( u(t))

Chapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12

Calculus and Differential Equations II

Exam in Systems Engineering/Process Control

Department of Mathematics IIT Guwahati

4 Second-Order Systems

MathQuest: Differential Equations

Section 5.4 (Systems of Linear Differential Equation); 9.5 Eigenvalues and Eigenvectors, cont d

Section 9.3 Phase Plane Portraits (for Planar Systems)

Nonlinear Autonomous Systems of Differential

STABILITY. Phase portraits and local stability

+ i. cos(t) + 2 sin(t) + c 2.

MATH 215/255 Solutions to Additional Practice Problems April dy dt

Chapter 1, Section 1.2, Example 9 (page 13) and Exercise 29 (page 15). Use the Uniqueness Tool. Select the option ẋ = x

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

MAT 22B - Lecture Notes

Understand the existence and uniqueness theorems and what they tell you about solutions to initial value problems.

Solutions to Dynamical Systems 2010 exam. Each question is worth 25 marks.

Differential Equations 2280 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 2015 at 12:50pm

Exam in Systems Engineering/Process Control

AC&ST AUTOMATIC CONTROL AND SYSTEM THEORY SYSTEMS AND MODELS. Claudio Melchiorri

Lesson 4: Non-fading Memory Nonlinearities

Models Involving Interactions between Predator and Prey Populations

Complex Dynamic Systems: Qualitative vs Quantitative analysis

Phase portraits in two dimensions

Intro. Computer Control Systems: F8

A plane autonomous system is a pair of simultaneous first-order differential equations,

Linearization of Differential Equation Models

Math 312 Lecture Notes Linear Two-dimensional Systems of Differential Equations

Appendix: A Computer-Generated Portrait Gallery

Copyright (c) 2006 Warren Weckesser

Exercises Automatic Control III 2015

Lecture 38. Almost Linear Systems

Math 266: Phase Plane Portrait

Math 232, Final Test, 20 March 2007

2.10 Saddles, Nodes, Foci and Centers

Solutions to Math 53 Math 53 Practice Final

Lecture 5. Outline: Limit Cycles. Definition and examples How to rule out limit cycles. Poincare-Bendixson theorem Hopf bifurcations Poincare maps

154 Chapter 9 Hints, Answers, and Solutions The particular trajectories are highlighted in the phase portraits below.

Subject: Optimal Control Assignment-1 (Related to Lecture notes 1-10)

Problem set 7 Math 207A, Fall 2011 Solutions

Math 331 Homework Assignment Chapter 7 Page 1 of 9

Math 3301 Homework Set Points ( ) ( ) I ll leave it to you to verify that the eigenvalues and eigenvectors for this matrix are, ( ) ( ) ( ) ( )

Math 215/255 Final Exam (Dec 2005)

Lecture Notes for Math 251: ODE and PDE. Lecture 27: 7.8 Repeated Eigenvalues

MA 527 first midterm review problems Hopefully final version as of October 2nd

Automatic Control Systems theory overview (discrete time systems)

Problem set 6 Math 207A, Fall 2011 Solutions. 1. A two-dimensional gradient system has the form

EE222 - Spring 16 - Lecture 2 Notes 1

A review of stability and dynamical behaviors of differential equations:

ODE, part 2. Dynamical systems, differential equations

In these chapter 2A notes write vectors in boldface to reduce the ambiguity of the notation.

Nonlinear Control Lecture 2:Phase Plane Analysis

Computer Problems for Methods of Solving Ordinary Differential Equations

Nonlinear dynamics & chaos BECS

Linear Planar Systems Math 246, Spring 2009, Professor David Levermore We now consider linear systems of the form

Return Difference Function and Closed-Loop Roots Single-Input/Single-Output Control Systems

Introduction to Dynamical Systems

Identification Methods for Structural Systems

Examples include: (a) the Lorenz system for climate and weather modeling (b) the Hodgkin-Huxley system for neuron modeling

OPTIMAL CONTROL. Sadegh Bolouki. Lecture slides for ECE 515. University of Illinois, Urbana-Champaign. Fall S. Bolouki (UIUC) 1 / 28

= F ( x; µ) (1) where x is a 2-dimensional vector, µ is a parameter, and F :

Modeling and Analysis of Dynamic Systems

Autonomous systems. Ordinary differential equations which do not contain the independent variable explicitly are said to be autonomous.

MCE693/793: Analysis and Control of Nonlinear Systems

7 Planar systems of linear ODE

Math 308 Final Exam Practice Problems

MS-E2133 Systems Analysis Laboratory II Assignment 2 Control of thermal power plant

Suppose that we have a specific single stage dynamic system governed by the following equation:

7 Two-dimensional bifurcations

Question: Total. Points:

Basic Procedures for Common Problems

CDS 101 Precourse Phase Plane Analysis and Stability

Part II Problems and Solutions

Linear Systems of ODE: Nullclines, Eigenvector lines and trajectories

Introduction to gradient descent

APPM 2360: Final Exam 10:30am 1:00pm, May 6, 2015.

EG4321/EG7040. Nonlinear Control. Dr. Matt Turner

State Regulator. Advanced Control. design of controllers using pole placement and LQ design rules

1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point

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

MODELING OF CONTROL 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

Homogeneous Constant Matrix Systems, Part II

4F3 - Predictive Control

Third In-Class Exam Solutions Math 246, Professor David Levermore Thursday, 3 December 2009 (1) [6] Given that 2 is an eigenvalue of the matrix

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

LECTURE 8: DYNAMICAL SYSTEMS 7

Computational Physics (6810): Session 8

CDS 101/110: Lecture 3.1 Linear Systems

Boulder School for Condensed Matter and Materials Physics. Laurette Tuckerman PMMH-ESPCI-CNRS

Nonlinear Control. Nonlinear Control Lecture # 18 Stability of Feedback Systems

Laplace Transform Part 1: Introduction (I&N Chap 13)

Uncertainty and Robustness for SISO Systems

CHAPTER 1 Basic Concepts of Control System. CHAPTER 6 Hydraulic Control System

APPPHYS 217 Tuesday 6 April 2010

Time Response of Systems

Nonlinear System Analysis

Transcription:

Contents lecture 6 (7) Automatic Control III Lecture 6 Linearization and phase portraits. Summary of lecture 5 Thomas Schön Division of Systems and Control Department of Information Technology Uppsala University.. General properties 3. Linearization and stationary points. Phase portraits Email: thomas.schon@it.uu.se, www: user.it.uu.se/~thosc Summary of lecture 5 (I/III) 3(7) Summary of lecture 5 (II/III) (7) H, H synthesis pros and cons: H and H synthesis: Make W u G wu, W S S, W T T small. H : Minimize ( W u G wu + W SS + W TT ) dω. H : Set an upper bound for W u G wu, W S S, W T T for all ω. Results in algebraic Riccati equations. (+) Directly handles the specifications on S, T and G wu (+) Let us know when certain specifications are impossible to achieve (via γ). (+) Easy to handle several different specifications (in the frequency domain) (-) Can be hard to control the behaviour in the time domain in detail. (-) Often results in complex controllers (number of states in the controller = number of states in G, W u, W S, W T ).

Summary of lecture 5 (III/III) 5(7) DC motor with a saturated control signal (I/V) 6(7) r ũ u y Σ K(s + ) s + s(s + ) Linear multivariable controller synthesis summary:. Perform an RGA analysis. Use simple SISO controllers of PID type if the RGA analysis indicates that it might be possible. 3. Otherwise make use of LQ, MPC or H /H synthesis. DC motor controlled using a lead controller. We want to control the motor angle. The saturation ũ ũ, u = sat(ũ) = ũ >, ũ <, renders the system nonlinear. DC motor with a saturated control signal (II/V) 7(7) DC motor with a saturated control signal (III/V) 8(7) Step responses for two different amplitudes of the reference r signal. Blue: Amplitude Red: Amplitude 5 (scaled with /5) Conclusion: The step response is amplitude dependent. If the system would have been linear the two step responses would have coincided....8.6.. 6 8 Red: Reference signal r. Blue: Output signal y. 9 8 7 6 5 3 6 8. 6 8 Both the ramp and the sine responses are roughly the same as for a linear system..6.5..3...

DC motor with a saturated control signal (IV/V) 9(7) Red: r. Blue: y. Green: y when r is a ramp (same as on the previous slide). DC motor with a saturated control signal (V/V) (7) Red: before the saturation (ũ). Blue: after the saturation (u). 8 3 6 6 8 Something happens here: There is no sine present in the response and the ramp error has increased... This violates the superposition principle and the frequency fidelity! 6 8 Linearization of a nonlinear system (7) We can approximate a nonlinear system ẋ = f (x, u), y = h(x, u), by linearizing the system around an equilibrium (stationary) point (x, u ). Intuitively this amounts to approximating the system by a flat hyperplane (straight line in the scalar case). Let x(t) = x(t) x, u(t) = u(t) u, y(t) = y(t) y. A Taylor expansion (only keeping the linear terms) results in A B {}}{{}}{ d dt x = f (x, u ) x y = h(x, u ) x } {{ } C x + f (x, u ) u x + f (x, u ) u } {{ } D u, u. Phase portraits for linear systems (7) Sign: Is the solution moving towards the origin or away from the origin (along the eigenvector)? Relative size: fast and slow eigenvectors, which is dominating the solution behaviour for t and t? Complex/real: Complex conjugated eigenvalues results in circles and spirals. Cases to consider:. Two distinct real valued eigenvalues implies two eigenvectors.. One multiple eigenvalue implies one eigenvector. 3. Complex eigenvalues.

Two distinct eigenvalues with the same sign 3(7) The solution is x(t) = c e λ t v + c e λ t v. Stable node: For eigenvalues λ < λ <. The first term dominates for small t, the second term dominates for large t. Unstable node: For eigenvalues < λ < λ. Also here, the first term dominates for small t, the second term dominates for large t. Two distinct eigenvalues with opposite sign (7) For eigenvalues λ < < λ (with corresponding eigenvectors v, v ) the solution is x(t) = c e λ t v + c e λ t v. Trajectories close to v will approach the origin. v is called the stable eigenvector. Trajectories close to v will move away from the origin. v is called the unstable eigenvector. A multiple eigenvalue 5(7) For multiple eigenvalues λ = λ. Two examples in 3D 6(7) Example of a generalization to 3D..8.8.6.6.....5.5.5.5.5...6.8 Stable node (unstable: change direction). Stable star node (unstable: change direction). Left: Focus + one real eigenvalue. Right: Three real eigenvalues.

A few concepts to summarize lecture 6 7(7) Equilibrium points: An equilibrium point is a point x, u where the system is in rest, i.e. f (x, u ) =. Also referred to as stationary points. Linearization: Find a Taylor expansion of the nonlinear system around an equilibrium point and only keep the linear parts. This means that we are approximating the system using a flat hyperplane. Phase plane: A two dimensional state space that is simple to visualize graphically. Phase portraits: A plot where one state variable is plotted against another state variable. Limit cycle: A limit cycle is a closed trajectory in phase space having the property that at least one other trajectory spirals into it either as time approaches ± infinity.