Lec03 11 PDF

EL2620 2011 EL2620 2011
Equilibrium Points for Linear Systems

x2 Phase-Plane Analysis for Nonlinear Systems
x1 Close to equilibrium points “nonlinear system” ≈ “linear system”
Theorem: Assume
Im λ
Re λ ẋ = f (x) = Ax + g(x),
with limkxk→0 kg(x)k/kxk = 0. If ż = Az has a focus, node, or
saddle point, then ẋ = f (x) has the same type of equilibrium at the
origin.
Remark: If the linearized system has a center, then the nonlinear

system has either a center or a focus.
Lecture 3 1 Lecture 3 2
EL2620 2011 EL2620 2011
How to Draw Phase Portraits

By hand:
Example: Phase-Locked Loop
1. Find equilibria A PLL tracks phase θin (t) of a signal sin (t) = A sin[ωt + θin (t)].
2. Sketch local behavior around equilibria
sin “θout ”
3. Sketch (ẋ1 , ẋ2 ) for some other points. Notice that Phase Filter VCO
Detector
dx2 ẋ2
=
dx1 ẋ1
4. Try to find possible periodic orbits θin e K θ̇out 1 θout
sin(·)
5. Guess solutions − 1 + sT s
By computer:
1. Matlab: dee or pplane
EL2620 2011 EL2620 2011
Phase-Plane Analysis of PLL Classification of Equilibria
Let (x1 , x2 ) = (θout , θ̇out ), K, T > 0, and θin (t) ≡ θin . Linearization gives the following characteristic equations:
n even:
ẋ1 (t) = x2 (t)
λ2 + T −1 λ + KT −1 = 0
−1 −1
ẋ2 (t) = −T x2 (t) + KT sin(θin − x1 (t)) K > (4T )−1 gives stable focus
Equilibria are (θin + nπ, 0) since 0 < K < (4T )−1 gives stable node
ẋ1 = 0 ⇒ x2 = 0 n odd:
ẋ2 = 0 ⇒ sin(θin − x1 ) = 0 ⇒ x1 = θin + nπ λ2 + T −1 λ − KT −1 = 0
Saddle points for all K, T >0
EL2620 2011 EL2620 2011
Phase-Plane for PLL Periodic Solutions

Example of an asymptotically stable periodic solution:
(K, T ) = (1/2, 1): focuses 2kπ, 0 , saddle points (2k + 1)π, 0
ẋ1 = x1 − x2 − x1 (x21 + x22 )
(1)
ẋ2 = x1 + x2 − x2 (x21 + x22 )
EL2620 2011 EL2620 2011
Periodic solution: Polar coordinates.

x1 = r cos θ ⇒ ẋ1 = cos θṙ − r sin θθ̇
x2 = r sin θ ⇒ ẋ2 = sin θṙ + r cos θθ̇ A system has a periodic solution if for some T >0
implies
ṙ

1

r cos θ r sin θ

ẋ1
x(t + T ) = x(t), ∀t ≥ 0
=
θ̇ r − sin θ cos θ ẋ2 A periodic orbit is the image of x in the phase portrait.
Now, from (1)
ẋ1 = r(1 − r2 ) cos θ − r sin θ • When does there exist a periodic solution?
2
ẋ2 = r(1 − r ) sin θ + r cos θ • When is it stable?
gives
ṙ = r(1 − r2 )
Note that x(t) ≡ const is by convention not regarded periodic
θ̇ = 1
Only r = 1 is a stable equilibrium!
EL2620 2011 EL2620 2011
Poincaré Map
Assume φt (x0 ) is a periodic solution with period T .
Let Σ ⊂ Rn be an n − 1-dim hyperplane transverse to f at x0 .
Flow Definition: The Poincaré map P : Σ → Σ is
The solution of ẋ = f (x) is sometimes denoted P (x) = φτ (x) (x)

φt (x0 ) where τ (x) is the time of first return.
to emphasize the dependence on the initial point x0 ∈ Rn Σ

φt (·) is called the flow.
P (x)
x φt (x)
EL2620 2011 EL2620 2011
Existence of Periodic Orbits Stable Periodic Orbit

∗ ∗ ∗
A point x such that P (x ) = x corresponds to a periodic orbit. The linearization of P around x∗ gives a matrix W such that
P (x) ≈ W x
P (x∗ ) = x∗ if x is close to x∗ .
• λj (W ) = 1 for some j
• If |λi (W )| < 1 for all i 6= j , then the corresponding periodic
orbit is asymptotically stable
• If |λi (W )| > 1 for some i, then the periodic orbit is unstable.

x∗ is called a fixed point of P .
EL2620 2011 EL2620 2011
Example—Stable Unit Circle

Rewrite (1) in polar coordinates: The Poincaré map is
ṙ = r(1 − r2 )
[1 + (r0−2 − 1)e−2·2π ]−1/2

P (r0 , θ0 ) =
θ̇ = 1 θ0 + 2π
Choose Σ = {(r, θ) : r > 0, θ = 2πk}. (r0 , θ0 ) = (1, 2πk) is a fixed point.
The solution is The periodic solution that corresponds to (r(t), θ(t)) = (1, t) is
asymptotically stable because
φt (r0 , θ0 ) = [1 + (r0−2 − 1)e−2t ]−1/2 , t + θ0 −4π
dP e 0
W = (1, 2πk) =
First return time from any point (r0 , θ0 ) ∈ Σ is d(r0 , θ0 ) 0 1
τ (r0 , θ0 ) = 2π. ⇒ Stable periodic orbit (as we already knew for this example)
EL2620 2011
Next Lecture
• Lyapunov methods for stability analysis
Lecture 3 17

Lec03 11 PDF

Uploaded by

Copyright:

Available Formats

Lec03 11 PDF

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Lec03 11 PDF

Uploaded by

Copyright:

Available Formats

EL2620 2011 EL2620 2011

Equilibrium Points for Linear Systems

Remark: If the linearized system has a center, then the nonlinear

EL2620 2011 EL2620 2011

How to Draw Phase Portraits

Phase-Plane Analysis of PLL Classification of Equilibria

EL2620 2011 EL2620 2011

Phase-Plane for PLL Periodic Solutions

Periodic solution: Polar coordinates.

EL2620 2011 EL2620 2011

The solution of ẋ = f (x) is sometimes denoted P (x) = φτ (x) (x)

to emphasize the dependence on the initial point x0 ∈ Rn Σ

Existence of Periodic Orbits Stable Periodic Orbit

• If |λi (W )| > 1 for some i, then the periodic orbit is unstable.

EL2620 2011 EL2620 2011

Example—Stable Unit Circle

You might also like