Transformadas de Laplace
Transformadas de Laplace
Transformadas de Laplace
TRANSFORMACIN DE LAPLACE
5.1 Introduccin
El clculo operacional y las transformadas integrales.
Las transformadas integrales lineales de Laplace, Fourier, Hankel, etc., presentan
una propiedad operacional bsica que proporciona una imagen algebraica de las expresiones
diferenciales. Aplicando las transformadas integrales a problemas con valores iniciales y
problemas con valores en la frontera, es posible obtener ecuaciones imagen, que en general
son mas simples de resolver que las originales. Si la transformacin inversa de la solucin
imagen es posible se obtiene la solucin buscada.
Los precursores del clculo operacional y las transformadas integrales fueron:
Laplace (1749-1827)
Cauchy (1789-1875)
Heaviside (1850-1925).
La transformacin T{F(t)} es lineal si para cada par de funciones F(t) y G(t) existen dos
constantes diferentes de cero C1 y C2, tal que
T{C1F(t) + C2G(t)} = T{C1F(t)} + T{C2G(t)} = C1T{F(t)} + C2T{G(t)},
(5.1)
T F t K t , s F t dt f s ,
b
(5.2)
(5.3)
e st
;
L 1 e 1 dt
0
s 0
por lo tanto, cuando s > 0, el valor de la integral converge y es
1
L 1 , s > 0.
s
b) Si, F(t) = e k t cuando t > 0 y en donde k es una constante
st
1 ( s k )t
L e e e dt
e
0
k s
por lo tanto, para convergencia es necesario que s > k y
1
L ek t
, s > k.
sk
kt
st k t
;
0
1
,
s2
L t 2
1
,
s3
L senkt
k
,
s k2
2
L cosh kt
s
s k2
2
cuando s > 0, para los dos primeros casos y s > | k | para el resto. Aunque hay otras formas
ms simples de obtener la imagen de funciones objeto, por ejemplo, utilizando la propiedad
de linealidad de las transformadas de Laplace pude ser til para obtener la transformada de
algunas funciones. Ejemplo para encontrar la transformada de Laplace de cos kt con ayuda
de la identidad de Euler
1
1 1
2s
1
1 1
L cos kt L eik t e i k t
;
2
2 2
2
2
2 s ik 2 s ik 2 s i k
s k .
Hay que notar la necesidad de definir el dominio de s, para que la integral de transformacin
de Laplace, la cual es una integral impropia, sea convergente.
5.3 Funciones Seccionalmente Continuas. Orden exponencial
Una funcin F(t) es seccionalmente continua sobre un intervalo limitado a < t < b, si
tal intervalo puede ser dividido en un numero finito de subintervalos, al interior de cada uno
de los cuales F es continua y tiene limite finito conforme t se aproxime desde el interior a
alguno de los puntos finales del subintervalo. La transformada de cada funcin de esta clase
sobre el intervalo (a,b), existe; y es la suma de las integrales de las funciones continuas
sobre los subintervalos.
La funcin escaln unitario
S k t 0 cuando 0 t k ,
1 cuando t k
es una funcin seccionalmente continua, que frecuentemente aparece en modelado de
procesos. La transformada de Laplace de esta funcin es
1
L S k t e st Sk t dt e st dt e st ;
k
0
k
s
e k s
.
s
Una funcin F(t) es de orden exponencial como t tiende a infinito, si alguna constante
existe tal que el producto
e t F t
este limitado para toda t mas grande que algn numero finito T. Por lo tanto F t no crece
mas rpido que Me t , donde M es alguna constante. Lo anterior es expresado diciendo que
F(t) es de orden exponencial, o que F(t) es de O e t . Si F(t) es seccionalmente continua y
de orden exponencial se asegura la convergencia absoluta de la integral de Laplace cuando
s (condicin suficiente mas no necesaria).
Una funcin que no cumple lo anterior es F t t
transformada existe y se obtiene como sigue:
1
1
2
L t 2 e st t 2 dt
0
s
1
2
2
,
s 2
e x dx
por lo tanto
12
L t
s
s 0.
Par la primera derivada de F(t), con la definicin de transformadas de Laplace y por medio
del mtodo de integracin por partes, formalmente se tiene que
L F ' t e st F ' t dt e st F t s e st F t dt
en donde f s L F t .
L F ' t sf s F 0 ,
(5.4)
(5.5)
L F (t ) s sL F t F 0 F ' 0 .
con lo cual se tiene la transformacin de una funcin de segundo orden
L F " t s 2 L t sF 0 F ' 0 ,
(5.6)
L F " t s 2 f ( s ) sF 0 F ' 0 .
(5.7)
en consecuencia
L sin kt
k
.
s k2
2
L F m 1 t 0 s m 1L t m m !
y por lo tanto
m!
s 0.
, m 1, 2,....
(5.9)
s m 1
Si el exponente no es entero, se puede recurrir a la definicin de la transformada de Laplace
L t m
k 1
L t k e st t k dt
x k e x dx
s0
(5.10)
0
s
La integral de la derecha representa la funcin gamma, o funcin factorial, con un
argumento k+1. Por lo tanto
k 1
(5.11)
L tk
k 1, s 0 .
s k 1
Notar que si k = m, se pueden igualar (5.9) y (5.11)
m 1 m !
m 1
s m 1
s
De donde
m 1 m !
Ejemplo. Clculo de L
Sea la funcin
G t F d
t
continua, y G(0) = 0. Tambin G(t) = F(t), excepto para aquellos valores de t para los
cuales F(t) es discontinua; por lo tanto G(t) es seccionalmente continua sobre cada intervalo
finito. Si la funcin G(t) es tambin de O(et), aplicando el teorema 1,
L G ' t sL G t L F t
Si se tiene una funcin f(s), que es la imagen de F(t), entonces existe un operador que
permite encontrar F(t) a partir de f(s), este operador es L1 define la transformada inversa de
Laplace, que es L1 f s F t . Como L es un operador lineal, su correspondiente
operador inverso L1 , cumple con la condicin
L1 Af s Bg s AL1 f s BL1 g s
AL1 f s BL1 g s AF t BG t
(5.12)
es decir, L1 tambin es una transformacin lineal y existe una relacin univoca entre ambos
operadores.
5.6 Teorema de sustitucin
f s a e s a t F t dt e st e at F t dt ,
s a.
(5.13)
(a)
sY 0 Y ' 0 2kY 0
s k 2
o
1
1
Y ' (0) kY (0)
.
(b)
s k
s k 2
En donde la ecuacin (b) es un arreglo ideal para encontrar la transformada inversa.
Aplicando L1 a la anterior se obtiene
Y(t) = (C1t + C2)et
y donde C1 = Y(0) k Y(0) y C2 = Y(0).
y s Y 0
s s
s 1
s 2 s 2 1
s 2 s 2 1
Luego
1 sY (0)
s2 s2 1
La transformada inversa de la anterior da
Y(t) = t + Y(0)cos t
y ( s)
5.8 La Convolucin
f s g s e s F d e svG v dv
s v
F G v d dv
f s g s F e
s v
G v dvd
Para la integral interna, sea t una nueva variable, tal que v t y dv dt , luego, este
producto es
f s g s F e stG t dtd
e st F u G t dud e st F G t d dt
F t * G t F G t d G F t d
t
La temperatura inicial de una placa de material homogneo limitado por los planos x = 0 y x
= l es T0, Se desea determinar el perfil de temperaturas en el slido despus de aislar la
cara x = 0 y de colocar la cara x = l a temperatura cero.
Este problema se describe matemticamente por medio de las siguientes ecuaciones
T x ,t
2T x ,t
, 0 x l, t 0 ,
a2
t
x2
(1)
T x ,0 T0 , 0 x l ,
(2)
T 0,t
0, t 0 ,
x
(3)
T l ,t 0 , t 0 .
(4)
e st
T x ,t
2T x ,t
dt e s t a 2
dt , 0 x l
0
t
x2
(5)
dx 2
dt
luego,
d 2 T x , s
s T x , s T0 a
, 0 xl
dx 2
(6)
d T 0 , s
0
dx
(7)
T l , s 0
(8)
T x , s C1e
s
x
a
C2 e
s
x
a
T0
s
(9)
s
s
C2
a
a
C1 C 2
T
s
x 0
a
s
Aplicando la condicin (8) en la ecuacin (10)
T
1
C 0
s
s
cosh
l
a
Luego, la solucin imagen es
T x , s C cosh
(10)
s
x
1 cosh
a
, 0 x l ,
T x , s T0
s
s
s cosh
l
(11)
De tablas
1
L1 1
s
y
n
cosh s x
2n 1x e
4 1
L
cos
1
2a
n 1 2n 1
s cosh s a
2 n 12 2t
4a2
2n 1
n 1
2a
(12)
(1)
U x ,0 0 , x 0 ,
(2)
U 0,t F0 , t 0 ,
(3)
lim U x ,t 0 .
(4)
u 0 , s
F0
,
s
(5)
(6)
lim U x , s 0 ,
(7)
u x , s C1e
sh
x
a
C2 e
sh
x
a
(8)
u x , s C1e
sh
x
a
(9)
F0
s
y la solucin imagen es
ux, s
sh
x
a
F0
e
s
(10)
4 a 2t
xe ht
U x ,t F0 1*
e
2a t 3
en donde
F0 x
2a
2 F0
x
2a t
x2
4 a 2
d
2 3
2 h x2
4 a 2 2
x
.
2a t
(1)
C
2C
D2
, 0 x , t 0
t
x2
(2)
C x ,0 0 , x 0
(3)
C x ,0 C 0 , 0 x
(4)
lim C x ,t 0 , t 0
(5)
lim C x ,t C 0 , t 0
(6)
C 0 ,t C 0 ,t , t 0
(7)
C 0 ,t
C 0 ,t
D1
D2
, t0
x
x
(8)
d 2 c x , s
, x0
dx 2
s c x , s C 0 D2
d 2 cx , s
, x0
dx 2
(9)
(10)
lim c x , s 0
(11)
lim c x , s
C0
s
(12)
c 0 ,s c 0 ,s
D1
(13)
dc 0 , s
dc 0 , s
D2
dx
dx
(14)
s
x
D1
c x , s De
s
x
D2
Be
s
x
D1
Ee
, x0
s
x
D2
C0
, x0
s
(15)
(16)
Aplicando las condiciones (11) y (12) a (15) y (16) respectivamente, se determina que
A 0 y E 0 , por lo cual se tiene que
c x , s Be
s
x
D1
, x0
(17)
y
c x , s De
s
x
D2
C0
, x0
s
(18)
C0
B
s
de donde
D B
C0
s
(19)
D1
s
s
B D2
D
D1
D2
o
B
D2
D
D
D1
(20)
D2
.
D1
Los valores de B y D se obtienen de (19) y (20)
en donde
1 C0
s
C0
D
s
Por lo tanto las soluciones imagen son:
c x , s
s
x
D1
C0 e
s
, x0
(21)
y
s
1 C e D2
0
c x , s C 0
s s
C x ,t
x
C0
erfc
2 Dt
, x 0
(22)
(23)
, x 0
C x ,t C 0 1
erfc
2 D t
2
C x ,t
(24)
Las ecuaciones (23) y (24) determinan las concentraciones de humedad dentro del poro.
Teorema 5.4 Si f(s) y g(s) son las transformadas de dos funciones F(t) y G(t), que son
seccionalmente continuas sobre cada intervalo 0 t T y de O e t , luego la transformada
de la convolucin F(t)*G(t) existe cuando s ; y es f(s)g(s). Por lo tanto la transformada
inversa de f(s)g(s) es dada por la formula
-1{f(s)g(s)} = F(t)*G(t),
(5.10)
(5.11)
Para probar el teorema sean h(s), g(s) y f(s) las funciones imagen de las funciones objeto
H(t), G(t) y F(t) y sea
h(s) = g(s)f(s)
Luego, a partir de (5.11)
h s L H t e st F G t ddt
f t
F G t u t
0
para t
para t
h s L H t e st F t G u t d dt
hs
e
e
st
st
Sea
hs
F t u t dtG d
F t dtG d
t
s
F d G d
e s e s F d G d
f s e s G d
0
f s e s G d
0
h s f s g s LQQD
Ejemplo:
Se desea encontrar la respuesta transitoria para un reactor tanque continuo agitado
operando con una variacin de la composicin de entrada con el tiempo. Asuma
comportamiento isotrmico y velocidad de consumo lineal para la especie A de acuerdo a:
k
A productos
Se puede asumir volumen y velocidad de flujo constante. Encuentre la respuesta
cuando la entrada lleva una perturbacin (a) con un cambio escalonado (se eleva la
concentracin permanentemente) y (2) una perturbacin de impulso (momentneamente
vara la concentracin).
es conveniente reacomodar las ecuaciones a una forma estndar, la cual produzca una
constante de tiempo del sistema explcita. As, se desea rearreglar para que la ecuacin tenga
la forma
dy
y f t
dt
Esto puede ser hecho al rescribir la ecuacin
dc
V A kV q c A qc A 0
dt
Despus, al dividir entre kV q y definiendo la constantes de tiempo como
V
kV q
por lo tanto, se encuentra
dc
A c A c A 0
(3)
dt
donde aparece la fraccin adimensional q q kV .
Si se aplica la trasformada de Laplace en este momento, se necesitar incluir la
condicin inicial, debido a la aparicin de la derivada dc A dt ; as,
dc
A sC A c A 0
dt
Es posible eliminar tanto exceso al escribir las variables de estmulo c A0 y respuesta c A en
trminos de la desviacin del estado estacionario inicial, para escribir
c A t c A t c A
(4)
c A0 t c A0 t c A0
c A 0 0
(5)
q , C A 0 t
2c
A0
A 0
V
C A t
c A0
c
A 0
t
En general, cualquier fraccin positiva o negativa de c A0 puede ser usada para la
magnitud de la funcin escaln de la entrada. Es claro que la composicin completa de
entrada es
c A0 t u t c A0 c A0 ;
c A0 u t c A0
(7)
por lo tanto,
C A s
c A 0
s s 1
(8)
c A 0
1
C A s
s s 1
(9)
1
1
s s 1
la cual se puede
0 ue
t
e at e a
d
a
a t
1
1 e at
a
c A t c A0 1 exp t
(10)
(11)
c A t
1 exp t
cA
c A c A y por supuesto que c A 2c A .
(12)
(13)
Esto significa que la constante de tiempo puede ser deducida cuando la respuesta para c A t
es alrededor del 63% de la distancia al estado estacionario final. La siguiente figura ilustra la
respuesta al escaln.
CA
2C
C A t
CA
C
(14)
El factor de peso c0 cuenta cuanto soluto A fue inyectado, en relacin a el flujo molar de
entrada de qc A0 .
As, suponiendo que la velocidad de flujo del reactor tanque agitado continuamente
es, q = 0.1 litros/segundo y c A0 1 mol/litro. Ahora se supone que se inyectan 100 cc de
soluto A 1 mol/litro de composicin c A0 dentro de la corriente que fluye en la entrada.
Modelando esta inyeccin como un impulso, ya que fue hecho rpidamente, relativo a la
constante de tiempo del sistema . Entonces, el flujo molar neto de entrada es
qc A0 t q c A0 0.1litros c A0 t
donde es claro que t tiene unidades recprocas de tiempo, seg-1. Ahora se ve que
as,
0.1
t c A0
c A0
q
c A0
q
dc A
c A c0t
dt
(15)
(16)
c0
exp t
(17)
Excepto por el factor de peso c0 , esta pudo haber sido obtenida directamente mediante la
diferenciacin de la respuesta escalonada. La respuesta esta ilustrada en la siguiente figura
cA
c
c A t
cA
c
Ejemplo:
El reactor tanque agitado continuo del ejemplo anterior esta sujeto a un cambio
escalonado en el sistema del flujo, mientras se mantiene un volumen constante para el
reactor. Encuentre la respuesta de la composicin de salida, c A t , a tales perturbaciones.
Las variables a derivar en este caso sern
qt q q t
(1)
c A t c A c A t
(2)
y c A0 t c A0 se mantiene a un valor constante. Insertando estas en el balance de materia para
el componente A, produce
V
dc A
q c A qc A kV c A q t c A t q t c A0 q t c A q c A t kVc A t
dt
(3)
El primer grupo de trminos en parntesis son iguales a cero, ya que este grupo define la
condicin de estado estacionario original, la cual en una forma simplificada es la misma que
en el ejemplo anterior,
c A c A 0
q t c A t 0
dc A
q kV c A q t c A0 q t c A
dt
estacionario requiere que c A c A0 ,
(4)
dc A
1 c A0 qt
c A
dt
q kV
(5)
Suponiendo que la velocidad de flujo toma un cambio escalonado con una magnitud igual al
50% del flujo original; esta sera expresada como