Transformada de Laplace

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 12

Transformada de Laplace

Si f (x, y) es una función de dos variables, la integral definida de la función f (x, y)


respecto a una de las variables produce una función de la otra variable, por ejemplo
Z 2
3x2 y 3 dx = 8y 3 .
0

Por esta razón a este tipo de integrales se les conoce como transformadas inte-
grales. Particularmente, en esta sección estudiamos una transformada integral
impropia muy importante llamada transformada de Laplace.
Sea f (x) una función definida para cada x ≥ 0. Se define la transformada de
Laplace de f (x) como
Z ∞
L [f (x)] = e−px f (x)dx (1)
0
Z b
= lı́m e−px f (x)dx,
b→∞ 0

si el lı́mite existe.
Cuando la integral impropia (1) converge, el resultado es una función de p. En
la práctica empleamos letras minúsculas para representar las funciones a que se
van a transformar (entradas) y mayúsculas para denotar sus transformadas de
Laplace(salidas), ejemplo

L[f (x)] = F (p), L[g(x)] = G(p), L[h(x)] = H(p).

Ejemplo 2. Evalúe L[1].

Solución:
Z ∞ Z b
−px
L[1] = e
(1)dx = lı́m e−px dx
0 b→∞
−px b
 0 −pb 
−e −e + 1 1
= lı́m = lı́m = ,
b→∞ p 0 b→∞ p p
si p > 0.

Ejemplo 3. Evalúe L[x].


2

Solución:
Z ∞ Z b
−px
L[x] = e xdx = lı́m xe−px dx (integramos por partes:u = x, dv = e−px dx)
0 b→∞ 0
 −px b Z b −px 
−xe e 1 1
= lı́m + dx = L[1] = 2 ,
b→∞ p 0 0 p p s
si p > 0.
Ejemplo 4. Evalúe L[e−2x ].
Solución:
∞ b
e−(p+2)x b
Z Z
−2x −px −2x
L[e ]= e e dx = lı́m e−(p+2)x dx = lı́m −
p+2 0

0 b→∞ 0 b→∞
−(p+2)b
e 1 1
= lı́m − + = ,
b→∞ p+2 p+2 p+2
si p+3 > 0, equivalentemente p > −3. ¿Qué podrá decir de L[eax ], a número real?

Ejemplo 5. Evalúe L[cos kx].


Solución:
Z ∞ Z b
−px
L[cos kx] = e cos kxdx = lı́m e−px cos kxdx, (u = cos kx, dv = e−px dx)
0 b→∞
Z 0
cos kx −px b k ∞ −px
= lı́m − e − e sen kxdx, (u = sen kx, dv = e−px dx)
b→∞ p 0 p 0
sen kx −px k ∞ −px
 Z 
1 k
= − lı́m − e + e cos kxdx
p p b→∞ p p 0
1 k2
= − 2 L[cos kx]
p p
De aquı́ obtenemos que
k2
 
1
1 + 2 L[cos kx] =
p p
 2
p + k2

1
L[cos kx] =
p2 p
p
L[cos kx] = , p > 0.
p2 + k 2
Observación 6. De aquı́ en adelante, ya no citaremos restricciones de p, se
supone que p está restringida de manera que la transformada de Laplace converja.
La siguiente tabla muestra las transformadas de Laplace de algunas funciones
elementales.
3

Función Transformada de Laplace


f (x) = 1 F (p) = p1
f (x) = x F (p) = p12
n!
f (x) = xn F (p) = pn+1
1
f (x) = eax F (p) = p−a
a
f (x) = sen ax F (p) = p2 +a2
p
f (x) = cos ax F (p) = p2 +a2
a
f (x) = senh ax F (p) = p2 −a2
p
f (x) = cosh ax F (p) = p2 −a2

De la definición de transformada de Laplace, es fácil demostrar que L es un


operador lineal, es decir, si f y g son funciones como en la definición anterior
y α es una constante,

i) L [f (x) + g(x)] = L [f (x)] + L [g(x)]

ii) L [αf (x)] = aL [f (x)].

La propiedad de linealidad de la trasformada de Laplace es importante ya que a


través de esta podemos determinar la transformada de combinaciones lineales de
las funciones elementales, Por ejemplo,
3 7 6
Ejemplo 7. ˆ L[3 + 7x + x3 ] = 3L[1] + 7L[x] = + 2+ 4
p p p
4 20
ˆ L[4e−3x − 10 sen 2x] = 4L[e−3x ] − 10L[sen 2x] = − 4
p+3 p +4
El siguiente teorema nos permite saber si dada una función existe o no su
transformada de Laplace.

Teorema. Si f es una función continua a pedazos para x ≥ 0 y |f (x)| ≤ M ecx


para cada x ≥ T , donde M , c > 0 y T > 0 son constantes, entonces L [f (x)] (p)
existe para p > c.

Cuando f (x) ≤ |f (x)| ≤ M ecx , decimos que f (x) es de orden exponencial.


Más aún

Teorema. Si f es una función continua a pedazos para x ≥ 0 y de orden expo-


nencial para x ≥ T , entonces

lı́m L [f (x)] (p) = lı́m F (p) = 0.


p→∞ p→∞


 1 si 0 ≤ x < 2,
Ejemplo 8. Evalúe L[f (x)] cuando f (x) = 0 si 2 ≤ x < 4,
2 si x ≥ 4

4

Solución: Como la función está definidad a trozos tenemos


Z ∞
L[f (t)] = e−px f (x)dx
Z0 2 Z 4 Z ∞
−px −px
= e (1)dx + e (0)x + e−px (2)dx
0 2 4
e−px 2 2e−px b
=− + lı́m −
p 0 b→∞ p 0

e−2p 1 2
=− + +
p p p
1
= (3 − e−2p ).
p

Transformada inversa de Laplace


Si L [f (x)] = F (p), decimos que f (x) es una transformada inversa de Laplace de
F (p) y la denotaremos por

L−1 [F (p)] = f (x).

En general, la transformada inversa de Laplace de F (p), no necesariamente es


única.
Naturalmente, tenemos las siguientes trasformadas inversas de funciones ele-
mentales:
Función: F (p) Transformada inversa de Laplace: L−1 [F (p)]
F (p) = p1 f (x) = 1
F (p) = p12 f (x) = x
n!
F (p) = pn+1 f (x) = xn
1
F (p) = p−a f (x) = eax
a
F (p) = p2 +a 2 f (x) = sen ax
p
F (p) = p2 +a2 f (x) = cos ax
a
F (p) = p2 −a 2 f (x) = senh ax
p
F (p) = p2 −a2 f (x) = cosh ax

Cuando evaluamos transformadas inversas con frecuencia sucede que la función de


p en consideración no corresponde exactamente a la forma de una transformada
de Laplace F (p) como aparece en la tabla. Por lo que será necesario adecuar la
función de p ya sea multiplicando o dividiendo por una constante apropiada como
lo ilustran los siguientes ejemplos:

Ejemplo 9. Determinar las siguiente transformadas inversas:


5
   
−1 3 −1 2
a) L b) L
p7 2
p +7
Solución:
  
−13 3 −1 6! 3
L 7
= L 6+1
= x6 .
p 6! p 6!
" √ #

 
2 2 7 2
L−1 2 = √ L−1 2 = √ sen( 7x).
p +7 7 p +7 7

Si las funciones son continuas, L−1 es un operador lineal, es decir


i) L−1 [F (p) + G(p)] = L−1 [F (p)] + L−1 [G(p)]

ii) L−1 [αF (p)] = αL−1 [F (p)], donde α es una constante.


 
−1 −2p + 6
Ejemplo 10. Determinar L
p2 − 4
Solución:
     
−1 −2p + 6 −1 p −1 2
L = −2L + 3L = −2 cosh 2x + 3 sinh 2x
p2 − 4 p2 − 4 p2 − 4
Aunque en teorı́a parece sencillo determinar la transformada inversa de una función,
en la práctica suele no serlo. Sin embargo, el método de fracciones parciales de-
sempeña un papel fundamental para dicha tarea.

Transformada de una derivada y solución a EDOs


lineales
Nuestro objetivo inmediato es utilizar la transformada de Laplace para resolver
ecuaciones diferenciales. Para ello necesitamos evaluar funciones como:
   00   (n) 
dy dy dy
L ,L 2
,...,L .
dx dx dxn
dy
Si y 0 (x) = es continua para x ≥ 0, entonces
dx
Z ∞
0
L[y (x)] = e−px y 0 (x)dx,
0

integrando por por partes tenemos


b
Z ∞
0
lı́m y(x)e−px 0 + p e−px y(x)dx

L[y (x)] =
b→∞ 0
= pL[y(x)] − y(0).
6

Como y 00 (x) = (y 0 (x))0 , entonces


0
L[y 00 (x)] = L[(y 0 (x)) ]
= pL[y 0 (x)] − y 0 (0)
= p [pL[y(x)] − y(0)] − y 0 (0)
= p2 L[y(x)] − py(0) − y 0 (0).
Continuando con este proceso, podemos generalizar el resultado anterior, como
lo indica el siguiente teorema.
Teorema. Si y, y 0 , . . . , y (n−1) son funciones en [0, ∞) de orden exponencial y si
y (n) (x) es continua por tramos en [0, ∞), entonces
L[y (n) (x)] = pn Y (p) − p(n−1) y(0) − p(n−2) y 0 (0) − ... − y (n−1) (0),
donde Y (p) = L[y(x)].
Este teorema hace que la tranformada de Laplace sea ideal para resolver proble-
mas de valor inicial en los que la ecuación diferencial tenga coeficientes constantes:
an y (n) + an−1 y (n−1) + · · · + a0 y = g(x) (11)
y(0) = y0 , y 0 (0) = y1 , . . . , y (n−1) (0) = yn−1
donde las ai , i = 0, 1, . . . , n y yk , k = 0, 1, . . . , n − 1 son constantes.
En efecto, aplicando la tranformada de Laplace a la ecuación (11):
an L[y (n) ] + an−1 L[y (n−1) ] + · · · + a0 L[y] = L[g(x)]
de donde obtenemos
an pn Y (p) − p(n−1) y(0) − p(n−2) y 0 (0) − ... − y (n−1) (0) +
 
(12)
an−1 p Y (p) − p(n−2) y(0) − p(n−2) y 0 (0) − ... − y (n−1) (0) + · · · + a0 Y (p) = G(p),
 n−1 

donde L[y(x)] = Y (p) y L[g(x)] = G(p).


Observe que la trasformada de Laplace de una ecuación diferencial lineal con
coeficientes constantes se transforma en una ecuación algebraica en Y (p).
Resolviendo (12) para determinar Y (p), obtenemos:
R(p) G(p)
Q(p)Y (p) = R(p) + G(p), equivalentemente Y (p) = + ,
Q(p) Q(p)
donde Q(p) = an pn + an−1 pn−1 + · · · + a1 p + a0 y R(p) es un polinomio en p de
grado menor o igual a n − 1, formado por diversos productos de los coeficientes
ai y las condiciones iniciales yk .
Por último, obtenemos la solución del problema de valor inicial
y(x) = L−1 [Y (p)].
Resumen del método para resolver un problema de valor inicial como (11):
7

1) Aplicar la trasformada de Laplace L, a la ecuación diferencial con condi-


ciones iniciales, para obtener una ecuación algebraica en Y (p).

2) Resolver la ecuación algebraica para determinar Y (p).

3) Aplicar la transformada inversa L−1 a Y (p) para determinar y(x).


Ejemplo 13. Resolver los siguientes problemas de valor inicial mediante la trans-
formada de Laplace
dy
1. − y = 1, y(0) = 0
dx
2. y 00 − 3y 0 + 2y = e−4x , y(0) = 1, y 0 (0) = 5

3. y 000 + 2y 00 − y 0 − 2y = sen 3x, y(0) = 0, y 0 (0) = 0, y 00 (0) = 1

Otras propiedades de la Transformada de Laplace


Traslación en el eje p
Otra propiedad muy útil de la transformada de Laplace tiene que ver con trasla-
ciones, es decir, si a es una número real
Z ∞
ax
L[e f (x)] = e−px eax f (x)dx
Z0 ∞
= e−(p−a)x f (x)dx
0
= F (p − a).

A veces es útil, usar el sı́mbolismo



ax
L[e f (x)] = L[f (x)] = F (p) .

p→p−a p→p−a

El resultado anterior es conocido como el primer teorema de traslación.


Ejemplo 14. Utilice el primer teorema de traslación hacer las siguientes evalua-
ciones:

a) L[e2x x3 ] b) L[e−x senh 5x] c) L[e3x cos 2x]

Para determinar la forma inversa del primer teorema de traslación, observe


que de F (p − a) primero se debe reconocer a F (p), para poder determinar f (x)
(f (x) = L−1 [F (p)]), y a continuación se multiplica f (x) por la función exponencial
eax , esto es,
L−1 [F (p − a)] = L−1 [F (p) s→s−a ] = eax f (x).

Ejemplo 15. Determinar las siguientes transformaciones inversas:


8
   
−1 2p + 5 −1 2p + 5
a) L b) L
(p − 3)2 2
p + 6p + 34

Ejemplo 16. Resolver el siguiente problema de valor inicial mediante la trans-


formada de Laplace

y 00 − 4y 0 + 4y = x3 e2x , y(0) = 0, y 0 (0) = 0

Traslación en el eje x
En ingenierı́a se presentan con frecuencia funciones que pueden estar “encendidas”
o “apagadas”. Por ello conviene definir una función especial que nos permita
trabajar en esos casos, dicha función se llama función escalón unitario y se define
como 
0 si 0 ≤ x < a,
U(x − a) =
1 si x ≥ a
Podemos entender la función U(x − a), que se mantiene apagada (vale 0) hasta
cierto momento x = a donde se enciende o activa (vale 1).
Algunas propiedades importantes de la función escalón unitario son: cuando
se multiplica por otra función definida para x ≥ 0, la función escalón unitario
“apaga” una parte de la gráfica de esa función. La función escalón unitario
también se puede usar para expresar en forma compacta funciones definidas a
pedazos.
Por ejemplo, si 
g(x) si 0 ≤ x < a,
f (x) =
h(x) si x ≥ a
equivale a
f (x) = g(x) − g(x)U(x − a) + h(x)U(x − a).

Ejemplo 17. Si 
 2 si 0 ≤ x < 2,
f (x) = −1 si 2 ≤ x < 3,
1 si x ≥ 3

equivale a
f (x) = 2 − 3U(x − 2) + 2U(x − 3).

Considere una función f (x) definida para x ≥ 0. La función definida por


tramos 
0 si 0 ≤ x < a,
f (x − a)U(x − a) =
f (x − a) si x ≥ a
básicamente es una traslación de la gráfica de f (x). Pues f (x−a)U(x−a) coincide
con la gráfica de f (x − a) para x ≥ a (que es la gráfica de f (x), x ≥ 0, desplazada
a unidades hacia la derecha sobre el eje x).
9

Teorema. Si a > 0 y f (x) es continua para x ≥ 0 y de orden exponencial,


entonces
L[U(x − a)f (x − a)] = e−ap L[f (x)] = e−ap [F (p)].
Al resultado anterior se le conoce como segundo teorema de traslación.
El teorema anterior en términos de la transformada inversa de Laplace, nos
dice que
L−1 [e−ap F (p)] = U(x − a)f (x − a).
Ejemplo 18. Use la transformada de Laplace para resolver el problema de valor
inicial dado:
a) y 00 + 4y = sen x u(x − 2π), y(0) = 1, y 0 (0) = 0.

 0 si 0 ≤ x < π,
b) y 00 + y = f (x), 0
y(0) = 0, y (0) = 1, donde f (x) = 1 si π ≤ x < 2π,
0 si x ≥ 2π.

Multiplicación de una función por xn


R∞
Supongamos que F (p) = L[f (t)] = 0
e−px f (x)dx, y que es posible intercambiar
la derivada y la integral, entonces
Z ∞
dF d
(p) = e−px f (x)dx
dp dp 0
Z ∞
∂ −px 
= e f (x) dx
∂p
Z0 ∞
= −xe−px f (x)dx
0
Z ∞
= − e−px (xf (x))dx
0
= −L [xf (x)] (p).
Es decir,
dF
L[xf (x)] = − (p). (19)
dp
Usando esta ecuación tenemos:
d
L[x2 f (x)] = L[x · xf (x)] = − L[xf (x)]
dp
 
d d
= − − L[f (x)]
dp dp
2
d
= L[f (x)]
dp2
d2 F
= (p)
dp2
En general, tenemos el siguiente teorema:
10

Teorema. Si F (p) = L[f (x)] y n entero positivo, entonces

dn
L[xn f (x)] = (−1)n F (p)
dpn

Ejemplo 20. Determinar L[x sen kx]

Ejemplo 21. Resolver el problema de valor inicial:


00
y + 16y = cos 4x, y(0) = 0, y 0 (0) = 1.

Convolución
Sean f y g funciones continuas a pedazos para x ≥ 0 y de orden exponencial.
Definimos el producto convolutivo o convolución de las funciones f y g, como
Z x
(f ∗ g)(x) = f (x − t)g(t)dt.
0

Si hacemos el cambio de variable u = x − t, en la definición de convolución, se


puede probar que (f ∗ g) = (g ∗ f ), es decir, la convolución de funciones es una
operación conmutativa.
La relación que existe entre la transforamada de Laplace y la convolución de
funciones, se muestra en el siguiente teorema.

Teorema. Si f y g son funciones continuas a pedazos para x ≥ 0 y de orden


exponencial, entonces

L[(f ∗ g)(x)](p) = L[f (x)](p)L[g(x)](p) = F (p)G(p).

Una manera equivalente de enunciar el teorema anterior es


Z x
−1
L [F (p)G(p)] = (f ∗ g)(x) = f (x − t)g(t)dt.
0

Nuevamente, si g(x) = 1 en el teorema anterior, la transformada de Laplace


de la integral de la función f se expresa como
Z x 
L[(f ∗ g)(x)](p) = L f (t)dt
0
F (p)
= .
p

Decimos que una función f es periódica con perı́odo T > 0, si para cada
x se cumple f (x + T ) = f (x). El siguiente teorema nos dice como obtener la
transformada de Laplace de una función periódica.
11

Teorema. Sea f una función continua a pedazos para x ≥ 0 y de orden exponen-


cial. Si f es periódica con perı́odo T , entonces
Z T
1
L[f (x)](p) = e−px f (x)dx.
1 − e−pT 0

Ejemplo 22. Resolver las siguientes ecuaciones diferenciales:


Rx
1. y 0 = 1 − sen x − 0 y(τ )dτ, y(0) = 0.
dy Rx
2. + 6y + 9 0 y(τ )dτ = 1, y(0) = 0.
dx

Delta de Dirac
Cuando se estudian sistemas mecánicos, eléctricos, etc., aparecen fuerzas externas
muy grandes en intervalos muy pequeños de tiempo , por ejemplo un golpe de
martillo en un sistema mecánico, un relámpago en un sistema eléctrico. La forma
de representar esta fuerza externa exterior es con la “función” delta de Dirac.
Si a y x0 son constantes positivas y x0 ≥ a, definimos

 0 0 ≤ x < x0 − a
1
δa (x − x0 ) = x0 − a ≤ x < x 0 + a
 2a
0 x ≥ x0 + a.

Note que Z ∞
δa (x − x0 ) dx = 1.
0

Debido a la propiedad anterior a δa (x − x0 ) se le conoce como la función impulso


unitario.
En la práctica conviene trabajar con otro tipo de impulso unitario, con una
“función” que aproxima a δa (x − x0 ), definida como

δ (x − x0 ) = lı́m δa (x − x0 ) .
a→0

A la expresión anterior, que no es una función, se le conoce como la “función”


delta de Dirac. Algunas de las propiedades más importantes de la “función” delta
de Dirac son

1. δ (x − x0 ) es infinita en x = x0 y cero si x 6= x0 .
R∞
2. 0 δ (x − x0 ) = 1.

e − e−ap
 ap 
−x0 p
3. L [δa (x − x0 )] = e .
2ap
4. L [δ (x − x0 )] = lı́m L [δa (x − x0 )] = e−x0 p .
a→0
12

5. Si x0 = 0, entonces L [δ (x)] = 1.
R∞
6. Si f es continua, 0 f (x)δ (x − x0 ) dt = f (x0 ).

7. L [f (x)δ (x − x0 )] = e−x0 p f (x0 )

Ejemplo 23. Resolver la ecuación diferencial y 00 + 2y 0 + 2y = δ(x − π), y(0) =


0, y 0 (0) = 0.

Concluimos este capı́tulo presentando una tabla con las propiedades princi-
pales de la transformada de Laplace.

Propiedades de la transformada de Laplace

L[αf (x) + βg(x)](p) = αF (p) + βG(p)


L[eax f (x)](p) = F (p − a)
L[U(x − a)f (x − a)](p) = e−ap F (p)
L[f 0 (x)](p) = pF (p) − f (0)
L[f 00 (x)](p) = p2 F (p) − pf (0) − f 0 (0)
L[f (n) (x)](p) = pn F (p) − p(n−1) f (0) − p(n−2) f 0 (0) − ... − f (n−1) (0)
Rx
L[ 0 f (t)dt](p) = F (p)
p
L[−xf (x)](p) = F 0 (p)
L[(−1)n xn f (x)](p) = F (n) (p)
R∞
L[ f (x)
x
](p) = p F (p)dp
Rx
L[ 0 f (x − t)g(t)dt](p) = F (p)G(p)
L [δ (x − x0 )] = e−x0 p
L [δ (x)] = 1
L [f (x)δ (x − x0 )] = e−x0 p f (x0 )

También podría gustarte