Lectures 1 - ME 340
Lectures 1 - ME 340
.
y0
c d
y
* Complex numbers are linear transformations represented by
matrices of the form
!
r cos q
r sin q
,
r sin q r cos q
by the corresponding rectangular form r cos q + jr sin q or by
the polar form r ejj , where the imaginary unit j satisfies the
condition that j j = 1. A complex number may be described
in terms of its real and imaginary parts r cos q and r sin q, or
in terms of its magnitude r and phase j.
* An eigenvalue l of the matrix
a
c
b
d
( x 0 , y0 ) = ax + by, cx + dy ,
(1)
y0 = 3
( x, y) = (5, 1)
2y = 3
( x 0 , y0 ) = (1, 2)
( x 0 , y0 )
x0
y0
:=
.
y0
c d
y
cx + dy
(2)
(3)
b
d
(4)
:=
.
(6)
g h
c d
ga + hc gb + hd
( x 0 , y0 ) = ( ax + by, cx + dy)
( x, y)
=
.
(7)
c d
0
0
If ad
with
1
2
(8)
b
d
:= ad
bc
(9)
(10)
equals
a
c
b
d
d/( ad bc)
c/( ad bc)
:=
b/( ad bc)
a/( ad bc)
(11)
The corresponding transformation reverses the changes to the coordinates, since the product
!
!
!
d/( ad bc)
b/( ad bc)
a b
1 0
=
(12)
c/( ad bc) a/( ad bc)
c d
0 1
equals the identity matrix, for which ( x 0 , y0 ) = ( x, y).
( x20
x10 )2 + (y20
= ( x2 cos j
y10 )2
y2 sin j
+ ( x2 sin j + y2 cos j
= ( x2
x1 )2 + ( y2
x1 sin j
y1 cos j)2
y1 )2 ,
(14)
y cos q = 0
(15)
y0 cos(q + j) = 0,
(16)
since
cos j sin(q + j)
(17)
(18)
Since the origin is transformed to the origin, the distance from each
such point to the origin is unchanged by the transformation.
j.
Exercises
1. Determine the point ( x, y) that is transformed to ( x 0 , y0 ) = (1, 5)
under the linear coordinate transformation
( x 0 , y0 ) = (4x
y, x + 3y).
y.
1
0
y0
3 0
1
2
1
2 1
x
y
0 2
0
0
1
3
1
0
1
2
1
2
Solutions
1. With ( x 0 , y0 ) = (1, 5) it follows that
1 = 4x
y, 5 = x + 3y,
1
1
y0
1 1
0
1
x
y
x
x+y
3. The equations
x 0 = x + 2y, y0 = 3x
1 0
1
( x + 2y0 ), y = (3x 0
7
7
y 0 ),
3y = 2 then
5x 0 + 11y0 = 14
4x
3y = 2
( x 0 , y0 ) = ( 5, 1)
( x, y) = ( 1, 2)
1
0
4
3
1
0
0
1
1
0
2
3
6. Multiplication yields
1
2
0
2
1
1
y, 1 = 2x + y
x
y
0
2
8. The determinant
equals 0 if a =
implies that 2x
1
2
a
3
0
1
1/2
0
= 3 + 2a
2
3
1
2
x
y
0
0
3y = 0.
0 2
0
Since
1
1
1
0
1
2
0
1
3
1
2
1
0
1
0
1
2
1
2
1
2
1
0
1
12
1
3
0
0
1
it follows that
1
2
1
12
1
3
!n
1
0
1
2
1
2n
0
1
3n
1
0
1
2
1
2
0
a
(20)
1
0
(22)
ab
= a
bx
by
abx
aby
x
y
. (25)
x
y
(26)
. (27)
x
y
= (ab)
x
y
= ja
bx
by
aby
abx
abx
aby
= j(ab)
x
y
= ( ab)
(24)
In particular, j j = 1, since application of the corresponding transformation twice in succession is equivalent to the transformation
( x 00 , y00 ) = ( y0 , x 0 ) = ( x, y) represented by the matrix product
!
!
!
0
1
0
1
1 0
=
.
(28)
1 0
1 0
0
1
Example 5. If matrix addition is defined so that
!
!
1
2
1 0
0
=
+
2 1
0 1
2
then this corresponds to the linear transformation
!
!
!
x0
1 0
x
0
=
+
0
y
0 1
y
2
2
0
2
0
(29)
x
y
(30)
10
x
y
x
y
+ j2
= (1 + j2)
x
y
11
(31)
=
(32)
2 1
2 1
4
3
or, equivalently,
(1 + j2) (1 + j2) = 1 1 + 1 j2 + j2 1 + j2 j2 =
since j2 j2 =
3 + j4,
(33)
4.
(34)
(35)
with r > 0 rotates every straight line through the origin by j about the
origin, and stretches the result by r.
Figure 4: The linear coordinate transformation in (35) rotates straight lines
through the origin by j and stretches
the result by r.
(r cos q, r sin q )
j
q
rr
y0
r sin j r cos j
r sin q
!
rr (cos j cos q sin j sin q )
=
rr (sin j cos q + cos j sin q )
!
rr cos( j + q )
=
,
(36)
rr sin( j + q )
where the last equality follows from the trigonometric addition formulae.
The positive factor r and the angle j are the magnitude and phase
of the complex number (35). The rectangular form
r cos j + jr sin j
(37)
(38)
Example 7. Application of the transformation represented by 1 + j2 followed by the transformation represented by 1 j2 is equivalent to a single
transformation represented by
(1
j2) (1 + j2) = 1 1 + 1 j2
j2 1
j2 j2 = 5.
(39)
It follows that application of the transformation represented by 1 + j2 followed by the transformation represented by 15 j 25 is equivalent to a single
transformation represented by
1
2
j
(1 + j2) = 1,
(40)
5
5
i.e., the identity transformation.
For an arbitrary complex number represented by a matrix
!
a
b
b a
(41)
b
a
!T
a
b
b
a
(42)
12
(a + jb) a + jb = (a + jb) (a
(43)
b
a
a/(a2 + b2 )
b/(a2 + b2 )
b/(a2 + b2 )
a/(a2 + b2 )
1
a jb
a + jb
a + jb
= 2
=
=
.
a + jb
a + b2
|a + jb|2
a + jb (a + jb)
(44)
(45)
Exercises
1. Compute
3
1
j2
(a + jb) (g + jd) = ab
bd + j(ad + bg)
b
a
g
d
d
g
4. Show that
b
a
g
d
d
g
5. Show that
(2
j3) (3 + j2).
13
j).
1 + j2.
j.
j)
1
.
j+1
1 j
2 + j2
1+j
.
2 j2
Solutions
1. By definition
j2
3
1
2
6
14
3. Here,
a
b
b
a
g
d
d
g
a
b
b
a
g
d
d
g
g
d
d
g
a
b
b
a
ag bd
ad + bg
ad
ag
bg
bd
ag bd
ad + bg
ad
ag
bg
bd
ag bd
ad + bg
ad
ag
bg
bd
4. Here,
and
.
5. Here,
(a + jb) (g + jd) = ag
bd
j(ad + bg),
(a + jb) (g + jd) = (a
jb) (g
whereas
= ag
jd)
j(ad + bg).
bd
(2
j3) (3 + j2) = 2 3 + 2 j2
j3 3
j3 j2 = 12
j5.
j) = j 1
j j = j + 1.
8. By substitution
(1
j2) (1
9. Here, j3 =
10. Here, |
j2)
2(1
j2) + 5 = 1
2 + j4 + 5 = 0.
j3.
1 + j2|2 = ( 1 + j2) ( 1
j4
1 + j2| =
17.
j) = (3 + j)/10.
12. Here,
(1
j)
1
= (1
j+1
j)
j
2
j.
15
j3.
p
15. Since cos p/3 = 1/2 and sin p/3 = 3/2,
!
!
!
1
p
1
p
x0
cos
sin
1
2
3
2
3
=
=
1
p
1
p
y0
1
2 sin 3
2 cos 3
4p
1+ 3
4
16. The difference between a complex number and its complex conjugate equals j2 multiplied by its imaginary part. Since
1 j
1
= (1
2 + j2
8
the imaginary part is
j)(2
j2) =
1
j ,
2
1
1
j2
= 1.
j
2
=
.
(47)
y0
1 2
y
x + 2y
In particular, points ( x, y) = (g, g) for any number g are transformed
to points ( x 0 , y0 ) = (g, g). It follows that the straight line y = x is
transformed to the straight line y0 = x 0 .
Similarly, points ( x, y) = (g, g) for any number g are transformed
to points ( x 0 , y0 ) = (3g, 3g). It follows that the straight line y = x is
transformed to the straight line y0 = x 0 .
A number l is a real eigenvalue of
a
c
b
d
(48)
16
= (a
u)(d
bc = 0,
u)
(49)
b
d
g
h
g
h
0
l
l
0
g
h
= l
(51)
2
1
1
2
l
1
1
2 l
g
h
(52)
(53)
g
h
= l
g
h
0
0
(55)
l)
(56)
(54)
l
1
1
2 l
l2
= (2
l)(2
4l + 3 = 0,
(57)
h = 0,
g + h = 0 ) h = g,
(58)
(59)
17
h = 0,
h=0)h=
(60)
g,
(61)
(62)
u
c
b
d
= (a
u)(d
u)
bc = 0.
(63)
b
d
gR
hR
gI
hI
gR
hR
gI
hI
The expression
gR + jg I
hR + jh I
(64)
lR
lI
lI
lR
(65)
(66)
(a
u )2 + b2 = 0
(67)
(68)
b2 . The equations
g I b,
(69)
bh I = gR b + g I a
(70)
bhR = gR a
h I b,
(71)
bg I + ah I = hR b + h I a
(72)
bgR + ahR = hR a
18
a
c
b
d
(74)
(75)
there exist angles j and y, and scaling factors r x and ry , such that the
transformation may be represented by the product
!
!
!
cos j
sin j
rx 0
cos y
sin y
,
(76)
sin j
cos j
0 ry
sin y
cos y
i.e., as a rotation by y, followed by a scaling of each of the coordinates by r x and ry , respectively, followed by a rotation by j. Here,
r x < 0 and ry < 0 correspond to reflection about the y and x axis,
respectively.
Example 11. Consider the linear coordinate transformation represented by
the matrix
!
2
1
p
.
(77)
0
3
This equals the product (76) provided that
r x cos j cos y
ry sin j sin y =
r x cos j sin y
ry sin j cos y = 1,
2,
(78)
(79)
ry
rx
, ry =
sin j sin y
sin j sin y
(80)
3.
(81)
(82)
19
then yields
r x cot j cot y
r x cot j
2,
ry =
(83)
ry cot y = 1,
r x cot y + ry cos j = 0,
r x + ry cot j cot y =
(84)
(85)
3.
(86)
r2x
rx
r2y
, cot y =
ry
r2x
r2y
(88)
p
r x ry
r x ry
2r x
3ry = r2y r2x cot j cot y = 2
= p
(89)
2
r x ry
3r x 2ry
or, in other words,
p
2 3r2x + 8r x ry
p
p
ry
2 3r2y = 0 ) r x = 3ry or p
3
(90)
2 2
r =
3 y
ry .
(91)
p
In the former case, cot j =
3 and cot y = 1, which imply that j =
p/6 + np and y = p/4 + mp, for arbitrary integers n and m and,
p
p
consequently, r x = ( 1)m+n 6 and ry = ( 1)m+n 2. In the latter
p
case, cot j = 1/ 3 and cot y = 1, which imply that j = p/3 + np
and y = p/4 + mp, for arbitrary integers n and m and, consequently,
p
p
r x = ( 1)m+n 2 and ry = ( 1)m+n 6.
The angles j and y, as well as the scaling factors r x and ry may be
obtained from the real eigenvalues and eigenvectors of
!
!
!
!
a b
a b
a c
a c
.
(93)
and
c d
b d
b d
c d
20
In both cases, the eigenvalues equal r2x and ry2 . The corresponding
eigenvectors are multiples of
!
!
cos j
sin j
,
(94)
sin j
cos j
and
cos y
sin y
sin y
cos y
(95)
respectively.
Example 12. In the case of the linear transformation represented by
!
2
1
p
,
0
3
(96)
the eigenvalues of
2
0
1
p
2
1
0
p
5
p
p !
3
3
(97)
(98)
(99)
respectively. Here,
cos j =
3h, sin j = h ) j =
1
3
2
0
1
p
p
+ np
6
4
2
2
4
(100)
(101)
(102)
(103)
respectively. Here,
sin y = g, cos y = g ) y =
p
+ mp
4
(104)
21
22
Exercises
1. Find the characteristic equation for the matrix
!
0 3
.
1
1
2. Compute all eigenvalues of the matrix
!
2 3/2
.
0
1
3. Verify that
1
0
is a real eigenvector of
1
0
5/2
4
1
3
G
F
D
Y
D
Y
s
t
equals
l
s
t
:=
:=
Gs+Dt
Fs+Yt
ls
lt
l
F
:= (G
l) (Y
l)
F D = 0.
23
Solutions
1. The characteristic equation equals
0
3
1 u
u
1
= 0 ) u2 + u
2 = 0. Since
p
b b2 4ac
,
2a
au2 + bu + c = 0 ) u =
it follows that the eigenvalues equal
p
1 9
= 2 and
2
3. The claim follows since
!
1 5/2
0
4
1
0
1.
1
0
3 = 0.
1
0
5
1
hR h I
hR h I
1 2
implies that
3gR + g I
5gR
2hR = 0,
gR + 3g I
3hR + h I = 0, 5g I
hR
2h I = 0,
3h I = 0
3
1
gR + g I , h I =
2
2
1
3
gR + g I .
2
2
24
2)2 = 0, so the
4u + 4 = (u
g + 3h = 2h ) h = g
c d
hR + jh I
equals
a
c
b
d
gR
hR
b
d
gR + jg I
hR + jh I
a
c
+j
Similarly,
(l R + jl I )
gI
hI
equals
gR
hR
gI
hI
lR
lI
+j
gR
hR
gI
hI
lI
lR
4h = 2g, 24g
10h = 2h ) h = 2g.
=
,
y0
24
10
2x
4x
i.e., points on the same straight line y0 = 2x 0 .
25
26
(106)
Figure 6: Eulers formula ejq =
cos q + j sin q implies that the complex
exponential ejq is a point on the unit circle in the complex plane with coordinates
given by the real and imaginary parts.
=
1
ejq
sin q
q
cos q
<
The product
1
2 cos(p/3) j2 sin(p/3)
e jp/3
=
=
.
4
2
2 ejp/3
2 ejp/3
jp/3
=1
(107)
(108)
d r(t)+jq (t)
dr (t)
dq (t)
e
=
+j
er(t)+jq (t) .
dt
dt
dt
(109)
jwt
= cos wt.
(110)
w sin wt,
(111)
Since
d (a+jb)t
e
= (a + jb) e(a+jb)t
dt
for constants a and b, it follows that
Z t
e(a+jb)t dt =
(112)
e(a+jb)t
.
a + jb
(113)
dt =
ejwt
e
j2w
jwt p/2w
=
0
1
,
w
(114)
which also results from integrating cos wt over the interval [0, p/2w ].
Exercises
1. Write the complex exponential e2+jp/3 in rectangular form.
2. Determine the magnitude of the complex exponential
3. Write the complex number 3
ejp/4 .
j3p/4 .
3 + j.
jpt/2
d jpt/2
e
.
dt
1/2+j)t .
27
Solutions
1. Here, r = e2 and q = p/3. The rectangular form is
p
p
p
1 2
3 2
2
2
e cos + je sin = e + j
e .
3
3
2
2
2. Since
1 = ejp ,
p
ej 4 =
1 ej 4 = ejp ej 4 = ej
5p
4
13 e
j arcsin p2
13
32 + 22 =
13
j=
j 3p
4
3p
= 3 ej 4 .
j arcsin p1
10 e
10
it follows that
3
j p3
e /2
10 e
6. Since
j arcsin p1
10
2e
j p3
p
= 2 10 e
r ejq
r
= e j( q
jj
h
he
j)
2e jp/2
ej2p/7
is 2.
7. Since
8. Here,
3 + j = 2 ej 6 ,
1
e
2
d t
e
dt
jt2
= (1
j p6
j2t) et
jt2
p +arcsin p1
3
10
28
9. Here,
j pt
2
d j pt
e 2 =e
dt
j pt
2
p j pt
e 2 = e
2
10. Since
j pt
2
p
p j p j pt
p
e 2 e 2 = ej 2
2
2
p
j .
2
d
C ejwt = C jw ejwt ,
dt
it follows that
d2
C ejwt = C (jw ) (jw ) ejwt =
dt2
11. The real part of the function e(
shown below.
1/2+j)t
Cw 2 ejwt =
is e
t/2 cos t
w 2 u ( t ).
whose graph is
1
e
t/2 cos t
2p
29
q)
cos q = p
a
a2
+ b2
, sin q = p
b
a2
+ b2
* Physical systems using electronic, mechanical, or fluid components may be modeled as dynamic systems, and designed
to output approximately the integral or derivative of a harmonic input.
30
31
e(t)
controller
u(t)
plant
o (t)
y(t)
sensor
o2 ( t )
i3 ( t )
o (t)
If the gain is denoted by K, then the value of the output o at time t equals
the value of the input i at time t multiplied by K. If the input signal i (t)
received by an amplifier equals e t sin t and the gain equals 5 then the
output signal o (t) returned by the component equals 5i (t) = 5e t sin t.
32
o (t)
c+
Rt
o (t)
If the initial value is denoted by c, then the value of the output o at time t
equals c plus the integral of the input i over the interval [0, t]. If i (t) = te t
and c = 0, then o (t) = 1 (1 + t)e t .
A signal of the form
A cos (wt
(115)
q)
with A > 0 is said to be harmonic with amplitude A, angular frequency w, and phase shift q. A harmonic signal with angular frequency2 w is periodic with period T = 2p/w.
A
A cos(wt
wt0
q
t0
q)
t
(116)
33
is a sum of two harmonic signals with the same angular frequency (and
therefore the same period) and phase shifts 0 and p/2, respectively, since
p
sin wt = cos wt
.
(117)
2
Every harmonic signal A cos (wt
From A cos(wt
that
Square each side, add, and take the square root to obtain the amplitude:
p
A = a2 + b2 .
Let the phase shift q be the unique angle in [ p, p ) that satisfies
the equations
cos q = p
a
a2
+ b2
and sin q = p
b
a2
+ b2
It follows that a signal of the form (116) can be written as a harmonic signal
with period 2p/w and with amplitude and phase shift given in terms of the
individual amplitudes a and b.
Example 20. The derivative of the signal
f (t) := a cos wt + b sin wt
equals
But,
and
d f (t)
=
dt
p
sin wt = cos wt +
= cos w t
2
p
cos wt = sin wt +
= sin w t
2
We conclude that
d f (t)
= w a cos w t
dt
(118)
T
4
(119)
+ b sin w t
T
4
T
4
T
4
(120)
(121)
(122)
T
4
(123)
3
(124)
t
b
cos wt
w
0
b
b
cos wt + .
w
w
a
sin wt
w
a
= sin wt
w
f (t ) dt =
But,
sin wt = cos wt
and
T
4
cos wt = sin wt
We conclude that
Z t
0
1
f (t ) dt = f
w
T
4
(125)
(126)
T
4
(127)
b
.
w
(128)
This equals the original signal f (t) amplified by 1/w, delayed by a quarter
of a period, and added to the constant signal b/w.
Vbc
C
!i
Vac
The voltage drop Vab across the capacitor in the direction of the
current i equals Q/C, where Q is the charge accumulated on the capacitor and C is the capacitance. Similarly, the voltage drop Vbc across
the resistor in the direction of the current equals Ri, where R is the
resistance.
34
(129)
Z t
0
i (t ) dt
(130)
equals the charge accumulated on the capacitor at time t. Differentiation on both sides of this relationship shows that
dQ
( t ) = i ( t ).
dt
(131)
q ),
(132)
dQ(t)
=
dt
wCB sin(wt
(133)
q)
and
wRCB sin(wt
q ).
(134)
q)
wRCB sin(wt
q ).
(135)
wRC cos q )
(136)
and
cos q = p
(137)
1 + (wRC )2
1 + (wRC )2
sin q = p
wRC
1 + (wRC )2
(138)
wRCB sin(wt
q)
wRCA sin wt
(139)
dVac
( t );
(140)
dt
the output voltage across the resistor is approximately proportional
to the derivative of the input voltage across the circuit.
Vbc (t) RC
35
q)
A
sin wt
wRC
1/w,
(141)
1
RC
Z t
0
Vac (t ) dt;
(142)
~Fe
~Fe
i
k
~Fs =
kxi
~Fd =
dvi
Let the unit vector i be parallel to the line of action of the suspension. Denote by x the displacement of the massless plate along i
relative to its position when the spring is undeformed. It follows that
x ( t ) = x (0) +
Z t
0
v(t ) dt,
(143)
36
37
(144)
The force ~Fs from the spring on the massless plate equals kxi,
where k is the stiffness. Similarly, the force ~Fd from the damper on
the massless plate equals dvi, where d is the damping coefficient.
Using Newtons 2nd law5 , and the fact that the plate is massless, we
conclude that
~Fe (t) = ~Fs (t) ~Fd (t)
(145)
equals the external force applied to the plate at time t.
Now suppose that ~Fe (t) = A cos(wt)i and assume that ~Fs (t) =
B cos(wt q )i for some unknown amplitude B and phase shift q.
Since ~Fs = kxi, it follows that
B
cos(wt
k
x (t) =
q ),
(146)
dx (t)
wB
=
sin(wt
dt
k
and
wdB
sin(wt
k
dv(t)i =
~Fd (t) =
(147)
q)
q )i.
(148)
B cos(wt
q) +
d
A = B cos q + w sin q , 0 =
k
or, in other words,
B= p
and
cos q =
1
1 + (wd/k)2
wdB
sin(wt
k
B sin q
q ).
d
w cos q
k
(150)
(151)
1 + (wd/k)2
sin q =
(149)
wd/k
1 + (wd/k )2
(152)
If we design the suspension in Fig. 14, such that d/k 1/w, then
B A and q p. In this case,
~Fd (t) =
wdB
sin(wt
k
wdA
q )i
sin(wt)i
k
(153)
d d~Fe
( t );
(154)
k dt
the force applied to the damper is proportional to the rate of change
of the external force applied to the plate.
If, instead, we design the suspension in Fig. 14, such that d/k
1/w, then wdB/k A and q p/2. In this case,
~Fd (t)
q )i
kA
sin(wt)i
wd
(155)
~Fs (t)
k
d
Z t
0
~Fe (t ) dt;
(156)
pa
g
pa
h
C
Dp
The pressure drop Pv across the valve equals rRq, where R is the
valve resistance, q is the fluid flow rate into the tank, and r is the liquid
38
density. Similarly, the pressure drop Pt between the bottom and top
of the liquid in the tank equals rgh, where g is the acceleration of
gravity and h is the head of the liquid.
By conservation of mass, it follows that
rq = rC
dh
,
dt
(157)
(158)
p0
B
+
cos(wt
rg rg
q ),
(159)
dh(t)
=
dt
and
Pv (t) = rRq(t) =
wCB
sin(wt
g
q)
(160)
wRCB
sin(wt
g
q ).
(161)
q)
wRCB
sin(wt
g
RC
A = B cos q + w
sin q , 0 = B sin q
g
or, in other words,
B= p
and
cos q = p
1
1 + (wRC/g)2
RC
w
cos q
g
q ).
(162)
(163)
(164)
1 + (wRC/g)2
sin q = p
wRC/g
1 + (wRC/g)2
(165)
If we design the reservor in Fig. 15, such that RC/g 1/w, then
B A and q 0. In this case,
Pv (t) =
wRCB
sin(wt
g
q)
wRCA
sin wt
g
RC dDp
( t );
g dt
(166)
(167)
39
p0 = B cos(wt
q)
gA
sin wt
wRC
(168)
p0 ) dt;
(169)
g
RC
p0
Z t
0
(Dp(t )
The time constant of the water reservoir equals the ratio RC/g
of the product of the valve resistance and tank cross section
with the acceleration of gravity. When the time constant is
much smaller than the period of the harmonic input signal
(RC/g 1/w), the reservoir acts as a differentiator. When
the time constant is much greater than the period of the harmonic input signal (RC/g
1/w), the reservoir acts as an
integrator.
Exercises
1. Find a unique angle q 2 [ p, p ), such that
1
cos q = p , sin q =
2
1
p .
2
1
p , sin q =
2
1
p .
2
sin 2t = A cos(wt
q ).
40
p/4) + 4 sin(t
2B sin(t
q)
q ).
8. Graph the signal f (t) := 3 cos pt sin pt, its derivative, and its
definite integral over the interval [0, t] in the same diagram.
9. Suppose that the input signal to each of the components below
is given by i (t) = e t for t
0 and 0 otherwise. In each case,
determine the output signal o (t) for t 0, and graph i (t) and o (t).
i(t)
i(t)
i(t)
o (t)
0.6
o (t)
0.6
0.6 +
Rt
o (t)
Solutions
p
1. Since cos is an even function, cos q = 1/ 2 implies that q is either
p/4 or p/4. Only the latter satisfies the condition that sin q =
p
1/ 2.
p
2. Since cos is an even function, cos q = 1/ 2 implies that q is
either 3p/4 or 3p/4. Only the latter satisfies the condition that
p
sin q = 1/ 2.
3. Since tan is a periodic function with period p, tan q = 1 if q is
either p/4 or 3p/4.
4. Since A cos(wt
w = 2 and
1.
2.
41
3/2,
p
3
3 3
3 sin(4t p/3) = sin 4t
cos 4t.
2
2
p
It follows that w = 4, a = 3 3/2, and b = 3/2.
6. Since
3 cos(t
4 sin(t
3
3
p/4) = p cos t + p sin t,
2
2
4
4
p/4) = p sin t p cos t,
2
2
and
B sin(wt + j) = B sin wt cos j + B cos wt sin j,
it follows that w = 1 and
1
7
p = B sin j, p = B cos j.
2
2
Squaring and adding gives 25 = B2 ) B = 5. The angle j lies in
p
the fourth quadrant and equals arcsin(1/5 2).
7. Since
cos t = cos(t
q + q ) = cos(t
q ) cos q
sin(t
q ) sin q,
it follows that
cos q = B, sin q = 2B,
p
i.e., B = 1/ 5 and q = arcsin(2/ 5).
8. Here,
and
d f (t)
=
dt
Z t
0
f (t ) dt =
3p sin pt
p cos pt
3
1
sin pt + cos pt
p
p
1
.
p
d f (t)/dt
5
p
Rt
0
for t
f (t ) dt
0 and 0 otherwise.
42
1
i(t)
o (t)
1
(t 0.6)
for t
1
i(t)
o (t)
1
For the integrator, o (t) = 1.6
for t
0 and 0 otherwise.
o (t)
1
i(t)
43
44
45
i(t)
i1 ( t ) +
i(t)
i1 ( t )
i2 ( t )
i3 ( t )
i(t)
i(t)
i3 ( t )
x (t)
c+
Here,
x (t) = Ki (t),
o (t) = c +
Z t
0
Rt
o (t)
Z t
0
x (t ) dt.
i (t ) dt.
Z t
o (t) = K c +
i (t ) dt ,
0
(170)
(171)
(172)
Here,
x (t) = c +
x (t)
Rt
c+
o (t)
o (t) = x (t
D ).
Z t
0
i (t ) dt,
(173)
Z t D
0
i (t ) dt.
(174)
Z t
0
D) dt = c +
i(t
Z t D
D
i (s) ds,
(175)
Z t
0
i (t ) dt
(176)
o (0) = c,
(177)
x (t)
Rt
c+
o (t)
Here,
x (t) = i(t)
ao (t),
o (t) = c +
Z t
0
x (t ) dt.
(178)
Z t
0
i(t )
ao (t ) dt,
(179)
46
o (0) = c,
(180)
the general form of a first-order, scalar, constant-coefficient, linear, nonhomogeneous differential equation.
Example 25. Consider the following block diagram
Figure 20: A double integrator with
feedback.
i(t) +
x (t)
y(t)
ay(t)
bo (t) and
c+
Rt
d+
Rt
0
o (t)
Z t
0
x (t ) dt,
o (t) = d +
Z t
0
y(t ) dt.
(181)
ay
bo,
do
=y
dt
(182)
o (0) = d,
do
(0) = c;
dt
(183)
the general form of a second-order, scalar, constant-coefficient, linear, nonhomogeneous differential equation.
47
Free-body diagram:
f (t)i
f (t)i
i
d
~Fs =
kxi
~Fd =
dvi
1
d
x (0) +
Rt
0
x (t)
48
49
x2
k2
x1
d2
d1
k1
Let the unit vector i again be parallel to the line of action of the
suspension. Denote by x1 and x2 the displacements of the lower and
upper massless plates, respectively, along i relative to their position
in the configuration where both springs are undeformed. It follows
that
x 1 (t) :=
dx1
( t ) = v1 ( t ),
dt
x 2 (t) :=
dx2
( t ) = v2 ( t ),
dt
(187)
where v1 (t) and v2 (t) are the velocities of the two plates along i at
time t.
Consider the free-body diagrams shown on the next page. If we
apply Newtons 2nd law to the lower plate, we obtain the equation of
force balance
x1 (t)) = k1 x1 (t) + d1 x 1 (t),
k 2 ( x2 ( t )
(188)
x1 (t) + d2 x 2 (t),
(189)
Z t
k 2 x2 ( t )
0
x1 ( t )
d1
k 1 x1 ( t )
dt
(190)
f (t)i
Free-body diagrams:
k 2 ( x2
x1 )i
x1 )i
k 2 ( x2
k1 x1 i
d2 v2 i
d1 v1 i
and
x2 ( t ) = x2 (0) +
Z t
f (t )
0
k 2 x2 ( t )
d2
x1 ( t )
dt
(191)
1
d2
x2 (0) +
Rt
0
x2 ( t )
+
k2
1
d1
x1 (0) +
Rt
0
x1 ( t )
k1
Now suppose that the upper plate has a nonzero mass m. Application of Newtons 2nd law then results in
m x2 (t) = f (t)
k 2 x2 ( t )
where
x2 (t) :=
x1 ( t )
d2 x2
( t ).
dt2
50
d2 x 2 (t)
(192)
(193)
If we let y(t) := x 2 (t) and y (t) = x2 (t), then this is equivalent to the
first-order equations
x 2 (t) = y(t),
my (t) = f (t)
k 2 x2 ( t )
x1 ( t )
(194)
d2 y ( t )
Z t
0
y(t ) dt
(195)
and
y ( t ) = y (0) +
Z t
f (t )
k 2 x2 ( t )
x1 ( t )
d2 y ( t )
dt
(196)
d2
f (t) +
1
m
y (0) +
Rt
0
y(t)
x2 (0) +
Rt
0
x2 ( t )
+
k2
1
d1
x1 (0) +
Rt
x1 ( t )
k1
Exercises
1. Determine i (t) and o (t) such that
i(t)
2 sin 2t
1+
Rt
0
o (t)
51
i(t)
2+
2 cos 3t
Rt
o (t)
c+
x (t)
Rt
0
4e
t sin 2t
x (t)
Rt
o (t)
x (t)
o (t)
Rt
0
do
( t ) + o ( t ) = i ( t ).
dt
do
(t) = i(t)
dt
o ( t ).
1,
do
d2 o
do
(0) = 2,
(t) + (t) + 4o (t) = i (t).
2
dt
dt
dt
Rt
0
y(t) +
1+
Rt
0
o (t)
52
1
3t + e
9
3t
1 , x2 ( t ) =
1
6t
18
3t
+1
Rt
1
2
x2 ( t )
+
2
Rt
x1 ( t )
Solutions
1. Here,
2 sin 2t = 2i (t), o (t) =
1+
Z t
0
2 sin 2t dt.
cos 2t.
2. Here,
2 cos 3t =
2+
Z t
0
Z t
0
i (t ) dt, 4e
4 cos 3t.
sin 2t = Kx (t)
x (t) =
1,
x (t) =
Z t
0
Z t
0
dt = t, o (t) = 0.
dt = t, o (t) = x (t
1) = t
1.
53
5. For 0 t 1,
x (t) = 0, o (t) =
For t
1,
x (t) = 1, o (t) =
Z t
0
Z t
0
x (t ) dt = 0.
x (t ) dt =
Z t
1
dt = t
1.
6. Here,
do
(t) = i(t)
dt
o ( t ), o (0) = 3 ) o ( t ) = 3 +
Z t
o (t ) dt
i(t )
3+
Rt
o (t)
7. Here,
do
i(t) o (t)
(t) =
, o (0) = 0 ) o ( t )
dt
2
Z t
i (t )
o (t )
2
dt.
Rt
1
2
o (t)
4o (t)
1+
Z t
0
y(t ) dt
y(t) ) y(t) = 2 +
Z t
0
i(t )
4o (t )
y(t ) dt
Rt
i(t) +
2+
Rt
0
y(t)
1+
o (t)
54
9. Here,
dy
(t) + 3y(t) = i (t), y(0) = 0
dt
do
(t) + 2o (t) = y(t), o (0) = 1
dt
This implies that
d2 o
do
dy
(t) + 2 (t) =
(t) = i(t)
dt
dt
dt2
and
3y(t) = i (t)
do
(0) = y (0)
dt
do
(t)
dt
6o (t)
2o (0) = 2
x1 ( t ) =
1
6
1
e
6
3t
Consequently,
Z t
0
2 x2 ( t )
x1 (t ) dt =
Z t
1
0
1
e
3
3t
dt = x1 (t)
and
Z t
1
0
2 x2 ( t )
x1 ( t )
dt =
Z t
1
0
1
+ e
3 6
3t
dt = x2 (t)
55
at
f (#) g (#) ( t ) : =
y (0) = c
f (#) e
Z t
0
a#
f (t ) g(t
( t ),
t ) dt
d (#) e
a#
=e
at
56
57
y (0) = c
(197)
where we have written y for the output and f for the input. The
solution to this initial-value problem is given by
at
y(t) = ce
If we let6
a#
f (#) e
Z t
0
(t) :=
f (t )e
a(t t )
dt.
(198)
6
Z t
f (t )e
a(t t )
dt
(199)
f (#) e
a#
(0) =
Z 0
f (t )eat dt = 0,
Z t
0
g(t, t) dt =
It follows that
d
f (#) e
dt
a#
Z t
g(t, t)
(t) =
dt + g(t, t).
a f (#) e
a#
(t) + f (t)
Finally,
dy
(t) =
dt
ace
at
a f (#) e
a#
(t) + f (t) =
ay(t) + f (t).
Z t
0
f (t ) dt.
(201)
at
(202)
Figure 27: The free response.
y(t)
c
ce
at
c/e
1/a
The output decays to a fraction 1/e of its original value after the characteristic time 1/a, known as the time constant of the system.
The unit step u(t) is a signal that equals 1 for t
wise.
0 and 0 other-
u(t)
Z t
t
1
y(t) = u(#) e a# (t) :=
e a(t t ) dt = e a(t t )
(203)
a
0
0
or, in other words,
1
1 at
e .
(204)
a a
The output converges to a steady-state value yss = 1/a in such a way
that y(t) yss decays to 0 exactly as the free response of a first-order system
with time constant 1/a.
y(t) =
58
u(t) u(t D)
D
1/D
1
D
Z D
0
a(t t )
u(t
dt =
at Z D
eat dt =
e at at
e
aD
(206)
0
eaD 1
e
aD
Using lHospitals rule, we recall that
y(t) =
lim
x !0
ex
1
x
at
= 1.
(207)
(208)
In the limit as the pulse width D ! 0, the output thus converges to the free
response with initial condition equal to 1. We call this the unit impulse
response. In particular, we have
u (#) u (# D )
at
a#
e
= lim
e
( t ).
(209)
D
D !0
We represent the right-hand side symbolically by the convolution
d(#) e a# (t).
(210)
59
60
The unit impulse d(t) acts in convolution with the exponentially decaying signal e at as multiplication by 1.
d(t)
D!0
)
t
i
d
f (t) +
1
d
x (0) +
Rt
0
x (t)
Z t
f (#)
f (t )
x (t) =
e k#/d (t) =
e
d
d
0
k (t t )/d
dt.
(211)
D
Then, for 0 t D,
x (t) =
Z t
1
0
dD
k(t t )/d
dt =
1
1
kD
1
e
kD
kt/d
k (t t )/d
t =t
(212)
t =0
(213)
D,
x (t) = x (D)e
k (t D)/d
(214)
But,
1
1 e kD/d ,
(215)
kD
which converges to 1/d as D ! 0. By substitution into (214) It follows that the unit impulse response equals
x (D) =
x (t) =
Suppose that f (t) = 2(D
1
e
d
kt/d
(216)
f (t)
2/D
61
Then, for 0 t D,
x (t) =
Z t
2( D
0
t)
dD2
k (t t )/d
dt = 2d
1 + k (D t )/d
e
k 2 D2
k (t t )/d
t =t
t =0
(217)
1 + k(D t)/d
k 2 D2
2d
1 + kD/d
e
k 2 D2
kt/d
(218)
D,
x (t) = x (D)e
k (t D)/d
(219)
But,
2d
kD/d
1
1
+
kD/d
e
(220)
k 2 D2
which also converges to 1/d as D ! 0 by lHospitals rule. It follows that the response to the given input f (t) converges to the unit
impulse response as D ! 0.
x (D) =
(a f (t ) + g(t )) dt = a
Z t
0
f (t ) dt +
Z t
0
g(t ) dt.
(221)
D ) + i2 ( t
D)
(222)
A system consisting of a single delay line is therefore linear. It is straightforward to show that the same conclusion applies to a system consisting of a
single amplifier.
Block diagrams consisting of linear components, summing, and
splitting junctions are linear.
Example 30. For a zero-initial-value integrator with feedback, linearity
follows from a study of the two block diagrams below. It also follows from
the explicit solution
Z t
a#
o ( t ) = i (#) e
(t) :=
i (t )e a(t t ) dt
(223)
0
62
i1 ( t ), i2 ( t ) +
Rt
0
o1 ( t ), o2 ( t )
a
ai1 (t) + i2 (t) +
Rt
Z t
y(t) = f (#)u(#) e a# (t) :=
f (t )u(t )e a(t t ) dt.
(224)
0
Z t
D
f (t
D)u(t
D)e
a(t t )
dt
(225)
Z t D
0
f (s)u(s)e
a(t D s)
f (#)u(#) e a# (t D),
ds,
63
(226)
(227)
i.e., the output signal obtained with f (t)u(t) as input, but delayed by D.
Thus, since
(228)
Z t
i (#) e a# (t) =
i (t )e a(t t ) dt,
(229)
0
Exercises
1. Use a variable substitution to show that
u (#) F (#) ( t ) = f ( t )
3. Compute the convolution
a#
b#
( t ).
when a 6= b.
4. Use the previous result to compute the convolution
64
for t
a#
sin w# e
(t)
1/a.
dy
(t) + 2y(t) = t,
dt
7. Find the solution for t
y (0) =
1.
dy
(t) + y(t) = e t ,
dt
y(0) = 1.
8. Find the free response, the unit step response, and the unit impulse response of the dynamic system shown below.
i(t) +
2+
Rt
o (t)
1
2
9. Find the free response, the unit step response, and the unit impulse response of the dynamic system shown below.
i(t) +
1+
Rt
o (t)
1
2
10. Find the response of the dynamic system below in the limit as
b ! and compare this to the unit impulse response. Graph the
input and output for several representative values of b.
be
bt
Rt
o (t)
65
Solutions
1. Here,
f (#) g (#) ( t ) =
f (t
2. Here,
s ) g(s ) ds =
u (#) F (#) ( t ) =
t ) dt
f (t ) g(t
t,
Z t
g (#) f (#) ( t )
t ) dt
F (t )u(t
t ) = 1 for 0 t t,
Z t
F (t ) dt = f (t)
3. Here,
Z t
a#
b#
(t) =
f (0) = f ( t ).
Z t
eat e b(t
t =t
e(a b)t + bt
a b
=
t =0
eat
a
t)
dt
e bt
.
b
since a 6= b.
4. Since
sin wt =
ejwt
it follows that
j2#
e
sin 2# sin 3# (t) =
e
j2
e
j2
jwt
j2#
ej3#
e
j2
But,
ej3t
ej2# ej3# (t) =
ej2# e
j3#
(t) =
ej2t
j3#
( t ).
ej2t
j
e
j5
j3t
ej3t e j2t
ej3# (t) =
j5
j2t
e
e j3t
e j2# e j3# (t) =
j
e
j2#
1 ej3t ej2t
ej2t e j3t
ej3t e
4
j
j5
j5
2
3
=
sin 3t + sin 2t
5
5
j2t
j2t
e
j
j3t
66
5. Since
ejwt
sin wt =
it follows that
sin w# e
a#
1 jw#
e e
j2
(t) =
jwt
e
j2
a#
1
e
j2
(t)
jw#
a#
But,
For t
ejw# e
jw#
a#
a#
(t) =
(t) =
ejwt e at
a + jw
jwt
e
jw
1/a, e at 0 and
jwt
1
e
a#
sin w# e
(t)
j2 a + jw
at
jwt
e
a
jw
a sin wt w cos wt
a2 + w 2
#e
2#
(t) =
Z t
0
#e
te
2#
2( t t )
it follows that
8. For t
y(t) =
t
2
Z t
(t) =
1
4
t
2
t
2
1 1
+ e
4 4
1
4
3
e
4
(t) =
It follows that
7. Since
dt =
2t
2t
(t t )
dt = te t ,
y ( t ) = (1 + t ) e t .
0,
o (t) = 2e
t/2
+ i (#) e
2( t t )
#/2
( t ).
t =t
t =0
(t)
67
t/2
For the unit step response, i (t) = u(t) and we ignore the contribution from the initial condition, in which case
#/2
o ( t ) = u (#) e
(t) = 2(1 e t/2 ).
For the unit impulse response, i (t) = d(t) and we ignore the
contribution from the initial condition, in which case
#/2
o ( t ) = d (#) e
(t) = e t/2 .
9. Here,
Z t
1+
o (t) =
2 i(t )
1
o (t )
2
do
(t) + o (t) = 2i (t), o (0) =
dt
It follows that, for t
dt
0,
o (t) =
+ 2i (#) e
(t)
e t.
For the unit step response, i (t) = u(t) and we ignore the contribution from the initial condition, in which case
#
o (t) = 2u(#) e
( t ) = 2(1 e t ).
For the unit impulse response, i (t) = d(t) and we ignore the
contribution from the initial condition, in which case
#
o (t) = 2d(#) e
(t) = 2e t .
10. For t
0 and b 6= 2,
o (t) = be
b#
2#
(t) = b
bt
e
b
2t
In the limit as b ! , o (t) converges to e 2t , i.e., the unit impulse response. The diagrams below show the input and output,
respectively, for three different values of b.
68
b!
be
bt
2t
o (t) = b e
bt
e
2 b
2t
69
L f (#) ( s ) : =
st
f (t)e
dt
1
s
L f (#
Moreover,
L
and
D) (s) = e
D ) u (#
Z
#
0
f (t ) dt (s) =
sD
L f (#) u (#) ( s ).
1
L f (#) ( s )
s
L o (#) ( s ) = L i (#) ( s ) H ( s ),
o ( t ) = i (#) h (#) ( t ).
70
L f (#) ( s ) : =
f (t)e st dt,
(230)
0
then the action of amplifiers, delay lines, zero-initial-value integrators, and zero-initial-value integrators with feedback in the time
domain correspond to multiplication with an appropriate function of
the complex number s in the frequency domain.
Example 32. The zero-initial-value integrator with negative feedback
may be represented in a block diagram in the time domain or the frequency
domain, as shown below.
Time domain:
i(t)
Rt
o (t)
a
Frequency domain:
L i (#) ( s ) +
1
s
L o (#) ( s )
Z t
0
i(t )
ao (t ) dt.
1
L o (#) ( s ) =
L i (#) ( s )
s
L o (#) ( s ) =
(231)
a L o (#) ( s )
1
L i (#) ( s )
s+a
(232)
(233)
We obtain the output signal o (t), for a given input i (t), by returning to the
time domain on both sides of this equation. It should come as no surprise
that, in this case, the answer is in the form of a convolution.
71
all t
0 and some M > 0, the Laplace transform L f (#) (s) is finite
for all complex numbers s with real part greater than a. In particular,
(234)
Example 33. The unit step f (t) = u(t) satisfies the inequality | f (t)|
Me at for all t
0 for M = 1 and a = 0. Its Laplace transform is finite
provided that the real part of s is positive and equals
Z
L u (#) ( s ) =
st
dt =
1
e
s
st
1
= .
s
(235)
The unit step is the unique signal whose Laplace transform equals 1/s.
Example 34. For every real or complex number a, the exponential function
f (t) = eat satisfies the inequality | f (t)| Me at for all t
0 for M = 1
and a equal to the real part of a. Its Laplace transform is finite provided that
the real part of s is greater than the real part of a, in which case
Z
h i
L ea# (s) =
e(a
s)t
dt =
1
a
e(a
s)t
=
0
1
s
(236)
The exponential function eat is the unique signal whose Laplace transform
equals 1/(s a).
If the input signal to the zero-initial-value integrator with negative
feedback in Example 32 is the unit step, i.e., i (t) = u(t), then
1
1
1 1
1
L o (#) ( s ) =
L i (#) ( s ) =
=
(237)
s+a
s(s + a)
a s
s+a
or, by linearity,
1
1 e at
(238)
a
0. This analysis is an example of the general observation:
o (t) =
for t
1
L u (#) ( s ) = .
s
72
Z
u (#) u (# D )
1 D
L
(s) =
e
D
D 0
By lHospitals rule,
1
lim
D !0
e
Ds
Ds
st
dt =
e
Ds
Ds
(239)
=1
(240)
Z t
0
73
f (t ) dt
0 and
(241)
F (t)e
st
1
F (t)e
s
dt =
st
+
0
1
s
Z
0
F0 (t)e
st
dt
(242)
#
0
f (t ) dt (s) =
1
L f (#) ( s ).
s
(243)
Time domain:
i(t)
o (t)
Rt
0
L i (#) ( s )
1
s
L o (#) ( s )
f (t ) dt =
Z t
0
u(t
t ) f (t ) dt =
u (#) f (#) ( t ).
(244)
(245)
i.e., convolution with the unit step in the time domain is equivalent to multiplication by the Laplace transform of the unit step in the frequency domain.
Z
0
st
Z Z t
0
t
0
st
t ) g(t ) dt
f (t
dt
t ) g(t ) dt dt
f (t
(246)
t=t
t=t
dt
t
dt
=
=
Z Z
0
Z
0
st0
s(t +t0 )
f (t0 ) dt0
st
g(t ) dt.
(247)
L o (#) ( s ) = L i (#) ( s )
1
.
( a)
(248)
74
Time domain:
i(t)
Rt
o (t)
75
a
Frequency domain:
L i (#) ( s )
L o (#) ( s )
1
s+a
L o (#) ( s ) = L i (#) ( s ) H ( s )
(249)
implies that
o (t) =
i (#) h (#) ( t )
(250)
o (t)
h(t)
or
L i (#) ( s )
H (s)
L o (#) ( s )
t ) = h(t) +
dh
(t)t + g(t
dt
t )t
(251)
u (#) u (# D )
h (#) ( t )
D
Z
1 D
dh
=
h(t) +
(t)t + g(t t )t dt
D 0
dt
= h(t) +
dh D
1
(t) +
dt
2
D
Z D
0
g(t
t )t dt
(252)
Figure 37: A convolution system represented in the time domain by convolution with the unit impulse response
h(t) or in the frequency domain by
multiplication with the transfer function
H ( s ).
D !0
u (#)
u (#
D
D)
h (#) ( t ) = h ( t )
(253)
It follows that h(t) is the unit impulse response of the system, and
H (s) is its Laplace transform.
Example 38. The transfer function of a zero-initial-value integrator is
1/s. But L u(#) (s) = 1/s. It follows that the unit impulse response of a
zero-initial-value integrator is the unit step u(t).
Similarly, the transfer function of a zero-initial-value integrator with
feedback is 1/(s + a). But L e a# (s) = 1/(s + a). It follows that the
unit impulse response of a zero-initial-value integrator with feedback is the
exponential function e at .
Example 39. Now consider the following equivalent block diagrams with
I (s) := L i (#) (s), X (s) = L x (#) (s), and O(s) := L o (#) (s).
b
i(t) +
x (t)
Rt
0
Rt
o (t)
I (s) +
X (s)
1
s
1
s
O(s)
It follows that
O(s) =
1
X ( s ),
s
X (s) =
1
I (s)
s
aX (s)
bO(s)
(254)
76
77
s2
(255)
1
I ( s ).
+ as + b
(256)
As suggested in the figure below, the system is thus equivalent to a convolution with transfer function
H (s) =
1
s2 + as + b
I (s)
(257)
O(s)
s2 +as+ b
I (s) +
1
s
O(s)
Let I (s) := L i (#) (s) and O(s) := L o (#) (s). It follows that
O(s) = bI (s) +
1
I (s)
s
aO(s)
(258)
b + 1s
bs + 1
I (s) =
I (s)
1 + as
s+a
(259)
As suggested in the figure below, the system is thus equivalent to a convolution with transfer function
H (s) =
bs + 1
1 ab
= b+
.
s+a
s+a
(260)
I (s)
Since
bs+1
s+a
H (s) = L bd(#) + (1
78
O(s)
a#
ab)e
( s ),
(261)
x2
k2
x1
d2
d1
k1
X1 ( s ) =
k 2 X2 ( s )
and
X1 ( s )
d1 s
k 1 X1 ( s )
k 2 X2 ( s ) X1 ( s )
d2 s
After solving for X1 (s) and X2 (s), we obtain
X2 ( s ) =
X1 ( s ) =
d1 d2
s2
F (s)
k2 F (s)
= H1 (s) F (s)
+ ( d2 k 1 + d1 k 2 + d2 k 2 ) s + k 1 k 2
(262)
(263)
(264)
79
and
X2 ( s ) =
( d2 s + k 1 + k 2 ) F ( s )
= H2 (s) F (s)
d1 d2 s2 + ( d2 k 1 + d1 k 2 + d2 k 2 ) s + k 1 k 2
(265)
f (t) +
1
d2
x2 (0) +
Rt
0
x2 ( t )
+
k2
1
d1
x1 (0) +
Rt
0
x1 ( t )
k1
F (s) +
1
d2
1
s
X2 ( s )
+
k2
1
d1
1
s
k1
X1 ( s )
80
F (s)
H1 (s)
H2 (s)
X1 ( s )
X2 ( s )
qi
pa
qo
pa
g
C
Ro
Assume that the mass flow rate out of the tank through a valve is
proportional to the pressure drop across the valve:
rqo =
Dp
,
Ro
(267)
81
(269)
CRo
.
g
(270)
If y(0) = y0 , then
y ( t ) = y0 e
t/T
1
q (#) e
C i
#/T
( t ).
(271)
If we consider qi (t) as the input and y(t) as the output, then for
y(0) = 0, the system is a convolution with transfer function
H (s) =
1 T
.
C sT + 1
(272)
L q i (#) ( s )
L y (#) ( s )
1 T
C sT +1
Example 41. The unit impulse response of the single-tank liquid level
system considered above is the signal whose Laplace transform equals the
transfer function, namely
y(t) =
1
e
C
t/T
(273)
which also equals the free response (in the absence of input) given an initial
condition y(0) = 1/C.
Example 42. The unit step response of the single-tank liquid level system
considered above is the signal whose Laplace transform equals
11 T
T 1
T
=
(274)
s C sT + 1
C s
sT + 1
or, in other words,
T
1 e t/T
(275)
C
for t 0. In particular, y (t) ! T/C = Ro /g as t ! corresponding to
a steady-state level of the liquid in the tank, as the amount of liquid flowing
into the tank equals that flowing out of the tank.
y(t) =
We replace the open inlet with a pump that raises the ambient
pressure by dp(t), and an inlet valve with resistance Ri connected to
the bottom of the tank, as shown in the figure below.
pa
g
pa
qi
dp
Ri
In this case,
Ro
d
(rCy) = rqi
dt
where
rqi =
dp
and
rqo =
so that
82
rqo ,
qo
(276)
rgy
Ri
(277)
rgy
,
Ro
(278)
dy
1
1
(t) + y(t) =
dp(t),
dt
T
rCRi
(279)
CRi Ro
.
g ( Ri + R o )
(280)
If y(0) = h0 , then
y ( t ) = y0 e
t/T
1
+
dp(#) e
rCRi
#/T
( t ).
(281)
If we consider dp(t) as the input and y(t) as the output, then for
y(0) = 0, the system is a convolution with transfer function
H (s) =
1
T
.
rCRi sT + 1
(282)
L dp(#) (s)
1
T
rCRi sT +1
L y (#) ( s )
pa
g
pa
qi
y1
dp
C1
R1
R2
C2
pa
pa
qo
y2
R3
across the pump and the liquid height in the first tank is a convolution with
transfer function
1
T1
H1 (s) =
,
(283)
rC1 R1 sT1 + 1
where
T1 =
C1 R1 R2
g ( R1 + R2 )
(284)
1
T2
,
C2 sT2 + 1
C2 R3
.
g
(285)
(286)
Finally, Eq. (278) implies that the relationship between the liquid height in
the first tank and the fluid flow rate into the second tank is a convolution
with transfer function
g
H1 2 (s) =
.
(287)
R2
The coupled liquid-level system may now be represented by the block
diagram below.
83
L dp(#) (s)
L y1 (#) ( s ) L q i (#) ( s )
T1
1
rC1 R1 sT1 +1
g
R2
T2
1
C2 sT2 +1
L y2 (#) ( s )
(288)
g
T1 T2
rC1 R1 C2 R2 (sT1 + 1)(sT2 + 1)
(289)
H (s) =
or, in other words,
H (s) =
Exercises
1. Consider the following block diagram in the time domain.
i(t) +
Rt
1
2
o (t)
Rt
0
y(t) +
Rt
0
o (t)
84
Rt
1
2
o (t)
+
2
Rt
0
y(t)
ea#
a
e b#
(s)
b
n!
L #n ( s ) = n +1 ,
s
dn
L #n f (# ) ( s ) = ( 1 ) n n L f (# ) ( s )
ds
and use the result to show that
n!
L #n ( s ) = n +1 ,
s
85
9. Show that
L f (#
D) (s) = e
D ) u (#
sD
L f (#) u (#) ( s ).
and illustrate the result using block diagrams in the time and
frequency domains.
10. Suppose that
L o (#) ( s ) =
1
L i (#) ( s )
s ( s + 2)
s+1
L o (#) ( s ) = 2
L i (#) ( s )
s +1
Draw a corresponding frequency-domain block diagram.
1 + sb
L o (#) ( s ) =
L i (#) ( s )
s+a
Solutions
1. A corresponding frequency-domain block diagram is given below.
L i (#) ( s ) +
Here,
1
2
1
s
1
L o (#) ( s ) =
L i (#) ( s )
2s
L o (#) ( s ) =
L o (#) ( s )
1
L i (#) ( s ).
1 + 2s
L o (#) ( s )
86
L i (#) ( s ) +
1
s
L y (#) ( s )
1
s
L o (#) ( s )
3
Here,
1
L o (#) ( s ) =
L y (#) ( s )
s
1
L y (#) ( s ) =
L i (#) ( s )
s
L o (#) ( s ) =
2L o (#) ( s )
3L y (#) ( s )
1
L i (#) ( s )
(s + 2)(s + 3)
L i (#) ( s ) +
L o (#) ( s )
1
1
s
2
1
s
L y (#) ( s )
Here,
1
L o (#) ( s ) =
L i (#) ( s ) 2 L o (#) ( s )
2s
2
L y (#) ( s ) =
L o (#) ( s ) L y (#) ( s )
s
L o (#) ( s ) =
s+2
L i (#) ( s )
2s(s + 3)
L y (#) ( s )
87
ea#
a
a ),
e b#
(s) =
b
L sin w# (s) =
1
s jw
L cos w# (s) =
1
s+jw
1
a)(s
(s
j2
Since,
1
s b
b)
jwt ) /j2,
1
s a
s2
w
+ w2
jwt ) /2,
1
s jw
1
s+jw
s
s2 + w 2
w
L sin w# (s) = L cos w# (s)
s
this agrees with (243).
6. Recall that L 1 (s) = 1/s, provided that the real part of s is greater
than 0. With f (t) = ntn 1 , it follows from (243) that
i
n h
L #n ( s ) = L #n 1 ( s )
s
and, consequently, that
7. Here,
nn 1 h n
L #n ( s ) =
L #
s s
( 1) n
(s) = =
dn
n d
L
f
(
#
)
(
s
)
=
(
1
)
dsn
dsn
= ( 1) n
Z t
0
n!
s n +1
f (t)e
st
dt
( t)n f (t)e
st
dt
or in other words,
dn
L #n f (# ) ( s ) = ( 1 ) n n L f (# ) ( s ).
ds
dn 1
n!
L #n ( s ) = ( 1 ) n n = n +1 .
ds s
s
h i
L #e# (s) =
d 1
1
=
.
ds s 1
( s 1)2
88
L f (#
D) (s) =
D ) u (#
sD
=e
D)u(t
f (t
Z
0
st
D)e
st
f (t )u(t )e
dt
dt,
Time domain:
i(t)
10. Since
1
1
=
s ( s + 2)
2
it follows that
L i (#) ( s )
o (t)
1
s
1
s+2
1
= L 1
2
L o (#) ( s )
Ds
2#
( s ),
1
2#
o (t) =
1 e
i (#) ( t ).
2
11. Let O(s) = L o (#) (s) and I (s) = L i (#) (s) so that
1
I (s)
s2
O(s) +
1
I ( s ).
s
Then,
Y (s) =
1
I (s)
s
O(s) ) O(s) =
1
Y ( s ) + O ( s ) + Y ( s ).
s
I (s) +
1
s
1
s
Y (s)
12. Since
h
1 + sb
1 ab
= b+
= L bd(#) + (1
s+a
s+a
it follows that
o (t) = bi (t) + (1
ab) e
a#
ab)e
a#
i (#) ( t ).
( s ),
O(s)
89
13. Since
2
s
1
3
= L 2e3#
e2# ](s),
e2t .
1
1)2
(s
1)2
90
1
=
(s
s2 + as + b
1
g1 )(s
g2 )
D at
e
= ae at + d(t).
dt
If limt#0 h(t) = 0, then sH (s) is the Laplace transform of the
derivative of h(t).
* If the function H (s) is a ratio of two polynomials in s, then
the corresponding inverse Laplace transform can be computed
using convolution integrals and distributional derivatives.
91
92
1
1
1
=
s g1 s g2
s2 + as + b
(290)
(291)
where
g1,2 :=
a2
4
Z t
Z t
g1 #
g2 #
e e
(t) =
eg1 t eg2 (t t ) dt = eg2 t
e(g1 g2 )t dt
(294)
0
eg1 t
g1
eg2 t
g2
(295)
(296)
if g1 = g2 = g.
Example 44. A double zero-initial-value integrator with negative feedback
is overdamped if a2 > 4b. In this case, g1 and g2 are real, negative, and
distinct.
h(t) = e
( a2 /4 b t)
at/2 sinhp
a2 /4 b
The unit impulse response shows an initial increase from zero and a
subsequent decrease back to zero without any overshoot.
93
h(t) = e
jx ) /j2.
2
at/2 sin(p b a /4 t)
2
b a /4
a2 /4
p
time constant 2/a, angular frequency b a2 /4, and no phase shift.
b
(298)
h(t) = te
at/2
(299)
at/2
94
4/a2
h (#) u (#) ( t )
t
Z t
h (#) u (#) ( t ) =
te at/2 u(t
0
=
or, in other words,
h (#) u (#) ( t ) =
2e
at/2 ( at
+ 2)
a2
4
2e
t ) dt
t
(300)
0
at/2 ( at + 2)
(301)
a2
at
= d(t)
ae
at
(303)
(304)
d(t) + jw e
j2
jwt
=w
ejwt + e
2
jwt
(305)
0,
D
sin wt = w cos wt
dt
(306)
i.e., the regular derivative of sin wt for all t. Since L sin w# (s) =
w/(s2 + w 2 ), it follows that the Laplace transform of cos wt equals
s
w
s
= 2
w s2 + w 2
s + w2
(307)
In contrast,
d(t) + jw ejwt + d(t)
2
jw e
jwt
= d(t)
ejwt
e
j2
jwt
(308)
0,
D
cos wt = d(t) w sin wt,
(309)
dt
i.e., the regular derivative of cos wt plus a unit impulse at t = 0 corresponding to the discontinuous jump from 0 to cos(0) = 1.
Example 48. The Laplace transform of the convolution
a#
a#
d(#) ae
e
(t) = e at ate
equals
at
(310)
s
1
s
=
s+as+a
( s + a )2
(311)
D
te
dt
(312)
0,
at
=e
at
ate
at
at .
s2 + s 2
s2 + s 2
=
s3 + 5s2 + 7s + 3
( s + 1)2 ( s + 3)
(313)
corresponding to the block diagram on the next page. We obtain the corresponding Inverse Laplace Transform by several successive convolution
integrals, as follows
1
1
#
#
L
(
t
)
=
e
e
(t) = te t .
(314)
(# + 1)2
1
(
t
)
=
#e
(# + 1)2 (# + 3)
Z t
3#
te t e 3(t
1 3t
=
e
+ 2te
4
(t)
t)
dt
e
(315)
95
+
1
I (s) +
1
s
1
s
1
s
+ O(s)
+
+
+
3
3e
3t
2te
+ 2e
4
d(t) + e
3e
3t
2te
4
+ 3e
(316)
and
3d(t) + 9e
3t
+ 2te
2e
+ 3d(t)
3e
9e
3t
+ 2te
4
5e
(317)
3t
+ 2te
1
2
e
4
t
3t
5e
+ 2te
1
3e 3t 2te t + 3e t
4
e t = e 3t te t
(318)
1
1
X ( s ) = O ( s ),
O(s) =
I (s) 3O(s) 2X (s)
(319)
s
s
from which we obtain
s2 O(s) + 3sO(s) + 2O(s) = sI (s).
(320)
s
s2 + 3s + 2
(321)
96
X (s)
I (s) +
97
1
s
1
s
O(s)
We obtain the corresponding unit impulse response by computing the distributional derivative of the unit impulse response corresponding to the
transfer function 1/ s2 + 3s + 2 .
Z t
o (t) = c +
i (t ) dt = cu(t) + i (#) u(#) (t)
(322)
0
c
I (s)
c + I (s)
+
=
.
s
s
s
(323)
Frequency domain:
Time domain:
c
i(t)
c+
Rt
o (t)
I (s) +
1
s
O(s)
Example 51. The figure on the next page shows two representations of a
nonzero-initial-value integrator with negative feedback. In this case,
O(s) =
1
c + I (s)
s
aO(s)
(324)
Time domain:
i(t)
c+
Rt
o (t)
a
Frequency domain:
c
I (s)
1
s
O(s)
c + I (s)
.
s+a
(325)
When I (s) = 0, the output equals the free response. When c = 0 and
I (s) = 1, the output equals the unit impulse response. The free response is
always a multiple of the unit impulse response and equals the unit impulse
response when c = 1.
Example 52. The figure on the next page shows two representations of a
double nonzero-initial-value integrator with negative feedback. In this case,
O(s) =
1
1
d + X (s) , X (s) =
c + I (s)
s
s
aX (s)
ds + c + da + I (s)
s2 + as + b
bO(s)
(326)
(327)
When I (s) = 0, the output equals the free response. When c = 0 and
I (s) = 1, the output equals the unit impulse response. The free response is a
linear combination of the unit impulse response of the transfer function
1
s2 + as + b
98
(328)
and its distributional derivative, and equals the unit impulse response when
c = 1 and d = 0.
Time domain:
b
i(t) +
Rt
c+
x (t)
d+
Rt
o (t)
Frequency domain:
b
c
I (s) +
99
1
s
X (s)
+
+
1
s
Exercises
1. Consider the mechanical suspension shown below.
f (t)i
xm
2 kg
i
3 N/m
0.1 Ns/m
O(s)
Use a free-body diagram and Newtons 2nd law to derive a differential equation governing x (t) and determine whether the system
is overdamped, underdamped, or critically damped. Graph the
corresponding unit impulse response.
2. Explain why the unit impulse response of an overdamped system
is a nonnegative function of time.
3. Use a frequency-domain block diagram representation of the
initial-value problem
do
(t) = e
dt
t/2
, o (0) = 1
u(t
1), o (0) = 1,
do
(0) = 3
dt
1.
7. Determine the initial values for which the free response of a double integrator with negative feedback equals the unit impulse
response.
8. Find the inverse Laplace transform of the function
4s3
4s2 + 4s + 5
4s2 + 4s + 5
=
.
2
+ 12s + 9s + 2
(2s + 1)2 (s + 2)
100
Solutions
1. Consider the following free-body diagram, where x = dx/dt.
f (t)i N
3xi N
0.1x i N
d2 x
dx
(t) + 0.1 (t) + 3x (t) = f (t)
dt
dt2
F (s)
1
2
1
s
1
s
X (s)
5
100
It follows that
F (s)
1
,
2
2 s + 0.05s + 1.5
which is underdamped. The graph of the corresponding unit
impulse response
p
1 t/40 sin( 2399/40t)
p
h(t) = e
2
2399/40
X (s) =
t
30
0.5
60
90
101
eg1 t
g1
eg2 t
.
g2
it follows that
h
L e
i
cos 2p# (s) =
g2 )t
s+1
,
(s + 1)2 + 4p 2
1
s +1
(s+1)2 +4p 2
1
s
O(s)
O(s) =
4. Since,
h
L #e
it follows that
#/2
(s) =
4
,
(2s + 1)2
1
4
(2s+1)2
1
s
O(s)
2
or, in other words,
O(s) =
5. Since
it follows that
L u (#)
4s2 + 4s + 5
.
4s3 + 12s2 + 9s + 2
u (#
1
1) ( s ) =
e
s
= 1,
102
2
3
1 e
s
X (s)
1
s
1
s
O(s)
3
or, in other words,
O(s) =
s2 + 6s + 1 e
s3 + 3s2 + 2s
1
1
0
+
+
3
X (s)
1
s
1
s
O(s)
2
It follows that
O(s) =
Since
L #e
o (t) = 5te
3s + 5
.
( s + 1)2
(s) =
1
,
( s + 1)2
D
te
dt
+3
= 2te
+ 3e t .
103
The free response (when I (s) = 0) equals the unit impulse response (when c = d = 0 and I (s) = 1), provided that
ds + c + da = 1
for all s, i.e., c = 1 and d = 0.
8. Since
L
it follows that
1
L
1
1
(t) = te
2
4
(2# + 1)
1
(t) =
(2# + 1)2 (# + 2)
1
#e
4
1
e
9
2t
t/2
#/2
2#
3t 2
e
18
(t)
t/2
4#3
4#2 + 4# + 5
D
D2
(t) = 5h(t) + 4 h(t) + 4 2 h(t)
2
dt
+ 12# + 9# + 2
dt
13 2t 6t 4 t/2
= e +
e
.
9
9
Since h(0) = dh
dt (0) = 0, both distributional derivatives are simply
regular derivatives.
9. Since
L
it follows that
1 + sb
L 1
(t) = e
s+a
1
(t) = e
#+a
at
+b
D
e
dt
at
at
= bd(t) + (1
ab)e
at
104
105
4
c
I (s) +
+
+
Here,
X (s) =
X (s)
1
s
1
c + I (s)
s
+
+
4O(s) , O(s) =
1
s
O(s)
1
d + X (s)
s
ds + c + I (s)
.
s2 + 4
(329)
(330)
Since
L sin(w#) (s) =
s2
it follows that
w
s
, L cos(w#) (s) = 2
,
2
+w
s + w2
c
I (s)
O(s) = dL cos 2# (s) + L sin 2# (s) + 2
.
2
s +4
(331)
(332)
c
sin 2t
2
(333)
d2 x
dx
dx
(t) + b (t) + kx (t) = 0, x (0) = x0 ,
(0) = v0
2
dt
dt
dt
(334)
106
Free-body diagram:
x
i
b
~Fs =
kxi ~Fd =
b dx
dt i
Z t
0
1
m
Z t
0
by(t ) + kx (t ) dt
(335)
bY (s) + kX (s)
s
(336)
translate into
X (s) =
v0
x0 + Y ( s )
, Y (s) =
s
1
m
b
m x0
,
b
k
ms + m
x0 s + v0 +
s2 +
Y (s) =
v0 s
s2 +
k
m x0
b
k
ms + m
107
(337)
1
s2 + 2
O1 (s)
O2 (s)
O2 (s)
(340)
O1 (s)
ds + c + I1 (s)
f s + e + I2 (s)
(341)
(342)
2
c
I1 (s) +
+
+
I2 (s)
+
+
d
1
s
1
s
+
+
+
+
1
s
1
s
O1 (s)
O2 (s)
e
2
When I1 (s) = I2 (s) = 0, the outputs o1 (t) and o2 (t) are harmonic
with angular frequency w, provided that
O1 (s) =
o11 s + o12
o s + o22
, O2 (s) = 212
2
2
s +w
s + w2
(343)
for some coefficients o11 , o12 , o21 , and o22 . Substitution into (342),
multiplication on both sides by s2 + w 2 , and the requirement that
coefficients of equal powers of s be identical on both sides of the
equation, shows that o11 = d, o12 = c, o21 = f , and o22 = e. In this
case,
!
!
1
O1 (s)
ds + c
= 2
(344)
s + w2
O2 (s)
fs + e
108
1
2
s +2
(345)
=
.
s + w2
1
s2 + 2
fs + e
fs + e
(346)
1
s2 + 2
= ( s2 + 3
u)(s2 + 1
u) = 0,
(347)
T
, we find
s t
and
s2 + 2
s2 + 2
( s2 + 3) s + t = 0 ) t = s
( s2 + 1) s + t = 0 ) t =
s,
(348)
(349)
respectively.
We conclude that the system of two coupled double integrators
p
with feedback have two natural frequencies: w = 3 and w = 1. If
the initial values are chosen so that f = d and e = c, associated with
the first eigenvector, then
O1 (s) = O2 (s) =
ds + c
,
s2 + 3
(350)
p
c
3t + p sin 3t
3
(351)
In this case, both outputs have identical amplitude and phase shifts.
Both outputs reach their maximal and minimal values simultaneously.
If, instead, the initial values are chosen so that f = d and e = c,
associated with the second eigenvector, then
O1 (s) =
O2 (s) =
ds + c
,
s2 + 1
(352)
(353)
109
In this case, the two outputs have identical amplitude, but phase
shifts that differ by p. When one of the outputs reaches its maximal
value, the other reaches its minimal value.
For arbitrary e and f ,
!
!
!
(d+ f )s+(c+e)
(d f )s+(c e)
ds + c
2
2
=
+
,
(354)
(d+ f )s+(c+e)
( f d)s+(e c)
fs + e
2
2
where the first column matrix is an eigenvector of
!
s2 + 2
1
1
s2 + 2
(355)
corresponding to the eigenvalue s2 + 3, and the second is an eigenvector corresponding to the eigenvalue s2 + 1. It follows by linearity
that the free response of this dynamical system is a combination of
p
a harmonic signal with angular frequency 3 and a harmonic signal with angular frequency 1. Such a combination is typically not
p
periodic, since 3 is irrational.
Each natural frequency is associated with a mode shape, corresponding to the relationship between the amplitudes and
phase shifts of the harmonic signals associated with each output. In the example, the mode shapes are 1 : 1 and 1 : 1,
corresponding to an in-phase and an out-of-phase mode, respectively.
d2 x1
( t ) + k 1 x1 ( t ) =
dt2
K x1 ( t )
x2 (t) , x1 (0) = d,
dx1
(0) = c
dt
(356)
m2
d2 x2
( t ) + k 2 x2 ( t ) =
dt2
K x2 ( t )
x1 ( t ) , x2 (0) = f ,
dx1
(0) = e
dt
(357)
(358)
(359)
110
x2
m2
K
x1
m1
k2
k1
k1 +K
m1
K
m2
K
m1
s2 + k2m+2K
X1 ( s )
X2 ( s )
ds + c
fs + e
(360)
x11 s + x12
x s + x22
, X2 (t) = 212
s2 + w 2
s + w2
(361)
for some coefficients x11 , x12 , x21 , and x22 . Substitution into (360),
multiplication on both sides by s2 + w 2 , and the requirement that
coefficients of equal powers of s be identical on both sides of the
equation, shows that x11 = d, x12 = c, x21 = f , and x22 = e. In this
case,
!
!
1
X1 ( s )
ds + c
= 2
(362)
s + w2
X2 ( s )
fs + e
must be a (formal) eigenvector of
s2 +
k1 +K
m1
K
m2
K
m1
s2 + k2m+2K
(363)
K/m1
s2 + (k2 + K )/m2
= 0,
111
(364)
It follows that v = l
(k1 + K )/m1
K/m2
K/m1
(k2 + K )/m2
= 0,
(365)
(366)
(367)
t2
k1 + K
k2 + K 2
4K2
d=
+
(368)
4
m1
m2
m1 m2
which is always positive.
It follows that l s2 is real and positive for both eigenvalues
of (363). If we denote these values by w12 and w22 , we conclude that
the mechanical system always has two natural frequencies. Since
the corresponding eigenvectors are also real, it follows that for each
modal oscillation, x1 (t) and x2 (t) are harmonic signals with different
amplitudes but identical phase shifts, if both components of the
eigenvector have the same sign (in-phase), or harmonic signals with
different amplitudes and phase shifts that differ by p (out-of-phase),
if the components of the eigenvector are of opposite sign.
For arbitrary e and f , the column matrix
!
ds + c
(369)
fs + e
may be expressed as a sum of an eigenvector corresponding to w1
and an eigenvector corresponding to w2 . It follows by linearity that
the free response of the system is a sum of modal oscillations with
angular frequencies w1 and w2 . Such a combination is periodic only
if w1 and w2 are rationally related, i.e., if there are integers m and n
such that mw1 = nw2 .
Exercises
1. Suppose that
5s + 6 kX (s)
.
s2 + 5
Determine k such that x (t) is harmonic with angular frequency 3.
X (s) =
112
X (s)
1
s
O(s)
c
and express x (t) in terms of the initial values c and d.
3. Show that the dynamical system shown below has no natural
frequencies.
c
1
s
X (s)
1
s
O(s)
9
4. Determine x0 and v0 such that the solution to the initial-value
problem
m
d2 x
dx
(t) + kx (t) = 0, x (0) = x0 ,
(0) = v0
dt
dt2
b
c
I (s) +
+
+
d
1
s
X (s)
+
+
1
s
O(s)
Let w denote its natural frequency. Show that the transfer function
of the corresponding system with zero initial values has poles at
jw.
113
3x2 (t) = 0
3x1 (t) = 0
114
x2
2 kg
4 N/m
x1
2 kg
4 N/m
i
10 N/m
11. Suppose that the free response of the mechanical system shown
below to some choice of initial conditions is harmonic in both
outputs with period approximately equal to 1.66 s. Sketch the configuration of the system when the first mass reaches its maximal
displacement in either direction, as well as when the first mass
is halfway between its two maximal displacements, and indicate
the elapsed time since the initial time at which both displacement
equal 0.
x2
4.1 kg
5.2 N/m
x1
0.5 kg
12.6 N/m
i
1.3 N/m
k 1 x1 ( t ) + K ( x2 ( t )
x1 (t)) + ad(t)
m2 x2 (t) =
k 2 x2 ( t )
x1 (t)) + bd(t)
K ( x2 ( t )
115
2.
Solutions
1. If x (t) is harmonic with angular frequency 3 then
X (s) =
ds + c
s2 + 9
1
c
s
s2
5s + 6
) k = 4.
+5+k
9O(s) , O(s) =
1
d + X (s)
s
cs 9d
.
s2 + 9
3d sin 3t and the natural frequency
X (s) =
It follows that x (t) = c cos 3t
is 3.
3. Here,
X (s) =
1
c
s
9X (s) , O(s) =
1
d + X (s)
s
c
ds + c + 9d
, O(s) =
s+9
s2 + 9s
r
r
k
m
k
x (t) = x0 cos
t + v0
sin
t
m
k
m
q
is harmonic with amplitude x02 + v20 m/k and phase shift q,
where
p
x0
v0 m/k
cos q = q
, sin q = q
.
x02 + v20 m/k
x02 + v20 m/k
p
p
It follows that x0 = 2 and v0 = 2k/m.
116
5. Here, w =
b. Furthermore,
O(s) =
ds + c + I (s)
s2 + b
s2
1
+b
( s2 + 3
u)(s2 + 1
T
satisfy
s t
s2 + 3
( s2 + 2 +
s2 + 3
( s2 + 2
u)
b.
2=0
3) s
2t = 0 ) t =
3) s
2t = 0 ) t =
p
2
and
p
1+ 3
s,
2
respectively.
7. If l is an eigenvalue then u = l satisfies the characteristic equation
s2 + a
g
This implies that v = l
u
s2
b
+d
= 0.
a
g
b
d
= 0,
v is a root of the
as a function of s. In the special case that v = w 2 and, consequently, l = s2 + w 2 , it follows that jw are two roots of the
determinant.
8. If they exist, the natural frequencies correspond to eigenvalues of
the matrix
!
s2 + b
a
a
s2 + b
117
1
1
d + c + 3X2 (s)
s
s
1
1
X2 ( s ) =
f + e + 3X1 (s)
s
s
X1 ( s ) =
4X1 (s)
12X2 (s)
3
2
s + 12
X1 ( s )
X2 ( s )
ds + c
fs + e
dx2
1
where c = dx
dt (0), d = x1 (0), e = dt (0), and f = x2 (0). The
natural frequencies w1 and w2 correspond to eigenvalues of the
matrix
!
s2 + 4
3
3
s2 + 12
s2
3
+ 12
= 0.
T
For the corresponding eigenvectors s t
,
and
s2 + 4
s +4
s2 + 3
s + 13
3t = 0 ) s = 3t
3t = 0 ) t =
3s
118
3
2
x1 (t) = 3 cos
3t
x2 (t) = cos
3t
2p
p
1
3p
2
3
Similarly, w2 corresponds to an out-of-phase mode shape 1 :
for which
p
x1 (t) = cos 13t
p
x2 (t) = 3 cos 13t
3,
x2 ( t ) =
3 cos
13t
1
x1 (t) = cos
p
1
13t
2p
3p
2
3
10. Here,
ds + c + 12 F1 (s) + 2 X2 (s)
s2 + 5
1
f s + e + 2 F2 (s) + 2 X1 (s)
X2 ( s ) =
s2 + 2
X1 ( s ) =
X1 ( s )
X2 ( s )
,
.
119
T
equivalently, that the column matrix 1 k
is an eigenvector
of
!
s2 + 7
2
2
s2 + 4
An eigenvalue l of this matrix is a solution u = l of the characteristic equation
s2 + 7 u
2
=0
2
2
s +4 u
i.e., l = s2 + 3 or l = s2 + 8. For the corresponding eigenvectors
T
,
s t
2
2
s +7
s +3
s 2t = 0 ) t = 2s
and
s +7
s +8
2t = 0 ) t =
1
s
2
1/2.
x1
x2
x2
t = 0 s:
t 0.42 s:
x1
120
x2
x2
t 0.83 s:
t 1.25 s:
x1
x1
ds + c0 + K ( X2 X1 )/m1
f s + e0 K ( X2 X1 )/m2
, X2 =
2
s + k1 /m1
s2 + k2 /m2
or, equivalently,
s2 +
k1 +K
m1
K
m2
K
m1
s2 + k2m+2K
X1
X2
ds + c0
f s + e0
and
K
m1
k2 +K
m2
ds + c0
f s + e0
is a corresponding eigenvector.
With the given numerical values for m1 , m2 , k1 , k2 , and K, the
eigenvalues are 8 and 3 with corresponding eigenvectors
!
!
2
1
,
1
2
121
2 f cos
p
p
p
2e0
e0
p sin 8t, x2 (t) = f cos 8t + p sin 8t
8
8
8t
p
p
p
c0
2c0
3t + p sin 3t, x2 (t) = 2d cos 3t + p sin 3t
3
3
If neither of these relations hold for the initial conditions, then the
response is a sum of the two modal oscillations, i.e.,
p
p
p
2B
D
p sin 8t + C cos 3t + p sin 3t
8
3
p
p
p
p
B
2D
x2 (t) = A cos 8t + p sin 8t + 2C cos 3t + p sin 3t
8
3
where
x1 ( t ) =
2A cos
8t
2A + C = d,
2B + D = c0 , A + 2C = f , B + 2D = e0
2
3
1
11
,B =
,C = ,D =
5
10
5
10
1
3
, B = 0, C = , D = 0
5
5
p
p
2
3
cos 8t + cos 3t,
5
5
p
p
1
6
x2 ( t ) =
cos 8t + cos 3t.
5
5
p
In particular, for 8t = p/2 + np, x1 : x2 = 1 : 2, i.e., according to
p
the second mode shape. Similarly, for 3t = p/2 + np, x1 : x2 =
2 : 1, i.e., according to the first mode shape.
x1 ( t ) =
2
, C = 0, D =
5
6
5
p
p
4
6
p sin 8t
p sin 3t,
5 8
5 3
p
p
2
12
p sin 3t
x2 (t) = p sin 8t
5 8
5 3
x1 ( t ) =
122
p
In particular, for 8t = np, x1 : x2 = 1 : 2, i.e., according to the
p
second mode shape. Similarly, for 3t = np, x1 : x2 = 2 : 1, i.e.,
according to the first mode shape.
123
q)
124
125
M
f (t)
t
|o (t)| =
Z t
0
t ) dt M
h ( t )i ( t
Z t
0
|h(t )| dt,
(370)
|h(t )| dt < .
(371)
H ( s ) = L h (#) ( s ) =
| H (s)|
st
h(t)e
dt
(372)
Z
0
|h(t)| dt < .
(373)
Rt
0
o (t)
126
1
s+a
(374)
at
(375)
1
<
a
(376)
Z
0
|e
at
| dt =
Example 56. Consider the sequential wiring of two identical zero-initialvalue integrators with negative feedback shown below.
Figure 67: A sequential wiring of two
identical zero-initial-value integrators
with negative feedback.
a
i(t) +
y(t) +
Rt
0
Rt
o (t)
1
( s + a )2
(377)
Z
0
|te
at
at
| dt =
(378)
1
<
a2
(379)
jwt
dt,
(380)
it follows that
o (t) ejwt
Z
0
h(t )e
jwt
dt = H (jw )ejwt
(381)
for large t. We call this the steady-state behavior of the output and
write oss (t).
127
The complex-valued function w 7! H (jw ) is known as the frequency response function. Since the magnitude of ejwt equals
1, the magnitude of the steady-state output oss (t) = H (jw )ejwt
equals the magnitude of the frequency-response function.
A pair of graphs of the magnitude (in dB10 ) and phase, respectively,
of H (jw ) versus w is called a Bode diagram. Each graph is a Bode
plot.
10
1
jw + 1
(382)
(383)
1 + w2
and phase
arctan w.
(384)
0
dB
p 1
1+ w 2
20
40
10
0
10
100
101
102
101
102
arctan w
p/4
p/2
10
10
100
128
1
=
(jw + 1)2
with magnitude
p
and phase
w2
1
+ 2jw + 1
(385)
(1
arctan
if w < 1,
arctan
(386)
w 2 )2 + 4w 2
2w
1 w2
2w
w2 1
(387)
(388)
0
dB
20
40
1
(1 w 2 )2 +4w 2
60
80
10
0
10
100
arctan 1 2ww2
arctan w2w
p
2 1
p/4
102
101
w<1
w>1
p/2
3p/4
p
10
10
100
101
102
d2 x
dx
(t) + b (t) + kx (t) = f (t)
dt
dt2
(389)
Free-body diagram:
f (t)i
f (t)i
i
d
kxi
b dx
dt i
g1 #
g2 #
h(t) = e e
(t)
(391)
where
g1,2 =
b2
2m
4mk
(392)
mw 2
1
1
=
k1
+ jbw + k
1
+ j2zW
W2
(393)
p
where wn := k/m, z := b/(2mwn ), and W = w/wn . As shown in
the Bode plot on the next page, its magnitude
1
p
k (1
1
W2 )2
+ 4z 2 W2
(394)
z2
129
(395)
130
2z 2
(396)
20
40
10
z 2 < 1/2
z 2 > 1/2
dB
1
(1 W2 )2 +4z 2 W2
10
101
i
k
f (t)
Let the unit vector i be parallel to the line of action of the suspension. Denote by x the displacement of the vehicle body of mass m
along i relative to the undeformed configuration of the spring, when
the lower end of the suspension is located at some reference level.
Suppose that displacements along i of the lower end of the suspension relative to the reference level are described by the function f (t).
Figure 72: A quarter-var vehicle suspension model that accounts for uneven
terrain through the function f (t).
d2 x
dx
df
m 2 (t) = b
(t)
(t)
dt
dt
dt
k x (t)
f (t)
(397)
1
m
Z t
0
Z t
0
b g(t )
Z t
0
y(t ) + k f (t )
g(t ) dt,
x (t )
(399)
dt.
(400)
x (0) + Y ( s )
f (0) + G ( s )
, F (s) =
,
s
s
(401)
and
where X (s) = L x (#) (s), Y (s) = L y(#) (s), F (s) = L f (#) (s), and
bs + k
F (s)
+ bs + k
ms2
(403)
i.e., the relationship between the displacement of the lower end of the
suspension and the displacement of the plate is a convolution with
transfer function
bs + k
H (s) =
.
(404)
2
ms + bs + k
The dynamic system governed by this transfer function is stable,
since the same is true for the transfer function 1/(ms2 + bs + k). The
magnitude of the steady-state response xss (t) to the complex exponential input f (t) = ejwt is given by the magnitude of the frequency
response function, i.e.,
p
p
k 2 + b2 w 2
1 + 4z 2 W2
= p
(405)
| H (jw )| = p
2
2
2
2
(k mw ) + b w
(1 W2 )2 + 4z 2 W2
p
where wn := k/m, z := b/(2mwn ), and W = w/wn . For W 1,
1, | H (jw )| 2z/W.
| H (jw )| 1. For W
As shown in the Bode plot below, there exists a unique global
maximum at the resonance frequency
rq
wn
w = w :=
1 + 8z 2 1
(406)
2z
131
132
| H (jw )| = s
(1 + 8z 2 )1/4
p
1 + 1 + 8z 2 + 1
1+
1+8z 2
4z 2
(407)
0
dB
10
20
10
1+4z 2 W2
(1 W2 )2 +4z 2 W2
10
101
d2 x
dx
(t) + b (t) + kx (t) = mdw 2 cos wt
dt
dt2
(408)
For design purposes, it may be desirable to select values of the stiffness and damping that ensure that x (t) remains below a critical
bound during steady-state operation across a range of frequencies.
The response x (t) equals the real part of the solution of the equation
d2 x
dx
M 2 (t) + b (t) + kx (t) = mdw 2 ejwt
(409)
dt
dt
Here, x (t) is the output of a double integrator with negative feedback
and gains b/M and k/M that results from the input mdw 2 ejwt /M. It
wt
i
b
follows that the magnitude of the steady-state response xss (t) equals
| H (jw )| = p
mdw 2
(k
Mw 2 )2 + b2 w 2
W2
md
p
M (1
W2 )2 + 4z 2 W2
(410)
p
where wn := k/M, z := b/(2Mwn ), and W = w/wn .
p
As shown in the Bode plot below, if z > 1/ 2, the magnitude
increases monotonically from 0 at w = 0 toward a limiting value of
p
md/M as w ! . In contrast, if z < 1/ 2, there exists a global
maximum at the resonance frequency
w = w := p
wn
1
(411)
2z 2
| H (jw )| =
md
1
p
M 2z 1
z2
(412)
The resonance peak in the corresponding Bode plot shifts to the right
from the natural frequency wn as the damping increases.
133
20
z 2 < 1/2
dB
z 2 > 1/2
20
40
10
134
W2
(1 W2 )2 +4z 2 W2
10
x1 , x2
i
z/2
f (t)
101
z/2
f (t)
1
k/2
1/2
1/2
k/2
1/
tions
1 dy1
(t) =
dt
k x1 ( t )
f (t)
z y1 ( t )
dy2
( t ) = x1 ( t )
dt
where
x1 ( t )
x2 ( t )
(413)
y2 ( t )
x2 ( t ) + z y1 ( t )
y2 ( t )
(414)
dx1
dx2
( t ), y2 ( t ) =
( t ).
(415)
dt
dt
In the case of zero initial displacements and velocities, this corresponds to the system of coupled integrators with feedback shown
below.
Here,
1
1
X1 (s) = Y1 (s), X2 (s) = Y2 (s)
(416)
s
s
y1 ( t ) =
F (s)
1
s
Y1 (s)
X1 ( s )
+
+
and
1
s
Y1 (s) =
k X1 ( s )
s
Y2 (s) =
1
s
X1 ( s )
Y2 (s)
1
s
X2 ( s )
+ z Y1 (s) Y2 (s) ,
X2 ( s )
F ( s ) + X1 ( s )
(417)
(418)
where X1 (s) = L x1 (#) (s), X2 (s) = L x2 (#) (s), Y1 (s) = L y1 (#) (s),
s4
+ z (1 + ) s3
135
(421)
(422)
are both stable. It follows that the ratio of the magnitudes of the
steady-state responses x1,ss (t) and x2,ss (t) to the complex exponential
input f (t) = ejwt is given by
p
| H1 (jw )|
(1 w 2 )2 + z 2 w 2
p
=
.
(423)
t :=
| H2 (jw )|
1 + z 2 w2
As seen in the diagram on the next page, this has a global minimum
at
rq
1
w = w :=
1 + 2z 2 1.
(424)
z
136
(1 w 2 )2 + z 2 w 2
1+ z 2 w 2
Exercises
1. Show that convolution with h(t) = u(t) is not stable.
2. Show that the convolution corresponding to the transfer function
H (s) = s/(s + 1) is stable.
3. Show that the convolution corresponding to the transfer function
H (s) = 1/(s2 + bs + 1) is stable for all b > 0.
4. Show that the convolution corresponding to the transfer function
H (s) = (s2 + 1)/(s4 + zs3 + 3s2 + zs + 1) is stable for all z > 0.
5. Suppose that oss (t) is the steady-state response of a stable convolution with transfer function H (s) to the input i (t) = ejwt . Find the
steady-state response to the input i (t) = e jwt .
6. Suppose that oss (t) is the steady-state response of a stable convolution with transfer function H (s) to the input i (t) = ejwt . Find the
steady-state response to the input i (t) = cos wt.
7. Describe the steady-state response of the stable convolution with
transfer function
s 1
H (s) = 2
s +s+4
to a general harmonic input, and specialize to the case that i (t) =
p
2 sin 3t.
8. Draw a Bode diagram for the transfer function H (s) = s/(s + 1).
9. Draw a Bode diagram for the transfer function H (s) = s/(s + 1)2 .
10. Show that the frequency response function corresponding to the
transfer function
1
H (s) =
ms2 + bs + k
has a resonance peak only if the dynamical system is underdamped.
11. Suppose that a convolution is described by the Bode diagram
shown below.
0
dB
20
40
60
1
100
101
102
100
101
102
cos 10t.
10
p
p/2
p
p/2
p
10
Solutions
1. Here,
o (t) =
where s = t
Z t
0
u ( t )i ( t
t ) dt =
Z t
0
i (s ) ds,
is not bounded.
Z t
0
u(s ) ds = t
137
2. Here,
h(t) =
and
Z t
0
h ( t )i ( t
D
e
dt
= d(t)
t ) dt = i (t)
h ( t )i ( t
t ) dt |i (t)| +
Z t
0
Z t
0
i(t
i(t
t ) dt
t ) dt .
h ( t )i ( t
t ) dt M + M
Z t
0
dt 2M.
for all t.
3. Since s2 + bs + 1 = (s g1 )(s g2 ) for some complex numbers g1
and g2 corresponding to the poles of H (s), it follows that h(t) is a
convolution of eg1 t and eg2 t . In particular,
Z t
0
|h(t )| dt <
1
2
For b 1, it follows g(b) g(0) b/2, i.e., that the real part
is negative for small, but positive b. But g(b) is purely imaginary
only if b = 0, since
g(b) = jw )
w 2 + jbw + 1 = 0
|h(t )| dt <
138
g (0) =
p
5 5
< 0.
20
g2 (0) + 1
=
4g2 (0) + 6
jzw 3
3w 2 + jzw + 1 = 0
h(t )e
jw (t t )
jwt
dt = e
then equals
jwt
Z t
0
h(t )ejwt dt
jwt
Z
0
Z
0
h(t )e
jwt
dt
ejwt + e jwt
2
the steady-state response to i (t) = cos wt equals
cos wt =
| H (jw )| cos(wt + q ).
139
q + arg H ( jw )),
p
3t = 2 cos( 3t
p/2)
implies that
p
oss (t) = 2 cos( 3t
p
p/2 + p/3) = 2 sin( 3t + p/3).
p
The amplitude of the steady-state response equals 2| H ( j 3)| = 2.
The absolute value of the difference in phase between the input
p
and the steady-state response equals | arg H ( j 3)| = p/3.
In general,
| H ( jw )| = p
p
(4
1 + w2
w 2 )2 + w 2
12.0.
140
0
dB
10
20
10
10
101
101
102
jw
jw + 1
w
1 + w2
and phase
p
arctan w
2
A corresponding Bode diagram is shown below.
0
dB
p w
1+ w 2
20
40
10
10
100
p/2
p
2
arctan w
p/4
0
10
10
100
jw
(jw + 1)2
101
102
141
has magnitude
w
1 + w2
and phase
p
2 arctan w
2
A corresponding Bode diagram is shown below.
0
dB
w
1+ w 2
20
40
10
10
100
101
102
p/2
p
2
p/4
2 arctan w
0
p/4
p/2
10
10
100
101
102
1
mw 2 + jbw + k
has magnitude
p
(k
mw 2 )2 + b2 w 2
(k
mw 2 )
b2
mw 2 )2 + b2 w 2
3/2
142
11. From the Bode diagram, the magnitude and phase of H (j3) are
approximately 0.56 and p. Similarly, the magnitude and phase of
H (j10) are approximately 0.1 and p/2. By linearity, it follows
that the steady-state response is approximately equal to
0.56 cos(3t + p )
0.1 cos(10t
p/2)
143
144
145
q (t)
`
i
A
Introduce a coordinate system such that the stationary end point of the
rod coincides with the origin, and A coincides with the z-axis, and assume
that this is inertial. Other than the stationary end point of the rod, all particles on the rod move along circles centered at the origin. It follows that these
particles cannot be isolated, and must be affected by forces.
Let q (t) denote the time-dependent angle of rotation of the rod relative
to the x-axis. A rotation corresponds to a positive angle, when the direction
of rotation would result in motion along the z-axis of a right-handed screw,
and is negative otherwise. The position vector to the particle at distance r
146
(425)
All particles of the mechanism may be described this way by considering all
values of r on the interval [0, `]. Since q (t) is the only generalized coordinate required to describe the positions of all the particles in the mechanism,
the mechanism has a single degree of freedom.
Step 3: Associate mass with each particle, keeping in mind
that the total mass in the mechanism must be finite.
Step 4: Compute the accelerations of each particle in the
mechanism, keeping in mind that the unit vectors i, j,
and k do not change with time.
Example 60. If the rod in Example 59 is homogeneous and has total mass
M, then each particle is associated with an infinitesimal mass dm = dr,
in terms of a constant density per unit length and an infinitesimal
length dr. In particular,
Z
rod
dm =
Z `
0
dr = ` = M,
(426)
i.e., = M/`.
Figure 80: An infinitesimal element
R of
mass dm, such that the integral rod dm
equals the total mass of the rod.
dm = dr
r
~a :=
d2~r
=
dt2
(427)
d~Fr
g dm j
d~Fg =
i
z
In the free-body diagram above, the force vector d~Fr denotes an unknown,
infinitesimal reaction force associated with internal interactions within the
rod and/or with the hinge. Newtons 2nd law now becomes
dm ~a = d~Fg + d~Fr
(428)
g dm j + d~Fr .
(429)
rod
~r d~Fr = 0,
(430)
147
i.e., that the total torque of the unknown reaction forces about the hinge
axis vanishes. Then, since
k ~r i =
(431)
the equation
Z
rod
k ~r i
Z
r dm q 2 (t) sin q (t) q(t) cos q (t) k ~r j
rod
Z
Z
~r d~Fr ,
=
g dm k ~r j + k
(432)
rod
rod
can be simplified to
M
q (t)
`
It follows that
Z `
0
r2 dr =
M
cos q (t) g
`
Z `
0
r dr.
2`
q + g cos q = 0
3
is the equation of motion of the rod.
(433)
(434)
(435)
(436)
body
dm = M
(437)
(438)
148
149
r x i + ry j + rz k
dm
~r
x
i
z
d~Fr
i
z
d~Fe
150
body
d~Fr = ~0
(439)
body
body
d~Fr +
body
d~Fe
(440)
can be simplified to
M x (t)i + y (t) j + z (t)k = ~Fe ,
where
~Fe :=
body
(441)
(442)
(443)
q (t)
R
(444)
151
and
~a =
(445)
cylinder
dm =
Z `/2 Z 2p Z R
`/2 0
r dr dj dz = pR2 ` = M,
(446)
i.e., = M/pR2 `.
j
j + q (t)
dm = r dr dj dz
r
M
r dr dj dz
r q 2 (t) cos( j + q (t)) + q(t) sin( j + q (t)) i
pR2
= d~Fr + d~Fe .
(447)
cylinder
r d~Fr = 0
(448)
i.e., the net reaction torque about the joint axis vanishes. Since
k (r i) =
(449)
d~Fe
d~Fr
it follows that
M
q (t)
pR2
Z `/2 Z 2p Z R
`/2 0
r3 dr dj dz = k
cylinder
+ k
can be simplified to
where
r d~Fr
cylinder
MR2
q (t) = Te
2
Te = k
cylinder
r d~Fe
r d~Fe .
(450)
(451)
(452)
x (t)
x (t)
~r = x (t) + r cos j
i + r sin j
j.
(455)
R
R
152
153
K
R
a
The expression x (t)/R for the angle of rotation ensures that the
disk rolls without slipping along the incline. The mechanism has a
single degree of freedom.
j y
x (t)i
k
z
x (t)/R
x (t)
x (t)
x (t)
x (t)
~a = x (t) r
cos
j
sin
j
i
R
R
R
R2
2
x (t)
x (t)
x (t)
x (t)
r
sin
j
+
cos
j
j.
(456)
R
R
R
R2
where x = dx (t)/dt and x = d2 x (t)/dt2 . Each particle may be
associated with the infinitesimal mass
dm =
M
r dr dj.
pR2
(457)
154
disk
~p d~Fr = 0,
(458)
i.e., that the net reaction torque about the point of contact with the inclined plane along a direction parallel to the axis of rolling vanishes.
Figure 90: The construction of the
vector ~p from the contact point to
the infinitesimal mass element, and a
superimposed free-body diagram.
d~Fs
~r
x (t)i
dm =
M
r dr dj
pR2
~p
d~Fr
R j
d~Fg
a
Since
x (t)
R r sin j
k
R
x (t)
~p j = r cos j
k,
R
~p i =
the contribution
d~Fg = g dm(sin ai
(459)
(460)
cos a j)
(461)
disk
~p d~Fg
Mg
pR2
Z 2p Z R
0
disk
r cos j
~p d~Fg =
x (t)
R
a + R sin a r dr dj k
MgR sin a k.
(462)
(463)
disk
(464)
in terms of the spring stiffness K, since this equals the net torque of
the force Kx (t)i applied at the center of the disk about the contact
point.
Finally,
Z
disk
~p dm ~a
M
=
pR3
Z 2p Z R
0
x (t)
r x (t)r cos j
R
0
x (t)
2
2
x (t) r + R + 2rR sin j
dr dj k
R
(465)
2
disk
~p dm ~a =
3MR x (t)
k.
2
(466)
(467)
it follows that
3MR x (t)
=
2
(468)
Exercises
1. Find the equation of motion for a small bead of mass m that is able
to translate along the straight line y = kx + in a two-dimensional
inertial coordinate system under the influence of a spring of undeformed length ` and stiffness K, connecting the bead to the origin.
y
y = kx +
m
K, `
x
155
2. Find the equation of motion for a small bead of mass m that is able
to translate along the parabolic curve y = bx2 in a vertical inertial
coordinate system under the influence of gravity.
y = bx2
x
3. Find the equation of motion for a small bead of mass m that is
able to translate along a vertical circular curve of radius R that
is spinning about a vertical axis through its center with angular
velocity w. Consider the influence of gravity.
w
m
R
4. Find the equation of motion for a small bead of mass m, suspended in a plane from a taut string that is wrapped around a
vertical pulley of radius R, and moving under the influence of
gravity.
R
g
5. Find the equations of motion for two small beads of mass m1 and
m2 that are able to translate along the x- and y-axis, respectively,
of an inertial coordinate system under the influence of a spring of
undeformed length ` and stiffness K, connecting the two beads.
156
y
m2
K, `
m1
x
m
r
g
R
157
m2
m1
Solutions
1. Introduce the basis vectors i and j such that the position vector to
the bead equals ~r = x (t)i + (kx (t) + ) j and
~a = x (t)i + k x (t) j.
y
j
p
The distance to the origin equals k~r k = x2 (t) + (kx (t) + )2 . It
follows that the force of the spring on the bead equals
~Fs =
K (k~r k
`)
~r
k~r k
158
~Fr
i + k j
j
~Fs
x
~Fr (i + k j) = 0
since ~r/x (t) = i + k j is tangent to the line. From Newtons 2nd
law,
m~a = ~Fs + ~Fr
or, after dot multiplication with ~r/x (t),
(1 + k 2 )m x (t) =
K (k~r k
`)
(1 + k 2 ) x (t) + k
.
k~r k
2. Introduce the basis vectors i and j such that the force of gravity on
the bead equals ~Fg = mg j and such that the position vector to the
bead equals ~r = x (t)i + bx2 (t) j for some constant b. It follows that
159
i + 2bx (t) j
~Fr
~Fg
x
i
From Newtons 2nd law,
m~a = ~Fg + ~Fr
or, after dot multiplication with ~r/x (t),
2mgbx (t).
R cos q (t)k
i
x
q (t)
j
y
y
wt
x
160
~Fr
m
~Fg
x
From Newtons 2nd law,
m~a = ~Fg + ~Fr
or, after dot multiplication with ~r/q (t),
mR2 q(t)
4. The distance between the bead and the point where the string
leaves the disk equals r0 + Rq (t) for some constant r0 and in terms
of the angle q (t) shown in the figure below. It follows that
+ R sin q (t)
r0 + Rq (t) cos q (t) j
and the acceleration ~a may be obtained by differentiating twice
with respect to t.
r0 + Rq (t)
R q (t)
q (t)
m
i
161
~Fr
q (t)
mg j
i
From Newtons 2nd law,
m~a = ~Fr
mg j
mg sin q (t).
5. Introduce the basis vectors i and j such that ~r1 = x (t)i and ~r2 =
y(t) j and, consequently, ~a1 = x (t)i and ~a2 = y (t) j.
y
j
m2
K, `
m1
x
p
The distance between the beads equals k~r2 ~r1 k = x2 (t) + y2 (t).
It follows that the force of the spring on the first bead equals
~F2,r
~F1,s
~F2,s
i
~F1,r
x
162
x (t)
2
2
m1 x (t) = K
x (t) + y (t) ` p
2
x ( t ) + y2 ( t )
q
y(t)
m2 y (t) = K
x 2 ( t ) + y2 ( t ) ` p
x 2 ( t ) + y2 ( t )
~r M = r x i + (y(t) + ry ) j
for some constants r x and ry and such that the spring is undeformed when y(t) = 0. Suppose that the rotating imbalance is
located at the top of its circular motion at t = 0 such that its displacement may be described by the position vector
~rm = (rc,x
rw 2 cos wt) j
M
wt r
y
rc,x i + (rc,y + y(t)) j
163
motor
j d~FM,r =
j ~Fm,r
i.e., that the vertical component of the reaction force on the particle
is equal and opposite to the vertical component of the reaction
force on the motor. Finally,
Z
motor
j d~Fs =
Ky(t),
motor
j d~Fd =
by (t)
since these are the net vertical force components acting on the
motor from the spring and damper, respectively.
~Fm,r
dM
d~FM,r
d~Fs
d~Fd
x
i
From Newtons 2nd law,
rw 2 cos wt) + ( M
m)y (t) =
Ky(t)
by (t).
164
cos q (t) j
R
q (t)
q (t)
In the free-body diagram below, d~Fr denotes an infinitesimal reaction force that ensures that the ring retains its shape and mass
distribution and is restricted to rotary motion about the horizontal
axis. Physical observations suggest that
Z
ring
k ~r d~Fr = 0
i.e., that the net torque of the reaction forces about the axis of
rotation vanishes. With dm = M dj/2p and d~Fg = dm g j, it
further follows that
Z
ring
k ~r d~Fg =
Finally, by assumption
Z
ring
k ~r d~Ft = t (t).
y
j
i
~r
d~Ft
k
dm
d~Fr
d~Fg
165
R cos q (t) j
~r2
m2
k
R
q (t)
m1
In the free-body diagrams below, ~Fr,1 and ~Fr,2 denotes reaction
forces that ensure that the beads remain at distance R from each
other, that the motion of the first bead remains in the vertical
plane, and the motion of the second bead remains along the axis.
~Fr,2
j
Kx (t)i
m2
k
m2 g j
~Fr,1
m1
m1 g j
~r2 ~Fr,1 = 0
166
i.e., the torque of the reaction force on the first bead about the axis
of rotation through the second bead vanishes. Moreover,
i ~Fr,1 + ~Fr,2 = 0
i.e., the horizontal component of ~Fr,2 is equal in magnitude and
opposite in direction to the horizontal component of ~Fr,1 .
From Newtons 2nd law,
m1~a1 =
m1 g j + ~Fr,1 , m2~a2 =
m2 g j
Kx (t)i + ~Fr,2
m1 gR sin q (t)
and
m1 Rq(t) cos q (t)
Kx (t).
167
168
Mq 2 (t)
~v ~v dm =
2`
rod
Z `
0
r2 dr =
M `2 q 2 (t)
.
6
(470)
169
d
T
T
= Fi
(471)
dt q i (t)
qi (t)
for i = 1, . . . , n. Here, q i denotes the derivative of qi (t) with respect
to time. Care should be taken when differentiating the kinetic energy
with respect to qi (t) and q i (t).
Example 63. For the rod in Example 59, an instantaneous infinitesimal
change dq to the value of the generalized coordinate q (t) results in an instantaneous infinitesimal change
d~r =
~r
dq = ( r sin q (t)i + r cos q (t) j) dq
q (t)
(472)
g dm j
d~Fg =
i
z
d~Fr
170
gravity
dW =
rod
Mg cos q (t)
`
d~Fg d~r =
`
0
r dr dq
(473)
Fq =
d
T
dt q (t)
T
= Fq
q (t)
(475)
(476)
With the expression for the kinetic energy in Example 62, we obtain
M `2 q 2 (t)
6
q (t)
M `2 q 2 (t)
d M `2 q (t)
M `2 q(t)
=
=
.
q (t)
6
dt
3
3
(477)
The equation of motion becomes
d
dt
M `2 q(t)
=
3
(478)
2`q(t)
+ g cos q (t) = 0
3
which is identical to the expression in Example 61.
(479)
(480)
(481)
171
r x i + ry j + rz k
body
dm = M
(482)
equals the total mass of the body. It follows that the total kinetic
energy equals
T=
1
2
body
~v ~v dm =
M 2
x (t) + y 2 (t) + z 2 (t) .
2
(483)
~r
~r
~r
dx +
dy +
dz = dx i + dy j + dz k
x (t)
y(t)
z(t)
(484)
to each particles position vector. The total infinitesimal work performed by the non-reaction forces then equals
Z
Z
~
~
dW =
d Fe d~r =
d Fe d~r
(485)
body
body
~Fe :=
body
(486)
(487)
172
173
and z(t) are the coefficients of dx, dy, and dz, respectively, in this
expression for the total infinitesimal work, i.e.,
Fx = Fe,x , Fy = Fe,y , Fz = Fe,z
(488)
Figure 93: A free-body diagram including known and unknown forces.
d~Fr
d~Fe
d
T
T
= Fx ,
dt x (t)
x (t)
d
T
T
= Fy ,
dt y (t)
y(t)
d
T
T
= Fz .
dt z (t)
z(t)
But,
d
dt
M
x (t) 2
x (t)
(489)
(490)
(491)
M 2
x (t) + y 2 (t) + z 2 (t) = M x (t)
2
(492)
and similarly for y(t) and z(t). The equations of motion become
M x (t) = Fe,x , My (t) = Fe,y , M z (t) = Fe,z
(493)
as expected.
The assumption that the unknown reaction forces do not contribute to the generalized forces implies that
Z
Z
~
~
d Fr d~r =
d Fr d~r = 0
(494)
body
body
body
d~Fr = ~0
(495)
which agrees with the physical observation that the net reaction force
should vanish.
174
q (t)
R
k
~v =
(496)
(497)
j + q (t)
dm = r dr dj dz
r
disk
dm =
Z `/2 Z 2p Z R
`/2 0
r dr dj dz = pR2 ` = M,
(498)
1
2
disk
~v ~v dm =
M q 2 (t)
2pR2 `
`/2 Z 2p Z R
`/2
r3 dr dj dz
MR2 q 2 (t)
.
4
(499)
(500)
d~Fe
d~Fr
An instantaneous infinitesimal change dq to the generalized coordinate q (t) results in the instantaneous infinitesimal change
d~r =
(501)
(502)
disk
d~Fe d~r =
disk
d~Fe (k ~r ) dq
(503)
disk
(504)
175
which equals the total (non-reaction) torque Te about the joint axis.
According to Lagranges equations
d
T
T
= Fq .
(506)
dt q (t)
q (t)
But,
d
dt
MR2 q 2 (t)
4
q (t)
MR2 q 2 (t)
MR2
=
q ( t ).
q (t)
4
2
(507)
(508)
as expected.
The assumption that the unknown reaction forces do not contribute to the generalized forces implies that
Z
~
~
~r d Fr dq = 0
d Fr d~r = k
(509)
disk
disk
disk
r d~Fr = 0
(510)
which agrees with the physical observation that the net reaction
torque about the joint axis should vanish.
x (t)
x (t)
~r = x (t) + r cos j
i + r sin j
j.
(513)
R
R
176
177
K
R
a
The expression x (t)/R for the angle of rotation ensures that the
disk rolls without slipping along the incline. The mechanism has a
single degree of freedom.
j y
x (t)i
k
z
x (t)/R
r
x (t)
r
~v = x (t) 1 + sin j
i x (t) cos j
R
R
R
x (t)
j
R
(514)
where x = dx (t)/dt. Each particle may be associated with the infinitesimal mass
M
dm =
r dr dj.
(515)
pR2
from which we obtain the total kinetic energy
T=
1
2
disk
~v ~v dm =
3M 2
x (t).
4
(516)
178
d~Fs
~r
x (t)i
dm =
M
r dr dj
pR2
~p
d~Fr
R j
d~Fg
a
An instantaneous infinitesimal change dx to the generalized coordinate x (t) results in the instantaneous infinitesimal change
r
x (t)
r
x (t)
d~r = dx 1 + sin j
i dx cos j
j
(517)
R
R
R
R
or, in other words,
dx
R
to each particles position vector. Recall the vector identity
d~r =
(k ~p)
dx
k (~p d~F ) .
R
(518)
(519)
(520)
Since
x (t)
R r sin j
k
R
x (t)
~p j = r cos j
k,
R
~p i =
and
d~Fg = g dm(sin a i
cos a j)
(521)
(522)
(523)
it follows that
~p d~Fg =
Mg
r cos j
pR2
x (t)
R
a + R sin a r dr dj
Z
dx
~
~
~p d Fg
dWg =
d Fg d~r = k
R
body
body
(524)
(525)
dx
~p d~Fs
dWs = k
= Kx (t) dx
R
body
(526)
(527)
in terms of the spring stiffness K, since this equals the net work of the
force Kx (t)i applied at the center of the disk. The generalized force
Fx corresponding to the generalized coordinate x (t) is the coefficient
of dx in the expression for the total infinitesimal work, i.e.,
Fx = Mg sin a
Kx (t).
d
T
T
= Fx
dt x (t)
x (t)
(528)
(529)
Kx (t)
(530)
which agrees with the equation of motion obtained using the Newtonian recipe.
The assumption that the unknown reaction forces do not contribute to the generalized forces implies that
Z
Z
dx
~
~
~p d F
d Fr d~r = k
=0
(531)
R
body
body
for arbitrary dx. This implies that
k
body
~p d~Fr = 0,
(532)
which agrees with the physical observation that the net reaction
torque about the contact point along a direction parallel to the axis of
rotation should vanish.
179
Exercises
1. Find the equation of motion for a small bead of mass m that is able
to translate along the straight line y = kx + in a two-dimensional
inertial coordinate system under the influence of a spring of undeformed length ` and stiffness K, connecting the bead to the origin.
y
y = kx +
m
K, `
x
2. Find the equation of motion for a small bead of mass m that is able
to translate along the parabolic curve y = bx2 in a vertical inertial
coordinate system under the influence of gravity.
y
y = bx2
x
3. Find the equation of motion for a small bead of mass m that is
able to translate along a vertical circular curve of radius R that
is spinning about a vertical axis through its center with angular
velocity w. Consider the influence of gravity.
w
m
R
4. Find the equation of motion for a small bead of mass m, suspended in a plane from a taut string that is wrapped around a
180
R
g
5. Find the equations of motion for two small beads of mass m1 and
m2 that are able to translate along the x- and y-axis, respectively,
of an inertial coordinate system under the influence of a spring of
undeformed length ` and stiffness K, connecting the two beads.
y
m2
K, `
m1
x
181
g
R
m2
m1
Solutions
1. Introduce the basis vectors i and j such that the position vector to
the bead equals ~r = x (t)i + (kx (t) + ) j and
~v = x (t)i + k x (t) j.
It follows that
T=
is the kinetic energy.
m
m x 2 (t)
~v ~v =
(1 + k 2 )
2
2
182
y
j
p
The distance to the origin equals k~r k = x2 (t) + (kx (t) + )2 . It
follows that the force of the spring on the bead equals
~Fs =
K (k~r k
`)
~r
k~r k
~Fr
y
j
~Fs
x
K (k~r k
`)
(1 + k 2 ) x (t) + k
dx.
k~r k
The generalized force Fx corresponding to the generalized coordinate x (t) is the coefficient of dx in this expression for the total
infinitesimal work, i.e.,
Fx =
K (k~r k
`)
(1 + k 2 ) x (t) + k
.
k~r k
d
T
T
= Fx .
dt x (t)
x (t)
183
But,
d
dt
m x 2 (t)
(1 + k 2 )
x (t)
2
m x 2 (t)
(1 + k 2 ) = m x (t)(1 + k 2 ).
x (t)
2
K (k~r k
`)
(1 + k 2 ) x (t) + k
.
k~r k
2. Introduce the basis vectors i and j such that the force of gravity on
the bead equals ~Fg = mg j and such that the position vector to the
bead equals ~r = x (t)i + bx2 (t) j for some constant b. It follows that
m
m x 2 (t)
~v ~v =
(1 + 4b2 x2 (t)).
2
2
y
The generalized force Fx corresponding to the generalized coordinate x (t) is the coefficient of dx in this expression for the total
infinitesimal work, i.e.,
Fx =
2mgbx (t).
184
~Fr
~Fg
x
i
According to Lagranges equations
d
T
T
= Fx .
dt x (t)
x (t)
But,
d
dt
m x 2 (t)
1 + 4b2 x2 (t)
x (t)
2
m x 2 (t)
1 + 4b2 x2 (t)
x (t)
2
2mgbx (t).
R cos q (t)k
mR2 2
q (t) + w 2 sin2 q (t) .
2
i
x
q (t)
j
y
y
wt
x
185
The generalized force Fq corresponding to the generalized coordinate q (t) is the coefficient of dq in this expression for the total
infinitesimal work, i.e.,
mgR sin q (t).
Fq =
k
~Fr
m
j
y
~Fg
x
According to Lagranges equations
d
T
T
= Fq .
dt q (t)
q (t)
But,
d
dt
mR2 2
q (t) + w 2 sin2 q (t)
q (t) 2
mR2 2
q (t) + w 2 sin2 q (t)
q (t) 2
= mR2 q(t) mR2 w 2 sin q (t) cos q (t).
186
4. The distance between the bead and the point where the string
leaves the disk equals r0 + Rq (t) for some constant r0 and in terms
of the angle q (t) shown in the figure below. It follows that
+ R sin q (t)
r0 + Rq (t) cos q (t) j
and the velocity ~v may be obtained by differentiating once with
respect to t. After simplification, the kinetic energy equals
T=
m
2
r0 + Rq (t) q 2 (t).
2
r0 + Rq (t)
R q (t)
q (t)
m
i
The generalized force Fq corresponding to the generalized coordinate q (t) is the coefficient of dq in this expression for the total
infinitesimal work, i.e.,
Fq =
~Fr
q (t)
mg j
187
d
T
T
= Fq .
dt q (t)
q (t)
But,
d
dt
m
2
r0 + Rq (t) q 2 (t)
q (t) 2
m
2
r0 + Rq (t) q 2 (t)
q (t) 2
2
= m r0 + Rq (t) q(t) + mR r0 + Rq (t) q 2 (t)
mg sin q (t).
5. Introduce the basis vectors i and j such that ~r1 = x (t)i and ~r2 =
y(t) j and, consequently, ~v1 = x (t)i and ~v2 = y (t) j. The total kinetic
energy then equals
T=
m1
m
m
m
~v ~v1 + 2 ~v2 ~v2 = 1 x 2 (t) + 2 y 2 (t).
2 1
2
2
2
y
j
m2
K, `
m1
x
p
The distance between the beads equals k~r2 ~r1 k = x2 (t) + y2 (t).
It follows that the force of the spring on the first bead equals
~F2,r
~F1,s
~F2,s
i
~F1,r
x
188
x (t) dx + y(t) dy
~
~
dW = F1,s d~r1 + F2,s d~r2 = K k~r2 ~r1 k `
.
k~r2 ~r1 k
The generalized forces Fx and Fy corresponding to the generalized
coordinates x (t) and y(t) are the coefficients of dx and dy in this
expression for the total infinitesimal work, i.e.,
x (t)
Fx = K k~r2 ~r1 k `
,
k~r2 ~r1 k
y(t)
Fy = K k~r2 ~r1 k `
.
k~r2 ~r1 k
According to Lagranges equations
d
T
T
= Fx ,
dt x (t)
x (t)
d
T
T
= Fy .
dt y (t)
y(t)
But,
d
dt
x (t)
d
dt
y (t)
m1 2
m2 2
x (t) +
y (t)
2
2
m1 2
m2 2
x (t) +
y (t) = m1 x (t)
x (t) 2
2
and
m1 2
m2 2
x (t) +
y (t)
2
2
m1 2
m2 2
x (t) +
y (t) = m2 y (t)
y(t) 2
2
x 2 ( t ) + y2 ( t )
` p
` p
x (t)
x 2 ( t ) + y2 ( t )
y(t)
x 2 ( t ) + y2 ( t )
,
.
189
~r M = r x i + (y(t) + ry ) j
for some constants r x and ry and such that the spring is undeformed when y(t) = 0. Suppose that the rotating imbalance is
located at the top of its circular motion at t = 0 such that its displacement may be described by the position vector
~rm = (rc,x
~v M = y (t) j, ~vm =
rw cos wt i + (y (t)
rw sin wt) j
m
2
y 2 (t) +
m 2
y (t) + r2 w 2
2
wt r
y
rc,x i + (rc,y + y(t)) j
190
center of the circular path of the rotating imbalance. An instantaneous infinitesimal change dy to the generalized coordinate y(t)
results in the instantaneous infinitesimal changes
d~r M = d~rm = dy j
to the position vectors to a particle on the motor and to the rotating imbalance. The total infinitesimal work performed by the
non-reaction forces then equals
dW =
motor
by (t) dy
Ky(t)
The generalized force Fy corresponding to the generalized coordinate y(t) is the coefficient of dy in this expression for the total
infinitesimal work, i.e.,
Fy =
Ky(t)
by (t).
~Fm,r
dM
d~FM,r
d~Fs
d~Fd
x
i
According to Lagranges equations
d
T
T
= Fy .
dt y (t)
y(t)
After simplification, the left-hand side equals
My (t)
mrw 2 cos wt =
Ky(t)
by (t)
191
cos q (t) j
MR2 q 2 (t)
2p
since dm = M dj/2p.
Z 2p
0
sin j
1) dj = MR2 q 2 (t).
y
j
i
x
R
q (t)
q (t)
In the free-body diagram below, d~Fr denotes an infinitesimal reaction force that ensures that the ring retains its shape and mass
distribution and is restricted to rotary motion about the horizontal
axis. An instantaneous infinitesimal change dq to the generalized
coordinate q (t) results in the instantaneous infinitesimal changes
d~r = R
sin j + q (t) + cos q (t) dq i
ring
=
=
d~Fg d~r
Z 2p
Mg
R cos j + q (t) + sin q (t) dj dq
2p
0
ring
d~Ft d~r = k
ring
~r d~Ft dq = t (t) dq
192
The generalized force Fq corresponding to the generalized coordinate q (t) is the coefficient of dq in the expression for the total
infinitesimal work, i.e.,
MgR sin q (t) + t (t).
Fq =
y
j
i
x
d~Ft
~r
k
dm
d~Fr
d~Fg
According to Lagranges equations
d
T
T
= Fq ,
dt q (t)
q (t)
i.e.
2MR2 q(t) =
R cos q (t) j
m1 + m2 2
m R2
x (t) + 1 q 2 (t) + m1 R cos q (t) x (t)q (t).
2
2
j
~r2
m2
k
R
q (t)
m1
193
~Fr,2
j
Kx (t)i
m2
k
m2 g j
~Fr,1
m1
m1 g j
Kx (t) dx
Kx (t), Fq =
m1 gR sin q (t).
d
T
T
= Fx ,
dt x (t)
x (t)
d
T
T
= Fq .
dt q (t)
q (t)
But,
T
x (t)
T
q (t)
T
x (t)
T
q (t)
194
(m1 + m2 ) x (t)
and
m1 gR sin q (t).
Kx (t)
195
Lesson 10 - Linearization
TAKE-AWAYS:
* Systems that violate the principle of superposition are nonlinear. Such systems may be linearized about an operating
point by substituting first-order Taylor expansions for of all
nonlinear terms. The linearization provides an approximate
description of the dynamics for as long as the system state is
close to the operating point.
* Equilibria correspond to constant signals for all system variables. Changes in the number of stability of equilibria are bifurcations.
* The linearization about an equilibrium may predict the
stability of the equilibrium. In particular, an equilibrium is
asymptotically stable if small perturbations decay to zero as
t ! , and unstable if some perturbations (no matter how
small) grow without bound as t ! .
* For lightly damped mechanical systems with more than one
degree of freedom, the linearization about asymptotically stable equilibria can be understood in terms of modal oscillations. These describe transient small oscillations about the
equilibrium in the nonlinear problem, as well as resonance
to harmonic excitation.
196
Let q (t) denote the angle between the ray from the center of the circle
to the bead and the vertical axis, as illustrated below, and define q (t) =
dq/dt and q(t) = d2 q/dt2 . Using either Newtons 2nd law or Lagranges
equations, we obtain the equation of motion
mRq(t)
mg sin q (t)
(533)
m
q (t)
x
bq (t).
m
y
y
wt
x
This equation is nonlinear due to the presence of the trigonometric functions sin q (t) and cos q (t). It is not true that aq (t) is a solution if q (t) is a
solution, in violation of the principle of superposition.
197
We linearize a nonlinear differential equation about an operating point by replacing every nonlinear term by a first-order Taylor
expansion in small deviations from the operating point.
Example 65. To linearize the nonlinear differential equation
mRq(t)
bq (t)
mg sin q (t)
(534)
about the operating point 0, assume that q (t) := #d(t) for some function
d(t) and 0 < # 1. It follows that
q (t) = #d(t), q(t) = #d(t),
(535)
bd(t)
(537)
(538)
Rw 2 )d(t) = 0,
(539)
since # 6= 0.
Analysis of the linearized problem gives insight into the original
nonlinear problem, as long as the linearization remains valid.
(541)
The free response of this dynamic system equals a combination of the unit
impulse response h(t) corresponding to the transfer function
H (s) =
mRs2
1
+ bs + m( g
Rw 2 )
g1 #
g2 #
h(t) = e e
( t ),
(542)
(543)
198
where
g1,2 =
b2
4m2 R( g
2mR
Rw 2 )
(544)
m
R
mg sin q (t)
bq (t)
(545)
(546)
199
3p
4
p
2
p
4
arccos( g/Rw 2 )
w2
g/R
p
4
p
2
To linearize the governing differential equation about an equilibrium operating point q , let q (t) = q + #d(t). It follows that
q (t) = #d(t), q(t) = #d(t)
(547)
sin q #d(t)
(548)
(549)
(550)
and
The linearized differential equation becomes
= mg sin q + #
mg cos q d(t)
bd(t)
(551)
200
201
Rw 2 cos 2q )d(t) = 0.
(552)
Rw 2 cos 2q =
g2
Rw 2 2 2 4
R w
g2
Rw 2
g2
R2 w 4
(553)
(554)
which is always positive when w 2 > g/R. Thus, when it exists, the
corresponding equilibrium configuration is asymptotically stable. In
the figure below, we modify the bifurcation diagram to use dashed
curve to represent unstable equilibria.
Figure 101: Solid curves represent stable equilibria. Dashed curves represent
unstable equilibria.
3p
4
p
2
p
4
arccos( g/Rw 2 )
w2
g/R
p
4
p
2
shown below. Suppose that motion takes place under the influence
of a uniform gravitational field with acceleration g, a spring of undeformed length ` and stiffness K connecting the second bead to a fixed
location along the axis, and a viscous damping force acting on the
second bead.
K, `
m2
m1
Let x (t) denote displacements of the second bead from the undeformed configuration of the spring, and let q (t) denote the angular
displacement of the first bead from a position immediately beneath
the second bead. Suppose that the units of mass, length, and time are
chosen so that m1 = 1, m2 = 1/, g = , and R = 1. The equations of
motion become
(1 + 1/) x (t)
and
Kx (t)
b x (t)
sin q (t),
(555)
(556)
(557)
It follows that
x (t) = #dx (t), x (t) = #dx (t), q (t) = #dq (t), q(t) = #dq (t).
(558)
sin q #d(t) = 1,
(559)
(560)
(561)
202
203
where the upper and lower signs refer to q = 0 and q = p, respectively. The linearized differential equations become
bdx (t)
(562)
(563)
(564)
(565)
dx (0)
dx (0)
1
s
1
s
D x (s)
+
dq (0)
1
s
+
+
1
s
Dq (s)
dq (0)
1+
Here, D x (s) := L dx (#) (s) and Dq (s) := L dq (#) (s) satisfy the
matrix equation
!
D x (s)
A(s)
= B ( s ),
(566)
Dq (s)
where
A(s) =
and
B(s) =
s2 + bs + K
(bs + K )
s2 (1 + )
(567)
(568)
It follows that
D x (s)
Dq (s)
=A
( s ) B ( s ),
(569)
where
A
(s) =
(570)
| A(s)| = s4 + bs3 + K (1 + ) s2 bs K.
(571)
The stability of either equilibrium configuration is clearly determined by the unit impulse response h(t) corresponding to the
transfer function
1
H (s) =
(572)
| A(s)|
since dx (t) and dq (t) are linear combinations of h(t) and its distributional derivatives. These decay exponentially as t ! provided
that all poles of H (s) have negative real part, and grow unbounded
exponentially fast if any pole has positive real part.
Exercises
1. Linearize the function f ( x ) = x/(1 + x2 ) about the operating point
x = 3.
p
2. Linearize the function g( x, q ) = 1 + x2 (q sin q ) about the
operating point ( x , q ) = (1, p/4).
3. The bead in the figure below slides along the straight line y =
kx + .
y
y = kx +
m
K, `
x
K (k~r k
`)
(1 + k 2 ) x (t) + k
,
k~r k
p
where k~r k = x2 (t) + (kx (t) + )2 . Find all equilibrium configurations, and draw a bifurcation diagram with ` along the horizontal axis for the case that = 1, k = 0.
204
R
g
mg sin q (t).
Find all physically realizable equilibrium configurations and describe the corresponding linearized dynamics.
6. Show that one of the roots of the polynomial
s2 + bs + K
bs + K
s2
s2
+1+
Solutions
1. Let x = 3 + #d for # 1. By Taylors theorem, to first order in #,
f (3 + #d) f (3) + # f 0 (3)d
205
1 x2
(1 + x 2 )2
3
10
2
(x
25
3)
g
g
g(1 + #dx , p/4 + #dq ) g(1, p/4) + #
(1, p/4)dx + (1, p/4)dq
x
q
But,
p
g(1, p/4) = p
2 2
and
g
x
( x, q ) = p
(q
x
1 + x2
p
g
( x, q ) = 1 + x2 (1
q
sin q ),
It follows that
p
g( x, q ) p
2 2
1+
p
p
4 2
1
(x
2
1) +
cos q )
p
4
x 2 + (kx + )2
if either (1 + k 2 ) x + k = 0 or
case,
x =
(1 + k 2 ) x + k
p
=0
x 2 + (kx + )2
206
q
2
1
`2
`2
1
2
4. Let
x (t) =
so that
k
+ #d(t)
1 + k2
2
1 + k2
1 + k2
md(t) + K 1
` d(t) = 0
2
p
As long as ` < ||/ 1 + k 2 , all solutions to this equation are
harmonic signals with angular frequency
v
s
u
uK
1 + k2
t
1
`
.
m
2
p
In contrast, for ` > ||/ 1 + k 2 , |d(t)| grows exponentially fast for
almost all initial conditions.
5. An equilibrium configuration corresponds to q (t) = q , q (t) =
q(t) = 0 for some constant q . Here, sin q = 0 implies that
q = 0 or p. A moments reflection shows that only the former is
a physically realizable equilibrium in which the string is taut and
able to provide a force that balances that of gravity.
To linearize the equation of motion about the equilibrium configuration, let q (t) = #d(t), such that
q (t) = #d(t), q(t) = #d(t).
207
g/r0 .
s2
i.e., s = jw or s = g, where
1
w2 =
and
+ K +
( 1
2
K = 0
+ K )2 + 4K
(1 + K )2 + 4K
2
are both positive. It follows that there exists one real positive root
when b = 0. As b increases, this can become negative only if
it crosses 0, but 0 is never a root for any b. It follows that one
root has positive real part when b > 0 and the corresponding
equilibrium configuration is unstable.
2
g =
1+
K +
7. Assume that no two roots coincide. Then, if s(b) denotes the dependence of a root on b, then s(0) satisfies the equation
s(0)4 + K + 1 + s(0)2 + K = 0
i.e., s(0) = jw, where
2
w =
K + 1 +
(K + 1 + )2
2
4K
1 + s2 (0)
2 1 + + K + 2s2 (0)
1 + + K
s 0 (0) =
1+ p
,
4
(1 + + K )2 4K
which is always negative, since
(1 + + K )2
4K
( 1 + + K )2 = 4
It follows that all roots have negative real part when b > 0 and the
corresponding equilibrium configuration is asymptotically stable.
208
1
s2 + 2
s2 + 1
1
T
For the corresponding eigenvectors s t
,
p
1 5
s=
t.
2
209