Systems Analysis and Control

Similar documents
Systems Analysis and Control

Systems Analysis and Control

Systems Analysis and Control

Systems Analysis and Control

Step Response Analysis. Frequency Response, Relation Between Model Descriptions

Systems Analysis and Control

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering Dynamics and Control II Fall 2007

Response to a pure sinusoid

Systems Analysis and Control

Systems Analysis and Control

Systems Analysis and Control

Plan of the Lecture. Goal: wrap up lead and lag control; start looking at frequency response as an alternative methodology for control systems design.

Systems Analysis and Control

Systems Analysis and Control

Systems Analysis and Control

Systems Analysis and Control

100 (s + 10) (s + 100) e 0.5s. s 100 (s + 10) (s + 100). G(s) =

Frequency methods for the analysis of feedback systems. Lecture 6. Loop analysis of feedback systems. Nyquist approach to study stability

Homework 7 - Solutions

Control Systems I Lecture 10: System Specifications

Review: transient and steady-state response; DC gain and the FVT Today s topic: system-modeling diagrams; prototype 2nd-order system

Frequency Response of Linear Time Invariant Systems

Course roadmap. Step response for 2nd-order system. Step response for 2nd-order system

Transient response via gain adjustment. Consider a unity feedback system, where G(s) = 2. The closed loop transfer function is. s 2 + 2ζωs + ω 2 n

EE C128 / ME C134 Fall 2014 HW 6.2 Solutions. HW 6.2 Solutions

ME 375 EXAM #1 Friday, March 13, 2015 SOLUTION

The requirements of a plant may be expressed in terms of (a) settling time (b) damping ratio (c) peak overshoot --- in time domain

EE C128 / ME C134 Fall 2014 HW 8 - Solutions. HW 8 - Solutions

r + - FINAL June 12, 2012 MAE 143B Linear Control Prof. M. Krstic

ROOT LOCUS. Consider the system. Root locus presents the poles of the closed-loop system when the gain K changes from 0 to. H(s) H ( s) = ( s)

MAS107 Control Theory Exam Solutions 2008

DESIGN USING TRANSFORMATION TECHNIQUE CLASSICAL METHOD

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

8.1.6 Quadratic pole response: resonance

Raktim Bhattacharya. . AERO 422: Active Controls for Aerospace Vehicles. Frequency Response-Design Method

EECS C128/ ME C134 Final Wed. Dec. 14, am. Closed book. One page, 2 sides of formula sheets. No calculators.

Root Locus Design Example #3

Radar Dish. Armature controlled dc motor. Inside. θ r input. Outside. θ D output. θ m. Gearbox. Control Transmitter. Control. θ D.

Modern Control Systems

Sinusoidal Forcing of a First-Order Process. / τ

16.30/31, Fall 2010 Recitation # 2

( ) Frequency Response Analysis. Sinusoidal Forcing of a First-Order Process. Chapter 13. ( ) sin ω () (

Frequency Response Analysis

Outline. Classical Control. Lecture 5

Studio Exercise Time Response & Frequency Response 1 st -Order Dynamic System RC Low-Pass Filter

Overview of Bode Plots Transfer function review Piece-wise linear approximations First-order terms Second-order terms (complex poles & zeros)

Asymptotic Bode Plot & Lead-Lag Compensator

Process Control & Instrumentation (CH 3040)

ELECTRONICS & COMMUNICATIONS DEP. 3rd YEAR, 2010/2011 CONTROL ENGINEERING SHEET 5 Lead-Lag Compensation Techniques

EECS C128/ ME C134 Final Thu. May 14, pm. Closed book. One page, 2 sides of formula sheets. No calculators.

12.7 Steady State Error

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

e iωt dt and explained why δ(ω) = 0 for ω 0 but δ(0) =. A crucial property of the delta function, however, is that

e iωt dt and explained why δ(ω) = 0 for ω 0 but δ(0) =. A crucial property of the delta function, however, is that

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

Automatic Control 2. Loop shaping. Prof. Alberto Bemporad. University of Trento. Academic year

Discrete Systems. Step response and pole locations. Mark Cannon. Hilary Term Lecture

INTRODUCTION TO DIGITAL CONTROL

Module 3F2: Systems and Control EXAMPLES PAPER 2 ROOT-LOCUS. Solutions

Dynamic System Response. Dynamic System Response K. Craig 1

AMME3500: System Dynamics & Control

Time Series Analysis. Solutions to problems in Chapter 7 IMM

FREQUENCY-RESPONSE DESIGN

CHAPTER 7 : BODE PLOTS AND GAIN ADJUSTMENTS COMPENSATION

Chapter 6 - Solved Problems

Homework 11 Solution - AME 30315, Spring 2015

Introduction to Feedback Control

4.5. Applications of Trigonometry to Waves. Introduction. Prerequisites. Learning Outcomes

APPLICATIONS FOR ROBOTICS

ECE 486 Control Systems

Frequency domain analysis

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

Systems Analysis and Control

EE 4343/ Control System Design Project LECTURE 10

Dynamic circuits: Frequency domain analysis

Topic # Feedback Control. State-Space Systems Closed-loop control using estimators and regulators. Dynamics output feedback

GATE EE Topic wise Questions SIGNALS & SYSTEMS

Exercises for lectures 13 Design using frequency methods

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering Dynamics and Control II Fall K(s +1)(s +2) G(s) =.

Chapter 12. Feedback Control Characteristics of Feedback 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

Professor Fearing EE C128 / ME C134 Problem Set 7 Solution Fall 2010 Jansen Sheng and Wenjie Chen, UC Berkeley

2.004 Dynamics and Control II Spring 2008

13.42 READING 6: SPECTRUM OF A RANDOM PROCESS 1. STATIONARY AND ERGODIC RANDOM PROCESSES

Control of Manufacturing Processes

Boise State University Department of Electrical Engineering ECE461 Control Systems. Control System Design in the Frequency Domain

Control Systems I. Lecture 9: The Nyquist condition

EE3CL4: Introduction to Linear Control Systems

EE 3CL4: Introduction to Control Systems Lab 4: Lead Compensation

16.30/31, Fall 2010 Recitation # 1

H(s) = s. a 2. H eq (z) = z z. G(s) a 2. G(s) A B. s 2 s(s + a) 2 s(s a) G(s) 1 a 1 a. } = (z s 1)( z. e ) ) (z. (z 1)(z e at )(z e at )

1 (s + 3)(s + 2)(s + a) G(s) = C(s) = K P + K I

Root Locus Design Example #4

Control of Manufacturing Processes

Parametric Functions and Vector Functions (BC Only)

Lecture 5: Frequency domain analysis: Nyquist, Bode Diagrams, second order systems, system types

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

UNIVERSITY OF BOLTON SCHOOL OF ENGINEERING BENG (HONS) IN BIOMEDICAL ENGINEERING SEMESTER 1 EXAMINATION 2017/2018 ADVANCED BIOMECHATRONIC SYSTEMS

Poles and Zeros and Transfer Functions

Chapter Eight. Transfer Functions. 8.1 Frequency Domain Modeling

Transcription:

Systems Analysis and Control Matthew M. Peet Illinois Institute of Technology Lecture 2: Drawing Bode Plots, Part 2

Overview In this Lecture, you will learn: Simple Plots Real Zeros Real Poles Complex Zeros Complex Poles M. Peet Lecture 2: Control Systems 2 / 33

Review Recall: Frequency Response Input: u(t) = M sin(ωt + φ) Output: Magnitude and Phase Shift 2.5 2 1.5 1 Linear Simulation Results y(t) = G(ıω) M sin(ωt + φ + G(ıω)) Amplitude.5.5 1 1.5 2 2 4 6 8 1 12 14 16 18 2 Time (sec) Frequency Response to sin ωt is given by G(ıω) M. Peet Lecture 2: Control Systems 3 / 33

Review Recall: Bode Plot Definition 1. The Bode Plot is a pair of log-log and semi-log plots: 1. Magnitude Plot: 2 log 1 G(ıω) vs. log 1 ω 2. Phase Plot: G(ıω) vs. log 1 ω Bite-Size Chucks: G(ıω) = i (ıωτ zi + 1) i (ıωτ pi + 1) 2 log G(ıω) = i 2 log ıωτ zi + 1 i 2 log ıωτ pi + 1 M. Peet Lecture 2: Control Systems 4 / 33

Plotting Simple Terms Plotting Normal Zeros Behaves as G 1 (s) = (τs + 1) { ıωτ + 1 1 ω << 1 τ = ıωτ ω >> 1 τ A constant at low ω. A zero at high ω. Break Point at ω = 1 τ. Magnitude: db ω << 1 τ 2 log ıωτ + 1 = 3.1 db ω = 1 τ 2 log ω ω >> 1 τ 4 35 3 25 2 15 1 5 1-2 1-1 1 1 1 1 2 4 35 3 25 2 15 ω = τ -1 1 5 1-2 1-1 1 1 1 1 2 M. Peet Lecture 2: Control Systems 5 / 33

Plotting Real Zeros Phase Geometry: (ıωτ + 1) = tan 1 ( ωτ 1 ω << 1 τ (ıωτ + 1) = 45 ω = 1 τ 9 ω >> 1 τ ) 1+ω 2 τ 2 1 Im(s) } ωτ } Re(s) 9 9 45 45 ω = τ -1 1-2 1-1 1 1 1 1 2 1-2 1-1 1 1 1 1 2 M. Peet Lecture 2: Control Systems 6 / 33

Plotting Real Zeros Phase 9 45 1-2 1-1 1 1 1 1 2 We can improve our sketch using a line (ıωτ + 1) ω <.1 τ = 45.1 (1 + log ωτ) τ < ω < 1 τ 9 ω > 1 τ 9 45 ω =.1 τ -1 ω = τ -1 ω = 1 τ -1 1-2 1-1 1 1 1 1 2 M. Peet Lecture 2: Control Systems 7 / 33

Plotting Real Poles Magnitude We can plot real zeros Now lets do the poles. Poles are the opposite of zeros G zero = G 1 pole 1 G pole (ıω) = ıωτ + 1 G zero (ıω) = ıωτ + 1 Magnitude: G pole (ıω) = 1 ıωτ + 1 2 log G pole (ıω) = 2 log ıωτ + 1 4 3 2 1-1 -2-3 -4 1-2 1-1 1 1 1 1 2 M. Peet Lecture 2: Control Systems 8 / 33

Plotting Real Poles Magnitude Our approximation is also a reflection of the zero 2 log 1 ıωτ + 1 = db ω << 1 τ 3.1 db ω = 1 τ 2 log ω ω >> 1 τ 5-5 -1-15 -2-25 -3 ω = τ -1-35 -4 1-2 1-1 1 1 1 1 2 5-5 ω = τ -1-1 -15-2 -25-3 -35-4 1-2 1-1 1 1 1 1 2 M. Peet Lecture 2: Control Systems 9 / 33

Plotting Real Poles Phase We can plot real zeros Now lets do the poles. Poles are the opposite of zeros G pole (ıω) = 1 ıωτ + 1 = G 1 zero(ıω) Phase: G pole (ıω) = tan 1 ωτ = G zero (ıω) 9 45-45 -9 1-2 1-1 1 1 1 1 2 M. Peet Lecture 2: Control Systems 1 / 33

Plotting Real Poles Phase Our approximation is also a reflection of the zero (ıωτ + 1) 1 ω <.1 τ = 45.1 (1 + log ωτ) τ < ω < 1 τ 9 ω > 1 τ -45-9 ω =.1 τ -1 ω = τ -1 ω = 1 τ -1 1-2 1-1 1 1 1 1 2 ω =.1 τ -1-45 ω = τ -1-9 ω = 1 τ -1 1-2 1-1 1 1 1 1 2 M. Peet Lecture 2: Control Systems 11 / 33

Plotting Real Poles and Zeros Lead-Lag We now have the pieces for a Lead-Lag Compensator. Example: Lead Compensator First Step: Standard Form G(s) = s + 1 s + 1 G(s) = s + 1 s + 1 = 1 s + 1 1 1 1 s + 1 Second Step: Plot the Constant c =.1 = 2dB. -6 9 2 1 45-45 -9-1 -2-3 -4 1-2 1-1 1 1 1 1 2 1 3-5 -6 M. Peet Lecture 2: Control Systems 12 / 33

Plotting Real Poles and Zeros Lead-Lag Third Step: Plot the Zero at ω = 1. Magnitude: db ω << 1 2 log ıω + 1 = 3.1 db ω = 1 2 log ω ω >> 1 4 3 2 1-1 ω = 1 ω = 1-2 -3 Phase: -4 1-2 1-1 1 1 1 1 2 (ıω + 1) ω <.1 = 45 (1 + log ω).1 < ω < 1 9 ω > 1 9 45-45 -9 ω = 1 zero 1-2 1-1 1 1 1 1 2 M. Peet Lecture 2: Control Systems 13 / 33

Plotting Real Poles and Zeros Lead-Lag Fourth Step: Plot the Pole at ω = 1. Magnitude: 2 log ıω 1 + 1 1 db ω << 1 = 3.1 db ω = 1 2 log ω ω >> 1 4 3 2 1-1 ω = 1 ω = 1-2 -3 Phase: ( ıω ) 1 1 + 1 ω < 1 = 45 ( ) 1 + log ω 1 1 < ω < 1 9 ω > 1-4 1-2 1-1 1 1 1 1 2 9 45 ω = 1 ω = 1-45 pole -9 1-2 1-1 1 1 1 1 2 M. Peet Lecture 2: Control Systems 14 / 33

Plotting Real Poles and Zeros Lead-Lag Add them all together 2 1 ω = 1-1 -2-3 ω = 1-4 -5-6 1-2 1-1 1 1 1 1 2 9 zero 45-45 ω = 1 ω = 1 total pole -9 1-2 1-1 1 1 1 1 2 M. Peet Lecture 2: Control Systems 15 / 33

Plotting Real Poles and Zeros Lead-Lag Compare with the True Plot 2 1 ω = 1-1 -2-3 True ω = 1-4 -5-6 1-2 1-1 1 1 1 1 2 9 45 True -45 ω = 1 ω = 1-9 1-2 1-1 1 1 1 1 2 M. Peet Lecture 2: Control Systems 16 / 33

Plotting Real Poles and Zeros Lead-Lag Example: Lag Compensator G(s) = 1 s + 1 1 s + 1 = s 1 + 1 s 1 + 1 Pole at ω = 1 Zero at ω = 1 Output Phase Lags behind the input. M. Peet Lecture 2: Control Systems 17 / 33

Plotting Real Poles and Zeros Lead-Lag Put them together to get a Lead-Lag Compensator. Zero at ω = 1 Pole at ω = 1 Pole at ω = 1 Zero at ω = 1 Magnitude: Changes in Slope. ω = 1: +2 ω = 1: -2 ω = 1: -2 ω = 1: +2 G(s) = s + 1 s + 1 s + 1 s + 1 = s + 1 s 4 3 2 1 1 + 1 s 1 + 1 s 1 + 1-1 1-2 1-1 1 1 1 1 2 1 3 1 4 1 5 M. Peet Lecture 2: Control Systems 18 / 33

Plotting Real Poles and Zeros Lead-Lag G(s) = s + 1 s + 1 s + 1 s + 1 Phase: Changes in Slope. ω =.1: +45 ω = 1: 45 ω = 1: 9 ω = 1: +9 ω = 1: +45 ω = 1: 45 1 9 45-45 -9-2 1-1 1 1 1 1 2 1 3 1 4 1 5 M. Peet Lecture 2: Control Systems 19 / 33

Complex Poles and Zeros Now we move on to the final topic: Complex Poles and Zeros Complex Zeros: G(s) = s 2 + 2ζω n s + ω 2 n First Step: Put in standard form ( ( s G(s) = ωn 2 ω n ) ) 2 + 2ζ s + 1 ω n So y ss = 1 Ignoring the constant G 1 (ıω) = ( ω ω n ) 2 + 2ζ ω ω n ı + 1 = 2ζı ( 1 ω ) 2 ω ω n ω n << 1 ω ω n = 1 ω ω n >> 1 M. Peet Lecture 2: Control Systems 2 / 33

Complex Poles and Zeros Magnitude: ( ) 2 ω 2 log G 1 (ıω) = 2 log + 2ζ sω ı + 1 ω n ω n ω db ω n << 1 ω = 2 log(2ζ) ω n = 1 ω 4(log ω log ω n ) ω n >> 1 8 6 4 2-2 1-2 1-1 1 1 1 1 2 M. Peet Lecture 2: Control Systems 21 / 33

Complex Poles and Zeros The Behavior near ω = ω n depends on ζ. When ζ =, 2 log G 1 (ıω) = 2 log(2ζ) = When ζ =.5, When ζ = 1, 2 log G 1 (ıω) = 2 log 1 = 2 log G 1 (ıω) = 2 log 2 = 6.2dB M. Peet Lecture 2: Control Systems 22 / 33

Complex Poles and Zeros Phase: ( ( ) ) 2 ω G 1 (ıω) = 1 + 2ζ ω ı ω n ω n ( ) = tan 1 2ζωωn ωn 2 ω 2 ω ω n << 1 = 9 ω ω n = 1 18 ω ω n >> 1 Re(s) Im(s) } 2 ζ ω/ω n < G(iω) } 1-(ω/ω n ) 2 18 135 9 45 1-2 1-1 1 1 1 1 2 M. Peet Lecture 2: Control Systems 23 / 33

Complex Poles and Zeros The Behavior near ω = ω n depends on ζ. For ω = ω n = 1, When ζ =, d G 1(ıω) = dω When ζ =.1, d G 1(ıω) dω When ζ =.5, d G 1(ıω) dω When ζ = 1, d G 1(ıω) dω = 1 = 2 = 1 M. Peet Lecture 2: Control Systems 24 / 33

Complex Poles and Zeros Comparison vs. True for ω n = 1, ζ =.5. 8 6 4 2-2 1-2 1-1 1 1 1 1 2 18 135 9 45 1-2 1-1 1 1 1 1 2 Good for magnitude, bad for phase M. Peet Lecture 2: Control Systems 25 / 33

Complex Poles and Zeros Comparison vs. True for ω n = 1, ζ =.1. 8 6 4 2-2 1-2 1-1 1 1 1 1 2 18 135 9 45 1-2 1-1 1 1 1 1 2 M. Peet Lecture 2: Control Systems 26 / 33

Complex Poles Treatment of Complex Poles is similar Complex Poles: 1 G(s) = ( ( ) ) 2 s ω n + 2ζ s ω n + 1 2-2 -4-6 -8 Magnitude: 1-2 1-1 1 1 1 1 2 ( ) 2 ω 2 log G 1 (ıω) = 2 log 1 + 2ζ sω ı ω n ω n ω db ω n << 1 ω = 2 log(2ζ) ω n = 1 ω 4(log ω log ω n ) ω n >> 1 M. Peet Lecture 2: Control Systems 27 / 33

Complex Poles The Behavior near ω = ω n depends on ζ. When ζ =, 2 log G 1 (ıω) = 2 log(2ζ) = + When ζ =.5, When ζ = 1, 2 log G 1 (ıω) = 2 log 1 = 2 log G 1 (ıω) = 2 log 2 = 6.2dB M. Peet Lecture 2: Control Systems 28 / 33

Complex Poles Phase: ( ( ω G 1 (ıω) = 1 ω n ω ω n << 1 = 9 ω ω n = 1 18 ω ω n >> 1 ) ) 2 + 2ζ ω ı ω n -45-9 -135-18 1-2 1-1 1 1 1 1 2 M. Peet Lecture 2: Control Systems 29 / 33

Complex Poles The Behavior near ω = ω n depends on ζ. For ω = ω n = 1, When ζ =, d G 1(ıω) = dω When ζ =.1, d G 1(ıω) dω When ζ =.5, d G 1(ıω) dω When ζ = 1, d G 1(ıω) dω = 1 = 2 = 1 M. Peet Lecture 2: Control Systems 3 / 33

Complex Poles Comparison vs. True for ω n = 1, ζ =.5. 2-2 -4-6 -8 1-2 1-1 1 1 1 1 2-45 -9-135 -18 1-2 1-1 1 1 1 1 2 Good for magnitude, bad for phase M. Peet Lecture 2: Control Systems 31 / 33

Complex Poles Comparison vs. True for ω n = 1, ζ =.1. 2-2 -4-6 -8 1-2 1-1 1 1 1 1 2-45 -9-135 -18 1-2 1-1 1 1 1 1 2 Phase improves, Magnitude gets worse. M. Peet Lecture 2: Control Systems 32 / 33

Summary What have we learned today? Simple Plots Real Zeros Real Poles Complex Zeros Complex Poles Next Lecture: Compensation in the Frequency Domain M. Peet Lecture 2: Control Systems 33 / 33