0% found this document useful (0 votes)
14 views3 pages

Lab # 10

The document discusses using the Routh stability criteria to analyze system stability in MATLAB. It provides the theory behind the Routh stability criteria for determining if a system is stable or unstable based on the signs of the coefficients of the characteristic equation. The lab objective is to construct a MATLAB program to generate the Routh matrix from input coefficients and check stability. The program operation and an example of running the program and checking stability are included.

Uploaded by

saadshahab622
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
0% found this document useful (0 votes)
14 views3 pages

Lab # 10

The document discusses using the Routh stability criteria to analyze system stability in MATLAB. It provides the theory behind the Routh stability criteria for determining if a system is stable or unstable based on the signs of the coefficients of the characteristic equation. The lab objective is to construct a MATLAB program to generate the Routh matrix from input coefficients and check stability. The program operation and an example of running the program and checking stability are included.

Uploaded by

saadshahab622
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/ 3

Lab # 10

Object: To Construct system stability by using Routh stability criteria in MATLAB.


Theory:
The theory of network synthesis states that any pole of the system lies on the right hand side
of the origin of the s plane, it makes the system unstable. A necessary and sufficient condition
for the system to be stable is that the real parts of the roots of the characteristic equation have
negative real parts. This ensures that the impulse response will decay exponentially with time.
On the basis of this condition, A. Hurwitz and E. J. Routh started investigating the necessary
and sufficient conditions of stability of a system.

The Routh Criterion: All the roots of the characteristic equation have negative real parts if
and only if the elements of the first column of the Routh table have the same sign. Otherwise,
the number of roots with positive real parts is equal to the number of changes of sign.
The Routh criterion is a method for determining continuous system stability, for systems with
an nth-order characteristic equation of the form:
𝑎𝑛 𝑠 𝑛 + 𝑎𝑛−1 𝑠 𝑛−1 + ⋯ + 𝑎1 𝑠 + 𝑎0 = 0
The criterion is applied using a Routh table defined as follows:
𝑠 𝑛 𝑎𝑛 𝑎𝑛−2 𝑎𝑛−4 ⋯
𝑠 𝑛−1 𝑎𝑛−1 𝑎𝑛−3 𝑎𝑛−5 ⋯
⋮ || 𝑏1 𝑏2 𝑏3 ⋯||
⋮ 𝑐1 𝑐2 𝑐3 ⋯
⋮ ⋯ ⋯ ⋯ ⋯
Procedure:
 Enter the program in the editor window.
 Execute the program.
 Enter the coefficients of characteristic equation in the command window.
 Finally save the program and run it.
 The Routh matrix is obtained in the command window.

Lab Practice:
%% firstly it is required to get first two row of routh matrix
e=input('Enter the coefficients of characteristic equation: ');
l=length(e);
m=mod(l,2);
if m==0
a=rand(1,(l/2));
b=rand(1,(l/2));
for i=1:(l/2)

65
a(i)=e((2*i)-1);
b(i)=e(2*i);
end
else e1=[e 0];
a=rand(1,((l+1)/2));
b=[rand(1,((l-1)/2)),0];
for i=1:((l+1)/2)
a(i)=e1((2*i)-1);
b(i)=e1(2*i);
end
end
%% now we generate the remaining rows of Routh matrix
l1=length(a);
c=zeros(l,l1);
c(1,:)=a;
c(2,:)=b;
for m=3:l
for n=1:l1-1
c(m,n)=-(1/c(m-1,1))*det([c((m-2),1) c((m-2),(n+1));c((m-1),1) c((m-
1),(n+1))]);
end
end
disp('The Routh Matrix is obtain as:')
disp(c)
%% now we check the stablity of system
if c(:,1)>0
disp('System is Stable')
else
disp('System is Unstable');
end

66
Output:
1. Enter the coefficients of characteristic equation: [1 1 3 1 6]

The Routh Matrix is obtain as:


1 3 6
1 1 0
2 6 0
−2 0 0
6 0 0
System is Unstable

2. Enter the coefficients of characteristic equation: [1 6 12 8]


The Routh Matrix is obtain as:
1 12
6 8
10.66 0
8 0
System is Stable

67

You might also like