Control Systems Design, SC4026. SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft

Similar documents
Introduction to Modern Control MT 2016

Introduction to Modern Control MT 2016

Freeman Dyson on describing the predictions of his model for meson-proton scattering to Enrico Fermi in 1953 [Dys04].

In the presence of viscous damping, a more generalized form of the Lagrange s equation of motion can be written as

Review: control, feedback, etc. Today s topic: state-space models of systems; linearization

CDS 101/110a: Lecture 2.1 Dynamic Behavior

Control Systems Design, SC4026. SC4026 Fall 2010, dr. A. Abate, DCSC, TU Delft

CDS 101/110a: Lecture 2.1 Dynamic Behavior

Control Systems Design, SC4026. SC4026 Fall 2010, dr. A. Abate, DCSC, TU Delft

Control Systems Design, SC4026. SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft

An Introduction for Scientists and Engineers SECOND EDITION

ELEC4631 s Lecture 2: Dynamic Control Systems 7 March Overview of dynamic control systems

System simulation using Matlab, state plane plots

Lagrange s Equations of Motion and the Generalized Inertia

EE 16B Midterm 2, March 21, Name: SID #: Discussion Section and TA: Lab Section and TA: Name of left neighbor: Name of right neighbor:

EE Homework 3 Due Date: 03 / 30 / Spring 2015

CDS 101: Lecture 2.1 System Modeling

Predictability: Does the Flap of a Butterfly s Wings in Brazil set off a Tornado

Fundamentals Physics. Chapter 15 Oscillations

Modeling and Experimentation: Compound Pendulum

Raktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Preliminaries

Chapter 2 SDOF Vibration Control 2.1 Transfer Function

Dr Ian R. Manchester Dr Ian R. Manchester AMME 3500 : Review

28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod)

Analog Signals and Systems and their properties

CDS 101: Lecture 5-1 Reachability and State Space Feedback

Robotics. Dynamics. University of Stuttgart Winter 2018/19

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

School of Engineering Faculty of Built Environment, Engineering, Technology & Design

System Modeling. Chapter Modeling Concepts

Definition of Reachability We begin by disregarding the output measurements of the system and focusing on the evolution of the state, given by

Numerics and Control of PDEs Lecture 1. IFCAM IISc Bangalore

CDS 101: Lecture 2.1 System Modeling. Lecture 1.1: Introduction Review from to last Feedback week and Control

Matlab-Based Tools for Analysis and Control of Inverted Pendula Systems

Linearization problem. The simplest example

Lecture A1 : Systems and system models

MODELING OF CONTROL SYSTEMS

Overview of the Seminar Topic

Robotics. Dynamics. Marc Toussaint U Stuttgart

EECS C128/ ME C134 Final Wed. Dec. 15, am. Closed book. Two pages of formula sheets. No calculators.

Non-Linear Response of Test Mass to External Forces and Arbitrary Motion of Suspension Point

Topic # Feedback Control Systems

General procedure for formulation of robot dynamics STEP 1 STEP 3. Module 9 : Robot Dynamics & controls

ECEN 420 LINEAR CONTROL SYSTEMS. Lecture 6 Mathematical Representation of Physical Systems II 1/67

D(s) G(s) A control system design definition

Module 03 Modeling of Dynamical Systems

Introduction to Controls

Automatic Control Systems. -Lecture Note 15-

The Inverted Pendulum

CALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems

Systems Engineering/Process Control L1

Lecture 25: Tue Nov 27, 2018

P321(b), Assignement 1

Swing-Up Problem of an Inverted Pendulum Energy Space Approach

Systems of Ordinary Differential Equations

Linear System Theory. Wonhee Kim Lecture 1. March 7, 2018

sc Control Systems Design Q.1, Sem.1, Ac. Yr. 2010/11

CALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems

STATE VARIABLE (SV) SYSTEMS

ME 132, Dynamic Systems and Feedback. Class Notes. Spring Instructor: Prof. A Packard

Reglerteknik, TNG028. Lecture 1. Anna Lombardi

State Space Representation

DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED.

Chapter 13 Lecture. Essential University Physics Richard Wolfson 2 nd Edition. Oscillatory Motion Pearson Education, Inc.

Lecture. Math Dylan Zwick. Spring Simple Mechanical Systems, and the Differential Equations

Designing Information Devices and Systems II Spring 2017 Murat Arcak and Michel Maharbiz Homework 9

Lecture 20 Aspects of Control

Lecture 9 Nonlinear Control Design. Course Outline. Exact linearization: example [one-link robot] Exact Feedback Linearization

SRV02-Series Rotary Experiment # 7. Rotary Inverted Pendulum. Student Handout

Today s goals So far Today 2.004

Transverse Linearization for Controlled Mechanical Systems with Several Passive Degrees of Freedom (Application to Orbital Stabilization)

PHYSICS 44 MECHANICS Homework Assignment II SOLUTION

Prüfung Regelungstechnik I (Control Systems I) Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam!

CL Digital Control

MAE143a: Signals & Systems (& Control) Final Exam (2011) solutions

Design and Comparison of Different Controllers to Stabilize a Rotary Inverted Pendulum

Analysis and Control of Multi-Robot Systems. Elements of Port-Hamiltonian Modeling

Outline. Classical Control. Lecture 1

Lecture 9 Nonlinear Control Design

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

ACM/CMS 107 Linear Analysis & Applications Fall 2016 Assignment 4: Linear ODEs and Control Theory Due: 5th December 2016

Massachusetts Institute of Technology. Department of Electrical Engineering and Computer Science : MULTIVARIABLE CONTROL SYSTEMS by A.

CDS 101/110a: Lecture 1.1 Introduction to Feedback & Control. CDS 101/110 Course Sequence

MEM04: Rotary Inverted Pendulum

DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED.

Contents. Dynamics and control of mechanical systems. Focus on

MEAM 510 Fall 2012 Bruce D. Kothmann

Dynamical Systems & Lyapunov Stability

Math 128A Spring 2003 Week 12 Solutions

Lecture 11 FIR Filters

CDS 101: Lecture 2.1 System Modeling

APPPHYS 217 Tuesday 6 April 2010

Neural Networks Lecture 10: Fault Detection and Isolation (FDI) Using Neural Networks

The Modeling of Single-dof Mechanical Systems

x(n + 1) = Ax(n) and y(n) = Cx(n) + 2v(n) and C = x(0) = ξ 1 ξ 2 Ex(0)x(0) = I

Consider a particle in 1D at position x(t), subject to a force F (x), so that mẍ = F (x). Define the kinetic energy to be.

Networks in systems biology

Imaginary. Axis. Real. Axis

Stabilization of a Chain of Exponential Integrators Using a Strict Lyapunov Function

Linear control of inverted pendulum

Transcription:

Control Systems Design, SC4026 SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft

Lecture 1 The concept of feedback The role of a controller What is a state? The concept of model in systems engineering: a few examples SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft 1

reasoning about a feedback system is difficult because the first system influences the second and the second system influences the first, leading to a circular argument. This makes reasoning based on cause and effect tricky, and it is necessary to analyze the system as a whole. A consequence of this is that the behavior of feedback systems is often counterintuitive, and it is therefore necessary to resort to formal methods to understand them. The concept of feedback Figure 1.1 illustrates in block diagram form the idea of feedback. We often use System 1 u System 2 y r System 1 u System 2 y (a) Closed loop (b) Open loop Figure 1.1: Open and closed loop systems. (a) The output of system 1 is used as the input of system 2 2, and the output of system 2 becomes the input of system CHAPTER1, 1. creating INTRODUCTION a closed loop system. (b) The interconnection between system 2 and system 1 is removed, and the system is said to be open loop. Figure 1.2: The centrifugal governor and the steam engine. The centrifugal governor on the left consists of a set of flyballs that spread apart as the speed of the engine increases. The SC4026 Fall 2009, dr. A. steam Abate, engine DCSC, on the TU right Delft uses a centrifugal governor (above and to the left of the flywheel) 2 to regulate its speed. (Credit: Machine a Vapeur Horizontale de Philip Taylor [1828].) the terms open loop and closed loop when referring to such systems. A system

The concept of feedback 16 CHAPTER 1. INTRODUCTION Figure 1.12: The wiring diagram of the growth-signaling circuitry of the mammalian cell [HW00]. The major pathways that are thought to play a role in cancer are indicated SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft 3 in the diagram. Lines represent interactions between genes and proteins in the cell. Lines ending in arrowheads indicate activation of the given gene or pathway; lines ending in a T-shaped head indicate repression. (Used with permission of Elsevier Ltd. and the authors.)

The concept of feedback 2.2. STATE SPACE MODELS 39 160 140 120 100 80 60 40 20 1845 1855 1865 1875 1885 1895 Hare Lynx 1905 1915 1925 1935 Figure 2.6: Predator versus prey. The photograph on the left shows a Canadian lynx and a snowshoe hare, the lynx s primary prey. The graph on the right shows the populations of hares and lynxes between 1845 and 1935 in a section of the Canadian Rockies [Mac37]. The data were collected on an annual basis over a period of 90 years. (Photograph copyright Tom and Pat Leeson.) discrete-time index (e.g., the month number), we can write H[k + 1] = H[k] + b r (u)h[k] al[k]h[k], SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft 4 L[k + 1] = L[k] + cl[k]h[k] d f L[k], (2.13) where b r (u) is the hare birth rate per unit period and as a function of the food

The concept of feedback SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft 5

The role of a controller Control: the use of algorithms and feedback in engineered systems 4 CHAPTER 1. INTRODUCTION noise external disturbances noise Actuators System Sensors Output Process Clock D/A Computer A/D Filter Controller operator input Figure 1.3: Components of a computer-controlled system. The upper dashed box represents SC4026 Fall 2009, the process dr. A. Abate, dynamics, DCSC, which TU Delft include the sensors and actuators in addition to the dynamical 6 system being controlled. Noise and external disturbances can perturb the dynamics of the process. The controller is shown in the lower dashed box. It consists of a filter and analog-todigital (A/D) and digital-to-analog (D/A) converters, as well as a computer that implements the control algorithm. A system clock controls the operation of the controller, synchronizing

CHAPTER 2. SYSTEM MODELING State-space Models: a First Example q c(q) q(t) = F m m q(t) = 1 m ( c( q) kq + u) k ############################ ss system with nonlinear damping. The position of the mass is denoted sponding to the rest position of the spring. The forces on the mass are pring with spring constant k and a damper with force dependent on the SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft 7 chanics! "

State-space Models: a First Example Block diagram for input-output relationship Input signal: u(t) ############################ Output signal: y(t) = q(t)! "!! SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft 8

! " State-space Models: a First Example Introduce state variables (integrator! outputs): x 1 (t) = q(t) and! x 2 (t) = q(t)! "!! SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft 9

Obtain system of first-order ODE: { ẋ1 (t) = x 2 (t) ẋ 2 (t) = 1 m ( c(x 2(t)) kx 1 (t) + u(t)) To find solution, need two initial conditions Note presence of linear & nonlinear parts SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft 10

State-space Models: a Second Example Predator-prey model (introduced before) Now we aim at obtaining a quantitative, abstract simplification of the actual dynamics State variables: time-dependent population level for the lynxes: l(t), t 0 and for the hares: h(t), t 0 Control Input: hare birth rate b(u), function of food Outputs: population levels l(t), h(t) SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft 11

cedure calls (RPCs). The server maintains a log of statistics of completed requests. The total number of requests being served, called RIS (RPCs in server), is also measured. The load on the server is controlled by a parameter called MaxUsers, Model parameters: Mortality rate d. Interaction rates a, c Dynamical model: { ḣ(t) = b(u)h(t) a l(t)h(t) l(t) = c l(t)h(t) d l(t) Simulation outputs of developed model: 40 CHAPTER 2. SYSTEM MODELING 250 200 Hares Lynxes Population 150 100 50 0 1850 1860 1870 1880 1890 1900 1910 1920 Year Figure 2.7: Discrete-time simulation of the predator prey model (2.13). Using the parameters a = c = 0.014, b r (u) = 0.6 and d = 0.7 in equation (2.13), the period and magnitude of the SC4026 Fall 2009, dr. A. Abate, lynx anddcsc, hare population TU Delft cycles approximately match the data in Figure 2.6. 12

State-space Models: a Third Example Control of inverted pendulum on moving cart (balance system, e.g. Segway) 36 CHAPTER 2. SYSTEM MODELING m θ l F M p (b) Saturn rocket (c) Cart pendulum system tems. (a) Segway Personal Transporter, (b) Saturn rocket and (c) art. Each of these examples uses forces at the bottom of the system SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft a generalization of the spring mass system we saw earlier. ics for a mechanical system in the general form (a) Segway (b) Saturn rock Figure 2.5: Balance systems. (a) Segway Pers inverted pendulum on a cart. Each13 of these exam to keep it upright.

State-space Models: a Third Example Dynamics can be derived via Lagrange equations States: position p and angle θ Kinetic energy: T M = 1 2 Mṗ2, T m = 1 2 m(ṗ2 + 2lṗ θ cos θ + l 2 θ2 ) Potential energy: V = mgl cos θ Overall state q = (p, θ), input u = F, output y = (p, θ) SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft 14

Lagrangian: L = T (q, q) V (q) = (T M + T m ) V Lagrange s Equations: d L dt q L q = [ F 0 ] Obtain [ (M + m) ml cos θ ml cos θ ml 2 ] [ p θ ] + [ ml sin θ θ2 mgl sin θ ] = [ F 0 ] Can synthetically write the dynamics as: M(q) q + K(q, q) = Bu SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft 15

Notice that we have built a nonlinear model The model will be linearized in Lecture 2 SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft 16

Models from Experiments How to construct state-space models? physics-based model (conservation laws, other physical laws, material properties, physical geometry and dimensions) models based on known interactions and properties (e.g.: energy-based models, stochiometric models) Models from experiments (data driven): use of transfer functions, possibly derived from experiments; measurement of model properties and use of fitting (connection to later part of class) SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft 17