Scribd
Upload
13 views

Ruthurwis

The document contains code to analyze the stability of a polynomial using the Routh-Hurwitz stability criterion. It takes a polynomial as input, constructs the Routh array, and calculates the number of right half plane poles to determine stability. It also has options to display the roots of the polynomial or plot the transfer function and step response.

Uploaded by

Hector J. Moscoso Roque
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as TXT, PDF, TXT or read online on Scribd
13 views

Ruthurwis

The document contains code to analyze the stability of a polynomial using the Routh-Hurwitz stability criterion. It takes a polynomial as input, constructs the Routh array, and calculates the number of right half plane poles to determine stability. It also has options to display the roots of the polynomial or plot the transfer function and step response.

Uploaded by

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

clc

clear
r=input('ingrse el polinomio : ');
num=input('ingrese el numerador : ');
m=length(r);
n=round(m/2);
q=1;
k=0;
for p = 1:m
if rem(p,2)==0
c_even(k)=r(p);
else
c_odd(q)=r(p);

k=k+1;
q=q+1;
end
end
a=zeros(m,n);

if m/2 ~= round(m/2)
c_even(n)=0;
end
a(1,:)=c_odd;
a(2,:)=c_even;
if a(2,1)==0
a(2,1)=0.01;
end
for i=3:m
for j=1:n-1
x=a(i-1,1);
if x==0
x=0.01;
end

a(i,j)=((a(i-1,1)*a(i-2,j+1))-(a(i-2,1)*a(i-1,j+1)))/x;

end
if a(i,:)==0
order=(m-i+1);
c=0;
d=1;
for j=1:n-1
a(i,j)=(order-c)*(a(i-1,d));
d=d+1;
c=c+2;
end
end
if a(i,1)==0
a(i,1)=0.01;
end
end
Right_poles=0;
for i=1:m-1
if sign(a(i,1))*sign(a(i+1,1))==-1
Right_poles=Right_poles+1;
end
end
fprintf('\n Tabla de Routh-Hurwitz:\n')
a
fprintf('\n numero de ceros Polos =%2.0f\n',Right_poles)

reply = input('�Necesita ra�ces del sistema? Y/N ', 's');


if reply=='y'||reply=='Y'
ROOTS=roots(r);
fprintf('\n Given Polynomials Coefficents Roots :\n')
ROOTS
else
end
tf(num,r)
step(num,r)

You might also like