Análisis de la T.
de Laplace en Python
Juan Felipe Plaza, Valeria Quinto
Pontificia Universidad Javeriana Cali
[email protected] - [email protected]
Enunciado del problema:
Enviar un correo electrónico con la solución de código en extensión .txt y un resumen en
Word o Latex de las conclusiones sobre el sistema que eligieron (por ejemplo, número de
polos, ceros, estabilidad, orden del sistema).
Sistema Seleccionado:
Antes de modificar el código base obtenemos de forma matemática los ceros y polos del
numerador y el denominador respectivamente:
Denominador:
Den=(s +1)( s+ 2)=0
s=−1
s=−2
Numerador:
Núm=4 S+3=0
S=−3/4=−0,75
Por lo tanto se espera que la gráfica en frecuencia tenga 2 polos en -2 y -1 y un cero en
-0,75.
Transformada inversa de Laplace:
Para calcular la transformada inversa de Laplace utilizamos fracciones parciales para
transformar nuestra ecuación de manera que obtenemos:
Utilizamos la siguiente propiedad de la transformada inversa de Laplace:
Y obtenemos la siente ecuación:
Al graficar la función en el dominio del tiempo obtenemos la siguiente gráfica:
De ella podemos concluir que en su estado dinámico el sistema tiene ganancia negativa y
los valores de la función varían de gran manera conforme pasa el tiempo, en su estado
estacionario el sistema se estabiliza en el valor cero debido a que los valores de euler en el
infinito tienden a cero.
Resultados en el código de Python:
Para el sistema seleccionado se tienen los coeficientes para el numerador: [4 , 3] y para el
denominador: [1 3 2]. Siendo de orden 1 y 2 respectivamente.
En la primera parte del código de Python toma los coeficientes y obtiene su función de
transferencia con la función Signal.TransferFunction y con La función Signal.Step se tiene la
respuesta escalonada del sistema en tiempo continuo.
En la segunda parte convierte el sistema en un sistema LTI con la función signal.lti , y
calcula la respuesta en frecuencia del sistema en tiempo continuo con la función
signal.freqresp
Finalmente se obtienen las gráficas de función de transferencia, el diagrama de Nyquist
para la magnitud y el mapa de Polos y Ceros para el análisis en Frecuencia.
Gráfica de Función de Transferencia.
Interpretando la gráfica obtenemos información del sistema, siendo este continuo de orden
2 (el grado del denominador),estrictamente propio (El grado del denominador es mayor que
el del numerador) y también posiblemente es un sistema estable al estar acotado.
Diagrama de Nyquist y Mapa de Polos y ceros.
Los polos y ceros son los esperados.
“Un sistema causal con función del sistema H(s) racional es estable si y sólo si todos los
polos de H(s) caen en la parte izquierda del plano S, es decir, todos los polos tienen parte
real negativa” [1]
Los polos son todos reales negativos, por lo que se encuentran en el lado izquierdo de la
gráfica, por lo tanto el sistema es estable y causal.
Resultados:
Como recapitulación de todo lo que obtuvimos a través de analizar la función y la gráfica
tenemos que:
- Presenta dos polos reales negativos, -2 y -1 y Un cero en -0.75
- Es un sistema Estrictamente Propio, de orden 2 y continua
- El sistema es porque todos sus polos se encuentran a la izquierda en el plano S.
Bibliografía:
[1]A. V. Oppenheim, «Análisis y Caracterización de los sistemas LTI usando la transformada de
Laplace,» de Señales y Sistemas, 1997, pp. 693 - 697.
[2] https://docs.scipy.org/doc/scipy/reference/generated/scipy.signal.freqresp.html
[3] https://docs.scipy.org/doc/scipy/reference/generated/scipy.signal.step.html
[4] https://docs.scipy.org/doc/scipy/reference/generated/scipy.signal.TransferFunction.html
[5] https://es.wikipedia.org/wiki/Diagrama_de_Nyquist#:~:text=El%20diagrama%20de%20Nyquist
%20es,de%20un%20sistema%20con%20realimentaci%C3%B3n.
[6] https://dademuch.com/2020/02/15/el-diagrama-de-nyquist/
