Lecture-9 Modeling in Time Domain
Lecture-9 Modeling in Time Domain
Lecture-9 Modeling in Time Domain
Note: I do not claim any originality in these lectures. The contents of this presentation are
mostly taken from the book of Ogata, Norman S Nise and various other internet sources. 1
Introduction
• Most classical control techniques were developed for linear constant
coefficient systems with one input and one output(perhaps a few
inputs and outputs).
• When nonlinearities and time variations are present, the very basis
for these classical techniques is removed.
dx
x2 dt
Velocity State (t=t1)
State (t=t1)
State
Vector
x1 Position
x
8
Example-1
• An electrical network is given in following figure, find a
state-space representation if the output is the current
through the resistor.
𝑑𝑣𝑐 (𝑡)
𝐶 = 𝑖𝑐
𝑑𝑡
𝑑𝑖𝐿 (𝑡)
𝐿 = 𝑣𝐿
𝑑𝑡
v c
• State Vector
iL
Example-1
• Step 3: Apply network theory, such as Kirchhoff’s voltage
and current laws, to obtain iC and vL in terms of the state
variables, vC and iL.
iL iR iC
𝑑𝑣𝑐 (𝑡)
iC iR iL 𝐶
𝑑𝑡
= 𝑖𝑐
𝑑𝑖𝐿 (𝑡)
dvC vC 𝐿 = 𝑣𝐿
C iL (1) 𝑑𝑡
dt R
Example-1
• Step 3: Apply network theory, such as Kirchhoff’s voltage
and current laws, to obtain iC and vL in terms of the state
variables, vC and iL.
diL
v(t ) L vR 𝑑𝑣𝑐 (𝑡)
dt 𝐶 = 𝑖𝑐
𝑑𝑡
𝑑𝑖𝐿 (𝑡)
diL 𝐿 = 𝑣𝐿
L vC v(t ) (2) 𝑑𝑡
dt
Example-1
dvC vC diL
C iL (1) L vC v(t ) (2)
dt R dt
Write equation (1) & (2) in standard form.
dvC 1 1
vC i L
dt RC C
State Equations
diL 1 1
vC v(t )
dt L L
Example-1
dvC 1 1 diL 1 1
vC i L vC v(t )
dt RC C dt L L
1 1
d vc RC v 0
C c 1 v(t )
1
dt iL 0 iL L
L
1 1
vc RC v 0
C c 1 v(t )
1
iL 0 iL L
L
Example-1
Step-4: The output is current through the resistor therefore, the
output equation is
1
iR vC
R
1 v c
iR 0
R iL
Example-1
1 1
vc RC v 0
C c 1 v(t )
1
iL 0 iL L
L Where,
x(t) --------------- State Vector
x(t ) Ax(t ) Bu(t ) A (nxn) ---------------- System Matrix
B (nxp) ----------------- Input Matrix
u(t) --------------- Input Vector
1 v c
iR 0
R iL Where,
y(t) -------------- Output Vector
C (qxn) ---------------- Output Matrix
y(t ) Cx(t ) Du (t ) D ----------------- Feed forward Matrix
Example-2
• Consider RLC Circuit shown below. Represent the system in Sate
Space and find (if L=1H, R=3Ω and C=0.5 F):
– State Vector
iL
– System Matrix
+ +
– Input Matrix & Input Vector Vc Vo
– Output Matrix & Output Vector - -
dvc diL
C u ( t ) iL L Ri L vc Vo Ri L
dt dt
• Choosing vc and iL as state variables
1 1 diL 1 R
dvc
iL u(t ) vc iL
dt C C dt L L
Example-2
dvc 1 1 diL 1 R
iL u(t ) vc iL
dt C C dt L L
1
vc 0 v 1
i 1
C
c
C u(t ) State Equation
L R iL
0
L L
Vo Ri L
vc
Vo 0 R Output Equation
iL
Example-3
• State Space Representation of Series RLC circuit
𝑑𝑞(𝑡) 𝑑 2 𝑞(𝑡) 1
𝑣 𝑡 = 𝑅+𝐿 2
+ 𝑞(𝑡)
𝑑𝑡 𝑑𝑡 𝐶
Example-3
𝑑𝑞(𝑡) 𝑑2 𝑞(𝑡) 1
𝑣 𝑡 = 𝑅+𝐿 2
+ 𝑞(𝑡)
𝑑𝑡 𝑑𝑡 𝐶
• We can convert above equation into two simultaneous, first-
order differential equations in terms of i(t) and q(t).
𝑑𝑞(𝑡)
= 𝑖(𝑡)
𝑑𝑡
• The second equation can be formed as
𝑑𝑖(𝑡) 𝑅 1 1
=− 𝑖 𝑡 − 𝑞 𝑡 + 𝑣 𝑡
𝑑𝑡 𝐿 𝐿𝐶 𝐿
Example-3
𝑑𝑞(𝑡) 𝑑𝑖(𝑡) 𝑅 1 1
= 𝑖(𝑡) =− 𝑖 𝑡 − 𝑞 𝑡 + 𝑣 𝑡
𝑑𝑡 𝑑𝑡 𝐿 𝐿𝐶 𝐿
q 0 1 q 0
i 1 R 1 v(t )
i
LC L L
• If capacitor Voltage 𝑣𝑐 𝑡 is the output then
1 1
𝑣𝑐 𝑡 = න 𝑖(𝑡) 𝑑𝑡 = 𝑞(𝑡)
𝐶 𝐶
1 q
vc 0
C i
Example-4: Home Work
• Find the state-space representation of the electrical
network shown in following Figure. The output is vo(t).
Example-5
• Represent the electrical network shown in following
figure in state space, where iR(t) is the output.
Example-5
Example-5
𝑑𝑣1 1 𝑑𝑖3 1
= 𝑖2 = 𝑣𝐿
𝑑𝑡 3 𝑑𝑡 2
𝑣𝐿 = 𝑣1 − 3(𝑖3 + 4𝑣1 )
𝑣𝐿 = −11𝑣1 − 3𝑖3
𝑑𝑣1 1 1 1
= − 𝑣1 − 𝑖3 + 𝑣𝑖
𝑑𝑡 9 3 9
Example-5
𝑑𝑣1 1 1 1 𝑑𝑖3 11 3
= − 𝑣1 − 𝑖3 + 𝑣𝑖 = − 𝑣1 − 𝑖3
𝑑𝑡 9 3 9 𝑑𝑡 2 2
• Writing state equations in vector matrix form
1 1
𝑣1ሶ − − 𝑣 1
= 9 3 1
+ 9 𝑣𝑖
𝑖3ሶ 11 3 𝑖3
− − 0
2 2
𝑖𝑅 = 4𝑣1 +𝑖3
𝑣1
𝑖𝑅 = 4 1 𝑖
3
Example-6
• Consider the following system
K
x(t)
M
f(t)
B
Differential equation of the system is:
d 2 x(t ) dx(t )
M 2
B Kx(t ) f (t )
dt dt
30
Example-6
• As we know
dx
v d 2x dv
dt dt 2 dt
• Choosing x and v as state variables
d 2 x(t ) dx(t )
dx M B Kx(t ) f (t )
v dt 2 dt
dt
dv B K 1
v x f (t )
dt M M M
x 0 1 x 0
K B 1 f (t )
v
v M M M
Example-6
x 0 1 x 0
K B 1 f (t )
v
v M M M
• If velocity v is the out of the system then output equation is given as
x
y (t ) 0 1
v
Example-7
• Find the state equations of following mechanical translational
system.
d 2 x1 dx1 d 2 x2
M1 2
D Kx1 Kx2 0 f (t ) M 2 2
Kx2 Kx1
dt dt dt
Example-7
• Now dx1 d 2 x1 dv1
v1
dt dt 2 dt
dx2 2
v2 d x2 dv2
2
dt dt dt
• Choosing x1, v1, x2, v2 as state variables
dx1
v1
dt
dv1
M1 Dv1 Kx1 Kx2 0
dt
dx2
v2
dt
dv2
f (t ) M 2 Kx2 Kx1
dt
Example-7
• In Standard form
dx1
v1
dt
dv1 D K K
v1 x1 x2
dt M1 M1 M1
dx2
v2
dt
dv2 K K 1
x2 x1 f (t )
dt M2 M2 M2
Example-7
dx1 dv1 D K K dx2
v1 v1 x1 x2 v2
dt dt M1 M1 M1 dt
dv2 K K 1
x2 x1 f (t )
dt M2 M2 M2
• In Vector-Matrix form
0 1 0 0 0
x1 K 1 x
0 0
D K
v1 M 1 M1 M1 v1
0 f (t )
x 2 0 0 0 1 x2
K 1
0 v 2
K
v2 M 0
2 M2 M2
Example-7
0 1 0 0 0
x1 K 1 x
0 0
D K
v1 M 1 M1 M1 v1
0 f (t )
x 2 0 0 0 1 x2
K 1
0 v 2
K
v2 M 0
2 M2 M2
• If x1 and v2 are the outputs of the system then
x1
1 0 0 0 v
1
y (t )
0 0 0 1
x2
v2
Example-8 (Home Work)
• Represent the translational mechanical system shown in
following figure in state space, where x3(t) is the output.
To download this lecture visit
http://imtiazhussainkalwar.weebly.com/
END OF LECTURE-9