Linear Control Systems Lecture # 23 Time-Varying Systems Solution of The State Equation
Linear Control Systems Lecture # 23 Time-Varying Systems Solution of The State Equation
Linear Control Systems Lecture # 23 Time-Varying Systems Solution of The State Equation
Lecture # 23
Time-Varying Systems
Solution of the State Equation
p. 1/2
x(t0 ) = x0
A( )x( ) d
t0
p. 2/2
Successive Approximation
x0 (t) = x0
x1 (t) = x0 +
x2 (t) = x0 +
= x0 +
A(1 )x0 d1
t0
Z t
t0
Z t
A(1 )x1 (1 ) d1
A(1 )x0 d1
t0
..
.
A(1 )
t0
xk (t) = x0 +
A(2 )x0 d2 d1
t0
t
t0
A(1 )xk1 (1 ) d1
p. 3/2
k1
X
j=0
Z t
A(
)x
d
kx1 (t) x0 (t)k =
1
0
1
t0
Z t
kA(1 )k kx0 k d1 T kx0 k
t0
p. 4/2
Z t
x
(
)]
d
1
1
1
0
1
1
t0
Z t
kA(1 )k kx1 (1 ) x0 (1 )k d1
t0
Z t
T kx0 k d1
t0
2
= T kx0 k(t t0 )
p. 5/2
Z t
x
(
)]
d
1
2
1
1
1
1
t0
Z t
kA(1 )k kx2 (1 ) x1 (1 )k d1
t0
Z t
2 T kx0 k(1 t0 ) d1
t0
= T kx0 k
2 (t t0 )2
2!
j (t t0 )j
j!
T kx0 k
j T j
j!
p. 6/2
xk (t) = x0 +
k1
X
j=0
k1
X
j=0
kx0 k + T kx0 k
k1
X
j=0
k1
X
j=0
j T j
j!
j T j
j!
eT as k
p. 7/2
A(1 )
t0
A(2 ) d2 d1 x0 +
t0
= 0 + A(t)x0 + A(t)
= A(t)[x0 +
= A(t)x(t)
A(2 ) d2 x0 +
t0
A(2 ) d2 x0 + ]
t0
p. 8/2
Z t
x(t) = I +
A(1 ) d1
t0
A(1 )
t0
A(2 ) d2 d1 + x0
t0
A(1 ) d1
t0
t
A(1 )
t0
A(2 ) d2 d1 +
t0
p. 9/2
p. 10/2
Properties of (t, ):
1.
(t, t) = I
2.
(t, ) = A(t)(t, )
3.
4.
t, ,
p. 11/2
5. (t, )
Proof of Property 3:
(t, ) = I +
+
A(1 ) d1
A(1 )
A(2 ) d2 d1 +
A(1 ) d1 = A( )
p. 12/2
A(1 )
= A( )
=
=
Z
Z
A(2 ) d2 d1
() d2 +
A(1 )
A(2 ) d2 d1
A(1 )A( ) d1
A(1 ) d1 A( )
(t, ) = I +
A(1 ) d1 + A( )
= (t, )A( )
p. 13/2
x( )
p. 14/2
Proof of Property 5:
(t, )(, t) = (t, t) = I
p. 15/2
Complete Solution
x(t0 ) = x0
Proof:
x(t0 ) = (t0 , t0 )x0 = x0
p. 16/2
+
(t, )B( )u( ) d
t0 t
= A(t)(t, t0 )x0 + B(t)u(t)
Z t
+
A(t)(t, )B( )u( ) d
x =
t0
p. 17/2
A d1 +
Z
Z
A d2 d1
A d3 d2 d1 +
A d1 = (t )A
A d2 d1 =
(t )2
2!
A d3 d2 d1 =
(t, ) =
X
(t )k
k=0
k!
A2
(t )3
3!
A3
Ak = e(t )A
p. 18/2
Special Case: n = 1
x = a(t)x
dx
= a(t) dt
x
Z
ln
x(t)
x(t0 )
dz
z
x(t)
x(t0 )
x(t) = e
Rt
t0
a( ) d
t0
a( ) d
t0
a( ) d
(t, t0 ) = e
Rt
t0
x(t0 )
a( ) d
p. 19/2
a1 (t)
A(t) =
(t, t0 ) =
a2 (t)
...
an (t)
1 (t, t0 )
2 (t, t0 )
...
n (t, t0 )
i (t, t0 ) = e
Rt
t0
ai ( ) d
p. 20/2
A(t) =
a22 (t) a23 (t)
a33 (t)
(t, t0 ) =
22 (t, t0 ) 23 (t, t0 )
33 (t, t0 )
Rt
t0
aii ( ) d
p. 21/2
Example
A=
"
2
t
0 1
11 (t, t0 ) = e2(tt 0 )
p. 22/2
e2(t ) e( t 0 ) x20 d
t0
e d = et (t 1) et0 (t0 1)
t0
x1 (t)
x2 (t)
"
e2(tt 0 ) 12 (t, t0 )
0
e(tt 0 )
#"
x10
x20
p. 23/2
Example:
x 1 = tx1 + tx2 ,
x2 (t) = e
x1 (t) = e
= e
2
2
1
(t
t
0)
2
2
2
1
(t
t
0)
2
Rt
t0
x1 (t0 ) +
x2 (t0 ) = e
Z
x1 (t0 ) + e
(t, t0 ) =
"
2
2
1
(t
t
0)
2
2
2
1
(t
)
2
x2 (t0 )
( 2 t20 )/2
x2 (t0 ) d
t0
2
2
1
(t
t
0)
2
1
(t2 t20 )
2
x 2 = tx2
1 2
2 (t
d x2 (t0 )
t0
1
(t2 t20 )
2
t20 )e
2
2
1
(t
t
0)
e2
p. 24/2