Topic 4 Describing Function PDF
Topic 4 Describing Function PDF
Part I
Prediction of Limit Cycles
Two Important Destinations of State Trajectories:
1. Equilibrium Points
2. Limit Cycles
V(e) e Stoichiometric
HEGO Sensor - Values
Switching Sensor
e
Basic Structure
Linear Dynamics
- G( s) Output e
V(e)
N
Time Invariant Nonlinearity
Large e:
Small Gain
Stability Analysis of the Closed-Loop System
K P e s Im
G( s)
Ts 1 Nyquist Plot of G(s)
Re
-1/N -1/N
large e small e ω increase
Re
-1/N -1/N
large e small e ω increase
Intersection of
G(jw) and -1/N(e)
u(t ) a sin(t )
3. In nonlinear systems, a periodic trajectory may contain
higher frequency harmonics, whose frequencies are
multiples of the base frequency (the frequency of the
limit cycle).
u(t ) a0 a1 sin(t ) b1 cos(t ) a2 sin(2t ) b2 cos(2t )
If the linear part G(s) is a low-pass filter, then the higher
frequency components of u(t) will be attenuated after passing
G(s).
Harmonic Approximation
Describing Function of M Linear Dynamics
A sin(t ) u x
Linearized N(A,ω) G(s)
System -
1. The linear part G(s) is low-pass: This will ensure that the base
frequency component is dominant in the closed-loop.
2. The nonlinear part is time invariant: This will ensure that it is
possible to approximate it by a linear time-invariant system.
3. The nonlinearity M is symmetric to the origin: This will ensure
that the DC component a0 = 0 (the signal has zero average).
M(x) M(x)
x x
Describing Function Analysis
Describing Function of M Linear Dynamics
A sin(t ) u x
N(A,ω) G(s)
-
1 G( j) N ( A, ) 0
i.e., for this value A, the closed-loop system has a pole at
jω
Two equations: real part and imaginary part
Two unknowns: A and ω
1 G( j) N ( A, ) 0
x ( x 1) x x 0
2
u x x M ( x)
2
x A sin t
M ( x) x 2 x A2 sin 2 (t ) A cos(t )
A3 A3
cos(t ) cos(3t )
4 4
Retain only the base frequency component
A3 A2 d ( A sin(t )) A2
u (t ) cos(t ) x
4 4 dt 4
The describing function:
A2
N ( A, ) j
4
Try to solve the equations:
1 G( j) N ( A, ) 0
A 2
1 j 0
( j ) ( j ) 1 4
2
j A2 4(1 2 ) j 4
1, A2
a0
c0 (t ) a1 cos(t ) b1 sin(t )
2
1
a0
M ( A sin(t ))d (t )
1
a1
M ( A sin(t )) cos(t )d (t )
1
b1
M ( A sin(t )) sin(t )d (t )
Describing Function:
1 b1
N ( A, ) [b1 ja1 ]
A A
M(e)
Relay (Switching Nonlinearity)
B
M is an odd function.
e
-B
2 4B N(A)
b1
B sin(t )d (t )
0
b1 4 B
N ( A)
A A
A
Saturation M is an odd function.
C(t)
M(e)
ka
a e t1 t
2
e
t1
t
Case 1: Aa It is in the linear range
N ( A) k
Case 2: Aa Saturation is taking effect
kA sin(t ), 0 t t1
c(t )
ka, t1 t / 2
2
b1
c(t ) sin(t )d (t )
0
/2
4
c(t ) sin(t )d (t )
0
4 1
t /2
kA sin (t )d (t ) ka sin(t )d (t )
2
0 t1
2kA a a2
t1 1 2
A A
At the saturation point
a
A sin(t1 ) a t1 sin 1
A
2kA 1 a a a2
b1 sin 1 2
A A A
Describing Function:
b1 2k 1 a a a2
N ( A, ) sin 1 2 N(A)
A A A A
k
a A
Dead Zone M is an odd function.
C(t)
M(e)
e t1 t
δ
2
e
t1
t
Case 1: A It is in the dead zone
N ( A) 0
Case 2: A Linear part is taking effect
0, 0 t t1
c(t )
k ( A sin(t ) ), t1 t / 2
4 /2
b1 k ( A sin(t ) ) sin(t )d (t )
t1
2kA 1 2
sin 1 2
2 A A A
Describing Function:
b1 2k 1 2
N ( A, ) sin 1 2
A 2 A A A
N(A)
δ A
Backlash M is NOT an odd function.
M(e)
k
e
δ
4 k
a1
A
1
kA
2
1 2 2 2
b1 sin 1 1 1 1
2 A A A
Describing Function:
1
N ( A, ) (b1 ja1 )
A
k
2
2 2 2 4 k
sin 1 1 1 1 1 j 1
2 A A A A A
|N(A)|
k
δ A
N ( A)
A
0o
-90 o
C(t)
M(e)
k
e t
δ t1 t2
e
t1
There is a phase delay!
t2
t
Describing Function Analysis
Part II
Main Issues
Time Invariant
Nonlinearity Linear Dynamics
u x
- M G( s)
e
r N(A) G(s) y
-
Nyquist Criterion:
Z NP
Z = the number of unstable zeros of 1+ L(s)
P = the number of unstable poles of L(s)
N = the number of clockwise encirclement of the Nyquist
plot of L around the critical point (-1,0).
For the closed-loop stability, Z = 0.
N P
The Nyquist plot must encircle (-1,0) counter-clockwise P times.
Suppose the open loop system L(s) is stable, P = 0
(-1,0)
(-1/k,0) (-1/k,0)
Closed-loop system
Is stable
Nyquist Plot of G(jω)
-1/N(A)
-1/N(A)
- G( s) Output e
V(e)
N
Time Invariant Nonlinearity
G( j) 1/ N ( A)
G( j) 1/ N ( A)
does not have a solution (A,ω) if their graphs do not intersect.
Im
Nyquist Plot of G(jω)
A increase
Re
ω increase
Plot of -1/N(A)
Re
Intersections of
G(jω) and -1/N(A) ω increase
(A1 ,ω1)
Plot of -1/N(A)
Intersections of G(jω) and -1/N(A): Prediction of two limit cycles:
One with magnitude A1 and frequency ω1
the other with magnitude A2 and frequency ω2
Stability Analysis of the Limit Cycles
Im
A increase
Nyquist Plot of G(jω)
(A2 ,ω2)
Re
ω increase
(A1 ,ω1)
Re
ω increase
(A1 ,ω1)
u 10s
M x
- s 2 2.1s 100
M(e)
Relay (Switching Nonlinearity)
B
e
2 4B
b1
B sin(t )d (t )
0
-B
N(A)
b1 4 B
N ( A)
A A
A
Nyquist plot of G(jω)
1
N ( A) (-4.76,0)
Magnitude: A 6.06B
10 j
G( j ) 4.76 j 0
100 j 2.1
2
Frequency: 10
1
N ( A)
(-4.76,0)
a e t1 t
2
e
t1
t
Describing Function:
b1 2k 1 a a a2
N ( A, ) sin 1 2
A A A A
N(A)
a A
Nyquist plot of G(jω)
1
N ( A) -1/k (-4.76,0)
1
4.76 k 0.21
k
1
N ( A) (-4.76,0) -1/k
Frequency: 10
N(A)
k
1
4.76 N ( A) 0.21 0.21
N ( A)
A0 A
Magnitude: A A0
Nyquist plot of G(jω)
1
N ( A)
(-4.76,0)
Linear Dynamics
Nonlinearity
u s
M x
- s 2 3s 2
Nyquist plot of G(jω)
(-1/3,0)
j
G( j )
2 2 j 3
At the intersection: 1
2, G( j 2)
3
δ
e t1 t
2
e
t1
t
Describing Function:
b1 2k 1 2
N ( A, ) sin 1 2
A 2 A A A
N(A)
δ A
1 Nyquist plot of G(jω)
N ( A) (-1/3,0)
-1/k
-1/k
1 1 k
N ( A) 3
N ( A) 3
3
δ A0 A
Magnitude: A A0
Case 3: Saturation Nonlinearity
M(e) C(t)
ka
a e t1 t
2
e
t1
t
Describing Function:
b1 2k 1 a a a2
N ( A, ) sin 1 2
A A A A
N(A)
a A
1 Nyquist plot of G(jω)
N ( A) (-1/3,0)
-1/k
-1/k
x ( x 1) x x 0
2
u x x M ( x)
2
x A sin t
M ( x) x 2 x A2 sin 2 (t ) A cos(t )
A3 A3
cos(t ) cos(3t )
4 4
Retain only the base frequency component
A3 A2 d ( A sin(t )) A2
u (t ) cos(t ) x
4 4 dt 4
The describing function:
A2
N ( A, ) j It is a function of
4 both A and ω
Equivalent System for Analysis on Limit Cycles:
A2
N ( A, ) G ( j ) j
4 ( j ) 2 j 1
A2 j
N ( A) G ( j )
4 ( j ) j 1
2
2 u s
A x
- 4 s2 s 1
j
1 j G( j ) 1
( j ) j 1
2
Nyquist plot of the equivalent system
G( j ) j G( j )
1 4
2
N ( A) A (-1,0)
j
Frequency: G ( j )
1 2 ( j )
1, G ( j1) 1
1 4
Magnitude: 1 2 1 A 2
N ( A) A