Control-Systems - FINAL -محاضرات 2019-2020
Control-Systems - FINAL -محاضرات 2019-2020
Control-Systems - FINAL -محاضرات 2019-2020
د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
التعرف على انظمة الديطرة الطوعية وطرق. يعتبر التعرف على مكونات دوائر الديطرة و أنواع المديطرات و استخداماتها و دوائرها التطبيقية:أهداف المادة
تحليلها وتمثيلها وتطبيقاتها العملية في المجاالت الهندسية
Week Syllabus
6th, 7th , 8th Transfer function, block diagram representation and reduction diagram.
17th , 18th Frequency domain analysis, phase margin and bode plot.
P, PI, PD, and PID Modes of Feedback Control, Realization of PID Controller
23rd, 24th
Using Active and Passive Elements.
25th, 26th , 27th Control system design using root locus method.
1
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
References:
Northern Technical University:
1.Feedback control system and synthesis: By J J .Dazzo
2. Modern control engineering: By OGATA. Engineering Technical College:
3. Automatic control engineering By: Francis Raven.
4. Linear control system By: B.S.Manke. Dept.: Medical Instruments Eng. Technology
5.Feedback and control system .SCHAUM’S OUTLINE
Stage: Fourth class
E-mail [email protected]
[email protected]
[email protected]
Course Coordinator
Course Objective
Course Description
Text Book
References
2
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Here, an input is applied to a controller and it produces an actuating signal or controlling signal.
This signal is given as an input to a plant or process which is to be controlled. So, the plant
produces an output, which is controlled. The traffic lights control system which we discussed
earlier is an example of an open loop control system.
In closed loop control systems, output is feedback to the input. So, the control action is
dependent on the desired output.
The following figure shows the block diagram of negative feedback closed loop control system.
The error detector produces an error signal, which is the difference between the input and
the feedback signal. This feedback signal is obtained from the block (feedback elements) by
3
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
considering the output of the overall system as an input to this block. Instead of the direct input,
the error signal is applied as an input to a controller.
So, the controller produces an actuating signal which controls the plant. In this combination, the
output of the control system is adjusted automatically till we get the desired response. Hence,
the closed loop control systems are also called the automatic control systems. Traffic lights
control system having sensor at the input is an example of a closed loop control system.
The differences between the open loop and the closed loop control systems are mentioned in the
following table.
Open Loop Control Systems Closed Loop Control Systems
Control action is independent of the desired Control action is dependent of the desired
output. output.
These are also called as non-feedback control These are also called as feedback control
systems. systems.
Inaccurate. Accurate.
Positive Feedback:
The positive feedback adds the reference input, R(s) and feedback output. The following figure
shows the block diagram of positive feedback control system.
4
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
The concept of transfer function will be discussed in later chapters. For the time being, consider
the transfer function of positive feedback control system is,
T=G/1−GH (Equation 1)
Where,
T is the transfer function or overall gain of positive feedback control system.
G is the open loop gain, which is function of frequency.
H is the gain of feedback path, which is function of frequency.
Negative Feedback
Negative feedback reduces the error between the reference input, R(s) and system output. The
following figure shows the block diagram of the negative feedback control system.
7
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Example
Consider the following electrical system as shown in the following figure. This circuit consists
of resistor, inductor and capacitor. All these electrical elements are connected in series. The
input voltage applied to this circuit is Vi and the voltage across the capacitor is the output
voltage Vo.
Here, we represented an LTI system with a block having transfer function inside it. And this
block has an input X(s) & output Y(s).
8
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Example
Previously, we got the differential equation of an electrical system as
d2Vo/dt2+(RL)dVo /dt +(1/LC)Vo = (1/LC)Vi
Apply Laplace transform on both sides.
S2 Vo(s)+(SRL)Vo(s)+(1/LC)Vo(s)=(1/LC)Vi(s)
⇒{S2 +(RL)S+1/LC}Vo(s)=(1/LC)Vi(s)
⇒Vo(s)/Vi(s)=1/[LCS2 +(RL)S+1/LC]
Where,
Vi(s) is the Laplace transform of the input voltage Vi
Vo(s) is the Laplace transform of the output voltage Vo
The above equation is a transfer function of the second order electrical system. The transfer
function model of this system is shown below.
Here, we show a second order electrical system with a block having the transfer function inside
it. And this block has an input Vi(s) & an output Vo(s).
Fm∝a
⇒Fm=Ma=Md2x/dt2
F=Fm=Md2x/dt2
Where,
F is the applied force
Fm is the opposing force due to mass
M is mass
a is acceleration
x is displacement
Spring
Spring is an element, which stores potential energy. If a force is applied on spring K, then it is
opposed by an opposing force due to elasticity of spring. This opposing force is proportional to
the displacement of the spring. Assume mass and friction are negligible.
10
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
F∝x
⇒Fk=Kx
F=Fk=Kx
Where,
F is the applied force
Fk is the opposing force due to elasticity of spring
K is spring constant
x is displacement
Dashpot
If a force is applied on dashpot B, then it is opposed by an opposing force due to friction of the
dashpot. This opposing force is proportional to the velocity of the body. Assume mass and
elasticity are negligible.
Fb∝ν
⇒Fb=Bν=Bdx/dt
F=Fb=Bdx/dt
Where,
Fb is the opposing force due to friction of dashpot
B is the frictional coefficient
v is velocity
x is displacement
11
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
If a torque is applied to a rotational mechanical system, then it is opposed by opposing torques
due to moment of inertia, elasticity and friction of the system. Since the applied torque and the
opposing torques are in opposite directions, the algebraic sum of torques acting on the system is
zero. Let us now see the torque opposed by these three elements individually.
Moment of Inertia
In translational mechanical system, mass stores kinetic energy. Similarly, in rotational
mechanical system, moment of inertia stores kinetic energy.
If a torque is applied on a body having moment of inertia J, then it is opposed by an opposing
torque due to the moment of inertia. This opposing torque is proportional to angular
acceleration of the body. Assume elasticity and friction are negligible.
Tj∝α
⇒Tj=Jα=Jd2θdt2
T=Tj=Jd2θdt2
Where,
T is the applied torque
Tj is the opposing torque due to moment of inertia
J is moment of inertia
α is angular acceleration
θ is angular displacement
Torsional Spring
In translational mechanical system, spring stores potential energy. Similarly, in rotational
mechanical system, torsional spring stores potential energy.
If a torque is applied on torsional spring K, then it is opposed by an opposing torque due to the
elasticity of torsional spring. This opposing torque is proportional to the angular displacement
of the torsional spring. Assume that the moment of inertia and friction are negligible.
12
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Tk∝θ
⇒Tk=Kθ
T=Tk=Kθ
Where,
T is the applied torque
Tk is the opposing torque due to elasticity of torsional spring
K is the torsional spring constant
θ is angular displacement
Dashpot
If a torque is applied on dashpot B, then it is opposed by an opposing torque due to
the rotational friction of the dashpot. This opposing torque is proportional to the angular
velocity of the body. Assume the moment of inertia and elasticity are negligible.
Tb∝ω
⇒Tb=Bω=Bdθdt
T=Tb=Bdθdt
Where,
Tb is the opposing torque due to the rotational friction of the dashpot
B is the rotational friction coefficient
ω is the angular velocity
θ is the angular displacement
13
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Electrical Analogies of Mechanical Systems
Two systems are said to be analogous to each other if the following two conditions are satisfied.
The two systems are physically different
Differential equation modelling of these two systems are same
Electrical systems and mechanical systems are two physically different systems. There are two
types of electrical analogies of translational mechanical systems. Those are force voltage
analogy and force current analogy.
Force Voltage Analogy
In force voltage analogy, the mathematical equations of translational mechanical system are
compared with mesh equations of the electrical system.
Consider the following translational mechanical system as shown in the following figure.
Consider the following electrical system as shown in the following figure. This circuit consists
of a resistor, an inductor and a capacitor. All these electrical elements are connected in a series.
The input voltage applied to this circuit is V volts and the current flowing through the circuit
is i Amps.
14
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Mesh equation for this circuit is
V=Ri+Ldi/dt+1/c∫idt (Equation 2)
Substitute, i =dq/dt in Equation 2.
V=Rdq/dt+Ld2q/dt2+q/C
⇒V=Ld2q/dt2+Rdq/dt+(1/c)q (Equation 3)
By comparing Equation 1 and Equation 3, we will get the analogous quantities of the
translational mechanical system and electrical system. The following table shows these
analogous quantities.
Translational Mechanical System Electrical System
Force(F) Voltage(V)
Mass(M) Inductance(L)
Displacement(x) Charge(q)
Velocity(v) Current(i)
Consider the following electrical system as shown in the following figure. This circuit consists
of current source, resistor, inductor and capacitor. All these electrical elements are connected in
parallel.
The above block diagram consists of two blocks having transfer functions G(s) and H(s). It is
also having one summing point and one take-off point. Arrows indicate the direction of the flow
of signals. Let us now discuss these elements one by one.
Block
The transfer function of a component is represented by a block. Block has single input and
single output.
The following figure shows a block having input X(s), output Y(s) and the transfer function
G(s).
16
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Transfer Function, G(s)=Y(s)/X(s)
⇒Y(s)=G(s)X(s)
Output of the block is obtained by multiplying transfer function of the block with input.
Summing Point
The summing point is represented with a circle having cross (X) inside it. It has two or more
inputs and single output. It produces the algebraic sum of the inputs. It also performs the
summation or subtraction or combination of summation and subtraction of the inputs based on
the polarity of the inputs. Let us see these three operations one by one.
The following figure shows the summing point with two inputs (A, B) and one output (Y).
Here, the inputs A and B have a positive sign. So, the summing point produces the output, Y
as sum of A and B.
i.e., Y = A + B.
The following figure shows the summing point with two inputs (A, B) and one output (Y).
Here, the inputs A and B are having opposite signs, i.e., A is having positive sign and B is
having negative sign. So, the summing point produces the output Y as the difference of A and
B.
Y = A + (-B) = A - B.
The following figure shows the summing point with three inputs (A, B, C) and one output (Y).
Here, the inputs A and B are having positive signs and C is having a negative sign. So, the
summing point produces the output Y as
17
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Y = A + B + (−C) = A + B − C.
Take-off Point
The take-off point is a point from which the same input signal can be passed through more than
one branch. That means with the help of take-off point, we can apply the same input to one or
more blocks, summing points.
In the following figure, the take-off point is used to connect the same input, R(s) to two more
blocks.
In the following figure, the take-off point is used to connect the output C(s), as one of the inputs
to the summing point.
18
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
By applying the Laplace transform to this circuit, will get the circuit in s-domain. The circuit is
as shown in the following figure.
19
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Equation 2 can be implemented with a block having transfer function, 1sC. The input and
output of this block are I(s) and Vo(s). The block diagram of Equation 2 is shown in the
following figure.
20
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
The overall block diagram of the series of RLC Circuit (s-domain) is shown in the following
figure.
Similarly, you can draw the block diagram of any electrical circuit or system just by following
this simple procedure.
Convert the time domain electrical circuit into an s-domain electrical circuit by applying
Laplace transform.
Write down the equations for the current passing through all series branch elements and voltage
across all shunt branches.
Draw the block diagrams for all the above equations individually.
Combine all these block diagrams properly in order to get the overall block diagram of the
electrical circuit (s-domain).
Where, Z(s)=G1(s)X(s)
⇒Y(s)=G2(s)[G1(s)X(s)]=G1(s)G2(s)X(s)
⇒Y(s)={G1(s)G2(s)}X(s)
Compare this equation with the standard form of the output equation, Y(s)=G(s)X(s).
Where, G(s)=G1(s)G2(s).
That means we can represent the series connection of two blocks with a single block. The
transfer function of this single block is the product of the transfer functions of those two blocks.
The equivalent block diagram is shown below.
Similarly, you can represent series connection of ‘n’ blocks with a single block. The transfer
function of this single block is the product of the transfer functions of all those ‘n’ blocks.
Parallel Connection
The blocks which are connected in parallel will have the same input. In the following figure,
two blocks having transfer functions G1(s) and G2(s) are connected in parallel. The outputs of
these two blocks are connected to the summing point.
22
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Compare this equation with the standard form of the output equation, Y(s)=G(s)X(s).
Where, G(s)=G1(s)+G2(s).
That means we can represent the parallel connection of two blocks with a single block. The
transfer function of this single block is the sum of the transfer functions of those two blocks.
The equivalent block diagram is shown below.
Similarly, you can represent parallel connection of ‘n’ blocks with a single block. The transfer
function of this single block is the algebraic sum of the transfer functions of all those ‘n’ blocks.
Feedback Connection
As we discussed in previous chapters, there are two types of feedback — positive feedback and
negative feedback. The following figure shows negative feedback control system. Here, two
blocks having transfer functions G(s) and H(s) form a closed loop.
23
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
This means we can represent the negative feedback connection of two blocks with a single
block. The transfer function of this single block is the closed loop transfer function of the
negative feedback. The equivalent block diagram is shown below.
Similarly, you can represent the positive feedback connection of two blocks with a single block.
The transfer function of this single block is the closed loop transfer function of the positive
feedback, i.e., G(s)/1−G(s)H(s)
Block Diagram Algebra for Summing Points
There are two possibilities of shifting summing points with respect to blocks −
Shifting summing point after the block
Shifting summing point before the block
Let us now see what kind of arrangements need to be done in the above two cases one by one.
Shifting Summing Point After the Block
Consider the block diagram shown in the following figure. Here, the summing point is present
before the block.
Summing point has two inputs R(s) and X(s). The output of it is {R(s)+X(s)}.
So, the input to the block G(s) is {R(s)+X(s)} and the output of it is –
Y(s)=G(s){R(s)+X(s)}
⇒Y(s)=G(s)R(s)+G(s)X(s) (Equation 1)
Now, shift the summing point after the block. This block diagram is shown in the following
figure.
24
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
25
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
26
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
27
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Here, X(s)=Y(s)=G(s)R(s)
When you shift the take-off point before the block, the output Y(s) will be same. But, there is
difference in X(s) value. So, in order to get same X(s) value, we require one more block G(s). It
is having the input R(s) and the output is X(s). This block diagram is shown in the following
figure.
28
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Example: Consider the block diagram shown in the following figure. Let us simplify (reduce)
this block diagram using the block diagram reduction rules.
Step 1 − Use Rule 1 for blocks G1 and G2. Use Rule 2 for blocks G3 and G4. The modified
block diagram is shown in the following figure.
29
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Step 2 − Use Rule 3 for blocks G1G2 and H1. Use Rule 4 for shifting take-off point after the
block G5. The modified block diagram is shown in the following figure.
Step 3 − Use Rule 1 for blocks (G3+G4) and G5. The modified block diagram is shown in the
following figure.
30
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Step 4 − Use Rule 3 for blocks (G3+G4)G5 and H3. The modified block diagram is shown in
the following figure.
Step 5 − Use Rule 1 for blocks connected in series. The modified block diagram is shown in the
following figure.
Step 6 − Use Rule 3 for blocks connected in feedback loop. The modified block diagram is
shown in the following figure. This is the simplified block diagram.
31
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Note − Follow these steps in order to calculate the transfer function of the block diagram having
multiple inputs.
Step 1 − Find the transfer function of block diagram by considering one input at a time and
make the remaining inputs as zero.
Step 2 − Repeat step 1 for remaining inputs.
Step 3 − Get the overall transfer function by adding all those transfer functions.
The block diagram reduction process takes more time for complicated systems. Because, we
have to draw the (partially simplified) block diagram after each step. So, to overcome this
drawback, use signal flow graphs (representation).
In the next two chapters, we will discuss about the concepts related to signal flow graphs, i.e.,
how to represent signal flow graph from a given block diagram and calculation of transfer
function just by using a gain formula without doing any reduction process.
The nodes present in this signal flow graph are y1, y2, y3 and y4.
y1 and y4 are the input node and output node respectively.
y2 and y3 are mixed nodes.
Branch
Branch is a line segment which joins two nodes. It has both gain and direction. For example,
there are four branches in the above signal flow graph. These branches have gains of a, b,
c and -d.
Construction of Signal Flow Graph
Let us construct a signal flow graph by considering the following algebraic equations −
y2=a12y1+a42y4
y3=a23y2+a53y5
y4=a34y3
y5=a45y4+a35y3
y6=a56y5
There will be six nodes (y1, y2, y3, y4, y5 and y6) and eight branches in this signal flow graph.
The gains of the branches are a12, a23, a34, a45, a56, a42, a53 and a35.
To get the overall signal flow graph, draw the signal flow graph for each equation, then
combine all these signal flow graphs and then follow the steps given below −
Step 1 − Signal flow graph for y2=a13y1+a42y4 is shown in the following figure.
33
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Step 2 − Signal flow graph for y3=a23y2+a53y5 is shown in the following figure.
Step 3 − Signal flow graph for y4=a34y3 is shown in the following figure.
Step 4 − Signal flow graph for y5=a45y4+a35y3 is shown in the following figure.
Step 5 − Signal flow graph for y6=a56y5 is shown in the following figure.
Step 6 − Signal flow graph of overall system is shown in the following figure.
34
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
35
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Represent the input signal R(s) and output signal C(s) of block diagram as input node R(s) and
output node C(s) of signal flow graph.
Just for reference, the remaining nodes (y1 to y9) are labelled in the block diagram. There are
nine nodes other than input and output nodes. That is four nodes for four summing points, four
nodes for four take-off points and one node for the variable between blocks G1 and G2.
The following figure shows the equivalent signal flow graph.
With the help of Mason’s gain formula (discussed in the next chapter), you can calculate the
transfer function of this signal flow graph. This is the advantage of signal flow graphs. Here, we
no need to simplify (reduce) the signal flow graphs for calculating the transfer function.
Mason's Gain Formula
Let us now discuss the Mason’s Gain Formula. Suppose there are ‘N’ forward paths in a signal
flow graph. The gain between the input and the output nodes of a signal flow graph is nothing
but the transfer function of the system. It can be calculated by using Mason’s gain formula.
Mason’s gain formula is
T=C(s)R(s)=ΣNi=1PiΔiΔ
Where,
C(s) is the output node
R(s) is the input node
T is the transfer function or gain between R(s) and C(s)
Pi is the ith forward path gain
Δ=1−(sum of all individual loop gains)
+(sum of gain products of all possible two non touching loops)
36
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
−(sum of gain products of all possible three non touching loops)+...
Δi is obtained from Δ by removing the loops which are touching the ith forward path.
Consider the following signal flow graph in order to understand the basic terminology involved
here.
Path
It is a traversal of branches from one node to any other node in the direction of branch arrows. It
should not traverse any node more than once.
Examples − y2→y3→y4→y5 and y5→y3→y2
Forward Path
The path that exists from the input node to the output node is known as forward path.
Examples − y1→y2→y3→y4→y5→y6 and y1→y2→y3→y5→y6.
Forward Path Gain
It is obtained by calculating the product of all branch gains of the forward path.
Examples − abcde is the forward path gain of y1→y2→y3→y4→y5→y6 and abge is the
forward path gain of y1→y2→y3→y5→y6.
Loop
The path that starts from one node and ends at the same node is known as loop. Hence, it is a
closed path.
Examples − y2→y3→y2 and y3→y5→y3.
Loop Gain
It is obtained by calculating the product of all branch gains of a loop.
Examples − bj is the loop gain of y2→y3→y2 and gh is the loop gain of y3→y5→y3.
Non-touching Loops
37
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
These are the loops, which should not have any common node.
Examples − The loops, y2→y3→y2 and y4→y5→y4 are non-touching.
Calculation of Transfer Function using Mason’s Gain Formula
Let us consider the same signal flow graph for finding transfer function.
39
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Here, both the transient and the steady states are indicated in the figure. The responses
corresponding to these states are known as transient and steady state responses.
Mathematically, we can write the time response c(t) as
c(t)=ctr(t)+css(t)
Where,
ctr(t) is the transient response
css(t) is the steady state response
Transient Response
After applying input to the control system, output takes certain time to reach steady state. So,
the output will be in transient state till it goes to a steady state. Therefore, the response of the
control system during the transient state is known as transient response.
The transient response will be zero for large values of ‘t’. Ideally, this value of ‘t’ is infinity and
practically, it is five times constant.
Mathematically, we can write it as
limt→∞ctr(t)=0
Steady state Response
The part of the time response that remains even after the transient response has zero value for
large values of ‘t’ is known as steady state response. This means, the transient response will be
zero even during the steady state.
Example
Let us find the transient and steady state terms of the time response of the control
system c(t)=10+5e−t
40
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Here, the second term 5e−t will be zero as t denotes infinity. So, this is the transient term. And
the first term 10 remains even as t approaches infinity. So, this is the steady state term.
Standard Test Signals
The standard test signals are impulse, step, ramp and parabolic. These signals are used to know
the performance of the control systems using time response of the output.
Unit Impulse Signal
A unit impulse signal, δ(t) is defined as
δ(t)=0 for t≠0
and ∫0+0−δ(t)dt=1
The following figure shows unit impulse signal.
So, the unit impulse signal exists only at ‘t’ is equal to zero. The area of this signal under small
interval of time around ‘t’ is equal to zero is one. The value of unit impulse signal is zero for all
other values of ‘t’.
Unit Step Signal
A unit step signal, u(t) is defined as
u(t)=1;t≥0
=0;t<0
Following figure shows unit step signal.
41
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
So, the unit step signal exists for all positive values of ‘t’ including zero. And its value is one
during this interval. The value of the unit step signal is zero for all negative values of ‘t’.
Unit Ramp Signal
A unit ramp signal, r(t) is defined as
r(t)=t;t≥0
=0;t<0
We can write unit ramp signal, r(t) in terms of unit step signal, u(t) as
r(t)=tu(t)
Following figure shows unit ramp signal.
So, the unit ramp signal exists for all positive values of ‘t’ including zero. And its value
increases linearly with respect to ‘t’ during this interval. The value of unit ramp signal is zero
for all negative values of ‘t’.
Unit Parabolic Signal
A unit parabolic signal, p(t) is defined as,
p(t)=t22;t≥0
=0;t<0
We can write unit parabolic signal, p(t) in terms of the unit step signal, u(t) as,
42
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
p(t)=t22u(t)
The following figure shows the unit parabolic signal.
So, the unit parabolic signal exists for all the positive values of ‘t’ including zero. And its value
increases non-linearly with respect to ‘t’ during this interval. The value of the unit parabolic
signal is zero for all the negative values of ‘t’.
Response of the First Order System
In this chapter, let us discuss the time response of the first order system. Consider the following
block diagram of the closed loop control system. Here, an open loop transfer function, 1sT is
connected with a unity negative feedback.
We know that the transfer function of the closed loop control system has unity negative
feedback as,
C(s)R(s)=G(s)1+G(s)
Substitute, G(s)=1sT in the above equation.
C(s)R(s)=1sT1+1sT=1sT+1
The power of s is one in the denominator term. Hence, the above transfer function is of the first
order and the system is said to be the first order system.
We can re-write the above equation as:
43
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
C(s)=(1sT+1)R(s)
Where,
C(s) is the Laplace transform of the output signal c(t),
R(s) is the Laplace transform of the input signal r(t), and
T is the time constant.
Follow these steps to get the response (output) of the first order system in the time domain.
Take the Laplace transform of the input signal r(t).
Consider the equation, C(s)=(1sT+1)R(s)
Substitute R(s) value in the above equation.
Do partial fractions of C(s) if required.
Apply inverse Laplace transform to C(s).
In the previous chapter, we have seen the standard test signals like impulse, step, ramp and
parabolic. Let us now find out the responses of the first order system for each input, one by one.
The name of the response is given as per the name of the input signal. For example, the
response of the system for an impulse input is called as impulse response.
Impulse Response of First Order System
Consider the unit impulse signal as an input to the first order system.
So, r(t)=δ(t)
Apply Laplace transform on both the sides.
R(s)=1
Consider the equation, C(s)=(1sT+1)R(s)
Substitute, R(s)=1 in the above equation.
C(s)=(1sT+1)(1)=1sT+1
Rearrange the above equation in one of the standard forms of Laplace transforms.
C(s)=1T( s+1T)⇒C(s)=1T(1s+1T)
Apply inverse Laplace transform on both sides.
c(t)=1Te(−tT)u(t)
The unit impulse response is shown in the following figure.
44
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
The unit impulse response, c(t) is an exponential decaying signal for positive values of ‘t’ and it
is zero for negative values of ‘t’.
Step Response of First Order System
Consider the unit step signal as an input to first order system.
So, r(t)=u(t)
Apply Laplace transform on both the sides.
R(s)=1s
Consider the equation, C(s)=(1sT+1)R(s)
Substitute, R(s)=1s in the above equation.
C(s)=(1sT+1)(1s)=1s(sT+1)
Do partial fractions of C(s).
C(s)=1s(sT+1)=As+BsT+1
⇒1s(sT+1)=A(sT+1)+Bss(sT+1)
On both the sides, the denominator term is the same. So, they will get cancelled by each other.
Hence, equate the numerator terms.
1=A(sT+1)+Bs
By equating the constant terms on both the sides, you will get A = 1.
Substitute, A = 1 and equate the coefficient of the s terms on both the sides.
0=T+B⇒B=−T
Substitute, A = 1 and B = −T in partial fraction expansion of C(s).
C(s)=1s−TsT+1=1s−TT(s+1T)
⇒C(s)=1s−1s+1T
Apply inverse Laplace transform on both the sides.
c(t)=(1−e−(tT))u(t)
The unit step response, c(t) has both the transient and the steady state terms.
45
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
The transient term in the unit step response is -
ctr(t)=−e−(tT)u(t)
The steady state term in the unit step response is -
css(t)=u(t)
The following figure shows the unit step response.
The value of the unit step response, c(t) is zero at t = 0 and for all negative values of t. It is
gradually increasing from zero value and finally reaches to one in steady state. So, the steady
state value depends on the magnitude of the input.
Ramp Response of First Order System
Consider the unit ramp signal as an input to the first order system.
So, r(t)=tu(t)
Apply Laplace transform on both the sides.
R(s)=1s2
Consider the equation, C(s)=(1sT+1)R(s)
Substitute, R(s)=1s2 in the above equation.
C(s)=(1sT+1)(1s2)=1s2(sT+1)
Do partial fractions of C(s).
C(s)=1s2(sT+1)=As2+Bs+CsT+1
⇒1s2(sT+1)=A(sT+1)+Bs(sT+1)+Cs2s2(sT+1)
On both the sides, the denominator term is the same. So, they will get cancelled by each other.
Hence, equate the numerator terms.
1=A(sT+1)+Bs(sT+1)+Cs2
By equating the constant terms on both the sides, you will get A = 1.
46
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Substitute, A = 1 and equate the coefficient of the s terms on both the sides.
0=T+B⇒B=−T
Similarly, substitute B = −T and equate the coefficient of s2 terms on both the sides. You will
get C=T2.
Substitute A = 1, B = −T and C=T2 in the partial fraction expansion of C(s).
C(s)=1s2−Ts+T2sT+1=1s2−Ts+T2T(s+1T)
⇒C(s)=1s2−Ts+Ts+1T
Apply inverse Laplace transform on both the sides.
c(t)=(t−T+Te−(tT))u(t)
The unit ramp response, c(t) has both the transient and the steady state terms.
The transient term in the unit ramp response is -
ctr(t)=Te−(tT)u(t)
The steady state term in the unit ramp response is -
css(t)=(t−T)u(t)
The following figure shows the unit ramp response.
The unit ramp response, c(t) follows the unit ramp input signal for all positive values of t. But,
there is a deviation of T units from the input signal.
Parabolic Response of First Order System
Consider the unit parabolic signal as an input to the first order system.
So, r(t)=t22u(t)
Apply Laplace transform on both the sides.
R(s)=1s3
Consider the equation, C(s)=(1sT+1)R(s)
Substitute R(s)=1s3 in the above equation.
47
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
C(s)=(1sT+1)(1s3)=1s3(sT+1)
Do partial fractions of C(s).
C(s)=1s3(sT+1)=As3+Bs2+Cs+DsT+1
After simplifying, you will get the values of A, B, C and D as 1, −T,T2and−T3respectively.
Substitute these values in the above partial fraction expansion of C(s).
C(s)=1s3−Ts2+T2s−T3sT+1⇒C(s)=1s3−Ts2+T2s−T2s+1T
Apply inverse Laplace transform on both the sides.
c(t)=(t22−Tt+T2−T2e−(tT))u(t)
The unit parabolic response, c(t) has both the transient and the steady state terms.
The transient term in the unit parabolic response is
Ctr(t)=−T2e−(tT)u(t)
The steady state term in the unit parabolic response is
Css(t)=(t22−Tt+T2)u(t)
From these responses, we can conclude that the first order control systems are not stable with
the ramp and parabolic inputs because these responses go on increasing even at infinite amount
of time. The first order control systems are stable with impulse and step inputs because these
responses have bounded output. But, the impulse response doesn’t have steady state term. So,
the step signal is widely used in the time domain for analyzing the control systems from their
responses.
We know that the transfer function of the closed loop control system having unity negative
feedback as
C(s)/R(s)=G(s)1+G(s)
48
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Substitute, G(s)=ω2ns(s+2δωn) in the above equation.
C(s)/R(s)=(ω2ns(s+2δωn))1+(ω2ns(s+2δωn))=ω2ns2+2δωns+ω2n
The power of ‘s’ is two in the denominator term. Hence, the above transfer function is of the
second order and the system is said to be the second order system.
The characteristic equation is -
s2+2δωns+ω2n=0
The roots of characteristic equation are -
s=−2ωδn±√(2δωn)2−4ω2n2=−2(δωn±ωn√δ2−1)2
⇒s=−δωn±ωn√δ2−1
The two roots are imaginary when δ = 0.
The two roots are real and equal when δ = 1.
The two roots are real but not equal when δ > 1.
The two roots are complex conjugate when 0 < δ < 1.
We can write C(s) equation as,
C(s)=(ω2ns2+2δωns+ω2n)R(s)
Where,
C(s) is the Laplace transform of the output signal, c(t)
R(s) is the Laplace transform of the input signal, r(t)
ωn is the natural frequency
δ is the damping ratio.
Follow these steps to get the response (output) of the second order system in the time domain.
Take Laplace transform of the input signal, r(t).
Consider the equation, C(s)=(ω2ns2+2δωns+ω2n)R(s)
Substitute R(s) value in the above equation.
Do partial fractions of C(s) if required.
Apply inverse Laplace transform to C(s).
Step Response of Second Order System
Consider the unit step signal as an input to the second order system.
Laplace transform of the unit step signal is,
R(s)=1/s
We know the transfer function of the second order closed loop control system is,
C(s)R(s)=ω2ns2+2δωns+ω2n
Case 1: δ = 0
Substitute, δ=0 in the transfer function.
C(s)R(s)=ω2ns2+ω2n
49
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
⇒C(s)=(ω2ns2+ω2n)R(s)
Substitute, R(s)=1s in the above equation.
C(s)=(ω2ns2+ω2n)(1s)=ω2ns(s2+ω2n)
Apply inverse Laplace transform on both the sides.
c(t)=(1−cos(ωnt))u(t)
So, the unit step response of the second order system when /delta=0 will be a continuous time
signal with constant amplitude and frequency.
Case 2: δ = 1
Substitute, /delta=1 in the transfer function.
C(s)R(s)=ω2ns2+2ωns+ω2n
⇒C(s)=(ω2n(s+ωn)2)R(s)
Substitute, R(s)=1s in the above equation.
C(s)=(ω2n(s+ωn)2)(1s)=ω2ns(s+ωn)2
Do partial fractions of C(s).
C(s)=ω2ns(s+ωn)2=As+Bs+ωn+C(s+ωn)2
After simplifying, you will get the values of A, B and C as 1,−1and−ωnrespectively. Substitute
these values in the above partial fraction expansion of C(s).
C(s)=1s−1s+ωn−ωn(s+ωn)2
Apply inverse Laplace transform on both the sides.
c(t)=(1−e−ωnt−ωnte−ωnt)u(t)
So, the unit step response of the second order system will try to reach the step input in steady
state.
Case 3: 0 < δ < 1
We can modify the denominator term of the transfer function as follows −
s2+2δωns+ω2n={s2+2(s)(δωn)+(δωn)2}+ω2n−(δωn)2
=(s+δωn)2+ω2n(1−δ2)
The transfer function becomes,
C(s)R(s)=ω2n(s+δωn)2+ω2n(1−δ2)
⇒C(s)=(ω2n(s+δωn)2+ω2n(1−δ2))R(s)
Substitute, R(s)=1s in the above equation.
C(s)=(ω2n(s+δωn)2+ω2n(1−δ2))(1s)=ω2ns((s+δωn)2+ω2n(1−δ2))
Do partial fractions of C(s).
C(s)=ω2ns((s+δωn)2+ω2n(1−δ2))=As+Bs+C(s+δωn)2+ω2n(1−δ2)
50
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
After simplifying, you will get the values of A, B and C as 1,−1and−2δωnrespectively.
Substitute these values in the above partial fraction expansion of C(s).
C(s)=1s−s+2δωn(s+δωn)2+ω2n(1−δ2)
C(s)=1s−s+δωn(s+δωn)2+ω2n(1−δ2)−δωn(s+δωn)2+ω2n(1−δ2)
C(s)=1s−(s+δωn)(s+δωn)2+(ωn√1−δ2)2−δ√1−δ2(ωn√1−δ2(s+δωn)2+(ωn√1−δ2)2)
Substitute, ωn√1−δ2 as ωd in the above equation.
C(s)=1s−(s+δωn)(s+δωn)2+ω2d−δ√1−δ2(ωd(s+δωn)2+ω2d)
Apply inverse Laplace transform on both the sides.
c(t)=(1−e−δωntcos(ωdt)−δ√1−δ2e−δωntsin(ωdt))u(t)
c(t)=(1−e−δωnt√1−δ2((√1−δ2)cos(ωdt)+δsin(ωdt)))u(t)
If √1−δ2=sin(θ), then ‘δ’ will be cos(θ). Substitute these values in the above equation.
c(t)=(1−e−δωnt√1−δ2(sin(θ)cos(ωdt)+cos(θ)sin(ωdt)))u(t)
⇒c(t)=(1−(e−δωnt√1−δ2)sin(ωdt+θ))u(t)
So, the unit step response of the second order system is having damped oscillations (decreasing
amplitude) when ‘δ’ lies between zero and one.
Case 4: δ > 1
We can modify the denominator term of the transfer function as follows −
s2+2δωns+ω2n={s2+2(s)(δωn)+(δωn)2}+ω2n−(δωn)2
=(s+δωn)2−ω2n(δ2−1)
The transfer function becomes,
C(s)R(s)=ω2n(s+δωn)2−ω2n(δ2−1)
⇒C(s)=(ω2n(s+δωn)2−ω2n(δ2−1))R(s)
Substitute, R(s)=1s in the above equation.
C(s)=(ω2n(s+δωn)2−(ωn√δ2−1)2)(1s)=ω2ns(s+δωn+ωn√δ2−1)(s+δωn−ωn√δ2−1)
Do partial fractions of C(s).
C(s)=ω2ns(s+δωn+ωn√δ2−1)(s+δωn−ωn√δ2−1)
=As+Bs+δωn+ωn√δ2−1+Cs+δωn−ωn√δ2−1
After simplifying, you will get the values of A, B and C as
1, 12(δ+√δ2−1)(√δ2−1)and −12(δ−√δ2−1)(√δ2−1) respectively. Substitute these values in
above partial fraction expansion of C(s).
C(s)=1s+12(δ+√δ2−1)(√δ2−1)(1s+δωn+ωn√δ2−1)−(12(δ−√δ2−1)(√δ2−1))(1s+δωn−ωn√δ2−1)
Apply inverse Laplace transform on both the sides.
c(t)=(1+(12(δ+√δ2−1)(√δ2−1))e−(δωn+ωn√δ2−1)t−(12(δ−√δ2−1)(√δ2−1))e−(δωn−ωn√δ2−1)t)
u(t)
51
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Since it is over damped, the unit step response of the second order system when δ > 1 will never
reach step input in the steady state.
Impulse Response of Second Order System
The impulse response of the second order system can be obtained by using any one of these two
methods.
Follow the procedure involved while deriving step response by considering the value of R(s) as
1 instead of 1s.
Do the differentiation of the step response.
The following table shows the impulse response of the second order system for 4 cases of the
damping ratio.
Condition of Impulse response for t ≥ 0
Damping ratio
δ=0 ωnsin(ωnt)
δ=1 ω2nte−ωnt
0<δ<1 (ωne−δωnt√1−δ2)sin(ωdt)
δ>1 (ωn2√δ2−1)(e−(δωn−ωn√δ2−1)t−e−(δωn+ωn√δ2−1)t)
Time Domain Specifications
In this chapter, let us discuss the time domain specifications of the second order system. The
step response of the second order system for the underdamped case is shown in the following
figure.
52
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
All the time domain specifications are represented in this figure. The response up to the settling
time is known as transient response and the response after the settling time is known as steady
state response.
Delay Time
It is the time required for the response to reach half of its final value from the zero instant. It is
denoted by td.
Consider the step response of the second order system for t ≥ 0, when ‘δ’ lies between zero and
one.
c(t)=1−(e−δωnt√1−δ2)sin(ωdt+θ)
The final value of the step response is one.
Therefore, at t=td, the value of the step response will be 0.5. Substitute, these values in the
above equation.
c(td)=0.5=1−(e−δωntd√1−δ2)sin(ωdtd+θ)
⇒(e−δωntd√1−δ2)sin(ωdtd+θ)=0.5
By using linear approximation, you will get the delay time td as
td=1+0.7δωn
Rise Time
It is the time required for the response to rise from 0% to 100% of its final value. This is
applicable for the under-damped systems. For the over-damped systems, consider the duration
from 10% to 90% of the final value. Rise time is denoted by tr.
53
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
At t = t1 = 0, c(t) = 0.
We know that the final value of the step response is one.
Therefore, at t=t2, the value of step response is one. Substitute, these values in the following
equation.
c(t)=1−(e−δωnt√1−δ2)sin(ωdt+θ)
c(t2)=1=1−(e−δωnt2√1−δ2)sin(ωdt2+θ)
⇒(e−δωnt2√1−δ2)sin(ωdt2+θ)=0
⇒sin(ωdt2+θ)=0
⇒ωdt2+θ=π
⇒t2=π−θωd
Substitute t1 and t2 values in the following equation of rise time,
tr=t2−t1
∴tr=π−θωd
From above equation, we can conclude that the rise time tr and the damped frequency ωd are
inversely proportional to each other.
Peak Time
It is the time required for the response to reach the peak value for the first time. It is denoted
by tp. At t=tp, the first derivate of the response is zero.
We know the step response of second order system for under-damped case is
c(t)=1−(e−δωnt√1−δ2)sin(ωdt+θ)
Differentiate c(t) with respect to ‘t’.
dc(t)dt=−(e−δωnt√1−δ2)ωdcos(ωdt+θ)−(−δωne−δωnt√1−δ2)sin(ωdt+θ)
Substitute, t=tp and dc(t)dt=0 in the above equation.
0=−(e−δωntp√1−δ2)[ωdcos(ωdtp+θ)−δωnsin(ωdtp+θ)]
⇒ωn√1−δ2cos(ωdtp+θ)−δωnsin(ωdtp+θ)=0
⇒√1−δ2cos(ωdtp+θ)−δsin(ωdtp+θ)=0
⇒sin(θ)cos(ωdtp+θ)−cos(θ)sin(ωdtp+θ)=0
⇒sin(θ−ωdtp−θ)=0
⇒sin(−ωdtp)=0⇒−sin(ωdtp)=0⇒sin(ωdtp)=0
⇒ωdtp=π
⇒tp=πωd
From the above equation, we can conclude that the peak time tp and the damped
frequency ωd are inversely proportional to each other.
Peak Overshoot
54
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Peak overshoot Mp is defined as the deviation of the response at peak time from the final value
of response. It is also called the maximum overshoot.
Mathematically, we can write it as
Mp=c(tp)−c(∞)
Where,
c(tp) is the peak value of the response.
c(∞) is the final (steady state) value of the response.
At t=tp, the response c(t) is -
c(tp)=1−(e−δωntp√1−δ2)sin(ωdtp+θ)
Substitute, tp=πωd in the right hand side of the above equation.
c(tP)=1−(e−δωn(πωd)√1−δ2)sin(ωd(πωd)+θ)
⇒c(tp)=1−(e−(δπ√1−δ2)√1−δ2)(−sin(θ))
We know that
sin(θ)=√1−δ2
So, we will get c(tp) as
c(tp)=1+e−(δπ√1−δ2)
Substitute the values of c(tp) and c(∞) in the peak overshoot equation.
Mp=1+e−(δπ√1−δ2)−1
⇒Mp=e−(δπ√1−δ2)
Percentage of peak overshoot % Mp can be calculated by using this formula.
%Mp=Mpc(∞)×100%
By substituting the values of Mp and c(∞) in above formula, we will get the Percentage of the
peak overshoot %Mp as
%Mp=(e−(δπ√1−δ2))×100%
From the above equation, we can conclude that the percentage of peak overshoot %Mp will
decrease if the damping ratio δ increases.
Settling time
It is the time required for the response to reach the steady state and stay within the specified
tolerance bands around the final value. In general, the tolerance bands are 2% and 5%. The
settling time is denoted by ts.
The settling time for 5% tolerance band is -
ts=3δωn=3τ
The settling time for 2% tolerance band is -
ts=4δωn=4τ
55
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
56
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Where,
R(s) is the Laplace transform of the reference Input signal r(t)
C(s) is the Laplace transform of the output signal c(t)
We know the transfer function of the unity negative feedback closed loop control system as
C(s)R(s)=G(s)1+G(s)
57
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
⇒C(s)=R(s)G(s)1+G(s)
The output of the summing point is -
E(s)=R(s)−C(s)
Substitute C(s) value in the above equation.
E(s)=R(s)−R(s)G(s)1+G(s)
⇒E(s)=R(s)+R(s)G(s)−R(s)G(s)1+G(s)
⇒E(s)=R(s)1+G(s)
Substitute E(s) value in the steady state error formula
ess=lims→0sR(s)1+G(s)
The following table shows the steady state errors and the error constants for standard input
signals like unit step, unit ramp & unit parabolic signals.
Input signal Steady state error ess Error constant
58
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
We can find the steady state errors only for the unity feedback systems. So, we have to convert
the non-unity feedback system into unity feedback system. For this, include one unity positive
feedback path and one unity negative feedback path in the above block diagram. The new block
diagram looks like as shown below.
59
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Simplify the above block diagram by keeping the unity negative feedback as it is. The following
is the simplified block diagram.
This block diagram resembles the block diagram of the unity negative feedback closed loop
control system. Here, the single block is having the transfer
function G(s)1+G(s)H(s)−G(s) instead of G(s). You can now calculate the steady state errors by
using steady state error formula given for the unity negative feedback systems.
Note − It is meaningless to find the steady state errors for unstable closed loop systems. So, we
have to calculate the steady state errors only for closed loop stable systems. This means we need
to check whether the control system is stable or not before finding the steady state errors. In the
next chapter, we will discuss the concepts-related stability.
60
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
This is the response of first order control system for unit step input. This response has the
values between 0 and 1. So, it is bounded output. We know that the unit step signal has the
value of one for all positive values of tincluding zero. So, it is bounded input. Therefore, the
first order control system is stable since both the input and the output are bounded.
Types of Systems based on Stability
We can classify the systems based on stability as follows.
Absolutely stable system
Conditionally stable system
Marginally stable system
Absolutely Stable System
If the system is stable for all the range of system component values, then it is known as
the absolutely stable system. The open loop control system is absolutely stable if all the poles of
the open loop transfer function present in left half of ‘s’ plane. Similarly, the closed loop control
system is absolutely stable if all the poles of the closed loop transfer function present in the left
half of the ‘s’ plane.
61
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
open loop control system is marginally stable if any two poles of the open loop transfer function
is present on the imaginary axis. Similarly, the closed loop control system is marginally stable if
any two poles of the closed loop transfer function is present on the imaginary axis.
So, to overcome this problem there we have the Routh array method. In this method, there is no
need to calculate the roots of the characteristic equation. First formulate the Routh table and
find the number of the sign changes in the first column of the Routh table. The number of sign
changes in the first column of the Routh table gives the number of roots of characteristic
equation that exist in the right half of the ‘s’ plane and the control system is unstable.
Follow this procedure for forming the Routh table.
Fill the first two rows of the Routh array with the coefficients of the characteristic polynomial
as mentioned in the table below. Start with the coefficient of sn and continue up to the
coefficient of s0.
Fill the remaining rows of the Routh array with the elements as mentioned in the table below.
Continue this process till you get the first column element of row s0 is an. Here, an is the
coefficient of s0 in the characteristic polynomial.
Note − If any row elements of the Routh table have some common factor, then you can divide
the row elements with that factor for the simplification will be easy.
The following table shows the Routh array of the nth order characteristic polynomial.
a0sn+a1sn−1+a2sn−2+...+an−1s1+ans0
sn a0 a2 a4 a6 ... ...
⋮ ⋮ ⋮ ⋮
s1 ⋮ ⋮
s0 an
Example
Let us find the stability of the control system having characteristic equation,
s4+3s3+3s2+2s+1=0
Step 1 − Verify the necessary condition for the Routh-Hurwitz stability.
All the coefficients of the characteristic polynomial, s4+3s3+3s2+2s+1 are positive. So, the
control system satisfies the necessary condition.
63
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Step 2 − Form the Routh array for the given characteristic polynomial.
s4 1 3 1
s3 3 2
s2 (3×3)−(2×1)3=73 (3×1)−(0×1)3=33=1
s1 (73×2)−(1×3)73=57
s0 1
Step 3 − Verify the sufficient condition for the Routh-Hurwitz stability.
All the elements of the first column of the Routh array are positive. There is no sign change in
the first column of the Routh array. So, the control system is stable.
Special Cases of Routh Array
We may come across two types of situations, while forming the Routh table. It is difficult to
complete the Routh table from these two situations.
The two special cases are −
The first element of any row of the Routh array is zero.
All the elements of any row of the Routh array are zero.
Let us now discuss how to overcome the difficulty in these two cases, one by one.
First Element of any row of the Routh array is zero
If any row of the Routh array contains only the first element as zero and at least one of the
remaining elements have non-zero value, then replace the first element with a small positive
integer, ϵ. And then continue the process of completing the Routh table. Now, find the number
of sign changes in the first column of the Routh table by substituting ϵ tends to zero.
Example
Let us find the stability of the control system having characteristic equation,
s4+2s3+s2+2s+1=0
Step 1 − Verify the necessary condition for the Routh-Hurwitz stability.
All the coefficients of the characteristic polynomial, s4+2s3+s2+2s+1 are positive. So, the
control system satisfied the necessary condition.
Step 2 − Form the Routh array for the given characteristic polynomial.
s4 1 1 1
s3 21 21
64
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
s2 (1×1)−(1×1)1=0 (1×1)−(0×1)1=1
s1
s0
The row s3 elements have 2 as the common factor. So, all these elements are divided by 2.
Special case (i) − Only the first element of row s2 is zero. So, replace it by ϵand continue the
process of completing the Routh table.
s4 1 1 1
s3 1 1
s2 ϵ 1
s1 (ϵ×1)−(1×1)ϵ=ϵ−1ϵ
s0 1
Step 3 − Verify the sufficient condition for the Routh-Hurwitz stability.
As ϵ tends to zero, the Routh table becomes like this.
s4 1 1 1
s3 1 1
s2 0 1
s1 -∞
s0 1
There are two sign changes in the first column of Routh table. Hence, the control system is
unstable.
All the Elements of any row of the Routh array are zero
In this case, follow these two steps −
Write the auxilary equation, A(s) of the row, which is just above the row of zeros.
Differentiate the auxiliary equation, A(s) with respect to s. Fill the row of zeros with these
coefficients.
65
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Example
Let us find the stability of the control system having characteristic equation,
s5+3s4+s3+3s2+s+3=0
Step 1 − Verify the necessary condition for the Routh-Hurwitz stability.
All the coefficients of the given characteristic polynomial are positive. So, the control system
satisfied the necessary condition.
Step 2 − Form the Routh array for the given characteristic polynomial.
s5 1 1 1
s4 31 31 31
s3 (1×1)−(1×1)1=0 (1×1)−(1×1)1=0
s2
s1
s0
The row s4 elements have the common factor of 3. So, all these elements are divided by 3.
Special case (ii) − All the elements of row s3 are zero. So, write the auxiliary equation, A(s) of
the row s4.
A(s)=s4+s2+1
Differentiate the above equation with respect to s.
dA(s)ds=4s3+2s
Place these coefficients in row s3.
s5 1 1 1
s4 1 1 1
s3 42 21
s2 (2×1)−(1×1)2=0.5 (2×1)−(0×1)2=1
s1 (0.5×1)−(1×2)0.5=−1.50.5=−3
s0 1
66
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
1∞+N(s)D(s)=0⇒N(s)D(s)=0⇒N(s)=0
If K=∞, then N(s)=0. It means the closed loop poles are equal to the open loop zeros when K is
infinity.
From above two cases, we can conclude that the root locus branches start at open loop poles and
end at open loop zeros.
The angle condition is the point at which the angle of the open loop transfer function is an odd
multiple of 1800.
Magnitude of G(s)H(s) is -
|G(s)H(s)|=√(−1)2+02=1
The magnitude condition is that the point (which satisfied the angle condition) at which the
magnitude of the open loop transfer function is one.
We know that the root locus branches start at the open loop poles and end at open loop zeros.
So, the number of root locus branches N is equal to the number of finite open loop poles P or
the number of finite open loop zeros Z, whichever is greater.
Mathematically, we can write the number of root locus branches N as
N=P if P≥Z
N=Z if P<Z
Rule 3 − Identify and draw the real axis root locus branches.
If the angle of the open loop transfer function at a point is an odd multiple of 180 0, then that
point is on the root locus. If odd number of the open loop poles and zeros exist to the left side of
a point on the real axis, then that point is on the root locus branch. Therefore, the branch of
points which satisfies this condition is the real axis of the root locus branch.
Rule 4 − Find the centroid and the angle of asymptotes.
If P=Z, then all the root locus branches start at finite open loop poles and end at finite open loop
zeros.
If P>Z , then Z number of root locus branches start at finite open loop poles and end at finite
open loop zeros and P−Z number of root locus branches start at finite open loop poles and end
at infinite open loop zeros.
If P<Z , then P number of root locus branches start at finite open loop poles and end at finite
open loop zeros and Z−P number of root locus branches start at infinite open loop poles and end
at finite open loop zeros.
So, some of the root locus branches approach infinity, when P≠Z. Asymptotes give the direction
of these root locus branches. The intersection point of asymptotes on the real axis is known
as centroid.
We can calculate the centroid α by using this formula,
α=∑Realpartoffiniteopenlooppoles−∑RealpartoffiniteopenloopzerosP−Z
The formula for the angle of asymptotes θ is
θ=(2q+1)1800P−Z
Where,
q=0,1,2,....,(P−Z)−1
Rule 5 − Find the intersection points of root locus branches with an imaginary axis.
We can calculate the point at which the root locus branch intersects the imaginary axis and the
value of K at that point by using the Routh array method and special case (ii).
If all elements of any row of the Routh array are zero, then the root locus branch intersects the
imaginary axis and vice-versa.
69
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Identify the row in such a way that if we make the first element as zero, then the elements of the
entire row are zero. Find the value of Kfor this combination.
Substitute this K value in the auxiliary equation. You will get the intersection point of the root
locus branch with an imaginary axis.
Rule 6 − Find Break-away and Break-in points.
If there exists a real axis root locus branch between two open loop poles, then there will be
a break-away point in between these two open loop poles.
If there exists a real axis root locus branch between two open loop zeros, then there will be
a break-in point in between these two open loop zeros.
Note − Break-away and break-in points exist only on the real axis root locus branches.
Follow these steps to find break-away and break-in points.
Write K in terms of s from the characteristic equation 1+G(s)H(s)=0.
Differentiate K with respect to s and make it equal to zero. Substitute these values of s in the
above equation.
The values of s for which the K value is positive are the break points.
Rule 7 − Find the angle of departure and the angle of arrival.
The Angle of departure and the angle of arrival can be calculated at complex conjugate open
loop poles and complex conjugate open loop zeros respectively.
The formula for the angle of departure ϕd is
ϕd=1800−ϕ
The formula for the angle of arrival ϕa is
ϕa=1800+ϕ
Where,
ϕ=∑ϕP−∑ϕZ
70
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Example
Let us now draw the root locus of the control system having open loop transfer
function, G(s)H(s)=Ks(s+1)(s+5)
Step 1 − The given open loop transfer function has three poles at s=0,s=−1and s=−5. It doesn’t
have any zero. Therefore, the number of root locus branches is equal to the number of poles of
the open loop transfer function.
N=P=3
The three poles are located are shown in the above figure. The line segment
between s=−1 and s=0 is one branch of root locus on real axis. And the other branch of the root
locus on the real axis is the line segment to the left of s=−5.
Step 2 − We will get the values of the centroid and the angle of asymptotes by using the given
formulae.
Centroid α=−2
The angle of asymptotes are θ=600,1800 and 3000.
The centroid and three asymptotes are shown in the following figure.
71
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Step 3 − Since two asymptotes have the angles of 600 and 3000, two root locus branches
intersect the imaginary axis. By using the Routh array method and special case(ii), the root
locus branches intersects the imaginary axis at j√5 and −j√5.
There will be one break-away point on the real axis root locus branch between the
poles s=−1 and s=0. By following the procedure given for the calculation of break-away point,
we will get it as s=−0.473.
The root locus diagram for the given control system is shown in the following figure.
72
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
In this way, you can draw the root locus diagram of any control system and observe the
movement of poles of the closed loop transfer function.
From the root locus diagrams, we can know the range of K values for different types of
damping.
Effects of Adding Open Loop Poles and Zeros on Root Locus
The root locus can be shifted in ‘s’ plane by adding the open loop poles and the open loop
zeros.
If we include a pole in the open loop transfer function, then some of root locus branches will
move towards right half of ‘s’ plane. Because of this, the damping ratio δ decreases. Which
implies, damped frequency ωd increases and the time domain specifications like delay time td,
rise time tr and peak time tp decrease. But, it effects the system stability.
If we include a zero in the open loop transfer function, then some of root locus branches will
move towards left half of ‘s’ plane. So, it will increase the control system stability. In this case,
the damping ratio δincreases. Which implies, damped frequency ωd decreases and the time
domain specifications like delay time td, rise time tr and peak time tp increase.
So, based on the requirement, we can include (add) the open loop poles or zeros to the transfer
function.
73
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Resonant Frequency
It is the frequency at which the magnitude of the frequency response has peak value for the first
time. It is denoted by ωr. At ω=ωr, the first derivate of the magnitude of T(jω) is zero.
Differentiate M with respect to u.
dMdu=−12[(1−u2)2+(2δu)2]−32[2(1−u2)(−2u)+2(2δu)(2δ)]
⇒dMdu=−12[(1−u2)2+(2δu)2]−32[4u(u2−1+2δ2)]
Substitute, u=ur and dMdu==0 in the above equation.
0=−12[(1−u2r)2+(2δur)2]−32[4ur(u2r−1+2δ2)]
⇒4ur(u2r−1+2δ2)=0
⇒u2r−1+2δ2=0
⇒u2r=1−2δ2
⇒ur=√1−2δ2
Substitute, ur=ωrωn in the above equation.
ωrωn=√1−2δ2
⇒ωr=ωn√1−2δ2
Resonant Peak
It is the peak (maximum) value of the magnitude of T(jω). It is denoted by Mr.
At u=ur, the Magnitude of T(jω) is -
Mr=1√(1−u2r)2+(2δur)2
75
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Substitute, ur=√1−2δ2 and 1−u2r=2δ2 in the above equation.
Mr=1√(2δ2)2+(2δ√1−2δ2)2
⇒Mr=12δ√1−δ2
Resonant peak in frequency response corresponds to the peak overshoot in the time domain
transient response for certain values of damping ratio δ. So, the resonant peak and peak
overshoot are correlated to each other.
Bandwidth
It is the range of frequencies over which, the magnitude of T(jω) drops to 70.7% from its zero
frequency value.
At ω=0, the value of u will be zero.
Substitute, u=0 in M.
M=1√(1−02)2+(2δ(0))2=1
Therefore, the magnitude of T(jω) is one at ω=0.
At 3-dB frequency, the magnitude of T(jω) will be 70.7% of magnitude of T(jω)at ω=0.
i.e., at ω=ωB,M=0.707(1)=1√2
⇒M=1√2=1√(1−u2b)2+(2δub)2
⇒2=(1−u2b)2+(2δ)2u2b
Let, u2b=x
⇒2=(1−x)2+(2δ)2x
⇒x2+(4δ2−2)x−1=0
⇒x=−(4δ2−2)±√(4δ2−2)2+42
Consider only the positive value of x.
x=1−2δ2+√(2δ2−1)2+1
⇒x=1−2δ2+√(2−4δ2+4δ4)
Substitute, x=u2b=ω2bω2n
ω2bω2n=1−2δ2+√(2−4δ2+4δ4)
⇒ωb=ωn√1−2δ2+√(2−4δ2+4δ4)
Bandwidth ωb in the frequency response is inversely proportional to the rise time tr in the time
domain transient response.
76
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
In both the plots, x-axis represents angular frequency (logarithmic scale). Whereas, yaxis
represents the magnitude (linear scale) of open loop transfer function in the magnitude plot and
the phase angle (linear scale) of the open loop transfer function in the phase plot.
The magnitude of the open loop transfer function in dB is -
M=20log|G(jω)H(jω)|
The phase angle of the open loop transfer function in degrees is -
ϕ=∠G(jω)H(jω)
Note − The base of logarithm is 10.
Basic of Bode Plots
The following table shows the slope, magnitude and the phase angle values of the terms present
in the open loop transfer function. This data is useful while drawing the Bode plots.
Type of term G(jω)H(jω) Slope(dB/dec) Magnitude (dB) Phase angle(degrees)
Constant K 0 20logK 0
The magnitude plot is a horizontal line, which is independent of frequency. The 0 dB line itself
is the magnitude plot when the value of K is one. For the positive values of K, the horizontal
line will shift 20logK dB above the 0 dB line. For the negative values of K, the horizontal line
will shift 20logK dB below the 0 dB line. The Zero degrees line itself is the phase plot for all
the positive values of K.
Consider the open loop transfer function G(s)H(s)=s.
Magnitude M=20logω dB
Phase angle ϕ=900
At ω=0.1 rad/sec, the magnitude is -20 dB.
At ω=1 rad/sec, the magnitude is 0 dB.
At ω=10 rad/sec, the magnitude is 20 dB.
78
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
The following figure shows the corresponding Bode plot.
The magnitude plot is a line, which is having a slope of 20 dB/dec. This line started
at ω=0.1 rad/sec having a magnitude of -20 dB and it continues on the same slope. It is touching
0 dB line at ω=1 rad/sec. In this case, the phase plot is 900 line.
Consider the open loop transfer function G(s)H(s)=1+sτ.
Magnitude M=20log√1+ω2τ2 dB
Phase angle ϕ=tan−1ωτ degrees
For ω<1τ , the magnitude is 0 dB and phase angle is 0 degrees.
For ω>1τ , the magnitude is 20logωτ dB and phase angle is 900.
The following figure shows the corresponding Bode plot.
79
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
The magnitude plot is having magnitude of 0 dB upto ω=1τ rad/sec. From ω=1τ rad/sec, it is
having a slope of 20 dB/dec. In this case, the phase plot is having phase angle of 0 degrees up
to ω=1τ rad/sec and from here, it is having phase angle of 900. This Bode plot is called
the asymptotic Bode plot.
As the magnitude and the phase plots are represented with straight lines, the Exact Bode plots
resemble the asymptotic Bode plots. The only difference is that the Exact Bode plots will have
simple curves instead of straight lines.
Similarly, you can draw the Bode plots for other terms of the open loop transfer function which
are given in the table.
Control Systems - Construction of Bode Plots
In this chapter, let us understand in detail how to construct (draw) Bode plots.
Rules for Construction of Bode Plots
Follow these rules while constructing a Bode plot.
Represent the open loop transfer function in the standard time constant form.
Substitute, s=jω in the above equation.
Find the corner frequencies and arrange them in ascending order.
80
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
th
Consider the starting frequency of the Bode plot as 1/10 of the minimum corner frequency or
0.1 rad/sec whichever is smaller value and draw the Bode plot up to 10 times maximum corner
frequency.
Draw the magnitude plots for each term and combine these plots properly.
Draw the phase plots for each term and combine these plots properly.
Note − The corner frequency is the frequency at which there is a change in the slope of the
magnitude plot.
Example: Consider the open loop transfer function of a closed loop control system
G(s)H(s)=10s(s+2)(s+5)
Let us convert this open loop transfer function into standard time constant form.
G(s)H(s)=10s2(s2+1)5(s5+1)
⇒G(s)H(s)=s(1+s2)(1+s5)
So, we can draw the Bode plot in semi log sheet using the rules mentioned earlier.
Stability Analysis using Bode Plots
From the Bode plots, we can say whether the control system is stable, marginally stable or
unstable based on the values of these parameters.
Gain cross over frequency and phase cross over frequency
Gain margin and phase margin
Phase Cross over Frequency
The frequency at which the phase plot is having the phase of -1800 is known as phase cross over
frequency. It is denoted by ωpc. The unit of phase cross over frequency is rad/sec.
Gain Cross over Frequency
The frequency at which the magnitude plot is having the magnitude of zero dB is known as gain
cross over frequency. It is denoted by ωgc. The unit of gain cross over frequency is rad/sec.
The stability of the control system based on the relation between the phase cross over frequency
and the gain cross over frequency is listed below.
If the phase cross over frequency ωpc is greater than the gain cross over frequency ωgc, then
the control system is stable.
If the phase cross over frequency ωpc is equal to the gain cross over frequency ωgc, then the
control system is marginally stable.
If the phase cross over frequency ωpc is less than the gain cross over frequency ωgc, then the
control system is unstable.
Gain Margin=== Gain margin GM is equal to negative of the magnitude in dB at phase cross
over frequency.
GM=20log(1Mpc)=20logMpc
81
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Where, Mpc is the magnitude at phase cross over frequency. The unit of gain margin (GM)
is dB.
Phase Margin The formula for phase margin PM is
PM=1800+ϕgc
Where, ϕgc is the phase angle at gain cross over frequency. The unit of phase margin is degrees.
The stability of the control system based on the relation between gain margin and phase margin
is listed below.If both the gain margin GM and the phase margin PM are positive, then the
control system is stable. If both the gain margin GM and the phase margin PM are equal to zero,
then the control system is marginally stable.
If the gain margin GM and / or the phase margin PM are/is negative, then the control system
is unstable.
This graph sheet consists of concentric circles and radial lines. The concentric circles and
the radial lines represent the magnitudes and phase angles respectively. These angles are
82
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
represented by positive values in anti-clock wise direction. Similarly, we can represent angles
with negative values in clockwise direction. For example, the angle 2700 in anti-clock wise
direction is equal to the angle −900 in clockwise direction.
Rules for Drawing Polar Plots
Follow these rules for plotting the polar plots.
Substitute, s=jω in the open loop transfer function.
Write the expressions for magnitude and the phase of G(jω)H(jω).
Find the starting magnitude and the phase of G(jω)H(jω) by substituting ω=0. So, the polar plot
starts with this magnitude and the phase angle.
Find the ending magnitude and the phase of G(jω)H(jω) by substituting ω=∞. So, the polar plot
ends with this magnitude and the phase angle.
Check whether the polar plot intersects the real axis, by making the imaginary term
of G(jω)H(jω) equal to zero and find the value(s) of ω.
Check whether the polar plot intersects the imaginary axis, by making real term
of G(jω)H(jω) equal to zero and find the value(s) of ω.
For drawing polar plot more clearly, find the magnitude and phase of G(jω)H(jω) by
considering the other value(s) of ω.
Example
Consider the open loop transfer function of a closed loop control system.
G(s)H(s)=5s(s+1)(s+2)
Let us draw the polar plot for this control system using the above rules.
Step 1 − Substitute, s=jω in the open loop transfer function.
G(jω)H(jω)=5jω(jω+1)(jω+2)
The magnitude of the open loop transfer function is
M=5ω(√ω2+1)(√ω2+4)
The phase angle of the open loop transfer function is
ϕ=−900−tan−1ω−tan−1ω2
Step 2 − The following table shows the magnitude and the phase angle of the open loop transfer
function at ω=0 rad/sec and ω=∞ rad/sec.
Frequency (rad/sec) Magnitude Phase angle(degrees)
0 ∞ -90 or 270
∞ 0 -270 or 90
83
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
0 0
So, the polar plot starts at (∞,−90 ) and ends at (0,−270 ). The first and the second terms within
the brackets indicate the magnitude and phase angle respectively.
Step 3 − Based on the starting and the ending polar co-ordinates, this polar plot will intersect
the negative real axis. The phase angle corresponding to the negative real axis is −180 0 or 1800.
So, by equating the phase angle of the open loop transfer function to either −180 0 or 1800, we
will get the ω value as √2.
By substituting ω=√2 in the magnitude of the open loop transfer function, we will get M=0.83.
Therefore, the polar plot intersects the negative real axis when ω=√2 and the polar coordinate is
(0.83,−1800).
So, we can draw the polar plot with the above information on the polar graph sheet.
Control Systems - Nyquist Plots
Nyquist plots are the continuation of polar plots for finding the stability of the closed loop
control systems by varying ω from −∞ to ∞. That means, Nyquist plots are used to draw the
complete frequency response of the open loop transfer function.
Nyquist Stability Criterion
The Nyquist stability criterion works on the principle of argument. It states that if there are P
poles and Z zeros are enclosed by the ‘s’ plane closed path, then the
corresponding G(s)H(s) plane must encircle the origin P−Z times. So, we can write the number
of encirclements N as,
N=P−Z
If the enclosed ‘s’ plane closed path contains only poles, then the direction of the encirclement
in the G(s)H(s) plane will be opposite to the direction of the enclosed closed path in the ‘s’
plane.
If the enclosed ‘s’ plane closed path contains only zeros, then the direction of the encirclement
in the G(s)H(s) plane will be in the same direction as that of the enclosed closed path in the ‘s’
plane.
Let us now apply the principle of argument to the entire right half of the ‘s’ plane by selecting it
as a closed path. This selected path is called the Nyquist contour.
We know that the closed loop control system is stable if all the poles of the closed loop transfer
function are in the left half of the ‘s’ plane. So, the poles of the closed loop transfer function are
nothing but the roots of the characteristic equation. As the order of the characteristic equation
increases, it is difficult to find the roots. So, let us correlate these roots of the characteristic
equation as follows.
The Poles of the characteristic equation are same as that of the poles of the open loop transfer
function.
84
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
The zeros of the characteristic equation are same as that of the poles of the closed loop transfer
function.
We know that the open loop control system is stable if there is no open loop pole in the the right
half of the ‘s’ plane.
i.e.,P=0⇒N=−Z
We know that the closed loop control system is stable if there is no closed loop pole in the right
half of the ‘s’ plane.
i.e.,Z=0⇒N=P
Nyquist stability criterion states the number of encirclements about the critical point (1+j0)
must be equal to the poles of characteristic equation, which is nothing but the poles of the open
loop transfer function in the right half of the ‘s’ plane. The shift in origin to (1+j0) gives the
characteristic equation plane.
Rules for Drawing Nyquist Plots
Follow these rules for plotting the Nyquist plots.
Locate the poles and zeros of open loop transfer function G(s)H(s) in ‘s’ plane.
Draw the polar plot by varying ω from zero to infinity. If pole or zero present at s = 0, then
varying ω from 0+ to infinity for drawing polar plot.
Draw the mirror image of above polar plot for values of ω ranging from −∞ to zero (0− if any
pole or zero present at s=0).
The number of infinite radius half circles will be equal to the number of poles or zeros at origin.
The infinite radius half circle will start at the point where the mirror image of the polar plot
ends. And this infinite radius half circle will end at the point where the polar plot starts.
After drawing the Nyquist plot, we can find the stability of the closed loop control system using
the Nyquist stability criterion. If the critical point (-1+j0) lies outside the encirclement, then the
closed loop control system is absolutely stable.
Stability Analysis using Nyquist Plots
From the Nyquist plots, we can identify whether the control system is stable, marginally stable
or unstable based on the values of these parameters.
Gain cross over frequency and phase cross over frequency
Gain margin and phase margin
Phase Cross over Frequency
The frequency at which the Nyquist plot intersects the negative real axis (phase angle is 180 0) is
known as the phase cross over frequency. It is denoted by ωpc.
Gain Cross over Frequency
85
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
The frequency at which the Nyquist plot is having the magnitude of one is known as the gain
cross over frequency. It is denoted by ωgc.
The stability of the control system based on the relation between phase cross over frequency
and gain cross over frequency is listed below.
If the phase cross over frequency ωpc is greater than the gain cross over frequency ωgc, then
the control system is stable.
If the phase cross over frequency ωpc is equal to the gain cross over frequency ωgc, then the
control system is marginally stable.
If phase cross over frequency ωpc is less than gain cross over frequency ωgc, then the control
system is unstable.
Gain Margin
The gain margin GM is equal to the reciprocal of the magnitude of the Nyquist plot at the phase
cross over frequency.
GM=1Mpc
Where, Mpc is the magnitude in normal scale at the phase cross over frequency.
Phase Margin
The phase margin PM is equal to the sum of 1800 and the phase angle at the gain cross over
frequency.
PM=1800+ϕgc
Where, ϕgc is the phase angle at the gain cross over frequency.
The stability of the control system based on the relation between the gain margin and the phase
margin is listed below.
If the gain margin GM is greater than one and the phase margin PM is positive, then the control
system is stable.
If the gain margin GM is equal to one and the phase margin PM is zero degrees, then the control
system is marginally stable.
If the gain margin GM is less than one and / or the phase margin PM is negative, then the
control system is unstable.
Control Systems - Compensators
There are three types of compensators — lag, lead and lag-lead compensators. These are most
commonly used.
Lag Compensator
The Lag Compensator is an electrical network which produces a sinusoidal output having the
phase lag when a sinusoidal input is applied. The lag compensator circuit in the ‘s’ domain is
shown in the following figure.
86
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Here, the capacitor is in series with the resistor R2 and the output is measured across this
combination.
The transfer function of this lag compensator is -
Vo(s)Vi(s)=1α(s+1τs+1ατ)
Where,
τ=R2C
α=R1+R2R2
From the above equation, α is always greater than one.
From the transfer function, we can conclude that the lag compensator has one pole
at s=−1ατ and one zero at s=−1τ . This means, the pole will be nearer to origin in the pole-zero
configuration of the lag compensator.
Substitute, s=jω in the transfer function.
Vo(jω)Vi(jω)=1α(jω+1τjω+1ατ)
Phase angle ϕ=tan−1ωτ−tan−1αωτ
We know that, the phase of the output sinusoidal signal is equal to the sum of the phase angles
of input sinusoidal signal and the transfer function.
So, in order to produce the phase lag at the output of this compensator, the phase angle of the
transfer function should be negative. This will happen when α>1.
Lead Compensator
The lead compensator is an electrical network which produces a sinusoidal output having phase
lead when a sinusoidal input is applied. The lead compensator circuit in the ‘s’ domain is shown
in the following figure.
87
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Here, the capacitor is parallel to the resistor R1 and the output is measured across resistor $R_2.
The transfer function of this lead compensator is -
Vo(s)Vi(s)=β(sτ+1βsτ+1)
Where,
τ=R1C
β=R2R1+R2
From the transfer function, we can conclude that the lead compensator has pole at s=−1β and
zero at s=−1βτ.
Substitute, s=jω in the transfer function.
Vo(jω)Vi(jω)=β(jωτ+1βjωτ+1)
Phase angle ϕ=tan−1ωτ−tan−1βωτ
We know that, the phase of the output sinusoidal signal is equal to the sum of the phase angles
of input sinusoidal signal and the transfer function.
So, in order to produce the phase lead at the output of this compensator, the phase angle of the
transfer function should be positive. This will happen when 0<β<1. Therefore, zero will be
nearer to origin in pole-zero configuration of the lead compensator.
Lag-Lead Compensator
Lag-Lead compensator is an electrical network which produces phase lag at one frequency
region and phase lead at other frequency region. It is a combination of both the lag and the lead
compensators. The lag-lead compensator circuit in the ‘s’ domain is shown in the following
figure.
88
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
This circuit looks like both the compensators are cascaded. So, the transfer function of this
circuit will be the product of transfer functions of the lead and the lag compensators.
Vo(s)Vi(s)=β(sτ1+1βsτ1+1)1α(s+1τ2s+1ατ2)
We know αβ=1.
⇒Vo(s)Vi(s)=(s+1τ1s+1βτ1)(s+1τ2s+1ατ2)
Where,
τ1=R1C1
τ2=R2C2
Control Systems - Controllers
The various types of controllers are used to improve the performance of control systems. In this
chapter, we will discuss the basic controllers such as the proportional, the derivative and the
integral controllers.
Proportional Controller
The proportional controller produces an output, which is proportional to error signal.
u(t)∝e(t)
⇒u(t)=KPe(t)
Apply Laplace transform on both the sides -
U(s)=KPE(s)
U(s)E(s)=KP
Therefore, the transfer function of the proportional controller is KP.
Where,
U(s) is the Laplace transform of the actuating signal u(t)
89
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
E(s) is the Laplace transform of the error signal e(t)
KP is the proportionality constant
The block diagram of the unity negative feedback closed loop control system along with the
proportional controller is shown in the following figure.
The proportional controller is used to change the transient response as per the requirement.
Derivative Controller
The derivative controller produces an output, which is derivative of the error signal.
u(t)=KDde(t)dt
Apply Laplace transform on both sides.
U(s)=KDsE(s)
U(s)E(s)=KDs
Therefore, the transfer function of the derivative controller is KDs.
Where, KD is the derivative constant.
The block diagram of the unity negative feedback closed loop control system along with the
derivative controller is shown in the following figure.
The derivative controller is used to make the unstable control system into a stable one.
Integral Controller
The integral controller produces an output, which is integral of the error signal.
u(t)=KI∫e(t)dt
90
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Apply Laplace transform on both the sides -
U(s)=KIE(s)s
U(s)E(s)=KIs
Therefore, the transfer function of the integral controller is KIs.
Where, KI is the integral constant.
The block diagram of the unity negative feedback closed loop control system along with the
integral controller is shown in the following figure.
91
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
The proportional derivative controller is used to improve the stability of control system without
affecting the steady state error.
Proportional Integral (PI) Controller
The proportional integral controller produces an output, which is the combination of outputs of
the proportional and integral controllers.
u(t)=KPe(t)+KI∫e(t)dt
Apply Laplace transform on both sides -
U(s)=(KP+KIs)E(s)
U(s)E(s)=KP+KIs
Therefore, the transfer function of proportional integral controller is KP+KIs.
The block diagram of the unity negative feedback closed loop control system along with the
proportional integral controller is shown in the following figure.
The proportional integral controller is used to decrease the steady state error without affecting
the stability of the control system.
Proportional Integral Derivative (PID) Controller
The proportional integral derivative controller produces an output, which is the combination of
the outputs of proportional, integral and derivative controllers.
u(t)=KPe(t)+KI∫e(t)dt+KDde(t)/dt
Apply Laplace transform on both sides -
U(s)=(KP+KIs+KDs)E(s)
U(s)E(s)=KP+KIs+KDs
Therefore, the transfer function of the proportional integral derivative controller
is KP+KIs+KDs.
The block diagram of the unity negative feedback closed loop control system along with the
proportional integral derivative controller is shown in the following figure.
92
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
The proportional integral derivative controller is used to improve the stability of the control
system and to decrease steady state error.
Control Systems - State Space Model
The state space model of Linear Time-Invariant (LTI) system can be represented as,
˙X=AX+BU
Y=CX+DU
The first and the second equations are known as state equation and output equation respectively.
Where,
X and ˙X are the state vector and the differential state vector respectively.
U and Y are input vector and output vector respectively.
A is the system matrix.
B and C are the input and the output matrices.
D is the feed-forward matrix.
Basic Concepts of State Space Model
The following basic terminology involved in this chapter.
State
It is a group of variables, which summarizes the history of the system in order to predict the
future values (outputs).
State Variable
The number of the state variables required is equal to the number of the storage elements
present in the system.
Examples − current flowing through inductor, voltage across capacitor
State Vector
It is a vector, which contains the state variables as elements.
In the earlier chapters, we have discussed two mathematical models of the control systems.
Those are the differential equation model and the transfer function model. The state space
93
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
model can be obtained from any one of these two mathematical models. Let us now discuss
these two methods one by one.
State Space Model from Differential Equation
Consider the following series of the RLC circuit. It is having an input voltage, vi(t) and the
current flowing through the circuit is i(t).
There are two storage elements (inductor and capacitor) in this circuit. So, the number of the
state variables is equal to two and these state variables are the current flowing through the
inductor, i(t) and the voltage across capacitor, vc(t).
From the circuit, the output voltage, v0(t) is equal to the voltage across capacitor, vc(t).
v0(t)=vc(t)
Apply KVL around the loop.
vi(t)=Ri(t)+Ldi(t)/dt+vc(t)
⇒di(t)/dt=−Ri(t)L−vc(t)L+vi(t)L
The voltage across the capacitor is -
vc(t)=1C∫i(t)dt
Differentiate the above equation with respect to time.
dvc(t)/dt=i(t)C
State vector, X=[i(t)vc(t)]
Differential state vector, ˙X=[di(t)/dt dvc(t)dt]
We can arrange the differential equations and output equation into the standard form of state
space model as,
˙X=[di(t)/dt dvc(t)/dt]=[−RL−1L1C0][i(t)vc(t)]+[1L0][vi(t)]
Y=[01][i(t)vc(t)]
Where,
A=[−RL−1L1C0],B=[1L0],C=[01]and D=[0]
State Space Model from Transfer Function
94
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Consider the two types of transfer functions based on the type of terms present in the
numerator.
Transfer function having constant term in Numerator.
Transfer function having polynomial function of ‘s’ in Numerator.
Transfer function having constant term in Numerator
Consider the following transfer function of a system
Y(s)U(s)=b0sn+an−1sn−1+...+a1s+a0
Rearrange, the above equation as
(sn+an−1sn−1+...+a0)Y(s)=b0U(s)
Apply inverse Laplace transform on both sides.
dny(t)dtn+an−1dn−1y(t)dtn−1+...+a1dy(t)dt+a0y(t)=b0u(t)
Let
y(t)=x1
dy(t)dt=x2=˙x1
d2y(t)dt2=x3=˙x2
.
.
.
dn−1y(t)dtn−1=xn=˙xn−1
dny(t)dtn=˙xn
and u(t)=u
Then,
˙xn+an−1xn+...+a1x2+a0x1=b0u
From the above equation, we can write the following state equation.
˙xn=−a0x1−a1x2−...−an−1xn+b0u
The output equation is -
y(t)=y=x1
The state space model is -
˙X=[˙x1˙x2⋮˙xn−1˙xn]
=[010…00001…00⋮⋮⋮…⋮⋮000…01−a0−a1−a2…−an−2−an−1][x1x2⋮xn−1xn]+[00⋮0b0][u]
Y=[10…00][x1x2⋮xn−1xn]
Here, D=[0].
Example
Find the state space model for the system having transfer function.
95
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Y(s)U(s)=1s2+s+1
Rearrange, the above equation as,
(s2+s+1)Y(s)=U(s)
Apply inverse Laplace transform on both the sides.
d2y(t)dt2+dy(t)dt+y(t)=u(t)
Let
y(t)=x1
dy(t)dt=x2=˙x1
and u(t)=u
Then, the state equation is
˙x2=−x1−x2+u
The output equation is
y(t)=y=x1
The state space model is
˙X=[˙x1˙x2]=[01−1−1][x1x2]+[01][u]
Y=[10][x1x2]
Transfer function having polynomial function of ‘s’ in Numerator
Consider the following transfer function of a system
Y(s)U(s)=bnsn+bn−1sn−1+...+b1s+b0sn+an−1sn−1+...+a1s+a0
⇒Y(s)U(s)=(1sn+an−1sn−1+...+a1s+a0)(bnsn+bn−1sn−1+...+b1s+b0)
The above equation is in the form of product of transfer functions of two blocks, which are
cascaded.
Y(s)U(s)=(V(s)U(s))(Y(s)V(s))
Here,
V(s)U(s)=1sn+an−1sn−1+...+a1s+a0
Rearrange, the above equation as
(sn+an−1sn−1+...+a0)V(s)=U(s)
Apply inverse Laplace transform on both the sides.
dnv(t)/dtn+an−1dn−1v(t)dtn−1+...+a1dv(t)dt+a0v(t)=u(t)
Let
v(t)=x1
dv((t)dt=x2=˙x1
d2v(t)dt2=x3=˙x2
.
.
96
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
.
dn−1v(t)dtn−1=xn=˙xn−1
dnv(t)dtn=˙xn
and u(t)=u
Then, the state equation is
˙xn=−a0x1−a1x2−...−an−1xn+u
Consider,
Y(s)V(s)=bnsn+bn−1sn−1+...+b1s+b0
Rearrange, the above equation as
Y(s)=(bnsn+bn−1sn−1+...+b1s+b0)V(s)
Apply inverse Laplace transform on both the sides.
y(t)=bndnv(t)dtn+bn−1dn−1v(t)dtn−1+...+b1dv(t)dt+b0v(t)
By substituting the state variables and y(t)=y in the above equation, will get the output equation
as,
y=bn˙xn+bn−1xn+...+b1x2+b0x1
Substitute, ˙xn value in the above equation.
y=bn(−a0x1−a1x2−...−an−1xn+u)+bn−1xn+...+b1x2+b0x1
y=(b0−bna0)x1+(b1−bna1)x2+...+(bn−1−bnan−1)xn+bnu
The state space model is
˙X=[˙x1˙x2⋮˙xn−1˙xn]
=[010…00001…00⋮⋮⋮…⋮⋮000…01−a0−a1−a2…−an−2−an−1][x1x2⋮xn−1xn]+[00⋮0b0][u]
Y=[b0−bna0b1−bna1...bn−2−bnan−2bn−1−bnan−1][x1x2⋮xn−1xn]
If bn=0, then,
Y=[b0b1...bn−2bn−1][x1x2⋮xn−1xn]
Control Systems - State Space Analysis
In the previous chapter, we learnt how to obtain the state space model from differential equation
and transfer function. In this chapter, let us discuss how to obtain transfer function from the
state space model.
Transfer Function from State Space Model
We know the state space model of a Linear Time-Invariant (LTI) system is -
˙X=AX+BU
Y=CX+DU
Apply Laplace Transform on both sides of the state equation.
sX(s)=AX(s)+BU(s)
⇒(sI−A)X(s)=BU(s)
97
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
⇒X(s)=(sI−A)−1BU(s)
Apply Laplace Transform on both sides of the output equation.
Y(s)=CX(s)+DU(s)
Substitute, X(s) value in the above equation.
⇒Y(s)=C(sI−A)−1BU(s)+DU(s)
⇒Y(s)=[C(sI−A)−1B+D]U(s)
⇒Y(s)U(s)=C(sI−A)−1B+D
The above equation represents the transfer function of the system. So, we can calculate the
transfer function of the system by using this formula for the system represented in the state
space model.
Note − When D=[0], the transfer function will be
Y(s)U(s)=C(sI−A)−1B
Example
Let us calculate the transfer function of the system represented in the state space model as,
˙X=[˙x1˙x2]=[−1−110][x1x2]+[10][u]
Y=[01][x1x2]
Here,
A=[−1−110],B=[10],C=[01]and D=[0]
The formula for the transfer function when D=[0] is -
Y(s)U(s)=C(sI−A)−1B
Substitute, A, B & C matrices in the above equation.
Y(s)U(s)=[01][s+11−1s]−1[10]
⇒Y(s)U(s)=[01][s−11s+1](s+1)s−1(−1)[10]
⇒Y(s)U(s)=[01][s1]s2+s+1=1s2+s+1
Therefore, the transfer function of the system for the given state space model is
Y(s)U(s)=1s2+s+1
State Transition Matrix and its Properties
If the system is having initial conditions, then it will produce an output. Since, this output is
present even in the absence of input, it is called zero input response x ZI R(t). Mathematically,
we can write it as,
xZIR(t)=eAtX(0)=L−1{[sI−A]−1X(0)}
From the above relation, we can write the state transition matrix ϕ(t) as
ϕ(t)=eAt=L−1[sI−A]−1
98
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
So, the zero input response can be obtained by multiplying the state transition matrix ϕ(t) with
the initial conditions matrix.
Following are the properties of the state transition matrix.
If t=0, then state transition matrix will be equal to an Identity matrix.
ϕ(0)=I
Inverse of state transition matrix will be same as that of state transition matrix just by replcing
‘t’ by ‘-t’.
ϕ−1(t)=ϕ(−t)
If t=t1+t2 , then the corresponding state transition matrix is equal to the multiplication of the
two state transition matrices at t=t1 and t=t2.
ϕ(t1+t2)=ϕ(t1)ϕ(t2)
Controllability and Observability
Let us now discuss controllability and observability of control system one by one.
Controllability
A control system is said to be controllable if the initial states of the control system are
transferred (changed) to some other desired states by a controlled input in finite duration of
time.
We can check the controllability of a control system by using Kalman’s test.
Write the matrix Qc in the following form.
Qc=[BABA2B...An−1B]
Find the determinant of matrix Qc and if it is not equal to zero, then the control system is
controllable.
Observability
A control system is said to be observable if it is able to determine the initial states of the control
system by observing the outputs in finite duration of time.
We can check the observability of a control system by using Kalman’s test.
Write the matrix Qo in following form.
Qo=[CTATCT(AT)2CT...(AT)n−1CT]
Find the determinant of matrix Qo and if it is not equal to zero, then the control system is
observable.
Example
99
زياد سعيد دمحم السليم. د: مدرش المادة- المرحلت الثالثت –قسم هندست تقنياث الحاسوبControl Engineering Fundamentals
-0202-0202 قسم هندست تقنياث االجهسة الطبيت- المرحلت الرابعتControl Systems الموصل/الكليت التقنيت الهندسيت
Let us verify the controllability and observability of a control system which is represented in the
state space model as,
˙x=[˙x1˙x2]=[−1−110][x1x2]+[10][u]
Y=[01][x1x2]
Here,
A=[−1−110],B=[10],[01],D=[0] and n=2
For n=2, the matrix Qc will be
Qc=[BAB]
We will get the product of matrices A and B as,
AB=[−11]
⇒Qc=[1−101]
|Qc|=1≠0
Since the determinant of matrix Qc is not equal to zero, the given control system is controllable.
For n=2, the matrix Qo will be -
Qo=[CTATCT]
Here,
AT=[−11−10] and CT=[01]
We will get the product of matrices AT and CT as
ATCT=[10]
⇒Qo=[0110]
⇒|Qo|=−1≠0
Since, the determinant of matrix Qo is not equal to zero, the given control system is observable.
Therefore, the given control system is both controllable and observable.
100