20 6 Transfer Functions PDF
20 6 Transfer Functions PDF
Introduction
In this Section we introduce the concept of a transfer function and then use this to obtain a Laplace
transform model of a linear engineering system. (A linear engineering system is one modelled by a
constant coefficient ordinary differential equation.)
We shall also see how to obtain the impulse response of a linear system and hence to construct the
general response by use of the convolution theorem.
' $
• be able to use the convolution theorem
HELM (2008): 63
Section 20.6: Transfer Functions
1. Transfer functions and linear systems
Linear engineering systems are those that can be modelled by linear differential equations. We
shall only consider those sytems that can be modelled by constant coefficient ordinary differential
equations.
Consider a system modelled by the second order differential equation.
d2 y dy
a 2
+ b + cy = f (t)
dt dt
in which a, b, c are given constants and f (t) is a given function. In this context f (t) is often called
the input signal or forcing function and the solution y(t) is often called the output signal.
We shall assume that the initial conditions are zero (in this case y(0) = 0, y 0 (0) = 0).
Now, taking the Laplace transform of the differential equation, gives:
(as2 + bs + c)Y (s) = F (s)
in which we have used y(0) = y 0 (0) = 0 and where we have designated L{y(t)} = Y (s) and
L{f (t)} = F (s).
We define the transfer function of a system to be the ratio of the Laplace transform of the output
signal to the Laplace transform of the input signal with the initial conditions as zero. The transfer
function (a function of s), is denoted by H(s). In this case
Y (s) 1
H(s) ≡ = 2
F (s) as + bs + c
Now, in the special case in which the input signal is the delta function, f (t) = δ(t), we have F (s) = 1
and so,
H(s) = Y (s)
We call the solution to the differential equation in this special case the unit impulse response
function and denote it by h(t)u(t) (we include the step function u(t) to emphasize its causality).
So
h(t)u(t) = L−1 {H(s)} when f (t) = δ(t)
Now, keeping this in mind and returning to the general case in which the input signal f (t) is not
necessarily the impulse function δ(t), we have:
Y (s) = H(s)F (s)
and so the solution for the output signal is, as usual, obtained by taking the inverse Laplace transform:
64 HELM (2008):
Workbook 20: The Laplace Transform
®
Key Point 17
Linear System Solution
The solution to a linear system, modelled by a constant coefficient ordinary differential equation,
is given by the convolution of the unit impulse response function h(t)u(t) with the input function
f (t).
This approach provides yet another method of solving a linear system as Example 6 illustrates.
Example 6
Find the impulse response function h(t) to a linear engineering system
modelled by the differential equation:
d2 y
2
+ 4y = e−t y(0) = 0 y 0 (0) = 0
dt
and hence solve the system.
Solution
Here
1 1
H(s) = = with a = 1, b = 0, c = 4
s2 +4 as2 + bs + c
This is obtained by replacing the forcing function e−t by the impulse function δ(t) and then taking
the Laplace transform. Using this:
1
h(t) = L−1 {H(s)} = L−1 { } = 12 sin 2t.u(t)
s2 +4
Then the output y(t) corresponding to the input e−t is given by the convolution of e−t and h(t).
That is,
Z t
−t
y(t) = (h ∗ e )(t) = 1
2
sin 2(t − x)e−x dx
0
1
sin 2t − 2 cos 2t + 2e−t
= 10
(Note: the last integral can be determined by integrating by parts (twice), or by use of a computer
algebra system such as Matlab.)
HELM (2008): 65
Section 20.6: Transfer Functions
Task
Use the transfer function approach to solve
dx
− 4x = sin t x(0) = 0.
dt
Your solution
Answer
You should find H(s) = 1/(s − 4) since the transfer function is the Laplace transform of the output
X(s) when the input is a delta function δ(t).
Now obtain an expression for the solution x(t) in terms of the convolution:
Your solution
Answer
You should obtain x(t) = (sin t ∗ h)(t) where
Z t
−1 −1 1
h(t) = L {H(s)} = L = e4t u(t) and x(t) = (sin x)e4(t−x) u(t − x) dx
s−4 0
Your solution
Answer
If t > 0 then u(t − x) = 1 and so
( t Z t )
Z t
sin x cos x
x(t) = sin x e4(t−x) dx = e4t − e−4x − − e−4x dx
0 4 0 0 4
Z t
4t sin t −4t 1 h cos x −4x it sin x −4x
= e − e + − e − e dt
4 4 4 0 0 4
1 4t
Therefore x(t) = − 41 sin t − 1
16
cos t + 16
e − 1
16
x(t)
1
−4 sin t − cos t + e4t
Hence x(t) =
17
66 HELM (2008):
Workbook 20: The Laplace Transform
®
HELM (2008): 67
Section 20.6: Transfer Functions
There are three basic components occurring in block diagrams which we now describe.
The first component, we have already met: the block relating the input to the output (Figure
23):
F(s) Y(s)
H(s)
—
X(s)
Here we have shown two incoming signals R(s), X(s) (but at a general summing point there may
be many incoming signals) and one outgoing signal (there should never be more than one outgoing
signal). The sign attached to the incoming signal defines whether the signal is adding to (+) or
subtracting from (−) the summing point. The outgoing signal is then calculated in an obvious way,
taking these signs into account.
The third component is a take-off point (Figure 25):
Y(s) Y(s)
Y (s)
Here the value of the signal Y (s) is found in such a way as not to affect the signal that is being
transmitted. (This situation can never be precisely realised in practice, but using sensitive measuring
devices it can be well approximated. As a simple example consider the problem of measuring the
temperature of a certain volume of liquid. The act of putting a thermometer in the liquid will usually
slightly affect the temperature we are trying to observe.)
68 HELM (2008):
Workbook 20: The Laplace Transform
®
An example of a block diagram is the so-called negative feedback loop, shown in Figure 26 (we are
using G(s) to denote the transfer function):
—
Y(s)
F(s) Y(s)
G(s)
1 + G(s)
dx
+ x = f (t)
dt
dy
3 − y = x(t)
dt
x(0) = 0 y(0) = 0
HELM (2008): 69
Section 20.6: Transfer Functions
In terms of Laplace transforms we have, as usual
sX(s) + X(s) = F (s) 3sY (s) − Y (s) = X(s)
so the transfer function for the first equation (G1 (s) say) satisfies
X(s) 1
G1 (s) ≡ =
F (s) s+1
whilst the transfer function for the second equation G2 (s) satisfies
Y (s) 1
G2 (s) ≡ =
X(s) 3s − 1
In pictorial terms this is shown in Figure 28:
Figure 28
So we have two transfer functions ‘in series’. To find how they combine we simply find an expression
connecting the final output Y (s) to the initial input F (s). Clearly
X(s) = G1 (s)F (s) and so Y (s) = G2 (s)X(s) = [G2 (s)G1 (s)] F (s)
So transfer functions in series are simply multiplied together. In this case the overall transfer function
H(s) is:
1
H(s) = G1 (s)G2 (s) =
(s + 1)(3s − 1)
Note that this result could be found directly from the differential equations used to model this system.
If we differentiate the second differential equation of the original pair we get:
d2 y dy dx
3 2 − =
dt dt dt
dx
Rearranging the first equation gives = f (t) − x
dt
d2 y dy
dy
Substituting gives: 3 2 − = f (t) − x = f (t) − 3 − y
dt dt dt
or
d2 y dy
3 2
+ 2 − y = f (t)
dt dt
which, on taking Laplace transforms, gives the s-relation (3s2 + 2s − 1)Y (s) = F (s) implying a
transfer function:
1 1
H(s) = 2 =
3s + 2s − 1 (s + 1)(3s − 1)
as obtained above.
70 HELM (2008):
Workbook 20: The Laplace Transform
®
Engineering Example 3
System response
Figure 29
Determine the system response v0 (t) when the input function is a unit step function when K = 2.5
and a = 0.
Solution
If the system has an overall transfer function H(s) then V0 (s) = H(s)V1 (s). But this particular
system is the negative feedback loop described earlier and so
K
G(s) K
H(s) = = 1 + as =
1 + G(s) K K + 1 + as
1+
1 + as
In this particular case
2.5 5
H(s) = =
3.5 + 0.5s 7+s
Thus the impulse response h(t) is
−1 −1 5
h(t) = L {H(s)} = L = 5e−7t u(t)
(7 + s)
and so the response to a step input u(t) is given by the convolution of h(t) with u(t)
Z t
v0 (t) = u(t − x)5e−7x u(t) dx
Z0 t
= 5e−7x dx t>0
0
t
5 −7x 5
= − e = − [e−7t − 1]
7 0 7
HELM (2008): 71
Section 20.6: Transfer Functions
Table of Laplace Transforms
Rule Causal function Laplace transform
1 f (t) F (s)
1
2 u(t)
s
n!
3 tn u(t)
sn+1
1
4 e−at u(t)
s+a
a
5 sin at . u(t)
s + a2
2
s
6 cos at . u(t)
s + a2
2
b
7 e−at sin bt . u(t)
(s + a)2 + b2
s+a
8 e−at cos bt u(t)
(s + a)2 + b2
72 HELM (2008):
Workbook 20: The Laplace Transform