EE3CL4: Response Plotting the Introduction to Linear Control - - PowerPoint PPT Presentation

ee3cl4
SMART_READER_LITE
LIVE PREVIEW

EE3CL4: Response Plotting the Introduction to Linear Control - - PowerPoint PPT Presentation

EE 3CL4, 8 1 / 77 Tim Davidson Transfer functions Frequency EE3CL4: Response Plotting the Introduction to Linear Control Systems freq. resp. Mapping Section 8: Frequency Domain Techniques Contours Nyquists criterion Ex: servo, P


slide-1
SLIDE 1

EE 3CL4, §8 1 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

EE3CL4: Introduction to Linear Control Systems

Section 8: Frequency Domain Techniques Tim Davidson

McMaster University

Winter 2020

slide-2
SLIDE 2

EE 3CL4, §8 2 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Outline

1

Transfer functions

2

Frequency Response

3

Plotting the freq. resp.

4

Mapping Contours

5

Nyquist’s criterion Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

6

Nyquist’s Stability Criterion as a Design Tool Relative Stability Gain margin and Phase margin Relationship to transient response

slide-3
SLIDE 3

EE 3CL4, §8 4 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Transfer Functions: A Quick Review

  • Consider a transfer function

G(s) = K

  • i(s + zi)
  • j(s + pj)
  • Zeros: −zi; Poles: −pj
  • Note that s + zi = s − (−zi),
  • This is the vector from −zi to s
  • Magnitude:

|G(s)| = |K|

  • i |s + zi|
  • j |s + pj| = |K|prod. dist’s from zeros to s
  • prod. dist’s from poles to s
  • Phase:

∠G(s) = ∠K + sum angles from zeros to s − sum angles from poles to s

slide-4
SLIDE 4

EE 3CL4, §8 6 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Frequency Response

  • For a stable, linear, time-invariant (LTI) system, the

steady state response to a sinusoidal input is a sinusoid of the same frequency but possibly different magnitude and different phase

  • Sinusoids are the eigenfunctions of convolution
  • If input is A cos(ω0t + θ)

and steady-state output is B cos(ω0t + φ), then the complex number B/Aej(φ−θ) is called the frequency response of the system at frequency ω0.

slide-5
SLIDE 5

EE 3CL4, §8 7 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Frequency Response, II

  • If a stable LTI system has a transfer function G(s),

then the frequency response at ω0 is G(s)|s=jω0

  • What if the system is unstable?
slide-6
SLIDE 6

EE 3CL4, §8 9 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Plotting the frequency response

  • For each ω, G(jω) is a complex number.
  • How should we plot it?
  • G(jω) =
  • G(jω)
  • ej∠G(jω)

Plot

  • G(jω)
  • versus ω, and ∠G(jω) versus ω
  • Plot 20 log10
  • G(jω)
  • versus log10(ω), and

∠G(jω) versus log10(ω)

  • G(jω) = Re
  • G(jω)
  • + j Im
  • G(jω)
  • Plot the curve
  • Re
  • G(jω)
  • , Im
  • G(jω)
  • n an “x–y” plot
  • Equiv. to curve
  • G(jω)
  • ej∠G(jω) as ω changes (polar plot)
slide-7
SLIDE 7

EE 3CL4, §8 10 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Polar plot, example 1

Let’s consider the example of an RC circuit

  • G(s) = V2(s)

V1(s) = 1 1+sRC

  • G(jω) =

1 1+jω/ω1 , where ω1 = 1/(RC).

  • G(jω) =

1 1+(ω/ω1)2 − j ω/ω1 1+(ω/ω1)2

  • G(jω) =

1

1+(ω/ω1)2 e−j atan(ω/ω1)

slide-8
SLIDE 8

EE 3CL4, §8 11 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Polar plot, example 1

  • G(jω) =

1 1+(ω/ω1)2 − j ω/ω1 1+(ω/ω1)2

  • G(jω) =

1

1+(ω/ω1)2 e−j atan(ω/ω1)

slide-9
SLIDE 9

EE 3CL4, §8 12 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Polar plot, example 2

Consider G(s) =

K s(sτ+1).

  • Poles at origin and s = −1/τ.
  • To use geometric insight to plot polar plot,

rewrite as G(s) =

K/τ s(s+1/τ)

  • Then
  • G(jω)
  • =

K/τ |jω| |jω+1/τ|

and ∠G(jω) = −∠(jω) − ∠(jω + 1/τ)

slide-10
SLIDE 10

EE 3CL4, §8 13 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Polar plot, ex. 2, G(s) =

K/τ s(s+1/τ)

  • When ω → 0+, |G(jω)| → ∞, ∠G(jω) → −90◦ from below

Tricky

  • To get a better feel, write G(jω) =

−Kω2τ ω2+ω4τ 2 − j ωK ω2+ω4τ 2

Hence, as ω → 0+, G(jω) → −Kτ − j∞

  • As ω increases, distances from poles to jω increase.

Hence |G(jω)| decreases

  • As ω increases, angle from pole at −1/τ increases.

Hence ∠G(jω) becomes more negative

slide-11
SLIDE 11

EE 3CL4, §8 14 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Polar plot, ex. 2, G(s) =

K/τ s(s+1/τ)

  • When ω = 1/τ, G(jω) = (K/τ)/
  • (1/τ)(

√ 2/τ)

  • e−j(90◦+45◦)

i.e., G(jω)|ω=1/τ = (Kτ/ √ 2)e−j135◦

  • As ω approaches +∞, both distances from poles get large.

Hence |G(jω)| → 0

  • As ω approaches +∞, angle from −1/τ approaches −90◦

from below. Hence ∠G(jω) approaches −180◦ from below

slide-12
SLIDE 12

EE 3CL4, §8 15 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Polar plot, ex. 2, G(s) =

K/τ s(s+1/τ) Summary

  • As ω → 0+, G(jω) → −Kτ − j∞
  • As ω increases,

|G(jω)| decreases, ∠G(jω) becomes more negative

  • When ω = 1/τ, G(jω) = (K/

√ 2)e−j135◦

  • As ω approaches +∞,

G(jω) approaches zero from angle −180◦

slide-13
SLIDE 13

EE 3CL4, §8 16 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Polar plot, ex. 2, G(s) =

K/τ s(s+1/τ)

slide-14
SLIDE 14

EE 3CL4, §8 17 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Bode Diagrams

  • Bode magnitude plot

20 log10 |G(jω)| against log10 ω

  • Bode phase plot

∠G(jω) against log10 ω

  • In 2CJ4 we developed rules to help sketch these plots
  • In this course we will use these sketches to design

controllers

slide-15
SLIDE 15

EE 3CL4, §8 18 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Sketching Bode Diagrams

  • Consider generic transfer function of LTI system

G(s) = K

i(s + zi) k(s2 + 2ζkωn,ks + ω2 n,k)

sN

j(s + pj) r(s2 + 2ζd,kωnd,rs + ω2 nd,r)

where zi and pj are real.

  • Unfortunately, not in the form that we are used to for

Bode diagrams

  • Divide numerator by

i zi

  • k ω2

n,k

  • Similarly for denominator
  • Then if ˜

K = K

i zi

  • k ω2

n,k/

  • j pj
  • r ω2

nd,r

  • ,

G(s) = ˜ K

i(1 + s/zi) k

  • 1 + 2ζk(s/ωn,k) + (s/ωn,k)2

sN

j(1 + s/pj) r

  • 1 + 2ζd,k(s/ωnd,r) + (s/ωnd,r)2
slide-16
SLIDE 16

EE 3CL4, §8 19 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Sketching Bode Diagrams, II

  • Now, frequency response can be written as:

G(jω) = ˜ K

i(1 + jω/zi)

(jω)N

j(1 + jω/pj)

×

  • k
  • 1 + 2ζk(jω/ωn,k) + (jω/ωn,k)2
  • r
  • 1 + 2ζd,k(jω/ωnd,r) + (jω/ωnd,r)2
  • Four key components:
  • Gain, ˜

K

  • Poles (or zeros) at origin
  • Poles and zeros on real axis
  • Poles and zeros in complex conjugate pairs
  • Each contributes to the Bode Diagram
slide-17
SLIDE 17

EE 3CL4, §8 20 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Bode Magnitude diagram

G(jω) = ˜ K

i(1 + jω/zi)

(jω)N

j(1 + jω/pj)

×

  • k
  • 1 + 2ζk(jω/ωn,k) + (jω/ωn,k)2
  • r
  • 1 + 2ζd,k(jω/ωnd,r) + (jω/ωnd,r)2
  • Bode Magnitude diagram:

20 log10 |G(jω)| against log10 ω

  • 20 log10 |G(jω)| is

Sum of 20 log10 of components of numerator − sum of 20 log10 of components of denominator

slide-18
SLIDE 18

EE 3CL4, §8 21 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Components for magnitude

G(jω) = ˜ K

i(1 + jω/zi)

(jω)N

j(1 + jω/pj)

×

  • k
  • 1 + 2ζk(jω/ωn,k) + (jω/ωn,k)2
  • r
  • 1 + 2ζd,k(jω/ωnd,r) + (jω/ωnd,r)2
  • Poles at origin: slope starts at −20N dB/dec
  • Gain | ˜

K| incorporated in position of that sloping line

  • First order component in numerator:

increase slope by 20 dB/dec at ω = zi

  • First order component in denominator:

decrease slope by 20 dB/dec at ω = pj

  • Second order components:

increase or decrease slope by 40 dB/dec at ω = ωn

slide-19
SLIDE 19

EE 3CL4, §8 22 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Bode Phase Diagram

G(jω) = ˜ K

i(1 + jω/zi)

(jω)N

j(1 + jω/pj)

×

  • k
  • 1 + 2ζk(jω/ωn,k) + (jω/ωn,k)2
  • r
  • 1 + 2ζd,k(jω/ωnd,r) + (jω/ωnd,r)2
  • Bode Phase Diagram

∠G(jω) against log10 ω

  • ∠G(jω) is

Sum of phases of components of numerator − sum of phases of components of denominator

slide-20
SLIDE 20

EE 3CL4, §8 23 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Components

G(jω) = ˜ K

i(1 + jω/zi)

(jω)N

j(1 + jω/pj)

×

  • k
  • 1 + 2ζk(jω/ωn,k) + (jω/ωn,k)2
  • r
  • 1 + 2ζd,k(jω/ωnd,r) + (jω/ωnd,r)2
  • Phase of ˜

K

  • Poles at origin: −N90◦
  • First order component in numerator:

linear phase change of +90◦ over ω ∈ [zi/10, 10zi]

  • First order component in denominator:

linear phase change of −90◦ over ω ∈ [pj/10, 10pj]

  • Second order components:

phase change of ±180◦ around ω = ωn

slide-21
SLIDE 21

EE 3CL4, §8 24 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Graphically

slide-22
SLIDE 22

EE 3CL4, §8 25 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Graphically

slide-23
SLIDE 23

EE 3CL4, §8 26 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Accuracy of Bode Sketches

Isolated first order pole (analogous for zero)

slide-24
SLIDE 24

EE 3CL4, §8 27 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Accuracy of Bode Sketches

Isolated complex conjugate pair of poles

slide-25
SLIDE 25

EE 3CL4, §8 28 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Accuracy of Bode Sketches

Isolated complex conjugate pair of poles

slide-26
SLIDE 26

EE 3CL4, §8 29 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Example

G(jω) = 5(1 + jω/10) jω(1 + jω/2)

  • 1 + 0.6(jω/50) + (jω/50)2
slide-27
SLIDE 27

EE 3CL4, §8 30 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Example

G(jω) = 5(1 + jω/10) jω(1 + jω/2)

  • 1 + 0.6(jω/50) + (jω/50)2
slide-28
SLIDE 28

EE 3CL4, §8 32 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Introduction

  • We have seen techniques that determine stability of a

system:

  • Routh-Hurwitz
  • root locus
  • However, both of them require a model for the plant
  • Today: frequency response techniques
  • Although they work best with a model
  • For an open-loop stable plant, they also work with

measurements

  • Key result: Nyquist’s stability criterion
  • Design implications: Bode techniques based on gain

margin and phase margin

slide-29
SLIDE 29

EE 3CL4, §8 33 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Characteristic equation

  • To determine the stability of the system we need to

examine the characteristic equation: F(s) = 1 + L(s) = 0 where L(s) = Gc(s)G(s)H(s).

  • The key result involves mapping a closed contour of

values of s to a closed contour of values of F(s).

  • We will investigate the idea of mappings first
slide-30
SLIDE 30

EE 3CL4, §8 34 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Simple example

  • Set F(s) = 2s + 1
  • Map the square in the "s-plane" to the contour in the

"F(s)-plane"

slide-31
SLIDE 31

EE 3CL4, §8 35 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Area enclosed

  • How might we define area enclosed by a closed contour?
  • We will be perfectly rigorous, but will go against

mathematical convention

  • Define area enclosed to be that to the right when the contour

is traversed clockwise

  • What you see when moving clockwise with eyes right
  • Sometimes we say that this area is the area “inside” the

clockwise contour

  • Notions of “enclosed” or “inside” will be applied to contours

in the s-plane

slide-32
SLIDE 32

EE 3CL4, §8 36 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Encirclement

  • In the F(s)-plane, we will be interested in the notion of

encirclement of the origin

  • A contour is said to encircle the origin in the clockwise

direction, if the contour completes a 360◦ revolution around the origin in the clockwise direction.

  • A contour is said to encircle the origin in the anti-clockwise

direction, if the contour completes a 360◦ revolution around the origin in the anti-clockwise direction.

  • We will say that an anti-clockwise encirclement is a

“negative” clockwise encirclement

slide-33
SLIDE 33

EE 3CL4, §8 37 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Example with rational F(s)

  • A mapping for F(s) =

s s+2

  • Note that s-plane contour encloses the zero of F(s)
  • How many times does the F(s)-plane contour encircle

the origin in the clockwise direction?

slide-34
SLIDE 34

EE 3CL4, §8 38 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Cauchy’s Theorem

  • Nyquist’s Criterion is based on Cauchy’s Theorem:
  • Consider a rational function F(s)
  • If the clockwise traversal of a contour Γs in the s-plane

encloses Z zeros and P poles of F(s) and does not go through any poles or zeros

  • then the corresponding contour in the F(s)-plane, ΓF

encircles the origin N = Z − P times in the clockwise direction

  • A sketch of the proof later.
  • First, some examples
slide-35
SLIDE 35

EE 3CL4, §8 39 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Example 1

  • A mapping for F(s) =

s s+1/2

  • s-plane contour encloses a zero and a pole
  • Theorem suggests no clockwise encirclements of origin
  • f F(s)-plane
  • This is what we have!
slide-36
SLIDE 36

EE 3CL4, §8 40 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Example 2

  • s-plane contour encloses 3 zeros and a pole
  • Theorem suggests 2 clockwise encirclements of the
  • rigin of the F(s)-plane
slide-37
SLIDE 37

EE 3CL4, §8 41 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Example 3

  • s-plane contour encloses one pole
  • Theorem suggests -1 clockwise encirclements of the
  • rigin of the F(s)-plane
  • That is, one anti-clockwise encirclement
slide-38
SLIDE 38

EE 3CL4, §8 42 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Informal Justification of Cauchy’s Theorem

  • Consider the case of F(s) = (s+z1)(s+z2)

(s+p1)(s+p2)

  • ∠F(s1) = φz1 + φz2 − φp1 − φp2
  • As the contour is traversed the nett contribution from

φz1 is 360 degrees

  • As contour is traversed, the nett contribution from other

angles is 0 degrees

  • Hence, as contour is traversed, ∠F(s) changes by 360
  • degrees. One encirclement!
slide-39
SLIDE 39

EE 3CL4, §8 43 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Informal Justification

  • Extending this to any number of poles and zeros inside

the contour

  • For a closed contour, the change in ∠F(s) is

360Z − 360P

  • Hence F(s) encircles origin Z − P times
slide-40
SLIDE 40

EE 3CL4, §8 45 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Cauchy’s Theorem (Review)

  • Consider a rational function F(s)
  • If the clockwise traversal of a contour Γs in the s-plane

encloses Z zeros and P poles of F(s) and does not go through any poles or zeros

  • then the corresponding contour in the F(s)-plane, ΓF

encircles the origin N = Z − P times in the clockwise direction

slide-41
SLIDE 41

EE 3CL4, §8 46 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Nyquist’s goal

  • Nyquist was concerned about testing for stability
  • How might one use Cauchy Theorem to examine this?
  • Perhaps choose F(s) = 1 + L(s), as this determines

stability

  • Which contour should we use?
slide-42
SLIDE 42

EE 3CL4, §8 47 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Nyquist’s contour

Actually, we have to be careful regarding poles and zeros on the jω-axis, including the origin. Standard approach is to indent contour so that it goes to the right of any such poles or zeros

slide-43
SLIDE 43

EE 3CL4, §8 48 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Modified Nyquist contour

Here’s an example for a model like that of the motor in the lab.

slide-44
SLIDE 44

EE 3CL4, §8 49 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Coarse Applic. of Cauchy

  • Recall that the zeros of F(s) = 1 + L(s) are the poles of

the closed loop

  • Let P denote the number of right half plane poles of

F(s)

  • The number of right half plane zeros of F(s) is N + P,

where N is the number of clockwise encirclements of the origin made by the image of Nyquist’s contour in the F(s) plane.

  • A little difficult to parse.
  • Perhaps we can apply Cauchy’s Theorem in a more

sophisticated way.

slide-45
SLIDE 45

EE 3CL4, §8 50 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Towards Nyquist’s Criterion

  • F(s) = 1 + L(s), where L(s) is the open loop transfer

function

  • Encirclement of the origin in F(s)-plane is the same as

encirclement of −1 in the L(s)-plane

  • This is more convenient, because L(s) is often

factorized, and hence we can easily determine P

  • Now that we are dealing with L(s), P is the number of

right-half plane poles of the open loop transfer function

  • If we handle the remainder of the components of

Cauchy’s theorem carefully we obtain:

slide-46
SLIDE 46

EE 3CL4, §8 51 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Nyquist’s Criterion: Simplified statement

  • Consider a unity feedback system with an open loop transfer

function L(s) = Gc(s)G(s)H(s),with no z’s or p’s on jω-axis

  • Let PL denote the number of poles of L(s) in RHP
  • Consider the Nyquist Contour in the s-plane
  • Let ΓL denote image of Nyquist Contour under L(s)
  • Let NL denote the number of clockwise encirclements that ΓL

makes of the point (−1, 0)

  • Nyquist’s Stability Criterion:

Number of closed-loop poles in RHP = NL + PL

  • Note that for a stable open loop, the closed-loop is stable if

the image of the Nyquist contour does not encircle (−1, 0).

slide-47
SLIDE 47

EE 3CL4, §8 52 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Ex: L(s) =

1000 (s+1)(s+10) (stable)

  • For 0 ≤ ω < ∞:
  • No zeros, two poles.
  • |L(0)| = 1000/(1 × 10) = 100; ∠L(0) = −0 − 0 = 0
  • Distances from poles to jω is increasing;

hence |L(jω)| is decreasing

  • Angles from poles to jω are increasing;

hence ∠L(jω) is decreasing

  • As ω → ∞, |L(jω)| → 0, ∠L(jω) → −180◦
  • Recall that L(−jω) = L(jω)∗
  • Remember to examine the r → ∞ part of the curve
slide-48
SLIDE 48

EE 3CL4, §8 53 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Ex: L(s) =

1000 (s+1)(s+10) (stable)

Note: No encirclements of (−1, 0) = ⇒ closed loop is stable

slide-49
SLIDE 49

EE 3CL4, §8 54 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Nyquist’s Criterion: Refined statement

  • Consider a unity feedback system with an open loop transfer

function L(s) = Gc(s)G(s)H(s),

  • Let PL denote the number of poles of L(s) in open RHP
  • Consider the modified Nyquist Contour in the s-plane

looping to the right of any poles or zeros on the jω-axis

  • Let ΓL denote image of mod. Nyquist Contour under L(s)
  • Let NL denote the number of clockwise encirclements that ΓL

makes of the point (−1, 0)

  • Nyquist’s Stability Criterion:

Number of closed-loop poles in open RHP = NL + PL

  • Now we can handle open-loop poles and zeros on jω-axis
slide-50
SLIDE 50

EE 3CL4, §8 55 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Example: Pole of L(s) at origin

  • Consider

L(s) = K s(τs + 1)

  • Like in servomotor
  • Problem with the original Nyquist contour
  • It goes through a pole!
  • Cauchy’s Theorem does not apply
  • Must modify Nyquist Contour to go around pole
  • Then Nyquist Criterion can be applied
slide-51
SLIDE 51

EE 3CL4, §8 56 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Example: Pole of L(s) at origin

Now three key aspects of the curve

  • Around the origin
  • Positive frequency axis;

remember negative freq. axis yields conjugate

  • At ∞
slide-52
SLIDE 52

EE 3CL4, §8 57 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Around the origin

  • L(s) =

K s(τs+1)

  • Around the origin, s = ǫejφ,

where φ goes from −90◦ to 90◦

  • In the L(s) plane: limǫ→0 L(ǫejφ)
  • This is: limǫ→0

K ǫejφ = limǫ→0 K ǫ e−jφ

slide-53
SLIDE 53

EE 3CL4, §8 58 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Up positive jω-axis

  • For 0 < ω < ∞, L(jω) =

K ω√ 1+ω2τ 2 e−j(90◦+atan(ωτ))

  • For small ω, L(jω) is large with phase −90◦

Actually, as we worked out in a previous lecture, as ω → 0+, L(jω) → −Kτ − j∞

  • For large ω, L(jω) is small with phase −180◦
  • For ω = 1/τ, L(jω) = Kτ/

√ 2 e−j135◦

slide-54
SLIDE 54

EE 3CL4, §8 59 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

For s = rejθ for large r

  • For s = rejθ with large r, and θ from +90◦ to −90◦,
  • limr→∞ L(rejθ) =

K τr 2 e−j2θ

  • How many encirclements of −1 in L(s) plane? None
  • Implies that closed loop is stable for all positive K
  • Consistent with what we know from root locus (Lab. 2)
slide-55
SLIDE 55

EE 3CL4, §8 60 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Example with open loop RHP pole, proportional control

  • Consider G(s) =

1 s(s−1)

  • Essentially the same as plant model for VTOL aircraft

example in root locus section

  • Consider prop. control, Gc(s) = K1, and H(s) = 1.
  • Hence, L(s) =

K1 s(s−1)

  • Observe that L(s) has a pole in RHP; hence PL = 1
slide-56
SLIDE 56

EE 3CL4, §8 61 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

  • Ex. with open loop RHP pole
  • L(s) =

K1 s(s−1). For s = jω and 0 < ω < ∞,

L(jω) = −K1 1 + ω2 + j K1 ω(1 + ω2) = K1 ω √ 1 + ω2 ∠

  • 90◦ + atan(ω)
  • For ω → 0+, L(jω) → −K1 + j∞
  • As ω increases, real and imag. parts decrease,
  • imag. part decreases faster
  • Equiv. magnitude decreases, phase increases
  • For ω → ∞, L(jω) is small with angle +180◦
  • Conjugate for −∞ < ω < 0
  • What about when s = ǫejθ for −90◦ ≤ θ ≤ 90◦?

L(s) = K1

ǫ ∠(−180◦ − θ)

slide-57
SLIDE 57

EE 3CL4, §8 62 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Example with open loop RHP pole

  • Recall PL = 1
  • Number clockwise encirclements of (−1, 0) is 1
  • Hence there are two closed loop poles in the RHP for all

positive values of K1

  • Consistent with root locus analysis
slide-58
SLIDE 58

EE 3CL4, §8 63 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Root locus of L(s) =

1 s(s−1)

slide-59
SLIDE 59

EE 3CL4, §8 64 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Example with open loop RHP pole, PD control

  • G(s) =

1 s(s−1) and H(s) = 1. L(s) = GC(s)G(s).

  • In the VTOL aircraft example, showed that closed-loop can

be stabilized by lead compensation, GC(s) = Kc(s+z)

(s+p)

  • It can also be stabilized by PD comp., GC(s) = K1(1 + K2s).

(Under the presumption that this can be realized. It can be realized when we have “velocity” feedback.)

  • Using the root locus, we can show that when K2 > 0 there is

a K1 > 0 that stabilizes the closed loop (see next page)

  • Can we see that in the Nyquist diagram?
  • Plot the Nyquist diagram of L(s) = Gc(s)G(s), where

G(s) =

K1 s(s−1) and Gc(s) = 1 + K2s

slide-60
SLIDE 60

EE 3CL4, §8 65 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Root locus analysis

Root locus of (1 + K2s)

1 s(s−1) for a given K2 > 0

  • Poles, zero and active sections of real axis
  • Complete root locus

Conclusion: For any given K2 > 0 there is a ¯ K1 > 0 such that closed loop is stable for all K1 > ¯

  • K1. We can find ¯

K1 using Routh-Hurwitz

slide-61
SLIDE 61

EE 3CL4, §8 66 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Nyquist diagram of (1 + K2s)

K1 s(s−1)

  • Recall that PL = 1
  • If K1K2 > 1, there is one anti-clockwise encirc. of −1
  • In that case, number closed-loop poles in RHP is

−1 + 1 = 0 and the closed loop is stable

  • Consistent with root locus analysis;

but gives ¯ K1 = 1/K2 directly

slide-62
SLIDE 62

EE 3CL4, §8 67 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

One more example

L(s) = K(s − 2) (s + 1)2 Open loop is stable, but has non-minimum phase (RHP) zero L(jω) = K √ ω2 + 4 ω2 + 1 ∠

  • 180◦ − atan(ω/2) − 2 atan(ω)
  • For small positive ω, L(jω) ≈ 2K∠180◦
  • For large positive ω, L(jω) ≈ K

ω ∠ − 90◦

  • In between, phase decreases monotonically, 180◦ → −90◦.

magnitude decreases monotonically (Bode mag dia.)

  • L(jω) =

2K

  • 2ω2−1+jω(5−ω2)
  • (1+ω2)2

; When ω = √ 5, L(jω) = K/2

  • When s = rejθ with r → ∞ and θ : 90◦ → −90◦,

L(s) → (K/r)e−jθ

slide-63
SLIDE 63

EE 3CL4, §8 68 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Nyquist plot of L(s)/K

  • Number of open loop RHP poles: 0
  • Number of clockwise encirclements of −1:

if K < 1/2: 0; if K > 1/2: 1

  • Hence closed loop is

stable for K < 1/2; unstable for K > 1/2

  • This is what we would expect from root locus
slide-64
SLIDE 64

EE 3CL4, §8 69 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Root locus of L(s) =

s−2 (s+1)2

slide-65
SLIDE 65

EE 3CL4, §8 71 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Nyquist’s Criterion (Review)

  • Consider a unity feedback system with an open loop transfer

function L(s) = Gc(s)G(s)H(s),

  • Let PL denote the number of poles of L(s) in open RHP
  • Consider the modified Nyquist Contour in the s-plane

(looping to the right of any poles or zeros on the jω-axis)

  • Let ΓL denote image of mod. Nyquist Contour under L(s)
  • Let NL denote the number of clockwise encirclements that ΓL

makes of the point (−1, 0)

  • Nyquist’s Stability Criterion:

Number of closed-loop poles in open RHP = NL + PL

slide-66
SLIDE 66

EE 3CL4, §8 72 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Relative Stability: Introductory Example

Consider L(s) = K s(τ1s + 1)(τ2s + 1) Nyquist Diagram:

slide-67
SLIDE 67

EE 3CL4, §8 73 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Zoom in

Since L(s) is minimum phase (no RHP zeros), we can zoom in For a given K,

  • how much extra gain would result in instability?

we will call this the gain margin

  • how much extra phase lag would result in instability?

we will call this the phase margin

slide-68
SLIDE 68

EE 3CL4, §8 74 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Formal definitions

  • Gain margin:

1 |L(jωx)|,

where ωx is the frequency at which ∠L(jω) reaches −180◦ amplifying the open-loop transfer function by this amount would result in a marginally stable closed loop

  • Phase margin: 180◦ + ∠L(jωc),

where ωc is the frequency at which |L(jω)| equals 1 adding this much phase lag would result in a marginally stable closed loop

  • These margins can be read from the Bode diagram
slide-69
SLIDE 69

EE 3CL4, §8 75 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Bode diagram

L(jω) = 1 jω(1 + jω)(1 + jω/5)

  • Gain margin ≈ 15 dB
  • Phase margin ≈ 43◦
slide-70
SLIDE 70

EE 3CL4, §8 76 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Phase margin and damping

  • Consider a second-order open loop of the form

L(s) =

ω2

n

s(s+2ζωn), with ζ < 1

  • Closed-loop poles s1, s2 = −ζωn ± jωn
  • 1 − ζ2
  • Let ωc be the frequency at which |L(jω)| = 1
  • Square and rearrange: ω4

c + 4ζ2ω2 nω2 c − ω4 n = 0;

Equivalently, ω2

c

ω2

n =

  • 4ζ4 + 1 − 2ζ2
  • By definition, φpm = 180◦ + ∠L(jωc)
  • Hence

φpm = atan

  • 2
  • (4 + 1/ζ4)1/2 − 2
  • Phase margin is an explicit function of damping ratio!
  • Approximation: for ζ < 0.7, ζ ≈ 0.01φpm, where φpm is

measured in degrees

slide-71
SLIDE 71

EE 3CL4, §8 77 / 77 Tim Davidson Transfer functions Frequency Response Plotting the

  • freq. resp.

Mapping Contours Nyquist’s criterion

Ex: servo, P control Ex: unst., P control Ex: unst., PD contr. Ex: RHP Z, P contr.

Nyquist’s Stability Criterion as a Design Tool

Relative Stability Gain margin and Phase margin Relationship to transient response

Previous example

L(jω) = 1 jω(1 + jω)(1 + jω/5)

  • Phase margin ≈ 43◦
  • Damping ratio ≈ 0.43