Mathmatical Model of Control Systems PDF
Mathmatical Model of Control Systems PDF
2.1
Introduction
Mathematical models of control systems are mathematical expressions which describe the
relationships among system inputs, outputs and other inner variables. Establishing the
mathematical model describing the control system is the foundation for analysis and design of
control systems. Systems can be described by differential equations including mechanical systems,
electrical systems, thermodynamic systems, hydraulic systems or chemical systems etc. The
response to the input (the output of the system) can be obtained by solving the differential
equations, and then the characteristic of the system can be analyzed. The mathematical model
should reflect the dynamics of a control system and be suitable for analysis of the system. Thus,
when we construct the model, we should simplify the problem to obtain the approximate model
which satisfies the requirements of accuracy.
Mathematical models of control systems can be established by theoretical analysis or practical
experiments. The theoretical analysis method is to analyze the system according to physics or
chemistry rules (such as Kirchhoffs voltage laws for electrical systems, Newtons laws for
mechanical systems and Law of Thermodynamics). The experimental method is to approximate the
system by the mathematical model according to the outputs of certain test input signals, which is
also called system identification. System identification has been developed into an independent
subject. In this chapter, the theoretical analysis method is mainly used to establish the mathematical
models of control system.
There are a number of forms for mathematical models, for example, the differential equations,
difference equations and state equations in time domain, the transfer functions and block diagram
models in the complex domain, and the frequency characteristics in the frequency domain. In this
chapter, we shall study the differential equation, transfer function and block diagram formulations.
2.2
on the right hand side and the output related terms on the left hand side in descending order.
The following are some examples for establishing differential equations.
EXAMPLE 2-1
Consider the R-C-L network shown in Figure 2-1. Write the differential equation between the input
voltage u , and the output voltage u c.
Figure 2-1
R-L-C network
Solution:
This is an electrical system. According to Kirchhoffs voltage law, we have,
di (t )
u c (t )
dt
du (t )
i (t ) c c
dt
u r (t ) Ri (t ) L
(2-1)
(2-2)
Eliminating the intermediate variable i in the two equations above leads to:
d 2 u c (t )
du (t )
LC
RC c u c (t ) u r (t )
2
dt
dt
(2-3)
Assume that R, L, C are constants. The equation above is a second-order linear time-invariant
differential equation.
EXAMPLE 2-2
Consider the spring-mass-damper system shown in
Figure 2-2, K is the spring constant, F is
damping ratio of the damper, m is the mass of the car.
The friction between the car and the floor is ignored.
Obtain the differential equations of the system,
where the force F (t ) is the input and the
displacement y (t ) of the car is the output.
Solution
This is a mechanical system. The force analysis of the
car is shown in Figure 2-3. By Newtons second law,
the horizontal motion can be described as
F (t ) f
By setting T
Figure 2-2
Figure 2-3
dy (t )
d 2 y (t )
Ky (t ) m
dt
dt 2
f
m
,
K
2 mK
spring-mass-damper system
(2-4)
T2
d 2 y (t )
dy (t )
F (t )
2T
y (t )
2
dt
K
dt
(2-5)
The Equation (2-5) is a second-order linear time-invariant differential equation when K, f and m
are constants.
EXAMPLE 2-3
Try to write the differential equation of the
armature-controlled DC motor shown in Figure
2-4. The armature voltage u a (t ) is the input. The
rotate speed of the DC motor (t ) is the output.
Solution
It is an electrical-mechanical system. The armature-controlled DC motor converts the electrical
energy into the rotational mechanical energy. The input armature voltage u a (t ) provides the
armature current ia (t ) in the armature loop. The current ia (t ) and the excitation magnetic flux
work interactively to generate the magnetic torque M m (t ) on the rotor of the motor. Thus, the
load is driven. The differential equation of the motor consists of the following three parts.
The voltage balance equation of the armature loop:
u a (t ) La
dia (t )
Ra ia (t ) E a
dt
(2-6)
where: E a is the back electromotive-force (EMF) voltage proportional to the motor speed. That is
E a C e (t )
(2-7)
where C e (V / rad / s ) is the coefficient of the back electromotive-force voltage.
The motor torque is:
M m (t ) C m (t )ia (t )
(2-8)
where: C m ( N .m / A) is the torque coefficient of the motor. M m (t ) is the motor torque which is
generated by the armature current.
The torque balance equation on the motor shaft is:
Jm
d (t )
dt
f m (t ) M m (t ) M c (t )
(2-9)
La J m
d 2 (t )
dt
Cmu a (t ) La
( La f m Ra J m )
d (t )
dt
( Ra f m CmCe ) (t )
dM c (t )
Ra M c (t )
dt
(2-10)
Tm
where: Tm
Kc
d (t )
(t ) K a u a (t ) K c M c (t )
dt
(2-11)
Cm
Ra J m
is the time constant of the motor. K a
,
Ra f m C m C e
Ra f m C m C e
Ra
are the transmission coefficients of the motor. When Tm K1 K 2 are all
Ra f m C m C e
(t )
d (t )
into Equation (2-11) leads to
dt
d 2 (t) d
Tm
Kaua (t) Kc Mc (t)
dt
dt2
(2-12)
EXAMPLE 2-4
Figure 2-5 shows the control system for the rotating speed of a motor. Obtain the differential
equation of the system.
Figure 2-5 The control system for the rotating speed of a motor.
Solution
In this system, the given voltage ur (t ) is the input. The rotating speed (t ) of the motor is the
output. The differential equations of each component in the signal flow sequence from the error
signal are:
Measuring component. Tachometer is the measuring component, which transfers the system
output (t ) to the voltage u f (t ) :
u f (t ) K t (t )
(2-13)
where K t is the transfer function of the tachometer, which can be seen as a constant.
Comparing component. Comparing element compares the feedback voltage u f (t ) and the
given voltage u r (t ) to generate the error voltage u e (t ) :
(2-14)
u e (t ) u r (t ) u f (t )
Amplifying component. In this system, amplifying components are voltage amplifier and power
amplifier. They amplify the voltage and power of the error voltage u e (t ) :
First amplifier
u1 (t )
R2
ue (t )
R1
(2-15a)
Second amplifier
u 2 (t )
R4
u1 (t )
R3
(2-15b)
Power amplifier
ua (t ) K 3 u2 (t )
(2-15c)
Thus
u a (t ) K u e (t )
where K
(2-15d)
R2 R4
K 3 is amplifying constant for the amplifiers.
R1R3
Actuator. DC motor is the actuator. It converts the armature voltage u a (t ) to the shaft angle
speed 1 (t ) . According to Equation (2-11) in Example 2-3, the differential equation of DC motor
is:
Tm
d1 (t )
1 (t ) K a u a (t ) K c M c (t )
dt
(2-16)
Reducer. Reducer is to reduce the speed and increase the moment. The differential equation of
the reducer is :
(t )
1
1 (t )
i
(2-17)
KK a
K
d (t ) i KK a K t
u r (t ) c M c (t )
(
) (t )
dt
iTm
iTm
iTM
(2-18)
It is obvious from the above mathematical models that different components or systems may
have the same mathematical model. For example, the mathematical models in Example 2-1 and 2-2
are second-order differential equations, while the mathematical models in example 2-3 and 2-4 are
first-order differential equations. Systems with similar mathematical models are similar systems.
Similar systems have the similar characteristics although they are different in physical formality.
The similar principle reveals the relationship between different physical phenomena. The computer
simulation of control systems are based on the similar system idea.
u1 K 1
d (i )
dt
(2-20)
where K1 is a constant. The magnetic flux of the coil is a nonlinear function of the current i shown
in Figure2-6(b). Substituting Equation (2-20) into Equation (2-19) leads to:
K1
d (i ) di
Ri u r
di dt
(2-22)
(a)
(b)
Figure 2-6 A solenoid coil and the curved shape of its magnetic flux (i )
If the voltage and the current of the coil varies vary in a small range around the equilibrium point
( u 0 ,i 0 ) and
series:
(i) is smooth over the range of i 0 , the magnetic flux could be expressed in Taylor
1 d 2
d
i
(i ) 2
0
2
2! dt i
di i0
0
where i i i0 . The higher order terms could be neglected if i is small enough in value.
Thus
0
where
d
di
i0
d
di
i0
0 C 1 i
0 C 1 i
Thus, the relation between the current and magnetic flux can be approximated by linear equation.
Rewrite u r i in the equation (2-21) into the increment equation around the equilibrium point.
u r u 0 u r
i i0 i
0 C1 i
By substituting the above three equations into equation (2-21) and eliminating the intermediate
variables, we obtain:
K 1C1
di
Ri u r
dt
(2-22)
Equation (2-22) is the linearized increment differential equation of the solenoid coil. It is
usually written in the following form
di
(2-23)
Ri u r
dt
where rather than the actual value, ur and i now represent the increments around the equilibrium
K 1C1
point.
Systems that can not be linearized will be studied in Chapter 7.
2.2.3
The aim of the differential equation is to analyze the dynamical characteristic by the
mathematic method quantificationally. Thus the differential equation needs to be solved. The
differential equation can be transformed into algebraic equation in complex domain, and then the
analytical solution of the differential equation will be obtained by taking inverse Laplace transform
of the solution of the algebraic equation.
EXAMPLE 2-6
The R-C passive network is shown in Figure2-7, ur (t ) U r 1(t ) , uc (0) u0 . Obtain the
uc (t ) when K is closed.
duc
, we
dt
have
RC
duc (t )
uc (t ) ur (t )
dt
2-24
RC[ sU c ( s ) u0 ] U c ( s)
Ur
s
Solving U c ( s ) , we have
U c (s)
u0
Ur
RC
U
Ur
u0 r
1
s( RCs 1) RCs 1
s s 1
s
RC
RC
2-25
Taking the inverse Laplace transform of the Equation (2-25), the analytical solution of the
differential equation is
uc (t ) U r (1 e
t
RC
) u0 e
t
RC
2-26
The first term of the above equation is the particular solution of the ur (t ) which is called as
the zero state response; the second term is the homogeneous solution of the u0 which is called
zero input response.
2.2.4
The solution of the differential equation is constituted by the general solution of the
homogeneous solution and special solution of the given signal. The general solution reflects the
law of the natural response. If the eigenroot is
1 2 , n
e t e t , e t are called as the dynamical mode of the differential equation, also called as
1
mode shape.
te t
t 2 e t , .
j ,
e t sin t e t cos t .
e ( j ) t e ( j ) t can be written as
Every mode can be viewed as the basic dynamical mode of the natural response of the linear
system which is the linear combination of the corresponding mode.
2.3
The differential equation is the mathematical model of control systems in the time domain. If the
external excitation and the initial condition are given, all the information of the output with time
can be obtained by solving the differential equation. Although this method is accurate, it is tedious
to rewrite the equation when the structure or the parameters of the system are changed.
The transfer function is the mathematical model in the S-plane, which is based on Laplace
transform. By using transfer function, we can study the dynamic response of the system, as well as
the effect of the structure or the parameter variation. The root locus method and the frequency
response method are based on the transfer function. It is the most important and fundamental
mathematical model for control systems.
2.3.1
Transfer function
dnc(t)
dn1c(t)
dc(t)
a
...a1
a0c(t)
n1
n
n1
dt
dt
dt
an
bm
d m r (t )
d m 1r (t )
dr (t )
b
... b1
b0 r (t )
m 1
m
m 1
dt
dt
dt
(2-27)
where c(t) is the output, r(t) is the input, and a n , a n 1 ,..., a 0 and bm , bm 1 ,..., b0 are the
coefficients decided by the system structure and parameters.
To obtain the transfer function of the linear system, we simply take the Laplace transform on
both sides of the equation and assume zero initial conditions. The result is
a s
n
an 1s n 1 .... a1s a0 C ( s )
bm s m b m 1 s m 1 ... b1s b0 R( s)
(2-28)
C ( s ) bm s m bm 1 s m 1 ... b1 s b0
G ( s)
R( s ) a n s n a n 1 s n 1 ... a1 s a 0
(2-29)
The transfer function is obtained under the assumption that all initial conditions are zero.
The assumption indicates that i) the input is connected to the system after the point t=0 , that is,
the input and its higher order derivatives were zero when t 0 ; ii) the system was in equilibrium
before the input is connected to the system, that is, the output and its higher order derivatives were
zero when t 0 . Most systems in practice satisfy this assumption. Under this assumption, the
calculation is simplified. This assumption is also essential to compare system responses under the
same condition.
EXAMPLE 2-7
Obtain the transfer function of the R-L-C network in Example 2-1.
Solution:
From Equation (2-3) of Example 2-1, the differential equation of this R-L-C network is
d 2 u c (t )
du (t )
LC
RC c u c (t ) u r (t )
2
dt
dt
When all initial conditions are zero, take the Laplace transform on both sides of the equation. The
transfer function of the network is consequently written
G (s)
U c (s)
U r (s)
LCs
1
RCs 1
In practice, loading and interaction between interconnected components or systems may occur.
If the loading of interconnected devices does occur, the engineer must account for this change in
the transfer function and use the corrected transfer function in the subsequent calculations.
2Properties of the transfer function
1The transfer function is a rational fraction of the complex variable s. It has the properties
of the complex function. In practice, the highest power (n) of s in the denominator is greater or
equal to that (m) of the numerator, that is nm, because of the inertia of the system and the limited
power of the power source.
2The transfer function is a property of a system itself, independent of the magnitude and
nature of the input or driving function.
3The transfer function is dependent of the differential equations.
4The inverse Laplace transform of the transfer function gives the impulse response of the
system. Therefore, the transfer function reflects the dynamic characteristics of the system. Since the
Laplace transform of the impulse function is 1 (R(s)= L (t ) 1 ), we have
C ( s)
L1 G ( s) L1
L1 C ( s) k (t )
R( s)
5The transfer function is related to the poles and zeros in s-plane.
When it is idle load, relation between angular displacement (t ) of a single regulation resistance
brush and output voltage can be expressed as:
u (t ) K 1 (t )
(2-30)
Here: K1 E / max is the output voltage corresponding to the unit angular displacement of the
brush. It is called regulation resistance transfer coefficient. E is electric power source voltage of the
regulation resistance. max is the maximum working angle of regulation resistance(rad).
G ( s)
U (s)
K1
( s )
(2-31)
Equation (2-30) shows that the transfer function of the regulation resistance is a constant decided
by electric power source voltage E and maximum working angle max . The transfer function of the
regulation resistance can be represented by the block diagram as what is shown in Figure 2-8 (d).
When angular deviation detector is composed by a pair of same regulation resistances, the
output voltage is:
u (t ) u1 (t ) u 2 (t ) K 1 [ 1 (t ) 2 (t )] K 1 (t )
Here: K1 is the transfer function of a single regulation resistance. (t ) 1 (t ) 2 (t ) is the
difference between the angular displacement of the two regulation resistances. It is called deviation
angle. Therefore, the transfer functions of deviation detector and single regulation resistance has
the same form when deviation angle is the input. As
G ( s)
U (s)
K1
( s )
E u (t )
RP RP
'
Figure 2-9
regulation resistance
u (t )
RP
(2-32)
'
u (t )
Rl
After re-organizing the equation, the output voltage of the regulation resistance can be written as
below:
u (t )
E
Rp
'
R
R
p (1 p )
'
Rp
Rp Rl
E (t )
R (t )
(t )
max [1 p
(1
)]
max
Rl max
So, because of the influence of the loaded resistance Rl , the relation between the output voltage
u (t ) and the angular displacement (t ) is not he linear one anymore. If the loaded resistance
Rl is very strong, for example, when Rl 10 R p , it can be obtained approximately:
u (t ) E (t ) / max K 1 (t )
2Differential synchro
Differential synchro consists of a transmitter and a receiver (called control transformer in another
way). Figure 2-10 shows the elementary diagram and the method of connection.
Figure 2-10
Figure 2-11
Ks
E
E
r c e
(2-34)
(2-35)
As to the detail about the working principle of the differential synchro and the derivation, please
refer to the books about electric machine control.
The differential synchro and the deviation angle detector consisted of the potentiometer have
the same block plan, as shown in Figure 2-11. The difference is that the voltage generated by
differential synchro is AC voltage.
3Tachometer
Figure 2-12 shows the mechanism of a tachometer. The rotor of the tachometer is connected with
the rotating shaft of the measured plant. The output voltage is proportional to the angle speed of the
rotor. Thus
u K t K t
d
dt
(2-36)
where, is rotating angle of the rotor, is rotating speed, u is output voltage, and K t is slope of
the output voltage. When the rotor changes its rotating direction, the polarity and phase position of
the output voltage will be changed.
Figure 2-12
Tachometers
U ( s ) K t ( s ) K t s( s )
(2-37)
G ( s)
U (s)
Kt
( s )
or
G(s)
U ( s)
Kt s
( s )
(2-38)
The above two equations are the transfer functions of the tachometer. The output of the first one is
the rotating speed (t ) , and the output of the second one is angle position (t ) .
It is obvious that a system may have different transfer functions for different inputs and
outputs.
The structure of tachometer can be found in Figure2-13.
Figure 2-13
4DC motor
The differential equation of the DC motor is given in Example 2-3:
Tm
d (t )
(t ) K a u a (t ) K c M c (t )
dt
Suppose the initial condition is zero. By applying Laplace transform to the above equation, we
have the following transfer function.
Ga ( s )
Ka
( s )
U a ( s ) Tm s 1
(2-39)
U a ( s ) s (Tm s 1)
(2-40)
polarity.
The torque-speed characteristic curve of the twophase asynchronous motor is nonlinear and
the slope is negative. Figure 2-15(b) shows the measured curves for different control voltage u a .
Since the asynchronous motor usually works around the speed of zero, the linear part of the low
speed characteristic is extended to the nonlinear part of the high speed characteristic, which is
presented by the dotted line shown in Figure 2-15(b). In addition, we can also use the following
linear equation.
M m f mm M s
where M m is the output torque of the motor,
(2-41)
f m dM m d m is the damping coefficient, the slope of the dotted line, and M s is the torque to
lock the rotor. From Figure 2-14(b), we have
M s CM ua
(2-42)
where C M M s E .
Figure 2-15 Second-phase asynchronous electric motor and its performance curve
d 2 m
d
fm m
(2-43)
2
dt
dt
where m is angular displacement of the motor rotor, J m is the moment of inertia.
Eliminating the intermediate variables M s and M m in Equation (2-40) and (2-41) followed
Mm Jm
by
G(s)
where
m (s)
Km
CM
U a ( s ) s ( J m s f m ) s (Tm s 1)
(2-44)
G ( s)
Km
m ( s)
U a ( s ) Tm s 1
(2-45)
The transfer function and structure (Figure 2-16) of the two-phase asynchronous motors have
the same form with the DC motor.
Equations (2-44) and (2-45) are the transfer functions of the two-phase asynchronous electric
motor in different forms. It is obvious that they are same as the transfer functions of the DC motor.
6Gear trains
A gear train is a set or system of gears arranged to
transfer rotational torque from one part of a mechanical
system to another. Since motors with high rotating
speed and low torque are used often in control systems,
the gear train is often adopted to reduce the speed and
increase the torque.
Consider the gear train shown in Figure 2-17. The
rotating speed and the number of teeth of the driving
gear are 1 and Z 1 . Those for the driven gear are
and Z 2 .
The transmission ratio of the gear train is:
1 Z 2
2 Z1
(2-46)
In control system, the gear train is usually used to reduce the speed, thus i 1 . The transfer
function of the gear train is:
G ( s)
2 ( s) 1
1 ( s ) i
(2-47)
Figure 2-18
The moment of inertia and the viscous friction coefficient to the rotor shaft are
1
J2
i2
1
f f1 2 f 2
i
J J1
For gear trains with idler gears, the moment of inertia and the viscous friction coefficient to the
rotor shaft are
1
1 2
J J1 ( )2 J 2 (
) J3
i1
i1 i2
(2-48)
1
1 2
f f1 ( ) 2 f 2 (
) f3
i1
i1 i2
(2-49)
It is evident from equation (2-48) and (2-49) that the moment of inertia and the viscous friction
coefficient of the gears which are close to the motor shaft (or the input shaft or the gear train) have
the dominant effect on the load of the motor. Therefore, reducing the moment of inertia and the
viscous friction coefficient of the gears which are close to the motor shaft will improve the
transient performance of the motor.
Typical link
Num
Name of the
Differential Equation
Transfer Function
Example
Regulation
resistance
amplifier
differential
synchro
Gain
c K r
First-order
factor
T c c r
1
Ts 1
CR circuit
AC/DC motor
Quadratic
factor
T 2 c 2T c c r
0 1
1
2 2
T s 2Ts 1
R-L-C circuit
Spring-mass-da
mper system
Integral
factor
c r
1
s
Derivative
factor
c r
Reciprocal
first-order
factor
c r r
s 1
Reciprocal
quadratic
factor
c 2 r 2 r r
2 s 2 2 s 1
Note that different components may have the same transfer function while a component may have
different transfer functions for different inputs and outputs.
Basic factors are elements of transfer functions. They frequently occur in arbitrary transfer
functions.
G (s)
K * (s z j )
j 1
(2-50)
(s p )
i
i 1
where, z1 , z 2 ,..., z m are the zeros of the transfer function, and p1 , p 2 ,..., p n are the poles.
The coefficients of the highest order terms of both the denominator and numerator are unit.
The transfer function is in the highest 1 form. The denominator and numerator in this form can be
factorized into the poles and zeros form.
2The lowest 1 form ( Basic factor form)
Rewrite Equation 2-26 into Equation 2-48.
G ( s) K
m1
m2
k 1
n1
l 1
n2
( k s 1) ( l2 s 2 2 l s 1)
s
(T s 1) (T
i 1
j 1
(2-51)
2
j
s 2T j s 1)
2
The coefficients of the lowest order terms of both the denominator and numerator are unit. The
transfer function is in the lowest 1 form. The denominator and numerator in this form can be
factorized into basic factors.
*
The K is called as gain. The relationship between the K and K is
m
K * z j
j 1
2-52
p
i 1
EXAMPLE 2-8
The transfer function of closed-loop system is ( s )
30( s 2)
s ( s 3)( s 2 2 s 2)
Obtain the K ;
Obtain the differential equation;
Draw the zero-pole diagram.
Solution:
*
Zero-pole diagram
K
( s)
30 2
10
3 2
C (s)
30( s 2)
30( s 2)
4
2
R ( s ) s ( s 3)( s 2s 2) s 5s 3 8s 2 6s
( s 4 5s 3 8s 2 6s )C ( s ) 30( s 2) R( s )
Take the Laplace transform, we can obtain the differential equation:
d 4c(t )
d 3c ( t )
d 2 c(t )
dc(t )
dr (t )
5
8
6
30
60r (t )
dt
dt
dt
dt
dt
The Zero-pole diagram is shown in Figure 2-19.
2.4
dia (t )
ua (t ) La dt Ra ia (t ) Ea
Ea Cem (t )
M m (t ) Cm (t )ia (t )
dm (t )
f mm (t ) M m (t ) M c (t )
Jm
dt
Solution: Write the sub-block diagram for each sub-equation which is shown in Figure 2-20(a).
Connect the sub-block diagram to construct the block diagram of the system which is shown in
Figure 2-20(b).
. 2-20(a)
2-20(b)
I a ( s)
1
U ( s ) E ( s ) L s R
a
a
a
a
Ea ( s) Ce
m ( s)
M ( s)
m
Cm
I
a (s)
m (s)
1
M ( s) M ( s) J s f
c
m
m
m
Figure 2-20(b)
Figure 2-20(a)
sub-block diagram
EXAMPLE 2-10
Obtain the corresponding block diagram for the structure diagram of function recorder control
system shown in Figure 2-21.
Figure 2-21
Solution: From 2.3.2, replacing the name of each component by its transfer function,
representing all the variables by the Laplace form. Thus the block diagram can be found in Figure
2-22.
C ( s ) G2 ( s )U ( s ) G2 ( s )G1 ( s ) R ( s )
holds true.
Thus,
G ( s)
C (s)
G2 ( s )G1 ( s )
U (s)
(2-53)
C ( s) G1 ( s) R( s) G2 ( s) R( s ) [G1 ( s) G2 ( s)]R( s)
Thus the transfer function of the paralleled system is:
G ( s ) G1 ( s ) G2 ( s )
(2-54)
Transformatio
n
Original Diagram
Equivalent Diagram
Algebra
Combining
blocks in
cascade
C ( s) G1 ( s )G2 ( s) R( s)
Combining
blocks in
Prallel
C ( s) [G1 ( s) G2 ( s)]R( s)
G ( s) R( s)
Eliminating a
feedback loop
C ( s)
Moving a
summing
point ahead of
a block
C ( s) R( s)G ( s) Q( s)
1 G ( s) H ( s)
[ R( s)
Q( s)
]G ( s )
G ( s)
Moving a
summing
point behind a
block
C ( s ) [ R ( s ) Q ( s )]G ( s )
R ( s )G ( s ) Q( s )G ( s )
Moving a
pickoff point
ahead of a
block
C ( s ) G ( s ) R( s )
Moving a
pickoff point
behind a block
Changing the
position of a
summing
point and a
pickoff point
R ( s ) R( s )G ( s )
1
G (s)
C ( s ) G ( s) R( s)
C ( s) R1 ( s) R2 ( s)
C ( s ) G ( s ) E ( s ) G ( s )[ R( s ) B( s )] G ( s )[ R ( s ) H ( s )C ( s )]
Thus
G ( s)
R( s )
1 G ( s) H ( s)
Therefore, the transfer function relating the output C ( s ) to the input R( s ) is
G ( s)
(s)
1 G ( s) H ( s)
C ( s)
(2-55)
When the transfer function of the feedback path is H ( s ) 1 , the system is a unity feedback system.
Consequently, the closed-loop transfer function is:
( s)
G( s)
1 G( s)
(2-56)
Table 2-2 summarizes the basic rules for moving the summing point and the pickoff point.
EXAMPLE 2-11
Reduce the block diagram shown in Figure 2-26. Obtain the closed-loop transfer function
( s)
C (s)
.
R( s)
Solution:
The block diagram reduction procedure is based on the utilization of Table 2-2, item 3, which
eliminates feedback loops. Therefore the other transformations are used to transform the diagram to
a form ready for eliminating feedback loops.
First, Moving the pickoff point a behind and the summing point b behind the block G2 ( s ) ,
G2 ( s)
in cascade and H 2 ( s ) in parallel,
G4 ( s)
( s)
C (s)
R( s )
It is interesting to consider if the pickoff point a can be moved before the summing point d directly.
We leave this problem to the readers.
2.5
A signal-flow graph (SFG) is a diagram that represents a set of simultaneous linear algebraic
equations. A signal flow graph contains essentially the same information as a block diagram. The
advantage of the signal flow graph is the availability of a flow graph gain formula, which provides
the relation between system variables without requiring any reduction procedure or manipulation of
the flow graph.
BranchRepresented by . A branch is a directed line segment joining two nodes. The arrow
head represents the direction of the signal transmission.
GainRepresented by the transfer function G . The gain is equivalent to the transfer function in
the block diagram. It represents the relation between the signals in the arrow head direction.
Figure 2-28(b) shows the corresponding signal flow graph of the block diagram shown in
Figure 2-28(a).
Before we discuss signal-flow graphs, we must define certain terms.
Input node or source. An input node or source is a node that has only outgoing branches. This
corresponds to an independent variable. Nodes R and N in Figure 2-28(b) are input nodes. They are
actually input signals.
Output node or sink. An output node or sink is a node that has only incoming branches. This
corresponds to a dependent variable. Node C in Figure 2-28 (b) is a output node. It is actually the
output signal.
Mixed node. A mixed node is a node that has both incoming and outgoing branches. Nodes E, P
and Q in Figure 2-28(b) are mixed nodes. The mixed mode is the equivalent of the pickoff point in
block diagrams.
Forward path. A forward path is a path from an input node (source) to an output node (sink) that
does not cross any nodes more than once. The paths REPQC and NPQC in Figure 2-28(b) are
forward paths.
Loop. A loop is a closed path that starts and ends at the same node and does not cross any other
nodes more than once. Path EPQE in Figure 2-28(b) is a loop.
Loop gain. The loop gain is the product of the branch transmittances of a loop.
Forward path gain. A forward path gain is the product of the branch transmittances of a forward
path.
Non-touching loops. Loops are non-touching if they do not possess any common nodes.
G ( s)
1 n
Pk k
k 1
(2-57)
1 La Lb Lc Ld Le L f
(2-58)
n = total number of forward paths between the output node and the input node.
Pk path gain or transmittance of kth forward path.
k cofactor of the kth forward path determinant of the graph with the loops touching the Kth
forward path removed, that is, the cofactor k is obtained from by removing the loops that
touch path Pk .
EXAMPLE 2-12
Consider the system shown in Figure 2-29. Obtain the closed-loop transfer function
Figure 2-29
C (s)
.
R(s)
Solution:
There are three forward paths and the forward path gains are P1 G1G2 G3 G4 , P2 G5G3G4
and P3 G1G6 . There are three loops and the loop gains are L1 G1G2G3G4 H 2 ,
L2 G1G6 H 2 and L3 G3 H1 . From Figure 2-27, we can see that there is only one pair of
nontouching loops, that is, the two loops are L2 and L3 .
Thus, the determinant is given by:
1 ( L1 L2 L3 ) ( L1 L3 )
1 G1G2G3G4 H 2 G1G6 H G3 H1 G1G3G6 H1 H 2
All loops are in touch with the forward path P1 and P2 . Thus 1 2 1 . Since the forward
path P3 does not touch the loop L3 , we obtain the cofactor 3 1 ( L3 ) 1 G3 H1 .
Therefore the closed-loop transfer function is given by:
C (s) 1
( P11 P2 2 P3 3 )
R( s)
G1G2G3G4 G3G4G5 G1G6 (1 G3 H1 )
Consider the system shown in Figure 2-30. Obtain the closed-loop transfer function
Figure 2-30
C ( s)
.
R( s)
Solution:
The following results are obtained by inspection of the block diagram.
There are two forward paths. The forward path gains and their cofactors are
p1 G1G2G3 ,
p2 H 4 ,
1 1,
2 1 G2G3 2
There are four individual loops. The gains of these loops are
L1 H 3 H 4 , L2 G1G2 G3 H 3 , L3 G2 G3 H 2 , L4 G1 H 1 .
There are two pairs of nontouching loops, that is, Loop L1 does not touch and loop L3 does not
touch L4 .
Using Equation 2-58, the closed-loop transfer function is written
R(s)
2.6
G1G2G3 H4(1G2G3H2
1H3H4 G1G2G3H3 G2G3H2 G1H1 G2G3H2H3H4 G1G2G3H1H2
G ( s) H ( s)
B( s )
G1 ( s )G2 ( s ) H ( s)
E ( s)
(2-59)
Note that the open loop transfer function is obtained by opening the closed-loop at the output of
H (s ) , or the summing point. Differing from the transfer function of open loop systems, it is for
closed-loop systems.
( s)
C ( s)
G1 ( s)G2 ( s)
N ( s ) 1 G1 ( s )G2 ( s) H ( s )
(2-60)
n ( s)
C (s)
G2 ( s )
N ( s ) 1 G1 ( s )G2 ( s ) H ( s )
(2-61)
C ( s)
G1 ( s)G2 ( s) R( s)
G2 ( s) N ( s )
1 G1 ( s)G2 ( s ) H ( s) 1 G1 ( s)G2 ( s) H ( s )
2.6.3 Transfer Functions between The Error Signal and The Input Signal
1The error transfer function subjected to the input signal
By setting the disturbance signal N ( s ) 0 , we obtain the error transfer function subjected to the
input.
e ( s)
E ( s)
1
R( s ) 1 G1 ( s )G2 ( s ) H ( s )
(2-62)
ne ( s )
E (s)
G2 ( s ) H ( s )
N ( s ) 1 G1 ( s )G2 ( s ) H ( s )
3The error subjected to both the input and the output signals
(2-62)
From Equation (2-63) and (2-64), the error subjected to both the input and the output signals is
E ( s)
G2 ( s ) H ( s) N ( s)
R( s )
1 G1 ( s)G2 ( s ) H ( s) 1 G1 ( s)G2 ( s) H ( s )
n ( s) ,
n (t ) (t )
Solution: The open-loop transfer function is
2
G (s)
s ( s 3)
The closed-loop transfer function is
Figure 2-32 block diagram of the system
2
C(S )
2
s( s 3)
( s)
2
R( s) 1
( s 1)( s 2)
s( s 3)
1
C (S )
s
s3
n (s)
( s 1)( s 2)
N (s) 1 2
s ( s 3)
The output is
C ( s ) ( s ) R ( s ) n ( s ) N ( s ))
2
1
s
( s 1)( s 2) s ( s 1)( s 2)
s2 2
1
3
3
s ( s 1)( s 2)
s s 1 s 2
c(t ) 1 3e t 3e2t
We should notice that the system is unit feedback when we solve e(t ) .
Summary
Mathematical modeling is the foundation of the control theory and the precondition for analyzing
and designing the control system. The contents of this chapter are mainly: four types of
mathematical models for control systems, three methods for obtaining transfer functions, and five
forms of the transfer functions of control systems.
This chapter focuses on the construction of the mathematical model of the system by the analysis
method. The process is shown in Figure 2-33
Figure 2-33 The process of constructing the mathematical model for the system
Exercises
E2-1
impedance method.
The mechanical system is shown in Figure 2-36(a), and the circuit system in Figure 2-36(b). Prove they are
similar systems, namely, they have the same form of mathematical model.
The diode shown in Fig.2-37 is a nonlinear component. The relationship between the current
voltage ud is id 10
14
(e
ud
0.026
id and
working point u (0) 2.39V , i (0) 2.19 10 A , try to find the linearization equation
for the id f (ud ) .
Figure 2-37
E2-5
dh
dt S
h and input Qr of a
1
Qr
S
h in terms of
is the length of
m.
rod,
is the angle, the mass is
Obtain the differential
equation of the system and linearized it.
E2-7 Obtain the image functions
Figure 2-38
Pendulum system
in Figure 2-39.
E2-8
(1)
X (s)
(2)
X ( s)
(3)
(4)
E2-9
e s
s 1
2
s 9
2
1
s ( s 2) 3 ( s 3)
s 1
X ( s)
s ( s 2 2 s 2)
X ( s)
Under the zero initial condition, a unit step response of the system can be represented as
c(t ) 1 2e 2t et
Determine the transfer function and impulsive response of the system.
E2-10
C (s)
2
2
R ( s ) s 3s 2
and the initial condition is c(0) 1 , c(0)
0.
r (t ) 1(t ) .
E2.11 Determine the transfer function
Qc ( s) .
Qr ( s)
Qc ( s) .
Qr ( s)
E2-14 The equation set of the system is shown as follows. Draw the block diagram and determine the closed-loop
transfer function
C (s)
.
R( s)
X 1 ( s ) G1 ( s ) R( s ) G1 ( s )[G7 ( s ) G8 ( s )]C ( s )
X ( s ) G ( s )[ X ( s ) G ( s ) X ( s )]
2
2
1
6
3
X 3 ( s ) [ X 2 ( s ) C ( s )G5 ( s )]G3 ( s )
C ( s ) G4 ( s ) X 3 ( s )
E2-15 Sketch the structure and signal diagram for the RC passive network shown in Figure2-43 and obtain the
transfer function
U c ( s)
.
U r ( s)
Figure 2-43
E2-16 Obtain the transfer functions
RC passive network
C (s)
for the systems which are shown in Figure 2-44 by the method of
R( s)
Figure 2-44
E2-17 A control system is shown in Figure 2-45. Draw the signal-flow diagram
Figure 2-45
E2-18 Draw the structure of system for the signal-flow which is shown in Figure 2-46.
Figure 2-47
C (s) C ( s)
,
.
R( s) E ( s )
G0
G1
G2
G3
H2
Figure 2-48
C (s) C ( s)
,
.
R( s) N ( s )
Figure 2-49
E2-22 A double pendulum system hung on the swivel without friction is shown in Figure 2-50. Suppose the
is small,
l is
Figure 2-50
double pendulums