Introduction To Nonlinear Control Lecture # 1 Preliminaries

Download as pdf or txt
Download as pdf or txt
You are on page 1of 62

Introduction to Nonlinear Control

Lecture # 1
Preliminaries

p. 1/6

Nonlinear State Model


x 1 = f1 (t, x1 , . . . , xn , u1 , . . . , up )
x 2 = f2 (t, x1 , . . . , xn , u1 , . . . , up )
..
..
.
.
x n = fn (t, x1 , . . . , xn , u1 , . . . , up )
x i denotes the derivative of xi with respect to the time
variable t
u1 , u2 , . . ., up are input variables
x1 , x2 , . . ., xn the state variables

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

Unforced state equation:


x = f (t, x)

It results from x = f (t, x, u) with u = (t, x)

p. 4/6

Time-Invariant System:
x = f (x, u)
y = h(x, u)

A time-invariant state model has a time-invariance property


with respect to shifting the initial time from t0 to t0 + a,
provided the input waveform is applied from t0 + a rather
than t0
Autonomous System:
x = f (x)

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

Existence and Uniqueness of Solutions


x = f (t, x)
f (t, x) is piecewise continuous in t and locally Lipschitz in
x over the domain of interest
f (t, x) is piecewise continuous in t on an interval J R if
for every subinterval [a, b] J , f is continuous in t for all
t [a, b], except, possibly, at a finite number of points
where f may have finite-jump discontinuities
f (t, x) is locally Lipschitz in x at a point x0 if there is a
neighborhood N (x0 , r) = {x Rn | kx x0 k < r}
where f (t, x) satisfies the Lipschitz condition
kf (t, x) f (t, y)k Lkx yk, L > 0

p. 8/6

A function f (t, x) is locally Lipschitz in x on a domain


(open and connected set) D Rn if it is locally Lipschitz at
every point x0 D
n=1:

|f (y) f (x)|
|y x|

(Finite Slope)

If f (x) is discontinuous or has infinite slope at x0 , it is not


locally Lipschitz at x0 . For example, f (x) = x1/3 is not
locally Lipschitz at x = 0. On the other hand, if f (x) is
continuous at x0 then |f (x)| k in a neighborhood of x0 .
Hence, f (x) satisfies the Lipschitz condition with L = k

p. 9/6

More generally, if f (t, x) and its partial derivatives fi /xj


are continuous for t J and x D , then f (t, x) is locally
Lipschitz in x on J D
Lemma: Let f (t, x) be piecewise continuous in t and
locally Lipschitz in x at x0 , for all t [t0 , t1 ]. Then, there is
> 0 such that the state equation x = f (t, x), with
x(t0 ) = x0 , has a unique solution over [t0 , t0 + ]
Without the local Lipschitz condition, we cannot ensure
uniqueness of the solution. For example, x = x1/3 with
x(0) = 0 has x(t) = (2t/3)3/2 and x(t) 0 as two
different solutions

p. 10/6

The lemma is a local result because it guarantees existence


and uniqueness of the solution over an interval [t0 , t0 + ],
but this interval might not include a given interval [t0 , t1 ].
Indeed the solution may cease to exist after some time
Example:
x = x2

f (x) = x2 is locally Lipschitz for all x


x(0) = 1 x(t) =

1
(t 1)

x(t) as t 1

the solution has a finite escape time at t = 1

p. 11/6

Global Existence and Uniqueness


A function f (t, x) is globally Lipschitz in x if
kf (t, x) f (t, y)k Lkx yk

for all x, y Rn with the same Lipschitz constant L


If f (t, x) and its partial derivatives fi /xj are continuous
for all x Rn , then f (t, x) is globally Lipschitz in x if and
only if the partial derivatives fi /xj are globally bounded,
uniformly in t
f (x) = x2 is locally Lipschitz for all x but not globally
Lipschitz because f (x) = 2x is not globally bounded

p. 12/6

Lemma: Let f (t, x) be piecewise continuous in t and


globally Lipschitz in x for all t [t0 , t1 ]. Then, x = f (t, x),
with x(t0 ) = x0 , has a unique solution over [t0 , t1 ]
The global Lipschitz condition holds for the linear system
x = A(t)x + g(t)

with kA(t)k L

but it is a restrictive condition for general nonlinear systems


Lemma: Let f (t, x) be piecewise continuous in t and
locally Lipschitz in x for all t t0 and all x in a domain D .
Let W be a compact subset of D , and suppose that every
solution of x = f (t, x), with x(t0 ) W lies entirely in W .
Then, there is a unique solution that is defined for all t t0

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

For the autonomous system x = f (x), the equilibrium


points are the real solutions of the equation
f (x) = 0

An equilibrium point could be isolated; that is, there are no


other equilibrium points in its vicinity, or there could be a
continuum of equilibrium points

p. 14/6

A linear system x = Ax can have an isolated equilibrium


point at x = 0 (if A is nonsingular) or a continuum of
equilibrium points in the null space of A (if A is singular)
It cannot have multiple isolated equilibrium points , for if xa
and xb are two equilibrium points, then by linearity any point
on the line xa + (1 )xb connecting xa and xb will be
an equilibrium point
A nonlinear state equation can have multiple isolated
equilibrium points

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)

Let x(t) = (x1 (t), x2 (t)) be a solution that starts at initial


state x0 = (x10 , x20 ). The locus in the x1 x2 plane of the
solution x(t) for all t 0 is a curve that passes through the
point x0 . This curve is called a trajectory or orbit
The x1 x2 plane is called the state plane or phase plane
The family of all trajectories is called the phase portrait
The vector field f (x) = (f1 (x), f2 (x)) is tangent to the
trajectory at point x because
dx2
dx1

f2 (x)
f1 (x)

p. 17/6

Vector Field diagram


Represent f (x) as a vector based at x; that is, assign to x
the directed line segment from x to x + f (x)
x2
x + fq (x) = (3, 2)



f (x)


*
q 


x = (1, 1)
x1

Repeat at every point in a grid covering the plane

p. 18/6

6
4

x2

2
0

2
4
6

x 1 = x2 ,

0
x1

x 2 = 10 sin x1

p. 19/6

Numerical Construction of the Phase Portrait:


Select a bounding box in the state plane
Select an initial point x0 and calculate the trajectory
through it by solving
x = f (x),

x(0) = x0

in forward time (with positive t) and in reverse time (with


negative t)
x = f (x),

x(0) = x0

Repeat the process interactively


Use Simulink or pplane

p. 20/6

Linear Second-Order Systems: x = Ax


Case 1: Real distinct eigenvalues: 1 6= 2 6= 0
#
"
1 0
z
x(t) = M z(t) = [v1 , v2 ]z(t) z =
0 2
v1 & v2 are the eigenvectors associated with 1 & 2
z 1 = 1 z1 ,
z1 (t) = z1 (0)e1 t ,
z2 =

2 /1
cz1
,

z 2 = 2 z2
z2 (t) = z2 (0)e2 t

c = z2 (0)/(z1 (0))2 /1

The shape of the trajectory depends on the sign of 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

e1 t and e2 t tend to zero as t


z2

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

Case 2: Complex eigenvalues: 1,2 = j


x(t) = M z(t) z 1 = z1 z2 ,
q
r = z12 + z22 ,
r = r,

= tan

z 2 = z1 + z2
 
1 z2
z1

r(t) = r0 et and (t) = 0 + t


< 0 r(t) 0 as t
> 0 r(t) as t
= 0 r(t) r0 t

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)

The eigenvalues of a matrix depend continuously on its


parameters
A node (with distinct eigenvalues), a saddle or a focus is
structurally stable because the qualitative behavior remains
the same under arbitrarily small perturbations in A
A stable node with multiple eigenvalues could become a
stable node or a stable focus under arbitrarily small
perturbations in A

p. 27/6

A center is not structurally stable


"
#
1
1
Eigenvalues = j
< 0 Stable Focus
> 0 Unstable Focus

p. 28/6

The qualitative behavior of a nonlinear system near an


equilibrium point can take one of the patterns we have seen
with linear systems. Correspondingly the equilibrium points
are classified as stable node, unstable node, saddle, stable
focus, unstable focus, or center
Can we determine the type of the equilibrium point of a
nonlinear system by linearization?

p. 29/6

Let p = (p1 , p2 ) be an equilibrium point of the system


x 1 = f1 (x1 , x2 ),

x 2 = f2 (x1 , x2 )

where f1 and f2 are continuously differentiable


Expand f1 and f2 in Taylor series about (p1 , p2 )
x 1 = f1 (p1 , p2 ) + a11 (x1 p1 ) + a12 (x2 p2 ) + H.O.T.
x 2 = f2 (p1 , p2 ) + a21 (x1 p1 ) + a22 (x2 p2 ) + H.O.T.
a11
a21


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

y 1 = x 1 = a11 y1 + a12 y2 + H.O.T.


y 2 = x 2 = a21 y1 + a22 y2 + H.O.T.

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

Eigenvalues of A Type of equilibrium point


of the nonlinear system
2 < 1 < 0
Stable Node
2 > 1 > 0
Unstable Node
2 < 0 < 1
Saddle
j, < 0
Stable Focus
j, > 0
Unstable Focus
j
Linearization Fails

p. 32/6

Multiple Equilibria

Example: Tunnel-Diode Circuit

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
"

"

0.5h (x1 ) 0.5


0.2
0.3
#

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

Eigenvalues : 3.57, 0.33


Eigenvalues : 1.77, 0.25
Eigenvalues : 1.33, 0.4

Q1 is a stable node; Q2 is a saddle; Q3 is a stable node

p. 36/6

x = 0.5 ( 17.76 x + 103.79 x2 229.62 x3 + 226.31 x4 83.72 x5 + y)


y = 0.2 ( x 1.5 y + 1.2)

1.5

0.5

0.5
0.5

0.5

1.5

p. 37/6

Oscillation: A system oscillates when it has a nontrivial


periodic solution
x(t + T ) = x(t), t 0

Linear (Harmonic) Oscillator:


"
#
0
z =
z
0
z1 (t) = r0 cos(t + 0 ),
q
r0 = z12 (0) + z22 (0),

z2 (t) = r0 sin(t + 0 )


1 z2 (0)
0 = tan
z1 (0)

p. 38/6

The linear oscillation is not practical because


It is not structurally stable. Infinitesimally small
perturbations may change the type of the equilibrium
point to a stable focus (decaying oscillation) or unstable
focus (growing oscillation)
The amplitude of oscillation depends on the initial
conditions
The same problems exist with oscillation of nonlinear
systems due to a center equilibrium point (e.g., pendulum
without friction)

p. 39/6

Limit Cycles

Example: Negative-Resistance Oscillator





L



C
CC

iC

CC

X
X



iL

i = h(v)
Resistive
Element

(a)

(b)

h(0) = 0,

h (0) < 0

h(v) as v , and h(v) as v

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

There is a unique equilibrium point at the origin


0
1
f

A=
=


x x=0
1 h (0)
2 + h (0) + 1 = 0

h (0) < 0 Unstable Focus or Unstable Node

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

Example: Van der Pol Oscillator:

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)

Stable Limit Cycle

x
1

(b)

Unstable Limit Cycle

p. 45/6

Stability of Equilibrium Points

x = f (x)
f is locally Lipschitz over a domain D Rn and x
D is
x) = 0
an equilibrium point; that is, f (

Characterize and study the stability of x

Without loss of generality, take x


= 0. If x
6= 0, let
y =xx

def

y = x = f (x) = f (y + x
) = g(y),

where g(0) = 0

p. 46/6

Definition: The equilibrium point x = 0 of x = f (x) is


stable if for each > 0 there is > 0 such that
kx(0)k < kx(t)k < , t 0

unstable if it is not stable


asymptotically stable if it is stable and for some c > 0
kx(0)k < c lim x(t) = 0
t

exponentially stable if for some c > 0, > 0, k 1


kx(0)k < c kx(t)k kkx(0)ket ,

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

The origin is stable if and only if xf (x) 0 in some


neighborhood of the origin
f(x)

Stable

f(x)

f(x)

Stable

Stable

p. 49/6

The origin is asymptotically stable if and only if xf (x) < 0


in some neighborhood of the origin
f(x)

f(x)

(a)

Asymptotically Stable

(b)

Globally Asymptotically Stable

p. 50/6

Definition: Let the origin be an asymptotically stable


equilibrium point of the system x = f (x), where f is a
locally Lipschitz function defined over a domain D Rn
( 0 D)

The region of attraction (also called region of


asymptotic stability, domain of attraction, or basin) is the
set of all points x0 in D such that the solution of
x = f (x),

x(0) = x0

is defined for all t 0 and converges to the origin as t


tends to infinity
The origin is said to be globally asymptotically stable if
the region of attraction is the whole space Rn

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

Example: Tunnel Diode Circuit


x = 0.5 ( 17.76 x + 103.79 x2 229.62 x3 + 226.31 x4 83.72 x5 + y)
y = 0.2 ( x 1.5 y + 1.2)

1.5

0.5

0.5
0.5

0.5

1.5

p. 53/6

Example: Pendulum Without Friction


x=y
y = 10 sin(x)

8
4

0
x

p. 54/6

Example: Pendulum With Friction


x=y
y = 10 sin(x) y

8
4

0
x

p. 55/6

Linear Time-Invariant Systems


x = Ax
x(t) = exp(At)x(0)
P 1 AP = J =

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

tk1 exp(i t)Rik

i=1 k=1

mi is the order of the Jordan block Ji


Re[i ] < 0 i Exponentially Stable
Re[i ] > 0 for some i Unstable
Re[i ] 0 i & mi > 1 for Re[i ] = 0 Unstable
Re[i ] 0 i & mi = 1 for Re[i ] = 0 Stable

If an n n matrix A has a repeated eigenvalue i of


algebraic multiplicity qi , then the Jordan blocks of i have
order one if and only if rank(A i I) = n qi

p. 57/6

Theorem: The equilibrium point x = 0 of x = Ax is stable if


and only if all eigenvalues of A satisfy Re[i ] 0 and for
every eigenvalue with Re[i ] = 0 and algebraic multiplicity
qi 2, rank(A i I) = n qi , where n is the dimension
of x. The equilibrium point x = 0 is globally exponentially
stable if and only if all eigenvalues of A satisfy Re[i ] < 0
When all eigenvalues of A satisfy Re[i ] < 0, A is called a
Hurwitz matrix

p. 58/6

Linearization
x = f (x),

f (0) = 0

f is continuously differentiable over D = {kxk < r}


J(x) =

f
x

(x)

h() = f (x) for 0 1

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

Set A = J(0) and add and subtract Ax


f (x) = [A + G(x)]x, where G(x) =

[J(x) J(0)] d

G(x) 0 as x 0

This suggests that in a small neighborhood of the origin we


can approximate the nonlinear system x = f (x) by its
linearization about the origin x = Ax

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

Stable if a = 0; Asymp stable if a < 0; Unstable if a > 0


When a < 0, the origin is not exponentially stable

p. 61/6

Example
x = x3

The origin is asymptotically stable


x(0)

x(t) = p
1 + 2tx2 (0)

x(t) does not satisfy |x(t)| ket |x(0)| because


x(t)| ket |x(0)|

Impossible because lim

e2t
1 + 2tx2 (0)
e2t
1+

2tx2 (0)

k2
=

p. 62/6

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy