EE C222/ME C237 - Spring'18 - Lecture 3 Notes: Murat Arcak January 24 2018
EE C222/ME C237 - Spring'18 - Lecture 3 Notes: Murat Arcak January 24 2018
EE C222/ME C237 - Spring'18 - Lecture 3 Notes: Murat Arcak January 24 2018
Invariant Sets
n( x )
M f (x)
( dc , ba )
saddle
%
x
Example 2:
ẋ1 = x1 + x2 − x1 ( x12 + x22 )
ẋ2 = −2x1 + x2 − x2 ( x12 + x22 )
x1
n( x ) =
x2 x2
f ( x ) · n( x ) = x12 + x1 x2 − x12 ( x12 + x22 ) − 2x1 x2 + x22 − x22 ( x12 + x22 )
= − x1 x2 + ( x12 + x22 ) − ( x12 + x22 )2
f (x)
1 1
− x1 x2 ≤ x12 + x22 (completion of squares) x1
2 2
Therefore, f ( x ) · n( x ) ≤ 23 r2 − r4 ≤ 0 if r2 ≥ 32 .
Two criteria:
x2
x1
q
3
Example 2 above: Br is positively invariant for r ≥ 2 but contains
the equilibrium x = 0.
" #
√
∂ f 1 1
= λ1,2 = 1 ∓ j 2 unstable focus.
∂x x=0
−2 1
A more general form of the PB Theorem states that, for time invari-
ant, planar systems, bounded trajectories converge to equilibria,
periodic orbits, or unions of equilibria connected by trajectories.
Corollary: No chaos for time invariant planar systems.
Index Theory
saddle -1
c0
The following corollary is useful for ruling out periodic orbits (like
Bendixson’s Theorem studied in the previous lecture):
Corollary: Inside any periodic orbit there must be at least one equi-
librium and the indices of the equilibria enclosed must add up to
+1.
Example (from last lecture):
ẋ1 = x2
ẋ2 = −δx2 + x1 − x13 + x12 x2 δ>0
The eigenvalues are real and have opposite signs, therefore (0, 0) is a
saddle: index = −1.
" #
∂f 0 1
∂x x =(∓1,0) = λ2 + (δ − 1)λ +2 = 0.
−2 1 − δ |{z}
>0
The eigenvalues are either real with the same sign (node) or complex
conjugates (focus or center), therefore (∓1, 0) each has index= +1.
Thus, the corollary above rules out the periodic orbit in the middle
plot below. It does not rule out the others, but does not prove their
existence either. Bendixson’s Criterion rules out neither of the three.
ee c222/me c237 - spring’18 - lecture 3 notes 5
x2
x1
√ √
x1 = − δ x1 = δ
x2
not possible
x1
√ √
x1 = − δ x1 = δ
x2
x1
√ √
x1 = − δ x1 = δ