Chapter 2
Chapter 2
d2 h i dh i
δr (t) −→ I ψ(t) + D ψ(t) = K δr (t) −→ ψ(t)
dt 2 dt
dh i h i
δr (t) −→ I r (t) + D r (t) = K δr (t) −→ r (t)
dt
Definition
Let F (s) be the Laplace transform of a function f (t), t > 0. The
following contour integral is called the inverse Laplace transform of
F (s):
Z c+j∞
−1 1
L [F (s)] ≡ f (t) = F (s)e st ds
2πj c−j∞
√
Where, j = −1 and c > σ0 .
7 5
1 F (s) = 5 F (s) =
s s2
+ 25
4 s
2 F (s) = 2 6 F (s) = 2
s s +9
2 12 2s
3 F (s) = 3 7 F (s) = 2 + 2
s s + 9 s + 81
6 5 2 1
4 F (s) = 3 + 2 + 8 F (s) = 2
s s s s +4
1
Nise (2011) page 37 to 44.
22 / 62 Professor Dr. Zobair Ibn Awal NAME 467: Control Engineering
Question
2
Given the differential equation, ddt y2 + 12 dy
dt + 32y = 32u(t), solve
for y (t) if all initial conditions are zero. Use the Laplace transform.
32
⇒Y (s) = 2 U(s)
s + 12s + 32
1
For a unit step input U(s) = . Hence,
s
32 32
Y (s) = =
s(s 2 + 12s + 32) s(s + 8)(s + 4)
32
K1 = =1
(s + 4)(s + 8)
s→0
32
K2 = = −2
s(s + 8)
s→−4
32
K3 = =1
s(s + 4)
s→−8
32 1 2 1
Y (s) = = − +
s(s + 8)(s + 4) s (s + 4) (s + 8)
Applying inverse Laplace transform results in:
Disadvantage
It requires systems to be at zero initial conditions.
Advantages
Indication on system stability
Convenience on system design
Convenience on understanding system dynamics due to
parameter changes
Convenience adding/removing system components
1 1 1
K1 = = =
s +2 0+2 2
s→0
1 1 1
K2 = = =−
s −2 2
s→−2
dψ 2 dψ
T + = K δr
dt 2 dt
dr
T + r = K δr
dt
ψ = Heading
dψ
r = Rate of turn =
dt
K = Rudder constant or Gain
T = Time constant
δr = Rudder command angle
Ψ(s) K /T
∴ G (s) = =
∆r (s) s(s + 1/T )
Aρ(dx − dy ) = qdt
Aρ(dx − dy ) = qdt
dx dy
⇒Aρ − =q
dt dt
dx dy q
⇒ − =
dt dt Aρ
dx dy P1 − P2
⇒ − = [From (2)]
dt dt RAρ
dx dy ky
⇒ − = [From (1)]
dt dt RA2 ρ
ky dy dx
⇒ 2 + =
RA ρ dt dt
ky dy dx
⇒L 2
+ =L
RA ρ dt dt
ky dy dx
⇒L + =L
RA2 ρ dt dt
k
⇒ Y (s) + sY (s) = sX (s)
RA2 ρ
k
⇒ + s Y (s) = sX (s)
RA2 ρ
Y (s) s
⇒ =
X (s) k
+s
RA2 ρ
s
∴ G (s) =
k
s+
RA2 ρ
d 2y dy
M 2
+B + Ky = Api
dt dt
d 2y
dy
⇒L M 2 + B + Ky = L [Api ]
dt dt
⇒Ms 2 Y (s) + BsY (s) + KY (s) = API (s)
Y (s) A
⇒ = 2
= G (s)
PI (s) Ms + Bs + K
PI (s) −→ G (s) −→ Y (s)
A
∴ G (s) =
Ms 2 + Bs + K
Solution: