Automatic Control ELTL 442
Automatic Control ELTL 442
Control
ELTL 442
01/02/1444
1 REGISTRATION
05/02/1444
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.
Based on some parameters, we can classify the control systems into the following ways.
The following figure shows the block diagram of the open loop control system.
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.
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 These are also called as feedback
systems. control systems.
Inaccurate. Accurate.
Types of Feedback:
1- Positive Feedback
The positive feedback adds the reference input, 𝑹(𝒔) and feedback output. The
following figure shows the block diagram of positive feedback control system.
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,
𝑮
𝑻= (Equation 1)
𝟏−𝑮𝑯
Where,
Where,
NOTE: The derivation of the above transfer function is present in later chapters.
Effects of Feedback:
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
𝝏𝑻
𝐏𝐞𝐫𝐜𝐞𝐧𝐭𝐚𝐠𝐞 𝐜𝐡𝐚𝐧𝐠𝐞 𝐢𝐧 𝐓
𝑺𝑻𝑮 = 𝑻
𝝏𝑮 = (Equation 3)
𝐏𝐞𝐫𝐜𝐞𝐧𝐭𝐚𝐠𝐞 𝐜𝐡𝐚𝐧𝐠𝐞 𝐢𝐧 𝐆
𝑻
𝟏 𝟏
𝑺𝑻𝑮 = (𝟏 + 𝑮𝑯) =
(𝟏 + 𝑮𝑯)𝟐 𝟏 + 𝑮𝑯
So, we got the sensitivity of the overall gain of closed loop control system as the
reciprocal of (𝟏 + 𝑮𝑯). So, Sensitivity may increase or decrease depending on the value of
(𝟏 + 𝑮𝑯).
• If the value of (𝟏 + 𝑮𝑯) 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 (𝟏 + 𝑮𝑯) 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.
Therefore, we have to properly choose the feedback in order to make the control
system stable.
To know the effect of feedback on noise, let us compare the transfer function relations
with and without feedback due to noise signal alone.
Consider an open loop control system with noise signal as shown below.
Consider a closed loop control system with noise signal as shown below.
In the closed loop control system, the gain due to noise signal is decreased by a factor
of (𝟏 + 𝑮𝒂 𝑮𝒃 𝑯) provided that the term (𝟏 + 𝑮𝒂 𝑮𝒃 𝑯) is greater than one.
𝒅𝒊
𝒗𝒊 = 𝑹𝒊 + 𝑳 + 𝒗𝒐
𝒅𝒕
𝒅𝒗𝒐 𝒅𝟐 𝒗𝒐
⇒ 𝒗𝒊 = 𝑹𝑪 + 𝑳𝑪 𝟐 + 𝒗𝒐
𝒅𝒕 𝒅𝒕
𝒅𝟐 𝒗𝒐 𝑹 𝒅𝒗 𝟏 𝟏
⇒ + ( ) 𝒐 + ( ) 𝒗𝒐 = ( ) 𝒗𝒊
𝒅𝒕 𝟐 𝑳 𝒅𝒕 𝑳𝑪 𝑳𝑪
If 𝒙(𝒕) and 𝒚(𝒕) are the input and output of an LTI system, then the corresponding Laplace
transforms are 𝑿(𝒔) and 𝒀(𝒔).
Therefore, the transfer function of LTI system is equal to the ratio of 𝒀(𝒔) and 𝑿(𝒔).
𝒀(𝒔)
𝒊. 𝒆. , 𝑻𝒓𝒂𝒏𝒔𝒇𝒆𝒓 𝑭𝒖𝒏𝒄𝒕𝒊𝒐𝒏 =
𝑿(𝒔)
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 𝑿(𝒔) & output 𝒀(𝒔).
𝒅𝟐 𝒗𝒐 𝑹 𝒅𝒗 𝟏 𝟏
+ ( ) 𝒐 + ( ) 𝒗𝒐 = ( ) 𝒗𝒊
𝒅𝒕 𝟐 𝑳 𝒅𝒕 𝑳𝑪 𝑳𝑪
𝒔𝑹 𝟏 𝟏
𝒔𝟐 𝑽𝒐 (𝒔) + ( ) 𝑽𝒐 (𝒔) + ( ) 𝑽𝒐 (𝒔) = ( ) 𝑽𝒊 (𝒔)
𝑳 𝑳𝑪 𝑳𝑪
𝑹 𝟏 𝟏
⇒ {𝒔𝟐 + ( ) 𝒔 + } 𝑽𝒐 (𝒔) = ( ) 𝑽𝒊 (𝒔)
𝑳 𝑳𝑪 𝑳𝑪
𝟏
𝑽𝒐 (𝒔) 𝑳𝑪
⇒ =
𝑽𝒊 (𝒔) 𝒔𝟐 + ( ) 𝒔 + 𝟏
𝑹
𝑳 𝑳𝑪
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 𝑽𝒊 (𝒔) & an output 𝑽𝒐 (𝒔).
𝒗𝒓 (𝒕) = 𝑹𝒊(𝒕)
𝒅 𝒗𝒄 (𝒕)
𝒊(𝒕) = 𝑪
𝒅𝒕
𝒅𝒗𝒄 (𝒕)
𝒗𝒓 (𝒕) = 𝑹𝑪
𝒅𝒕
𝒅𝒗𝒄 (𝒕)
𝒆(𝒕) = 𝑹𝑪 + 𝒗𝒄 (𝒕)
𝒅𝒕
We note that the differential equation of the system is from the first derivative,
because the system is first order.
𝑬(𝒔) 𝑹𝑪𝑺 + 𝟏
=
𝑽𝑪 (𝒔) 𝟏
𝑽𝑪 (𝒔) 𝟏
𝑮(𝒔) = =
𝑬(𝒔) 𝑹𝑪𝑺 + 𝟏
We note that the largest exponent in the denominator of the transfer function is equal
to one, because the system is first order.
1- Differential Equation:
2- Transfer Function:
𝑽𝑪 (𝒔) 𝟏
𝑮(𝒔) = =
𝑬(𝒔) 𝑹. 𝑪. 𝑺 + 𝟏
𝑽𝒄 (𝒔) 𝟏 𝟏
𝑮(𝒔) = = =
𝑬(𝒔) (𝟏 × 𝟏𝟎𝟑 )(𝟏𝟎𝟎 × 𝟏𝟎−𝟔 ). 𝑺 + 𝟏 𝟎. 𝟏𝑺 + 𝟏
𝒗𝒓 (𝒕) = 𝑹𝒊(𝒕)
𝒅𝒊(𝒕)
𝒗𝒍 (𝒕) = 𝑳
𝒅𝒕
𝒅𝒗𝒄 (𝒕)
𝒊(𝒕) = 𝑪
𝒅𝒕
𝒅𝒗𝒄 (𝒕)
𝒗𝒓 (𝒕) = 𝑹𝑪
𝒅𝒕
𝒅 𝒗𝒄 (𝒕) 𝒅𝟐 𝒗𝒄 (𝒕)
𝒆(𝒕) = 𝑹𝑪 + 𝑳𝑪 + 𝒗𝒄 (𝒕)
𝒅𝒕 𝒅𝒕𝟐
𝒅𝟐 𝒗𝒄 (𝒕) 𝒅 𝒗𝒄 (𝒕)
𝒆(𝒕) = 𝑳𝑪 + 𝑹𝑪 + 𝒗𝒄 (𝒕)
𝒅𝒕𝟐 𝒅𝒕
𝒆(𝒕) = 𝑳𝑪 𝒗′′ ′
𝒄 (𝒕) + 𝑹𝑪 𝒗𝒄 (𝒕) + 𝒗𝒄 (𝒕)
We note that the differential equation of the system is from the second derivative,
because the system is second order.
𝑽𝑪 (𝒔) 𝟏
𝑮(𝒔) = =
𝑬(𝒔) 𝑳𝑪𝑺𝟐 + 𝑹𝑪𝑺 + 𝟏
We note that the largest exponent in the denominator of the transfer function is equal
to two, because the system is second order.
Example:
1- Differential Equation:
𝒆(𝒕) = 𝑳𝑪 𝒗′′ ′
𝒄 (𝒕) + 𝑹. 𝑪 𝒗𝒄 (𝒕) + 𝒗𝒄 (𝒕)
𝒆(𝒕) = 𝟎. 𝟏𝒗′′
𝒄 (𝒕) + 𝟎. 𝟐𝒗𝒄̀ (𝒕) + 𝒗𝒄 (𝒕)
2- Transfer Function:
𝑽𝑪 (𝒔) 𝟏
𝑮(𝒔) = =
𝑬(𝒔) 𝑳𝑪𝑺𝟐 + 𝑹𝑪𝑺 + 𝟏
𝑽𝑪 (𝒔) 𝟏 𝟏
𝑮(𝒔) = = =
𝑬(𝒔) (𝟐𝟎)(𝟓 × 𝟏𝟎−𝟑 ). 𝑺𝟐 + (𝟒𝟎)(𝟓 × 𝟏𝟎−𝟑 ). 𝑺 + 𝟏 𝟎. 𝟏𝑺𝟐 + 𝟎. 𝟐𝑺 + 𝟏
The basic elements of a block diagram are a block, the summing point and the take-off
point. Let us consider the block diagram of a closed loop control system as shown in the
following figure to identify these elements.
The above block diagram consists of two blocks having transfer functions 𝑮(𝒔) and
𝑯(𝒔). 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 𝑿(𝒔), output 𝒀(𝒔) and the transfer
function 𝑮(𝒔).
𝒀(𝒔)
Transfer Function, 𝑮(𝒔) =
𝑿(𝒔)
⇒ 𝒀(𝒔) = 𝑮(𝒔)𝑿(𝒔)
Output of the block is obtained by multiplying transfer function of the block with input.
The summing point is represented with a circle having cross (𝑿) 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.
𝒊. 𝒆. , 𝒀 = 𝑨 + 𝑩.
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.
𝒀 = 𝑨 + (−𝑩) = 𝑨 − 𝑩.
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
𝒀 = 𝑨 + 𝑩 + (−𝑪) = 𝑨 + 𝑩 − 𝑪.
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, 𝑹(𝒔) to two
more blocks.
In the following figure, the take-off point is used to connect the output 𝑪(𝒔), as one of
the inputs to the summing point.
In this section, let us represent an electrical system with a block diagram. Electrical
systems contain mainly three basic elements — resistor, inductor and capacitor.
Consider a series of RLC circuit as shown in the following figure. Where, 𝑽𝒊 (𝒕) and
𝑽𝒐 (𝒕) are the input and output voltages. Let 𝒊(𝒕) be the current passing through the circuit.
This circuit is in time domain.
𝑽𝒊 (𝒔) − 𝑽𝒐 (𝒔)
𝑰(𝒔) =
𝑹 + 𝒔𝑳
𝟏
⇒ 𝑰(𝒔) = { } {𝑽𝒊 (𝒔) − 𝑽𝒐 (𝒔)} (Equation 1)
𝑹+𝒔𝑳
𝟏
𝑽𝒐 (𝒔)= ( ) 𝑰(𝒔) (Equation 2)
𝒔𝑪
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).
𝟏
Equation 1 can be implemented with a block having the transfer function, . The
𝑹+𝒔𝑳
input and output of this block are {𝑽𝒊 (𝒔) − 𝑽𝒐 (𝒔)} and 𝑰(𝒔). We require a summing point to
get {𝑽𝒊 (𝒔) − 𝑽𝒐 (𝒔)}. The block diagram of Equation 1 is shown in the following figure.
𝟏
Equation 2 can be implemented with a block having transfer function, . The input
𝒔𝑪
and output of this block are 𝑰(𝒔) and 𝑽𝒐 (𝒔). 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.
• 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).
1- Series Connection:
Series connection is also called cascade connection. In the following figure, two blocks
having transfer functions 𝑮𝟏 (𝒔). and 𝑮𝟐 (𝒔). are connected in series.
Compare this equation with the standard form of the output equation, 𝒀(𝒔) = 𝑮(𝒔)𝑿(𝒔).
Where, 𝑮(𝒔) = 𝑮𝟏 (𝒔)𝑮𝟐 (𝒔).
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.
Compare this equation with the standard form of the output equation, 𝒀(𝒔) = 𝑮(𝒔)𝑿(𝒔).
Where, 𝑮(𝒔) = 𝑮𝟏 (𝒔) + 𝑮𝟐 (𝒔).
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.
3- 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.
𝒀(𝒔) 𝑮(𝒔)
⇒ = 𝟏 + 𝑮(𝒔)𝑯(𝒔)
𝑿(𝒔) 𝟏 + 𝑮(𝒔)𝑯(𝒔)
𝑮(𝒔)
Therefore, the negative feedback closed loop transfer function is
𝟏+𝑮(𝒔)𝑯(𝒔)
This means we can represent the negative feedback connection of two blocks with a
single block. The transfer function of this single block is the closed loop transfer function of
the negative feedback. The equivalent block diagram is shown below.
Similarly, you can represent the positive feedback connection of two blocks with a
single block. The transfer function of this single block is the closed loop transfer function of
𝑮(𝒔)
the positive feedback, i.e.,
𝟏−𝑮(𝒔)𝑯(𝒔)
There are two possibilities of shifting summing points with respect to blocks:
Consider the block diagram shown in the following figure. Here, the summing point is
present before the block.
Summing point has two inputs 𝑹(𝒔) and 𝑿(𝒔). The output of it is {𝑹(𝒔) + 𝑿(𝒔)}.
So, the input to the block 𝑮(𝒔) is {𝑹(𝒔) + 𝑿(𝒔)} 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 ‘𝑮(𝒔)𝑹(𝒔)’ 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 𝑮(𝒔). It is
having the input 𝑿(𝒔) and the output of this block is given as input to summing point instead
of 𝑿(𝒔). This block diagram is shown in the following figure.
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 ‘𝑮(𝒔)𝑹(𝒔)′ 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 more block . It is
G(s)
having the input 𝑿(𝒔) and the output of this block is given as input to summing point instead
of 𝑿(𝒔). This block diagram is shown in the following figure.
There are two possibilities of shifting the take-off points with respect to blocks:
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.
Consider the block diagram shown in the following figure. Here, the take-off point is
present after the block.
Follow these rules for simplifying (reducing) the block diagram, which is having many
blocks, summing points and take-off points.
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.
Step 2 − Use Rule 3 for blocks 𝑮𝟏 𝑮𝟐 and 𝑯𝟏 . Use Rule 4 for shifting take-off point after the
block 𝑮𝟓 . The modified block diagram is shown in the following figure.
Step 3 − Use Rule 1 for blocks (𝑮𝟑 + 𝑮𝟒 ) and 𝑮𝟓 . The modified block diagram is shown in
the following figure.
Step 4 − Use Rule 3 for blocks (𝑮𝟑 + 𝑮𝟒 )𝑮𝟓 and 𝑯𝟑 . 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.
𝒀(𝒔) 𝑮𝟏 𝑮𝟐 𝑮𝟓 𝟐 (𝑮𝟑 + 𝑮𝟒 )
=
𝑹(𝒔) (𝟏 + 𝑮𝟏 𝑮𝟐 𝑯𝟏 ){𝟏 + (𝑮𝟑 + 𝑮𝟒 )𝑮𝟓 𝑯𝟑 }𝑮𝟓 − 𝑮𝟏 𝑮𝟐 𝑮𝟓 (𝑮𝟑 + 𝑮𝟒 )𝑯𝟐
Note − Follow these steps in order to calculate the transfer function of the block diagram
having multiple inputs.
• Step 1 − Find the transfer function of block diagram by considering one input at a
time and make the remaining inputs as zero.
• Step 2 − Repeat step 1 for remaining inputs.
• Step 3 − Get the overall transfer function by adding all those transfer functions.
The block diagram reduction process takes more time for complicated systems.
Because, we have to draw the (partially simplified) block diagram after each step. So, to
overcome this drawback, use signal flow graphs (representation).
Example1.
Example2.
Example3.
Example5.
Nodes and branches are the basic elements of signal flow graph.
Node is a point which represents either a variable or a signal. There are three types of nodes:
Let us consider the following signal flow graph to identify these nodes.
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.
Let us construct a signal flow graph by considering the following algebraic equations:
𝒚𝟐 = 𝒂𝟏𝟐 𝒚𝟏 + 𝒂𝟒𝟐 𝒚𝟒
𝒚𝟑 = 𝒂𝟐𝟑 𝒚𝟐 + 𝒂𝟓𝟑 𝒚𝟓
𝒚𝟒 = 𝒂𝟑𝟒 𝒚𝟑
𝒚𝟔 = 𝒂𝟓𝟔 𝒚𝟓
There will be six nodes (𝒚𝟏 , 𝒚𝟐 , 𝒚𝟑 , 𝒚𝟒 , 𝒚𝟓 and 𝒚𝟔 ) and eight branches in this signal
flow graph. The gains of the branches are 𝒂𝟏𝟐 , 𝒂𝟐𝟑 , 𝒂𝟑𝟒 , 𝒂𝟒𝟓 , 𝒂𝟓𝟔 , 𝒂𝟒𝟐 , 𝒂𝟓𝟑 and 𝒂𝟑𝟓 .
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 𝒚𝟐 = 𝒂𝟏𝟐 𝒚𝟏 + 𝒂𝟒𝟐 𝒚𝟒 is shown in the following figure.
Step 2 − Signal flow graph for 𝒚𝟑 = 𝒂𝟐𝟑 𝒚𝟐 + 𝒂𝟓𝟑 𝒚𝟓 is shown in the following figure.
Step 3 − Signal flow graph for 𝒚𝟒 = 𝒂𝟑𝟒 𝒚𝟑 is shown in the following figure.
Step 4 − Signal flow graph for 𝒚𝟓 = 𝒂𝟒𝟓 𝒚𝟒 + 𝒂𝟑𝟓 𝒚𝟑 is shown in the following figure.
Step 5 − Signal flow graph for 𝒚𝟔 = 𝒂𝟓𝟔 𝒚𝟓 is shown in the following figure.
Step 6 − Signal flow graph of overall system is shown in the following figure.
Follow these steps for converting a block diagram into its equivalent signal flow graph.
• Represent all the signals, variables, summing points and take-off points of block
diagram as nodes in signal flow graph.
• Represent the blocks of block diagram as branches 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
(but there is no block in between), 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 𝑹(𝒔) and output signal 𝑪(𝒔) of block diagram as input
node 𝑹(𝒔) and output node 𝑪(𝒔) of signal flow graph.
Just for reference, the remaining nodes (𝒚𝟏 to 𝒚𝟗 ) 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 𝑮𝟏 and 𝑮𝟐 .
With the help of Mason’s gain formula (discussed in the next chapter), you can
calculate the transfer function of this signal flow graph. This is the advantage of signal flow
graphs. Here, we no need to simplify (reduce) the signal flow graphs for calculating the
transfer function.
𝑪(𝒔) ∑𝑵𝒊=𝟏 𝑷𝒊 ∆𝒊
𝑻= =
𝑹(𝒔) ∆
Where,
Consider the following signal flow graph in order to understand the basic terminology
involved here.
• Path:
It is a traversal of branches from one node to any other node in the direction of branch
arrows. It should not traverse any node more than once.
The path that exists from the input node to the output node is known as forward path.
Examples: 𝒚𝟏 →𝒚𝟐 →𝒚𝟑 →𝒚𝟒 →𝒚𝟓 →𝒚𝟔 and 𝒚𝟏 →𝒚𝟐 →𝒚𝟑 →𝒚𝟓 →𝒚𝟔 .
• Forward Path Gain:
It is obtained by calculating the product of all branch gains of the forward path.
Examples: 𝒂𝒃𝒄𝒅𝒆 is the forward path gain of 𝒚𝟏 →𝒚𝟐 →𝒚𝟑 →𝒚𝟒 →𝒚𝟓 →𝒚𝟔 and 𝒂𝒃𝒈𝒆 is the
forward path gain of 𝒚𝟏 →𝒚𝟐 →𝒚𝟑 →𝒚𝟓 →𝒚𝟔 .
• 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: 𝒃𝒋 is the loop gain of 𝒚𝟐 →𝒚𝟑 →𝒚𝟐 and 𝒈𝒉 is the loop gain of 𝒚𝟑 →𝒚𝟓 →𝒚𝟑 .
These are the loops, which should not have any common node.
Examples: The loops, 𝒚𝟐 →𝒚𝟑 →𝒚𝟐 and 𝒚𝟒 →𝒚𝟓 →𝒚𝟒 are non-touching.
Let us consider the same signal flow graph for finding transfer function.
Higher number of (more than two) non-touching loops are not present in this signal flow
graph. We know,
𝑪(𝒔) ∑𝟐𝒊=𝟏 𝑷𝒊 ∆𝒊 𝑷𝟏 ∆𝟏 + 𝑷𝟐 ∆𝟐
𝑻= = =
𝑹(𝒔) ∆ ∆
If the output of control system for an input varies with respect to time, then it is called
the time response of the control system. The time response consists of two parts.
• 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,
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 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 𝒄(𝒕) = 𝟏𝟎 + 𝟓𝒆−𝒕
Here, the second term 𝟓𝒆−𝒕 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.
The standard test signals are impulse, step, ramp and parabolic. These signals are used to
know the performance of the control systems using time response of the output.
𝜹(𝒕) = 𝟎 for 𝒕 ≠ 𝟎
𝟎+
and ∫𝟎− 𝜹(𝒕) 𝒅𝒕 = 𝟏
𝒖(𝒕) = 𝟏; 𝒕 ≥ 𝟎
𝒖(𝒕) = 𝟎; 𝒕 < 𝟎
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’.
𝒓(𝒕) = 𝒕; 𝒕 ≥ 𝟎
𝒓(𝒕) = 𝟎; 𝒕 < 𝟎
We can write unit ramp signal, 𝒓(𝒕) in terms of unit step signal, 𝒖(𝒕) as
𝒓(𝒕) = 𝒕𝒖(𝒕)
𝒕𝟐
𝒑(𝒕) = ;𝒕 ≥ 𝟎
𝟐
𝒑(𝒕) = 𝟎 ; 𝒕 < 𝟎
We can write unit parabolic signal, 𝒑(𝒕) in terms of the unit step signal, 𝒖(𝒕) as,
𝒕𝟐
𝒑(𝒕) = 𝒖(𝒕)
𝟐
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,
𝑪(𝒔) 𝑮(𝒔)
=
𝑹(𝒔) 𝟏 + 𝑮(𝒔)
𝟏
Substitute, 𝑮(𝒔) = in the above equation.
𝝉𝑺
𝟏
𝑪(𝒔) 𝟏
= 𝑺𝝉 =
𝑹(𝒔) 𝟏 + 𝟏 𝝉𝑺 + 𝟏
𝝉𝑺
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.
𝟏
𝑪(𝒔) = ( ) 𝑹(𝒔)
𝝉𝑺 + 𝟏
Where,
Follow these steps to get the response (output) of the first order system in the time domain.
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.
Consider the unit impulse signal as an input to the first order system. So, 𝒓(𝒕) = 𝜹(𝒕)
Rearrange the above equation in one of the standard forms of Laplace transforms.
𝟏 𝟏 𝟏
𝑪(𝒔) = ⇒ 𝑪(𝒔) = ( )
𝟏 𝝉 𝑺+𝟏
𝝉 (𝑺 + )
𝝉 𝝉
𝟏 (− 𝒕 )
𝒄(𝒕) = 𝒆 𝝉 𝒖(𝒕)
𝝉
The unit impulse response is shown in the following figure (Exponential Function).
Consider the unit step signal as an input to first order system. So, 𝒓(𝒕) = 𝒖(𝒕)
𝟏
𝑹(𝒔) =
𝑺
𝟏
Consider the equation, 𝑪(𝒔) = ( ) 𝑹(𝒔)
𝝉𝑺+𝟏
𝟏
Substitute, 𝑹(𝒔) = in the above equation.
𝑺
𝟏 𝟏 𝟏
𝑪(𝒔) = ( )( ) =
𝝉𝑺 + 𝟏 𝑺 𝑺(𝝉𝑺 + 𝟏)
𝟏 𝑨 𝑩 𝑨(𝝉𝑺 + 𝟏) + 𝑩𝑺
𝑪(𝒔) = = + =
𝑺(𝝉𝑺 + 𝟏) 𝑺 𝝉𝑺 + 𝟏 𝑺(𝝉𝑺 + 𝟏)
On both the sides, the denominator term is the same. So, they will get cancelled by
each other. Hence, equate the numerator terms.
𝟏 = 𝑨(𝝉𝑺 + 𝟏) + 𝑩𝑺
By equating the constant terms on both the sides, you will get 𝑨 = 𝟏.
Substitute, 𝑨 = 𝟏 and equate the coefficient of the S terms on both the sides.
𝟎 = 𝝉 + 𝑩 ⇒ 𝑩 = −𝝉
The value of the unit step response, c(t) is zero at 𝒕 = 𝟎 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.
Consider the unit ramp signal as an input to the first order system. So, 𝒓(𝒕) = 𝒕𝒖(𝒕)
𝟏
𝑹(𝒔) =
𝑺𝟐
𝟏
Consider the equation, 𝑪(𝒔) = ( ) 𝑹(𝒔)
𝝉𝑺+𝟏
𝟏
Substitute, 𝑹(𝒔) = in the above equation.
𝑺𝟐
𝟏 𝟏 𝟏
𝑪(𝒔) = ( ) ( 𝟐) = 𝟐
𝝉𝑺 + 𝟏 𝑺 𝑺 (𝝉𝑺 + 𝟏)
On both the sides, the denominator term is the same. So, they will get cancelled by
each other. Hence, equate the numerator terms.
By equating the constant terms on both the sides, you will get 𝑨 = 𝟏.
Substitute, 𝑨 = 𝟏 and equate the coefficient of the S terms on both the sides.
𝟎 = 𝝉 + 𝑩 ⇒ 𝑩 = −𝝉
Similarly, substitute 𝑩 = −𝝉 and equate the coefficient of 𝑺𝟐 terms on both the sides. You
will get 𝑪 = 𝝉𝟐.
Substitute 𝑨 = 𝟏, 𝑩 = −𝝉 and 𝑪 = 𝝉𝟐 in the partial fraction expansion of 𝑪(𝒔).
𝟏 𝝉 𝝉𝟐 𝟏 𝝉 𝝉𝟐 𝟏 𝝉 𝝉
𝑪(𝒔) = 𝟐 − + = 𝟐− + = 𝟐− +
𝑺 𝑺 𝝉𝑺 + 𝟏 𝑺 𝑺 𝝉 (𝑺 + 𝟏) 𝑺 𝑺 𝑺+𝟏
𝝉 𝝉
The unit ramp response, 𝒄(𝒕) has both the transient and the steady state terms.
𝒕𝟐
Consider the unit parabolic signal as an input to the first order system. So, 𝒓(𝒕) = 𝒖(𝒕)
𝟐
𝟏
𝑹(𝒔) =
𝑺𝟑
𝟏
Consider the equation, 𝑪(𝒔) = ( ) 𝑹(𝒔)
𝝉𝑺+𝟏
𝟏
Substitute 𝑹(𝒔) = in the above equation.
𝑺𝟑
𝟏 𝟏 𝟏
𝑪(𝒔) = ( ) ( 𝟑) = 𝟑
𝝉𝑺 + 𝟏 𝑺 𝑺 (𝝉𝑺 + 𝟏)
After simplifying, you will get the values of 𝑨, 𝑩, 𝑪 and 𝑫 as 1, −𝝉, 𝝉𝟐 and
−𝝉𝟑 respectively. Substitute these values in the above partial fraction expansion of 𝑪(𝒔).
𝟏 𝝉 𝝉𝟐 𝝉𝟑 𝟏 𝝉 𝝉𝟐 𝝉𝟐
𝑪(𝒔) = 𝟑 − 𝟐 + − = − + −
𝑺 𝑺 𝑺 𝝉𝑺 + 𝟏 𝑺𝟑 𝑺𝟐 𝑺 𝑺 + 𝟏
𝝉
Apply inverse Laplace transform on both the sides.
𝒕𝟐 𝒕
−( )
𝒄(𝒕) = ( − 𝝉𝒕 + 𝝉𝟐 − 𝝉𝟐 𝒆 𝝉 ) 𝒖(𝒕)
𝟐
The unit parabolic response, 𝒄(𝒕) has both the transient and the steady state terms.
𝒕𝟐
𝒄𝒔𝒔 (𝒕) = ( − 𝝉𝒕 + 𝝉𝟐 ) 𝒖(𝒕)
𝟐
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.
𝑲=𝟓
𝑽𝒊𝒏 = 𝒓 = 𝟐 𝑽
𝑹𝑪 = 𝝉 = (𝟓𝟎)(𝟏𝟎 × 𝟏𝟎−𝟑 ) = 𝟎. 𝟓 𝒔
𝒄(𝒕)𝒆−𝒕⁄𝝉 = 𝟎
• Transfer Function:
𝑪(𝒔) 𝑽𝒐 (𝒔) 𝑲 𝑲 5×2 10
𝑮(𝒔) = = = = = =
𝑹(𝒔) 𝑽𝒊𝒏 (𝒔) 𝝉𝑺 + 𝟏 𝑹𝑪𝑺 + 𝟏 (50)(10 × 10−3 )𝑆 + 1 0.5𝑆 + 1
• Differential Equation:
𝑲 𝒓(𝒕) = 𝝉𝒄̀ (𝒕) + 𝒄(𝒕)
We know that the transfer function of the closed loop control system having unity negative
feedback as
𝑪(𝒔) 𝑮(𝒔)
=
𝑹(𝒔) 𝟏 + 𝑮(𝒔)
𝝎𝟐𝒏
Substitute, 𝑮(𝒔) = in the above equation.
𝑺(𝑺+𝟐𝜹𝝎𝒏 )
𝝎𝟐𝒏
𝑪(𝒔) ( ) 𝝎𝟐𝒏
𝑺(𝑺 + 𝟐𝜻𝝎𝒏 )
= = 𝟐
𝑹(𝒔) 𝝎𝟐𝒏 𝑺 + 𝟐𝜻𝝎𝒏 𝑺 + 𝝎𝟐𝒏
𝟏+( )
𝑺(𝑺 + 𝟐𝜻𝝎𝒏 )
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.
𝑺𝟐 + 𝟐𝜻𝝎𝒏 𝑺 + 𝝎𝟐𝒏 = 𝟎
Where,
Follow these steps to get the response (output) of the second order system in the time
domain.
Consider the unit step signal as an input to the second order system.
𝟏
𝑹(𝒔) =
𝑺
We know the transfer function of the second order closed loop control system is,
𝑪(𝒔) 𝝎𝟐𝒏
=
𝑹(𝒔) 𝑺𝟐 + 𝟐𝜻𝝎𝒏 𝑺 + 𝝎𝟐𝒏
• Case 1: 𝜻 = 𝟎
Substitute, 𝜻 = 𝟎 in the transfer function.
𝑪(𝒔) 𝝎𝟐𝒏 𝝎𝟐𝒏
= ⇒ 𝑪(𝒔) = ( 𝟐 ) 𝑹(𝒔)
𝑹(𝒔) 𝑺𝟐 + 𝝎𝟐𝒏 𝑺 + 𝝎𝟐𝒏
𝟏
Substitute, 𝑹(𝒔) = in the above equation.
𝑺
𝝎𝟐𝒏 𝟏 𝝎𝟐𝒏
𝑪(𝒔) = ( 𝟐 )( ) =
𝑺 + 𝝎𝟐𝒏 𝑺 𝑺(𝑺𝟐 + 𝝎𝟐𝒏 )
So, the unit step response of the second order system when 𝜻= 𝟎 will be a continuous
time signal with constant amplitude and frequency.
• Case 2: 𝜻 = 𝟏
Substitute, 𝜻= 𝟏 in the transfer function.
𝑪(𝒔) 𝝎𝟐𝒏 𝝎𝟐𝒏
= ⇒ 𝑪(𝒔) = ( ) 𝑹(𝒔)
𝑹(𝒔) 𝑺𝟐 + 𝟐𝝎𝒏 𝑺 + 𝝎𝟐𝒏 (𝑺 + 𝝎𝒏 )𝟐
𝟏
Substitute, 𝑹(𝒔) = in the above equation.
𝑺
𝝎𝟐𝒏 𝟏 𝝎𝟐𝒏
𝑪(𝒔) = ( )( ) =
(𝑺 + 𝝎𝒏 )𝟐 𝑺 𝑺(𝑺 + 𝝎𝒏 )𝟐
After simplifying, you will get the values of 𝑨, 𝑩 and 𝑪 as 𝟏, −𝟏 and −𝝎𝒏 respectively.
Substitute these values in the above partial fraction expansion of 𝑪(𝒔).
𝟏 𝟏 𝝎𝒏
𝑪(𝒔) = − −
𝑺 𝑺 + 𝝎𝒏 (𝑺 + 𝝎𝒏 )𝟐
So, the unit step response of the second order system will try to reach the step input in
steady state.
𝟏
Substitute, 𝑹(𝒔) = in the above equation.
𝑺
𝝎𝟐𝒏 𝟏 𝝎𝟐𝒏
𝑪(𝒔) = ( )( ) =
(𝑺 + 𝜻𝝎𝒏 )𝟐 + 𝝎𝟐𝒏 (𝟏 − 𝜻𝟐 ) 𝑺 𝑺((𝑺 + 𝜻𝝎𝒏 )𝟐 + 𝝎𝟐𝒏 (𝟏 − 𝜻𝟐 ))
After simplifying, you will get the values of 𝑨, 𝑩 and 𝑪 as 𝟏, −𝟏 and −𝟐𝜻𝝎𝒏 respectively.
Substitute these values in the above partial fraction expansion of 𝑪(𝒔).
𝟏 𝑺 + 𝟐𝜻𝝎𝒏
𝑪(𝒔) = −
𝑺 (𝑺 + 𝜻𝝎𝒏 )𝟐 + 𝝎𝟐𝒏 (𝟏 − 𝜻𝟐 )
𝟏 (𝑺 + 𝜻𝝎𝒏 ) 𝜻 𝝎𝒏 √(𝟏 − 𝜻𝟐 )
𝑪(𝒔) = − 𝟐− ( 𝟐)
𝑺 √(𝟏 − 𝜻𝟐 ) (𝑺 + 𝜻𝝎 )𝟐 + (𝝎 √(𝟏 − 𝜻𝟐 ))
(𝑺 + 𝜻𝝎𝒏 )𝟐 + (𝝎𝒏 √(𝟏 − 𝜻𝟐 )) 𝒏 𝒏
𝟏 (𝑺 + 𝜻𝝎𝒏 ) 𝜻 𝝎𝒅
𝑪(𝒔) = − − ( )
𝑺 (𝑺 + 𝜻𝝎𝒏 )𝟐 + (𝝎𝒅 )𝟐 √(𝟏 − 𝜻𝟐 ) (𝑺 + 𝜻𝝎𝒏 )𝟐 + (𝝎𝒅 )𝟐
𝜻
𝒄(𝒕) = (𝟏 − 𝒆−𝜻𝝎𝒏 𝒕 𝐜𝐨𝐬(𝝎𝒅 𝒕) − 𝒆−𝜻𝝎𝒏 𝒕 𝐬𝐢𝐧(𝝎𝒅 𝒕)) 𝒖(𝒕)
√(𝟏 − 𝜻𝟐 )
𝒆−𝜻𝝎𝒏 𝒕
𝒄(𝒕) = (𝟏 − ((√(𝟏 − 𝜻𝟐 )) 𝐜𝐨𝐬(𝝎𝒅 𝒕) + 𝜻𝐬𝐢𝐧(𝝎𝒅 𝒕))) 𝒖(𝒕)
√(𝟏 − 𝜻𝟐 )
If √(𝟏 − 𝜻𝟐 ) = 𝐬𝐢𝐧(𝜽), then ‘𝜻’ will be 𝒄𝒐𝒔(𝜽). Substitute these values in the above equation.
𝒆−𝜻𝝎𝒏 𝒕
𝒄(𝒕) = (𝟏 − (𝐬𝐢𝐧(𝜽) 𝐜𝐨𝐬(𝝎𝒅 𝒕) + 𝒄𝒐𝒔(𝜽) 𝐬𝐢𝐧(𝝎𝒅 𝒕))) 𝒖(𝒕)
√(𝟏 − 𝜻𝟐 )
𝒆−𝜻𝝎𝒏 𝒕
= (𝟏 − ( ) 𝐬𝐢𝐧(𝝎𝒅 𝒕 + 𝜽)) 𝒖(𝒕)
√(𝟏 − 𝜻𝟐 )
So, the unit step response of the second order system is having damped oscillations
(decreasing amplitude) when ‘𝜻’ lies between zero and one.
• Case 4: 𝜻 > 𝟏
𝟏
Substitute, 𝑹(𝒔) = in the above equation.
𝑺
𝝎𝟐𝒏 𝟏
𝑪(𝒔) = ( 𝟐 ) ( )
𝟐 𝟐 𝑺
(𝑺 + 𝜻𝝎𝒏 ) − (𝝎𝒏 √(𝜻 − 𝟏))
𝝎𝟐𝒏
=
𝑺 (𝑺 + 𝜻𝝎𝒏 + 𝝎𝒏 √(𝜻𝟐 − 𝟏)) (𝑺 + 𝜻𝝎𝒏 − 𝝎𝒏 √(𝜻𝟐 − 𝟏))
𝟏 𝟏 𝟏
𝑪(𝒔) = +( )( )
𝑺 𝟐 (𝜻 + √(𝜻𝟐 − 𝟏)) (√(𝜻𝟐 − 𝟏)) 𝑺 + 𝜻𝝎𝒏 + 𝝎𝒏 √(𝜻𝟐 − 𝟏)
−𝟏 𝟏
−( )( )
𝟐 (𝜻 − √(𝜻𝟐 − 𝟏)) (√(𝜻𝟐 − 𝟏)) 𝑺 + 𝜻𝝎𝒏 − 𝝎𝒏 √(𝜻𝟐 − 𝟏)
Since it is over damped, the unit step response of the second order system when 𝜻 > 𝟏
will never reach step input in the steady state.
The impulse response of the second order system can be obtained by using any one of
these two methods.
• Follow the procedure involved while deriving step response by considering the value
𝟏
of 𝑹(𝒔) as 1 instead of 𝑺.
The following table shows the impulse response of the second order system for 4 cases of the
damping ratio.
𝜻=𝟎 𝝎𝒏 𝐬𝐢𝐧(𝝎𝒏 𝒕)
𝝎𝒏 𝒆−𝜻𝝎𝒏 𝒕
𝟎<𝜻<𝟏 ( ) 𝐬𝐢𝐧(𝝎𝒅 𝒕)
√(𝟏 − 𝜻𝟐 )
𝟐 𝟐
𝝎𝒏 −(𝜻𝝎𝒏 −𝝎𝒏 √(𝜻 −𝟏)𝒕) −(𝜻𝝎𝒏 +𝝎𝒏√(𝜻 −𝟏)𝒕)
𝜻>𝟏 (𝒆 −𝒆 )
𝟐
𝟐√(𝜻 − 𝟏)
( )
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 𝒕𝒅 .
Consider the step response of the second order system for 𝒕 ≥ 𝟎, when ‘𝜻’ lies between
zero and one.
𝒆−𝜻𝝎𝒏 𝒕
𝒄(𝒕) = 𝟏 − ( ) 𝐬𝐢𝐧(𝝎𝒅 𝒕 + 𝜽)
√𝟏 − 𝜻𝟐
Therefore, at 𝒕 = 𝒕𝒅 , the value of the step response will be 0.5. Substitute, these values
in the above equation.
𝒆−𝜻𝝎𝒏 𝒕𝒅 𝒆−𝜻𝝎𝒏 𝒕𝒅
𝒄(𝒕𝒅 ) = 𝟎. 𝟓 = 𝟏 − ( ) 𝐬𝐢𝐧(𝝎𝒅 𝒕𝒅 + 𝜽) ⇒ ( ) 𝐬𝐢𝐧(𝝎𝒅 𝒕𝒅 + 𝜽) = 𝟎. 𝟓
√𝟏 − 𝜻𝟐 √𝟏 − 𝜻𝟐
𝟏 + 𝟎. 𝟕𝜻
𝒕𝒅 =
𝝎𝒏
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 𝒕𝒓 .
At 𝒕 = 𝒕𝟏 = 𝟎, 𝒄(𝒕) = 𝟎.
Therefore, at 𝒕 = 𝒕𝟐 , the value of step response is one. Substitute, these values in the
following equation.
𝒆−𝜻𝝎𝒏 𝒕
𝒄(𝒕) = 𝟏 − ( ) 𝐬𝐢𝐧(𝝎𝒅 𝒕 + 𝜽)
√𝟏 − 𝜻𝟐
𝒆−𝜻𝝎𝒏 𝒕𝟐
𝒄(𝒕𝟐 ) = 𝟏 = 𝟏 − ( ) 𝐬𝐢𝐧(𝝎𝒅 𝒕𝟐 + 𝜽)
√𝟏 − 𝜻𝟐
𝒆−𝜻𝝎𝒏 𝒕𝟐
⇒ ( ) 𝐬𝐢𝐧(𝝎𝒅 𝒕𝟐 + 𝜽) = 𝟎
√𝟏 − 𝜻𝟐
⇒ 𝐬𝐢𝐧(𝝎𝒅 𝒕𝟐 + 𝜽) = 𝟎
⇒ 𝝎𝒅 𝒕𝟐 + 𝜽 = 𝝅
𝝅−𝜽
⇒ 𝒕𝟐 =
𝝎𝒅
𝒕𝒓 = 𝒕𝟐 − 𝒕𝟏
𝝅−𝜽
∴ 𝒕𝒓 = −𝟎
𝝎𝒅
−𝟏 𝝅 −𝟏 𝝅
𝝅 − 𝜽 𝝅 − [(𝐜𝐨𝐬 𝜻) (𝟏𝟖𝟎∘ )] 𝝅 − [(𝐜𝐨𝐬 𝜻) (𝟏𝟖𝟎∘ )]
𝒕𝒓 = = =
𝝎𝒅 √𝝎𝟐𝒏 − 𝜶𝟐 𝝎𝒏 √𝟏 − 𝜻𝟐
From above equation, we can conclude that the rise time 𝒕𝒓 and the damped
frequency 𝝎𝒅 are inversely proportional to each other.
It is the time required for the response to reach the peak value for the first time. It is
denoted by 𝒕𝒑 . At 𝒕 = 𝒕𝒑 , the first derivate of the response is zero.
We know the step response of second order system for under-damped case is:
𝒆−𝜻𝝎𝒏 𝒕
𝒄(𝒕) = 𝟏 − ( ) 𝐬𝐢𝐧(𝝎𝒅 𝒕 + 𝜽)
√𝟏 − 𝜻𝟐
𝒅𝒄(𝒕)
Substitute, 𝒕 = 𝒕𝒑 and = 𝟎 in the above equation.
𝒅𝒕
𝒆−𝜻𝝎𝒏 𝒕𝒑
𝟎 = −( ) [𝝎𝒅 𝐜𝐨𝐬(𝝎𝒅 𝒕𝒑 + 𝜽) − 𝜻𝝎𝒏 𝒆−𝜻𝝎𝒏 𝒕𝒑 𝐬𝐢𝐧(𝝎𝒅 𝒕𝒑 + 𝜽)]
√𝟏 − 𝜻𝟐
⇒ √𝟏 − 𝜻𝟐 𝐜𝐨𝐬(𝝎𝒅 𝒕𝒑 + 𝜽) −𝜻𝐬𝐢𝐧(𝝎𝒅 𝒕𝒑 + 𝜽) = 𝟎
⇒ 𝐬𝐢𝐧(𝜽 − 𝝎𝒅 𝒕𝒑 − 𝜽) = 𝟎
⇒ 𝝎𝒅 𝒕𝒑 = 𝝅
𝝅 𝝅 𝝅
𝒕𝒑 = = =
𝝎𝒅 √𝝎𝟐𝒏 − 𝜶𝟐 𝝎𝒏 √𝟏 − 𝜻𝟐
From the above equation, we can conclude that the peak time 𝒕𝒑 and the damped
frequency 𝝎𝒅 are inversely proportional to each other.
Peak overshoot 𝑴𝒑 is defined as the deviation of the response at peak time from the final
value of response. It is also called the maximum overshoot.
𝑴𝒑 = 𝒄(𝒕𝒑 ) − 𝒄(∞)
Where,
𝝅
Substitute, 𝒕𝒑 = in the right hand side of the above equation.
𝝎𝒅
𝝅 𝜻𝝅
−𝜻𝝎𝒏 ( ) −( )
𝒆 𝝎𝒅 𝝅 𝒆 √𝟏−𝜻𝟐
𝒄(𝒕𝒑 ) = 𝟏 − ( ) 𝐬𝐢𝐧 (𝝎𝒅 ( ) + 𝜽) = 𝟏 − (− 𝐬𝐢𝐧(𝜽))
√𝟏 − 𝜻𝟐 𝝎𝒅 √𝟏 − 𝜻𝟐
( )
We know that
𝐬𝐢𝐧(𝜽) = √𝟏 − 𝜻𝟐
Substitute the values of 𝒄(𝒕𝒑 ) and 𝒄(∞) in the peak overshoot equation.
𝜻𝝅 𝜻𝝅
−( ) −( )
𝑴𝒑 = 𝟏 + 𝒆 √𝟏−𝜻𝟐 −𝟏=𝒆 √𝟏−𝜻𝟐
By substituting the values of 𝑴𝒑 and 𝒄(∞) in above formula, we will get the Percentage of
the peak overshoot %𝑴𝒑 as
𝜻𝝅
−( )
%𝑴𝒑 = (𝒆 √𝟏−𝜻𝟐 ) × 𝟏𝟎𝟎% = (𝒆−𝜶𝒕𝒑 ) × 𝟏𝟎𝟎% = (𝒆−𝜻𝝎𝒏 𝒕𝒑 ) × 𝟏𝟎𝟎%
∴ 𝑴𝒑 = %𝑴𝒑 ⋅ 𝑪𝒔𝒔
From the above equation, we can conclude that the percentage of peak
overshoot %𝑴𝒑 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 𝟐% and
𝟓%. The settling time is denoted by 𝒕𝒔 .
𝟑
𝒕𝒔 = = 𝟑𝝉
𝜻𝝎𝒏
𝟒
𝒕𝒔 = = 𝟒𝝉
𝜻𝝎𝒏
𝟏
Where, 𝝉 is the time constant and is equal to .
𝜻𝝎𝒏
• Both the settling time 𝒕𝒔 and the time constant 𝝉 are inversely proportional to the
damping ratio 𝜻.
• Both the settling time 𝒕𝒔 and the time constant 𝝉 are independent of the system gain.
That means even the system gain changes, the settling time 𝒕𝒔 and time constant 𝝉 will
never change.
system.
We know that the standard form of the transfer function of the second order closed
loop control system as:
𝝎𝟐𝒏
𝑺𝟐 + 𝟐𝜻𝝎𝒏 𝑺 + 𝝎𝟐𝒏
By equating these two transfer functions, we will get the un-damped natural
frequency 𝝎𝒏 as 𝟐 𝒓𝒂𝒅/𝒔𝒆𝒄 and the damping ratio 𝜻 as 𝟎. 𝟓.
We know the formula for damped frequency 𝝎𝒅 as
𝝎𝒅 = 𝝎𝒏 √𝟏 − 𝜻𝟐
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.
Time
domain Formula Substitution of values in Formula Final value
specification
𝟏 + 𝟎. 𝟕𝛅 𝟏 + 𝟎. 𝟕(𝟎. 𝟓)
Delay time 𝒕𝒅 = 𝒕𝒅 = 𝒕𝒅 = 𝟎. 𝟔𝟕𝟓 𝒔𝒆𝒄
𝝎𝒏 𝟐
𝝅−𝜽 𝝅
𝝅 − (𝟑)
Rise time 𝒕𝒓 = 𝒕𝒓 = 𝒕𝒓 = 𝟏. 𝟐𝟎𝟕 𝒔𝒆𝒄
𝝎𝒅 𝟏. 𝟕𝟑𝟐
𝝅 𝝅
Peak time 𝒕𝒑 = 𝒕𝒑 = 𝒕𝒑 = 𝟏. 𝟖𝟏𝟑 𝒔𝒆𝒄
𝝎𝒅 𝟏. 𝟕𝟑𝟐
𝜻𝝅 𝟎.𝟓𝝅
% Peak −( ) −( )
%𝑴𝒑 = (𝒆 √𝟏−𝜻𝟐 ) × 𝟏𝟎𝟎% %𝑴𝒑 = (𝒆 √𝟏−(𝟎.𝟓)𝟐 ) × 𝟏𝟎𝟎% %𝑴𝒑 = 𝟏𝟔. 𝟑𝟐%
overshoot
Settling time
for 2% 𝟒 𝟒
𝒕𝒔 = 𝒕𝒔 = 𝒕𝒔 = 𝟒 𝒔𝒆𝒄
tolerance 𝜻𝝎𝒏 (𝟎. 𝟓)(𝟐)
band
𝒄(𝒕) = 𝒄′ (𝒕) = 𝟎
𝟎, 𝒕<0
𝒓(𝒕) = {
𝟏, 𝒕≥𝟎
• Differential Equation:
𝑲. 𝝎𝒏 𝟐 𝒓(𝒕) = 𝒄′′ (𝒕) + 𝟐𝜶𝒄′ (𝒕) + 𝝎𝒏 𝟐 𝒄(𝒕)
• Transfer Function:
𝑪(𝒔) 𝑽𝒐 (𝒔) 𝑲𝝎𝒏 𝟐 𝑲𝝎𝒏 𝟐 5×9 45
𝑮(𝒔) = = = 𝟐 = = =
𝑹(𝒔) 𝑽𝒊𝒏 (𝒔) 𝑺 + 𝟐𝜶 𝑺 + 𝝎𝒏 𝟐 𝑺𝟐 + 𝟐𝜻𝝎𝒏 𝑺 + 𝝎𝒏 𝟐 𝑆 2 + 2 𝑆 + 9 𝑆 2 + 2 𝑆 + 9
• Damping Coefficient (𝜶):
𝟐𝜶 = 𝟐
2
𝜶= =1
2
𝝎𝒏 = √9 = 3 𝑟𝑎𝑑/𝑠𝑒𝑐
45 45
𝑲𝒔 = 2
= =5
𝜔𝑛 9
𝝅 𝝅 𝝅
𝒕𝑷 = = = = 1.111 𝑠𝑒𝑐
𝝎𝒅 𝝎𝒏 √𝟏 − 𝜻𝟐 1 2
3√1 − (3)
𝟑𝟐. 𝟗
𝑴𝒑 = × 5 =1.645 𝑉
𝟏𝟎𝟎
𝑴𝑷 = 𝑪𝑷 − 𝑪𝒔𝒔 = 6.65 − 5 = 1.65 𝑉
𝒖(𝒕) ∝ 𝒆(𝒕)
⇒ 𝒖(𝒕) = 𝑲𝑷 𝒆(𝒕)
𝑼(𝒔) = 𝑲𝑷 𝑬(𝒔)
𝑼(𝒔)
= 𝑲𝑷
𝑬(𝒔)
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.
The derivative controller produces an output, which is derivative of the error signal.
𝒅𝒆(𝒕)
𝒖(𝒕) = 𝑲𝑫
𝒅𝒕
𝑼(𝒔) = 𝑲𝑫 𝑺 𝑬(𝒔)
𝑼(𝒔)
= 𝑲𝑫 𝑺
𝑬(𝒔)
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.
The integral controller produces an output, which is integral of the error signal.
𝒖(𝒕) = 𝑲𝑰 ∫ 𝒆(𝒕) 𝒅𝒕
𝑲𝑰 𝑬(𝒔)
𝑼(𝒔) =
𝑺
𝑼(𝒔) 𝑲𝑰
=
𝑬(𝒔) 𝑺
𝑲𝑰
Therefore, the transfer function of the integral controller is .
𝑺
The block diagram of the unity negative feedback closed loop control system along
with the integral controller is shown in the following figure.
𝒅𝒆(𝒕)
𝒖(𝒕) = 𝑲𝑷 𝒆(𝒕) + 𝑲𝑫
𝒅𝒕
𝑼(𝒔)
= 𝑲𝑷 + 𝑲𝑫 𝑺
𝑬(𝒔)
The block diagram of the unity negative feedback closed loop control system along
with the proportional derivative controller is shown in the following figure.
𝑲𝑰
𝑼(𝒔) = (𝑲𝑷 + ) 𝑬(𝒔)
𝑺
𝑼(𝒔) 𝑲𝑰
= 𝑲𝑷 +
𝑬(𝒔) 𝑺
𝑲𝑰
Therefore, the transfer function of proportional integral controller is 𝑲𝑷 + .
𝑺
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.
𝒅𝒆(𝒕)
𝒖(𝒕) = 𝑲𝑷 𝒆(𝒕) + 𝑲𝑰 ∫ 𝒆(𝒕) 𝒅𝒕 + 𝑲𝑫
𝒅𝒕
𝑼(𝒔) 𝑲𝑰
= 𝑲𝑷 + + 𝑲𝑫 𝑺
𝑬(𝒔) 𝑺
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.
Example:
• Block Diagram:
• Differential Equation:
𝒅𝒆(𝒕)
𝒖(𝒕) = 𝑲𝑷 𝒆(𝒕) + 𝑲𝑰 ∫ 𝒆(𝒕) 𝒅𝒕 + 𝑲𝑫
𝒅𝒕
𝑑𝑒(𝑡)
⇒ 10 𝑒(𝑡) + 0.5 ∫ 𝑒(𝑡) 𝑑𝑡 + 0.01
𝑑𝑡
• Transfer Function:
𝑼(𝒔) 𝑲𝑰
= 𝑲𝑷 + + 𝑲𝑫 𝑺
𝑬(𝒔) 𝑺
0.5
⇒ 10 + + 0.01 𝑆
𝑆
• Input and Output Signal:
Let us discuss how to find steady state errors for unity feedback and non-unity
feedback control systems one by one.
Consider the following block diagram of closed loop control system, which is having
unity negative feedback.
Where,
We know the transfer function of the unity negative feedback closed loop control system as
The following table shows the steady state errors and the error constants for standard
input signals like unit step, unit ramp & unit parabolic signals.
𝟏 𝑲𝒑 =𝐥𝐢𝐦 𝑮(𝒔)
unit step signal 𝒔→𝟎
𝟏 + 𝑲𝒑
𝟏 𝑲𝒗 =𝐥𝐢𝐦 𝑺𝑮(𝒔)
unit ramp signal 𝒔→𝟎
𝑲𝒗
𝟏 𝑲𝒂 =𝐥𝐢𝐦 𝑺𝟐 𝑮(𝒔)
unit parabolic signal 𝒔→𝟎
𝑲𝒂
Where, 𝑲𝒑 , 𝑲𝒗 and 𝑲𝒂 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 𝒓(𝒕) = (𝟓 + 𝟐𝒕 + ) 𝒖(𝒕) of
𝟐
𝟓(𝑺+𝟒)
unity negative feedback control system with 𝑮(𝒔) =
𝑺𝟐 (𝑺+𝟏)(𝑺+𝟐𝟎)
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.
𝒕𝟐 𝑲𝒂 =𝐥𝐢𝐦 𝑺𝟐 𝑮(𝒔) = 𝟏 𝟏
𝒓𝟑 (𝒕) = 𝒖(𝒕) 𝒔→𝟎 𝒆𝒔𝒔𝟑 = =𝟏
𝟐 𝑲𝒂
We will get the overall steady state error, by adding the above three steady state errors.
⇒ 𝒆𝒔𝒔 = 𝟎 + 𝟎 + 𝟏 = 𝟏
Therefore, we got the steady state error 𝒆𝒔𝒔 as 1 for this example.
Consider the following block diagram of closed loop control system, which is having
non-unity negative feedback.
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.
This block diagram resembles the block diagram of the unity negative feedback closed
loop control system. Here, the single block is having the transfer
𝑮(𝒔)
function instead of 𝑮(𝒔). 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.
• Input signal:
1
𝑹(𝒔) =
𝑆
1
𝑮(𝒔) =
𝑆+1
• Transfer Function:
1 1 1
𝑪(𝒔) 𝑮(𝒔) 𝑆 + 1 𝑆 + 1 𝑆 + 1
𝑻(𝒔) = = = = =
𝑹(𝒔) 𝟏 + 𝑮(𝒔) 1 + ( 1 ) . (1) (𝑆 + 1 + 1) . (1) 𝑆 + 2
𝑆+1 𝑆+1 𝑆+1
1 𝑆+1 1 1 1
= × = × =
𝑆+1 𝑆+2 1 𝑆+2 𝑆+2
1 1
⇒= (1 − )
𝑆 𝑆+2
1 1 1 1
𝒆𝒔𝒔 = 𝐥𝐢𝐦 𝑺 𝑬(𝒔) = lim (𝑆 (1 − )) = lim (1 − )=1− = 0.5
𝒔→𝟎 𝑆→0 𝑆 𝑆+2 𝑆→0 𝑆+2 0+2
• Input signal:
1
𝑹(𝒔) =
𝑆
𝑲𝑷 (𝒔) = 10
1
𝑮(𝒔) =
𝑆+1
• Transfer Function:
1 10
𝑪(𝒔) 𝑮(𝒔) 𝑲𝑷 (𝒔)𝑮(𝒔) (10). ( )
𝑻(𝒔) = = = = 𝑆 + 1 = 𝑆 +1
𝑹(𝒔) 𝟏 + 𝑮(𝒔) 𝟏 + 𝑲𝑷 (𝒔)𝑮(𝒔) 1 + (10). ( 1 ) . (1) 1 + ( 10 ) . (1)
𝑆+1 𝑆+1
10 10
𝑆+1 10 𝑆+1 10
= = 𝑆+1 = × =
𝑆 + 1 + 10 𝑆 + 11 𝑆 + 1 𝑆 + 11 𝑆 + 11
𝑆+1 𝑆+1
1 10
⇒= (1 − )
𝑆 𝑆 + 11
1 10 10 10
𝒆𝒔𝒔 = 𝐥𝐢𝐦 𝑺 𝑬(𝒔) = lim (𝑆 (1 − )) = lim (1 − )=1− = 0.0.091
𝒔→𝟎 𝑆→0 𝑆 𝑆 + 11 𝑆→0 𝑆 + 11 0 + 11
• Input signal:
1
𝑹(𝒔) =
𝑆
2
𝑲𝑰 (𝒔) =
𝑆
1
𝑮(𝒔) =
𝑆+1
• Transfer Function:
2 1 2
𝑪(𝒔) 𝑮(𝒔) 𝑲𝑰 (𝒔)𝑮(𝒔) ( ).( ) 𝑆 (𝑆 + 1)
𝑻(𝒔) = = = = 𝑆 𝑆+1 =
𝑹(𝒔) 𝟏 + 𝑮(𝒔) 𝟏 + 𝑲𝑰 (𝒔)𝑲𝑭𝑶 (𝒔) 1 + (2) . ( 1 ) . (1) 1 + 2
𝑆 𝑆+1 𝑆 (𝑆 + 1)
2
𝑆 (𝑆 + 1) 2 𝑆 (𝑆 + 1) 2 2
= = × = = 2
𝑆 (𝑆 + 1) + 2 𝑆 (𝑆 + 1) 𝑆 (𝑆 + 1 ) + 2 𝑆 (𝑆 + 1) + 2 𝑆 + 𝑆 + 2
𝑆 (𝑆 + 1)
1 2
⇒= (1 − 2 )
𝑆 𝑆 +𝑆+2
1 2 2 2
𝒆𝒔𝒔 = 𝐥𝐢𝐦 𝑺 𝑬(𝒔) = lim (𝑆 (1 − 2 )) = lim (1 − 2 )=1− =0
𝒔→𝟎 𝑆→0 𝑆 𝑆 +𝑆+2 𝑆→0 𝑆 +𝑆+2 0+0+2
• Input signal:
1
𝑹(𝒔) =
𝑆
10 𝑆 2 + 3 𝑆 + 10
𝑲(𝒔) = 3 + +1𝑆 =
𝑆 𝑆
1
𝑮(𝒔) =
𝑆+1
• Transfer Function:
𝑆 2 + 3 𝑆 + 10 1 𝑆 2 + 3 𝑆 + 10
𝑪(𝒔) 𝑮(𝒔) 𝑲(𝒔)𝑮(𝒔) (
) . ( )
𝑆 𝑆+1 𝑆(𝑆 + 1)
𝑻(𝒔) = = = = 2 =
𝑹(𝒔) 𝟏 + 𝑮(𝒔) 𝟏 + 𝑲(𝒔)𝑮(𝒔) 𝑆 + 3 𝑆 + 10 1 𝑆 2 + 3 𝑆 + 10
1+( ).( ) . (1) 1 +
𝑆 𝑆+1 𝑆(𝑆 + 1)
𝑆 2 + 3 𝑆 + 10 𝑆 2 + 3 𝑆 + 10 𝑆 2 + 3 𝑆 + 10
𝑆(𝑆 + 1) 𝑆(𝑆 + 1) 𝑆(𝑆 + 1)
= = 2 =
𝑆(𝑆 + 1) + 𝑆 2 + 3 𝑆 + 10 𝑆 + 𝑆 + 𝑆 2 + 3 𝑆 + 10 2 𝑆 2 + 4 𝑆 + 10
𝑆(𝑆 + 1) 𝑆(𝑆 + 1) 𝑆(𝑆 + 1)
𝑆 2 + 3 𝑆 + 10 𝑆(𝑆 + 1) 𝑆 2 + 3 𝑆 + 10
= × =
𝑆(𝑆 + 1) 2 𝑆 2 + 4 𝑆 + 10 2 𝑆 2 + 4 𝑆 + 10
1 𝑆2 + 3 𝑆 + 10
⇒= (1 − 2 )
𝑆 2 𝑆 + 4 𝑆 + 10
1 𝑆 2 + 3 𝑆 + 10 𝑆 2 + 3 𝑆 + 10 0 + 0 + 10
= lim (𝑆 (1 − 2
)) = lim (1 − 2
)=1− =0
𝑆→0 𝑆 2 𝑆 + 4 𝑆 + 10 𝑆→0 2 𝑆 + 4 𝑆 + 10 0 + 0 + 10
• Input signal:
1
𝑹(𝒔) =
𝑆
𝑲𝑷 (𝒔) = 10
1
𝑮(𝒔) =
𝑆+2
2
𝑯(𝒔) =
𝑆+4
• Transfer Function:
1 10
𝑪(𝒔) 𝑮(𝒔) 𝑲𝑷 (𝒔)𝑮(𝒔) (10). ( )
𝑻(𝒔) = = = = 𝑆 + 2 = 𝑆 +2
𝑹(𝒔) 𝟏 + 𝑮(𝒔) 𝟏 + 𝑲𝑷 (𝒔)𝑮(𝒔)𝑯(𝒔) 1 2 20
1 + (10). (𝑆 + 2) . (𝑆 + 4) 1 +
(𝑆 + 2)(𝑆 + 4)
10
𝑆 +2 10 (𝑆 + 2)(𝑆 + 4) 10(𝑆 + 4) 10 𝑆 + 40
= = × = = 2
(𝑆 + 2)(𝑆 + 4) + 20 𝑆 + 2 (𝑆 + 2)(𝑆 + 4) + 20 (𝑆 + 2)(𝑆 + 4) + 20 𝑆 + 6 𝑆 + 28
(𝑆 + 2)(𝑆 + 4)
1 10 𝑆 + 40
⇒= (1 − 2 )
𝑆 𝑆 + 6 𝑆 + 28
1 10 𝑆 + 40 10 𝑆 + 40 0 + 40
𝒆𝒔𝒔 = 𝐥𝐢𝐦 𝑺 𝑬(𝒔) = lim (𝑆 (1 − 2 )) = lim (1 − 2 )=1− = −0.428
𝒔→𝟎 𝑆→0 𝑆 𝑆 + 6 𝑆 + 28 𝑆→0 𝑆 + 6 𝑆 + 28 0 + 0 + 28
• Input signal:
𝟏
𝑹(𝒔) =
𝑺
1
𝑮(𝒔) =
𝑆+4
5
𝑯(𝒔) =
𝑆 + 10
• Transfer Function:
1 𝐾
𝑪(𝒔) 𝑮(𝒔) 𝑲(𝒔)𝑮(𝒔) (𝐾). ( )
𝑻(𝒔) = = = = 𝑆 + 4 = 𝑆 + 4
𝑹(𝒔) 𝟏 + 𝑮(𝒔) 𝟏 + 𝑲𝑷 (𝒔)𝑮(𝒔)𝑯(𝒔) 1 + (𝐾). ( 1 ) . ( 5 ) 1 + 5𝐾
𝑆+4 𝑆 + 10 (𝑆 + 4)(𝑆 + 10)
𝐾
𝑆 + 4 𝐾 (𝑆 + 4)(𝑆 + 10) 𝐾(𝑆 + 10)
= = × =
(𝑆 + 4)(𝑆 + 10) + 5 𝐾 𝑆 + 4 (𝑆 + 4)(𝑆 + 10) + 5 𝐾 (𝑆 + 4)(𝑆 + 10) + 5 𝐾
(𝑆 + 4)(𝑆 + 10)
𝐾 𝑆 + 10 𝐾
=
𝑆 2 + 14 𝑆 + 40 + 5 𝐾
1 𝐾 𝑆 + 10 𝐾
⇒= (1 − 2 )
𝑆 𝑆 + 14 𝑆 + 40 + 5 𝐾
1 𝐾 𝑆 + 10 𝐾
= lim (𝑆 (1 − 2 ))
𝑆→0 𝑆 𝑆 + 14 𝑆 + 40 + 5 𝐾
𝐾 𝑆 + 10 𝐾 10 𝐾
= lim (1 − ) = 1 −
𝑆→0 𝑆 2 + 14 𝑆 + 40 + 5 𝐾 40 + 5 𝐾
What is Stability?
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.
If the system is stable for all the range of system component values, then it is known as
the absolutely stable system. The open loop control system is absolutely stable if all the
poles of the open loop transfer function present in left half of ‘S’ plane. Similarly, the closed
loop control system is absolutely stable if all the poles of the closed loop transfer function
present in the left half of the ‘S’ plane.
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.
Routh-Hurwitz stability criterion is having one necessary condition and one sufficient
condition for stability. If any control system doesn’t satisfy the necessary condition, then we
can say that the control system is unstable. But, if the control system satisfies the necessary
condition, then it may or may not be stable. So, the sufficient condition is helpful for
knowing whether the control system is stable or not.
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.
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.
If all the roots of the characteristic equation exist to the left half of the ‘S’ plane, then
the control system is stable. If at least one root of the characteristic equation exists to the
right half of the ‘S’ plane, then the control system is unstable. So, we have to find the roots
of the characteristic equation to know whether the control system is stable or unstable. But, it
is difficult to find the roots of the characteristic equation as order increases.
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 𝑺𝒏 and
continue up to the coefficient of 𝑺𝟎 .
• 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 𝑺𝟎 is 𝒂𝒏 .
Here, 𝒂𝒏 is the coefficient of 𝑺𝟎 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.
𝒂𝟏 𝒂𝟐 − 𝒂𝟑 𝒂𝟎 𝒂𝟏 𝒂𝟒 − 𝒂𝟓 𝒂𝟎 𝒂𝟏 𝒂𝟔 − 𝒂𝟕 𝒂𝟎
𝒏−𝟐 𝒃𝟏 = 𝒃𝟐 = 𝒃𝟑 =
𝑺 ... ... ...
𝒂𝟏 𝒂𝟏 𝒂𝟏
𝒃𝟏 𝒂𝟑 − 𝒃𝟐 𝒂𝟏 𝒃𝟏 𝒂𝟓 − 𝒃𝟑 𝒂𝟏
𝑺𝒏−𝟑 𝒄𝟏 = 𝒄𝟐 = ⋮
𝒃𝟏 𝒃𝟏
⋮ ⋮ ⋮ ⋮
𝑺𝟏 ⋮ ⋮
𝑺𝟎 𝒂𝒏
Example:
Let us find the stability of the control system having characteristic equation,
All the coefficients of the characteristic polynomial, 𝑺𝟒 + 𝟑𝑺𝟑 + 𝟑𝑺𝟐 + 𝟐𝑺+ 𝟏 are positive.
So, the control system satisfies the necessary condition.
Step 2 − Form the Routh array for the given characteristic polynomial.
𝑺𝟒 𝒂𝟎 = 𝟏 𝒂𝟐 = 𝟑 𝒂𝟒 = 𝟏
𝑺𝟑 𝒂𝟏 = 𝟑 𝒂𝟑 = 𝟐 𝒂𝟓 = 𝟎
𝒂𝟏 𝒂𝟔 − 𝒂𝟕 𝒂𝟎
𝒃𝟑 =
(𝟑 × 𝟑) − (𝟐 × 𝟏) (𝟑 × 𝟏) − (𝟎 × 𝟏) 𝒂𝟏
𝒃𝟏 = 𝒃𝟐 =
𝟑 𝟑 (𝟑 × 𝟎) − (𝟎 × 𝟏)
𝑺𝟐 =
𝟕 𝟑 𝟑
= = =𝟏
𝟑 𝟑 𝟎
= =𝟎
𝟑
𝑺𝟎 1
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.
We may come across two types of situations, while forming the Routh table. It is
difficult to complete the Routh table from these two situations.
Example:
Let us find the stability of the control system having characteristic equation,
All the coefficients of the characteristic polynomial, 𝑺𝟒 + 𝟐𝑺𝟑 +𝑺𝟐 + 𝟐𝑺+ 𝟏 are positive. So,
the control system satisfied the necessary condition.
Step 2 − Form the Routh array for the given characteristic polynomial.
𝟐 𝟐
𝑺𝟑 𝒂𝟏 = =𝟏 𝒂𝟑 = =𝟏 𝒂𝟓 = 𝟎
𝟐 𝟐
(𝟏 × 𝟏) − (𝟏 × 𝟏) (𝟏 × 𝟏) − (𝟎 × 𝟏)
𝑺𝟐 𝒃𝟏 = =𝟎 𝒃𝟐 = =𝟏
𝟏 𝟏
𝑺𝟏
𝑺𝟎
The row 𝑺𝟑 elements have 2 as the common factor. So, all these elements are divided by 2.
• Special case (i) − Only the first element of row 𝑺𝟐 is zero. So, replace it by 𝝐 and
continue the process of completing the Routh table.
𝑺𝟒 1 1 1
𝑺𝟑 1 1
𝑺𝟐 𝝐 1
(𝝐× 𝟏) − (𝟏 × 𝟏) 𝝐 − 𝟏
𝑺𝟏 𝒄𝟏 = =
𝝐 𝝐
𝑺𝟎 1
𝑺𝟒 1 1 1
𝑺𝟑 1 1
𝑺𝟐 0 1
𝑺𝟏 −∞
𝑺𝟎 1
There are two sign changes in the first column of Routh table. Hence, the control system
is unstable.
ENG. KHALID H ALSAYED & ENG. FLAIH R ALHARBY 97
2- All the Elements of any row of the Routh array are zero:
• Write the auxiliary equation, 𝑨(𝒔) of the row, which is just above the row of zeros.
• Differentiate the auxiliary equation, 𝑨(𝒔) 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,
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.
𝑺𝟓 𝒂𝟎 = 𝟏 𝒂𝟐 = 𝟏 𝒂𝟒 = 𝟏
𝟑 𝟑 𝟑
𝑺𝟒 𝒂𝟏 = =𝟏 𝒂𝟑 = =𝟏 𝒂𝟓 = =𝟏
𝟑 𝟑 𝟑
(𝟏 × 𝟏) − (𝟏 × 𝟏) (𝟏 × 𝟏) − (𝟏 × 𝟏)
𝑺𝟑 𝒃𝟏 = =𝟎 𝒃𝟐 = =𝟎
𝟏 𝟏
𝑺𝟐
𝑺𝟏
𝑺𝟎
The row 𝑺𝟒 elements have the common factor of 3. So, all these elements are divided by 3.
• Special case (ii) − All the elements of row 𝑺𝟑 are zero. So, write the auxiliary
equation, 𝑨(𝒔) of the row 𝑺𝟒 .
𝑨(𝒔) = 𝑺𝟒 + 𝑺𝟐 + 𝟏
𝒅𝑨(𝒔)
= 𝟒𝑺𝟑 + 𝟐𝑺
𝒅𝒔
𝑺𝟓 𝒂𝟎 = 𝟏 𝒂𝟐 = 𝟏 𝒂𝟒 = 𝟏
𝟑 𝟑 𝟑
𝑺𝟒 𝒂𝟎 = =𝟏 𝒂𝟐 = =𝟏 𝒂𝟒 = =𝟏
𝟑 𝟑 𝟑
𝟒 𝟐
𝑺𝟑 𝒂𝟏 = =𝟐 𝒂𝟑 = =𝟏 𝒂𝟓 = 𝟎
𝟐 𝟐
(𝟐 × 𝟏) − (𝟏 × 𝟏) (𝟐 × 𝟏) − (𝟎 × 𝟏)
𝑺𝟐 𝒃𝟏 = = 𝟎. 𝟓 𝒃𝟐 = =𝟏
𝟐 𝟐
(𝟎. 𝟓 × 𝟏) − (𝟏 × 𝟐) −𝟏. 𝟓
𝑺𝟏 𝒄𝟏 = =
𝟎. 𝟓 𝟎. 𝟓
𝑺𝟎 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.
𝑿̇ = 𝑨𝑿 + 𝑩𝑼
𝒀 = 𝑪𝑿 + 𝑫𝑼
The first and the second equations are known as state equation and output equation
respectively.
Where,
• 𝑿 and 𝑿̇ are the state vector and the differential state vector respectively.
• 𝑼 and 𝒀 are input vector and output vector respectively.
• 𝑨 is the system matrix.
• 𝑩 and 𝑪 are the input and the output matrices.
• 𝑫 is the feed-forward matrix.
• State: 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.
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 (inductor and capacitor) in this circuit. So, the number
of the state variables is equal to two and these state variables are the current flowing through
the inductor, 𝒊(𝒕) and the voltage across capacitor, 𝒗𝑪 (𝒕).
From the circuit, the output voltage, 𝒗𝟎 (𝒕) is equal to the voltage across capacitor, 𝒗𝑪 (𝒕).
𝒗𝟎 (𝒕) = 𝒗𝑪 (𝒕)
𝒅𝒊(𝒕)
𝒗𝒊 (𝒕) = 𝑹𝒊(𝒕) + 𝑳 + 𝒗𝑪 (𝒕)
𝒅𝒕
𝟏
𝒗𝑪 (𝒕) = ∫ 𝒊(𝒕) 𝒅𝒕
𝑪
𝒊(𝒕)
State vector, 𝑿 = [ ]
𝒗𝒄 (𝒕)
𝒅𝒊(𝒕)
𝒅
Differential state vector, 𝑿̇ = [𝒅𝒗 𝒕(𝒕)]
𝑪
𝒅𝒕
We can arrange the differential equations and output equation into the standard form
of state space model as,
ENG. KHALID H ALSAYED & ENG. FLAIH R ALHARBY 101
𝒅𝒊(𝒕) 𝑹 𝟏
− − 𝒊(𝒕) 𝟏
𝒅 𝑳 𝑳
𝑿̇ = 𝒕 =[ ][ ] + [𝑳] [𝒗𝒊 (𝒕)]
𝒅𝒗𝑪 (𝒕) 𝟏 𝒗𝒄 (𝒕)
𝟎 𝟎
[ 𝒅𝒕 ] 𝑪
𝒊(𝒕)
𝒀 = [ 𝟎 𝟏] [ ] + [𝟎][𝒖(𝒕)]
𝒗𝒄 (𝒕)
Where,
𝑹 𝟏
−− 𝟏
𝑨=[ 𝑳 𝑳 ] , 𝑩 = [𝑳] , 𝑪 = [𝟎 𝟏], 𝑫 = [𝟎]
𝟏
𝟎 𝟎
𝑪
Consider the two types of transfer functions based on the type of terms present in the
numerator. We will study one type.
𝒀(𝒔) 𝒃𝟎
= 𝒏
𝑼(𝒔) 𝑺 + 𝒂𝒏−𝟏 𝑺𝒏−𝟏 + ⋯ + 𝒂𝟏 𝑺 + 𝒂𝟎
Let
𝒚(𝒕) = 𝒙𝟏
𝒅𝟐 𝒚(𝒕)
= 𝒙𝟑 = 𝒙̇ 𝟐
𝒅𝒕𝟐
𝒅𝒏−𝟏 𝒚(𝒕)
= 𝒙𝒏 = 𝒙̇ 𝒏−𝟏
𝒅𝒕𝒏−𝟏
𝒅𝒏 𝒚(𝒕)
= 𝒙̇ 𝒏
𝒅𝒕𝒏
and 𝒖(𝒕) = 𝒖
Then,
𝒙̇ 𝒏 + 𝒂𝒏−𝟏 𝒙𝒏 + ⋯ + 𝒂𝟏 𝒙𝟐 + 𝒂𝟎 𝒙𝟏 = 𝒃𝟎 𝒖
From the above equation, we can write the following state equation.
𝒙̇ 𝒏 = −𝒂𝟎 𝒙𝟏 − 𝒂𝟏 𝒙𝟐 − ⋯ − 𝒂𝒏−𝟏 𝒙𝒏 + 𝒃𝟎 𝒖
𝒚(𝒕) = 𝒚 = 𝒙𝟏
𝒙̇ 𝟏 𝟎 𝟏 𝟎 ⋯ 𝟎 𝟎 𝒙𝟏 𝟎
𝒙̇ 𝟐 𝟎 𝟎 𝟏 ⋯ 𝟎 𝟎 𝒙 𝟐 𝟎
̇𝑿 = ⋮ = ⋮ ⋮ ⋮ ⋯ ⋮ ⋮ ⋮ + ⋮ [𝒖]
𝒙̇ 𝒏−𝟏 𝟎 𝟎 𝟎 ⋯ 𝟎 𝟏 𝒙𝒏−𝟏 𝟎
[ 𝒙̇ 𝒏 ] [−𝒂𝟎 −𝒂𝟏 −𝒂𝟐 ⋯ −𝒂𝒏−𝟐 −𝒂𝒏−𝟏 ] [ 𝒙𝒏 ] [𝒃𝟎 ]
𝒙𝟏
𝒙𝟐
𝒀 = [𝟏 𝟎 ⋯ 𝟎 𝟎] ⋮ + [𝟎][𝒖]
𝒙𝒏−𝟏
[ 𝒙𝒏 ]
𝟎 𝟏 𝟎 ⋯ 𝟎 𝟎 𝟎
𝟎 𝟎 𝟏 ⋯ 𝟎 𝟎 𝟎
𝑨= ⋮ ⋮ ⋮ ⋯ ⋮ ⋮ , 𝑩 = ⋮ , 𝑪 = [𝟏 𝟎 ⋯ 𝟎 𝟎 ], 𝑫 = [𝟎 ]
𝟎 𝟎 𝟎 ⋯ 𝟎 𝟏 𝟎
[−𝒂𝟎 −𝒂𝟏 −𝒂𝟐 ⋯ −𝒂𝒏−𝟐 −𝒂𝒏−𝟏 ] [𝒃 𝟎 ]
Example:
Find the state space model for the system having transfer function.
𝒀(𝒔) 𝟏
= 𝟐
𝑼(𝒔) 𝑺 + 𝑺 + 𝟏
𝒅𝟐 𝒚(𝒕) 𝒅𝒚(𝒕)
+ + 𝒚(𝒕) = 𝒖(𝒕)
𝒅𝒕𝟐 𝒅𝒕
Let
𝒚(𝒕) = 𝒙𝟏
𝒅𝒚(𝒕)
= 𝒚̇ 𝟏 = 𝒙𝟐 = 𝒙̇ 𝟏
𝒅𝒕
𝒅𝟐 𝒚(𝒕)
= 𝒚̇ 𝟐 = 𝒙𝟑 = 𝒙̇ 𝟐
𝒅𝒕𝟐
and 𝒖(𝒕) = 𝒖
𝒙𝟑 + 𝒙𝟐 + 𝒙𝟏 = 𝒖
𝒙̇ 𝟐 = −𝟏𝒙𝟏 − 𝟏𝒙𝟐 + 𝟏𝒖
𝒚(𝒕) = 𝒚 = 𝒙𝟏
𝒚 = 𝑪𝒙 + 𝑫𝒖
𝒙𝟏 𝒙𝟐 𝒖𝟏
𝒙̇ = 𝑨 = [𝒙 𝒙𝟐 ] 𝒙 + 𝑩 = [𝒖𝟐 ] 𝒖
𝟏
𝒚 = 𝑪 = [𝒙𝟏 𝒙𝟐 ]
𝒙̇ 𝟎 𝟏 𝒙𝟏 𝟎
𝑿̇ = [ 𝟏 ] = [ ] [𝒙 ] + [ ] [𝒖]
𝒙̇ 𝟐 −𝟏 −𝟏 𝟐 𝟏
𝒙𝟏
𝒀 = [𝟏 𝟎] [𝒙 ] + [𝟎][𝒖]
𝟐
Where,
𝟏 𝟏 𝟎
𝑨=[ ] , 𝑩 = [ ] , 𝑪 = [𝟏 𝟎], 𝑫 = [𝟎]
−𝟏 −𝟏 𝟏
Example:
Find the state space model for the system having transfer function.
𝒀(𝒔) 𝑺𝟐 + 𝟑𝑺 + 𝟖 𝟏
= 𝟑 𝟐
=( 𝟑 𝟐
) (𝑺𝟐 + 𝟑𝑺 + 𝟖)
𝑼(𝒔) 𝑺 + 𝟔𝑺 + 𝟏𝟎𝑺 + 𝟓 𝑺 + 𝟔𝑺 + 𝟏𝟎𝑺 + 𝟓
Let
𝒚(𝒕) = 𝒙𝟏
𝒅𝒚(𝒕)
= 𝒚̇ 𝟏 = 𝒙𝟐 = 𝒙̇ 𝟏
𝒅𝒕
𝒅𝟐 𝒚(𝒕)
= 𝒚̇ 𝟐 = 𝒙𝟑 = 𝒙̇ 𝟐
𝒅𝒕𝟐
𝒅𝟑 𝒚(𝒕)
= 𝒚̇ 𝟑 = 𝒙𝟒 = 𝒙̇ 𝟑
𝒅𝒕𝟑
𝒅𝒖(𝒕)
= 𝒙𝟐
𝒅𝒕
𝒅𝟐 𝒖(𝒕)
= 𝒙𝟑
𝒅𝒕𝟐
𝒙𝟒 + 𝒙𝟑 + 𝒙𝟐 + 𝒙𝟏 = (𝒙𝟑 + 𝒙𝟐 + 𝒙𝟏 )𝒖
𝒚(𝒕) = 𝒚 = 𝒙𝟏
𝒖(𝒕) = 𝒖 = 𝒙𝟏
𝒙̇ = 𝑨𝒙 + 𝑩𝒖
𝒚 = 𝑪𝒙 + 𝑫𝒖
𝒙𝟏 𝒙𝟐 𝒖𝟏
𝒙̇ = 𝑨 = [𝒙 𝒙𝟐 ] 𝒙 + 𝑩 = [𝒖𝟐 ] 𝒖
𝟏
𝒚 = 𝑪 = [𝒙𝟏 𝒙𝟐 ]
𝒙̇ 𝟏 𝟎 𝟏 𝟎 𝒙𝟏 𝟎
̇𝑿 = [𝒙̇ 𝟐 ] = [ 𝟎 𝟎 𝟏 ] [𝒙𝟐 ] + [𝟎] [𝒖]
𝒙̇ 𝟑 −𝟓 −𝟏𝟎 −𝟔 𝒙𝟑 𝟏
𝒙𝟏
𝒀 = [𝟖 𝟑 𝟏] [𝒙𝟐 ] + [𝟎][𝒖]
𝒙𝟑
Where,
𝟎 𝟏 𝟎 𝟎
𝑨=[ 𝟎 𝟎 𝟏 ] , 𝑩 = [𝟎] , 𝑪 = [𝟖 𝟑 𝟏], 𝑫 = [𝟎]
−𝟓 −𝟏𝟎 −𝟔 𝟏
𝒙̇ 𝟏 = 𝒖 − 𝒙𝟐 𝒙̇ 𝟐 = −𝒙𝟐 − 𝟐𝒙𝟏 − 𝒖 𝒚 = 𝒙𝟐 + 𝒙𝟏
𝑿̇ = 𝑨𝑿 + 𝑩𝑼
𝒀 = 𝑪𝑿 + 𝑫𝑼
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.
Let us calculate the transfer function of the system represented in the state space model as,
𝒙̇ −𝟏 −𝟏 𝒙𝟏 𝟏
𝑿̇ = [ 𝟏 ] = [ ] [𝒙 ] + [ ] [𝒖]
𝒙̇ 𝟐 𝟏 𝟎 𝟐 𝟎
𝒙𝟏
𝒀 = [𝟎 𝟏] [𝒙 ] + [𝟎][𝒖]
𝟐
Here,
−𝟏 −𝟏 𝟏
𝑨=[ ] , 𝑩 = [ ] , 𝑪 = [𝟎 𝟏], 𝑫 = [𝟎]
𝟏 𝟎 𝟎
𝒔 + 𝟏 𝟏 −𝟏 𝟏
⇒ [𝟎 ]
𝟏 [ ] [ ]
−𝟏 𝒔 𝟎
𝟏 𝒔 𝟏 𝟏
⇒ [ 𝟎 𝟏] [ [ ]] [ ]
𝒔𝟐 + 𝒔 + 𝟏 −𝟏 𝒔 + 𝟏 𝟎
𝒔 𝟏
⇒ [𝟎 𝒔 𝟐 + 𝒔 + 𝟏 𝒔 𝟐 + 𝒔 + 𝟏] [ 𝟏]
𝟏] [ −𝟏 𝒔+𝟏 𝟎
𝒔𝟐 + 𝒔 + 𝟏 𝒔𝟐 + 𝒔 + 𝟏
𝟏 𝒔+𝟏 𝟏 𝟏
⇒[ 𝟐 ] [ ] =
𝒔 +𝒔+𝟏 𝒔𝟐 + 𝒔 + 𝟏 𝟎 𝒔𝟐 + 𝒔 + 𝟏
Therefore, the transfer function of the system for the given state space model is
𝒀(𝒔) 𝟏
=[ 𝟐 ]
𝑼(𝒔) 𝒔 +𝒔+𝟏
Let us calculate the transfer function of the system represented in the state space model as,
𝒙̇ 𝟎 𝟏 𝒙𝟏 𝟏
𝑿̇ = [ 𝟏 ] = [ ] [𝒙 ] + [ ] [𝒖]
𝒙̇ 𝟐 −𝟐 −𝟐 𝟐 𝟏
𝒙𝟏
𝒀 = [𝟐 𝟑] [𝒙 ] + [𝟎][𝒖]
𝟐
Here,
𝟎 𝟏 𝟏
𝑨=[ ] , 𝑩 = [ ] , 𝑪 = [𝟐 𝟑], 𝑫 = [𝟎]
−𝟐 −𝟐 𝟏
−𝟏
𝒀(𝒔) 𝒔 𝟎 𝟎 𝟏 𝟏
= [𝟐 𝟑] [[ ]−[ ]] [ ]
𝑼(𝒔) 𝟎 𝒔 −𝟐 −𝟐 𝟏
𝒔 −𝟏 −𝟏 𝟏
⇒ [𝟐 𝟑] [ ] [ ]
𝟐 𝒔+𝟐 𝟏
𝟏 𝒔+𝟐 𝟏 𝟏
⇒ [𝟐 𝟑] [ 𝟐 [ ]] [ ]
𝒔 + 𝟐𝒔 + 𝟐 −𝟐 𝒔 𝟏
𝒔+𝟐 𝟏
⇒ [𝟐 𝒔𝟐 + 𝟐𝒔 + 𝟐 𝒔𝟐 + 𝟐𝒔 + 𝟐] [𝟏]
𝟑] [ −𝟐 𝒔 𝟏
𝒔𝟐 + 𝟐𝒔 + 𝟐 𝟐
𝒔 + 𝟐𝒔 + 𝟐
𝒔+𝟐 𝟏
+
⇒ [𝟐 𝒔𝟐 + 𝟐𝒔 + 𝟐 𝒔𝟐 + 𝟐𝒔 + 𝟐]
𝟑] [ −𝟐 𝒔
𝟐
+ 𝟐
𝒔 + 𝟐𝒔 + 𝟐 𝒔 + 𝟐𝒔 + 𝟐
𝒔+𝟑 𝒔−𝟐
⇒ [𝟐 𝟑] [ ]
𝒔𝟐 + 𝟐𝒔 + 𝟐 𝒔𝟐 + 𝟐𝒔 + 𝟐
𝟐𝒔 + 𝟔 𝟑𝒔 − 𝟔 𝟓𝒔
⇒[ 𝟐 ] = 𝟐
𝒔 + 𝟐𝒔 + 𝟐 𝒔𝟐 + 𝟐𝒔 + 𝟐 𝒔 + 𝟐𝒔 + 𝟐
Therefore, the transfer function of the system for the given state space model is
𝒀(𝒔) 𝟓𝒔
=[ 𝟐 ]
𝑼(𝒔) 𝒔 + 𝟐𝒔 + 𝟐
Let us now discuss controllability and observability of control system one by one.
1. 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.
• Find the determinant of matrix 𝑸𝑪 and if it is not equal to zero, then the control system
is controllable.
2. 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.
⋮
[𝑪𝑨𝒏−𝟏 ]
• Find the determinant of matrix 𝑸𝑶 and if it is not equal to zero, then the control
system is observable.
Example:
𝒙̇ −𝟏 −𝟏 𝒙𝟏 𝟏
𝑿̇ = [ 𝟏 ] = [ ] [𝒙 ] + [ ] [𝒖]
𝒙̇ 𝟐 𝟏 𝟎 𝟐 𝟎
Here,
−𝟏 −𝟏 𝟏
𝑨=[ ] , 𝑩 = [ ] , 𝑪 = [𝟎 𝟏], 𝑫 = [𝟎] 𝒂𝒏𝒅 𝒏 = 𝟐
𝟏 𝟎 𝟎
−𝟏 −𝟏 𝟏 −𝟏
𝑨𝑩 = [ ][ ] = [ ]
𝟏 𝟎 𝟎 𝟏
𝟏 −𝟏
⇒ 𝑸𝑪 = [ ] , (𝟏 × 𝟏) − (−𝟏 × 𝟎) = 𝟏
𝟎 𝟏
|𝑸𝑪 | = 𝟏 ≠ 𝟎
Since the determinant of matrix 𝑸𝑪 is not equal to zero, the given control system is
controllable.
For 𝒏 = 𝟐, the matrix 𝑸𝑶 will be -
𝑸𝑶 = [𝑪𝑻 𝑨𝑻 𝑪𝑻 ]
Here,
−𝟏 𝟏 𝟎
𝑨𝑻 = [ ] 𝒂𝒏𝒅 𝑪𝑻 = [ ]
−𝟏 𝟎 𝟏
𝟎 𝟏
⇒ 𝑸𝑶 = [ ] , (𝟎 × 𝟎) − (𝟏 × 𝟏) = −𝟏
𝟏 𝟎
⇒ |𝑸𝑶 | = −𝟏 ≠ 𝟎
Since, the determinant of matrix 𝑸𝑶 is not equal to zero, the given control system is
observable.
𝒙𝟏
𝒀 = [𝟏 𝟐 𝟏] [𝒙𝟐 ] + [𝟎][𝒖]
𝒙𝟑
Controllable? Observable?
𝑸𝑪 = [𝑩 𝑨𝑩 𝑨𝟐 𝑩], 𝝆(𝑸𝑪 ) = 𝒏
𝒙̇ = 𝑨𝒙 + 𝑩𝒖
𝒚 = 𝑪𝒙 + 𝑫𝒖 ,𝑫 = 𝟎
➢ If the controllability matrix has full rank, then the system is controllable.
𝟏 𝟎 𝟎
𝑸𝑪 = [𝑩 𝑨𝑩 𝑨 𝑩] = [𝟎 𝟐 𝟎 −𝟏] , ∆(𝑸𝑪 ) ≥ 𝟏, 𝝆(𝑸𝑪 ) = 𝒏(𝑨) = 𝟑
𝟎 −𝟏 𝟑
➢ If the observability matrix has full rank, then the system is observable.
𝑪 𝟏 𝟐 𝟏
𝑸𝑶 = [ 𝑪𝑨 ] = [−𝟏 −𝟐 −𝟏] , ∆(𝑸𝑶 ) ≤ 𝟎 , 𝝆(𝑸𝑶 ) ≠ 𝒏(𝑨) = 𝟏
𝑪𝑨𝟐 𝟏 𝟐 𝟏
% evaluate (sI-A)^(-1)
% get denominator
det(C*((s*eye(size(A,1))-A)^(-1))*B)
𝒙𝟏
𝒀 = [𝟏 𝟎 𝟏] [𝒙𝟐 ] + [𝟎][𝒖]
𝒙𝟑
Controllable? Observable?
𝑸𝑪 = [𝑩 𝑨𝑩 𝑨𝟐 𝑩], 𝝆(𝑸𝑪 ) = 𝒏
𝒙̇ = 𝑨𝒙 + 𝑩𝒖
𝒚 = 𝑪𝒙 + 𝑫𝒖 ,𝑫 = 𝟎
➢ If the controllability matrix has full rank, then the system is controllable.
𝟎 𝟏 𝟏 𝟎 𝟎 𝟎
𝑸𝑪 = [𝑩 𝑨𝑩 𝟐]
𝑨 𝑩 = [𝟏 𝟎 𝟎 𝟎 𝟐 𝟎] , ∆(𝑸𝑪 ) ≥ 𝟏, 𝝆(𝑸𝑪 ) = 𝒏(𝑨) = 𝟑
𝟎 𝟎 𝟐 𝟎 −𝟐 𝟎
➢ If the observability matrix has full rank, then the system is observable.
𝑪 𝟏 𝟎 𝟏
𝑸𝑶 = [ 𝑪𝑨 ] = [𝟎 𝟑 −𝟏] , ∆(𝑸𝑶 ) ≥ 𝟏 , 𝝆(𝑸𝑶 ) ≠ 𝒏(𝑨) = 𝟑
𝑪𝑨𝟐 𝟎 −𝟐 𝟒
% evaluate (sI-A)^(-1)
% get denominator
det(C*((s*eye(size(A,1))-A)^(-1))*B)
Homework:
𝒙𝟏
𝒀 = [𝟑 𝟎 𝟎] [𝒙𝟐 ] + [𝟎][𝒖]
𝒙𝟑
Controllable? Observable?
01/02/1444
1 REGISTRATION
05/02/1444
08/02/1444
2 • Laplace Transforms.
12/02/1444
15/02/1444
3 1 • Inverse Laplace Transforms. LAB # 1
19/02/1444 (7 points)
29/02/1444
5 • Response of the First Order System.
03/03/1444 LAB # 2
2
06/03/1444 (7 points)
6 • Response of the Second Order System
10/03/1444
13/03/1444
7 • Proportional Controller (P).
17/03/1444
20/03/1444
8 • Integral Controller (I).
24/03/1444
LAB # 3
27/03/1444 3
9 • Proportional–Integral Controller (PI). (7 points)
02/04/1444
05/04/1444 • Proportional–Integral–Derivative
10 Controller (PID).
09/04/1444
12/04/1444 FINAL EXAMINATION
11 LEC # 4
16/04/1444 (13 points)
Notes:
Laboratory
LAB # 1: 7 points, Continuous evaluation 1
LAB # 2: 7 points, Continuous evaluation 2
LAB # 3: 7 points, Continuous evaluation 3
LAB # 4: 13 points, Final Examination.
15
𝐹(𝑠) =
𝑠
2. 𝒇(𝒕) = 𝟓 + 𝟒𝒆−𝟐𝒕
5 4 5 4 1 5𝑠 + 10 + 4𝑠 9𝑠 + 10
𝐹(𝑠) = + = +( × )= = 2
𝑠 𝑠+2 𝑠 1 𝑠+2 𝑠 2 + 2𝑠 𝑠 + 2𝑠
3. 𝒇(𝒕) = 𝟓 + 𝟒𝒆+𝟐𝒕
5 4 5𝑠 − 10 + 4𝑠 9𝑠 − 10
𝐹(𝑠) = + = = 2
𝑠 𝑠−2 𝑠 2 − 2𝑠 𝑠 − 2𝑠
4. 𝒇(𝒕) = 𝟏𝒕 − 𝟐𝒆−𝟏𝒕
1 2 𝑠 + 1 − 2𝑠 2 −2𝑠 2 + 𝑠 + 1
𝐹(𝑠) = 2 − = 3 =
𝑠 𝑠+1 𝑠 + 1𝑠 2 𝑠 3 + 1𝑠 2
5. 𝒙(𝒕) = 𝟐𝟎 𝐬𝐢𝐧 𝟒𝒕
4 80 80
𝑋(𝑠) = 20 × = 2 = 2
𝑠2 +42 𝑠 + (4 × 4) 𝑠 + 16
6. 𝒙(𝒕) = 𝐬𝐢𝐧 𝟒𝒕
4 4 4
𝑋(𝑠) = 1 × = 2 = 2
𝑠2 +42 𝑠 + (4 × 4) 𝑠 + 16
7. 𝒚(𝒕) = 𝟐𝒕 + 𝟏𝐜𝐨𝐬 𝟑𝒕
2 𝑠 2 𝑠 2𝑠 2 + 18 + 𝑠 3 𝑠 3 + 2𝑠 2 + 18
𝑌(𝑠) = + = + = =
𝑠 2 𝑠 2 + 32 𝑠 2 𝑠 2 + 9 𝑠 4 + 9𝑠 2 𝑠 4 + 9𝑠 2
8. 𝒚(𝒕) = 𝟐𝒕 + 𝟐𝒄𝒐𝒔 𝟑𝒕
2 2 𝑠 2 2𝑠 2𝑠 2 + 18 + 2𝑠 3 2𝑠 3 + 2𝑠 2 + 18
𝑌(𝑠) = + ( × ) = + = =
𝑠2 1 𝑠 2 + 32 𝑠2 𝑠2 + 9 𝑠 4 + 9𝑠 2 𝑠 4 + 9𝑠 2
9. 𝒉(𝒕) = 𝟏𝟎𝟎 + 𝟏𝟒𝒕 + 𝟖 𝐜𝐨𝐬 𝟏𝒕
100 14 𝑠 100 14 8𝑠
𝐻(𝑠) = + 2 + (8 × 2 ) = + +
𝑠 𝑠 𝑠 + 12 𝑠 𝑠2 𝑠2 + 1
10. 𝒇(𝒕) = 𝒕𝟒
4! 4 × 3 × 2 × 1 24
𝐹(𝑠) = = = 5
𝑠 4+1 𝑠5 𝑠
11. 𝒇(𝒕) = 𝟐𝒕𝟒
4! 4 × 3 × 2 × 1 2 × 24 48
𝐹(𝑠) = = = = 5
𝑠 4+1 𝑠5 𝑠5 𝑠
8! 8 × 7 × 6 × 5 × 4 × 3 × 2 × 1 40320
𝐹(𝑠) = = =
𝑠 8+1 𝑠9 𝑠9
1
𝐹(𝑠) =
(𝑠 + 5)2
3×1 3
𝐹(𝑠) = =
(𝑠 − 5)2 (𝑠 − 5)2
4! 4×3×2×1 24
𝐹(𝑠) = = =
(𝑠 + 5)4+1 (𝑠 + 5)5 (𝑠 + 5)5
3! 3×2×1 6
𝐹(𝑠) = = =
(𝑠 + 5)3+1 (𝑠 + 5)4 (𝑠 + 5)4
6! 6×5×4×3×2×1 720
𝐹(𝑠) = = =
(𝑠 + 5)6+1 (𝑠 + 5)7 (𝑠 + 5)7
2 2
𝐹(𝑠) = =
(𝑠 + 6)2 + 22 (𝑠 + 6)2 + 4
2 8
𝐹(𝑠) = 4 × =
(𝑠 + 6)2 + 22 (𝑠 + 6)2 + 4
𝑠+3 𝑠+3
𝐹(𝑠) = =
(𝑠 + 3)2 + 52 (𝑠 + 3)2 + 25
𝑠−3 2(𝑠 − 3) 2𝑠 − 6
𝐹(𝑠) = 2 × = =
(𝑠 − 3)2 + 52 (𝑠 − 3)2 + 52 (𝑠 − 3)2 + 25
syms f t s
% STEP FUNCTION:
% f=10
% RAMP FUNCTION:
% f=10*t
% EXPONENTIAL FUNCTION:
% f=exp(-2*t);
% SINUSOIDAL FUNCTION:
% f=sin(10*t)
% f=cos(10*t)
% f(t)= t^3
% f(t)= t^3*exp(-5*t)
% f(t)= exp(-3*t)*sin(2*t)
% f(t)= exp(-3*t)*cos(2*t)
F=laplace(f,t,s)
pretty(F)
ℒ −1 [𝐹(𝑠)] = 𝑓(𝑡) = 𝑡 3
𝒔
4. 𝑭(𝒔) =
𝒔𝟐 +𝟑𝟔
𝟏𝟓 𝟓
1. 𝑭(𝒔) = 𝒔𝟐 +𝟐𝟓 = 𝟑 × 𝒔𝟐 +𝟓𝟐
syms f s t
% F(s)=1/(s+2)
% F(s)= (s-3)/((s-3)^2+81)
% F(s)=6/s^4
% F(s)= s/(s^2+36)
% F(s)= 15/(s^2+25)
% F(s)= 2/(2*s+4)
% F(s)= 15/((s+4)^2+9)
% F(s)= 96/((s+5)^5)
f(t)=ilaplace(F)
pretty(f)
(𝑠 + 3) (𝑠 + 1) (𝑠 + 3) (𝑠 + 3) −1 + 3 2
𝐴1 = |(𝑠 + 1) ( )| = |( )( )| =| | = = =2
(𝑠 + 1)(𝑠 + 2) 𝑠=−1 1 (𝑠 + 1)(𝑠 + 2) 𝑠=−1 (𝑠 + 2) 𝑠=−1 −1 + 2 1
(𝑠 + 3) (𝑠 + 2) (𝑠 + 3) (𝑠 + 3) −2 + 3 1
𝐴2 = |(𝑠 + 2) ( )| = |( )( )| =| | = = = −1
(𝑠 + 1)(𝑠 + 2) 𝑠=−2 1 (𝑠 + 1)(𝑠 + 2) 𝑠=−2 (𝑠 + 1) 𝑠=−2 −2 + 1 −1
𝟐 −𝟏
∴ 𝑭(𝒔) = +
(𝒔 + 𝟏) (𝒔 + 𝟐)
𝒔 𝒔 𝟏 𝟐 𝑨 𝑨
2. 𝑭(𝒔) = 𝒔𝟐 +𝟏𝟏𝒔+𝟑𝟎 = (𝒔+𝟓)(𝒔+𝟔) = (𝒔+𝟓) + (𝒔+𝟔) 𝒔𝟏 = −𝟓, 𝒔𝟐 = −𝟔
𝑠 (𝑠 + 5) 𝑠 𝑠 −5 −5
𝐴1 = |(𝑠 + 5) ( )| = |( )( )| =| | = = = −5
(𝑠 + 5)(𝑠 + 6) 𝑠=−5 1 (𝑠 + 5)(𝑠 + 6) 𝑠=−5 (𝑠 + 6) 𝑠=−5 −5 + 6 1
𝑠 (𝑠 + 6) 𝑠 𝑠 −6 −6
𝐴2 = |(𝑠 + 6) ( )| = |( )( )| =| | = = =6
(𝑠 + 5)(𝑠 + 6) 𝑠=−6 1 (𝑠 + 5)(𝑠 + 6) 𝑠=−6 (𝑠 + 5) 𝑠=−6 −6 + 5 −1
−𝟓 𝟔
∴ 𝑭(𝒔) = +
(𝒔 + 𝟓) (𝒔 + 𝟔)
𝟐𝟎𝟎 𝑨𝟏 𝟐𝑨
3. 𝑭(𝒔) = 𝒔(𝒔+𝟏𝟎) = + (𝒔+𝟏𝟎) 𝒔=𝟎
𝒔
𝟐𝟎 −𝟐𝟎
∴ 𝑭(𝒔) = +
𝒔 (𝒔 + 𝟏𝟎)
10 (𝑠 + 3) 10 10 10 10 5
𝐴1 = |(𝑠 + 3) ( )| = |( )( )| =| | = = =
(𝑠 + 3)(𝑠 + 9) 𝑠=−3 1 (𝑠 + 3)(𝑠 + 9) 𝑠=−3 (𝑠 + 9) 𝑠=−3 −3 + 9 6 3
10 (𝑠 + 9) 10 10 10 10 5
𝐴2 = |(𝑠 + 9) ( )| = |( )( )| =| | = = =
(𝑠 + 3)(𝑠 + 9) 𝑠=−9 1 (𝑠 + 3)(𝑠 + 6) 𝑠=−9 (𝑠 + 3) 𝑠=−9 −9 + 3 −6 −3
𝟓 𝟓
∴ 𝑭(𝒔) = 𝟑 − 𝟑
(𝒔 + 𝟑) (𝒔 + 𝟗)
5 5
ℒ −1 [𝐹(𝑠)] = 𝑓(𝑡) = 𝑒 −3𝑡 − 𝑒 −9𝑡
3 3
(𝒔+𝟑) (𝒔+𝟑) 𝟏 𝑨 𝟐𝑨
5. 𝑭(𝒔) = 𝒔𝟐 +𝟖𝒔+𝟏𝟐 = (𝒔+𝟐)(𝒔+𝟔) = (𝒔+𝟐) + (𝒔+𝟔)
(𝑠 + 3) (𝑠 + 2) (𝑠 + 3) (𝑠 + 3) −2 + 3 1
𝐴1 = |(𝑠 + 2) ( )| = |( )( )| =| | = =
(𝑠 + 2)(𝑠 + 6) 𝑠=−2 1 (𝑠 + 2)(𝑠 + 6) 𝑠=−2 (𝑠 + 6) 𝑠=−2 −2 + 6 4
(𝑠 + 3) (𝑠 + 6) (𝑠 + 3) (𝑠 + 3) −6 + 3 −3 3
𝐴2 = |(𝑠 + 6) ( )| = |( )( )| =| | = = =
(𝑠 + 2)(𝑠 + 6) 𝑠=−6 1 (𝑠 + 2)(𝑠 + 6) 𝑠=−6 (𝑠 + 2) 𝑠=−6 −6 + 2 −3 4
𝟏 𝟑
∴ 𝑭(𝒔) = 𝟒 + 𝟒
(𝒔 + 𝟐) (𝒔 + 𝟔)
1 3
ℒ −1 [𝐹(𝑠)] = 𝑓(𝑡) = 𝑒 −2𝑡 + 𝑒 −6𝑡
4 4
(𝒔+𝟏) (𝒔+𝟏) 𝟏 𝑨 𝑨
𝟐
6. 𝑭(𝒔) = 𝒔𝟐 +𝟏𝟐𝒔+𝟑𝟓 = (𝒔+𝟓)(𝒔+𝟕) = (𝒔+𝟓) + (𝒔+𝟕)
(𝑠 + 1) (𝑠 + 5) (𝑠 + 1) (𝑠 + 1) −5 + 1 −4
𝐴1 = |(𝑠 + 5) ( )| = |( )( )| =| | = = = −2
(𝑠 + 5)(𝑠 + 7) 𝑠=−5 1 (𝑠 + 5)(𝑠 + 7) 𝑠=−5 (𝑠 + 7) 𝑠=−5 −5 + 7 2
(𝑠 + 1) (𝑠 + 7) (𝑠 + 1) (𝑠 + 1) −7 + 1 −6
𝐴2 = |(𝑠 + 7) ( )| = |( )( )| =| | = = =3
(𝑠 + 5)(𝑠 + 7) 𝑠=−7 1 (𝑠 + 5)(𝑠 + 7) 𝑠=−7 (𝑠 + 5) 𝑠=−7 −7 + 5 −2
−𝟐 𝟑
∴ 𝑭(𝒔) = +
(𝒔 + 𝟓) (𝒔 + 𝟕)
12 (𝑠) 12 12 12 12
𝐴1 = |(𝑠) ( )| = |( ) ( )| =| | = =
𝑠(𝑠 + 1)(𝑠 + 4) 𝑠=0 1 𝑠(𝑠 + 1)(𝑠 + 4) 𝑠=0 (𝑠 + 1)(𝑠 + 4) 𝑠=0 (0 + 1)(0 + 4) 4
=3
12 (𝑠 + 1) 12 12 12 12
𝐴2 = |(𝑠 + 1) ( )| = |( )( )| =| | = = = −4
𝑠(𝑠 + 1)(𝑠 + 4) 𝑠=−1 1 𝑠(𝑠 + 1)(𝑠 + 4) 𝑠=−1 𝑠(𝑠 + 4) 𝑠=−1 −1 × (−1 + 4) −3
12 (𝑠 + 4) 12 12 12 12
𝐴3 = |(𝑠 + 4) ( )| = |( )( )| =| | = = =1
𝑠(𝑠 + 1)(𝑠 + 4) 𝑠=−4 1 𝑠(𝑠 + 1)(𝑠 + 4) 𝑠=−4 𝑠(𝑠 + 1) 𝑠=−4 −4 × (−4 + 1) 12
𝟑 𝟒 𝟏
∴ 𝑭(𝒔) = − +
𝒔 (𝒔 + 𝟏) (𝒔 + 𝟒)
𝒄(𝒕)𝒆−𝒕⁄𝝉 = 𝟎
𝟎, 𝑡<0
𝒓(𝒕) = {
𝟏, 𝑡≥0
𝑲 𝒓(𝒕) = 𝝉𝒄̀ (𝒕) + 𝒄(𝒕)
𝒄(𝒕)𝒆−𝒕⁄𝝉 = 𝟎
𝟎, 𝑡<0
𝒓(𝒕) = {
𝟐, 𝑡≥0
𝑲 𝒓(𝒕) = 𝝉𝒄̀ (𝒕) + 𝒄(𝒕)
𝒄(𝒕)𝒆−𝒕⁄𝝉 = 𝟎
𝟎, 𝑡<0
𝒓(𝒕) = {
𝟑, 𝑡≥0
(
𝑲 𝒓(𝒕) = 𝝉𝒄̀ 𝒕) + 𝒄(𝒕)
𝒄(𝒕)𝒆−𝒕⁄𝝉 = 𝟎
𝟎, 𝑡<0
𝒓(𝒕) = {
𝟒, 𝑡≥0
𝑲 𝒓(𝒕) = 𝝉𝒄̀ (𝒕) + 𝒄(𝒕)
𝒄(𝒕)𝒆−𝒕⁄𝝉 = 𝟎
𝟎, 𝑡<0
𝒓(𝒕) = {
𝟓, 𝑡≥0
𝑲 𝒓(𝒕) = 𝝉𝒄̀ (𝒕) + 𝒄(𝒕)
n=5
d=[2 1]
G_s=tf(n, d)
step (G_s)
figure(1);
grid on
grid minor
xlabel ('$$Time$$', 'interpreter', 'latex');
ylabel ('$$Vc(t)$$', 'interpreter', 'latex');
title ('$$Step Response$$', 'interpreter', 'latex');
𝒄(𝒕) = 𝒄′ (𝒕) = 𝟎
𝟎, 𝒕<0
𝒓(𝒕) = {
𝟏, 𝒕≥𝟎
𝒄(𝒕) = 𝒄′ (𝒕) = 𝟎
𝟎, 𝒕<0
𝒓(𝒕) = {
𝟐, 𝒕≥𝟎
𝒄(𝒕) = 𝒄′ (𝒕) = 𝟎
𝟎, 𝒕<0
𝒓(𝒕) = {
𝟑, 𝒕≥𝟎
𝒄(𝒕) = 𝒄′ (𝒕) = 𝟎
𝟎, 𝒕<0
𝒓(𝒕) = {
𝟒, 𝒕≥𝟎
𝒄(𝒕) = 𝒄′ (𝒕) = 𝟎
𝟎, 𝒕<0
𝒓(𝒕) = {
𝟓, 𝒕≥𝟎
n=18
d=[1 6 9]
G_s=tf(n, d)
step (G_s)
figure(1);
grid on
grid minor
xlabel ('$$Time$$', 'interpreter', 'latex');
ylabel ('$$Vc(t)$$', 'interpreter', 'latex');
title ('$$Step Response$$', 'interpreter', 'latex');
Library Browser:
• Step
• Sum
• Gain
• Integrator
• Derivative
• Add
• Transfer fcn
• Mux
• Scope
Library Browser:
• Step
• Sum
• Gain
• Integrator
• Add
• Transfer fcn
• Mux
• Scope
Library Browser:
• Step
• Sum
• Gain
• Integrator
• Add
• Transfer fcn
• Mux
• Scope
Library Browser:
• Step
• Sum
• Gain
• Integrator
• Derivative
• Add
• Transfer fcn
• Mux
• Scope