1234567898

3.1

∞ ∞
a) 1e[ − bt
]
sin ωt = ∫ e − bt
sin ωt e dt = ∫ sin ωt e − ( s + b )t dt
− st

0 0

[− (s + b) sin ωt − ω cos ωt ]
∞

= e − ( s + b ) t 
 ( s + b ) 2 + ω2 0
ω
=
( s + b) 2 + ω2
∞ ∞
b) [
1 e − bt
]
cos ωt = ∫ e − bt
cos ωt e dt = ∫ cos ωt e − ( s + b )t dt
− st

0 0

[− (s + b) cos ωt + ω sin ωt ]
∞

= e − ( s + b ) t 
 ( s + b) 2 + ω2 0
s+b
=
( s + b) 2 + ω2

3.2

a) The Laplace transform provided is

4
Y ( s) =
s + 3s + 4 s 2 + 6 s + 4
4 3

We also know that only sin ωt is an input, where ω = 2 . Then

ω 2 2
X ( s) = = =
s +ω
2 2
s2 + 2 ( ) 2
s +2
2

Since Y(s) = D-1(s) X(s) where D(s) is the characteristic polynomial (when
all initial conditions are zero),

Solution Manual for Process Dynamics and Control, 2nd edition,

Copyright © 2004 by Dale E. Seborg, Thomas F. Edgar and Duncan A. Mellichamp.

3-1
 2 2 2
Y (s) =
( s + 3s + 2) ( s + 2)
2 2

and the original ode was

d2y dy
2
+ 3 + 2 y = 2 2 sin 2t with y ′(0) = y (0) = 0
dt dt

b) This is a unique result.

c) The solution arguments can be found from

2 2 2
Y (s) =
( s + 1)( s + 2) + ( s 2 + 2)

which in partial fraction form is

α1 α a s + a2
Y (s) = + 2 + 12
s +1 s + 2 s +2

Thus the solution will contain four functions of time

e-t , e-2t , sin 2 t , cos 2 t

3.3

a) Pulse width is obtained when x(t) = 0

Since x(t) = h – at

tω : h − atω = 0 or tω = h/a

b)
h

slope = -a slope = a
x(t) x(t)

slope = -a

x(t) = hS(t) – atS(t) + a(t -tω) S(t-tω)

3-2
 h a ae − stω h e − stω − 1
c) X ( s) = − 2 + 2 = +
s s s s s2

d) Area under pulse = h tω/2

3.4

a) f(t) = 5 S(t) – 4 S(t-2) – S(t- 6)

F (s ) = = ( 5 - 4e -2s - e -6s )
1
s

b)

x(t)

x1 a
a
x4
tr 2tr 3tr

-a -a

x2 x3

x(t) = x1(t) + x2(t) + x3(t) + x4(t)

= at – a(t − tr)S(t − tr) − a(t −2tr)S(t − 2tr) + a(t − 3tr)S(t − 3tr)

following Eq. 3-101. Thus

X(s) =
a
s 2
[
1 − e −tr s − e − 2tr s + e −3tr s ]
by utilizing the Real Translation Theorem Eq. 3-104.

3-3
3.5

55 55
T(t) = 20 S(t) + t S(t) – (t-30) S(t-30)
30 30
T (s) =
20 55 1
+
s 30 s 2
−
55 1 −30 s 20 55 1
30 s 2
e = +
s 30 s 2
1 − e −30 s ( )

3.6

s ( s + 1) α α α
a) X ( s) = = 1 + 2 + 3
( s + 2)( s + 3)( s + 4) s + 2 s + 3 s + 4

s ( s + 1)
α1 = =1
( s + 3)( s + 4) s = −2

s ( s + 1)
α2 = = −6
( s + 2)( s + 4) s = −3

s ( s + 1)
α3 = =6
( s + 2)( s + 3) s = −4

1 6 6
X (s) = − + and x(t ) = e−2t − 6e −3t + 6e −4t
s+2 s+3 s+4

s +1 s +1
b) X (s) = =
( s + 2)( s + 3)( s + 4) ( s + 2)( s + 3)( s + 2 j )( s − 2 j )
2

α1 α α + jβ 3 α 3 − j β 3
X ( s) = + 2 + 3 +
s+2 s+3 s+2j s+2j

s +1 1
α1 = =−
( s + 3)( s 2 + 4) s = −2
8
s +1 2
α2 = =
( s + 2)( s 2 + 4) s = −3
13
s +1 1− 2 j − 3 + 11 j
α 3 + jβ 3 = = =
( s + 2)( s + 3)( s − 2 j ) s = −2 j
− 40 − 8 j 208

3-4
 1 2  −3   11 
x(t ) = − e − 2t + e −3t + 2  cos 2t + 2  sin 2t
8 13  208   208 

1 2 3 11
= − e − 2t + e −3t − cos 2t + sin 2t
8 13 104 104

s+4 α α2
c) X ( s) = = 1 + (1)
( s + 1) 2
s + 1 ( s + 1) 2

α 2 = ( s + 4) s =−1 = 3

In Eq. 1, substitute any s≠-1 to determine α1. Arbitrarily using s=0, Eq. 1
gives

4 α1 3
= + or α1 = 1
12 1 12

1 3
X (s) = + and x( t ) = e− t + 3te − t
s + 1 ( s + 1) 2

1 1 1
d) X (s) = 2 = =
s + s +1  1  3 (s + b ) + ω
2 2 2

s +  +
 2 4
1 3
where b = and ω =
2 2

t
1 2 −2 3
x(t ) = e −bt sin ω t = e sin t
ω 3 2

s +1
e) X(s) = e −0.5 s
s ( s + 2)( s + 3)

To invert, we first ignore the time delay term. Using the Heaviside
expansion with the partial fraction expansion,

s +1 A B C
Xˆ ( s ) = = + +
s ( s + 2)( s + 3) s s + 2 s + 3

Multiply by s and let s → 0

3-5
 1 1
A= =
(2)(3) 6
Multiply by (s+2) and let s→ −2

− 2 +1 −1 1
B= = =
(−2)(−2 + 3) (−2)(1) 2

Multiply by (s+3) and let s→-3

− 3 +1 −2 2
C= = =−
(−3)(−3 + 2) (−3)(−1) 3

Then

1 6 1 2 −2 3
Xˆ ( s ) = + +
s s+2 s+3

1 1 − 2 t 2 −3t
xˆ (t ) = + e − e
6 2 3

Imposing shift theorem

1 1 − 2 (t −0.5) 2 −3( t − 0.5)

x(t ) = xˆ (t − 0.5) = + e − e
6 2 3

for t ≥ 0.5

3.7

6( s + 1) 6 α α
a) Y (s) = = 2 = 1 + 22
s ( s + 1) s
2
s s

6
α2 = s2 =6 α1 = 0
s2 s =0

6
Y (s) = 2
s

12( s + 2) α1 α 2 s + α 3
b) Y (s) = = + 2
s ( s 2 + 9) s s +9

3-6
 Multiplying both sides by s(s2+9)

12( s + 2) = α 1 ( s 2 + 9) + (α 2 s + α 3 )( s ) or
12 s + 24 = (α 1 + α 2 ) s + α 3 s + 9α 12

Equating coefficients of like powers of s,

s2: α1 + α2 = 0
s1: α3 = 12
0
s : 9α1 = 24

Solving simultaneously,

8 −8
α1 = , α2 = , α 3 = 12
3 3
 8 
 − s + 12 
8 1  3 
Y (s) = +
3s s2 + 9

( s + 2)( s + 3) α α α
c) Y (s) = = 1 + 2 + 3
( s + 4)( s + 5)( s + 6) s + 4 s + 5 s + 6

( s + 2)( s + 3)
α1 = =1
( s + 5)( s + 6) s = −4

( s + 2)( s + 3)
α2 = = −6
( s + 4)( s + 6) s = −5

( s + 2)( s + 3)
α3 = =6
( s + 4)( s + 5) s = −6

1 6 6
Y (s) = − +
s+4 s+5 s+6

1 1
Y (s) = =
[(s + 1) ]
d)
2 2
+ 1 ( s + 2) ( s + 2s + 2) 2 ( s + 2)
2

α1 s + α 2 α s + α4 α
= + 2 3 + 5
s + 2s + 2 ( s + 2s + 2)
2 2
s+2

3-7
 Multiplying both sides by ( s 2 + 2s + 2) 2 ( s + 2) gives

1 = α1s4 + 4α1s3 + 6α1s2 +4α1s + α2s3 +4α2s2 +6α2s +4α2 + α3s2 +2α3s +
α4s + 2α4 + α5s4 + 4α5s3 + 8α5s2 + 8α5s + 4α5

Equating coefficients of like power of s,

s4 : α1 + α5 = 0

s3 : 4α1 + α2 + 4α5 = 0

s2 : 6α1 + 4α2 + α3 + 8α5 = 0

s1 : 4α1 + 6α2 + 2α3 + α4 + 8α5 = 0

s0 : 4α2 + 2α4 + 4α5 = 1

Solving simultaneously:

α1 = -1/4 α2 = 0 α3=-1/2 α4=0 α5 = ¼

− 1 / 4s − 1 / 2s 1/ 4
Y (s) = + 2 +
s + 2 s + 2 ( s + 2 s + 2)
2 2
s+2

3.8

a) From Eq. 3-100

t  1
1  ∫ f (t * )dt *  = F ( s )
0  s
 t −τ  1
we know that 1  ∫ e dτ = 21 e − τ =[ ] 1
s ( s + 1)
0  s

∴ Laplace transforming yields

2
s2X(s) + 3X(s) + 2X(s) =
s ( s + 1)
2
or (s2 + 3s + 1) X(s) =
s ( s + 1)

3-8
 2
X(s) =
s ( s + 1) 2 ( s + 2)

and x(t) = 1 − 2te-t − e-2t

b) Applying the final Value Theorem

2
lim x(t ) = lim sX ( s ) = lim =2
t →∞ s →0 s →0 ( s + 1) ( s + 2)
2

[ Note that Final Value Theorem is applicable here]

3.9

6( s + 2) 6( s + 2)
a) X (s) = =
( s + 9 s + 20)( s + 4) ( s + 4)( s + 5)( s + 4)
2

 6s ( s + 2) 
x(0) = lim  =0
s →∞ ( s + 5)( s + 4) 2 
 
 6s ( s + 2) 
x(∞) = lim =0
s →0 ( s + 5)( s + 4) 2 
 

x(t) is converging (or bounded) because [sX(s)] does not have a limit at
s = −4, and s = −5 only, i.e., it has a limit for all real values of s ≥ 0.

x(t) is smooth because the denominator of [sX(s)] is a product of real

factors only. See Fig. S3.9a.

10s 2 − 3 10s 2 − 3
b) X (s) = 2 =
( s − 6s + 10)( s + 2) ( s − 3 + 2 j ) ( s − 3 − 2 j )( s + 2)

 10s 3 − 3s 
x(0) = lim  2  = 10
s →∞ ( s − 6 s + 10)( s + 2)
 

Application of final value theorem is not valid because [sX(s)] does not
have a limit for some real s ≥ 0, i.e., at s = 3±2j. For the same reason, x(t)
is diverging (unbounded).

x(t) is oscillatory because the denominator of [sX(s)] includes complex

factors. See Fig. S3.9b.

3-9
 16s + 5 16 s + 5
c) X (s) = =
( s + 9) ( s + 3 j ) ( s − 3 j )
2

16s 2 + 5s 
x(0) = lim  2  = 16
s →∞
 ( s + 9) 

Application of final value theorem is not valid because [sX(s)] does not
have a limit for real s = 0. This implies that x(t) is not diverging, since
divergence occurs only if [sX(s)] does not have a limit for some real value
of s>0.

x(t) is oscillatory because the denominator of [sX(s)] is a product of

complex factors. Since x(t) is oscillatory, it is not converging either. See
Fig. S3.9c
0.4

0.35

0.3

0.25

0.2
x(t)

0.15

0.1

0.05

-0.05
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
Time

Figure S3.9a. Simulation of X(s) for case a)

8000

6000

4000

2000

0
x(t)

-2000

-4000

-6000

-8000

-10000

-12000
0 0.5 1 1.5 2 2.5
Time

Figure S3.9b. Simulation of X(s) for case b)

3-10
 20

15

10

x(t)
0

-5

-10

-15

-20
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
Time

Figure S3.9c. Simulation of X(s) for case c)

The Simulink block diagram is shown below. An impulse input should be

used to obtain the function’s behavior. In this case note that the impulse
input is simulated by a rectangular pulse input of very short duration. (At
time t = 0 and t =0.001 with changes of magnitude 1000 and –1000
respectively). The MATLAB command impulse might also be used.

Figure S3.9d. Simulink block diagram for cases a), b) and c).

3-11
3.10

a)
2 2 A B C
i) Y(s) = = 2 = 2+ +
s ( s + 4 s ) s ( s + 4) s
2
s s+4

∴ y(t) will contain terms of form: constant, t, e-4t

2 2 A B C
ii) Y(s) = = = + +
s ( s + 4 s + 3) s ( s + 1)( s + 3) s s + 1 s + 3
2

∴ y(t) will contain terms of form: constant, e-t, e-3t

2 2 A B C
iii) Y (s) = = = + +
s ( s + 4 s + 4) s ( s + 2)
2 2
s ( s + 2) 2
s+2

∴ y(t) will contain terms of form: constant, e-2t , te-2t

2
iv) Y (s) =
s ( s + 4 s + 8)
2

s 2 + 4 s + 8 = ( s 2 + 4 s + 4) + (8 − 4) = ( s + 2) 2 + 2 2
2
Y (s) =
s[(s + 2) 2 + 2 2 ]

∴ y(t) will contain terms of form: constant, e-2t sin2t, e-2tcos2t

2( s + 1) 2( s + 1) A Bs C
b) Y (s) = = = + 2 + 2
s ( s + 4) s ( s + 2 ) s s + 2
2 2 2 2
s + 22
2( s + 1) 1
A = lim 2 =
s →0 ( s + 4) 2

2(s+1) = A(s2+4) + Bs(s) + Cs

2s+2 = As2 + 4A + Bs2 + Cs

Equating coefficients on like powers of s

1
s2 : 0=A+B → B = −A = −
2
s1 : 2= C → C=2
1
s0 : 2 = 4A → A=
2

3-12
 1 2 − (1 2) s 2
∴ Y(s) = + 2 + 2
s s +2 2
s + 22

1 1 2
y(t) = − cos 2t + sin 2t
2 2 2

1
y(t) = (1 − cos 2t ) + sin 2t
2

3.11

Since convergent and oscillatory behavior does not depend on initial

dx 2 (0) dx(0)
conditions, assume = = x(0) = 0
dt 2 dt
a) Laplace transform of the equation gives

3
s 3 X ( s ) + 2 s 2 X ( s ) + 2 sX ( s ) + X ( s ) =
s
3 3
X (s) = =
s ( s + 2 s 2 + 2 s + 1)
3
1 3 1 3
s ( s + 1)( s + + j )( s + − j)
2 2 2 2
Denominator of [sX(s)] contains complex factors so that x(t) is oscillatory,
and denominator vanishes at real values of s= −1 and -½ which are all <0
so that x(t) is convergent. See Fig. S3.11a.
2
b) s 2 X ( s) − X (s) =
s −1
2 2
X (s) = =
( s − 1)( s − 1) ( s − 1) 2 ( s + 1)
2

The denominator contains no complex factors; x(t) is not oscillatory.

The denominator vanishes at s=1 ≥0; x(t) is divergent. See Fig. S3.11b.

1
c) s 3 X ( s) + X (s) =
s +1
2

1 1
X (s) = =
( s + 1)( s 3 + 1)
2
1 3 1 3
( s + j )( s − j )( s + 1)( s − + j )( s − − j)
2 2 2 2

The denominator contains complex factors; x(t) is oscillatory.

The denominator vanishes at real s = 0, ½; x(t) is not convergent. See Fig.
S3.11c.

3-13
 4
d) s 2 X ( s ) + sX ( s ) =
s
4 4
X (s) = = 2
s ( s + s ) s ( s + 1)
2

The denominator of [sX(s)] contains no complex factors; x(t) is not

oscillatory.
The denominator of [sX(s)] vanishes at s = 0; x(t) is not convergent. See
Fig. S3.11d.

3.5

2.5

2
x(t)

1.5

0.5

-0.5
0 1 2 3 4 5 6 7 8 9 10
time

Figure S3.11a. Simulation of X(s) for case a)

700

600

500

400
x(t)

300

200

100

0
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
time

Figure S3.11b. Simulation of X(s) for case b)

3-14
 80

60

40

x(t)
20

-20

-40
0 1 2 3 4 5 6 7 8 9 10
time

Figure S3.11c. Simulation of X(s) for case c)

18

16

14

12

10
x(t)

0
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5

time

Figure S3.11d. Simulation of X(s) for case d)

3.12

Since the time function in the solution is not a function of initial

conditions, we Laplace Transform with

dx(0)
x ( 0) = =0
dt

τ1τ2s2X(s) + (τ1+τ2)sX(s) + X(s) = KU(s)

3-15
 K
X ( s) = U (s)
τ1 τ 2 s + ( τ1 + τ 2 ) s + 1
2

Factoring denominator

K
X ( s) = U (s)
(τ1 s + 1)(τ 2 s + 1)

a
a) If u(t) = a S(t) then U(s)=
s

Ka
X a (s) = τ1 ≠ τ 2
s (τ1 s + 1)(τ 2 s + 1)

xa(t) = fa( S(t), e-t/τ1, e– t/τ2)

bτ
b) If u(t) = be-t/τ then U(s) =
τs +1

Kbτ
X b ( s) = τ ≠ τ1 ≠ τ 2
(τs + 1)(τ1 s + 1)(τ 2 s + 1)

xb(t) = fb(e-t/τ , e-t/τ1, e– t/τ2)

τc
c) If u(t) =ce-t/τ where τ = τ1 , then U(s) =
τ1 s + 1
Kcτ
X c ( s) =
(τ1 s + 1) 2 (τ 2 s + 1)

xc(t) = fc(e– t/τ1, t e– t/τ1, e– t/τ2)

dω
d) If u(t) = d sin ωt then U(s) =
s + ω2
2

Kd
X d (s) =
( s + ω )(τ1 s + 1)(τ 2 s + 1)
2 2

xd(t) = fd(e– t/τ1, e– t/τ2, sin ωt, cos ωt)

3-16
3.13

dx3 d 2 x(0) dx(0)

a) 3
+ 4 x = et with = = x(0) = 0
dt dt 2 dt

Laplace transform of the equation,

1
s 3 X(s) + 4 X(s) =
s −1

1 1
X (s) = =
( s − 1)( s + 4) ( s − 1)( s + 1.59)( s − 0.79 + 1.37 j )( s − 0.79 − 1.37 j )
3

α1 α2 α 3 + jβ 3 α 3 − jβ 3
= + + +
s − 1 s + 1.59 s − 0.79 + 1.37 j s − 0.79 − 1.37 j

1 1
α1 = =
( s + 4)
3
s =1
5
1 1
α2 = =−
( s − 1)( s − 0.79 + 1.37 j )( s − 0.79 − 1.37 j ) s = −1.59
19.6
1
α 3 + jβ 3 = = −0.74 − 0.59 j
( s − 1)( s + 1.59)( s − 0.79 − 1.37 j ) s = 0.79 −1.37 j

1 1
−
− 0.074 − 0.059 j − 0.074 + 0.059 j
X(s) = 5 + 19.6 + +
s − 1 s + 1.59 s − 0.79 + 1.37 j s − 0.79 − 1.37 j

1 1 −1.59t
x(t ) = e t − e − 2e 0.79t (0.074 cos 1.37t + 0.059 sin 1.37t )
5 19.6
dx
b) − 12 x = sin 3t with x(0) = 0
dt

3
sX (s) − 12 X(s) =
s +92

3 3
X (s) = =
( s + 9)( s − 12) ( s + 3 j )( s − 3 j )( s − 12)
2

α 1 + j β1 α 1 − j β 1 α3
= + +
s+3j s−3j s − 12

3-17
 3 3 1 4
α 1 + j β1 = = =− − j
( s − 3 j )( s − 12) s = −3 j
− 18 + 72 j 102 102

3 1
α3 = =
( s + 9 ) s =12 51
2

1 4 1 4 1
− − j − + j
X (s) = 102 102 + 102 102 + 51
s+3j s −3j s − 12

1 1
x(t ) = − (cos 3t + 4 sin 3t ) + e12 t
51 51

d2x dx dx(0)
c) 2
+ 6 + 25 x = e −t with = x(0) = 0
dt dt dt

1 1
s 2 X ( s ) + 6 sX ( s ) + 25X ( s ) + X ( s ) = or X( s ) =
s +1 ( s + 1 )( s + 6 s + 25 )
2

1 α α + β2 j α 2 − β2 j
X ( s) = = 1 + 2 +
( s + 1)( s + 3 + 4 j )( s + 3 − 4 j ) s + 1 s + 3 + 4 j s + 3 − 4 j
1 1
α1 = =
( s + 6 s + 25)
2
s = −1
20

1 1 1
α 2 + jβ 2 = =− − j
( s + 1)( s + 3 − 4 j ) s = −3− 4 j
40 80

1 1 1 1 1
− − j − − j
X ( s ) = 20 + 40 80 + 40 80
s +1 s + 3 + 4 j s + 3− 4 j

1 −t 1 1
x(t ) = e − e −3t ( cos 4t + sin 4t )
20 20 40

d) Laplace transforming (assuming initial conditions = 0, since they do not

affect results)

sY1(s) + Y2(s) = X1(s) (1)

sY2(s) – 2Y1(s) + 3 Y2(s) = X2(s) (2)

3-18
From (2),

(s+3) Y2(s) = X2(s) + 2Y1(s)

1 2
Y2 ( s ) = X 2 (s) + Y1 ( s )
s+3 s+3

Substitute in Eq.1

1 2
sY1(s) + X 2 (s) + Y1 ( s ) = X1(s)
s+3 s+3

We neglect X2(s) since it is equal to zero.

[s(s + 3) + 2]Y1 (s) = ( s + 3) X 1 (s)

( s 2 + 3s + 2)Y1 ( s ) = ( s + 3) X 1 ( s )

s+3 s+3
Y1 ( s ) = X 1 (s) = X 1 (s)
s + 3s + 2
2
( s + 1)( s + 2)

1
Now if x1(t) = e-t then X1(s) =
s +1
s+3 A B C
∴ Y1 ( s ) = = + +
( s + 1) ( s + 2) ( s + 1) ( s + 1)
2 2
s+2

so that y1(t) will contain e-t/τ , te-t/τ, e–2t functions of time.

For Y2(s)

2 A B C
Y2 ( s ) = = + +
( s + 1) ( s + 2) ( s + 1) ( s + 1)
2 2
s+2

so that y2(t) will contain the same functions of time as y1(t) (although
different coefficients).

3-19
3.14

d 2 y (t ) dy (t ) d ( x − 2)
2
+3 + y (t ) = 4 − x(t − 2)
dt dt dt

Taking the Laplace transform and assuming zero initial conditions,

s2Y(s) + 3sY(s) + Y(s) = 4 e-2ssX(s) −e-2sX(s)

Rearranging,

Y (s) − (1 − 4s )e −2 s
= G(s) = 2 (1)
X ( s) s + 3s + 1

a) The standard form of the denominator is : τ2s2 + 2ζτs + 1

From (1) , τ = 1 , ζ = 1.5

Thus the system will exhibit overdamped and non-oscillatory response.

b) Steady-state gain

K = lim G ( s ) = −1 (from (1))

s →0

c) For a step change in x

1.5 − (1 − 4 s )e −2 s 1.5
X(s) = and Y(s) = 2
s ( s + 3s + 1) s

Therefore yˆ (t ) = −1.5 + 1.5e-1.5t cosh(1.11t) + 7.38e-1.5t sinh(1.11t)

Using MATLAB-Simulink, y(t)= yˆ (t − 2) is shown in Fig. S3.14
1.5

0.5

-0.5

-1

-1.5
0 5 10 15 20 25 30

Figure S3.14. Output variable for a step change in x of magnitude 1.5

3-20
3.15

f (t ) = hS (t ) − hS (t − 1 / h)

+ 4 x = h[S (t ) − S (t − 1 / h)]
dx
, x(0)=0
dt

Take Laplace transform,

 1 e−s / h 
sX ( s ) + 4 X ( s ) = h − 
s s 
1 α α 
X ( s ) = h(1 − e − s / h ) = h(1 − e − s / h )  1 + 2 
s ( s + 4)  s s + 4
1 1 1 1
α1 = = , α2 = =−
s + 4 s =0 4 s s =−4 4

h 1 1 
X ( s) = (1 − e − s / h )  − 
4 s s + 4

h 1 e −s / h 1 e −s / h 
=  − − + 
4 s s s + 4 s + 4

0 t <0
h
x(t ) = (1 − e − 4t ) 0 < t < 1/h
4

4
[
h − 4( t −1 / h )
e − e −4t ] t > 1/h

1
h=1
h=10
0.9 h=100

0.8

0.7

0.6
x(t)

0.5

0.4

0.3

0.2

0.1

0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
time

Figure S3.15. Solution for values h= 1, 10 and 100

3-21
3.16

a) Laplace transforming

[s Y (s) − sy(0) − y′(0)]+ 6[sY (s) − y(0)] + 9Y (s) = s s+ 1

2
2

s
(s2 + 6s + 9)Y(s) − s(1) − 2 –(6)(1)=
s +1
2

s
(s2 + 6s + 9)Y(s) = +s+8
s +1
2

2 s + s 3 + s + 8s 2 + 8
(s + 6s + 9)Y(s) =
s2 +1

s 3 + 8s 2 + 2 s + 8
Y(s) =
( s + 3) 2 ( s 2 + 1)

To find y(t) we have to expand Y(s) into its partial fractions

A B Cs D
Y (s) = + + 2 + 2
( s + 3) 2
s + 3 s +1 s +1

y(t) = Ate-3t + Be-3t + C cost + D sint

s +1
b) Y(s)=
s ( s + 4 s + 8)
2

42
Since < 8 we know we will have complex factors.
4

∴ complete square in denominator

s2 + 4s + 8 = s2 + 4s + 4 + 8−4

= s2 + 4s + 4 + 4 = (s+2)2 + (2)2 { b = 2 , ω=2}

∴ Partial fraction expansion gives

3-22
 A B ( s + 2) C s +1
Y (s) = + 2 + 2 =
s s + 4 s + 8 s + 4 s + 8 s ( s + 4 s + 8)
2

Multiply by s and let s→0

A=1/8

Multiply by s(s2+4s+8)

A(s2+4s+8) + B(s+2)s + Cs = s + 1

As2 + 4As + 8A + Bs2 + 2Bs + Cs = s + 1

1
s2 : A + B = 0 → B = −A = −
8
1 1 3
s1 : 4A + 2B + C = 1 → C = 1 + 2  − 4  =
8 8 4
1
s0 : 8A = 1 →A= (This checks with above result)
8

1 / 8 (− 1 / 8)( s + 2) 3/ 4
Y (s) = + +
s ( s + 2) + 2
2 2
( s + 2) 2 + 2 2

1 1  3
y(t) =   −   e-2t cos 2t +   e-2t sin 2t
8 8 8

3.17

dC
V + qC = qCi
dt

Since V and q are constant, we can Laplace Transform

sVC(s) + qC(s) = q Ci(s)

Note that c(t = 0) = 0

Also, ci(t) = 0 , t ≤0
ci(t) = ci , t>0

3-23
Laplace transforming the input function, a constant,

ci
Ci ( s ) =
s
so that
ci qci
sVC(s) + qC(s) = q or C(s) =
s ( sV + q ) s

Dividing numerator and denominator by q

ci
C(s) =
V 
 s +1 s
q 

Use Transform pair #3 in Table 3.1 to invert (τ =V/q)

 − t 
V

c(t) = ci 1 − e q 

 
 

Using MATLAB, the concentration response is shown in Fig. S3.17.

(Consider V = 2 m3, Ci=50 Kg/m3 and q = 0.4 m3/min)

50

45

40

35

30
c(t)

25

20

15

10

0
0 5 10 15 20 25 30

Time

Figure S3.17. Concentration response of the reactor effluent stream.

3-24
3.18

KAω
a) If Y(s) =
s( s 2 + ω2 )
Aω
and input U(s) = = 1 {A sin ωt}
(s + ω2 )
2

then the differential equation had to be

dy
= Ku (t ) with y(0) = 0
dt

KAω α α s α ω
b) Y(s) = = 1+ 2 2 2 + 2 3 2
s( s + ω ) s s + ω
2 2
s +ω

KAω KA
α1 = =
s + ω2
2
s →0 ω

Find α2 and α3 by equating coefficients

KAω= α1(s2+ω2) + α2s2+α3ωs

KAω = α1s2 + α1ω2 + α2s2 + α3ωs

− KA
s2 : 0 = α1 + α2 → α2 = −α1 =
ω

s: 0 = α3 ω → α3 = 0

KAω KA / ω ( KA / ω) s
∴ Y(s) = = − 2
s( s + ω )
2 2
s s + ω2
KA
y(t) = (1 − cos ωt )
ω

3-25
 c)

u(t)

-A

2KA/ω

y(t)

Time

i) We see that y(t) follows behind u(t) by 1/4 cycle = 2π/4= π/2 rad.
which is constant for all ω

ii) The amplitudes of the two sinusoidal quantities are:

y : KA/ω
u: A

Thus their ratio is K/ω, which is a function of frequency.

3-26