routh_t

Arguments

p

a real polynomial

h

a real SISO transfer system

k

a real polynomial or a scalar

kp

a scalar: proportional controller constant

normalized

a boolean (%t (default value) or %f)

r

a matrix or a list: Routh array elements

num

a scalar: the number of sign changes

Description

r=routh_t(p) computes Routh's table of the polynomial p.

r=routh_t(h,k) computes Routh's table of denominator of the system described by transfer matrix SISO h with the feedback by the gain k.

If k=poly(0,'k') we will have a polynomial or a rational matrix with dummy variable k, formal expression of the Routh table.

If normalized=%f we will have a polynomial matrix with non normalized elements. In the other case, we will have a rational and normalized matrix.

The second output argument num returns the number of sign changes in the first column of Routh's table. Note that, this argument value will only have sense when the table is normalized.

Hint: If h=1/p, then routh_t(h, kp) is equivalent to routh_t(p+kp) .

Examples

s=%s;
P=5*s^3-10*s^2+7*s+20;
routh_t(P)

// Transfer function with formal feedback, normalized case
routh_t((1+s)/P,poly(0,'k'))

// Transfer function with formal feedback, non normalized case
routh_t((1+s)/P,poly(0,'k'),%f)

// One of the coefficients in the polynomial equals zero
P1=2*s^3-24*s+32;
routh_t(P1)

// A row full of zeros
P2=s^4-6*s^3+10*s^2-6*s+9;
routh_t(P2)

// The number of roots in the rhp could be retrieved as the second output argument
P3=5*s^3-10*s^2+7*s;
[r,num]=routh_t(1/P3,20)
if num==0
   disp("System is stable")
else
   mprintf("There is %g sign changes in entries of first column.\nTherefore, system is unstable.", num)
end

//

See also

• roots — racines d'un polynôme
• spec — valeurs propres, et vecteurs propres d'une matrice ou d'un faisceau de matrices
• evans — Evans root locus
• kpure — continuous SISO system limit feedback gain
• krac2 — continuous SISO system limit feedback gain

Bibliography

http://controls.engin.umich.edu/wiki/index.php/RouthStability

http://www.jdotec.net/s3i/TD_Info/Routh/Routh.pdf

Comments on the Routh-Hurwitz criterion, Shamash, Y.,Automatic Control, IEEE T.A.C Volume 25, Issue 1, Feb 1980 Page(s): 132 - 133

History