1

INFORME DE LA PRÁCTICA DE LABORATORIO N°3:

“LA TRANSFORMADA DE LAPLACE”
 Adicionalmente se ajusta el grosor de la gráfica para
mejor visualización con el comando line width.
INTRODUCCIÓN
En este informe se describe el desarrollo del laboratorio Se obtiene el siguiente resultado:
la transformada de Laplace, donde se hace uso de la
herramienta de software Matlab para encontrar las exp(-3 t) 3 exp(-t)
transformadas de Laplace y su inversa de diferentes tipos --------- - --------- + 1
de ejercicios y adicionalmente su representación en el 4 4
dominio del tiempo. El uso de dicha herramienta tiene
como fin adquirir conocimientos en el uso de software
Se obtiene la siguiente gráfica:
para la solución de problemas de sistemas de control.

I. OBJETIVOS

Aprender a utilizar la herramienta computacional de

Matlab para encontrar la transformada de la Laplace y su
transformada inversa de diferentes funciones.

II. DESARROLLO DE LA PRÁCTICA

1. Encontrar el voltaje en el capacitor en el tiempo con
la transformada inversa de Laplace y graficar la señal
resultante.

0.5s2  2s  3
Vc (s)  Figura 1. Señal resultante representando el voltaje del
s(s  1)(s  3) capacitor
Se realiza el siguiente código en el script:

2. Encontrar la función de transferencia del siguiente

syms s circuito, aplicar transformada inversa de Laplace y
f=(0.5*(s^2)+2*s+3)/(s*(s+1)*(s+3)); graficar la señal resultante.
F=ilaplace(f);
pretty(F)
set(ezplot(F),'LineWidth',5);

Primero se declara la variable S como simbólica mediante

el comando syms, después se utiliza el comando ilaplace
para obtener la inversa de Laplace, el comando pretty es
para visualizar mejor la respuesta y el comando ezplot se
utiliza para graficar expresiones simbólicas.
Figura 2. Circuito con amplificador operacional
 2

Resultando la función de transferencia de la siguiente exp(-2 s) 2 2 exp(-s) 1

manera:
----------- - --------- + -
Vout 400K
 2 s s s
Vin s  2Ks  2M Se obtiene la siguiente gráfica:
Se realiza el siguiente código de script para encontrar la
transformada inversa de Laplace:

syms s
f=(400e3)/((s^2)+2e3*s+2e6);
F=ilaplace(f);
pretty(F)

Se obtiene el siguiente resultado:

sin(1000 t) exp(-1000 t) 400

3. Dada la siguiente señal utilizando la función

heaviside, encuentre la transformada de Laplace Figura 4. Grafica de v(t) mediante Matlab
utilizando la integral de Laplace y la función escalón.
Seguido, se calcula la transformada de Laplace
recordando que:

F (s)    f (t )   e st f (t )dt
0

Figura 3. Función v(t) representada en escalones unitarios

Se realiza el siguiente código de script utilizando el
comando heaviside para representar escalones unitarios Se realiza el siguiente código de script:
para obtener la gráfica de v(t) y obtener su transformada
de Laplace: syms t s infinito
f1=1;
syms t s f2=-1;
he=[heaviside(t)-2*heaviside(t- f3=1;
1)+2*heaviside(t-2)]; A=int(exp(-s*t)*f1,t,0,1);
set(ezplot(he,[-1,4]),'LineWidth',5); B=int(exp(-s*t)*f2,t,1,2);
F=laplace(he); C=int(exp(-s*t)*f3,t,2,infinito);
pretty(F) pretty(simplify(A+B+C));

Al representar la función a tramos se realiza la integración

mediante el comando int con sus respectivos límites de
2 exp(-s) - exp(-2 s) 2 + exp(-infinito s) - 1
integración, obteniendo el siguiente resultado:
- ------------------------------------------------ Como se aprecia se obtiene el mismo resultado por medio
s de las integrales que con el comando Laplace

Se obtiene el siguiente resultado:

 3

4. Dada la siguiente señal, encuentra la expresión del

tiempo para V(t) en términos de la función heaviside
y su transformada de Laplace

Figura 5. Función v(t)

Se realiza el siguiente código de script utilizando el Figura 6. Grafica de v(t) mediante Matlab
comando heaviside para representar escalones unitarios
para obtener la gráfica de v(t) y obtener su transformada Seguido, se calcula la transformada de Laplace
de Laplace: recordando que:

syms t s F (s)    f (t )   e st f (t )dt
he=[(2*t+1)*(heaviside(t)-heaviside(t- 0
1))+(3)*(heaviside(t-1)-heaviside(t-2))+(- Se realiza el siguiente código de script:
t+3)*(heaviside(t-2)-heaviside(t-3))]; syms t s
set(ezplot(he,[-1,4]),'LineWidth',5); f1=2*t+1;
F=simplify(laplace(he)); f2=3;
pretty(F) f3=-t+3;
A=int(exp(-s*t)*f1,t,0,1);
B=int(exp(-s*t)*f2,t,1,2);
Se obtiene el siguiente resultado:
C=int(exp(-s*t)*f3,t,2,3);
pretty(simplify(A+B+C))
s - 2 exp(-s) - exp(-2 s) + exp(-3 s) - s exp(-2 s) 2 + 2
-----------------------------------------------------------
Al representar la función a tramos se realiza la integración
2
mediante el comando int con sus respectivos límites de
s integración, obteniendo el siguiente resultado:

Se obtiene la siguiente grafica s - 2 exp(-s) - exp(-2 s) + exp(-3 s) - s exp(-2 s) 2 + 2

----------------------------------------------------------------
2
s

5. Calcular las transformadas y las transformadas

inversas de los siguientes ejercicios
 4


a) e4t cos(2 t)  24 (s + 5)
---------------------
b) 2cos(3t)
2
5 ((s + 5) + 160000)
c) f (t)  4.8e5t cos(400t ) Resultado inciso d:
3s  4 / / sqrt(5) t \ \
d) F(s) 
s 2  3s  1 | sqrt(5) sinh| --------- | |
/ 3t\ | / sqrt(5) t \ \ 2 /|
2(s 5)
e) F(s)  exp| - --- | | cosh| --------- | - ------------------------- | 3
(s 1)2 \ 2 / \ \ 2 / 15 /

Se realiza el siguiente código de script utilizando el

comando Laplace para obtener la transformada de Resultado inciso e:
Laplace e ilaplace para obtener la transformada inversa
de Laplace: 2 exp(-t) + 8 t exp(-t)
syms t s
a=exp(-4*t)*cos(2*t); III. CONCLUSIONES
A=laplace(a);
pretty(A)  Se identifican los comandos respectivos para la
escritura de manera simbólica y su correcta
b=2*cos(3*t); manipulación.
B=laplace(b);  Se reconoce la importancia del uso del comando
pretty(B) heaviside para representar funciones a tramos
para su posterior manipulación como aplicar la
c=4.8*exp(-5*t)*cos(400*t); transformada de Laplace.
C=laplace(c);  Es posible contrastar la transformada de Laplace
pretty(C) obtenida por el análisis integral como la obtenida
por los comandos en Matlab, concluyendo, como
d=(3*s+4)/((s^2)+3*s+1); es de esperar, que se llega al mismo resultado.
D=ilaplace(d);
pretty(D)
BIBLIOGRAFIA
e=2*(s+5)/((s+1)^2);
E=ilaplace(e);  [Link]
pretty(E) s/curso/[Link]
 White Robert: Computational Mathematics:
Models, Methods, and Analysis with MATLAB
and MPI. Chapman and Hall/CRC; 1 edition, 2003.
Resultado inciso a:  Andrew Knight Basics of MATLAB and Beyond.
s+4 Chapman and Hall/CRC; 1 edition, 1999.
------------
2
(s + 4) + 4
Resultado inciso b:
2s
------
2
s +9
Resultado inciso c:
5