Control Systems Quick Guide
Control Systems Quick Guide
Here, the control system is represented by a single block. Since, the output is controlled by varying input, the control system got this name. We will vary
this input with some mechanism. In the next section on open loop and closed loop control systems, we will study in detail about the blocks inside the
control system and how to vary this input in order to get the desired response.
Traffic lights control system is an example of control system. Here, a sequence of input signal is applied to this control system and the output is one
of the three lights that will be on for some duration of time. During this time, the other two lights will be off. Based on the traffic study at a particular
junction, the on and off times of the lights can be determined. Accordingly, the input signal controls the output. So, the traffic lights control system
operates on time basis.
Control Systems can be classified as continuous time control systems and discrete time control systems based on the type of the signal used.
In continuous time control systems, all the signals are continuous in time. But, in discrete time control systems, there exists one or more
discrete time signals.
Control Systems can be classified as SISO control systems and MIMO control systems based on the number of inputs and outputs present.
SISO SingleInputandSingleOutput control systems have one input and one output. Whereas, MIMO MultipleInputsandMultipleOutputs control
systems have more than one input and more than one output.
In open loop control systems, output is not fed-back to the input. So, the control action is independent of the desired output.
The following figure shows the block diagram of the open loop control system.
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.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 1/90
17/01/2018 Control Systems Quick Guide
In closed loop control systems, output is fed back 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 feedbackelements by 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.
Control action is independent of the desired output. Control action is dependent of the desired output.
These are also called as non-feedback control systems. These are also called as feedback control systems.
Inaccurate. Accurate.
Types of Feedback
There are two types of feedback −
Positive feedback
Negative feedback
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.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 2/90
17/01/2018 Control Systems Quick Guide
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,
G
T= 1 − GH
Equation1
Where,
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.
Where,
Effects of Feedback
Let us now understand the effects of feedback.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 3/90
17/01/2018 Control Systems Quick Guide
From Equation 2, we can say that the overall gain of negative feedback closed loop control system is the ratio of 'G' and 1 + GH. So, the overall gain
may increase or decrease depending on the value of 1 + GH.
If the value of 1 + GH is less than 1, then the overall gain increases. In this case, 'GH' value is negative because the gain of the feedback path is
negative.
If the value of 1 + GH is greater than 1, then the overall gain decreases. In this case, 'GH' value is positive because the gain of the feedback path is
positive.
In general, 'G' and 'H' are functions of frequency. So, the feedback will increase the overall gain of the system in one frequency range and decrease in the
other frequency range.
Sensitivity of the overall gain of negative feedback closed loop control system (T) to the variation in open loop gain (G) is defined as
∂T
T Percentage change in T
S TG = ∂G
= Equation3
Percentage change in G
G
∂T
∂G
=
∂
∂G ( )
G
1 + GH
=
( 1 + GH ) .1 − G ( H )
( 1 + GH ) 2
=
1
( 1 + GH ) 2
Equation5
1 1
S TG = (1 + GH) =
(1 + GH) 2 1 + GH
So, we got the sensitivity of the overall gain of closed loop control system as the reciprocal of 1 + GH. So, Sensitivity may increase or decrease
depending on the value of 1 + GH.
If the value of 1 + GH is less than 1, then sensitivity increases. In this case, 'GH' value is negative because the gain of feedback path is negative.
If the value of 1 + GH is greater than 1, then sensitivity decreases. In this case, 'GH' value is positive because the gain of feedback path is positive.
In general, 'G' and 'H' are functions of frequency. So, feedback will increase the sensitivity of the system gain in one frequency range and decrease in the
other frequency range. Therefore, we have to choose the values of 'GH' in such a way that the system is insensitive or less sensitive to parameter
variations.
A system is said to be stable, if its output is under control. Otherwise, it is said to be unstable.
In Equation 2, if the denominator value is zero i. e. , GH = − 1, then the output of the control system will be infinite. So, the control system
becomes unstable.
Therefore, we have to properly choose the feedback in order to make the control system stable.
Consider an open loop control system with noise signal as shown below.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 4/90
17/01/2018 Control Systems Quick Guide
Consider a closed loop control system with noise signal as shown below.
C(s) Gb
N(s)
= 1 + G aG bH
Equation8
In the closed loop control system, the gain due to noise signal is decreased by a factor of (1 + G aG bH) provided that the term (1 + G aG bH) is greater than
one.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 5/90
17/01/2018 Control Systems Quick Guide
Get the differential equation in terms of input and output by eliminating the intermediate variables.
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 v i and the voltage across the capacitor is the output voltage v o.
di
v i = Ri + L + vo
dt
dv o
Substitute, the current passing through capacitor i = c dt
in the above equation.
dv o d 2v o
⇒ v i = RC + LC + vo
dt dt 2
d 2v o
⇒
dt 2
+ () ( ) ( )
R dv o
L dt
+
1
LC
vo =
LC
1
vi
If x(t) and y(t) are the input and output of an LTI system, then the corresponding Laplace transforms are X(s) and Y(s).
Therefore, the transfer function of LTI system is equal to the ratio of Y(s) and X(s).
Y(s)
i. e. , Transfer Function =
X(s)
The transfer function model of an LTI system is shown in the following figure.
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).
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 6/90
17/01/2018 Control Systems Quick Guide
Example
d 2v o
dt 2
+ () ( ) ( )
R dv o
L dt
+
LC
1
vo =
1
LC
vi
s 2V o(s) + ( )
sR
L
V o(s) + ( )
1
LC
V o(s) = ( ) 1
LC
V i(s)
⇒ { ()
s2 +
R
L
s+
1
LC } V o(s) = ( )
LC
1
V i(s)
1
V o(s) LC
⇒ =
V i(s)
s2 + ()R
L
s+
1
LC
Where,
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 V i(s) & an output V o(s).
If a force is applied to a translational mechanical system, then it is opposed by opposing forces due to mass, elasticity and friction of the system. Since
the applied force and the opposing forces are in opposite directions, the algebraic sum of the forces acting on the system is zero. Let us now see the force
opposed by these three elements individually.
Mass
Mass is the property of a body, which stores kinetic energy. If a force is applied on a body having mass M, then it is opposed by an opposing force due
to mass. This opposing force is proportional to the acceleration of the body. Assume elasticity and friction are negligible.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 7/90
17/01/2018 Control Systems Quick Guide
Fm ∝ a
d 2x
⇒ F m = Ma = M
dt 2
d 2x
F = Fm = M
dt 2
Where,
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.
F∝ x
⇒ F k = Kx
F = F k = Kx
Where,
K is spring constant
x is displacement
Dashpot
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 8/90
17/01/2018 Control Systems Quick Guide
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 ∝ ν
dx
⇒ F b = Bν = B
dt
dx
F = Fb = B
dt
Where,
v is velocity
x is displacement
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 ∝ α
d 2θ
⇒ T j = Jα = J
dt 2
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 9/90
17/01/2018 Control Systems Quick Guide
d 2θ
T = Tj = J
dt 2
Where,
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.
Tk ∝ θ
⇒ T k = Kθ
T = T k = Kθ
Where,
θ 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.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 10/90
17/01/2018 Control Systems Quick Guide
Tb ∝ ω
dθ
⇒ T b = Bω = B
dt
dθ
T = Tb = B
dt
Where,
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.
Consider the following translational mechanical system as shown in the following figure.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 11/90
17/01/2018 Control Systems Quick Guide
F = Fm + Fb + Fk
d 2x dx
⇒F=M + B dt + Kx Equation1
dt 2
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.
dq
Substitute, i = in Equation 2.
dt
dq d 2q q
V=R +L 2 +
dt dt C
⇒V=L
d 2q
dt 2
dq
+ R dt + ()
1
c
q Equation3
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.
ForceF VoltageV
MassM InductanceL
Spring ConstantK 1
Reciprocal of Capacitance ( c )
Displacementx Chargeq
Velocityv Currenti
Similarly, there is torque voltage analogy for rotational mechanical systems. Let us now discuss about this analogy.
In this analogy, the mathematical equations of rotational mechanical system are compared with mesh equations of the electrical system.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 12/90
17/01/2018 Control Systems Quick Guide
T = Tj + Tb + Tk
d 2θ dθ
⇒T=J + B dt + kθ Equation4
dt 2
By comparing Equation 4 and Equation 3, we will get the analogous quantities of rotational mechanical system and electrical system. The following table
shows these analogous quantities.
TorqueT VoltageV
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.
dΨ
Substitute, V = in Equation 5.
dt
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 13/90
17/01/2018 Control Systems Quick Guide
i=
1 dΨ
R dt
+ ()
1
L
Ψ+C
d 2Ψ
dt 2
⇒i=C
d 2Ψ
dt 2
+ () ()
1
R
dΨ
dt
+
1
L
Ψ Equation6
By comparing Equation 1 and Equation 6, we will get the analogous quantities of the translational mechanical system and electrical system. The
following table shows these analogous quantities.
ForceF Currenti
MassM CapacitanceC
Frictional coefficientB 1
Reciprocal of Resistance( R )
Spring constantK 1
Reciprocal of Inductance( L )
Velocityv VoltageV
Similarly, there is a torque current analogy for rotational mechanical systems. Let us now discuss this analogy.
In this analogy, the mathematical equations of the rotational mechanical system are compared with the nodal mesh equations of the electrical
system.
By comparing Equation 4 and Equation 6, we will get the analogous quantities of rotational mechanical system and electrical system. The following table
shows these analogous quantities.
TorqueT Currenti
In this chapter, we discussed the electrical analogies of the mechanical systems. These analogies are helpful to study and analyze the non-electrical
system like mechanical system from analogous electrical system.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 14/90
17/01/2018 Control Systems Quick Guide
The above block diagram consists of two blocks having transfer functions Gs and Hs. 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 Xs, output Ys and the transfer function Gs.
Y(s)
Transfer Function, G(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.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 15/90
17/01/2018 Control Systems Quick Guide
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
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, Rs to two more blocks.
In the following figure, the take-off point is used to connect the output Cs, as one of the inputs to the summing point.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 16/90
17/01/2018 Control Systems Quick Guide
Consider a series of RLC circuit as shown in the following figure. Where, Vit and Vot are the input and output voltages. Let it be the current passing
through the circuit. This circuit is in time domain.
By applying the Laplace transform to this circuit, will get the circuit in s-domain. The circuit is as shown in the following figure.
V i(s) − V o(s)
I(s) =
R + sL
⇒ I(s) = { }{ 1
R + sL
V i(s) − V o(s) } Equation1
V o(s) = ( )
sC
1
I(s) Equation2
Let us now draw the block diagrams for these two equations individually. And then combine those block diagrams properly in order to get the overall
block diagram of series of RLC Circuit s − domain.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 17/90
17/01/2018 Control Systems Quick Guide
1
Equation 1 can be implemented with a block having the transfer function, R + sL { }
. The input and output of this block are V i(s) − V o(s) and I(s). We
{ }
require a summing point to get V i(s) − V o(s) . The block diagram of Equation 1 is shown in the following figure.
1
Equation 2 can be implemented with a block having transfer function, sC
. The input and output of this block are I(s) and V o(s). The block diagram of
Equation 2 is shown in the following figure.
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.
Series Connection
Series connection is also called cascade connection. In the following figure, two blocks having transfer functions G 1(s) and G 2(s) are connected in
series.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 18/90
17/01/2018 Control Systems Quick Guide
Y(s) = G 2(s)Z(s)
Compare this equation with the standard form of the output equation, Y(s) = G(s)X(s). Where, G(s) = G 1(s)G 2(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 G 1(s) and G 2(s) are
connected in parallel. The outputs of these two blocks are connected to the summing point.
Compare this equation with the standard form of the output equation, Y(s) = G(s)X(s).
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 19/90
17/01/2018 Control Systems Quick Guide
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.
Y(s) = E(s)G(s)
Y(s) G(s)
⇒ =
X(s) 1 + G(s)H(s)
G(s)
Therefore, the negative feedback closed loop transfer function is 1+G(s)H(s)
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
G(s)
loop transfer function of the positive feedback, i.e., 1−G(s)H(s)
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 20/90
17/01/2018 Control Systems Quick Guide
There are two possibilities of shifting summing points with respect to blocks −
Let us now see what kind of arrangements need to be done in the above two cases one by one.
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 –
Now, shift the summing point after the block. This block diagram is shown in the following figure.
The first term ‘G(s)R(s) ′ is same in both the equations. But, there is difference in the second term. In order to get the second term also same, we require
one more block G(s). It is having the input X(s) and the output of this block is given as input to summing point instead of X(s). This block diagram is
shown in the following figure.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 21/90
17/01/2018 Control Systems Quick Guide
Consider the block diagram shown in the following figure. Here, the summing point is present after the block.
Now, shift the summing point before the block. This block diagram is shown in the following figure.
The first term ‘G(s)R(s) ′ is same in both equations. But, there is difference in the second term. In order to get the second term also same, we require one
1
more block . It is having the input X(s) and the output of this block is given as input to summing point instead of X(s). This block diagram is shown
G(s)
in the following figure.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 22/90
17/01/2018 Control Systems Quick Guide
Let us now see what kind of arrangements are to be done in the above two cases, one by one.
Consider the block diagram shown in the following figure. In this case, the take-off point is present before the block.
When you shift the take-off point after the block, the output Y(s) will be same. But, there is difference in X(s) value. So, in order to get the same X(s) value,
1
we require one more block . It is having the input Y(s) and the output is X(s). This block diagram is shown in the following figure.
G(s)
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 23/90
17/01/2018 Control Systems Quick Guide
Consider the block diagram shown in the following figure. Here, the take-off point is present after the block.
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.
Rule 3 − Check for the blocks connected in feedback loop and simplify.
Rule 4 − If there is difficulty with take-off point while simplifying, shift it towards right.
Rule 5 − If there is difficulty with summing point while simplifying, shift it towards left.
Rule 6 − Repeat the above steps till you get the simplified form, i.e., single block.
Note − The transfer function present in this single block is the transfer function of the overall block diagram.
Example
Consider the block diagram shown in the following figure. Let us simplify reduce this block diagram using the block diagram reduction rules.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 24/90
17/01/2018 Control Systems Quick Guide
Step 1 − Use Rule 1 for blocks G 1 and G 2. Use Rule 2 for blocks G 3 and G 4. The modified block diagram is shown in the following figure.
Step 2 − Use Rule 3 for blocks G 1G 2 and H 1. Use Rule 4 for shifting take-off point after the block G 5. The modified block diagram is shown in the
following figure.
Step 3 − Use Rule 1 for blocks (G 3 + G 4) and G 5. The modified block diagram is shown in the following figure.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 25/90
17/01/2018 Control Systems Quick Guide
Step 4 − Use Rule 3 for blocks (G 3 + G 4)G 5 and H 3. 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.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 26/90
17/01/2018 Control Systems Quick Guide
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 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 partiallysimplified 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.
Node
Node is a point which represents either a variable or a signal. There are three types of nodes — input node, output node and mixed node.
Mixed Node − It is a node, which has both incoming and outgoing branches.
Example
Let us consider the following signal flow graph to identify these nodes.
The nodes present in this signal flow graph are y1, y2, y3 and y4.
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.
y 2 = a 12y 1 + a 42y 4
y 3 = a 23y 2 + a 53y 5
y 4 = a 34y 3
y 5 = a 45y 4 + a 35y 3
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 27/90
17/01/2018 Control Systems Quick Guide
y 6 = a 56y 5
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 y 2 = a 13y 1 + a 42y 4 is shown in the following figure.
Step 2 − Signal flow graph for y 3 = a 23y 2 + a 53y 5 is shown in the following figure.
Step 3 − Signal flow graph for y 4 = a 34y 3 is shown in the following figure.
Step 4 − Signal flow graph for y 5 = a 45y 4 + a 35y 3 is shown in the following figure.
Step 5 − Signal flow graph for y 6 = a 56y 5 is shown in the following figure.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 28/90
17/01/2018 Control Systems Quick Guide
Step 6 − Signal flow graph of overall system is shown in the following figure.
Represent all the signals, variables, summing points and take-off points of block diagram as nodes in signal flow graph.
Represent the transfer functions inside the blocks of block diagram as gains of the branches in signal flow graph.
Connect the nodes as per the block diagram. If there is connection between two nodes butthereisnoblockinbetween, then represent the gain of the
branch as one. For example, between summing points, between summing point and takeoff point, between input and summing point, between
take-off point and output.
Example
Let us convert the following block diagram into its equivalent signal flow graph.
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 G 1 and G 2.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 29/90
17/01/2018 Control Systems Quick Guide
With the help of Mason’s gain formula discussedinthenextchapter, 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.
C(s) ΣN PΔ
i=1 i i
T= =
R(s) Δ
Where,
Δ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.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 30/90
17/01/2018 Control Systems Quick Guide
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 − y 2 → y 3 → y 4 → y 5 and y 5 → y 3 → y 2
Forward Path
The path that exists from the input node to the output node is known as forward path.
Examples − y 1 → y 2 → y 3 → y 4 → y 5 → y 6 and y 1 → y 2 → y 3 → y 5 → y 6.
It is obtained by calculating the product of all branch gains of the forward path.
Examples − abcde is the forward path gain of y 1 → y 2 → y 3 → y 4 → y 5 → y 6 and abge is the forward path gain of y 1 → y 2 → y 3 → y 5 → y 6.
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 − y 2 → y 3 → y 2 and y 3 → y 5 → y 3.
Loop Gain
Non-touching Loops
These are the loops, which should not have any common node.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 31/90
17/01/2018 Control Systems Quick Guide
Higher number of morethantwo non-touching loops are not present in this signal flow graph.
We know,
So, Δ 1 = 1.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 32/90
17/01/2018 Control Systems Quick Guide
2
C(s) Σ i = 1P iΔ i
T= =
R(s) Δ
C(s) P 1Δ 1 + P 2Δ 2
T= =
R(s) Δ
Transient response
Steady state response
The response of control system in time domain is shown in the following figure.
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.
Where,
Transient Response
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 33/90
17/01/2018 Control Systems Quick Guide
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.
lim c tr(t) = 0
t→∞
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
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.
δ(t) = 0 for t ≠ 0
0+
and ∫ 0 − δ(t)dt = 1
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’.
u(t) = 1; t ≥ 0
= 0; t < 0
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 34/90
17/01/2018 Control Systems Quick Guide
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’.
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)
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’.
t2
p(t) = ;t ≥ 0
2
= 0; t < 0
We can write unit parabolic signal, p(t) in terms of the unit step signal, u(t) as,
t2
p(t) = u(t)
2
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 35/90
17/01/2018 Control Systems Quick Guide
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’.
We know that the transfer function of the closed loop control system has unity negative feedback as,
C(s) G(s)
=
R(s) 1 + G(s)
1
Substitute, G(s) = in the above equation.
sT
1
C(s) sT 1
= =
R(s) 1 sT + 1
1+ sT
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.
C(s) = ( )
1
sT + 1
R(s)
Where,
Follow these steps to get the response output of the first order system in the time domain.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 36/90
17/01/2018 Control Systems Quick Guide
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.
R(s) = 1
C(s) = ( ) 1
sT + 1
(1) =
1
sT + 1
Rearrange the above equation in one of the standard forms of Laplace transforms.
C(s) =
( )
T s+
1
1
T
⇒ C(s) =
1
T
( )
s+
1
1
T
c(t) =
1
T
e ( ) u(t)
t
−T
The unit impulse response, ct is an exponential decaying signal for positive values of ‘t’ and it is zero for negative values of ‘t’.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 37/90
17/01/2018 Control Systems Quick Guide
1
R(s) =
s
1
Substitute, R(s) = in the above equation.
s
C(s) = ( )( )1
sT + 1
1
s
=
1
s(sT + 1)
1 A B
C(s) = = +
s(sT + 1) s sT + 1
1 A(sT + 1) + Bs
⇒ =
s(sT + 1) s(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
1 T 1 T
C(s) = − = −
s sT + 1 s
T s+( ) 1
T
1 1
⇒ C(s) = −
s 1
s+ T
(
c(t) = 1 − e − ( ) )u(t)
t
T
The unit step response, ct has both the transient and the steady state terms.
c tr(t) = − e − ( ) u(t)
t
T
c ss(t) = u(t)
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 38/90
17/01/2018 Control Systems Quick Guide
The value of the unit step response, ct 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.
1
R(s) =
s2
1
Substitute, R(s) = in the above equation.
s2
C(s) = ( )( 1
sT + 1 s
1
2 ) = 2
1
s (sT + 1)
1 A B C
C(s) = = + +
2
s (sT + 1) s 2 s sT + 1
1 A(sT + 1) + Bs(sT + 1) + Cs 2
⇒ 2
=
s (sT + 1) s 2(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) + Cs 2
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
Similarly, substitute B = −T and equate the coefficient of s 2 terms on both the sides. You will get C = T 2.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 39/90
17/01/2018 Control Systems Quick Guide
1 T T2 1 T T2
C(s) = − + = 2 − +
2 s sT + 1 s
( )
s s 1
T s+ T
1 T T
⇒ C(s) = − +
s2 s 1
s+ T
(
c(t) = t − T + Te − ( ) )u(t)
t
T
The unit ramp response, ct has both the transient and the steady state terms.
c tr(t) = Te − ( ) u(t)
t
T
c ss(t) = (t − T)u(t)
The unit ramp response, ct follows the unit ramp input signal for all positive values of t. But, there is a deviation of T units from the input signal.
t2
So, r(t) = u(t)
2
1
R(s) =
s3
1
Substitute R(s) = in the above equation.
s3
C(s) = ( )( 1
sT + 1
1
s3 ) =
1
s 3(sT + 1)
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 40/90
17/01/2018 Control Systems Quick Guide
1 A B C D
C(s) = = + + +
s 3(sT + 1) s3 s2 s sT + 1
After simplifying, you will get the values of A, B, C and D as 1, − T, T 2 and − T 3 respectively. Substitute these values in the above partial fraction
expansion of Cs.
1 T T2 T3 1 T T2 T2
C(s) = − + − ⇒ C(s) = − + −
s3 s2 s sT + 1 s3 s2 s
s+ T
1
c(t) =
( t2
2
− Tt + T 2 − T 2e − ( ) )u(t)
t
T
The unit parabolic response, ct has both the transient and the steady state terms.
C tr(t) = − T 2e − ( ) u(t)
T
t
C ss(t) =
( t2
2
− Tt + T 2 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) G(s)
=
R(s) 1 + G(s)
ω 2n
Substitute, G(s) = in the above equation.
s ( s + 2δω n )
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 41/90
17/01/2018 Control Systems Quick Guide
( )
2
ωn
s ( s + 2δω n )
C(s) ω 2n
= =
R(s) s 2 + 2δω ns + ω 2n
1+
( ) ω 2n
s ( s + 2δω n )
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.
s 2 + 2δω ns + ω 2n = 0
− 2ωδ n ±
√(2δω ) n
2 − 4ω 2n
√
− 2(δω n ± ω n δ 2 − 1)
s= =
2 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.
C(s) =
( ω 2n
s 2 + 2δω ns + ω 2n ) R(s)
Where,
Follow these steps to get the response output of the second order system in the time domain.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 42/90
17/01/2018 Control Systems Quick Guide
1
R(s) =
s
We know the transfer function of the second order closed loop control system is,
C(s) ω 2n
=
R(s) s 2 + 2δω ns + ω 2n
Case 1: δ = 0
Substitute, δ = 0 in the transfer function.
2
C(s) ωn
=
R(s) s 2 + ω 2n
⇒ C(s) =
( ) ω 2n
s2 + ωn
2
R(s)
1
Substitute, R(s) = s
in the above equation.
( )( )
2 2
ωn 1 ωn
C(s) = =
s 2 + ω 2n s s(s 2 + ω 2n)
(
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
C(s) ω 2n
=
R(s) s 2 + 2ω ns + ω 2n
( )
2
ωn
⇒ C(s) = R(s)
(s + ω n) 2
1
Substitute, R(s) = in the above equation.
s
C(s) =
( )( ) ω 2n
(s + ω n) 2
1
s
=
ω 2n
s(s + ω n) 2
ω 2n A B C
C(s) = = + +
s(s + ω n) 2 s s + ω n (s + ω ) 2
n
After simplifying, you will get the values of A, B and C as 1, − 1 and − ω n respectively. Substitute these values in the above partial fraction expansion of
C(s).
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 43/90
17/01/2018 Control Systems Quick Guide
1 1 ωn
C(s) = − −
s s + ωn (s + ω n) 2
So, the unit step response of the second order system will try to reach the step input in steady state.
s 2 + 2δω ns + ω n =
2
{s 2
}
+ 2(s)(δω n) + (δω n) 2 + ω n − (δω n) 2
2
= (s + δω n) 2 + ω 2n(1 − δ 2)
C(s) ω 2n
=
R(s) (s + δω n) 2 + ω 2n(1 − δ 2)
⇒ C(s) =
( ω 2n
(s + δω n) 2 + ω 2n(1 − δ 2) ) R(s)
1
Substitute, R(s) = s
in the above equation.
C(s) =
( ω 2n
(s + δω n) 2 + ω 2n(1 − δ 2) )( ) 1
s
=
(
ω 2n
s (s + δω n) 2 + ω 2n(1 − δ 2) )
Do partial fractions of C(s).
ω 2n A Bs + C
C(s) = = +
s (s + δω ) 2 + ω 2(1 − δ 2)
(
s (s + δω n) 2 + ω 2n(1 − δ 2) ) n n
After simplifying, you will get the values of A, B and C as 1, − 1 and − 2δω n respectively. Substitute these values in the above partial fraction expansion of
Cs.
1 s + 2δω n
C(s) = −
s (s + δω n) 2 + ω 2n(1 − δ 2)
1 s + δω n δω n
C(s) = − −
s (s + δω n) 2 + ω 2n(1 − δ 2) (s + δω n) 2 + ω 2n(1 − δ 2)
C(s) =
1
s
−
( 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.
√
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 44/90
17/01/2018 Control Systems Quick Guide
C(s) =
1
s
−
(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 − δ 2
e − δω 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 decreasingamplitude when ‘δ’ lies between zero and one.
Case 4: δ > 1
s 2 + 2δω ns + ω n =
2
{s 2
}
+ 2(s)(δω n) + (δω n) 2 + ω n − (δω n) 2
2
(
= s + δω n )2 − ω 2n (δ2 − 1 )
The transfer function becomes,
C(s) ω 2n
=
R(s) (s + δω n) 2 − ω 2n(δ 2 − 1)
⇒ C(s) =
( ω 2n
(s + δω n) 2 − ω 2n(δ 2 − 1) ) R(s)
1
Substitute, R(s) = s
in the above equation.
( )( )
2 2
ωn 1 ωn
C(s) = =
s
√
( s + δω n ) 2 − ( ω n δ 2 − 1 ) 2 √ √
s ( s + δω n + ω n δ 2 − 1 ) ( s + δω n − ω n δ 2 − 1 )
2
ωn
C(s) =
s(s + δω n + ω n δ 2 − 1)(s + δω n − ω n δ 2 − 1)
√ √
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 45/90
17/01/2018 Control Systems Quick Guide
A B C
= + +
s s + δω + ω δ 2 − 1 s + δω − ω δ 2 − 1
n n n n √ √
1 −1
After simplifying, you will get the values of A, B and C as 1, and respectively. Substitute these values in above
2(δ+ √ δ2 − 1 ) ( √ δ2 − 1 ) 2(δ− √ δ2 − 1 ) ( √δ 2 − 1 )
partial fraction expansion of C(s).
C(s) =
1
+
s 2(δ +
√ δ2
1
− 1)( √ δ2 (
− 1) s + δω n + ω n
1
√ δ2 −1 )(−
2(δ − √ δ2
1
− 1)( √ δ2 − 1) )( 1
√
s + δω n − ω n δ 2 − 1 )
Apply inverse Laplace transform on both the sides.
c(t) =
( (
1+
2(δ+ √
1
δ2 − 1 ) ( √ δ2 − 1 ) )
e − ( δω n + ω n√δ
2−1)t
−
( 2(δ− √
1
δ2 − 1 ) ( √ δ2 − 1 ) ) e − ( δω n − ω n√δ
2−1)t
) u(t)
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.
The following table shows the impulse response of the second order system for 4 cases of the damping ratio.
δ=1 ω 2nte − ω nt
0<δ<1
( )
ω ne − δω nt
√1 − δ 2
sin(ω dt)
δ>1
( )(
ωn
2 δ2 − 1
√
e − ( δω n − ω n√δ
2
−1)t − e − ( δω n + ω n√δ
2
−1)t
)
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 46/90
17/01/2018 Control Systems Quick Guide
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 t d.
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 + θ)
Therefore, at t = t d, the value of the step response will be 0.5. Substitute, these values in the above equation.
(√ )
c(t d) = 0.5 = 1 −
e − δω nt d
1 − δ2
sin(ω dt d + θ)
⇒
(√ )
e − δω nt d
1 − δ2
sin(ω dt d + θ) = 0.5
1 + 0.7δ
td =
ω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.
At t = t1 = 0, ct = 0.
Therefore, at t = t 2, the value of step response is one. Substitute, these values in the following equation.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 47/90
17/01/2018 Control Systems Quick Guide
(√ )
c(t) = 1 −
e − δω nt
1 − δ2
sin(ω dt + θ)
(√ )
c(t 2) = 1 = 1 −
e − δω nt 2
1 − δ2
sin(ω dt 2 + θ)
⇒
(√ ) e − δω nt 2
1 − δ2
sin(ω dt 2 + θ) = 0
⇒ sin(ω dt 2 + θ) = 0
⇒ ω dt 2 + θ = π
π−θ
⇒ t2 =
ωd
tr = t2 − t1
π−θ
∴ tr =
ωd
From above equation, we can conclude that the rise time t r 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 t p. At t = t p, 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 + θ)
dc(t)
dt
= −
(√ )
e − δω nt
1 − δ2
ω dcos(ω dt + θ) −
( − δω ne − δω nt
√1 − δ 2 ) sin(ω dt + θ)
dc ( t )
Substitute, t = t p and = 0 in the above equation.
dt
0= −
(√ )[
e − δω nt p
1 − δ2
ω dcos(ω dt p + θ) − δω nsin(ω dt p + θ) ]
⇒ ω n 1 − δ 2cos(ω dt p + θ) − δω nsin(ω dt p + θ) = 0
√
⇒ √1 − δ2cos(ωdtp + θ) − δsin(ωdtp + θ) = 0
⇒ sin(θ)cos(ω dt p + θ) − cos(θ)sin(ω dt p + θ) = 0
⇒ sin(θ − ω dt p − θ) = 0
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 48/90
17/01/2018 Control Systems Quick Guide
⇒ ω dt p = π
π
⇒ tp =
ωd
From the above equation, we can conclude that the peak time t p and the damped frequency ω d are inversely proportional to each other.
Peak Overshoot
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.
M p = c(t p) − c(∞)
Where,
At t = t p, the response ct is -
c(t p) = 1 −
(√ )
e − δω nt p
1 − δ2
sin(ω dt p + θ)
π
Substitute, t p = in the right hand side of the above equation.
ωd
( )
( )
π
e − δω n
(( ) )
ωd
π
c(t P) = 1 − sin ω d +θ
ωd
√1 − δ2
( )
( )
δπ
e− √1 − δ 2
⇒ c(t p) = 1 − ( − sin(θ))
√1 − δ 2
We know that
sin(θ) = √1 − δ 2
So, we will get c(t p) as
c(t p) = 1 + e − ( ) δπ
√1 − δ 2
Substitute the values of c(t p) and c(∞) in the peak overshoot equation.
Mp = 1 + e − ( )δπ
√1 − δ 2 −1
⇒ Mp = e − ( ) δπ
√1 − δ 2
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 49/90
17/01/2018 Control Systems Quick Guide
Mp
%M p = × 100%
c(∞)
By substituting the values of M p and c(∞) in above formula, we will get the Percentage of the peak overshoot %M p as
( ( ))
%M p = e −
δπ
√1 − δ 2 × 100%
From the above equation, we can conclude that the percentage of peak overshoot %M p 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 t s.
3
ts = = 3τ
δω n
4
ts = = 4τ
δω n
1
Where, τ is the time constant and is equal to .
δω n
Both the settling time t s and the time constant τ are inversely proportional to the damping ratio δ.
Both the settling time t s and the time constant τ are independent of the system gain. That means even the system gain changes, the settling time t s
and time constant τ will never change.
Example
4
Let us now find the time domain specifications of a control system having the closed loop transfer function when the unit step signal is applied
s 2 + 2s + 4
as an input to this control system.
We know that the standard form of the transfer function of the second order closed loop control system as
ω 2n
s 2 + 2δω ns + ω 2n
By equating these two transfer functions, we will get the un-damped natural frequency ω n as 2 rad/sec and the damping ratio δ as 0.5.
√
ωd = ωn 1 − δ2
⇒ ω d = 2 1 − (0.5) 2
√
⇒ ω d = 1.732 rad / sec
θ = cos − 1δ
π
⇒ θ = cos − 1(0.5) = rad
3
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 50/90
17/01/2018 Control Systems Quick Guide
Substitute the above necessary values in the formula of each time domain specification and simplify in order to get the values of time domain
specifications for given transfer function.
The following table shows the formulae of time domain specifications, substitution of necessary values and the final values.
( ( )) ( ( ))
% Peak overshoot % M p=16.32%
δπ 0.5π
%M p = e − √ 1 − δ2 × 100% %M p = e − √1 − ( 0.5 ) 2 × 100%
Where,
Let us discuss how to find steady state errors for unity feedback and non-unity feedback control systems one by one.
Where,
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 51/90
17/01/2018 Control Systems Quick Guide
We know the transfer function of the unity negative feedback closed loop control system as
C(s) G(s)
=
R(s) 1 + G(s)
R(s)G(s)
⇒ C(s) =
1 + G(s)
R(s)G(s)
E(s) = R(s) −
1 + G(s)
R(s)
⇒ E(s) =
1 + G(s)
sR(s)
e ss = lim
s→0 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.
1
unit parabolic signal K a = lim s → 0s 2G(s)
Ka
Where, K p, K v and K a are position error constant, velocity error constant and acceleration error constant respectively.
Note − If any of the above input signals has the amplitude other than unity, then multiply corresponding steady state error with that amplitude.
Note − We can’t define the steady state error for the unit impulse signal because, it exists only at origin. So, we can’t compare the impulse response with
the unit impulse input as t denotes infinity.
Example
(
Let us find the steady state error for an input signal r(t) = 5 + 2t +
t2
2 )u(t) of unity negative feedback control system with G(s) =
5(s+4)
s 2 ( s + 1 ) ( s + 20 )
The given input signal is a combination of three signals step, ramp and parabolic. The following table shows the error constants and steady state error
values for these three signals.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 52/90
17/01/2018 Control Systems Quick Guide
t2 K a = lim s → 0s 2G(s) = 1 1
r 3(t) = u(t) e ss3 = =1
2 ka
We will get the overall steady state error, by adding the above three steady state errors.
⇒ e ss = 0 + 0 + 1 = 1
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.
Simplify the above block diagram by keeping the unity negative feedback as it is. The following is the simplified block diagram.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 53/90
17/01/2018 Control Systems Quick Guide
This block diagram resembles the block diagram of the unity negative feedback closed loop control system. Here, the single block is having the transfer
G(s)
function 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.
What is Stability?
A system is said to be stable, if its output is under control. Otherwise, it is said to be unstable. A stable system produces a bounded output for a given
bounded input.
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 t including zero. So, it is bounded input. Therefore, the first order control system is
stable since both the input and the output are bounded.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 54/90
17/01/2018 Control Systems Quick Guide
stable if all the poles of the closed loop transfer function present in the left half of the ‘s’ plane.
If the system is stable for a certain range of system component values, then it is known as conditionally stable system.
If the system is stable by producing an output signal with constant amplitude and constant frequency of oscillations for bounded input, then it is known
as marginally stable system. The 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.
The necessary condition is that the coefficients of the characteristic polynomial should be positive. This implies that all the roots of the characteristic
equation should have negative real parts.
a 0s n + a 1s n − 1 + a 2s n − 2 + . . . + a n − 1 s 1 + a ns 0 = 0
Note that, there should not be any term missing in the nth order characteristic equation. This means that the nth order characteristic equation should
not have any coefficient that is of zero value.
The sufficient condition is that all the elements of the first column of the Routh array should have the same sign. This means that all the elements of the
first column of the Routh array should be either positive or negative.
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.
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 s n and continue up to the coefficient of s 0.
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 s 0 is a n. Here, a n is the coefficient of s 0 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.
a 0s n + a 1s n − 1 + a 2s n − 2 + . . . + a n − 1 s 1 + a ns 0
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 55/90
17/01/2018 Control Systems Quick Guide
sn a0 a2 a4 a6 ... ...
sn − 1 a1 a3 a5 a7 ... ...
sn − 3 b 1a 3 − b 2a 1 b 1a 55 − b 3a 1 ⋮
c1 = b1
c2 = b1
⋮ ⋮ ⋮ ⋮
s1 ⋮ ⋮
s0 an
Example
Let us find the stability of the control system having characteristic equation,
s 4 + 3s 3 + 3s 2 + 2s + 1 = 0
All the coefficients of the characteristic polynomial, s 4 + 3s 3 + 3s 2 + 2s + 1 are positive. So, the control system satisfies the necessary condition.
Step 2 − Form the Routh array for the given characteristic polynomial.
s4 1 3 1
s3 3 2
s2 (3×3) − (2×1)
=
7 (3×1) − (0×1)
=
3
=1
3 3 3 3
s1
( )
7
3
×2 − (1×3)
5
7
= 7
3
s0 1
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 56/90
17/01/2018 Control Systems Quick Guide
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.
Let us now discuss how to overcome the difficulty in these two cases, one by one.
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,
s 4 + 2s 3 + s 2 + 2s + 1 = 0
All the coefficients of the characteristic polynomial, s 4 + 2s 3 + s 2 + 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
s2 (1×1) − (1×1)
=0
(1×1) − (0×1)
=1
1 1
s1
s0
The row s 3 elements have 2 as the common factor. So, all these elements are divided by 2.
Special case i − Only the first element of row s 2 is zero. So, replace it by ϵ and continue the process of completing the Routh table.
s4 1 1 1
s3 1 1
s2 ϵ 1
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 57/90
17/01/2018 Control Systems Quick Guide
s1 (ϵ×1) − (1×1)
=
ϵ−1
ϵ ϵ
s0 1
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
Write the auxilary equation, As of the row, which is just above the row of zeros.
Differentiate the auxiliary equation, As with respect to s. Fill the row of zeros with these coefficients.
Example
Let us find the stability of the control system having characteristic equation,
s 5 + 3s 4 + s 3 + 3s 2 + s + 3 = 0
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)
=0
(1×1) − (1×1)
=0
1 1
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 58/90
17/01/2018 Control Systems Quick Guide
s2
s1
s0
The row s 4 elements have the common factor of 3. So, all these elements are divided by 3.
Special case ii − All the elements of row s 3 are zero. So, write the auxiliary equation, As of the row s 4.
A(s) = s 4 + s 2 + 1
dA(s)
= 4s 3 + 2s
ds
s5 1 1 1
s4 1 1 1
s3 42 21
s2 (2×1) − (1×1)
= 0.5
(2×1) − (0×1)
=1
2 2
s1 ( 0.5 × 1 ) − ( 1 × 2 )
=
− 1.5
= −3
0.5 0.5
s0 1
There are two sign changes in the first column of Routh table. Hence, the control system is unstable.
In the Routh-Hurwitz stability criterion, we can know whether the closed loop poles are in on left half of the ‘s’ plane or on the right half of the ‘s’ plane
or on an imaginary axis. So, we can’t find the nature of the control system. To overcome this limitation, there is a technique known as the root locus. We
will discuss this technique in the next two chapters.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 59/90
17/01/2018 Control Systems Quick Guide
The Root locus is the locus of the roots of the characteristic equation by varying system gain K from zero to infinity.
We know that, the characteristic equation of the closed loop control system is
1 + G(s)H(s) = 0
N(s)
G(s)H(s) = K
D(s)
Where,
Ns represents the numerator term having factored nth order polynomial of ‘s’.
Ds represents the denominator term having factored mth order polynomial of ‘s’.
N(s)
1+k =0
D(s)
⇒ D(s) + KN(s) = 0
Case 1 − K = 0
If K = 0, then D(s) = 0.
That means, the closed loop poles are equal to open loop poles when K is zero.
Case 2 − K = ∞
K ( 1
K
+
N(s)
D(s) ) =0⇒
1
K
+
N(s)
D(s)
=0
1 N(s) N(s)
+ =0⇒ = 0 ⇒ N(s) = 0
∞ D(s) D(s)
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 points on the root locus branches satisfy the angle condition. So, the angle condition is used to know whether the point exist on root locus branch or
not. We can find the value of K for the points on the root locus branches by using magnitude condition. So, we can use the magnitude condition for the
points, and this satisfies the angle condition.
1 + G(s)H(s) = 0
⇒ G(s)H(s) = − 1 + j0
∠G(s)H(s) = tan − 1 ( )
−1
0
= (2n + 1)π
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
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 60/90
17/01/2018 Control Systems Quick Guide
The magnitude condition is that the point whichsatisfiedtheanglecondition at which the magnitude of the open loop transfer function is one.
Rule 1 − Locate the open loop poles and zeros in the ‘s’ plane.
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.
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 1800, 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.
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.
(2q + 1)180 0
θ=
P−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.
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 K for 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.
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.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 61/90
17/01/2018 Control Systems Quick Guide
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.
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.
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.
ϕ d = 180 0 − ϕ
ϕ a = 180 0 + ϕ
Where,
ϕ= ∑ ϕP − ∑ ϕZ
Example
K
Let us now draw the root locus of the control system having open loop transfer function, G(s)H(s) = s(s+1) (s+5)
Step 1 − The given open loop transfer function has three poles at s = 0, s = − 1 and 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
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 62/90
17/01/2018 Control Systems Quick Guide
The centroid and three asymptotes are shown in the following figure.
Step 3 − Since two asymptotes have the angles of 60 0 and 300 0, two root locus branches intersect the imaginary axis. By using the Routh array method
and special caseii, 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.
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.
The root locus can be shifted in ‘s’ plane by adding the open loop poles and the open loop zeros.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 63/90
17/01/2018 Control Systems Quick Guide
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 t d, rise time t r
and peak time t p 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 t d, rise time t r and peak time t p increase.
So, based on the requirement, we can include add the open loop poles or zeros to the transfer function.
If a sinusoidal signal is applied as an input to a Linear Time-Invariant LTI system, then it produces the steady state output, which is also a sinusoidal
signal. The input and output sinusoidal signals have the same frequency, but different amplitudes and phase angles.
G(s) = G(jω)
The amplitude of the output sinusoidal signal is obtained by multiplying the amplitude of the input sinusoidal signal and the magnitude of G(jω)
at ω = ω 0.
The phase of the output sinusoidal signal is obtained by adding the phase of the input sinusoidal signal and the phase of G(jω) at ω = ω 0.
Where,
ω 0 = 2πf 0
Here, f 0 is the frequency of the input sinusoidal signal. Similarly, you can follow the same procedure for closed loop control system.
Consider the transfer function of the second order closed loop control system as,
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 64/90
17/01/2018 Control Systems Quick Guide
C(s) ω 2n
T(s) = =
R(s) 2
s 2 + 2δω ns + ω n
2
ωn
T(jω) =
(jω) 2 + 2δω n(jω) + ω 2n
ω 2n ω 2n
⇒ T(jω) = =
− ω 2 + 2jδωω n + ω 2n
ω 2n 1 −
( ω2
ωn
2
+
2jδω
ωn )
1
⇒ T(jω) =
( )()
1−
ω2
2
ωn
+j
2δω
ωn
ω
Let, = u Substitute this value in the above equation.
ωn
1
T(jω) =
(1 − u 2) + j(2δu)
Magnitude of T(jω) is -
1
M = | T(jω) | =
2 2
√(1 − u ) + (2δu) 2
Phase of T(jω) is -
∠T(jω) = − tan − 1
( ) 2δu
1 − u2
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.
−3
dM 1
⇒
du
= −
2 [
(1 − u 2) 2 + (2δu) 2 ] [4u(u
2 2 − 1 + 2δ 2) ]
dM
Substitute, u = u r and du
== 0 in the above equation.
3
1
0= −
2 [
(1 − u 2r ) 2 + (2δu r) 2 ] [4u (u
−2
r
2
r
− 1 + 2δ 2) ]
2
⇒ 4u r(u r − 1 + 2δ 2) = 0
2
⇒ u r − 1 + 2δ 2 = 0
⇒ u 2r = 1 − 2δ 2
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 65/90
17/01/2018 Control Systems Quick Guide
⇒ ur = √1 − 2δ2
ωr
Substitute, u r = in the above equation.
ωn
ω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 M r.
1
Mr =
2
√(1 − u r )2 + (2δu r) 2
2
Substitute, u r = √1 − 2δ2 and 1 − ur = 2δ2 in the above equation.
1
Mr =
√(2δ ) 2 2
+ (2δ 1 − 2δ 2) 2
√
1
⇒ Mr =
2δ 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.
Substitute, u = 0 in M.
1
M= =1
√(1 − 02) 2 + (2δ(0)) 2
Therefore, the magnitude of T(jω) is one at ω = 0.
1
i.e., at ω = ω B, M = 0.707(1) =
√2
1 1
⇒M= =
√2 √(1 − u 2b) 2 + (2δu b) 2
2 2
⇒ 2 = (1 − u b) 2 + (2δ) 2u b
2
Let, u b = x
⇒ 2 = (1 − x) 2 + (2δ) 2x
⇒ x 2 + (4δ 2 − 2)x − 1 = 0
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 66/90
17/01/2018 Control Systems Quick Guide
− (4δ 2 − 2) ± √(4δ2 − 2) 2 + 4
⇒x=
2
x = 1 − 2δ 2 + √(2δ2 − 1) 2 + 1
⇒ x = 1 − 2δ 2 + √(2 − 4δ2 + 4δ4)
2
ωb
Substitute, x = u 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 t r in the time domain transient response.
Magnitude plot
Phase plot
In both the plots, x-axis represents angular frequency logarithmicscale. Whereas, yaxis represents the magnitude linearscale of open loop transfer function
in the magnitude plot and the phase angle linearscale of the open loop transfer function in the phase plot.
M = 20 log | G(jω)H(jω) |
ϕ = ∠G(jω)H(jω)
Constant K 0 20logK 0
Zero at jω 20 20logω 90
origin
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 67/90
17/01/2018 Control Systems Quick Guide
origin 1
jω
Simple 1 + jωr 20 1 1
0 for ω < 0 for ω <
r r
zero
1 1
20 logωr for ω > 90 for ω >
r r
Simple 1 − 20 1 1
0 for ω < 0 for ω <
1 + jωr r r
pole
1 1
− 20 logωr for ω > − 90 or 270 for ω >
r r
Second
( )
40 40 log ω n for ω < ω n 0 for ω < ω n
order ω2 2jδω
ω 2n 1 − +
derivative ωn
2 ωn 20 log (2δω 2n) for ω = ω n 90 for ω = ω n
term
40 log ω for ω > ω n 180 for ω > ω n
( )
order
ω2
− 20 log (2δω 2n) for ω = ω n
2jδω
integral 2
ωn 1 − + ω − 90 for ω = ω n
2
ωn n
term
− 40 log ω for ω > ω n − 180 for ω > ω n
Magnitude M = 20 logK dB
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 68/90
17/01/2018 Control Systems Quick Guide
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 20 logK dB above the 0 dB line. For the negative values of K, the horizontal line will shift 20 logK dB
below the 0 dB line. The Zero degrees line itself is the phase plot for all the positive values of K.
Magnitude M = 20logω dB
Phase angle ϕ = 90 0
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 69/90
17/01/2018 Control Systems Quick Guide
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.
√
Magnitude M = 20 log 1 + ω 2τ 2 dB
1
For ω < , the magnitude is 0 dB and phase angle is 0 degrees.
τ
1
For ω > τ
, the magnitude is 20 logωτ dB and phase angle is 900.
1 1
The magnitude plot is having magnitude of 0 dB upto ω = rad/sec. From ω = rad/sec, it is having a slope of 20 dB/dec. In this case, the phase plot is
τ τ
1
having phase angle of 0 degrees up to ω = τ
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.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 70/90
17/01/2018 Control Systems Quick Guide
Similarly, you can draw the Bode plots for other terms of the open loop transfer function which are given in the table.
Represent the open loop transfer function in the standard time constant form.
Consider the starting frequency of the Bode plot as 1/10th of the minimum corner frequency or 0.1 rad/sec whichever is smaller value and draw
the Bode plot upto 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
10s
G(s)H(s) =
(s + 2)(s + 5)
Let us convert this open loop transfer function into standard time constant form.
10s
G(s)H(s) =
2 ( )( )
s
2
+1 5
s
5
+1
s
⇒ G(s)H(s) =
( )( )
1+
s
2
1+
s
5
So, we can draw the Bode plot in semi log sheet using the rules mentioned earlier.
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.
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.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 71/90
17/01/2018 Control Systems Quick Guide
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
( )
1
M pc
= 20logM pc
Where, M pc is the magnitude at phase cross over frequency. The unit of gain margin GM is dB.
Phase Margin
PM = 180 0 + ϕ 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.
The Polar plot is a plot, which can be drawn between the magnitude and the phase angle of G(jω)H(jω) by varying ω from zero to ∞. The polar graph
sheet is shown in the following figure.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 72/90
17/01/2018 Control Systems Quick Guide
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 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.
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 values 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 values of ω.
For drawing polar plot more clearly, find the magnitude and phase of G(jω)H(jω) by considering the other values of ω.
Example
Consider the open loop transfer function of a closed loop control system.
5
G(s)H(s) =
s(s + 1)(s + 2)
Let us draw the polar plot for this control system using the above rules.
5
G(jω)H(jω) =
jω(jω + 1)(jω + 2)
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 73/90
17/01/2018 Control Systems Quick Guide
5
M=
√ √
ω( ω + 1)( ω 2 + 4)
2
ω
ϕ = − 90 0 − 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.
0 ∞ -90 or 270
∞ 0 -270 or 90
So, the polar plot starts at (∞,−900) and ends at (0,−2700). 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 −1800 or 1800. So, by equating the phase angle of the open loop transfer function to either −1800 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.
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.
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
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 74/90
17/01/2018 Control Systems Quick Guide
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.
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.
The frequency at which the Nyquist plot intersects the negative real axis (phase angle is 1800) is known as the phase cross over frequency. It is
denoted by ω pc.
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.
1
GM =
M pc
Where, M pc 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 = 180 0 + ϕ gc
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.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 75/90
17/01/2018 Control Systems Quick Guide
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.
Here, the capacitor is in series with the resistor R 2 and the output is measured across this combination.
( )
1
V o(s) 1 s+ τ
=
V i(s) α 1
s+ ατ
Where,
τ = R 2C
R1 + R2
α=
R2
1 1
From the transfer function, we can conclude that the lag compensator has one pole at s = − and one zero at s = − . This means, the pole will be
ατ τ
nearer to origin in the pole-zero configuration of the lag compensator.
( )
1
V o(jω) 1 jω + τ
=
V i(jω) α 1
jω + ατ
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
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 76/90
17/01/2018 Control Systems Quick Guide
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.
Here, the capacitor is parallel to the resistor R 1 and the output is measured across resistor $R_2.
( )
V o(s) sτ + 1
=β
V i(s) βsτ + 1
Where,
τ = R 1C
R2
β=
R1 + R2
1 1
From the transfer function, we can conclude that the lead compensator has pole at s = − β
and zero at s = − βτ
.
( )
V o(jω) jωτ + 1
=β
V i(jω) βjωτ + 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.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 77/90
17/01/2018 Control Systems Quick Guide
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.
( )
1
( )
V o(s) sτ 1 + 1 s+
1 τ2
=β
V i(s) βsτ 1 + 1 α 1
s+ ατ 2
We know αβ = 1.
( )( )
1 1
V o(s) s+ s+
τ1 τ2
⇒ =
V i(s) 1 1
s+ s+
βτ 1 ατ 2
Where,
τ 1 = R 1C 1
τ 2 = R 2C 2
Proportional Controller
The proportional controller produces an output, which is proportional to error signal.
u(t) ∝ e(t)
⇒ u(t) = K Pe(t)
U(s) = K PE(s)
U(s)
= KP
E(s)
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 78/90
17/01/2018 Control Systems Quick Guide
Where,
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.
de(t)
u(t) = K D
dt
U(s) = K DsE(s)
U(s)
= K Ds
E(s)
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) = K I∫ e(t)dt
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 79/90
17/01/2018 Control Systems Quick Guide
K IE(s)
U(s) =
s
U(s) KI
=
E(s) s
KI
Therefore, the transfer function of the integral controller is .
s
The block diagram of the unity negative feedback closed loop control system along with the integral controller is shown in the following figure.
de(t)
u(t) = K Pe(t) + K D
dt
U(s) = (K P + K Ds)E(s)
U(s)
= K P + K Ds
E(s)
The block diagram of the unity negative feedback closed loop control system along with the proportional derivative controller is shown in the following
figure.
The proportional derivative controller is used to improve the stability of control system without affecting the steady state error.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 80/90
17/01/2018 Control Systems Quick Guide
( )
U(s) = K P +
KI
s
E(s)
U(s) KI
= KP +
E(s) s
KI
Therefore, the transfer function of proportional integral controller is K P + .
s
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.
de(t)
u(t) = K Pe(t) + K I∫ e(t)dt + K D
dt
(
U(s) = K P +
KI
s
+ K Ds E(s)
)
U(s) KI
= KP + + K Ds
E(s) s
KI
Therefore, the transfer function of the proportional integral derivative controller is K P + + K Ds.
s
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.
The proportional integral derivative controller is used to improve the stability of the control system and to decrease steady state error.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 81/90
17/01/2018 Control Systems Quick Guide
The state space model of Linear Time-Invariant LTI system can be represented as,
Ẋ = AX + BU
Y = CX + DU
The first and the second equations are known as state equation and output equation respectively.
Where,
X and Ẋ are the state vector and the differential state vector respectively.
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.
State Vector
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 model can be obtained from any one of these two mathematical models. Let us now discuss these two methods one by
one.
There are two storage elements inductorandcapacitor 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, v c(t).
From the circuit, the output voltage, v 0(t) is equal to the voltage across capacitor, v c(t).
v 0(t) = v c(t)
di(t)
v i(t) = Ri(t) + L + v c(t)
dt
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 82/90
17/01/2018 Control Systems Quick Guide
1
C∫
v c(t) = i(t)dt
dv c(t) i(t)
=
dt C
State vector, X =
[ ]
i(t)
v c(t)
[ ]
di ( t )
dt
Differential state vector, Ẋ = dv c ( t )
dt
We can arrange the differential equations and output equation into the standard form of state space model as,
[ ] [ ][ []
di ( t ) R 1
]
−L −L 1
dt i(t)
Ẋ = dv c ( t ) = 1 v c(t) +
L
[ vi(t) ]
0 0
dt C
Y = [0 1]
[ ]
i(t)
v c(t)
Where,
[ ] []
R 1
−L −L 1
A= 1
, B= L , C = [0 1 ] and D = [ 0 ]
0 0
C
Y(s) b0
=
U(s) s n + a n−1
n − 1s + . . . + a 1s + a 0
(s n + a n − 1s n − 1 + . . . + a 0)Y(s) = b 0U(s)
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 83/90
17/01/2018 Control Systems Quick Guide
Let
y(t) = x 1
dy(t)
= x 2 = ẋ 1
dt
d 2y(t)
= x 3 = ẋ 2
dt 2
d n − 1y(t)
= x n = ẋ n − 1
dt n − 1
d ny(t)
= ẋ n
dt n
and u(t) = u
Then,
ẋ n + a n − 1x n + . . . + a 1x 2 + a 0x 1 = b 0u
From the above equation, we can write the following state equation.
ẋ n = − a 0x 1 − a 1x 2 − . . . − a n − 1x n + b 0u
y(t) = y = x 1
[]
ẋ 1
ẋ 2
Ẋ = ⋮
ẋ n − 1
ẋ n
[ ][ ] [ ]
0 1 0 … 0 0 x1 0
0 0 1 … 0 0 x2 0
= ⋮ ⋮ ⋮ … ⋮ ⋮ ⋮ + ⋮ [u ]
0 0 0 … 0 1 xn − 1 0
− a0 − a1 − a2 … − an − 2 − an − 1 xn b0
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 84/90
17/01/2018 Control Systems Quick Guide
[]
x1
x2
Y = [1 0 … 0 0] ⋮
xn − 1
xn
Here, D = [0].
Example
Find the state space model for the system having transfer function.
Y(s) 1
= 2
U(s) s +s+1
(s 2 + s + 1)Y(s) = U(s)
d 2y(t) dy(t)
+ + y(t) = u(t)
dt 2 dt
Let
y(t) = x 1
dy(t)
= x 2 = ẋ 1
dt
and u(t) = u
ẋ 2 = − x 1 − x 2 + u
y(t) = y = x 1
Ẋ =
[]ẋ 1
ẋ 2
=
[ 0
−1
1
−1 ][ ] [ ]
x1
x2
+
0
1
[u]
Y = [1 0]
[]
x1
x2
Y(s) b ns n + b n − 1 s n − 1 + . . . + b 1s + b 0
= n
U(s) s + a n − 1 s n − 1 + . . . + a 1s + a 0
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 85/90
17/01/2018 Control Systems Quick Guide
⇒
Y(s)
U(s)
=
( sn + an − 1 sn − 1
1
+ . . . + a 1s + a 0 )(b ns n + b n − 1s n − 1 + . . . + b 1s + b 0)
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) 1
=
U(s) s n + a n−1
n − 1s + . . . + a 1s + a 0
(s n + a n − 1s n − 1 + . . . + a 0)V(s) = U(s)
Let
v(t) = x 1
dv((t)
= x 2 = ẋ 1
dt
d 2v(t)
= x 3 = ẋ 2
dt 2
d n − 1v(t)
= x n = ẋ n − 1
dt n − 1
d nv(t)
= ẋ n
dt n
and u(t) = u
ẋ n = − a 0x 1 − a 1x 2 − . . . − a n − 1x n + u
Consider,
Y(s)
= b ns n + b n − 1 s n − 1 + . . . + b 1s + b 0
V(s)
Y(s) = (b ns n + b n − 1s n − 1 + . . . + b 1s + b 0)V(s)
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 86/90
17/01/2018 Control Systems Quick Guide
By substituting the state variables and y(t) = y in the above equation, will get the output equation as,
y = b nẋ n + b n − 1x n + . . . + b 1x 2 + b 0x 1
y = b n( − a 0x 1 − a 1x 2 − . . . − a n − 1x n + u) + b n − 1x n + . . . + b 1x 2 + b 0x 1
[]
ẋ 1
ẋ 2
Ẋ = ⋮
ẋ n − 1
ẋ n
[ ][ ] [ ]
0 1 0 … 0 0 x1 0
0 0 1 … 0 0 x2 0
= ⋮ ⋮ ⋮ … ⋮ ⋮ ⋮ + ⋮ [u ]
0 0 0 … 0 1 xn − 1 0
− a0 − a1 − a2 … − an − 2 − an − 1 xn b0
[]
x1
x2
Y = [b 0 − b na 0 b 1 − b na 1 ... b n − 2 − b na n − 2 b n − 1 − b na n − 1 ] ⋮
xn − 1
xn
If b n = 0, then,
[]
x1
x2
Y = [b 0 b1 ... bn − 2 b n − 1] ⋮
xn − 1
xn
Ẋ = AX + BU
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 87/90
17/01/2018 Control Systems Quick Guide
Y = CX + DU
Y(s)
⇒ = C(sI − A) − 1B + D
U(s)
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.
Y(s)
= C(sI − A) − 1B
U(s)
Example
Let us calculate the transfer function of the system represented in the state space model as,
[] [ ][ ] [ ]
ẋ 1 −1 −1 x1 1
Ẋ = = + [u]
ẋ 2 1 0 x2 0
Y = [0 1]
[]
x1
x2
Here,
A=
[ −1
1
−1
0
,
] B=
[]
1
0
, C = [0 1] and D = [0]
Y(s)
= C(sI − A) − 1B
U(s)
Y(s)
U(s)
= [0 1]
[ s+1
−1
1
s ] []
−1 1
0
s −1
1 s+1 [ ]
⇒
Y(s)
U(s) [
= 0 1]
1
(s + 1)s − 1( − 1) 0 []
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 88/90
17/01/2018 Control Systems Quick Guide
⇒
Y(s)
=
[0 1]
[] s
1
=
1
U(s) 2
s +s+1 2
s +s+1
Therefore, the transfer function of the system for the given state space model is
Y(s) 1
=
U(s) s2 +s+1
{
x ZIR(t) = e AtX(0) = L − 1 [sI − A] − 1X(0) }
From the above relation, we can write the state transition matrix ϕ(t) as
ϕ(t) = e At = L − 1[sI − A] − 1
So, the zero input response can be obtained by multiplying the state transition matrix ϕ(t) with the initial conditions 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 = t 1 + t 2 , then the corresponding state transition matrix is equal to the multiplication of the two state transition matrices at t = t 1 and t = t 2.
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.
Qc = B [ AB A 2B ... A n − 1B ]
Find the determinant of matrix Q c 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.
Qo = CT[ A TC T (A T) 2C T ... (A T) n − 1C T ]
Find the determinant of matrix Q o and if it is not equal to zero, then the control system is observable.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 89/90
17/01/2018 Control Systems Quick Guide
Example
Let us verify the controllability and observability of a control system which is represented in the state space model as,
ẋ =
[] ẋ 1
ẋ 2 = [ −1
1
−1
0 ][ ] [ ]
x1
x2 +
1
0 [u]
Y = [0 1]
[] x1
x2
Here,
A= [ −1
1
−1
0 , ] B= [] 1
0 , [0 1 ], D = [0] and n=2
Q c = [B AB]
AB = [ ]
−1
1
⇒ Qc = [ ]
1
0
−1
1
| Qc | = 1 ≠ 0
Since the determinant of matrix Q c is not equal to zero, the given control system is controllable.
Qo = CT [ A TC T ]
Here,
AT = [ ]
−1
−1
1
0 and CT = []
0
1
A TC T = [] 1
0
⇒ Qo = [ ]
0 1
1 0
⇒ | Qo | = − 1 ≠0
Since, the determinant of matrix Q o is not equal to zero, the given control system is observable.
https://www.tutorialspoint.com/cgi-bin/printpage.cgi 90/90