1

C
h
a
p
t
e
r

1
3
Frequency Response Analysis
Sinusoidal Forcing of a First-Order Process
For a first-order transfer function with gain K and time constant ,
the response to a general sinusoidal input, is:

( )
sin = x t A t
( )
( )
/
2 2
cos sin (5-25)
1

= +
+
t
KA
y t e t t
Note that y(t) and x(t) are in deviation form. The long-time
response, y
l
(t), can be written as:
( ) ( )
2 2
sin for (13-1)
1
= +
+
A
KA
y t t t
where:
( )
1
tan

=
2
C
h
a
p
t
e
r

1
3
Figure 13.1 Attenuation and time shift between input and output
sine waves (K= 1). The phase angle of the output signal is given
by , where is the (period) shift and P
is the period of oscillation.

Time shift / 360 =

For ( ) sin , sin as where : = = +

and tan
1

= =
+
KA
A
1. The output signal is a sinusoid that has the same frequency, ,
as the [Link], x(t) =Asint.
2. The amplitude of the output signal, , is a function of the
frequency and the input amplitude, A:

A
2 2

(13-2)
1
=
+
KA
A
3. The output has a phase shift, , relative to the input. The
amount of phase shift depends on .
4
C
h
a
p
t
e
r

1
3
Dividing both sides of (13-2) by the input signal amplitude A
yields the amplitude ratio (AR)
2 2

AR (13-3a)
1
= =
+
A K
A
which can, in turn, be divided by the process gain to yield the
normalized amplitude ratio (AR
N
)
N
2 2
1
AR (13-3b)
1
=
+
5
C
h
a
p
t
e
r

1
3
Shortcut Method for Finding
the Frequency Response
The shortcut method consists of the following steps:
Step 1. Set s=j in G(s) to obtain .
Step 2. Rationalize G(j); We want to express it in the form.
G(j)=R + jI
where R and I are functions of . Simplify G(j) by
multiplying the numerator and denominator by the
complex conjugate of the denominator.
Step 3. The amplitude ratio and phase angle of G(s) are given
by:
( )
G j
2 2
1
AR
tan ( / )
R I
R I

= +
=
Memorize
6
C
h
a
p
t
e
r

1
3
Example 13.1
Find the frequency response of a first-order system, with
( )
1
(13-16)
1
G s
s
=
+
Solution
First, substitute in the transfer function s j =
( )
1 1
(13-17)
1 1
G j
j j
= =
+ +
Then multiply both numerator and denominator by the complex
conjugate of the denominator, that is, 1 j +
( )
( )( )
( )
2 2
2 2 2 2
1 1

1 1
1

1
(13-18)
1 1
j j
G j
j j
j R jI
+ +
= =
+ +
+

= + = +
+ +
7
C
h
a
p
t
e
r

1
3
2 2
1
(13-19a)
1
R =
+
where:
2 2

(13-19b)
1
I

=
+
From Step 3 of the Shortcut Method,
2 2
2 2
2 2 2 2
1
AR
1 1

| | | |
= + = +
| |
+ + \ . \ .
R I
( )
( )
2 2
2
2 2
2 2
1
1
AR (13-20a)
1
1
+
= =
+
+
Also,
( ) ( )
1 1 1
tan tan tan (13-20b)

| |
= = =
|
\ .
I
R
or
8
C
h
a
p
t
e
r

1
3
Complex Transfer Functions
Consider a complex transfer G(s),
( )
( ) ( ) ( )
( ) ( ) ( )
1 2 3
(13-22) =
"
"
a b c
G s G s G s
G s
G s G s G s
( )
( ) ( ) ( )
( ) ( ) ( )
1 2 3

(13-23)

=
"
"
a b c
G j G j G j
G j
G j G j G j
Substitute s=j,
From complex variable theory, we can express the magnitude and
angle of as follows:
( )
G j
( )
( ) ( ) ( )
( ) ( ) ( )
1 2 3

(13-24a)

=
"
"
a b c
G j G j G j
G j
G j G j G j
( ) ( ) ( ) ( )
( ) ( ) ( )
1 2 3

[ ] (13-24b)
a b c
G j G j G j G j
G j G j G j
= + + +
+ + +
"
"
9
C
h
a
p
t
e
r

1
3
Bode Diagrams
A special graph, called the Bode diagram or Bode plot,
provides a convenient display of the frequency response
characteristics of a transfer function model. It consists of
plots of AR and as a function of .
Ordinarily, is expressed in units of radians/time.

Bode Plot of A First-order System

10
C
h
a
p
t
e
r

1
3
Figure 13.2 Bode diagram for a first-order process.
11
C
h
a
p
t
e
r

1
3
Note that the asymptotes intersect at , known as the
break frequency or corner frequency. Here the value of AR
N
from (13-21) is:
1/
b
= =
( )
N
1
AR 0.707 (13-30)
1 1
b
= = =
+
Some books and software defined AR differently, in terms of
decibels. The amplitude ratio in decibels AR
d
is defined as
d
AR 20 log AR (13-33) =
12
C
h
a
p
t
e
r

1
3
Integrating Elements
The transfer function for an integrating element was given in
Chapter 5:
( )
( )
( )
(5-34)
Y s
K
G s
U s s
= =
( )
AR (13-34)

K K
G j
j
= = =
( ) ( )
90 (13-35) G j K = = =
D
Second-Order Process
A general transfer function that describes any underdamped,
critically damped, or overdamped second-order system is
( )
2 2
(13-40)
2 1
K
G s
s s
=
+ +
13
C
h
a
p
t
e
r

1
3
Substituting and rearranging yields: s j =
( )
( )
2
2
2 2
AR (13-41a)
1 2
K
=
+
1
2 2
2
tan (13-41b)
1

(
=
(

Figure 13.3 Bode diagrams for second-order processes.
14
C
h
a
p
t
e
r

1
3
Time Delay
Its frequency response characteristics can be obtained by
substituting ,
s j =
( )

(13-53)
j
G j e

=
which can be written in rational form by substitution of the
Euler identity,
( )

cos sin (13-54)

= =
j
G j e j
From (13-54)
( )
( )
2 2
1
AR cos sin 1 (13-55)
sin
tan
cos
G j
G j

= = + =
| |
= =
|
\ .
or
(13-56) =
15
C
h
a
p
t
e
r

1
3
Figure 13.6 Bode diagram for a time delay, .
s
e

16
C
h
a
p
t
e
r

1
3
Figure 13.7 Phase angle plots for and for the 1/1 and 2/2
Pad approximations (G
1
is 1/1; G
2
is 2/2).
s
e

17
C
h
a
p
t
e
r

1
3
Process Zeros
Consider a process zero term,
( )
( 1) = + G s K s
Substituting s=j gives
( )
( 1) = + G j K j
Thus:
( )
( ) ( )
2 2
1
AR 1
tan

= = +
= = +
G j K
G j
Note: In general, a multiplicative constant (e.g., K) changes
the AR by a factor of K without affecting .

18
C
h
a
p
t
e
r

1
3
Frequency Response Characteristics of
Feedback Controllers
Proportional Controller. Consider a proportional controller with
positive gain
( )
(13-57)
c c
G s K =
In this case , which is independent of .
Therefore,
( )

c c
G j K =
AR (13-58)
c c
K =
and
0 (13-59)
c
=
D
19
C
h
a
p
t
e
r

1
3
Proportional-Integral Controller. A proportional-integral (PI)
controller has the transfer function (cf. Eq. 8-9),
( )
1 1
1 (13-60)

I
c c c
I I
s
G s K K
s s
| | | |
+
= + =
| |
\ . \ .
Substitute s=j:
( )
1 1 1
1 1

| | | | | |
+
= + = =
| | |

\ . \ . \ .
I
c c c c
I I I
j
G j K K K j
j j
Thus, the amplitude ratio and phase angle are:
( )
( )
( )
2
2
1
1
AR 1 (13-62)

I
c c c c
I
I
G j K K
+
= = + =
( ) ( ) ( )
1 1
tan 1/ tan 90 (13-63)
c c I I
G j

= = =
D
20
C
h
a
p
t
e
r

1
3
Figure 13.9 Bode plot of a PI controller,
( )
10 1
2
10
c
s
G s
s
+
| |
=
|
\ .
21
C
h
a
p
t
e
r

1
3
Ideal Proportional-Derivative Controller. For the ideal
proportional-derivative (PD) controller (cf. Eq. 8-11)
( ) ( )
1 (13-64)
c c D
G s K s = +
The frequency response characteristics are similar to those of a
LHP zero:
( )
2
AR 1 (13-65)
c c D
K = +
( )
1
tan (13-66)
D

=
Proportional-Derivative Controller with Filter. The PD
controller is most often realized by the transfer function
( )
1
(13-67)
1
D
c c
D
s
G s K
s
| |
+
=
|
+
\ .
22
C
h
a
p
t
e
r

1
3
Figure 13.10 Bode
plots of an ideal PD
controller and a PD
controller with
derivative filter.
Idea:
With Derivative
Filter:
G s
( ) ( )
2 4 1
c
s = +
( )
4 1
2
0.4 1
c
s
G s
s
+
| |
=
|
+
\ .
23
C
h
a
p
t
e
r

1
3
PID Controller Forms
Parallel PID Controller. The simplest form in Ch. 8 is
( )
1
1
1

c c D
G s K s
s
| |
= + +
|
\ .
Series PID Controller. The simplest version of the series PID
controller is
( ) ( )
1
1
1
1 (13-73)

c c D
s
G s K s
s
| |
+
= +
|
\ .
Series PID Controller with a Derivative Filter.
( )
1
1
1 1
1
D
c c
D
s s
G s K
s s
| || |
+ +
=
| |
+
\ .\ .
24
C
h
a
p
t
e
r

1
3
Figure 13.11 Bode
plots of ideal parallel
PID controller and
series PID controller
with derivative filter
( = 1).
Idea parallel:
Series with
Derivative Filter:
( )
1
2 1 4
10
c
G s s
s
| |
= + +
|
\ .
( )
10 1 4 1
2
10 0.4 1
c
s s
G s
s s
+ +
| || |
=
| |
+
\ .\ .
25
C
h
a
p
t
e
r

1
3
Nyquist Diagrams
Consider the transfer function
( )
1
(13-76)
2 1
G s
s
=
+
with
( )
( )
2
1
AR (13-77a)
2 1
G j = =
+
and
( ) ( )
1
tan 2 (13-77b) G j

= =
26
C
h
a
p
t
e
r

1
3
Figure 13.12 The Nyquist diagram for G(s) = 1/(2s + 1)
plotting and
( )
( )
Re G j
( )
( )
Im . G j
27
C
h
a
p
t
e
r

1
3
Figure 13.13 The Nyquist diagram for the transfer
function in Example 13.5:
6
5(8 1)
( )
(20 1)(4 1)
s
s e
G s
s s

+
=
+ +