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.

a d 2 y d t 2 + b d y d t + c y = f ( t )

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 ).

Now, taking the Laplace transform of the differential equation, gives:

( a s 2 + b s + c ) Y ( s ) = F ( s )

in which we have used y ( 0 ) = y ( 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

H ( s ) Y ( s ) F ( s ) = 1 a s 2 + b s + 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:

y ( t ) = L 1 { Y ( s ) } = L 1 { H ( s ) F ( s ) } = ( h f ) ( t )

using the convolution theorem.

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:

d 2 y d t 2 + 4 y = e t y ( 0 ) = 0 y ( 0 ) = 0

and hence solve the system.

Solution

Here

H ( s ) = 1 s 2 + 4 ( = 1 a s 2 + b s + c  with a = 1 , b = 0 , c = 4 )

This is obtained by replacing the forcing function e t by the impulse function δ ( t ) and then taking the Laplace transform. Using this:

h ( t ) = L 1 { H ( s ) } = L 1 { 1 s 2 + 4 } = 1 2 sin 2 t . u ( t )

Then the output y ( t ) corresponding to the input e t is given by the convolution of e t and h ( t ) . That is,

y ( t ) = ( h e t ) ( t ) = 0 t 1 2 sin 2 ( t x ) e x d x = 1 10 sin 2 t 2 cos 2 t + 2 e t

(Note: the last integral can be determined by integrating by parts (twice), or by use of a computer algebra system such as Matlab.)

Task!

Use the transfer function approach to solve

d x d t 4 x = sin t x ( 0 ) = 0.

Begin by finding the transfer function H ( s ) :

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:

You should obtain x ( t ) = ( sin t h ) ( t ) where

h ( t ) = L 1 { H ( s ) } = L 1 1 s 4 = e 4 t u ( t )  and x ( t ) = 0 t ( sin x ) e 4 ( t x ) u ( t x ) d x

Now complete the evaluation of this integral:

If t > 0 then u ( t x ) = 1 and so

x ( t ) = 0 t sin x e 4 ( t x ) d x = e 4 t sin x 4 e 4 x 0 t 0 t cos x 4 e 4 x d x = e 4 t sin t 4 e 4 t + 1 4 cos x 4 e 4 x 0 t 0 t sin x 4 e 4 x d t

Therefore x ( t ) = 1 4 sin t 1 16 cos t + 1 16 e 4 t 1 16 x ( t )

Hence x ( t ) = 1 17 4 sin t cos t + e 4 t