Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd Logo
Upload

Welcome to Scribd!

  • 10 views

    Ruthurwis

    Uploaded by

    Hector J. Moscoso Roque
    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.

    Copyright:

    © All Rights Reserved
    Download as TXT, PDF, TXT or read online from Scribd

    Ruthurwis

    Uploaded by

    Hector J. Moscoso Roque
    10 views2 pages
    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.

    Original Title

    ruthurwis

    Copyright

    © © All Rights Reserved

    Available Formats

    TXT, PDF, TXT or read online from Scribd

    Did you find this document useful?

    Is this content inappropriate?

    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.

    Copyright:

    © All Rights Reserved
    Download as TXT, PDF, TXT or read online from Scribd
    Download as txt, pdf, or txt
    10 views2 pages

    Ruthurwis

    Uploaded by

    Hector J. Moscoso Roque
    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.

    Copyright:

    © All Rights Reserved
    Download as TXT, PDF, TXT or read online from Scribd
    Download as txt, pdf, or txt
    of 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