Scribd
Upload
25 views4 pages

Lab No. 4-System Stabilty Using Routh Table, Root Locus

The document discusses the Routh-Hurwitz criterion and root locus analysis in MATLAB. It provides an example of using the Routh-Hurwitz criterion to check the stability of a transfer function with a variable gain K. The poles are plotted for different K values, showing the system is stable for K between 6 and 20. It also provides an example of using root locus analysis to plot the root locus and find the value of K where the damping ratio is 0.707.

Uploaded by

Ashno Khan
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
25 views4 pages

Lab No. 4-System Stabilty Using Routh Table, Root Locus

The document discusses the Routh-Hurwitz criterion and root locus analysis in MATLAB. It provides an example of using the Routh-Hurwitz criterion to check the stability of a transfer function with a variable gain K. The poles are plotted for different K values, showing the system is stable for K between 6 and 20. It also provides an example of using root locus analysis to plot the root locus and find the value of K where the damping ratio is 0.707.

Uploaded by

Ashno Khan
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 4

Lab No.

04
Section:A Routh-Hurwitz Criterion in Matlab
General Description:
This is a very useful method which gives us the information about the stability of a system
without the need to solve the system and find the closed loop system poles. Using this method
we can tell how many closed loop system-poles are located in right half plane, in the left half
plane and on the jw axis.
The Routh Hurwitz criterion declares that the number of roots of polynomial that are in right half
plane is equal to the number of sign changes in first column.
Problem:
If Transfer function of a system is given as:

i) K=[1:1:20]
ii) For what value of K system will be stable.
Source Code:

k=[6:0.5:20]; % Declaring the value of k

n=length(k); %Calculating the length of k

for i=1:n % using the for loop which executes n times

q=[1 10 10 5*k(i)]; % the denominator of a given T.F

p(:,i)=roots(q); %Calculating the poles for different Value of K

end %Ending of loop

plot(real(p),imag(p),'x') %Plotting the Graph


grid
System Response:

Fig 4.1: Location of Poles of the Closed Loop System

Poles are located in the left half plane, so the system is stable.
Section:B Root Locus in Matlab
General Description:
In control theory and stability theory, root locus analysis is a graphical method for examining
how the roots of a system change with variation of a certain system parameter, commonly a gain
within a feedback system.
The root locus technique can be used to analyze and design the effect of loop gain upon the
system's transient response and stability. This technique provides very useful information about
system parameters. As one says if you can measure it , you can control it so this technique is
very useful in this regard.
Root Locus in Matlab:
Syntax
rlocus
rlocus(sys)
rlocus(sys1,sys2,...)
Description:
rlocus computes the Evans root locus of a SISO open-loop model. The root locus gives the
closed-loop pole trajectories as a function of the feedback gain (assuming negative feedback).
Root loci are used to study the effects of varying feedback gains on closed-loop pole locations.
In turn, these locations provide indirect information on the time and frequency responses.
rlocus(sys) calculates and plots the root locus of the open-loop SISO model sys.
Problem 01:

Plot the root locus


ii) Find the Value where damping ratio is 0.707.

Source Code:
%Problem #01
k=[1:10:1000] % The value of K where damping ratio is 0.7
ng=[1] % declaring the feed-forward neumenator.

dg=[1 8 10 1] % declaring the feed-forward denuminator.


nh=[1] % declaring the feed-back neumenator.

dh=[1] % declaring the feed-back denuminator.

[numc,denc]=feedback(ng,dg,nh,dh,-1) %Solving for unity feed back gain


sys=tf(numc,denc)
%finding the loci
rlocus(sys,k)

You might also like