syslin

Programa en Scilab

s=%s

num=s^2+2*s+3;

den=(s+1)^3;

g=syslin('c',num/den);

gf=tf2ss(g);

se=pfss(gf)


Solucion:

 se  =
 
 
       se(1)
 
                2     
      3 + 2s + s      
    --------------    
               2   3  
    1 + 3s + 3s + s

Como vemos no nos ha solucionado nada debido al polo multiple, por lo que lo resolveremos como lo hacemos normalmente pero utilizando Scilab. Es decir solucionaremos las siguienter ecuaciones, para obtener los coeficientes:

La descomposicion en fracciones simples nos quedaria:

Esta ecuacion la programaremos con Scilab de la siguiente manera:

 

Programa en Scilab

s=%s

num=s^2+2*s+3;

den=(s+1)^3;

g=syslin('c',num/den);

rd=roots(den);

[n d k]=factors(g);

a(3)=horner(g*d(1)^3,rd(1))/2;

a(2)=horner(derivat(g*d(1)^3),rd(1));

a(1)=horner(derivat(derivat(g*d(1)^3)),rd(1))

 

 Solucion de los coeficientes:

 a  =
 
    2.  
    0   
    1.