Automatic Control
Automatic Control
Automatic Control
(LECTURE NOTE)
By
1.0 INTRODUCTION
The functional diagrams shown below will help clarify the difference
between open and closed loop control systems.
Desired Output
Comparison Controller Process Out put
Response
Measurement
Desired Output
Comparison Controller Process Out put
Response
Measurement
An advantage of the closed loop control system is the fact that the
use of feedback makes the system response relatively insensitive to
external disturbances and internal variations in system parameters. It is
thus possible to use relatively inaccurate and inexpensive components to
obtain the accurate control of a given plant, whereas doing so is impossible
in the open-loop case. From the point of view of stability, the open-loop
control system is easier to build because system stability is not a major
problem. On the other hand, stability is a major problem in the closed-loop
control system, which may tend to overcorrect errors and thereby can
Switch Input
Temperature
Power Source
Oven
Current
Heating Element
Power Source
Fig. 3: Open loop Temperature – Oven Control System
1.2 Stability
Support
Stable unstable
𝑂𝑢𝑡𝑝𝑢𝑡
𝑇𝑟𝑎𝑛𝑠𝑓𝑒𝑟 𝐹𝑢𝑛𝑐𝑡𝑖𝑜𝑛 = =1
𝑖𝑛𝑝𝑢𝑡
2.0 Modeling
Before we attempt to design a control device for any system, we
have to test whether the system stabilizes or not. There are basically two
ways to do this.
i. Perturbed the actual system or subject it to actual input signal and
determine how it responds. Now use trial and error control until a
suitable control is achieved.
ii. Propose a physical and mathematical model of the system, using
suitable assumptions to an idealization by experimenting with the
real system. Determine the parameters of the model and then use
the model to select a proper controlling scheme. When an
appropriate model has been formed with respect to any type of input
signal, then the response to any type of input signal can easily be
determine by the computational method alone. However, at the
experimental stage, the response of the system selected signal gives
some good hints as to the stability of the system. The main type of
standard signals used are;
Model of physical systems are usually derived using physical laws such as;
In most cases the model takes the form of near ordinary differential
equations. One major example of a control system is the servo of the
angular position of heavy masses especially where absolute precision is
important.
Motor Platform
𝜃𝑖(𝑡)
Bearing
i. The platform and the gun is a simple rigid body to make one mass
ii. That the friction at the bearing is wet and hence a linear damper
iii. That the error [𝜀 = 𝜃𝑖(𝑡) − 𝜃𝑜(𝑡) ]is small and therefore linear
analysis applies.
iv. The presence of external disturbance, e.g. wind and some internal
losses are negligible.
v. That the shaft connecting the motor to the platform is rigid
𝜃𝑜
𝜃𝑖(𝑡) Motor C
Damper
C(s)+ 𝜃𝑜(𝑡)
𝜃𝑖(𝑡)
Amplifier Motor
R(s)-
Damper
Measurement
Where; 𝜀 = (𝜃𝑖 − 𝜃𝑜 )
Hence; 𝑇𝑚 = 𝑘(𝜃𝑖 − 𝜃𝑜 )
𝑑𝜃𝑜
Torque from the damper 𝑇𝐷 = 𝐶
𝑑𝑡
𝑑 2 𝜃𝑜
Σ𝑚
̅̅̅ = 𝐼
𝑑𝑡 2
𝑑𝜃𝑜 𝑑 2 𝜃𝑜
Hence; 𝑘(𝜃𝑖 − 𝜃𝑜 ) − 𝐶 =𝐼
𝑑𝑡 𝑑𝑡 2
Rearranging we have;
𝑑 2 𝜃𝑜 𝑑𝜃𝑜
𝐼 +𝐶 + 𝑘𝜃𝑜 = 𝑘𝜃𝑖
𝑑𝑡 2 𝑑𝑡
There are basically four main approaches that can be used to analyze
control system;
There are three main criteria that are used to assess the performance
of a control system, they are;
i. Adequate stability
ii. Fast response
iii. Adequate steady state accuracy
use of which the system model is developed. When the feed-back loop is
closed, the system might be
i. Unstable
ii. Have too large steady state error
iii. Have too small damping error
iv. Have too small natural frequency
iii. The tolerance band may be specified and setting time 𝑇𝑠 given. It is
usually tabulated as shown below;
Tolerance (%) 1 2 3 4 5 6 10 15 20
Optimum (b) 0.82 0.80 0.75 0.72 0.69 0.63 0.57 0.52 0.46
Setting Time (𝑤𝑛 𝑡𝑚 ) 4.0 3.5 3.2 3.0 2.8 2.6 2.3 2.0 1.8
𝑘 𝑘 𝐶
We know that 𝑤𝑛 = √ which implies that 𝑤𝑛 2 = then = 2𝑏𝑤𝑛
𝐼 𝐼 𝐼
Hence;
𝜃𝑜̈ + 2𝑏𝑤𝑛 𝜃𝑜̇ + 𝑤𝑛 2 𝜃𝑜 = 𝑤𝑛 2 𝜃𝑖 ……………….. (ii)
The solution to the differential equation (ii) is divided into two parts
i.e. the homogenous solution and the particular integral solution.
𝜃𝑜(𝑡) = 𝜃𝑜ℎ(𝑡) + 𝜃𝑜𝑝(𝑡) ……………………………….. (iii)
To get the homogenous solution we set the right-hand side to zero.
Therefore;
𝜃𝑜̈ + 2𝑏𝑤𝑛 𝜃𝑜̇ + 𝑤𝑛 2 𝜃𝑜 = 0 ………………………… (iv)
Consider a solution of the form;
𝜃𝑜 = 𝐴𝑒 𝑚𝑡 …………………………………………….. (v)
Substituting equation (v) into (iv) to get
(𝑚2 + 2𝑏𝑤𝑛 𝑚 + 𝑤𝑛 2 )𝐴𝑒 𝑚𝑡 = 0
At 𝑡 = 0 and A = 1 the solution becomes;
𝑚1,2 = −𝑏𝑤𝑛 ± 𝑤𝑛 √𝑏 2 − 1 ……………………… (vi)
In this equation there are three cases to consider;
When
i. 𝑏>1
ii. 𝑏=1
iii. 𝑏<1
𝜃𝑖 𝑏<1
𝑏=1
𝑏>1
From the graph above we see that the most interesting case is
when 𝑏 < 1, we now proceed to rewrite the equation for the case, and the
equation we are going to rewrite is (xvi).
Now boundary conditions;
𝜃𝑜 = 0 When 𝑡 = 0, i.e., 𝜃𝑜̇ (0) = 0 ……… (xvii)
From equation (xvi) differentiating we have;
𝜃𝑜̇ (𝑡) = 𝑒 (−𝑏𝑤𝑛 )𝑡 [−𝐴 𝑤𝑑 sin 𝑤𝑑 𝑡 + 𝐵 𝑤𝑑 cos 𝑤𝑑 𝑡] − 𝑏𝑤𝑛 𝑒 (−𝑏𝑤𝑛 )𝑡 [Acos 𝑤𝑑 𝑡 +
𝐵 sin 𝑤𝑑 𝑡] ..…. (xviii)
𝑡 −𝑏𝑤𝑛
𝜃𝑜(𝑡) = 𝜃𝑖 [ 1 − 𝑒 (−𝑏𝑤𝑛 ) [cos 𝑤𝑑 𝑡 + sin 𝑤𝑑 𝑡] ……….. (xxi)
𝑤𝑑
𝑤𝑑 𝑏𝑤𝑛
cCos𝛼 =
𝑤𝑑
Substituting in xxi;
𝑡
𝜃𝑜(𝑡) = 𝜃𝑖 [ 1 − 𝑒 (−𝑏𝑤𝑛 ) [𝑐Sin𝛼. cos 𝑤𝑑 𝑡 − 𝑐Cos𝛼 sin 𝑤𝑑 𝑡]
𝑡
𝜃𝑜(𝑡) = 𝜃𝑖 [ 1 − 𝑒 (−𝑏𝑤𝑛 ) . 𝑐[Sin 𝑤𝑑 𝑡 + 𝛼] ……… (xxii)
𝑏2 𝑤𝑛 2
𝑏2 𝑤𝑛 2 𝑏2 𝑐2 = 𝑎𝑛𝑑
Where; 𝑐 = √1 + = √1 + ( 𝑤𝑑 2 )
𝑤𝑑 2 1−𝑏2 2
𝑤𝑑 = 𝑤𝑛 √ 1 − 𝑏
1
This implies that 𝑐 = ………….. (xxiii)
√1−𝑏2
𝑤𝑑 √1−𝑏2
𝑡𝑎𝑛𝛼 = = …………. (xxiv)
𝑏𝑤𝑛 𝑏
𝑡
𝑒 (−𝑏𝑤𝑛)
Hence; 𝜃𝑜 = 𝜃𝑖 [ 1 − [Sin 𝑤𝑑 𝑡 + 𝛼] ………… (xxv)
√1−𝑏2
The first term in the bracket i.e., 1 is constant in all time, but the second
term vanishes for larger values of time. It may be noted that actually the first
term is the particular solution while the second term is the homogenous.
Because of the time behavior of the second term the first term is often called
steady state solution and the second term is called transient solution. From
𝜃𝑜 (𝑡)
the graph of , we note that b ≥ 1 the response is slow but it has no
𝑤𝑛 (𝑡)
overshoot on the other hand when b < 1 the response is very fast but there is
an overshoot of the steady state value. Thus, a compromise has always got to
be struck between fast response and stability of the response. From this
discussion a number of criteria can be deduced for the desirable performance
of a dynamic system.
𝜃𝑖 = 1 𝑟𝑎𝑑
Now, 𝜃𝑜̈ + 2𝑏𝑤𝑛 𝜃𝑜̇ + 𝑤𝑛 2 𝜃𝑜 = 𝑤𝑛 2 𝜃𝑖
𝑘 20
We know that 𝑤𝑛 = √ = √ = 14.14 𝑟𝑎𝑑/𝑠
𝐼 0.1
𝑤𝑛 2 = 14.142 ≅ 200
2𝑏𝑤𝑛 = 2 × 0.5 × 14.14
= 14.14
Putting these values in 𝜃𝑜̈ + 2𝑏𝑤𝑛 𝜃𝑜̇ + 𝑤𝑛 2 𝜃𝑜 = 𝑤𝑛 2 𝜃𝑖 we have;
𝜃𝑜̈ + 14.14𝜃𝑜̇ + 200𝜃𝑜 = 200 × 1
For transient solution;
𝑒 −𝑏𝑤𝑛𝑡
𝜃𝑜 (𝑡) = sin (𝑤𝑑 𝑡 + 𝛼)
√1−𝑏2
√1 − 𝑏 2 = √1 − (0.5)2 = 0.866
∴ 𝜃𝑜 (𝑡) = 1.15[𝑒 −7.07𝑡 𝑠𝑖𝑛(𝑤𝑑 𝑡 + 𝛼)]
But 𝑤𝑑 = 𝑤𝑛 √1 − 𝑏 2 = 14.14 × 0.866 = 12.25 𝑟𝑎𝑑/𝑠
𝑤𝑑 1−𝑏2
𝑇𝑎𝑛𝛼 = =√ = 1.732
𝑏𝑤𝑛 𝑏
𝑘
But 𝑤𝑛 = √ 𝑘 = 𝑤𝑛 2 𝐼 = 1. 42 × 0.1 = 0.196𝑁/𝑚
𝐼
𝐶
For C = 2𝑏𝑤𝑛 𝐶 = 2𝑏𝑤𝑛 𝐼 = 2 × 0.69 × 1.4 × 0.1
𝐼
𝐶 = 0.1932
6.0 Complex Domain Analysis
Example 5: Let 𝑓(𝑡) = 𝑒 𝑎𝑡 . Compute 𝐹(𝑠) = ℒ(𝑓(𝑡; 𝑠)) directly. Give the
region in the complex s-plane.
∞ ∞ 1 ∞
𝐹(𝑠) = ℒ(𝑓(𝑡; 𝑠)) = ∫0 𝑒 𝑎𝑡 . 𝑒 −𝑠𝑡 𝑑𝑡 = ∫0 𝑒 (𝑎−𝑠)𝑡 𝑑𝑡 = 𝑒 (𝑎−𝑠)𝑡 ⌉
𝑎−𝑠 0
1 1 1
(0 − 𝑒 0 ) = (−1) =
𝑎−𝑠 𝑎−𝑠 𝑠−𝑎
50
𝐵= | ∴ 𝐵 = 16.67
𝑠(𝑠+0.5) 𝑠=−2
50
𝐶= | ∴ 𝐶 = −66.67
𝑠(𝑠+2) 𝑠=−0.5
50 50 16.67 −66.67
𝐸𝑜 (𝑡) = = + +
𝑠(𝑠+2)(𝑠+0.5) 𝑠 𝑠+2 𝑠+0.5
6+𝑆 2 +3𝑆 𝐾1 2 𝐾 3 𝐾
Step III: = + (𝑆+1) + (𝑆+2)
𝑆(𝑆+1)(𝑆+2) 𝑆
6+𝑆 2 +3𝑆 6
And 𝐾1 = (𝑆+1)(𝑆+2)⌋ = =3
𝑆=0 2
3 4 2
Then: 𝑥(𝑠) = − (𝑆+1) + (𝑆+2)
5
I (s ) =
10
(s + 1)(s + 2)
10
Solution: 𝐼(𝑠) = (𝑆+1)(𝑆+2)
10 1 𝐾 2 𝐾
(𝑆+1)(𝑆+2)
= (𝑆+1) + (𝑆+2)
10 1 𝐾
(𝑆+1)(𝑆+2)
= (𝑆+1)
10 10
𝐾1 = (𝑆+2)| = (−1+2)
= 10
𝑆 = −1
Now for 𝐾2
10 2 𝐾
(𝑆+1)(𝑆+2)
= (𝑆+2)
10 10
𝐾2 = (𝑆+1)| = (−2+1)
= −10
𝑆 = −2
10 10
Therefore: 𝐼(𝑠) = (𝑆+1) − (𝑆+2)
Where R(s) is the Laplace transform of the input to the system and
C(s) is the Laplace transform of the output of the system, where we assume
that all initial conditions involved are zero. It follows that the output C(s) can
be written as the product of G(s) and R(s):
𝐶(𝑠) = 𝑅(𝑠) 𝐺(𝑠).
leading away from the block represents the output. Such arrows are referred
to as signals.
R(s) C(s)
Input Transfer Function G(s) Output
A C = A±B
G1
±
B G2
Summing Point
Take Off Point
𝐶(𝑠) 𝐺(𝑠)
= …6
𝑅(𝑠) 1+𝐻(𝑠).𝐺(𝑠)
𝐶(𝑠) 𝐺(𝑠)
= Closed loop transfer function =
𝑅(𝑠) 1+𝐻(𝑠).𝐺(𝑠)
Put the values of C(s) from equation (1) into equation (6)
𝐸(𝑠) 𝐺(𝑠) 𝐺(𝑠)
=
𝑅(𝑠) 1+𝐻(𝑠).𝐺(𝑠)
𝐸(𝑠) 1
= …8
𝑅(𝑠) 1+𝐻(𝑠).𝐺(𝑠)
𝐸(𝑠) 1
= Error Ratio =
𝑅(𝑠) 1+𝐻(𝑠).𝐺(𝑠)
Put the values of C(s) from equation (1) into equation (2)
𝐵(𝑠) = 𝐻(𝑠). 𝐸(𝑠). 𝐺(𝑠) … 10
Put the values of E(s) from equation (3) into equation (10)
𝐵(𝑠) = 𝐻(𝑠). 𝐺(𝑠)[𝑅(𝑠) − 𝐵(𝑠)]
𝐵(𝑠) 𝐻(𝑠).𝐺(𝑠)
= … 11
𝑅(𝑠) 1+𝐻(𝑠).𝐺(𝑠)
𝐵(𝑠) 𝐻(𝑠).𝐺(𝑠)
= Primary feedback ratio =
𝑅(𝑠) 1+𝐻(𝑠).𝐺(𝑠)
𝐶(𝑠) 𝐺(𝑠)
= For positive Feedback
𝑅(𝑠) 1−𝐻(𝑠).𝐺(𝑠)
R2(s)
R1(s) C(s)
E(s) G1(s) ± G2(s)
±
B(s)
H(s)
Thus,
𝐶1 (𝑠)
= 𝐺1 (𝑠)
𝑅(𝑠)
𝐶(𝑠)
= 𝐺2 (𝑠)
𝐶1 (𝑠)
𝐶1 (𝑠) 𝐶(𝑠) 𝐶(𝑠)
∴ . (𝑠)
= = 𝐺1 (𝑠). 𝐺2 (𝑠)
𝑅(𝑠) 𝐶1 𝑅(𝑠)
G1(s)
G3(s)
C(s)
R(s) C(s)
R(s) G(s) G(s)
≡
G(s)
≡ 1
𝐺(𝑠)
+ + X + X+Z + X + Z-Y
X X-Y X–Y+Z
+ + −
−
Y Z ≡ Z Y
Example 11: Derive the transfer function of the given figure below using block
reduction technique.
R(s) + + + C(s)
G1(s) G2(s)
− − −
H1(s) H2(s)
H3(s)
Solution
Step I: There are two internal closed loops in the system, reduce them by using
equation 6.
𝐺1 𝐺2
R(s) + C(s)
1 + 𝐻1 . 𝐺1 1 + 𝐻2 . 𝐺2
−
H3(s)
𝐺1 𝐺2
R(s) + C(s)
(1 + 𝐻1 . 𝐺1 )(1 + 𝐻2 . 𝐺2 )
−
H3(s)
Step III: In step II the block is a closed loop again, use equation 6.
R(s) 𝐺1 𝐺2 C(s)
(1 + 𝐻1 . 𝐺1 + 𝐻2 . 𝐺2 + 𝐺1 𝐺2 𝐻3 . +𝐺1 𝐺2 𝐻1 𝐻2 )
𝐶(𝑠) 𝐺1 𝐺2
∴ = (1+𝐻
𝑅(𝑠) 1 .𝐺1 +𝐻2 .𝐺2 +𝐺1 𝐺2 𝐻3 .+𝐺1 𝐺2 𝐻1 𝐻2 )
Example 12: Find the overall transfer function for the system shown in the figure
1
below. The pick of point is beyond the block .
𝑠𝐶2
Solution:
𝑠𝐶2
− 1
Vi(s) + 1 + 1 1 Vo(s)
𝑅1 + 𝑠𝐶1 𝑅2 𝑠𝐶2
− -
Step I:
𝑠𝐶2
− 1
Vi(s) + 1 + 1 Vo(s)
𝑅1 + 𝑠𝐶1 𝑠𝐶2 𝑅2
− -
Step II:
𝑠𝐶2
− 1
Vi(s) + 1 1 Vo(s)
𝑅1 + 𝑠𝐶1 1 + 𝑠𝐶2 𝑅2
−
Step III:
𝑠𝐶2
− 1
Vi(s) + 1 1 Vo(s)
𝑅1 + 𝑠𝐶1 1 + 𝑠𝐶2 𝑅2
−
1 + 𝑠𝐶2 𝑅2
Step IV:
𝑠𝐶2
− 1
Vi(s) + 1 Vo(s)
𝑅1 + 𝑠𝐶1 + 𝑅2 𝐶2 𝐶1 𝑠 2
−
1 + 𝑠𝐶2 𝑅2
Step V:
Vi(s) + 1 1 Vo(s)
𝑅1 𝑠𝐶1 + 𝑅2 𝐶2 𝐶1 𝑠 2 + 𝑠𝐶2
−
1 + 𝑠𝐶2 𝑅2
Step VI:
Vi(s) + 1 Vo(s)
𝑅1 𝑠𝐶1 + 𝑅1 𝑅2 𝐶2 𝐶1 𝑠 2 + 𝑅1 𝑠𝐶2
−
1 + 𝑠𝐶2 𝑅2
Step VII:
Vi(s) 1 Vo(s)
2
1 + 𝑅1 𝑅2 𝐶2 𝐶1 𝑠 + (𝑅1 𝐶1 + 𝑅1 𝐶2 + 𝑅2 𝐶2 )𝑠
𝑉𝑜 (𝑠) 1
= 2
𝑉𝑖 (𝑠) 1 + 𝑅1 𝑅2 𝐶2 𝐶1 𝑠 + (𝑅1 𝐶1 + 𝑅1 𝐶2 + 𝑅2 𝐶2 )𝑠
𝐶(𝑠)
Example 13: Determine the ratio the system shown in the figure below
𝑅(𝑠)
𝐻2
-
R(s) + + + C(s)
𝐺1 𝐺2 𝐺3
− −
𝐻1
Example 14: The figure below shows a block diagram of control system, determine
i. Opened loop transfer function of the system
ii. Closed loop transfer function of the system
iii. Characteristic equation of the system
iv. Reduce the system block diagram to single closed loop system
R(s) E(s) C(s)
50/S 30/S
20
20
Thus;
R(s) E(s) T(s) C(s)
50/S 30/S+600
(50)(30) 1500
Now; 𝐷𝑇𝐺 = =
𝑆(𝑆+600) 𝑆(𝑆+600)
50 30
Open loop transfer function = Loop Gain = × ×5
𝑆 𝑆+600
7500 7500
= =
𝑆(𝑆+600) 𝑆 2 +600𝑆
𝐷𝑇𝐺
Closed loop transfer function=
1+𝐿𝑜𝑜𝑝 𝐺𝑎𝑖𝑛
1500
𝑆(𝑆+600) 1500
= 7500 =
1+ 𝑆 2 +600𝑆+7500
𝑆(𝑆+600)
𝑆 2 + 600𝑆 + 7500 = 0
The reduced system diagram is as follows:
R(s) 1500 C(s)
𝑆 2 + 600𝑆 + 7500
Example15: Find the transfer function of the system shown in the figure below
1
3
𝑠
-1
-6 1
𝑠
-1
3 1
𝑠
-1
Step I:
1
3
𝑠+1
1
R(s) -6 C(s)
𝑠+2
3 1
𝑠+3
Step II 3
𝑠+1
−6
R(s) 𝑠+2 C(s)
3
𝑠+3
Step III
6
R(s) 𝑠3 + 6𝑠 2 + 11𝑠 + 6 C(s)
𝐶(𝑠) 6
∴ =
𝑅(𝑠) 𝑠 3 +6𝑠 2 +11𝑠+6
Row
0 𝑺𝒏 𝒂𝟎 𝒂𝟐 𝒂𝟒 𝒂𝟔
1 𝑺𝒏−𝟏 𝒂𝟏 𝒂𝟑 𝒂𝟓 𝒂𝟕
2 𝑺𝒏−𝟐 𝒃𝟏 𝒃𝟐 𝒃𝟑 𝒃𝟒
3 𝑺𝒏−𝟑 𝒄𝟏 𝒄𝟐 𝒄𝟑 𝒄𝟒
. . . . . .
. . . . . .
n-1 𝑺𝟏 𝒚𝟏 𝒚
n 𝑺𝟎 𝒛𝟏
Where the coefficients 𝑎0 , 𝑎1 … … . 𝑎𝑛 are all of the same sign and none is
zero. For this polynomial the Routh Array is where each element in the
third row onwards is obtained from the elements in the previous rows in
the following manner.
𝒂𝟏 𝒂𝟐 −𝒂𝟎 𝒂𝟑 𝒂𝟏 𝒂𝟒 −𝒂𝟎 𝒂𝟓
Here: 𝒃𝟏 = 𝒃𝟐 =
𝒂𝟏 𝒂𝟏
𝒃𝟏 𝒂𝟑 −𝒂𝟏 𝒃𝟐 𝒃𝟏 𝒂𝟓 −𝒂𝟏 𝒃𝟑
Similarly: 𝒄𝟏 = 𝒄𝟐 =
𝒃𝟏 𝒃𝟏
S 4 + 2S 3 + 6 S 2 + 4S + 1 = 0
Row
0 S4 1 6 1
1 S3 2 4
2 S2 4 1
3 S1 3.5
4 S0 1
𝒂𝟏 𝒂𝟐 −𝒂𝟎 𝒂𝟑 𝒂𝟏 𝒂𝟒 −𝒂𝟎 𝒂𝟓
Here: 𝒃𝟏 = 𝒃𝟐 =
𝒂𝟏 𝒂𝟏
(𝟐×𝟔)−(𝟏×𝟒) 𝟐×𝟏−𝟏×𝟎
𝒃𝟏 = =𝟒 𝒃𝟐 = =𝟏
𝟐 𝟐
𝒃𝟏 𝒂𝟑 −𝒂𝟏 𝒃𝟐
Similarly: 𝒄𝟏 =
𝒃𝟏
𝟒×𝟒−𝟐×𝟏
𝒄𝟏 = = 𝟑. 𝟓
𝟒
𝟑.𝟓×𝟏−𝟒×𝟎
Also; 𝒅𝟏 = =𝟏
𝟑.𝟓
Since all the coefficients in the column are of the same sign (positive), the
given equation has no roots with positive real parts. Hence, the system is
stable.
Solution:
Row
0 S5 1 6 8
1 S4 7 42 56
1 6 8
2 S3 0 0 0
3 S2
4 S1
5 S0
𝒂𝟏 𝒂𝟐 −𝒂𝟎 𝒂𝟑 𝒂𝟏 𝒂𝟒 −𝒂𝟎 𝒂𝟓
𝒃𝟏 = 𝒃𝟐 =
𝒂𝟏 𝒂𝟏
𝒃𝟏 𝒂𝟑 −𝒂𝟏 𝒃𝟐 𝒃𝟏 𝒂𝟓 −𝒂𝟏 𝒃𝟑
Similarly: 𝒄𝟏 = 𝒄𝟐 =
𝒃𝟏 𝒃𝟏
We are faced with the problem of zeros in third row, thus from the new
polynomials using the entries in the row above zeros. The polynomials will
start with power of s in that row, and continue by skipping every other
power of s, i.e.
𝐹(𝑠) = 𝑺𝟒 + 𝟔𝑺𝟐 + 𝟖
𝑑𝐹(𝑆)
Then; = 4𝑺𝟑 + 𝟏𝟐𝑺 + 𝟎
𝑑𝑠
Hence, the row with all zeros in the Routh table is replaced with these
coefficients. Therefore:
2 𝑺𝟑 𝟒 𝟏𝟐 𝟎
1 3 0
3 𝑺𝟐 3 8 0
4 𝑺𝟏 𝟏 0 0
𝟑
5 𝑺𝟎 8 0 0
Since there are no sign changes in the first column, the system is said to be
stable.
Example 14: Check the stability of the system having the following
characteristic equation
S 4 + 2 S 3 + 3 S 2 + 4S + 5 = 0
Solution:
Row
0 S4 1 3 5
1 S3 2 4
2 S2 1 5
3 S1 -6
4 S0 5
From the above result there are two number of sign changes, thus the
system is unstable.