Introduction To Nonlinear Control Lecture # 1 Preliminaries
Introduction To Nonlinear Control Lecture # 1 Preliminaries
Introduction To Nonlinear Control Lecture # 1 Preliminaries
Lecture # 1
Preliminaries
p. 1/6
p. 2/6
x=
x1
x2
..
.
..
.
xn
f1 (t, x, u)
u1
f2 (t, x, u)
..
, u =
.
, f (t, x, u) =
..
..
u
p
fn (t, x, u)
x = f (t, x, u)
p. 3/6
x = f (t, x, u)
y = h(t, x, u)
x is the state, u is the input
y is the output (q -dimensional vector)
Special Cases:
Linear systems: x = A(t)x + B(t)u
y = C(t)x + D(t)u
p. 4/6
Time-Invariant System:
x = f (x, u)
y = h(x, u)
p. 5/6
Linearization
A common engineering practice in analyzing a nonlinear
system is to linearize it about some nominal operating point
and analyze the resulting linear model
What are the limitations of linearization?
Since linearization is an approximation in the
neighborhood of an operating point, it can only predict
the local behavior of the nonlinear system in the
vicinity of that point. It cannot predict the nonlocal or
global behavior
There are essentially nonlinear phenomena that can
take place only in the presence of nonlinearity
p. 6/6
Nonlinear Phenomena
Finite escape time
Multiple isolated equilibrium points
Limit cycles
Subharmonic, harmonic, or almost-periodic oscillations
Chaos
Multiple modes of behavior
p. 7/6
p. 8/6
|f (y) f (x)|
|y x|
(Finite Slope)
p. 9/6
p. 10/6
1
(t 1)
x(t) as t 1
p. 11/6
p. 12/6
with kA(t)k L
p. 13/6
Equilibrium Points
A point x is an equilibrium point of x = f (t, x) if
x(t0 ) = x x(t) x , t t0
p. 14/6
p. 15/6
Pendulum Equation
ml = mg sin kl
x1 = ,
x2 =
x 1 = x2 ,
mg
x 2 =
g
l
sin x1
k
m
x2
Equilibrium Points:
0 = x2 ,
0=
g
l
sin x1
k
m
x2
(n, 0) for n = 0, 1, 2, . . .
p. 16/6
Second-Order Systems
x 1 = f1 (x1 , x2 ) = f1 (x)
x 2 = f2 (x1 , x2 ) = f2 (x)
f2 (x)
f1 (x)
p. 17/6
f (x)
*
q
x = (1, 1)
x1
p. 18/6
6
4
x2
2
0
2
4
6
x 1 = x2 ,
0
x1
x 2 = 10 sin x1
p. 19/6
x(0) = x0
x(0) = x0
p. 20/6
2 /1
cz1
,
z 2 = 2 z2
z2 (t) = z2 (0)e2 t
c = z2 (0)/(z1 (0))2 /1
p. 21/6
2 < 1 < 0
z2 =
/
cz1 2 1
with
2
1
> 1 and
dz2
dz1
=c
2
1
[(2 /1 )1]
z1
z1
Stable Node
p. 22/6
x2
v2
(a)
Stable Node
2 < 1 < 0
x2
v2
v1
v1
x1
x1
(b)
Unstable Node
2 > 1 > 0
p. 23/6
2 < 0 < 1
e1 t , while e2 t 0 as t
z2 =
2 /1
cz1
,
2 /1 < 0
x2
z2
v2
v1
x1
z1
(a)
(b)
Saddle
p. 24/6
= tan
z 2 = z1 + z2
1 z2
z1
p. 25/6
z2
(a)
z2
z1
=0
Unstable Focus
x2
(b)
x1
(c)
z1
>0
Stable Focus
z1
<0
(a)
(b)
x2
Center
(c)
x2
x1
p. 26/6
Effect of Perturbations
A A + A
(A arbitrarily small)
p. 27/6
p. 28/6
p. 29/6
x 2 = f2 (x1 , x2 )
f1 (x1 , x2 )
=
,
x1
x=p
f2 (x1 , x2 )
=
,
x1
x=p
a12
a22
f1 (x1 , x2 )
=
x2
x=p
f2 (x1 , x2 )
=
x2
x=p
p. 30/6
f1 (p1 , p2 ) = f2 (p1 , p2 ) = 0
y1 = x1 p1
y2 = x2 p2
A=
a11 a12
a21 a22
y Ay
f
f1
1
x1 x2
f2
f2
x1
x2
x=p
f
=
x x=p
p. 31/6
p. 32/6
Multiple Equilibria
iL
X
X
R
E
P
P
P
P
P
P
P
P
P
P
L
+ vL
s
CC
i,mA
iC
CC
iR
i=h(v)
0.5
vC C
vR
0.5
(a)
0.5
1 v,V
(b)
iC = C
dvC
dt
vL = L
diL
dt
x1 = vC , x2 = iL
p. 33/6
x 1 =
1
C
[h(x1 ) + x2 ] ,
x 2 =
1
L
[x1 Rx2 + E]
Equilibrium Points:
0 = h(x1 ) + x2 ,
0 = x1 Rx2 + E
i
R
1.2
1
h(x1 ) =
E
R
1
R
0.8
x1
0.6
Q
1
Q
2
0.4
Q
3
0.2
0
0
0.5
p. 34/6
x 1 = 0.5[h(x1 ) + x2 ]
x 2 = 0.2(x1 1.5x2 + 1.2)
h(x1 ) = 17.76x1 103.79x21 +229.62x31 226.31x41 +83.72x51
i
R
1
Q1 = (0.063, 0.758)
Q2 = (0.285, 0.61)
Q3 = (0.884, 0.21)
0.8
0.6
Q
1
Q
2
0.4
Q
3
0.2
0
0
0.5
1 v
p. 35/6
f
x
"
"
3.598 0.5
A1 =
,
0.2 0.3
"
#
1.82 0.5
A2 =
,
0.2 0.3
"
#
1.427 0.5
A3 =
,
0.2 0.3
p. 36/6
1.5
0.5
0.5
0.5
0.5
1.5
p. 37/6
z2 (t) = r0 sin(t + 0 )
1 z2 (0)
0 = tan
z1 (0)
p. 38/6
p. 39/6
Limit Cycles
C
CC
iC
CC
X
X
iL
i = h(v)
Resistive
Element
(a)
(b)
h(0) = 0,
h (0) < 0
p. 40/6
iC + iL + i = 0
C
dv
dt
CL
1
L
d2 v
dt2
v(s) ds + h(v) = 0
+ v + Lh (v)
= t/ CL,
v =
dv
d
dv
dt
=0
v
+ h (v)v + v = 0
L/C
p. 41/6
State model: x1 = v,
x 1 = x2 ,
x2 = v
x 2 = x1 h (x1 )x2
0
1
f
A=
=
x x=0
1 h (0)
2 + h (0) + 1 = 0
p. 42/6
Energy Analysis:
2
1
1
1
1
2
2
2
E = 2 CvC + 2 LiL = Cx1 + 2 L h(x1 ) x2
E = Cx1 h(x1 )
p. 43/6
h(v) = v + 31 v 3
x 1 = x2
x 2 = x1 + (1 x21 )x2
4
x2
0
x
1
2
2
(a)
(b)
= 0.2
=1
p. 44/6
(a)
x
1
(b)
p. 45/6
x = f (x)
f is locally Lipschitz over a domain D Rn and x
D is
x) = 0
an equilibrium point; that is, f (
def
y = x = f (x) = f (y + x
) = g(y),
where g(0) = 0
p. 46/6
t0
p. 47/6
First-Order Systems (n = 1)
The behavior of x(t) in the neighborhood of the origin can
be determined by examining the sign of f (x)
The requirement for stability is violated if xf (x) > 0 on
either side of the origin
f(x)
f(x)
Unstable
f(x)
Unstable
Unstable
p. 48/6
Stable
f(x)
f(x)
Stable
Stable
p. 49/6
f(x)
(a)
Asymptotically Stable
(b)
p. 50/6
x(0) = x0
p. 51/6
Second-Order Systems (n = 2)
Type of equilibrium point
Center
Stable Node
Stable Focus
Unstable Node
Unstable Focus
Saddle
Stability Property
p. 52/6
1.5
0.5
0.5
0.5
0.5
1.5
p. 53/6
8
4
0
x
p. 54/6
8
4
0
x
p. 55/6
i 1
0
i
.
..
Ji = .
..
..
.
0 ...
block diag[J1 , J2 , . . . , Jr ]
0 ... ... 0
1
0 ... 0
..
...
.
...
0
...
1
. . . . . . 0 i m m
i
p. 56/6
exp(At) = P exp(Jt)P 1 =
r X
mi
X
i=1 k=1
p. 57/6
p. 58/6
Linearization
x = f (x),
f (0) = 0
f
x
(x)
h () = J(x)x
Z 1
h(1) h(0) =
h () d, h(0) = f (0) = 0
0
f (x) =
J(x) d x
0
p. 59/6
f (x) =
J(x) d x
0
[J(x) J(0)] d
G(x) 0 as x 0
p. 60/6
Theorem:
The origin is exponentially stable if and only if
Re[i ] < 0 for all eigenvalues of A
The origin is unstable if Re[i ] > 0 for some i
Linearization fails when Re[i ] 0 for all i, with
Re[i ] = 0 for some i
Example
x = ax3
f
2
A=
= 3ax x=0 = 0
x x=0
p. 61/6
Example
x = x3
x(t) = p
1 + 2tx2 (0)
e2t
1 + 2tx2 (0)
e2t
1+
2tx2 (0)
k2
=
p. 62/6