Transformadas Integrales

Ingeniería en Diseño Electrónico

XXxx
XXxx

Tarea 1: La transformada de Laplace:

Importancia y aplicaciones
Módulo 1: Transformada de Laplace

XXxx

XXxx
Actividad I
PARTE 1
Una empresa se encuentra desarrollando un prototipo robótico el cual incluye un
mecanismo que se puede modelar como un sistema masa-resorte-amortiguador. Con
las características físicas del sistema con 𝑚𝑚 = 1.0 𝑘𝑘𝑘𝑘, 𝑘𝑘 = 0.1, y 𝑏𝑏 = 0.2. Se obtiene la
siguiente ecuación diferencial:
𝑑𝑑 2 𝑥𝑥 𝑑𝑑𝑑𝑑
𝑚𝑚 2
+ 𝑏𝑏 + 𝑘𝑘𝑘𝑘 = 𝑢𝑢 (1)
𝑑𝑑𝑡𝑡 𝑑𝑑𝑑𝑑

Utilice la transformada de Laplace para resolver la ecuación diferencial y determinar

cuál es la función que representa la posición del sistema
• Método Analítico
Aplicando la transformada de Laplace a la ecuación 1,
𝑑𝑑 2 𝑥𝑥 𝑑𝑑𝑑𝑑
ℒ �𝑚𝑚 ∙ 2 + 𝑏𝑏 + 𝑘𝑘 ∙ 𝑥𝑥� = ℒ{𝑢𝑢}
𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑

Utilizando las propiedades de la transformada

𝑑𝑑 2 𝑥𝑥 𝑑𝑑𝑑𝑑
𝑚𝑚 ∙ ℒ � 2 � + 𝑏𝑏 ∙ ℒ � � + 𝑘𝑘 ∙ ℒ{𝑥𝑥} = ℒ{𝑢𝑢}
𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑
2
𝑚𝑚 ∙ [𝑠𝑠 𝑋𝑋(𝑠𝑠) − 𝑠𝑠𝑠𝑠(0) − 𝑥𝑥̇ (0)] + 𝑏𝑏 ∙ [𝑠𝑠 ∙ 𝑋𝑋(𝑠𝑠) − 𝑥𝑥(0)] + 𝑘𝑘 ∙ 𝑋𝑋(𝑠𝑠) = 𝑈𝑈(𝑠𝑠)
Si las condiciones iniciales 𝑥𝑥(0) = 0 y 𝑥𝑥̇ (0) = 0, la ecuación anterior quedaría de la
siguiente manera:
𝑚𝑚 ∙ [𝑠𝑠 2 ∙ 𝑋𝑋(𝑠𝑠)] + 𝑏𝑏 ∙ [𝑠𝑠 ∙ 𝑋𝑋(𝑠𝑠)] + 𝑘𝑘 ∙ 𝑋𝑋(𝑠𝑠) = 𝑈𝑈(𝑠𝑠)
Reduciendo y simplificando esta ecuación
𝑋𝑋(𝑠𝑠) ∙ (𝑚𝑚 ∙ 𝑠𝑠 2 + 𝑏𝑏 ∙ 𝑠𝑠 + 𝑘𝑘) = 𝑈𝑈(𝑠𝑠)
1
Despejando 𝑋𝑋(𝑠𝑠), si 𝑈𝑈(𝑠𝑠) = y dividiendo el denominador entre 𝑚𝑚
𝑠𝑠

1
𝑋𝑋(𝑠𝑠) =
𝑏𝑏 𝑘𝑘 (2)
𝑠𝑠 �𝑠𝑠 2 + ∙ 𝑠𝑠 + �
𝑚𝑚 𝑚𝑚

Sustituyendo los valores dados al inicio del problema, en la ecuación 2

1
𝑋𝑋(𝑠𝑠) =
𝑠𝑠(𝑠𝑠 2 + 0.2 ∙ 𝑠𝑠 + 0.1) (3)

pág. 1
Para aplicar la transformada inversa de Laplace, primero es necesario descomponer la
ecuación 3 en fracciones parciales. Para ello:
1 𝐴𝐴 𝐵𝐵𝑠𝑠 + 𝐶𝐶
= +
𝑠𝑠(𝑠𝑠 2 + 0.2 ∙ 𝑠𝑠 + 0.1) 𝑠𝑠 (𝑠𝑠 2 + 0.2 ∙ 𝑠𝑠 + 0.1) (4)

Resolviendo las fracciones parciales

1 = 𝐴𝐴 ∙ (𝑠𝑠 2 + 0.2 ∙ 𝑠𝑠 + 0.1) + [(𝐵𝐵𝑠𝑠 + 𝐶𝐶) ∙ 𝑠𝑠]
1 = (𝐴𝐴 + 𝐵𝐵) ∙ 𝑠𝑠 2 + (0.2 ∙ 𝐴𝐴 + 𝐶𝐶) ∙ 𝑠𝑠 + 0.1 ∙ 𝐴𝐴
El sistema de ecuaciones quedaría
(𝐴𝐴 + 𝐵𝐵) = 0
(0.2 ∙ 𝐴𝐴 + 𝐶𝐶) = 0
0.1 ∙ 𝐴𝐴 = 1
Así encontramos que 𝐴𝐴 = 10; 𝐵𝐵 = −10 y 𝐶𝐶 = −2
Sustituyendo estos valores encontrados en la ecuación 4.
1 10 10𝑠𝑠 + 2
𝑋𝑋(𝑠𝑠) = = − 2
𝑠𝑠(𝑠𝑠 2 + 0.2 ∙ 𝑠𝑠 + 0.1) 𝑠𝑠 (𝑠𝑠 + 0.2 ∙ 𝑠𝑠 + 0.1) (5)

Aplicando transformada inversa de Laplace a la ecuación 5

10 10𝑠𝑠 + 2
ℒ −1 {𝑋𝑋(𝑠𝑠)} = ℒ −1 � − 2 �
𝑠𝑠 (𝑠𝑠 + 0.2 ∙ 𝑠𝑠 + 0.1)
1 10𝑠𝑠 2
ℒ −1 {𝑋𝑋(𝑠𝑠)} = 10 ∙ ℒ −1 � � − �ℒ −1 � 2 � + ℒ −1 � 2 ��
𝑠𝑠 (𝑠𝑠 + 0.2 ∙ 𝑠𝑠 + 0.1) (𝑠𝑠 + 0.2 ∙ 𝑠𝑠 + 0.1)
Para aplicar las formulas de las tablas que se encuentran en el apéndice es necesario
a completar los términos en la transformada inversa, así se tiene que:
1 𝑠𝑠 + 0.1 10 −1 0.3
ℒ −1 {𝑋𝑋(𝑠𝑠)} = 10 ∙ ℒ −1 � � − �10 ∙ ℒ −1 � � + ∙ ℒ � ��
𝑠𝑠 (𝑠𝑠 + 0.1)2 + (0.3)2 3 (𝑠𝑠 + 0.1)2 + (0.3)2
Simplificando
1 𝑠𝑠 + 0.1 1 −1 0.3
ℒ −1 {𝑋𝑋(𝑠𝑠)} = 10 ∙ ℒ −1 � � − 10 ∙ �ℒ −1 � � + ∙ ℒ � ��
𝑠𝑠 (𝑠𝑠 + 0.1)2 + (0.3)2 3 (𝑠𝑠 + 0.1)2 + (0.3)2

pág. 2
Aplicando la transformada inversa de Laplace
1
𝑥𝑥(𝑡𝑡) = 10 − �10 ∙ 𝑒𝑒 −0.1∙𝑡𝑡 ∙ �cos(0.3 ∙ 𝑡𝑡) + ∙ sin (0.3 ∙ 𝑡𝑡)�� (6)
3
• Método Numérico
Para la resolución de este problema ahora se planteo el uso del software de Matlab
donde se programó el script que se muestra en la tabla 1.
% Borrar todo y cerrar ventanas
clc
clear all
close all

% variables simbolicas
syms s t

% Parametros del sistema

m = 1;
b = 0.2;
k = 0.1;
F = 1;

% Resolucion de la ecuacion diferencial

solucion = dsolve('1*D2x + 0.2*Dx + 0.1*x = 1', 'x(0) = 0', 'Dx(0) = 0', 't');
disp("La solucion de la ecuacion diferencial es: ")
simplify(solucion)
figure(1)
fplot(solucion, [0, 50])
title("Ecuacion diferencial metodo numerico")

% Resolucion por Laplace

num = F/s;
den = m*s^2 + b*s + k;
Xs = num / den;
x = ilaplace(Xs);
disp("La solucion de la ecuacion diferencial es: ")
simplify(x)
figure(2)
fplot(x, [0, 50])
title("Ecuacion diferencial por Laplace")

Tabla 1. Script de Matlab que calcula la respuesta de un sistema Resorte-Masa-Amortiguador

Como resultado también podemos calcular la solución de la transformada de Laplace

de manera simbólica, tal y como se muestra en la figura 1. Donde se aprecia el mismo
resultado obtenida de manera analítica.

Figura 1. Resultado obtenido de manera numérica

pág. 3
Por ultimo se muestra la respuesta del sistema el cual corresponde a uno de tipo
subamortiguado. Ambas respuestas ya sea usando la instrucción dsolve o iLaplace
son de tipo numéricas lo cual da por resultado que la respuesta sea similar.

Figura 2. Respuesta numérica a problema de Resorte Masa Amortiguador

PARTE 2
Resuelve el siguiente circuito RLC en serie

Figura 3 Circuito RLC

• Método Analítico
Para resolver es te circuito de acuerdo a la ley de mallas de Kirchoff, donde menciona
que la suma algébrica de todas las caídas de voltaje en una malla, bucle o circuito
cerrado como el que se muestra en la figura 1 se puede expresar de la siguiente
manera:

� 𝑉𝑉𝑖𝑖 = 0

pág. 4
Rescribiendo la ecuación de malla de Kirchoff
𝑉𝑉𝑅𝑅 + 𝑉𝑉𝐶𝐶 + 𝑉𝑉𝐼𝐼 − 𝑣𝑣(𝑡𝑡) = 0
Sustituyendo las caídas de voltaje en cada elemento
𝑞𝑞 𝑑𝑑𝑑𝑑
𝑖𝑖𝑖𝑖 +
+ 𝐿𝐿 − 𝑣𝑣(𝑡𝑡) = 0 (7)
𝐶𝐶 𝑑𝑑𝑑𝑑
Si la corriente está definida como el cambio de carga de un elemento a través del
tiempo,
𝑑𝑑𝑑𝑑
𝑖𝑖 =
𝑑𝑑𝑑𝑑
Si introducimos esta relación a la ecuación (7), obtenemos
𝑑𝑑 2 𝑞𝑞 𝑑𝑑𝑑𝑑 1
𝐿𝐿 ∙ 2
+ 𝑅𝑅 ∙ + ∙ 𝑞𝑞 = 𝑣𝑣(𝑡𝑡) (8)
𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑 𝐶𝐶

Aplicando la transformada de Laplace a la ecuación 8,

𝑑𝑑 2 𝑞𝑞 𝑑𝑑𝑞𝑞 1
ℒ �𝐿𝐿 ∙ 2
+ 𝑅𝑅 + ∙ 𝑞𝑞� = ℒ{𝑣𝑣(𝑡𝑡)}
𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑 𝐶𝐶

Utilizando las propiedades de la transformada

𝑑𝑑2 𝑞𝑞 𝑑𝑑𝑞𝑞 1
𝐿𝐿 ∙ ℒ � � + 𝑅𝑅 ∙ ℒ � � + ∙ ℒ{𝑞𝑞} = ℒ{𝑣𝑣(𝑡𝑡)}
𝑑𝑑𝑑𝑑 2 𝑑𝑑𝑑𝑑 𝐶𝐶
1
𝐿𝐿 ∙ [𝑠𝑠 2 𝑄𝑄(𝑠𝑠) − 𝑠𝑠𝑞𝑞(0) − 𝑞𝑞̇ (0)] + 𝑅𝑅 ∙ [𝑠𝑠 ∙ 𝑄𝑄(𝑠𝑠) − 𝑞𝑞(0)] + ∙ 𝑄𝑄(𝑠𝑠) = 𝑉𝑉(𝑠𝑠)
𝐶𝐶
Si las condiciones iniciales 𝑞𝑞(0) = 0 y 𝑞𝑞̇ (0) = 0, la ecuación anterior quedaría de la
siguiente manera:
1
𝐿𝐿 ∙ [𝑠𝑠 2 ∙ 𝑄𝑄(𝑠𝑠)] + 𝑅𝑅 ∙ [𝑠𝑠 ∙ 𝑄𝑄(𝑠𝑠)] + ∙ 𝑄𝑄(𝑠𝑠) = 𝑉𝑉(𝑠𝑠)
𝐶𝐶
Reduciendo y simplificando esta ecuación
1
𝑄𝑄(𝑠𝑠) ∙ �𝐿𝐿 ∙ 𝑠𝑠 2 + 𝑅𝑅 ∙ 𝑠𝑠 + � = 𝑉𝑉(𝑠𝑠)
𝐶𝐶
Despejando 𝑄𝑄(𝑠𝑠), y dividiendo entre 𝐿𝐿
1
𝑄𝑄(𝑠𝑠) = ∙ 𝑉𝑉(𝑠𝑠)
𝑅𝑅 1 (9)
𝑠𝑠 2 + ∙ 𝑠𝑠 +
𝐿𝐿 𝐿𝐿𝐿𝐿

pág. 5
Si para 𝑣𝑣(𝑡𝑡) = 24 ∙ sin (10𝑡𝑡), a partir de las tablas en el apéndice, su transformada de
Laplace estaría dada por
10
𝑉𝑉(𝑠𝑠) = 48 ∙ � �
𝑠𝑠 2 + 102

Sustituyendo los valores dados donde 𝐿𝐿 = 0.5ℎ, 𝑅𝑅 = 6Ω, 𝐶𝐶 = 0.02𝑓𝑓 y 𝑉𝑉(𝑠𝑠)

1 480
𝑄𝑄(𝑠𝑠) = ∙ 2
𝑠𝑠 2 + 12 ∙ 𝑠𝑠 + 100 𝑠𝑠 + 100

480
𝑄𝑄(𝑠𝑠) =
(𝑠𝑠 2 + 100)(𝑠𝑠 2 + 12 ∙ 𝑠𝑠 + 100) (10)

Para aplicar la transformada inversa de Laplace, primero es necesario descomponer la

ecuación 10 en fracciones parciales. Para ello:
480 𝐴𝐴𝑠𝑠 + 𝐵𝐵 𝐶𝐶𝑠𝑠 + 𝐷𝐷
= +
(𝑠𝑠 2 + 100)[(𝑠𝑠 + 6)2 + 64] (𝑠𝑠 2 + 100) [(𝑠𝑠 + 6)2 + 64] (11)

Resolviendo las fracciones parciales

480 = (𝐴𝐴𝑠𝑠 + 𝐵𝐵) ∙ (𝑠𝑠 2 + 12𝑠𝑠 + 100) + (𝐶𝐶𝐶𝐶 + 𝐷𝐷) ∙ (𝑠𝑠 2 + 100)
480 = 𝐴𝐴𝑠𝑠 3 + 12𝐴𝐴𝑠𝑠 2 + 100𝐴𝐴𝐴𝐴 + 𝐵𝐵𝑠𝑠 2 + 12𝐵𝐵𝐵𝐵 + 100𝐵𝐵 + 𝐶𝐶𝑠𝑠 3 + 100𝐶𝐶𝐶𝐶 + 𝐷𝐷𝑠𝑠 2 + 100𝐷𝐷
El sistema de ecuaciones quedaría
(𝐴𝐴 + 𝐶𝐶) = 0
(12𝐴𝐴 + 𝐵𝐵 + 𝐷𝐷) = 0
100𝐴𝐴 + 12𝐵𝐵 + 100𝐶𝐶 = 0
100𝐵𝐵 + 100𝐷𝐷 = 480
Así encontramos que 𝐴𝐴 = −0.4; 𝐵𝐵 = 0; 𝐶𝐶 = 0.4; 𝐷𝐷 = 4.8
Sustituyendo estos valores encontrados en la ecuación 11.
480 −0.4 ∙ 𝑠𝑠 0.4 ∙ 𝑠𝑠 + 4.8 (12)
= +
(𝑠𝑠 2 + 100)[(𝑠𝑠 + 6)2 + 64] (𝑠𝑠 2 + 100) [(𝑠𝑠 + 6)2 + 64]

pág. 6
Aplicando transformada inversa de Laplace a la ecuación 12
−0.4 ∙ 𝑠𝑠 0.4 ∙ 𝑠𝑠 + 4.8
ℒ −1 {𝑄𝑄(𝑠𝑠)} = ℒ −1 � 2 + �
(𝑠𝑠 + 100) [(𝑠𝑠 + 6)2 + 64]
𝑠𝑠 𝑠𝑠
ℒ −1 {𝑄𝑄(𝑠𝑠)} = −0.4 ∙ ℒ −1 � � + 0.4 ∙ ℒ −1 � �
(𝑠𝑠 2
+ 100) [(𝑠𝑠 + 6)2 + 64]
1
+ 4.8ℒ −1 � �
[(𝑠𝑠 + 6)2 + 64]
Para aplicar las fórmulas de las tablas que se encuentran en el apéndice es necesario
a completar los términos en la transformada inversa, así se tiene que:
𝑠𝑠 𝑠𝑠 + 6
ℒ −1 {𝑄𝑄(𝑠𝑠)} = −0.4 ∙ ℒ −1 � � + 0.4 ∙ ℒ −1
� �
(𝑠𝑠 2 + 100) [(𝑠𝑠 + 6)2 + 64]
2.4 −1 8
+ ℒ � �
8 [(𝑠𝑠 + 6)2 + 64]
Aplicando la transformada inversa de Laplace
𝑞𝑞(𝑡𝑡) = −0.4 ∙ cos(10𝑡𝑡) + 0.1 ∙ 𝑒𝑒 −6𝑡𝑡 ∙ [4 ∙ cos(8𝑡𝑡) + 3 ∙ sin (8𝑡𝑡)] (13)
Para conocer la corriente del circuito solo es necesario derivar la ecuación 13
𝑑𝑑𝑑𝑑 (14)
𝑖𝑖(𝑡𝑡) = = 4 ∙ sin(10𝑡𝑡) − 5𝑒𝑒 −6𝑡𝑡 ∙ sin (8𝑡𝑡)
𝑑𝑑𝑑𝑑

• Método Numérico
Para la resolución de este problema se hizo del uso del software de Matlab donde se
programó el script que se muestra en la tabla 1.
% Borrar todo y cerrar ventanas
clc
clear all
close all

% variables simbolicas
syms s t

% Parametros del sistema

L = 0.5;
R = 6;
C = 0.02;
vt = 24*sin(10*t);

% Resolucion de la ecuacion diferencial

solucion = dsolve('0.5*D2x + 6*Dx + 50*x = 24*sin(10*t)', 'x(0) = 0', 'Dx(0) = 0', 't');
disp("La solucion de la ecuacion diferencial es: ")
simplify(solucion)
figure(1)
fplot(solucion, [0, 20])
title("Ecuacion diferencial metodo numerico")

pág. 7
% Resolucion por Laplace
num = 48*(10/(s^2+100));
den = L*s^2 + R*s + 1/C;
Xs = num / den;
x = ilaplace(Xs);
disp("La solucion de la ecuacion diferencial es: ")
simplify(x)
figure(2)
fplot(x, [0, 20])
title("Ecuacion diferencial por Laplace")

Tabla 2. Script de Matlab que calcula la respuesta de un circuito RLC

Como resultado también podemos calcular la solución de la transformada de Laplace

de manera simbólica, tal y como se muestra en la figura 5. Donde se aprecia el mismo
resultado obtenida de manera analítica.

Figura 4 Resultado obtenido de manera numérica

Por último, se muestra la respuesta del sistema el cual corresponde a uno de tipo
oscilante. Ambas respuestas ya sea usando la instrucción dsolve o iLaplace son de
tipo numéricas lo cual da por resultado que la respuesta sea similar.

Figura 5. Respuesta numérica a problema de circuito RLC

pág. 8
Conclusiones
La transformada de Laplace y su inversa son herramientas fundamentales en el análisis
de sistemas dinámicos en ingeniería, permitiendo convertir ecuaciones diferenciales
complejas en algebraicas más manejables. Su importancia radica en la capacidad de
analizar y resolver sistemas de control, como el resorte-masa-amortiguador y los
circuitos RLC, de manera eficiente.
El cálculo analítico proporciona soluciones exactas y comprensión teórica profunda,
mientras que el uso de MATLAB para cálculos numéricos ofrece rapidez y capacidad
de manejar sistemas más complejos, aunque con aproximaciones. La combinación de
ambos métodos asegura soluciones precisas y prácticas en ingeniería.

Bibliografía

• Vázquez N., J. (2023). Historia de la transformada de Laplace. Telecomunicado

(2021). ENTIENDE con 7 CLAVES la TRANSFORMADA DE LAPLACE [para
PRINCIPIANTES]. [Vídeo]. YouTube. Recuperado de:
https://www.youtube.com/watch?v=0WXQbYDH4yI
• El traductor de ingeniería (2015). Transformada de Laplace - Parte 1 de 3 | El
Traductor. [Vídeo]. YouTube. Recuperado
de: https://www.youtube.com/watch?v=McSN9g7DbYA
• Vázquez N., J. (2023). Resolviendo ecuaciones diferenciales con Matlab y la
transformada de Laplace. [Video]. Anahuac Online.
• Vázquez N., J. (2023). Transformada inversa de Laplace con Matlab. [Video].
Anahuac Online.

pág. 9
Apéndice 1

pág. 10