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

Similar documents
Dynamic circuits: Frequency domain analysis

Core Concepts Review. Orthogonality of Complex Sinusoids Consider two (possibly non-harmonic) complex sinusoids

Lecture 16 FREQUENCY RESPONSE OF SIMPLE CIRCUITS

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

EE221 Circuits II. Chapter 14 Frequency Response

Systems Analysis and Control

MAS107 Control Theory Exam Solutions 2008

Exercise 1 (A Non-minimum Phase System)

EE221 Circuits II. Chapter 14 Frequency Response

Exercise 1 (A Non-minimum Phase System)

Solution: K m = R 1 = 10. From the original circuit, Z L1 = jωl 1 = j10 Ω. For the scaled circuit, L 1 = jk m ωl 1 = j10 10 = j100 Ω, Z L

H(s) = 2(s+10)(s+100) (s+1)(s+1000)

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

Signals and Systems. Lecture 11 Wednesday 22 nd November 2017 DR TANIA STATHAKI

Frequency Response part 2 (I&N Chap 12)

Lab 3: Poles, Zeros, and Time/Frequency Domain Response

Designing Information Devices and Systems II Fall 2018 Elad Alon and Miki Lustig Discussion 5A

ECE382/ME482 Spring 2005 Homework 6 Solution April 17, (s/2 + 1) s(2s + 1)[(s/8) 2 + (s/20) + 1]

Overview of Sampling Topics

Response to a pure sinusoid

Chapter 8: Converter Transfer Functions

If you need more room, use the backs of the pages and indicate that you have done so.

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

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

Frequency Response Analysis

8.1.6 Quadratic pole response: resonance

Control Systems I Lecture 10: System Specifications

Introduction to Feedback Control

Frequency domain analysis

First and Second Order Circuits. Claudio Talarico, Gonzaga University Spring 2015

R-L-C Circuits and Resonant Circuits

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

J. McNames Portland State University ECE 223 Sampling Ver

'XNH8QLYHUVLW\ (GPXQG73UDWW-U6FKRRORI(QJLQHHULQJ. EGR 224 Spring Test II. Michael R. Gustafson II

'XNH8QLYHUVLW\ (GPXQG73UDWW-U6FKRRORI(QJLQHHULQJ. ECE 110 Fall Test II. Michael R. Gustafson II

CHAPTER 7 : BODE PLOTS AND GAIN ADJUSTMENTS COMPENSATION

Sinusoidal steady-state analysis

Refinements to Incremental Transistor Model

Homework 7 - Solutions

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)

'XNH8QLYHUVLW\ (GPXQG73UDWW-U6FKRRORI(QJLQHHULQJ. EGR 224 Spring Test II. Michael R. Gustafson II

Frequency response. Pavel Máša - XE31EO2. XE31EO2 Lecture11. Pavel Máša - XE31EO2 - Frequency response

Lecture 4: R-L-C Circuits and Resonant Circuits

R. W. Erickson. Department of Electrical, Computer, and Energy Engineering University of Colorado, Boulder

Class 13 Frequency domain analysis

STABILITY. Have looked at modeling dynamic systems using differential equations. and used the Laplace transform to help find step and impulse

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

Systems Analysis and Control

FREQUENCY-RESPONSE DESIGN

Overview of Discrete-Time Fourier Transform Topics Handy Equations Handy Limits Orthogonality Defined orthogonal

Step Response Analysis. Frequency Response, Relation Between Model Descriptions

MAE 143B - Homework 9

Fast Fourier Transform Discrete-time windowing Discrete Fourier Transform Relationship to DTFT Relationship to DTFS Zero padding

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

APPLICATIONS FOR ROBOTICS

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

Lecture 9 Infinite Impulse Response Filters

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

Fast Fourier Transform Discrete-time windowing Discrete Fourier Transform Relationship to DTFT Relationship to DTFS Zero padding

MAE 143A: Homework 7 Solutions

Laplace Transform Analysis of Signals and Systems

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

Time Response Analysis (Part II)

Systems Analysis and Control

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

ME 375 Final Examination Thursday, May 7, 2015 SOLUTION

ECSE 4962 Control Systems Design. A Brief Tutorial on Control Design

Sinusoidal Steady-State Analysis

INTRODUCTION TO DIGITAL CONTROL

Problem Weight Score Total 100

EE 4343/ Control System Design Project LECTURE 10

Steady State Frequency Response Using Bode Plots

Control for. Maarten Steinbuch Dept. Mechanical Engineering Control Systems Technology Group TU/e

(b) A unity feedback system is characterized by the transfer function. Design a suitable compensator to meet the following specifications:

VALLIAMMAI ENGINEERING COLLEGE SRM Nagar, Kattankulathur

Frequency Response. Re ve jφ e jωt ( ) where v is the amplitude and φ is the phase of the sinusoidal signal v(t). ve jφ

2.010 Fall 2000 Solution of Homework Assignment 1

Systems Analysis and Control

Transfer func+ons, block diagram algebra, and Bode plots. by Ania- Ariadna Bae+ca CDS Caltech 11/05/15

The Frequency-response Design Method

Frequency Response Techniques

(a) Find the transfer function of the amplifier. Ans.: G(s) =

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

Test II Michael R. Gustafson II

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

Delhi Noida Bhopal Hyderabad Jaipur Lucknow Indore Pune Bhubaneswar Kolkata Patna Web: Ph:

Homework Assignment 3

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

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

MAE 143B - Homework 7

BIOEN 302, Section 3: AC electronics

Frequency Response DR. GYURCSEK ISTVÁN

Transfer function and linearization

Outline. Classical Control. Lecture 1

2.161 Signal Processing: Continuous and Discrete Fall 2008

Compensator Design to Improve Transient Performance Using Root Locus

Chapter 2. Classical Control System Design. Dutch Institute of Systems and Control

Poles and Zeros and Transfer Functions

Assignment 4 Solutions Continuous-Time Fourier Transform

PROBLEMS OF CHAPTER 4: INTRODUCTION TO PASSIVE FILTERS.

Transcription:

Overview of Bode Plots Transfer function review Piece-wise linear approximations First-order terms Second-order terms (complex poles & zeros) J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 1

Transfer Function Review x(t) H(s) y(t) Recall that if H(s) is known and x(t) =A cos(ωt + φ), then we can find the steady-state solution for y(t): y ss (t) =A H(jω) cos (ωt + φ + H(jω)) J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 2

x(t) Bode Plots H(s) y(t) Bode plots are standard method of plotting the magnitude and phase of H(s) Both plots use a logarithmic scale for the x-axis Frequency is in units of radians/second (rad/s) The phase is plotted on a linear scale in degrees Magnitude is plotted on a linear scale in decibels H db (jω) 2 log 1 H(jω) J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 3

Decibel Scales It is important to become adept at translating between amplitude, H(jω), and decibels, H db (jω). Amplitude ( H(jω) ) Decibels (2 log 1 H(jω) ) 1 2 log 1 1 = 1 2 log 1 1 = 1 2 log 1 1 = 1 2 log 1 1 =.1 2 log 1.1 =.1 2 log 1.1 =.1 2 log 1.1 = 1 1 2 log 1 2 = -6.26 2 2 log 1 2 = 1 1 2 2 log 1 2 = J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 4

Example 1: Bode Plots 2 nf 1kΩ 1 kω v s (t) R L + v o (t) - 1. Find the transfer function of the circuit shown above. 2. Generate the bode plot. J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 5

Example 1: Workspace J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 6

Example 1: Bode Plot Active Lowpass RC Filter H(jω) (db) H(jω) (degrees) 2 1 1 1 1 2 1 3 1 4 1 5 18 16 14 12 1 1 1 1 2 1 3 Frequency (rad/s) 1 4 1 5 J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 7

w = logspace(1,5,5); H = -5e3./(j*w + 5e3); subplot(2,1,1); h = semilogx(w,2*log1(abs(h))); set(h, LineWidth,1.4); ylabel( H(j\omega) (db) ); title( Active Lowpass RC Filter ); set(gca, Box, Off ); grid on; set(gca, YLim,[-5 25]); Example 1: MATLAB Code subplot(2,1,2); h = semilogx(w,angle(h)*18/pi); set(h, LineWidth,1.4); ylabel( \angle H(j\omega) (degrees) ); set(gca, Box, Off ); grid on; set(gca, YLim,[85 185]); xlabel( Frequency (rad/s) ); J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 8

Example 2: Bode Plots 2 μf 1kΩ 1 μf 1kΩ v s (t) + R L - 1. Find the transfer function of the circuit shown above. 2. Generate the bode plot. J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 9

Example 2: Workspace J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 1

H(jω) (degrees) H(jω) (db) 8 6 4 2 Example 2: Bode Plot Active Lead/Lag RC Filter 2 1 1 1 2 1 3 1 4 1 5 16 165 17 175 18 (rad/s) 1 1 1 2 1 3 Frequency 1 4 1 5 J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 11

w = logspace(1,5,5); H = -2*(j*w+5)./(j*w + 1); Example 2: MATLAB Code subplot(2,1,1); h = semilogx(w,2*log1(abs(h))); set(h, LineWidth,1.4); ylabel( H(j\omega) (db) ); title( Active Lead/Lag RC Filter ); set(gca, Box, Off ); grid on; set(gca, YLim,[-2 8]); subplot(2,1,2); h = semilogx(w,angle(h)*18/pi); set(h, LineWidth,1.4); ylabel( \angle H(j\omega) (degrees) ); set(gca, Box, Off ); grid on; set(gca, YLim,[-18-16]); xlabel( Frequency (rad/s) ); J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 12

Bode Plot Approximations Until recently (late 198 s) bode plots were drawn by hand There were many rules-of-thumb, tables, and template plots to help Today engineers primarily use MATLAB, or the equivalent Why discuss the old method of plotting by hand? It is still important to understand how the poles, zeros, and gain influence the Bode plot These ideas are used for transfer function synthesis, analog circuit design, and control systems We will discuss simplified methods of generating Bode plots Based on asymptotic approximations J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 13

Alternate Transfer Function Expressions There are many equivalent expressions for transfer functions. H(s) = N(s) D(s) = b ms m + b m 1 s m 1 + + b 1 s + b a n s n + a n 1 s n 1 + + a 1 s + a = b m s ±l (s z 1)(s z 2 )...(s z m ) a n (s p 1 )(s p 2 )...(s p n ) ( )( ) ( ) 1 s = ks ±l z 1 1 s z 2... 1 s z m ( )( ) ( ) 1 s p 1 1 s p 2... 1 s p n This last expression is called standard form The first step in making bode plots is to convert H(s) to standard form J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 14

Magnitude Components Consider the expression for the transfer function magnitude: H db (jω) = 2log 1 H(jω) = 2log 1 k s±l (1 s z 1 )...(1 s (1 s p 1 )...(1 s p n ) z m ) s=jω jω 1 ±l z = 2log 1 k jω 1... 1 jω z m 1 jω p 1... 1 jω p n = 2log 1 k ±l2 log 1 ω +2 log 1 1 jω z 1 + +2log 1 1 jω z m 2 log 1 1 jω p 1 2 log 1 1 jω p n J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 15

Magnitude Components Comments H db (ω) = 2log 1 k ±l2 log 1 ω +2 log 1 1 jω z 1 + +2log 1 1 jω z m 2 log 1 1 jω p 1 2 log 1 1 jω p n Thus, H db (ω) can be written as a sum of simple functions This is similar like using basis functions {δ(t),u(t),& r(t)} to write an expression for a piecewise linear signal We will use this approach to generate our piecewise linear approximations of the bode plot Note that there are four types of components in this expression Constant Linear term Zeros Poles J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 16

Magnitude Components: Constant H(jω) (db) 4 2-2 ω (rad/sec) -4 The constant term, 2 log 1 k, is a straight line on the Bode plot. J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 17

Magnitude Components: Linear Term H(jω) (db) 4 2-2 ω (rad/sec) -4 The linear term, ±l 2 log 1 ω, is a line on the magnitude plot with a slope equal to ±l 2 db per decade. The x-axis intercept occurs at ω =1rad/s. Plot the bode magnitude plots for H(s) =s, 1 s,s2, 1 s 2. J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 18

Magnitude Components: Real Zeros Consider two limiting conditions for a term containing a zero, 2 log 1 1 jω z First condition: ω z lim 2 log ω z 1 1 jω z = Thus, if ω z 1, then2 log 1 Second condition: ω z 1 jω z ω lim z 2 log 1 1 jω z. =2log1 jω z =2log 1 ω 2 log 1 z ω Thus, if z 1, then this term is linear (on a log scale) with a slope of 2 db per decade and an x-axis intercept at ω = z. J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 19

Magnitude Components: Real Zeros Continued H(jω) (db) 4 2-2 ω (rad/sec) -4 Our piecewise approximation joins these two linear asymptotic approximations at ω = z. Plot the piecewise approximation of the term 2 log 1 1 jω z. J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 2

Magnitude Components: Real Zeros Continued 2 5 Bode Magnitude Real Zero: z = ±1 rad/s 4 3 Mag (db) 2 1 1 1 2 1 1 1 Frequency (rad/sec) 1 1 1 2 The approximation is least accurate at ω = z. The true magnitude is 3 db higher than the approximation at this corner frequency. J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 21

Magnitude Components: Real Poles H(jω) (db) 4 2-2 ω (rad/sec) -4 Consider two limiting conditions for a term containing a pole, 2 log 1 1 jω p This is just the negative of the expression for a zero The piecewise approximation is the mirror image of that for a zero J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 22

Magnitude Components: Real Poles Continued 1 Bode Magnitude Real Pole: p = 1 rad/s 1 Mag (db) 2 3 4 5 1 2 1 1 1 Frequency (rad/sec) 1 1 1 2 The approximation is least accurate at ω = p. The true magnitude is 3 db less than the approximation at this corner frequency. J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 23

Complex Poles & Zeros Complex poles and zeros require special attention Will discuss later You will not be expected to plot approximations with complex poles or zeros on exams There are essentially 3 steps to generating piecewise linear approximations of bode plots 1. Convert to standard form 2. Plot the components 3. Graphically add the components together J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 24

Example 3: Magnitude Components Draw the piecewise approximation of the bode magnitude plot for H(s) = (s + 1)(s + 1)2 1s 2 (s + 1) J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 25

Example 3: Workspace H(jω) 6 (db) 4 2-2 -4-6 1 1 1 1 2 1 3 1 4 1 5 ω (rad/sec) J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 26

Example 3: Solution Bode Magnitude Example 1 6 4 Mag (db) 2 2 4 1 1 1 1 1 1 2 1 3 1 4 1 5 Frequency (rad/sec) J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 27

Example 4: Magnitude Components Draw the piecewise approximation of the bode magnitude plot for H(s) = 1 11 s(s + 1) (s + 1)(s + 1)(s +1, ) 2 J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 28

Example 4: Workspace H(jω) 6 (db) 4 2-2 -4-6 1 1 1 1 2 1 3 1 4 1 5 ω (rad/sec) J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 29

Example 4: Solution Mag (db) 8 7 6 5 4 3 2 1 1 2 3 Bode Magnitude Example 2 1 1 2 1 4 1 6 Frequency (rad/sec) J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 3

Phase Components H(jω) = k (jω)l (1 jω z 1 )...(1 jω z m ) (1 jω p 1 )...(1 jω p n ) Each of these terms can be expressed in polar form: a + jb = Ae jθ. Note that (jω) l = ω l j l = ω l (e j π 2 ) l = ω l e j π 2 l H(jω) = k e jη kπ (ωl e jl π 2 )N 1 e jθ 1...N m e jθ m D 1 e jφ 1...Dn e jφ n = k ω l N 1...N m D 1...D n exp ( j(η k π + l π 2 + θ 1 + + θ n φ 1 φ n ) ) where η k = { k 1 k< J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 31

Phase Components: Comments H(jω) = η k π + l π 2 + θ 1 + + θ n φ 1 φ n ( = η k π + l π 2 + 1 jω ) ( + + 1 jω ) z 1 z m ( 1 jω ) ( 1 jω ) p 1 p n Thus the phase of H(jω) is also a linear sum of the phases due to each component We will consider each of the four components in turn Constant Linear term Zeros Poles J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 32

H(jω) 18 9 Phase Components: Constant (deg) 9 ω (rad/sec) 18 The complex angle of the constant term, k, is either if k> or 18 if k<. J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 33

H(jω) 18 9 Phase Components: Linear Term (deg) 9 ω (rad/sec) 18 The linear term, (jω) l = j l = l 9, is a constant multiple of 9 Plot the bode phase plots for H(s) =s, 1 s,s2, 1 s 2. J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 34

Phase Components: Real Zeros Consider three limiting conditions for a term containing a zero, 1 jω z First condition: ω z lim ( 1 jω ) ω z z = ). ω Thus, if z 1, then ( 1 jω z Second condition: ω = z ( ) 1 jω ω= z z = (1 jη z )= η z 45 where η z = sign(z) Third condition: ω z ω lim z ( 1 jω z Thus, if ω z 1, then ( 1 jω z ) = ( jω z ) ηz 9. ) = ηz 9 J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 35

Phase Components: Real Zeros Continued H(jω) 18 9 (deg) 9 ω (rad/sec) 18 Our piecewise approximation joins these three linear asymptotic approximations at ω =1 1 z and ω =1 z. Plot the piecewise approximation of the term ( 1 jω ) z. Assume that z is in the left-hand plane (i.e. Re{z} < ) J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 36

Phase Components: Real Zeros in Left Plane Phase (deg) 1 9 8 7 6 5 4 3 2 1 1 1 2 1 1 1 Frequency (rad/sec) 1 1 1 2 The approximation is least accurate at ω =.1 z and ω =1 z. J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 37

Phase Components: Real Zeros Continued 2 H(jω) 18 9 (deg) 9 ω (rad/sec) 18 If the zero is in the right half plane (i.e. Re{z} > ), then the phase approaches 9 asymptotically. Plot the piecewise approximation of the term ( 1 jω ) z. Assume that z is in the right half plane. J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 38

Phase Components: Real Zeros in Right Plane Phase (deg) 1 1 2 3 4 5 6 7 8 9 1 1 2 1 1 1 Frequency (rad/sec) 1 1 1 2 The approximation is least accurate at ω =.1 z and ω =1 z. J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 39

Phase Components: Real Poles Real poles in the left half plane have the same phase as real zeros in the right half plane. We will only discuss poles in the left half plane because only these systems are stable. First condition: ω p lim ω p ( 1 jω p ) = Second condition: ω = p ( 1 jω p = (1 sign(p) j) = (1 + j) = 45 ) ω= p Third condition: ω p ( lim 1 jω ω p p ) = ( jω p ) = j = 9 J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 4

Phase Components: Real Poles H(jω) 18 9 (deg) 9 ω (rad/sec) 18 Plot the piecewise approximation of the term that p is in the left-hand plane (i.e. Re{p} < ) ( 1 jω p ). Assume J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 41

Phase Components: Real Poles in Left Plane Phase (deg) 1 1 2 3 4 5 6 7 8 9 1 1 2 1 1 1 Frequency (rad/sec) 1 1 1 2 The approximation is least accurate at ω =.1 p and ω =1 p. J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 42

Example 5: Phase Components H(jω) 27 18 9 (deg) 9 18-27 1 1 1 1 2 1 3 1 4 1 5 (rad/sec) Draw the piecewise approximation of the bode phase plot for H(s) = (s + 1)(s + 1)2 1s 2 (s + 1) J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 43

Example 5: Solution 5 Phase (deg) 5 1 15 1 1 1 1 1 1 2 1 3 1 4 1 5 Frequency (rad/sec) J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 44

Example 6: Phase Components H(jω) 27 18 9 (deg) 9 18-27 1 1 1 1 2 1 3 1 4 1 5 (rad/sec) Draw the piecewise approximation of the bode phase plot for H(s) = 1 11 s(s + 1) (s + 1)(s + 1)(s +1, ) 2 J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 45

Example 6: Solution 1 5 Phase (deg) 5 1 15 1 1 2 1 4 1 6 Frequency (rad/sec) J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 46

Example 7: Magnitude H(jω) 6 (db) 4 2-2 -4-6 1 1 1 1 2 1 3 1 4 1 5 ω (rad/sec) Plot the magnitude of H(s) =1 s +1 s + 1 J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 47

Example 7: Phase H(jω) 27 18 9 (deg) 9 18-27 1 1 1 1 2 1 3 1 4 1 5 (rad/sec) Plot the phase of H(s) =1 s +1 s + 1 J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 48

Example 7: Solution 4 Bode Example 1 Mag (db) 2 2 1 1 1 1 1 1 1 2 1 3 1 4 1 5 Phase (deg) 5 1 1 1 1 1 1 2 1 3 1 4 1 5 Frequency (rad/sec) J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 49

Example 8: Magnitude H(jω) 6 (db) 4 2-2 -4-6 1 1 1 1 2 1 3 1 4 1 5 ω (rad/sec) Plot the magnitude of H(s) = 1 s 1 s + 1 J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 5

Example 8: Phase H(jω) 27 18 9 (deg) 9 18-27 1 1 1 1 2 1 3 1 4 1 5 (rad/sec) Plot the phase of H(s) = 1 s 1 s + 1 J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 51

Example 8: Solution 4 Bode Example 2 Mag (db) 2 2 1 1 1 1 1 1 2 1 3 1 4 1 5 Phase (deg) 5 1 15 1 1 1 1 1 1 2 1 3 1 4 1 5 Frequency (rad/sec) J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 52

1 mh 1 mf Example 9: Circuit Example + v s (t) 11 Ω v o (t) - Draw the straight-line approximations of the transfer function for the circuit shown above. Hint: Recall from one of the Transfer Functions Examples R L H(s) = s s 2 + R L s + 1 LC J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 53

Example 9: Magnitude H(jω) 6 (db) 4 2-2 -4-6 1 1 1 1 2 1 3 1 4 1 5 ω (rad/sec) J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 54

Example 9: Phase H(jω) 27 18 9 (deg) 9 18-27 1 1 1 1 2 1 3 1 4 1 5 (rad/sec) J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 55

Example 9: Solution 2 Bode Example 3 Mag (db) 2 4 1 1 1 1 1 1 1 2 1 3 1 4 Phase (deg) 5 5 1 1 1 1 1 1 Frequency 1 2 (rad/sec) 1 3 1 4 J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 56

Complex Poles Complex poles can be expressed in the following form: C(s) = ω 2 n s 2 +2ζω n s + ω 2 n = 1 1+2ζ s ω n + ( ) 2 s ω n ω n is called the undamped natural frequency ζ (zeta) is called the damping ratio The poles are p 1,2 =( ζ ± ζ 2 1) ω n If ζ 1, the poles are real If <ζ<1, the poles are complex If ζ =, the poles are imaginary: p 1,2 = ±jω n If ζ<, the poles are in the right half plane (Re{p} > ) and the system is unstable J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 57

Complex Poles Continued The transfer function C(s) can also be expressed in the following form C(s) = 1 1+2ζ s ω n + ( ) 2 = s ω n 1 1+ s Qω n + ( ) 2 s ω n where Q 1 2ζ The meaning of Q, thequality factor, will become clear in the following slides. 1 C(jω)= 1 ( ) 2 ω ω n + jω Qω n J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 58

Complex Poles Magnitude 2 log 1 C(jω) = 2log 1 1 ( ) 2 ω + jω ω n Qω n 1 ( ) 2 = 2 log 1 1 ω2 ωn 2 + ( ω Qω n ) 2 For ω ω n, 2 log 1 C(jω) 2 log 1 1 =db For ω ω n, 2 log 1 C(jω) 2 log 1 ω 2 ω 2 n = 4 log 1 ω ω n db At these extremes, the behavior is identical to two real poles. J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 59

For ω = ω n, Complex Poles Magnitude Continued ( ) 2 2 log 1 C(jω) = 2 log 1 1 ω2 ωn 2 + ( ω Qω n 2 log 1 C(jω n ) = 2 log 1 1 Q =2log 1 Q = Q db ) 2 J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 6

Mag (db) 4 3 2 1 1 2 3 4 5 Complex Poles Magnitude Complex Poles Q =.1 Q =.5 Q =.77 Q = 1 Q = 2 Q = 1 Q = 1 6 1 1 1 1 1 Frequency (rad/sec) J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 61

For ω ω n, For ω ω n, For ω = ω n, Complex Poles Phase 1 C(jω) = ( ) 2 ω 1 ω n + jω Qω n C(jω) 1 ω Qω n C(jω) 1 = = Qω n ω = 1= 18 C(jω) 1 Qj = 1 j = j = 9 At the extremes, the behavior is identical to two real poles. At other values of ω near ω n, the behavior is more complicated. J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 62

Phase (deg) 2 4 6 8 1 12 14 16 18 Complex Poles Phase Complex Poles Q =.1 Q =.5 Q =.77 Q = 1 Q = 2 Q = 1 Q = 1 1 2 1 1 1 1 1 1 2 Frequency (rad/sec) J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 63

Complex Zeros Left half plane: Inverted magnitude of complex poles Inverted phase of complex poles Right half plane: Inverted magnitude of complex poles Same phase of complex poles This is the same relationship real zeros had to real poles J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 64

Complex Zeros Magnitude 8 Complex Zeros 6 4 Mag (db) 2 2 4 Q =.1 Q =.5 Q =.77 Q = 1 Q = 2 Q = 1 Q = 1 1 1 1 1 1 Frequency (rad/sec) J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 65

Complex Zeros Phase 18 Complex Zeros 16 14 Phase (deg) 12 1 8 6 4 2 Q =.1 Q =.5 Q =.77 Q = 1 Q = 2 Q = 1 Q = 1 1 2 1 1 1 1 1 1 2 Frequency (rad/sec) J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 66

Complex Poles Maximum What is the frequency at which C(jω) is maximized? C(jω) = C(jω) = 1 ( ) ( ) 2 1+ jω Qω n + jω ω n 1 ( ) 2 ( 1 ω2 + ω 2 n ) 2 ω Qω n For high values of Q, the maximum of C(jω) > 1 This is called peaking The largest Q before the onset of peaking is Q = 1 2.77 This curve is said to be maximally flat This is also called a Butterworth response In this case, C(jω n ) = 3 db and ω n is the cutoff frequency J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 67

Complex Poles Maximum Continued If Q>.77, the maximum magnitude and frequency are as follows: ω r = ω n 1 1 Q 2Q 2 C(jω r ) = 1 1 4Q 2 ω r is called the resonant frequency or the damped natural frequency As Q, ω r ω n For sufficiently large Q (say Q>5) ω r ω n C(jω r ) Q Peaked responses are useful in the synthesis of high-order filters Complex zeros (in the left half plane) have the inverted magnitude and phase of complex poles J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 68

5 mh R Complex Poles Example + v s (t) 2 nf v o (t) - Generate the bode plot for the circuit shown above. H(s) = 1 LC s 2 + R L s + 1 LC = ω 2 n s 2 +2ζω n s + ω 2 n = ω 2 n s 2 + ω n Q s + ω 2 n J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 69

Complex Poles Example Continued 1 1 ω n = =1k rad/s LC ζ = R R LC = CL = R.1 2L 2 1 L L Q = LC R = 1 C R = 5 R R =5Ω ζ =.5 Q = 1 Very Light Damping R =5Ω ζ =.5 Q =1 Light Damping R = 77 Ω ζ =1.41 Q =.77 Strong Damping R =1kΩ ζ =1 Q =.5 Critical Damping R =5kΩ ζ =5 Q =.1 Over Damping J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 7

Complex Poles Example Continued 2 4 Resonance Example Mag (db) 2 2 4 6 R = 5 Ω R = 5 Ω 8 R = 77 Ω R = 1 kω 1 R = 5 kω R = 5 kω 12 1 2 1 3 1 4 Frequency (rad/sec) 1 5 1 6 J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 71

Complex Poles Example Continued 3 Resonance Example 2 4 Phase (deg) 6 8 1 12 14 16 18 R = 5 Ω R = 5 Ω R = 77 Ω R = 1 kω R = 5 kω R = 5 kω 1 2 1 3 1 4 1 5 1 6 Frequency (rad/sec) J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 72

Complex Poles Example Continued 4 Output (V) 2 1.8 1.6 1.4 1.2 1.8.6.4.2 Resonance Example R = 5 Ω R = 5 Ω R = 77 Ω R = 1 kω R = 5 kω R = 5 kω.1.2.3.4.5.6.7.8.9.1 Time (sec) J. McNames Portland State University ECE 222 Bode Plots Ver. 1.19 73