Lecture 3-Transfer Functions v2.0
Lecture 3-Transfer Functions v2.0
Contents
Objective ........................................................................................................................................... 2
What is a model? Why do we need one? ............................................................................. 2
Linear Time Invariant (LTI) systems ................................................................................... 4
TIME INVARIANT.......................................................................................................................... 4
LINEARITY ....................................................................................................................................... 5
Recap on Laplace Transforms ................................................................................................. 6
Transfer function of a first order ODE ................................................................................. 6
Capacitor........................................................................................................................................ 11
First Order Resistor Capacitor (RC) network ................................................................ 11
Transfer function for a second order ODE ...................................................................... 15
Second order RLC circuit ........................................................................................................ 17
Higher Order Systems .............................................................................................................. 19
Poles and zeros ........................................................................................................................... 19
Remarks on behaviour ............................................................................................................ 20
Summary........................................................................................................................................ 21
Exercises ........................................................................... Error! Bookmark not defined.
_____________________________________________________________________________________
Page 1 of 21
Objective
_____________________________________________________________________________________
Page 2 of 21
2. Unmodeled dynamics (to make simpler model).
There is always a tradeoff between simplicity and accuracy. It’s often possible
to improve the accuracy of a mathematical model by increasing its complexity.
Simplification often means ignoring some inherent physical properties, e.g.,
ignore nonlinearities in linear, lumped-parameter models.
In general, it’s desirable to start with a simplified model to get a “general feel,”
increasing complexity only if the controlled system does not meet
performance requirements.
KEY POINT: “All models are wrong, but some are useful”
TIME INVARIANT
The first system property that we look at is that of time invariance.
_____________________________________________________________________________________
Page 4 of 21
We can test a system for this property using ideas from the figure.
LINEARITY
The second property that we look at is linearity.
𝑥3 (𝑡) = 𝛼𝑥1 (𝑡) + 𝛽𝑥2 (𝑡) → 𝑦3 (𝑡) = 𝛼𝑦1 (𝑡) + 𝛽𝑦2 (𝑡) ,
_____________________________________________________________________________________
Page 5 of 21
Input 𝑥3 (𝑡) = 𝛼𝑥1 (𝑡) + 𝛽𝑥2 (𝑡) to the system and measure 𝑦3 (𝑡)
If 𝑦3 (𝑡) = 𝛼𝑦1 (𝑡) + 𝛽𝑦2 (𝑡)for all possible α and β values, and 𝑥1 (𝑡) and
𝑥2 (𝑡), then the system is linear.
𝑑𝑦(𝑡)
𝑇 + 𝑦(𝑡) = 𝐾𝑢(𝑡)
𝑑𝑡
Taking Laplace transforms (ignore initial conditions) gives:
Rearranging gives:
K
y ( s) u ( s) y ( s) G(s)u (s)
Ts 1
Where K
G( s)
Ts 1
_____________________________________________________________________________________
Page 6 of 21
From the above expression, the Laplace transform of the output signal 𝑦̅(𝑠) is
given by the transfer function representation of the system 𝐺(𝑠) multiplied
by the Laplace transform of the input signal 𝑢
̅ (𝑠). The transfer function 𝐺(𝑠)
provides us with a simple representation of process dynamics and provides a
direct relationship between inputs and outputs which is needed for control
design and analysis.
FORCE BALANCE: The normal way to model these scenarios is to consider the
force as being made up of different (parallel) components, here
𝑓1 (𝑡) and 𝑓2 (𝑡).
_____________________________________________________________________________________
Page 7 of 21
𝑓1 (𝑡) represents the friction force/drag. This is assumed to be proportional
to speed.
1/𝑓𝑣 𝐾
𝐺(𝑠) = ≡
𝑀
𝑠 + 1 𝑇𝑠 + 1
𝑓𝑣
_____________________________________________________________________________________
Page 8 of 21
Where
1 𝑀
𝐾= and 𝑇 =
𝑓𝑣 𝑓𝑣
1/𝑓𝑣
𝑀
𝑠+1
𝑓𝑣
𝑓1 (𝑡) = 𝑘𝑥(𝑡);
𝑑𝑥(𝑡)
𝑓2 (𝑡) = 𝑓𝑣 ;
𝑑𝑡
𝑓(𝑡) = 𝑓1 (𝑡) + 𝑓2 (𝑡)
The simple model for a spring-damper system is a 1st order ODE with
parameters the friction coefficient 𝑓𝑣 and the spring stiffness 𝑘.
_____________________________________________________________________________________
Page 9 of 21
𝑑𝑥(𝑡)
𝑓(𝑡) = 𝑓1 (𝑡) + 𝑓2 (𝑡) = 𝑘𝑥(𝑡) + 𝑓𝑣
𝑑𝑡
1/𝑘
𝑓𝑣
𝑠+1
𝑘
1st order models of simple electrical networks tend to arise for series
arrangements of components and thus that is the focus here. For series
arrangements we use Kirchhoff’s voltage laws (KVL).
_____________________________________________________________________________________
Page 10 of 21
Capacitor
A capacitor is a storage device. It stores charge. When current (charge/sec)
flows through a capacitor, this charge 𝑞(𝑡) is stored or released from storage
according to the direction of flow. The voltage 𝑉𝐶 (𝑡)across a capacitor is
linked to the total charge 𝑞(𝑡) stored.
2. Application of KVL.
𝑑𝑞(𝑡) 1
𝑉(𝑡) = 𝑉𝑟 (𝑡) + 𝑉𝐶 (𝑡) = 𝑅 + 𝑞(𝑡)
𝑑𝑡 𝐶
_____________________________________________________________________________________
Page 11 of 21
Taking Laplace transforms assuming zero initial conditions
1
𝑉̅ (𝑠) = 𝑅𝑠𝑞̅(𝑠) + 𝑞̅(𝑠)
𝐶
1
𝑞̅(𝑠) = 𝑉̅ (𝑠) ⇒ 𝑞̅(𝑠) = 𝐺(𝑠)𝑉̅ (𝑠)
1
𝑅𝑠 +
𝐶
1
𝐺(𝑠) =
1
𝑅𝑠 +
𝐶
Multiply the numerator and denominator by 𝐶 yields
𝐶 𝐾
𝐺(𝑠) = ≡
𝑅𝐶𝑠 + 1 𝑇𝑠 + 1
Where
𝐾 = 𝐶 and 𝑇 = 𝑅𝐶
𝐶
𝑅𝐶𝑠 + 1
_____________________________________________________________________________________
Page 12 of 21
1
𝑉̅𝐶 (𝑠) = 𝑉̅ (𝑠) ⇒ 𝑉̅ (𝑠) = 𝐺(𝑠)𝑉̅𝐶 (𝑠)
𝑅𝐶𝑠 + 1
Therefore
1 𝐾
𝐺(𝑠) = ≡
𝑅𝐶𝑠 + 1 𝑇𝑠 + 1
Where
𝐾 = 1 and 𝑇 = 𝑅𝐶
1
𝑅𝐶𝑠 + 1
A tank will store fluid, thus it is worth considering how much it can store and
how this storage is linked to flow rate into and out of the tank. Assume
constant cross-sectional area A, stored volume V, tank height h, flow
resistance R, fluid density 𝝆, gravitational force g=9.81m/s2 , flow in Fin and
flow out of Fout. Assume that flow is created by differential pressure between
the two ends of the pipe.
P1 P2 Fout R
_____________________________________________________________________________________
Page 13 of 21
V Ah (t )
Fout (t ) gh(t ) / R
dV dh(t )
A Fin (t ) Fout (t )
dt dt
P1 P 2 Fout (t ) R
P1 P 2 gh(t )
𝑃1 − 𝑃2 𝜌𝑔ℎ(𝑡)
𝐹𝑜𝑢𝑡 (𝑡) = =
𝑅 𝑅
𝑑ℎ(𝑡)
𝐴 = 𝐹𝑖𝑛 (𝑡) − 𝐹𝑜𝑢𝑡 (𝑡)
𝑑𝑡
𝑑ℎ(𝑡) 𝜌𝑔ℎ(𝑡)
𝐴 + = 𝐹𝑖𝑛 (𝑡)
𝑑𝑡 𝑅
Taking Laplace transforms assuming zero initial conditions
𝜌𝑔
𝐴𝑠ℎ̅(𝑠) + ℎ̅(𝑠) = 𝐹̅𝑖𝑛 (𝑠)
𝑅
1
ℎ̅(𝑠) = ̅
𝜌𝑔 𝐹𝑖𝑛 (𝑠)
𝐴𝑠 +
𝑅
𝑅
Multiply the numerator and denominator by yields
𝜌𝑔
𝑅/𝜌𝑔
ℎ̅(𝑠) = 𝐹̅𝑖𝑛 (𝑠) ⇒ ℎ̅(𝑠) = 𝐺(𝑠)𝐹̅𝑖𝑛 (𝑠)
𝐴𝑅
𝑠+1
𝜌𝑔
Where
𝑅/𝜌𝑔 𝐾
𝐺(𝑠) = ≡
𝐴𝑅
𝑠 + 1 𝑇𝑠 + 1
𝜌𝑔
And
_____________________________________________________________________________________
Page 14 of 21
𝐴𝑅 𝑅
𝑇= and 𝐾 =
𝜌𝑔 𝜌𝑔
𝑅/𝜌𝑔
𝐴𝑅
𝑠+1
𝜌𝑔
𝑑2 𝑑
𝑎 2 𝑦(𝑡) + 𝑏 𝑦(𝑡) + 𝑐𝑦(𝑡) = 𝐾𝑢(𝑡)
𝑑𝑡 𝑑𝑡
Assume all initial conditions are zero, then taking Laplace transforms gives,
or
and so,
𝛼
𝑦̅(𝑠) = 𝑢̅(𝑠)
𝑎𝑠 2 + 𝑏𝑠 + 𝑐
𝑦̅(𝑠) = 𝐺(𝑠)𝑢̅(𝑠)
The function
𝛼
𝐺(𝑠) =
𝑎𝑠 2 + 𝑏𝑠 + 𝑐
is called the transfer function from 𝑢̅(𝑠) (the input) to 𝑦̅(𝑠) (the output). Clearly
the same technique will work for higher order linear ordinary differential
equations with constant coefficients. For such systems, the transfer function
can be regarded as a placeholder for the coefficients of the differential equation.
𝑑𝑣(𝑡) 𝑑 2 𝑥(𝑡)
∑ 𝐹 = 𝑀𝑎 = 𝑀 =𝑀
𝑑𝑡 𝑑𝑡 2
The spring force is proportional to the displacement of the mass, 𝑥(𝑡), and the
viscous damping force is proportional to the velocity of the mass, 𝑀 . Both
forces oppose the motion of the mass and are, therefore in the negative 𝑥 -
direction. Note also that 𝑥(𝑡) = 0, corresponds to the position of the mass
when the spring is unstretched.
and these are all drawn together using Newton’s second law of motion:
𝑑 2 𝑥(𝑡) 𝑑𝑥(𝑡)
𝑀 + 𝑓𝑣 + 𝐾𝑥(𝑡) = 𝑓(𝑡)
𝑑𝑡 2 𝑑𝑡
_____________________________________________________________________________________
Page 16 of 21
Taking Laplace transforms assuming zero initial conditions
𝑥̅ (𝑠)(𝑀𝑠 2 + 𝑓𝑣 𝑠 + 𝑘) = 𝑓(̅ 𝑠)
1
𝑥̅ (𝑠) = 𝑓 (̅ 𝑠) ⇒ 𝑥̅ (𝑠) = 𝐺(𝑠)𝑓(̅ 𝑠)
𝑀𝑠 2 + 𝑓𝑣 𝑠 + 𝑘
1
𝐺(𝑠) =
𝑀𝑠 2 + 𝑓𝑣 𝑠 + 𝑘
1
𝑀𝑠 2 + 𝑓𝑣 𝑠 + 𝑘
In any problem, the designer must first decide what the input and output
should be. In this network, several variables could have been chosen to be the
output. For example, the inductor voltage, the capacitor voltage, the resistor
voltage, or the current. The problem statement, however, is clear in this case:
We are to treat the capacitor voltage as the output and the applied voltage as
the input.
_____________________________________________________________________________________
Page 17 of 21
By writing down the voltage quantities across the inductor, capacitor and
resistor
1 1
𝑉𝑅 (𝑡) = 𝑖(𝑡)𝑅 𝑉𝐶 (𝑡) = ∫ 𝑖(𝜏)𝑑𝜏
𝐶 0
𝑑𝑖(𝑡)
𝑉𝐿 (𝑡) = 𝐿
𝑑𝑡
Using Kirchhoff’s voltage law
𝑑𝑖(𝑡) 1 1
𝑉(𝑡) = 𝐿 + 𝑖(𝑡)𝑅 + ∫ 𝑖(𝜏)𝑑𝜏
𝑑𝑡 𝐶 0
𝑑𝑞(𝑡)
Changing variables from current to charge using 𝑖(𝑡) =
𝑑𝑡
𝑑 2 𝑞(𝑡) 𝑑𝑞(𝑡) 1
𝑉(𝑡) = 𝐿 + 𝑅 + 𝑞(𝑡)
𝑑𝑡 2 𝑑𝑡 𝐶
From the voltage-charge relationship for a capacitor in
1/𝐿𝐶
𝑅 1
𝑠2 + 𝑠 +
𝐿 𝐿𝐶
Clearly the same technique will work for higher order linear ordinary
differential equations with constant coefficients.
𝑑𝑛 𝑥 𝑑 𝑛−1 𝑥 𝑑𝑥 𝑑 𝑛−1 𝑢 𝑑𝑢
𝑎𝑛 𝑛 + 𝑎𝑛−1 𝑛−1 + ⋯ + 𝑎1 + 𝑎𝑜 𝑥 = 𝑏𝑛−1 𝑛−1 + ⋯ + 𝑏1 + 𝑏𝑜 𝑢
𝑑𝑡 𝑑𝑡 𝑑𝑡 𝑑𝑡 𝑑𝑡
Taking Laplace transforms assuming zero initial conditions
4𝑠 2 − 8𝑠 − 60 4𝑠 2 − 8𝑠 − 60
𝐺(𝑠) = 3 =
𝑠 + 2𝑠 2 + 2𝑠 𝑠(𝑠 2 + 2𝑠 + 2)
Imag(s)
Real(s)
Remarks on behaviour
The transfer function can be generalised as
𝑥̅ (𝑠) = 𝐺(𝑠)𝑢̅(𝑠)
_____________________________________________________________________________________
Page 20 of 21
Let 𝑢(𝑡) have a particular behaviour, then from Laplace one can imply 𝑥(𝑡)
includes this behaviour. As will be seen, this is invaluable where 𝑢̅(𝑠) is not
known explicitly, such as in a feedback loop.
1 𝐺(𝑠)
𝑢(𝑡) = 1 𝑢̅(𝑠) = 𝑥(𝑡) = ℒ −1 { }
𝑠 𝑠
1 𝐺(𝑠)
𝑢(𝑡) = 𝑡 𝑢̅(𝑠) = 2 𝑥(𝑡) = ℒ −1 { 2 }
𝑠 𝑠
𝜔 𝐺(𝑠)𝜔
𝑢(𝑡) = sin(𝜔𝑡) 𝑢̅(𝑠) = 2 𝑥(𝑡) = ℒ −1 { 2 }
𝑠 + 𝜔2 𝑠 + 𝜔2
Summary
Transfer functions facilitates
_____________________________________________________________________________________
Page 21 of 21