GEE338 LCS Chapter - 1 Math Model For LTI System
GEE338 LCS Chapter - 1 Math Model For LTI System
GEE338 LCS Chapter - 1 Math Model For LTI System
Chapter 1
Mathematical Modeling of Linear Time
Invariant (LTI) Systems
Prof. Fawzy Ibrahim and Eng. Medhat Toubar
Electronics and Communication Department
Misr International University (MIU)
1 of 70 GEE338 Ch.1 Math Models Prof Fawzy Ibrahim
Chapter Contents
1.0 Introduction
1.0.1 Control Systems Terminologies
1.0.2 Control systems Classifications
1.0.3 Control systems Applications
1.2 Mathematical Model
1.3 System Representation using Differential Equations and
Transfer Function for:
1.3.1 Electrical and Electronic Elements
1.3.2 Translation Mechanical Elements
1.3.3 Rotational Mechanical Elements
1.3.4 Electro Mechanical Elements
1.4 Block Diagram Models
1.4.1 Components of a block diagram for a LTI system
1.4.2 Block Diagram Reduction Technique
1.5 Signal Flow Graph Models
1.5.1 Definitions 1.5.2 Flow graphs of control systems
1.5.3 Signal Flow Graph Algebra
1.5.4 Signal Flow Graph Reduction Technique
1.5.5 Block Diagram to Signal Flow Graph
Appendix A The Laplace Transform
2 of 70 GEE338 Ch.1 Math Models Prof Fawzy Ibrahim
1.0 Introduction
1.0.1 Control Systems Terminologies
System: A system is a combination of components that act together and perform a
certain objective.
Linear Time-Invariant Systems (LTI systems) :
• LTI systems are a class of systems used in signals and systems that are
both linear and time-invariant.
• Linear systems are systems whose outputs for a linear combination of inputs
are the same as a linear combination of individual responses to those inputs. In
other words the superposition theorem is applied.
• Time-invariant systems are systems where the output does not depend
on when an input was applied. In other words if the input is delayed by t the
output is also delayed by t .
Control System: An interconnection of components forming a system
configuration that will provide a desired response. It can be described by a
block diagram as shown in Fig. 1.1.
• Manual Control
closed-loop control
Control Action open-loop control
(feedback control)
(a)
Program setting Washed clothes
Laundry
(Input) Machine (Output)
(b)
Fig. 1.5 (a) An open-loop system; (b) An open-loop system example.
9 of 70 GEE338 Ch.1 Math Models Prof Fawzy Ibrahim
1.0 Introduction
1.0.2 Control systems Classifications
2. Closed-loop control systems
• It is the system in which the control action is dependent on the output or the
input depends on the output as shown in Fig. 1.6.
• The input transducer converts the form of the input to the form used by the
controller.
• An output transducer or sensor measures the output response and converts it
into the form used by the controller. For example, if the controller uses electric
signal to operate the valve of a temperature control system: the input position
can be converted to voltage by a potentiometer (Input transducer) and the
output temperature can be converted to voltage by a thermistor (Output
transducer or Sensor).
• Closed-loop control systems. Feedback control systems are often referred to as
closed-loop control systems.
• In practice, the terms feedback control and closed-loop control are used
interchangeably.
• In a closed-loop control system the actuating error signal, which is the
difference between the input signal and the feedback signal , is fed to the
controller so as to reduce the error and bring the output of the system to a
desired value.
• Examples:- Refrigerator, Iron
10 of 70 GEE338 Ch.1 Math Models Prof Fawzy Ibrahim
1.0 Introduction
1.0.2 Control systems Classifications
2. Closed-loop control systems
(a)
Open-loop control
Closed-loop control systems
systems
• Simple construction, • High accuracy
ease of maintenance, • Not sensitive to disturbance
and less expensive.
• Controllable transient
Advantages • There is no stability
response
concern.
• Convenient when output • Controllable steady state
is hard to measure. error
• Disturbances and
changes in calibration • More Complex, and More
cause errors, and the Expensive.
output may be different • Possibility of instability.
Disadvantages
from what is desired.
• Recalibration is • Need for output
necessary from time to measurement.
time
13 of 70 GEE338 Ch.1 Math Models Prof Fawzy Ibrahim
1.0 Introduction
1.0.3 Control systems Applications
• Control systems have played a vital role in the advance of engineering and
science.
• Their Applications include:
- Robotic systems
- Aircraft Autopilot.
- Autonomous Cars.
- Missile guidance systems.
- Industrial processes.
- Transportations.
- Medicine.
- Health.
- Military.
- … etc.
C (s) 32 32
G (s)
s 12s 32 ( s 4)(s 8)
2
R( s)
16 of 70 GEE338 Ch.1 Math Models Prof Fawzy Ibrahim
1.1 Mathematical Model
Example 1.1 Solution
b) The impulse response function g(t) can be described as follows:
• Knowing the system transfer function, G(s) and the input, R(s) then output, C(s)
can be determined: C ( s) G ( s) R( s)
• Since the Laplace transform of the unit-impulse is unity the response of the
system on impulse is: If r (t ) (t ) and R ( s ) 1 C (s) G (s)
• Impulse response function g(t) is the inverse Laplace transform of G(s) as:
• Using partial fraction, we get: g (t ) L G ( s )
-1
32 32 A B
G(s) 2
s 12s 32 ( s 4)( s 8) ( s 4) ( s 8)
32 32
A ( s 4)G ( s) s 4 8
( s 8) s 4
4
8 8
G ( s)
32 32 ( s 4) ( s 8)
B ( s 8)G( s) s 8 8
( s 4) s 8 4
• Use the Laplace Transform tables
in Appendix A, g(t) is:
g (t ) 8(e 4t e 8t )u (t )
17 of 70 GEE338 Ch.1 Math Models Prof Fawzy Ibrahim
1.1 Mathematical Model
Example 1.1 Solution
c) The frequency response G(j) is defined as G(s) as s j:
32 32
G ( j ) G ( s ) s j G ( ) ( )
( s 4)( s 8) s j
( j 4)( j 8)
where |G()| is called the amplitude response and () is the phase response.
d) The Block diagram is shown in the following figure:
(a)
(b)
(c)
Fig. 1.6 Block diagram: (a) in time domain; (b) in s-domain and (c) in frequency domain
18 of 70 GEE338 Ch.1 Math Models Prof Fawzy Ibrahim
1.1 Mathematical Model
Example 1.2: A Linear Time Invariant (LTI) control system described by the
following transfer function G(s): 2
G (s)
For this system derive: ( s 1)( s 2)
a) The differential equation b) The impulse response g(t).
c) The frequency response G(j) d) draw its block diagram.
Solution
a) The system D.E. is calculated as follows:
C (s) 2 2
Then G (s) 2
R ( s ) ( s 1)( s 2) s 3s 2
s 2C ( s) 3sC ( s) 2C ( s) 2 R( s)
• Use the Laplace Transform tables in Appendix A, take the inverse Laplace
Transform for both sides, we have:
d 2c(t ) dc(t )
2
3 2c(t ) 2r (t )
dt dt
(a)
(b)
(c)
Fig. 1.7 Block diagram: (a) in time domain; (b) in s-domain and (c) in frequency domain
20 of 70 GEE338 Ch.1 Math Models Prof Fawzy Ibrahim
1.2 System Representation using Differential Equations
and Transfer Function
• The approach to dynamic system problems can be as follows:
1. Define the system and its components.
2. Formulate the mathematical model and list the necessary assumptions
3. Write the differential equations describing the model.
4. Solve the equations for the desired output variables.
5. Examine the solution and the assumptions.
6. If necessary reanalyze or redesign the system.
• The general form of the differential equation for LTI-System is given by:
d n c(t ) d n 1c(t )
an n
an 1 n 1
a0 c(t )
dt dt
d m r (t ) d m 1r (t ) (1.1)
bm m
bm 1 m 1
b0 r (t )
dt dt
• Where r(t) is the system input or excitation and c(t) is its output or response.
• The transfer function, G(s) of a LTI system is defined as the ratio of the
Laplace transform of the output (response function), C(s) and input (driving
function), R(s) under assumption that all initial conditions are zero.
C ( s) bm s m bm1s m1 b0
G( s) (1.2c)
R( s) an s n an1s n1 a0
bm s m bm 1s m 1 b0
R (s ) C (s )
an s n an 1s n 1 a0
Capacitor
Resistor Inductor
24 of 70 GEE338 Ch.1 Math Models Prof Fawzy Ibrahim
1.2 System Representation using DE and TF
1.2.1 DE and TF for Electrical and Electronic Elements
Example 1.3: For the series RLC circuit shown Write:
(a) the intego- differential equations
b) the transfer function of the following circuit, G(s) = Vs(s)/V(s).
Solution
a) Summing the voltages around the loop, assuming zero initial
conditions, yields the integro-differential equation as:
t
di(t ) 1
L Ri (t ) i ( )d v(t ) RLC circuit
dt C0
Changing variables from current to charge using i(t) = dq/dt, yields:
d 2 q (t ) dq(t ) 1
L R q (t ) v (t )
dt 2 dt C
the voltage-charge relationship for a capacitor q(t ) Cvc (t )
So, the system differential equation is:
d 2 vC (t ) dvC (t )
LC RC vC (t ) v(t )
dt 2 dt
b) Taking the Laplace transform assuming zero initial conditions, rearranging
terms, and simplifying yields
( s 2 LC sRC 1)VC ( s ) V ( s )
25 of 70 GEE338 Ch.1 Math Models Prof Fawzy Ibrahim
1.2 System Representation using DE and TF
1.2.1 DE and TF for Electrical and Electronic Elements
Example 1.3 Solution:
Solving for the transfer function, G(s) = Vc(s)/V(s) we obtain:
VC ( s ) 1 / LC
G (s)
V (s) R 1
s2 s
Another solution; L LC
• writing a mesh equation using the impedances as, we get Block diagram
V ( s ) ( R sL 1 / sC ) I ( s )
But the voltage across the capacitor is:
Vc ( s) I ( s) x1 / sC
Substitute we get the above results.
VC ( s ) 1 1 / LC
G (s)
V (s) LCs 2 RCs 1 s 2 R s 1
So, L LC
( s LC sRC 1)VC ( s ) V ( s )
2
Take the inverse Laplace transform, we get the differential equation as:
d 2 vC (t ) dvC (t )
LC 2
RC vC (t ) v(t )
dt dt
26 of 70 GEE338 Ch.1 Math Models Prof Fawzy Ibrahim
2.2 System Representation using DE and TF
1.2.1 DE and TF for Electrical and Electronic Elements
Example 1.3: Write the integodifferential equations and the transfer function of the
following circuit, G(s) = I2(s)/V(s). d iL 1
v L (t) L i L v L (t)dt
dt L
d vC 1
iC (t) C vC iC (t)dt
dt C
d
i 1 R1 L ( i1 i 2 ) v(t)
Time Domain dt
d t
1
i2 R2 L ( i 2 i1 )
c
i 2 dt 0
dt 0
Take Laplace Transform of both
sides then find G(s) = I2 (s)/V(s)
Laplace Domain
Block diagram
27 of 70 GEE338 Ch.1 Math Models Prof Fawzy Ibrahim
1.2 System Representation using DE and TF
2.2.1 DE and TF for Electrical and Electronic Elements
1. Write the equations that describe the dynamic behavior for each component.
2. Take Laplace transform of these equations, assuming zero initial conditions.
3. Represent each Laplace-transformed equation individually in block form.
4. Assembly the elements into a complete block diagram.
Example 1.4:
Find the transfer function, G(s) = V3(s)/V(s), for the circuit shown below.
Solution
Using KCL at each:
At node (2)
V2 ( s ) V1 ( s ) V2 ( s ) 0 V2 ( s ) V3 ( s )
0 (1)
R1 sL1 R3
At node (3)
V3 ( s ) V1 ( s ) V3 ( s ) V2 ( s ) V3 ( s ) 0
0 (2)
sL2 R2 sL3
Put L1 = L2 = L3 = 1 H and
R1 = R2 = R3 = 1
Eq. (1) becomes:
1 (1)
V1 ( s ) 1 1 V2 ( s ) V3 ( s ) 0
s
28 of 70 GEE338 Ch.1 Math Models Prof Fawzy Ibrahim
2.2 System Representation using DE and TF
1.2.1 DE and TF for Electrical and Electronic Elements
Example 1.4 Solution :
Find the transfer function, G(s) = V3(s)/V(s), for the circuit shown below.
Eq. (2) becomes:
V1 ( s ) 1 1
V2 ( s ) 1V3 ( s) 0 (2)
s s s
From Eq. (2) we get:
V1 ( s ) s 2
V2 ( s ) V3 ( s ) 0
s s
Substitute in Eq. (1) we get:
2s 1 V1 ( s) 2s 1
V1 ( s) V3 ( s) V3 ( s) 0 Block Diagram
s s s
Put V1(s) = V(s), then the transfer function, G(s) = V3(s)/V(s) is given by:
V3 ( s) (2s 1) / s 2 (2s 1)
G( s)
V ( s) (2s 3) / s s(2s 3)
Block diagram
ct r t
R2
G s
C ( s) R
R1 2
R( s ) R1 C s Z
G ( s) 2
Fig. 1.9 Inverting amplifier . Rs Z1
30 of 70 GEE338 Ch.1 Math Models Prof Fawzy Ibrahim
1.2 System Representation using DE and TF
1.2.1 DE and TF for Electrical and Electronic Elements
Operational Amplifiers (Op Amps) in system control
- Op Amps can be utilized for:
- Signal amplification in the sensor circuits
- Filters used for compensation purposes
- Modeling of the “real world” systems
- Lead or lag networks - Design of controllers Block diagram
C ( s) R2
c t
R2
1
r t G s
1
R1 R( s) R1
Fig. 1.10 Non-inverting amplifier.
31 of 70 GEE338 Ch.1 Math Models Prof Fawzy Ibrahim
1.2 System Representation using DE and TF
1.2.1 DE and TF for Electrical and Electronic Elements
Operational Amplifiers (Op Amps) in system control
Integrator
C s Z R2 R C s 1 RC R2 s 1 / R1C1
2 x 1 1 1 1 x x
R s Z1 R2C2 s 1 R1 R1 R2C2 s 1 / R2C2
C1 s 1/ R2C 2
x Where is
C2 s 1 /( ) R1C1
Fig. 1.12 Op Amp as a lead compensator ( < 1) as a lag compensator ( > 1) .
33 of 70 GEE338 Ch.1 Math Models Prof Fawzy Ibrahim
1.2 System Representation using DE and TF
1.2.1 DE and TF for Electrical and Electronic Elements
Operational Amplifiers (Op Amps) in system control
Z1
R1C1s 1R3
R1 R3 C1s 1 R2C2 s 1R4
Z2
R2 R4 C2 s 1
C s Z R
G (s) 2 x 6
Fig. 1.13 Op Amp as lead-lag network. R s Z1 R5
Forces Ma G(s)
V (s)
s
F ( s ) Ms 2 f v s K
Solution
1. Draw the free body diagram for the masses M1 and M2 as shown:
Forces Ma
2 N1 T2 N 2
1 N 2 T1 N1
bm s m bm 1s m 1 b0
R(s ) C (s )
an s n an 1s n 1 a0
System block diagram
41 of 70 GEE338 Ch.1 Math Models Prof Fawzy Ibrahim
1.3 Block Diagram Models
1.3.1 Components of a block diagram for a LTI system
G1G2G3
T ( s)
1 G1G2 H1 G2G3 H 2 G1G2G3
0.5
(0.15)
T (s) s(0.2s 1)( s 25) 0.075
0.075
015x0.5 s(0.2s 1)( s 25) 0.075 0.2s 3 6s 2 25s 0.075
1
s ( 0.2 s 1)( s 25)
56 of 70 GEE338 Ch.1 Math Models Prof Fawzy Ibrahim
1.4 Signal Flow Graph Model
1.4.1 Definitions
• Node: a point representing a signal or variable.
• Branch: unidirectional line segment joining two nodes.
• Path: a branch or a continuous sequence of branches that can be traversed
from one node to another node.
• Loop: a closed path that originates and terminates on the same node and
along the path no node is met twice.
• Nontouching loops: two loops are said to be nontouching if they do not have
a common node.
C (s) k p k k
T (s)
R(s)
• Where
∆ = 1 – Σ(All different loop gains)
+ Σ(Gain products of all combinations of two non-touching loops)
- Σ(Gain products of all combinations of three non-touching loops)
+ …
• Pk : The paths connecting the input R(s) and the output C(s).
• ∆k : is ∆ with the loops touching the kth path removed.
L1 G2 H 2 , L2 H 3G3 ,
L3 G6 H 6 , L4 G7 H 7
1 1 L3 L4 1 (G6 H 6 G7 H 7 )
2 1 L1 L2 1 (G2 H 2 G3 H 3 )
6. The TF of the system is
C ( s ) P11 P2 2
T (s)
R( s)
G1G2G3G4 1 (G6 H 6 G7 H 7 ) G5G6G7G8 1 (G2 H 2 G3 H 3 )
1 G2 H 2 G3 H 3 G6 H 6 G7 H 7 G2 H 2G6 H 6 G2 H 2G7 H 7 G3 H 3G6 H 6 G3 H 3G7 H 7
Where
∆ = 1 – Σ(All different loop gains)
+ Σ(Gain products of all combinations of two non-touching loops)
- Σ(Gain products of all combinations of three non-touching loops)
64 of 70 GEE338 Ch.1 Math Models Prof Fawzy Ibrahim
1.4 Signal Flow Graph Models
1.4.5 Block Diagram to Signal Flow Graph
Example 1.11 Solution
Pk : The paths connecting the input R(s) and the output C(s).
∆k : is ∆ with the loops touching the kth path removed
1. Calculate forward path transfer function Pk for each forward path k as:
1 1 1 100
P1 100 5 2 and P2 100
500
ss s s s
2. Calculate all loop TF’s as:
2 3 1 1 500 1 100
L1 , L2 , L3 100 5 1 2 and L4 100 1
s s s s s s s
3. L1 and L2 are nontouching loops at a time.
4. Calculate Δ as:
2 3 500 100 6
1 L1 L2 L3 L4 ( L1 L2 ) 1 2
105 506
2 1 2
s s s s s s s
5. Calculate Δk for each forward path as:
1 1 and 2 1
6. The TF of the system is calculated as:
500 100
x1 x1
C ( s ) P11 P2 2 s 2
s 100 s 500 100( s 5)
T (s) 2 2
R( s) 1
105 506
2 s 105s 506 s 105s 506
s s
65 of 70 GEE338 Ch.1 Math Models Prof Fawzy Ibrahim
1.4 Signal Flow Graph Models
1.4.5 Block Diagram to Signal Flow Graph
Example 1.11 Solution
b) The equivalent block diagram is shown below:
2. Eliminating the feedback loops: (1/s and 2), (1/s and 3) and Parallel (s and 5)
we get:
100( s 5)
( s 2)( s 3) 100( s 5) 100( s 5)
T (s) 2
100( s 5) ( s 2)( s 3) 100( s 5) ( s 5s 506)
1
( s 2)( s 3)
2
x1
C ( s ) P11 s 2
0 . 7 s 1 2( s 2 )
T (s) 3
R( s) 2 2 s ( s 2) s 4 . 7 s 2
4. 4 s 4
1
( s 2)( s 0.7 s 1)
2