Story Transcript
CAPÍTULO
6 La transformada de Laplace
6.2 Definición de la transformada de Laplace 6.2.1
Definición y primeras observaciones
En la gran mayoría de los sistemas de interés para la física y la ingeniería es posible (al menos en principio) predecir su comportamiento futuro partiendo de condiciones dadas en un determinado tiempo, el cual podemos desde luego suponer que es t D 0. Sólo en muy contados ejemplos es factible predecir el comportamiento pasado del sistema. En lo que sigue nos ocuparemos solamente de la parte de las funciones f .t/ definida para valores t 0, sin darle importancia a lo que sucede para t < 0. Con esta aclaración podemos enunciar la siguiente definición de la transformada de Laplace, en ocasiones denominada unilateral: La transformada de Laplace (TL) de una función f .t/ se define mediante:
Lf f .t/g D F .s/ D
Z
1
e
st
f .t/ dt;
0
para todos los valores s 2 R para los cuales la integral impropia anterior sea convergente. Si Lf f .t/g D F .s/, diremos también que f .t/ es la transformada inversa de Laplace de F .s/, lo que denotaremos como: 1 f .t/ D L fF .s/g: A esta conexión entre f .t/ y F .s/ se le suele representar mediante el esquema: f .t/
1. canek.azc.uam.mx: 24/ 9/ 2010
1
! F .s/:
2
Ecuaciones diferenciales ordinarias
Observaciones: 1. Recordemos que una integral impropia como la anterior se define como un límite de integrales sobre intervalos finitos, es decir, Z Z R
1
e
0
st
f .t/ dt D lím
e
R!1
st
f .t/ dt;
0
si el límite existe. El resultado de la integral será en general una función que depende de s. Al calcular una TL no siempre se especifica el rango de valores s para los cuales la integral existe; no obstante en algunos casos es recomendable determinarlo.
2. En ocasiones a esta transformada se le llama unilateral. Podemos ampliar la definición dada la TL bilateral. Ésta es Z 1 Lf f .t/g D F .s/ D e st f .t/ dt: 1
No obstante, el manejo y aplicación de ésta tiene otras complicaciones debidas en parte a problemas de convergencia que caen fuera del alcance e interés de este libro.
6.2.2
Cálculo de la TL
Usaremos la definición para calcular la TL de algunas funciones. Ejemplo 6.2.1 Calcular la TL de la función escalón unitario ( 0; si t < 0I u.t/ D 1; si t 0: H
Usando la definición de la TL: Z 1 U.s/ D e
st
D
1 s
u.t/dt D lím
R!1
0
lím .e
sR
R!1
R
Z
e
st
1/ D
dt D lím
R!1
0
1 s
lím
R!1
"
1 e sR
# 1 st R e D s 0 1 :
Observamos ahora que el límite anterior existe, sólo si s > 0. En este caso, lím e
sR
R!1
consecuencia, U.s/ D
1 s
. 1/ D
D lím
R!1
1 D 0, en e sR
1 : s
De esta manera, hemos hallado nuestra primera fórmula de TL: u.t/
1 ; con s > 0: s
!
(6.1)
Ejemplo 6.2.2 Calcular la TL de la función f .t/ D e at , para t 0 & a 2 R. H
Aquí 1
F .s/ D
Z
D
1
e
st
0
a
s
e
at
dt D lím
R!1
h
lím e .a
R!1
s/R
Z
R
e
0
i
1 :
.a s/t
dt D lím
R!1
"
1 a
s
e .a
R# s/t D 0
(6.2)
6.2 Definición de la transformada de Laplace
3
Al estudiar el límite anterior, concluimos su existencia siempre y cuando a que lím e .a s/R D 0, por lo que:
s < 0. En este caso, tenemos
R!1
1
F .s/ D
a
s
1
. 1/ D
s
; con s > a:
a
También podemos escribir: e at
!
1 s
a
; con s > a:
(6.3)
Ejemplo 6.2.3 Calcular la TL de la función f .t/ D H
Partiendo de la definición de la TL: Z 1 Z 1 st Lf f .t/g D e f .t/ dt D e 0 .2 s/t
(
e 2t ; 4;
st 2t
e dt C
0
si t < 1I si 1 t:
Z
1
e
st
1
4 dt D
1 sR e e2 s e s C 4 lím e D D R!1 2 s 0 2 s s s 2 s s 1 e 4e D C , para s > 0 y s ¤ 2: s 2 s
Para s D 2, se reduce la integral: Z 1 F .2/ D e
2t
0
D Entonces:
Z
0
f .t/ dt D
1
dt C 4
Z
1
e
Z
2t
1
1
1
e
2t 2t
0
dt D 1
e
dt C
1 2
Z
lím
R!1
4e s e2 s C ; s 2 s Lf f .t/g D F .s/ D 1 C 2e 2;
s
Z
1
e
.2 s/t
0
C
4e s
dt C 4
s
Z
1
e
st
1
dt D
D
1
e
2t
4 dt D R! 2t 2e D 1 C 2e 1
2
:
1
si s > 0 y s ¤ 2I si s D 2:
Propiedad de linealidad de la TL
Una propiedad importante de la TL es la siguiente: Z 1 Z st Lf af .t/ ˙ bg.t/g D Œaf .t/ ˙ bg.t/e dt D a 0
1
f .t/e
0
st
dt ˙ b
D aLf f .t/g ˙ b Lf g.t/g D aF .s/ ˙ bG.s/:
Z
0
1
g.t/e
st
dt D
Es decir:
Lf af .t/ ˙ bg.t/g D aLf f .t/g ˙ Lf g.t/g: Una fórmula recursiva para la TL
(6.4)
Hay fórmulas de integración que reducen el resultado de una integral a una integral más sencilla, por lo común mediante una integración por partes. Vamos a calcular Lf t n g (para ello usaremos en el desarrollo):
4
Ecuaciones diferenciales ordinarias
u D tn st
dv D e n
Lf t g D
Z
1
e
st n
t dt D
0
D
n 1 R lím s R!1 e Rs
Z
) du D nt n dt
n
st
t e
dtD lím
R!1
0
n 0 C s
"
s t ne s
st
e
st n 1
t
: R# 0
1
Z
dtI
st
e
) vD
1
1
dt:
Z
0
1
st
e
s
nt n
1
dt D
0
Rn D 0. Por lo tanto: R!1 e sR
Si aplicamos n-veces la regla de L0 Hôpital, entonces, lím
Lf t n g D
n s
Z
1
e
st n 1
t
dt;
si s > 0:
0
Registramos lo anterior como una fórmula recursiva: n ˚ Lf t n g D L t n 1 ; para n D 1; 2; 3; s
(6.5)
Aplicando esta fórmula podemos obtener los siguientes resultados: 1 ˚ s ˚ 2 2 ˚ 1 L t D L t D s ˚ 3 3 ˚ 2 L t D L t D s ˚ 4 4 ˚ 3 L t D L t D s Generalizando, vemos lo siguiente:
1 s 2 s 3 s 4 s
1 1 1 D 2: s s s 1 21 2 D 3 : s s 21 321 3 D : s s4 321 4321 D : 4 s s5
L t 1 D L t 0 D Lf 1g D ˚
Lf t n g D
nŠ ; s nC1
para
Ejemplo 6.2.4 Calcular la TL de la funciones f1 .t/ D t 3 C 3t, f2 .t/ D 2t H
Usando las propiedades anteriores:
(6.6)
n D 1; 2; 3; 5 & f3 .t/ D 2e t C 3e 2t
e
3t
.
3Š 1Š 6 3 6 C 3s 2 C 3 D C D : s4 s2 s4 s2 s4 1 1 2 5s Lf 2t 5g D 2Lf t g 5Lf u.t/g D 2 2 5 D 2 : s˚ s s ˚ ˚ ˚ L 2e t C 3e 2t e 3t D 2L e t C 3L e 2t L e 3t D 3 1 2.s 2/.s C 3/ C 3.s 1/.s C 3/ .s 1/.s 2/ 2 C D D D s 1 s 2 sC3 .s 1/.s 2/.s C 3/ 2.s 2 C s 6/ C 3.s 2 C 2s 3/ .s 2 3s C 2/ 4s 2 C 11s 23 D D : .s 1/.s 2/.s C 3/ .s 1/.s 2/.s C 3/
L t 3 C 3t D L t 3 C 3Lf t g D ˚
˚
De manera similar es posible calcular la TL de cualquier función que sea combinación lineal de potencias enteras de t y exponenciales e at , para a 2 R. Hasta el momento, a partir de la definición y primeras propiedades de la TL tenemos una cantidad limitada de funciones cuyas TL podemos calcular. Conforme avancemos en las siguientes secciones se ampliará considerablemente la clase de funciones para las cuales podemos calcular su TL.
6.2 Definición de la transformada de Laplace
6.2.3
5
Más ejemplos de cálculo de la TL
Al estudiar ED lineales homogéneas con coeficientes constantes, consideramos conveniente introducir los números complejos de la forma z D ˛ C iˇ y las correspondientes fórmulas de Euler como sigue: e i D cos C i sen
&
i
e
D cos
i sen :
Las propiedades de la integral son válidas para constantes complejas, es decir, que si z D ˛ C iˇ, entonces: Z Z Z Z zf .t/ dt D .˛ C iˇ/f .t/ dt D ˛ f .t/ dt C iˇ f .t/ dt D Z Z D .˛ C iˇ/ f .t/ dt D z f .t/ dt: Por otro lado, la fórmula
L e zt D ˚
1 z
a
es válida con z complejo. Ya que, si z D ˛ C iˇ, en la deducción de la fórmula (véase ejemplo 6.2.2), requeriríamos considerar: lím e .z
s/R
D lím e iˇR e .˛
R!1
s/R
R!1
D0
cuando ˛
s < 0:
Así, para que la fórmula anterior sea válida es preciso que su parte real cumpla Re.z/ D ˛ < s. Utilicemos ahora las consideraciones previas para resolver los ejemplos siguientes: Ejemplo 6.2.5 Calcular la TL de las funciones f .t/ D cos kt
&
f .t/ D sen kt
para k 2 R:
H Si usamos la observación previa sobre linealidad de la integral con números complejos y la fórmula de Euler: Z 1 Z 1 ˚ ikt st i k t L e D e e dt D e st .cos kt C i sen kt/ dt D 0 0 Z 1 Z 1 st D e cos kt dt C i e st sen kt dt D Lf cos ktg C i Lf sen ktg: 0
0
También, de acuerdo a lo dicho antes, sobre la transformada de una función con exponente complejo, ˚
L e
ikt
D
1 s
ik
D
1 s
ik
s C ik s C ik
D
s C ik s k D 2 Ci 2 : 2 2 2 s Ck s Ck s C k2
Por lo tanto:
Lf cos ktg C i Lf sen kt g D
s k Ci 2 : s2 C k2 s C k2
De aquí, al igualar las partes reales e imaginarias respectivamente, determinamos que cos kt sen kt
s ; C k2 k ! 2 ; s C k2 !
s2
con s > 0; con s > 0:
La validez de estas fórmulas se logra puesto que Re.i k/ D 0 < s.
6
Ecuaciones diferenciales ordinarias Se puede generalizar la fórmula (6.5) para Lf t n g al caso de exponentes no enteros, introduciendo la función Gamma .z/, estudiada por Euler en el siglo XVIII: .z/ D
Z
1 x z 1
e
x
dx;
0
que se puede definir para números complejos z, con Re.z/ > 0, aunque aquí restringimos nuestra atención a .r /, con r > 0. La propiedad que caracteriza a .r / es la de generalizar el factorial de enteros positivos, ya que integrando por partes: .r C 1/ D
Z
R! x
1
e
x r
x dxD lím
e
R!1
0
Z
x r
0
1
e
x
rx
r 1
0
dx D r
Z
1
e
x r 1
x
0
dx D r.r /:
Hemos usado para resolver por partes la primera integral:
˚
u D xr
dv D e
x
) du D r x r
dx
1 x
) vD e
dxI
:
De este modo, la fórmula recursiva obtenida (6.7)
.r C 1/ D r.r / junto con la observación de que .1/ D 1 nos permiten concluir que .2/ D .1 C 1/ D 1.1/ D 1I .3/ D 2.2/ D 2 1 D 2ŠI .4/ D 3.3/ D 3 2 1 D 3ŠI y en general
(6.8)
.n C 1/ D nŠ n D 1; 2; 3; Ejemplo 6.2.6 Calcular Lf t r g, con r > 0. H
Aplicando la definición de TL:
Lf t r g D
Z
0
1
e
st r
t dtD
Z
0
1
e
x
x r s
dx s
D
1 s r C1
Z
0
1
e
x r
x dx D
Donde hemos usado para resolver la primera integral el cambio de variable x D st.
.r C 1/ : s r C1
De aquí que para la función f .t/ D t r , con r > 0:
Lf t r g D
.r C 1/ : s r C1
Puesto que, usando (6.8)
Lf t n g D
.n C 1/ nŠ D nC1 : s nC1 s
Vemos que la fórmula (6.9) generaliza la anterior (6.6).
(6.9)
6.2 Definición de la transformada de Laplace
7
p 1 D . 2 Usando la definición de : Z 1 1 D e 2 0
Vamos a demostrar H
Haciendo uD
se tiene
x
1
x2
1
dx D
Z
1
1
0
dx:
x2
p x ) x D u2 ) dx D 2u du;
Z 1 u2 Z 1 e 1 2u du D 2e D u 2 0 0
Ahora
x
e
u2
du:
Z 1 Z 1 Z 1Z 1 1 u2 v2 D 2e du 2e dv D 4 e 2 0 0 0 0 Usando coordenadas polares (u D r cos I v D r sen ): 2
Z Z 1 2 1 D4 e 2 0 0 2
r2
r drd D 4
de donde
Z
2
0
1 e 2
r2
1 Z d D 4
0
0
p 1 D : 2
2
.u2 Cv 2 /
dudv:
2 1 d D 2 D ; 0C 2 0
A partir del resultado anterior, usando (6.7): 1 1 3 1 D D C1 D 2 2 2 2 3 3 5 3 D D C1 D 2 2 2 2 :: : n o 3 Ejemplo 6.2.7 Calcular L t 2 . H
1p : 2 3p : 4
Usamos (6.9): n
L t
3 2
o
D
5 2 5
s2
3p r 3 4 p D D : 4 s5 s5
1 p . t
1
Ejemplo 6.2.8 Calcular L H
L p
t
n DL t
1 2
o
D
s
1 C1 2 1 C1 2
D
1 2 1
s2
r p D p D : s s
Podemos aplicar, entonces, la fórmula para los exponentes r negativos que cumplan 1 < r < 0. 1 Por supuesto, para exponentes r que no sean enteros o múltiplos impares de , el cálculo de .r C 1/ puede 2 ser muy complicado y requerir de alguna aproximación numérica.
8
Ecuaciones diferenciales ordinarias
Ejercicios 6.2.1 Definición. Soluciones en la página 9 Utilizar la definición para hallar la transformada de Laplace de cada una de las siguientes funciones: 1
1. f .t/ D e 5 t . 2. f .t/ D e t 3. f .t/ D 6
2
.
t 2.
4. f .t/ D t 8 C e t . ( 3t; si t < 1I 5. f .t/ D 0; si t 1: ( t 2 C 3t 2; 6. f .t/ D 0; 7. f .t/ D t cos at; 8. cosh kt D
ek t C e 2
si 1 t 2I si t … Œ1; 2 :
a constante. kt
&
senh kt D
ek t
e 2
kt
.
6.2 Definición de la transformada de Laplace
9
Ejercicios 6.2.1 Definición. Página 8 1.
5 5s
1
;
para s >
2
e ; para s > 1. s 1 6 2 3. ; para s > 0. s s3 2.
4.
7s 2 C 9s 1 . s 2 .s 1/
1 . 5
5. 3
„
1 e s
s
1 e s2
s
C
« 1 . s2
e s e 2s .s 2/ C 3 .s C 2/. 3 s s s 2 a2 7. . .s 2 C a2 /2 s k 8. Lf cosh.k t /g D 2 ; Lf senh.k t /g D 2 : 2 s k s k2 6.