chapter 5: frequency domain analysis of feedback systems

62
Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering, Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020 Chapter 5: Frequency Domain Analysis of Feedback Systems In this chapter we review absolute stability of a linear feedback system with a nonlinear block through circle and Popov criterion. Then by quasi-linear approximation of nonlinear feedback systems into a linear and a nonlinear block, existence, stability and frequency and amplitude of limit cycles are analyzed, by using describing function analysis. Nonlinear Control

Upload: others

Post on 13-Mar-2022

5 views

Category:

Documents


0 download

TRANSCRIPT

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Chapter 5: Frequency Domain Analysis of Feedback Systems

In this chapter we review absolute stability of a linear feedback system with a nonlinear block through circle and Popov criterion. Then by quasi-linear approximation of nonlinear feedback systems into a linear and a nonlinear block, existence, stability and frequency and amplitude of limit cycles are analyzed, by using describing function analysis.

Nonlinear Control

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

WelcomeTo Your Prospect Skills

On Nonlinear System Analysis and

Nonlinear Controller Design . . .

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

About ARAS

ARAS Research group originated in 1997 and is proud of its 22+ years of brilliant background, and its contributions to

the advancement of academic education and research in the field of Dynamical System Analysis and Control in the

robotics application.ARAS are well represented by the industrial engineers, researchers, and scientific figures graduated

from this group, and numerous industrial and R&D projects being conducted in this group. The main asset of our

research group is its human resources devoted all their time and effort to the advancement of science and technology.

One of our main objectives is to use these potentials to extend our educational and industrial collaborations at both

national and international levels. In order to accomplish that, our mission is to enhance the breadth and enrich the

quality of our education and research in a dynamic environment.

Get more global exposure Self confidence is the first key

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Contents

In this chapter we review absolute stability of a linear feedback system with a nonlinear block. Circle and Popov criterion in single variable, andmultivariable systems are described in detail. Then by quasi-linear approximation of nonlinear feedback systems into a linear and a nonlinear block, existence, stability and frequency and amplitude of limit cycles are analyzed, by using describing function analysis.

4

Absolute StabilityIntroduction, definitions, sector nonlinearity, Lure’s problem, Multivariable and single variable circle and Popov criteria.1

Describing Function MethodIllustrating example, assumptions and definitions, computing describing functions for common nonlinearities.2

Describing Function MethodReview of Nyquist criterion, Existence and stability of limit cycles, examples.3

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Frequency Domain Analysis of Feedback Systems5

β€’ Introduction

Feedback Connection representation of many nonlinear

physical systems

β€’ Separation of linear and nonlinear elements.

– Linear + one nonlinear block

– Genuinely nonlinear but separable

β€’ Assume π‘Ÿ = 0, stability of unforced system

– Linear Control theory and frequency response.

– Extension of Nyquist stability criteria

β€’ Absolute Stability

– The origin is globally uniformly asymptotically Stable.

– The Circle and Popov criteria

Feedback system 5.1

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Frequency Domain Analysis of Feedback Systems6

β€’ Absolute Stability

Consider the feedback connected system by closed-loop non-autonomous system

ሢπ‘₯ = 𝐴π‘₯ + 𝐡𝑒 (7.1)

𝑦 = 𝐢π‘₯ + 𝐷𝑒 (7.2)

𝑒 = βˆ’πœ“(𝑑, 𝑦) (7.3)

β€’ where (𝐴, 𝐡) is controllable, (𝐴, 𝐢) is observable,

β€’ And πœ“ is a memoryless, possibly time varying nonlinearity

– Piecewise continuous in 𝑑 and locally Lipschitz in 𝑦

β€’ The transfer matrix of the system:

𝐺 𝑠 = 𝐢 𝑠𝐼 βˆ’ 𝐴 βˆ’1𝐡 + 𝐷 (7.4)

β€’ For all nonlinearities, origin is the eq. point.

β€’ Lure’s problem: Study the stability for sector type nonlinearity.

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Absolute Stability7

β€’ Definitions

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Absolute Stability8

β€’ Definitions

Absolute Stability

β€’ The stability is examined using two Ly. Functions

– Weighted norm Lyapunov function:

– Lure type Lyapunov function:

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Absolute Stability9

β€’ Multivariable Circle Criterion

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Absolute Stability10

β€’ Circle Criterion:

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Absolute Stability11

β€’ Circle Criterion:

Graphical Representation

β€’ Disk 𝐷(𝛼 , 𝛽) is shown graphically

β€’ Extension of Nyquist criteria by replacement of βˆ’1/π‘˜ to the disk 𝐷(𝛼 , 𝛽).

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Absolute Stability12

β€’ Circle Criterion:

Example 1: Consider the system

β€’ Sector Nonlinearity [𝛼, 𝛽]

– Case (a): 𝛼 < 0 < 𝛽

– For absolute stability the Nyquist must lie inside circle:

– 1st choice:

𝐷(𝛼 , 𝛽) = [βˆ’1/4 ,+1/4] = [βˆ’0.25,+0.25]

β€’ This is not the largest sector.

– 2nd choice: dotted circle:

– D(Ξ± ,Ξ²) = [-1/4.4 , +1/1.4]

= [-0.227, +0.714]

β€’ Better than the 1st choice.

4-1.4

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Absolute Stability13

β€’ Circle Criterion:

Example 1: (Cont.)

β€’ Sector Nonlinearity [𝛼, 𝛽]

– Case (b): 0 = 𝛼 < 𝛽

– For absolute stability the Nyquist must lie to the right of line Re(s) -

0.875:

– 𝐷(𝛼 , 𝛽) = [0 ,+1/0.875]

= [0 , 1.117]

– This is not the largest sector.

-0.875

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Absolute Stability14

β€’ Circle Criterion: Example 2: Consider the system

β€’ Linear system + saturation

β€’ Saturation lies in the sector [0 , 1]:

β€’ According to Ex 1 case b, the system is absolutely stable.

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Absolute Stability15

β€’ Circle Criterion:

Example 3: Let

β€’ Not Hurwitz

– Nyquist plot must not enter the disk, but encircle it once in c.c.w.

direction.

– Plot Nyquist.

– Plot circle:

Center: (0,βˆ’3.2) Radius: 0.1688

– The disk is 𝐷( βˆ’3.36 ,βˆ’3.03)

– The stable sector is

[0.297, 0.330]

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Absolute Stability16

β€’ Circle Criterion:

Example 4: Let

β€’ Nonlinearity is saturation belonging to sector [0, 1]

– Not Hurwitz, Nyquist plot must not enter the

disk, but encircle it once in c.c.w. direction.

– Plot Nyquist, and the circle:

Center ∢ (0, βˆ’1),Radius ∢ 0.45

– The disk is 𝐷(βˆ’0.55,βˆ’1.45)

– The stable sector is [0.690, 1.818]

– The system is not globally asymptotically

stable, and it is only absolutely stable in a

finite domain.

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Absolute Stability17

β€’ Popov Criteria:

Consider the feedback connected system by closed-loop non-autonomous system

ሢπ‘₯ = 𝐴π‘₯ + 𝐡𝑒 (7.13)

𝑦 = 𝐢π‘₯ (7.14)

𝑒𝑖 = βˆ’πœ“π‘– 𝑦𝑖 , 1 ≀ 𝑖 ≀ 𝑝 (7.15)

β€’ where (𝐴, 𝐡) is controllable, (𝐴, 𝐢) is observable,

β€’ And πœ“π‘– is a memoryless nonlinearity belonging to sector [0, π‘˜π‘–]

β€’ The transfer matrix of the system is strictly proper:

𝐺 𝑠 = 𝐢 𝑠𝐼 βˆ’ 𝐴 βˆ’1𝐡

β€’ With the Lure-type Lyapunov function absolute stability of the system is proved.

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Absolute Stability18

β€’ Multivariable Popov Criteria:

β€’ For 𝑀 + 𝐼 + 𝑠Γ 𝐺 𝑠 be strictly positive real 𝐺 𝑠 must be Hurwitz.

For scalar case the following theorem may be stated.

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Absolute Stability19

β€’ Popov Criteria:

β€’ Applicable only for SISO and stable systems.

β€’ Only Sufficient Condition.

β€’ Graphical representation of condition (7.19)

Theorem 7.4: Consider scalar nonlinearity πœ™ and the feedback system (5.1) that satisfies the conditions:

β€’ The matrix A is Hurwitz and the pair (A,B) is controllable.β€’ The nonlinearity πœ“ belongs to the sector [0, π‘˜]β€’ There exists a positive number 𝛾 such that

Then the origin is globally asymptotically stable.

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Absolute Stability20

Popov Criteria:

Example: Consider the system

β€’ Let πœ“ = β„Ž 𝑦 βˆ’ 𝛼𝑦, 𝛼 > 0 to make the system matrix Hurwitz.

– Assume the sector nonlinearity β„Ž ∈ 𝛼, 𝛽 , 𝛽 > 𝛼 > 0

– Then πœ“ ∈ 0, π‘˜ , k = 𝛽 βˆ’ 𝛼. Then, the Popov condition is

– For all positive values of 𝛼 and π‘˜ this is satisfied by 𝛾 > 1.

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Absolute Stability21

Popov Criteria:

Example: (Cont.)

– Even for π‘˜ = ∞, the popov condition is satisfied for all πœ”

– Hence, the system is absolutely stable for all β„Ž in the sector 𝛼,∞ ,

where 𝛼 can be arbitrary small.

– The Popov plot is shown for 𝛼 = 1.

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Contents

In this chapter we review absolute stability of a linear feedback system with a nonlinear block. Circle and Popov criterion in single variable, andmultivariable systems are described in detail. Then by quasi-linear approximation of nonlinear feedback systems into a linear and a nonlinear block, existence, stability and frequency and amplitude of limit cycles are analyzed, by using describing function analysis.

22

Absolute StabilityIntroduction, definitions, sector nonlinearity, Lure’s problem, Multivariable and single variable circle and Popov criteria.1

Describing Function MethodIllustrating example, assumptions and definitions, computing describing functions for common nonlinearities.2

Describing Function MethodReview of Nyquist criterion, Existence and stability of limit cycles, examples.3

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Frequency Domain Analysis of Feedback Systems23

β€’ Describing Function Method

Illustrating Example: Van der Pol

β€’ System Dynamics: ሷπ‘₯ + 𝛼 π‘₯2 βˆ’ 1 ሢπ‘₯ + π‘₯ = 0

β€’ Assume, a disk type limit cycle π‘₯ 𝑑 = 𝐴𝑠𝑖𝑛(πœ”π‘‘)

– With amplitude 𝐴 and frequency πœ”.

Separate the linear and nonlinear elements, and Integrate it into

feedback system (5.2)

ሷπ‘₯ βˆ’ 𝛼 ሢπ‘₯ + π‘₯ = 𝛼(βˆ’π‘₯2 ሢπ‘₯)

– Consider the nonlinear term 𝑀 π‘₯ = βˆ’π‘₯2 ሢπ‘₯ in the feedback.

𝐺(𝑠)𝑀(β‹…)π‘₯(𝑑)𝑀(𝑑)βˆ’π‘₯(𝑑)π‘Ÿ = 0

βˆ’+

Nonlinearity Linear SystemFeedback System 5.2

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Method24

Illustrating Example: Van der Pol (Cont.)

β€’ Then the linear system is

𝐺 𝑠 =π‘₯

𝑒=

𝛼

𝑠2βˆ’π›Ό 𝑠+1

β€’ Approximate the nonlinear element

– For a limit cycle π‘₯ 𝑑 β‰ˆ 𝐴 sin(πœ”π‘‘) β‡’ ሢπ‘₯ = π΄πœ” cos πœ”π‘‘

𝑀 = βˆ’π‘₯2 ሢπ‘₯ = βˆ’π΄2 sin2(πœ”π‘‘) π΄πœ” cos πœ”π‘‘

=βˆ’1

2𝐴3πœ” 1 βˆ’ cos 2πœ”π‘‘ cos πœ”π‘‘ =

βˆ’1

4𝐴3πœ” cos(πœ”π‘‘) βˆ’ cos(3πœ”π‘‘)

β€’ Approximate the output

𝑀 β‰ˆβˆ’π΄3

4πœ” cos πœ”π‘‘ =

𝐴2

4

𝑑

π‘‘π‘‘βˆ’π΄ sin πœ”π‘‘ =

𝐴2

4(βˆ’π‘ π‘₯)

β€’ Input/Output Quasi-Linear approximation

𝑀

π‘₯β‰ˆπ΄2

4βˆ’π‘ 

– Amplitude dependent

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Method25

Quasi-linear block diagram

β€’ The frequency response estimate:

𝑀 = 𝑁 𝐴,πœ” βˆ’π‘₯ β‡’ 𝑁 𝐴,πœ” =𝐴2

4π‘—πœ”

β€’ Existence of limit cycle: βˆƒ π‘₯ 𝑑 β‰ˆ 𝐴 sin πœ”π‘‘ , βˆ‹

π‘₯ 𝑑 β‰ˆ 𝐴 sin πœ”π‘‘ = 𝐺 π‘—πœ” 𝑁(𝐴, πœ”)(βˆ’π‘₯) β‡’ 1 + 𝐺 π‘—πœ” 𝑁 𝐴,πœ” π‘₯

– For a nontrivial limit cycle: π‘₯ 𝑑 β‰  0 β†’ 1 + 𝐺 π‘—πœ” 𝑁 𝐴,πœ” = 0

– Compare Nyquist criteria: 1 + 𝐺 π‘—πœ” β‹… π‘˜ = 0

𝛼

𝑠2 βˆ’ 𝛼𝑠 + 1

π‘₯(𝑑)𝑀(𝑑)βˆ’π‘₯(𝑑)π‘Ÿ = 0

βˆ’+

𝐴2

4𝑝

Quasi-linear Approximation 𝐺(𝑠)

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Method26

Illustrating Example: Van der Pol (Cont.)

β€’ Graphical Solution

– Intersect 𝐺(π‘—πœ”) and βˆ’1/𝑁(𝐴,πœ”)

β€’ Analytical Solution (for 𝑝 = π‘—πœ”)

2( )

1G j

j

2 2

1 4 4

( , )

j

N A A j A

2

4, 1 2

jj A

A

πœ” = 1

𝐺(π‘—πœ”)

βˆ’1/𝑁(𝐴,πœ”)

1,2p

2

2 2 2

4 0 2

11 1/ 64 ( 4)

A A

A

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Method

Illustrating Example: Van der Pol (Cont.)

β€’ Graphical Solution

The approx. limit cycle

parameters does not

depend on 𝛼

πœ” = 1

βˆ’1/𝑁(𝐴, πœ”)

𝐺(π‘—πœ”), 𝛼 = 2

𝐺(π‘—πœ”), 𝛼 = 1

𝐺(π‘—πœ”), 𝛼 = 4

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Method

Illustrating Example: Van der Pol (Cont.)

β€’ The approximate L.C. parameters does not depend on 𝛼

β€’ The REAL L.C. parameters depends on 𝛼

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Method29

β€’ Assumptions and Definitions

Systems Under D.F. Analysis

β€’ Separation of linear and nonlinear elements.

– Linear + one nonlinear block

– Genuinely nonlinear but separable

β€’ only one time invariant nonlinearity

β€’ The 1st harmonics of nonlinearity pass through the linear part

(𝐺(𝑠) is low pass).

𝐺(𝑠)𝑀(β‹…)π‘₯(𝑑)𝑀(𝑑)βˆ’π‘₯(𝑑)π‘Ÿ = 0

βˆ’+

Nonlinearity Linear SystemFeedback System 5.2

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Method30

β€’ Assumptions and Definitions

Nonlinear function

β€’ For a Limit cycle the output is periodic.

β€’ Use Fourier Series to represent the output.

β€’ For odd nonlinearities π‘Ž0 = 0.

β€’ By 3rd assumption consider just the first harmonics.

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Method31

β€’ Assumptions and Definitions

Nonlinear function (cont.)

– In which

– In complex representation

β€’ Describing Function

– Extension of input/output transfer function:

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Method32

β€’ Computing Describing Functions

Common nonlinearities

β€’ Hardening Spring

– Fourier series for input: :

– Nonlinearity is odd: π‘Ž0 = π‘Ž1 = 0, and

– Therefore,

– Describing Function is:

– Note: D.F. is real and does not depend on πœ”

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Method33

β€’ Computing Describing Functions

Common nonlinearities

β€’ Saturation:

– Consider the input:

– Saturated output:

– In which,

– Nonlinearity is odd: π‘Ž0 = π‘Ž1 = 0, and

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Method34

β€’ Computing Describing Functions

Common nonlinearities

β€’ Saturation:

– Describing Function:

– D.F. Features:

– Graphical Representation

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Method35

β€’ Computing Describing Functions

Common nonlinearities

β€’ Relay: Special case of saturation

– Direct determination:

– Describing Function:

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Method36

β€’ Computing Describing Functions

Common nonlinearities

β€’ Piecewise Linear

– Describing Function

1( , )b

N aa

( )N a

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Method37

β€’ Computing Describing Functions

Common nonlinearities

β€’ Piecewise Linear

– Describing Function: Graphical Representation

β€’ Saturation

– Special case of piecewise linear with

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Method38

β€’ Computing Describing Functions

Common nonlinearities

β€’ Dead-Zone

– Nonlinear function

– In which,

– Nonlinearity is odd: π‘Ž0 = π‘Ž1 = 0, and

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Method39

β€’ Computing Describing Functions

Common nonlinearities

β€’ Dead-Zone

– Describing Function

– Real, A dependent

– Not depending on πœ”

– No phase shift

– Graphical

Representation

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Method40

β€’ Computing Describing Functions

Common nonlinearities

β€’ Backlash

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Method41

β€’ Computing Describing Functions

β€’ Backlash

Nonlinear Function:

In which

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Method42

β€’ Computing Describing Functions

β€’ Backlash

– Describing Function:

– Therefore,

– D.F. Graphical Representation

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Contents

In this chapter we review absolute stability of a linear feedback system with a nonlinear block. Circle and Popov criterion in single variable, andmultivariable systems are described in detail. Then by quasi-linear approximation of nonlinear feedback systems into a linear and a nonlinear block, existence, stability and frequency and amplitude of limit cycles are analyzed, by using describing function analysis.

43

Absolute StabilityIntroduction, definitions, sector nonlinearity, Lure’s problem, Multivariable and single variable circle and Popov criteria.1

Describing Function MethodIllustrating example, assumptions and definitions, computing describing functions for common nonlinearities.2

Describing Function AnalysisReview of Nyquist criterion, Existence and stability of limit cycles, examples.3

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Method44

β€’ Review of Nyquist Criterion

Feedback system

β€’ Loop gain 𝐿(𝑠) = 𝐢(𝑠)𝑃(𝑠)

– Characteristic Equation 1 + 𝐿(𝑠) = 0 or 𝐿(𝑠) = βˆ’1

– Draw the Nyquist plot of 𝐿(π‘—πœ”) for βˆ’ < πœ” < +

– Count the number of c.w. encirclement of the Nyquist plot around point (βˆ’1,0). 𝑁

– Compute the number of unstable poles of the loop gain 𝐿(𝑠). 𝑃

– The number of unstable roots of the characteristic equation (or the number of unstable poles of the closed-loop transfer function) is denoted by 𝑍 and is found from: 𝑍 = 𝑁 + 𝑃

𝑃(𝑠)𝐢(𝑠)+ -

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Analysis45

β€’ Review of Nyquist Criterion

Feedback system

β€’ Loop gain 𝐿(𝑠) = 𝐢(𝑠)𝑃(𝑠)

– Characteristic Equation 1 + π‘˜ 𝐿(𝑠) = 0 or 𝐿(𝑠) = βˆ’1 / π‘˜

– Draw the Nyquist plot of 𝐿(π‘—πœ”) for βˆ’ < πœ” < +

– Count the number of c.w. encirclement of the Nyquist plot around

point (βˆ’1 / π‘˜ , 0). 𝑁

– Compute the number of unstable poles of the loop gain 𝐿(𝑠). 𝑃

– The number of unstable roots of the characteristic equation (or the number of unstable poles of the closed-loop transfer function) is denoted by 𝑍 and is found from: 𝑍 = 𝑁 + 𝑃

𝐿(𝑠)π‘˜+

-

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Analysis46

β€’ Existence of Limit Cycle

Feedback Systemπ‘₯ = βˆ’π‘¦

𝑀 = 𝑁 𝐴,πœ” π‘₯𝑦 = 𝐺 π‘—πœ” 𝑀

β‡’ 𝑦 = 𝐺 π‘—πœ” 𝑁 𝐴,πœ” (βˆ’π‘¦)

β€’ Characteristics Equation: (for 𝑦 0 )

𝐺 π‘—πœ” 𝑁 𝐴,πœ” + 1 = 0 OR 𝐺 π‘—πœ” = βˆ’1

𝑁 𝐴,πœ”

β€’ Solve it Graphically

β€’ Amplitude dependent describing function

– Intersect 𝐺(π‘—πœ”) and βˆ’1/𝑁(𝐴)

– If intersection occurs L.C. exists

– Amplitude of L.C. 𝐴 (on βˆ’1/𝑁(𝐴))

– Frequency of L.C. πœ” (on 𝐺(π‘—πœ”))

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Analysis47

β€’ Existence of Limit Cycle

Graphical Solution

β€’ General Describing Function (𝑁(𝐴,πœ”))

– Plot 𝐺(π‘—πœ”)

– Plot a set of βˆ’1/𝑁(𝐴,πœ”) for different constant frequencies.

– Analyze the intersection points

– If an intersection point in which

The frequencies match L.C. exists

– Amplitude of L.C. 𝐴 (on -1/𝑁(𝐴,πœ”))

– Frequency of L.C. the matched πœ”

(on 𝐺(π‘—πœ”) and on -1/𝑁(𝐴,πœ”))

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Analysis48

β€’ Stability of Limit Cycle

Feedback System

β€’ Characteristics Equation: (for 𝑦 0 )

𝐺 π‘—πœ” 𝑁 𝐴,πœ” + 1 = 0 OR 𝐺 π‘—πœ” = βˆ’1

𝑁 𝐴,πœ”

Compare it with Nyquist criteria:

𝐿 𝑠 π‘˜ + 1 = 0 or 𝐿(𝑠) = βˆ’1 / π‘˜

β€’ Extended Nyquist Criteria

– Replace the point βˆ’1/π‘˜ by βˆ’1/𝑁(𝐴,πœ”) and count encirclements

– Apply the same Nyquist criteria to determine stability

– Note: stability decay of amplitude

– Instability increase of amplitude

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Unstable Region

Describing Function Analysis49

β€’ Stability of Limit Cycle

Example:

β€’ Consider 𝐺(𝑠) is stable 𝑃 = 0

β€’ Suppose two L.C. 𝐿1, 𝐿2β€’ Stability analysis of 𝐿1

– Perturb 𝐿1by increasing its amplitude

β€’ 𝐿1β€² is inside 𝐺(π‘—πœ”) (unstable)

β€’ The amplitude increases

β€’ 𝐿1β€² moves away from 𝐿1 Unstable Limit Cycle

– Perturb 𝐿1by decreasing its amplitude

β€’ 𝐿1" is outside 𝐺(π‘—πœ”) (stable)

β€’ The amplitude decreases

β€’ 𝐿1" moves away from 𝐿1 Unstable Limit Cycle

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Analysis50

'

2L

"

2L

β€’ Stability of Limit Cycle

Example: (cont.)

β€’ Stability analysis of 𝐿2

– Perturb 𝐿2 by increasing its amplitude

β€’ 𝐿2β€² is outside G(jΟ‰) (stable)

β€’ The amplitude decreases

β€’ 𝐿2β€² moves toward 𝐿2 Stable Limit Cycle

– Perturb 𝐿2 by decreasing its amplitude

β€’ 𝐿2β€²β€² is inside G(jΟ‰) (unstable)

β€’ The amplitude increases

β€’ 𝐿2β€²β€² moves toward 𝐿2 Stable Limit Cycle

Unstable Region

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

β€’ Example 1:

Consider 𝐺 𝑠 =1

𝑠(𝑠+1)(𝑠+2)with unit

β€’ Nyquist Plot

𝐺 π‘—πœ” = βˆ’3πœ” βˆ’ 𝑗(2 βˆ’ πœ”2)

9πœ”3 + πœ” 2 βˆ’ πœ”2 2

– @ point A

β„‘ 𝐺 π‘—πœ” = 0 β†’ 2 βˆ’ πœ”2 = 0 β†’ πœ” = 2

𝐴 = 𝐺 𝑗 2 = βˆ’3 2

18 2 + 0= βˆ’

1

6

(a) For unit saturation (π‘˜ = 1)

– 𝑁(𝐴) < 1 No Crossing, No Limit Cycle

Describing Function Analysis51

a) Saturation

b) Relay

βˆ’1/𝑁(𝐴)

𝐺(π‘—πœ”)

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

β€’ Example 1: (Cont.)

(b) For unit Relay (M=1)

– βˆ’ < βˆ’1/𝑁(𝐴) < 0 Crossing @ point A

– Amplitude:βˆ’1

𝑁 𝐴= βˆ’

1

6β†’

𝑁 𝐴 =4

πœ‹π΄= 6 β†’ 𝐴 =

2

3πœ‹

– Frequency: πœ” = 2 and Limit Cycle is stable.

Describing Function Analysis52

βˆ’1/𝑁(𝐴)

𝐺(π‘—πœ”)

Unstable Region

𝐴

𝐿′

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Analysis53

β€’ Example 1: (Cont.)

Simulation verifications:

– Frequency :

– 2𝑇 = 10 β†’ 𝑇 = 5 and πœ” = 2πœ‹/5 = 1.26 (estimated by 1.414)

– Amplitude 𝐴 = 2/3πœ‹ = 0.21 (red circle)

2𝑇 = 10 sec

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Analysis54

β€’ Example 2:

Consider 𝐺 𝑠 =βˆ’π‘ 

𝑠2+0.8𝑠+8with

β€’ Nyquist Plot

– @ point A

(a) For unit saturation (π‘˜ = 1)

– βˆ’1

𝑁 𝐴= βˆ’1.25 β†’ 𝑁 𝐴 = 0.8

– 𝐴 1.455 @ πœ” = 2.83 Stable Limit Cycle

a) Saturation, 1, 1

b) Dead-zone, 1, 0.5

k a

k

2 2

2 2 2

0.8 (8 )( )

0.64 (8 )

jG j

2( ( )) 0 (8 ) 0 2 2G j

( 2) 1.25A G j

βˆ’1/𝑁(𝐴)

𝐺(π‘—πœ”)

βˆ’1 Unstable Region

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Describing Function Analysis55

β€’ Example 2:

Simulation verifications:

β€’ Frequency :

– 7𝑇 = 16 𝑇 = 2.285 and πœ” = 2πœ‹/𝑇 = 2.75 (estimated by 2.83)

– Amplitude 𝐴 = 1.45 (red ellipse)

7𝑇 = 16sec

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

56

Describing Function Analysis

β€’ Example 2: (Cont.)

β€’ (b) For Dead-zone 𝛿 = 1, π‘˜ = 0.5

β€’ 0 < 𝑁 𝐴 < 0.5 𝑁 𝐴 0.8

β€’ No Crossing No Limit Cycle

βˆ’2βˆ’1/𝑁(𝐴)

𝐺(π‘—πœ”)

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

57

Describing Function Analysis

β€’ Example 3:

Consider Raleigh’s Equation ሷπ‘₯ + π‘₯ = πœ–( ሢπ‘₯ βˆ’1

3ሢπ‘₯3)

β€’ Separation ሷπ‘₯ βˆ’ πœ– ሢπ‘₯ + π‘₯ = βˆ’πœ–(1

3ሢπ‘₯3)

– Consider 𝑦 = ሢπ‘₯ and πœ™ 𝑦 =1

3𝑦3

– 𝑠2 βˆ’ πœ–π‘  + 1 π‘₯ = βˆ’πœ–πœ™ 𝑦 = βˆ’πœ–πœ™ 𝑠π‘₯ = βˆ’πœ–π‘ πœ™ π‘₯

– Hence,

β€’ Linear element: 𝐺 𝑠 =π‘₯

πœ™=

πœ–π‘ 

𝑠2βˆ’πœ–π‘ +1

β€’ Nonlinear function: πœ™ 𝑦 =1

3𝑦3

𝐺(𝑠)πœ™(β‹…)π‘₯(𝑑)πœ™(𝑑)βˆ’π‘₯(𝑑)π‘Ÿ = 0

βˆ’+

Nonlinearity Linear System

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

58

Describing Function Analysis

β€’ Example 3:

Raleigh’s Equation (Cont.)

β€’ Describing Function πœ™ 𝑦 =1

3𝑦3

β€’ 𝑁 𝐴 =𝑏1

𝐴=

2

3πœ‹A0πœ‹π΄ sin πœƒ 3 sin πœƒ π‘‘πœƒ =

1

4𝐴2

β€’ Nyquist Plot

𝐺 𝑠 =πœ–π‘ 

𝑠2 βˆ’ πœ–π‘  + 1β†’ 𝐺 π‘—πœ” =

πœ–π‘—πœ” 1 βˆ’ πœ”2 + πœ–π‘—πœ”

1 βˆ’ πœ”2 2 + πœ–2πœ”2

– @ point A

β„‘ 𝐺 π‘—πœ” = 0 β†’ 1 βˆ’ πœ”2 = 0 β†’ πœ” = 1

– 𝐴 = 𝐺 𝑗1 = βˆ’πœ–2

πœ–2= βˆ’1

– Regardless of πœ–

𝐴

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

59

Describing Function Analysis

β€’ Example 3:

Raleigh’s Equation (Cont.)

β€’ Intersect βˆ’1 / 𝑁(𝐴) with βˆ’1

β€’ 𝑁 𝐴 =1

4𝐴2 = 1 β†’ 𝐴 = 2

β€’ A Limit cycle with

– Amplitude = 2

– Frequency = 1

β€’ Stability analysis

– G(s) is unstable 𝑃 = 2

– For stable region 𝑁 = βˆ’2

– direction of βˆ’1/𝑁(π‘Ž)

β€’ The limit cycle is STABLE

-40

-30

-20

-10

0

Magnitu

de (

dB

)

10-1

100

101

-270

-225

-180

-135

-90

Phase (

deg)

Bode Diagram

Frequency (rad/sec)

βˆ’1/𝑁(𝐴)

βˆ’1

N=-2: Stable Region

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

60

Describing Function Analysis

β€’ Example 3: Raleigh’s Equation (Cont.)

Simulation verifications:β€’ Frequency :

– 2𝑇 = 13.34 πœ” = 2πœ‹/6.67 = 0.95 (estimated by 1)

– Amplitude 𝐴 = 2/3πœ‹ = 0.21 (red circle)

πŸπ‘» = πŸπŸ‘. πŸ‘πŸ’ 𝒔𝒆𝒄

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Chapter 5: Frequency Domain Analysis of Feedback Systems

To read more and see the course videos visit our course website:

http://aras.kntu.ac.ir/arascourses/nonlinear-control/

Thank You

Nonlinear Control

Nonlinear Control K. N. Toosi University of Technology, Faculty of Electrical Engineering,

Prof. Hamid D. Taghirad Department of Systems and Control, Advanced Robotics and Automated Systems October 21, 2020

Hamid D. Taghirad has received his B.Sc. degree in mechanical engineering

from Sharif University of Technology, Tehran, Iran, in 1989, his M.Sc. in mechanical

engineering in 1993, and his Ph.D. in electrical engineering in 1997, both

from McGill University, Montreal, Canada. He is currently the University Vice-

Chancellor for Global strategies and International Affairs, Professor and the Director

of the Advanced Robotics and Automated System (ARAS), Department of Systems

and Control, Faculty of Electrical Engineering, K. N. Toosi University of Technology,

Tehran, Iran. He is a senior member of IEEE, and Editorial board of International

Journal of Robotics: Theory and Application, and International Journal of Advanced

Robotic Systems. His research interest is robust and nonlinear control applied to

robotic systems. His publications include five books, and more than 250 papers in

international Journals and conference proceedings.

About Hamid D. Taghirad

Hamid D. TaghiradProfessor