Scribd
Upload
17 views2 pages

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
17 views2 pages

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

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy