chapter 5: frequency domain analysis of feedback systems
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