Tema Laplace
Tema Laplace
Tema Laplace
2011-2012
Tema 7. La Transformada de Laplace
Depto. Teorı́a de la Señal y Comunicaciones e IT
Universidad de Valladolid
1 Introducción
En los temas anteriores hemos estudiado los aspectos básicos del análisis de Fourier. Esta
herramienta nos permite describir las señales de interés en términos de una superposición
(discreta o continua, dependiendo de la periodicidad o no de la señal) de exponenciales
complejas (en tiempo discreto o continuo, dependiendo de la naturaleza discreta o continua
de la señal). Esta descripción tiene dos ventajas:
1. Puesto que las exponenciales complejas son autofunciones de cualquier sistema LTI,
la caracterización de éstos se reduce al producto de la entrada con la respuesta en
frecuencia, operación conceptualmente mucho más sencilla que la convolución.
1
1 j ω+π
X(ω) = Π(ω)− log
2 2π ω−π
Puesto que x(t) en este caso es de energı́a finita, su transformada de Fourier es otra señal
que también tiene energı́a finita, aunque vemos que esta última presenta dos singularidades
en ω = ±π. Ahora bien, ¿qué ocurrirá si consideramos la señal t·x(t)? De las propiedades
d
de la transformada de Fourier deducimos que su transformada será j dω X(ω). Puesto que
X(ω) contiene funciones escalón, al derivar obtendremos deltas de Dirac, correspondientes
al hecho de que t·x(t) = π1 sin(πt) u(t) no tiene energı́a finita:
Z ∞ Z ∞ N Z n+1 N
2 1 2 1 X
2
X 1
|t·x(t)| dt = 2 sin (πt)dt = 2 lim sin (πt)dt = lim 2
= ∞.
∞ π 0 π N →∞ N →∞ 2π
n=0 n n=0
(4)
De hecho, la transformada de Fourier de t·x(t) puede calcularse de forma sencilla:
1 j 1/2π 1/2π
F {t·x(t)} (ω) = F sin(πt) u(t) (ω) = (δ(ω + π) − δ(ω − π)) + − . (5)
π 2 ω+π ω−π
Ahora tenemos dos polos en ±π (términos 1/(ω ± π) que ahora no son de cuadrado
integrable) correspondientes a las derivadas de la parte imaginaria de X(ω), y además
aparecen las deltas de Dirac correspondientes a las derivadas de las funciones escalón.
2
y por tanto no se puede operar con ellas como si lo fueran. En el presente tema intro-
duciremos un nuevo tipo de transformada para señales no periódicas en tiempo continuo,
la transformada de Laplace. Recordemos que la transformada de Fourier transformaba
la variable real t en la variable real ω. La transformada de Laplace transforma la variable
real t en la variable compleja s = σ + jω, siendo ésta la principal diferencia entre ambas.
Como veremos a lo largo del tema, esta nueva herramienta tiene dos ventajas sobre la
transformada de Fourier:
Con respecto al segundo punto, conviene recordar aquı́ que las exponenciales complejas
son autofunciones de cualquier sistema LTI. Es decir, ∀s ∈ C:
Z ∞ Z ∞
st st
e ∗h(t) = exp(s(t − τ ))h(τ )dτ = e exp(−sτ )h(τ )dτ = H(s)est . (6)
−∞ −∞
Z ∞
análisis: X(s) = L {x(t)} (s) = x(t) exp(−st)dt. (7)
−∞
Como puede verse, si restringimos la variable compleja s al eje imaginario, s = jω, la ex-
presión se reduce a la transformada de Fourier. Para entender el porqué de esta definición,
vamos a retomar el ejemplo anterior. Supongamos que σ = <{s} > 0, y calculemos la
transformada de Laplace de x(t):
Z ∞
sin(πt) j s − jπ
X(s) = L {x(t)} (s) = exp(−st)dt = log , (8)
0 πt 2π s + jπ
3
Figura 2: Módulo al cuadrado del integrando que define la transformada de Laplace de
x(t), | π1 sin(πt) exp(−σt)|2 , para σ > 0 (izquierda) y σ < 0 (derecha).
1 ∞ 1 ∞ exp(jπt) − exp(−jπt)
Z Z
L {t·x(t)} (s) = sin(πt) exp(−st)dt = exp(−st)dt
π 0 π 0 2j
Z ∞ Z ∞
1 1
= exp((jπ − s)t)dt − exp((−jπ − s)t)dt
2jπ 0 2jπ 0
∞ ∞
1 1
= exp((jπ − s)t) − exp((−jπ − s)t)
2jπ(jπ − s) t=0 2jπ(−jπ − s) t=0
−1 1 1
= + = 2 , (9)
2jπ(jπ − s) 2jπ(−jπ − s) s + π2
que ahora no contiene distribuciones delta de Dirac, sino que es una función relativamente
sencilla en la variable compleja s (con polos en s = ±π).
Por último, notemos que la transformada de Laplace no existirá para cualquier posible
valor de la variable s. En particular, si en el ejemplo anterior tomamos σ = <{s} < 0
las exponenciales exp(−σt) son crecientes (ver Figura 2), y por tanto las integrales no
existen. Este comportamiento es completamente general: como veremos a continuación, la
transformada de Laplace de una cierta señal x(t) estará definida sólo para un subconjunto
de valores s ∈ ROC ⊆ C, aquéllos contenidos en la región de convergencia (ROC).
4
2.1 Regiones de convergencia (ROC)
La ROC de la transformada de Laplace de una señal x(t) se define como el conjunto de
números complejos s ∈ C para los cuales x(t) exp(−st) es absolutamente integrable2 :
Z ∞ Z ∞
∆
ROC = ∀s ∈ C tal que: |x(t) exp(−st)| dt = |x(t)| exp(−σt)dt < ∞ . (10)
−∞ −∞
ya que en este caso la integral sólo converge si σ es estrictamente positiva. Notemos que
la función 1s tiene un polo en s = 0 y por tanto este valor no puede estar contenido en la
ROC, ya que corresponde a una integral no convergente.
2
La condición de “absolutamente integrable” es más restrictiva que la de integrable. A veces se define la
ROC como los s ∈ C para los que x(t) exp(−st) es simplemente integrable, pero esto conduce a conceptos
matemáticos que van mucho más allá de lo que se pretende en esta asignatura. En general, y salvo casos
muy excepcionales, ambas definiciones dan lugar a exactamente la misma ROC.
5
1
Figura 3: Integrando de la transformada de Laplace de la función racional 1+t2
truncada,
para σ > 0 (izquierda), σ = 0 (centro) y σ < 0 (derecha).
ya que en este caso es el lı́mite inferior el que es infinito, y por tanto la convergencia
requiere −(α + σ) > 0 ⇒ σ < −α. Comparando este ejemplo con el ejemplo anterior,
vemos que en ambos casos la expresión funcional de la transformada de Laplace es idéntica,
y ambas transformadas se distinguen tan solo por sus respectivas regiones de convergencia.
Esta consideración da lugar a una importante conclusión:
1
Función racional truncada. Supongamos la señal x(t) = 1+t 2 u(t). Su ROC vendrá
dada por los valores de σ tal que la integral que define la transformada de Laplace es
absolutamente convergente:
Z ∞ Z ∞
1 1
1 + t2 exp(−st)u(t) dt = exp(−σt)dt < ∞. (15)
−∞
0 1 + t2
6
1
Figura 4: Integrando de la transformada de Laplace de la función racional 1+t2
sin
truncar, para σ > 0 (izquierda), σ = 0 (centro) y σ < 0 (derecha).
Incluso ası́, la ROC es en este caso ROC = {∀s ∈ C : <{s} > 0}, excluyendo el eje jω:
1
Función racional. Por último, supongamos la señal x(t) = 1+t 2 , esta vez sin truncar a
los valores mayores que 0. La condición para que s esté contenida en la ROC es ahora:
Z ∞ Z ∞
1 1
1 + t2 exp(−st) dt = exp(−σt)dt < ∞. (17)
−∞
−∞ 1 + t2
Si representamos el integrando de la ecuación anterior para diferentes valores de σ, como se
ha hecho en la Figura 4, podemos comprobar que la integral no es convergente ni para σ > 0
ni para σ < 0, puesto que en ambos casos tenemos colas que crecen de forma exponencial.
Para σ = 0 la integral puede calcularse de forma explı́cita como en la ecuación (16) y el
resultado es π. Puesto que la ROC debe ser un conjunto abierto, pero la transformada de
Laplace sólo es convergente para <{s} = 0, esta función no admite una transformada
de Laplace.
1
Ejercicio: Demostrar que la función racional x(t) = 1+t 2 no admite transformada de
Laplace pero sı́ que admite transformada de Fourier. Calcular esta última.
Para finalizar esta sección, la Figura 5 muestra las ROC (si éstas existen) de cada uno de
los ejemplos propuestos hasta este momento.
3
El que la ROC se defina como un abierto tiene que ver con quePX(s) es analı́tica en todos los puntos
de ésta: si s0 ∈ ROC existe un disco de radio ρ > 0 donde X(s) = ∞ n
n=0 an (s − s0 ) , ∀s : |s − s0 | < ρ. Si
tomamos un s0 en la frontera de la ROC, cualquier disco de radio ρ > 0 centrado en s0 contendrá valores
fuera de la ROC, donde la serie de potencias no será convergente: es decir, X(s) no es analı́tica en la
frontera de la ROC, y por eso ésta no incluye su frontera.
7
Figura 5: Representación gráfica de las zonas donde la integral que define la transformada
de Laplace es absolutamente convergente. Las lineas sólidas indican que la integral es ab-
solutamente convergente para los valores en la frontera, mientras que las lineas punteadas
indican lo contrario. Las cruces rojas son los polos en la transformada de Laplace. Las
ROC respectivas son entonces las zonas sombreadas. Para la función racional no truncada
no existe una ROC, y por tanto esta función no admite transformada de Laplace.
La ROC está formada por bandas paralelas al eje jω. Como ya hemos comentado,
esta propiedad se deduce fácilmente teniendo en cuenta que la convergencia de la trans-
formada de Laplace depende sólo de σ = <{s}.
puesto que s0 ∈ ROC. Pero, puesto que |X(s0 )| < ∞, s0 no puede ser un polo de X(s).
8
acotar el integrando de la siguiente forma:
Z ∞ Z t0
|x(t) exp(−st)| dt = |x(t)| · exp(−σt)dt
−∞ −t0
Z t0 Z t0
≤ |x(t)| · exp(|σ|t0 )dt = exp(|σ|t0 ) |x(t)| , (19)
−t0 −t0
que al ser x(t) absolutamente integrable tendrá un valor finito para todos los valores finitos
de σ. No obstante, si σ → ∞ la integral deja de estar acotada:
que será finita por pertenecer <{s} = σ0 a la ROC. Por tanto, si σ0 ∈ ROC, cualquier
valor σ > σ0 está también contenido en ella y la ROC se extenderá hacia la derecha. No
obstante puede haber un polo en +∞.
Si x(t) es bilateral la ROC es una banda paralela al eje jω, o bien no existe
transformada de Laplace. Una señal bilateral x(t) podrá escribirse como la suma de dos
señales unilaterales:
1. x1 (t) = x(t)u(t), cuya transformada X1 (s) convergerá para <{s} > σ1 .
2. x2 (t) = x(t)u(−t), cuya transformada X2 (s) convergerá para <{s} < σ2 .
Para que la transformada de x(t) sea convergente:
Z ∞ Z 0 Z ∞
|x(t)| exp(−σt)dt = |x1 (t)| exp(−σt)dt + |x2 (t)| exp(−σt)dt < ∞, (21)
−∞ −∞ 0
ası́ que X(s) será convergente para aquellos valores de s que pertenezcan simultáneamente
a las ROC de X1 (s) y X2 (s), es decir, a la banda σ1 < <{s} < σ2 . Podrı́a ocurrir que
no existiese tal banda si σ1 > σ2 , y entonces no existe la transformada de Laplace. En la
siguiente sección veremos algunas consideraciones al respecto.
Ejercicio: Para x(t) = sin(t), escribir esta señal como una señal a derechas x1 (t) más
una señal a izquierdas x2 (t). Calcular σ1 y σ2 y demostrar que esta señal no admite
transformada de Laplace.
9
2.3 Definiciones alternativas de la transformada de Laplace
Como se vio en el ejercicio inmediatamente anterior, la señal bilateral x(t) = sin(t) no ad-
mite transformada de Laplace, a pesar de ser una función sencilla. De hecho hay multitud
de señales bilaterales sencillas que no tienen transformada de Laplace. Para superar esta
dificultad, en algunas disciplinas se usa una definición alternativa a la de la ecuación (7)
para la transformada de Laplace:
Z ∞
Lunilateral {x(t)}(s) = x(t) exp(−st)dt, (22)
0
ya que <{s} = 0 pertenece a la ROC. Es decir, la integral anterior está bien definida y es
finita para cualquier valor de ω, por lo que podemos concluir:
Si para una señal x(t) <{s} = 0 está contenido en la ROC de su transformada de Laplace,
entonces x(t) tiene transformada de Fourier X̃(ω) y ésta es una función en el sentido clásico
(no incluye deltas de Dirac). Además dicha transformada de Fourier se puede calcular a
partir de la transformada de Laplace de x(t), X(ω):
10
Supongamos ahora que <{s} = σ0 pertenece a la ROC, y calculemos la transformada de
Laplace de x(t) en s = σ0 + jω:
Z ∞
L{x(t)}(σ0 + jω) = x(t) exp(−(σ0 + jω)t)dt
−∞
Z ∞
= (x(t) exp(−σ0 t)) exp(−jωt)dt = F{x(t) exp(−σ0 t)}(ω). (25)
−∞
Este resultado permite encontrar una fórmula de inversión para la transformada de Laplace,
como se muestra a continuación.
X(σ0 + jω) = F{x(t) exp(−σ0 t)}(ω) ⇒ x(t) = exp(σ0 t) · F−1 {X(σ0 + jω)}(t)
1 ∞ 1 ∞
Z Z
= exp(σ0 t) X(σ0 + jω) exp(jωt)dω = X(σ0 + jω) exp((σ0 + jω)t)dω. (27)
2π −∞ 2π −∞
La integral previa puede verse como una integral paramétrica en el plano complejo: σ0 +
jω, ω ∈ (−∞, ∞) define una recta vertical (paralela al eje jω) en el plano complejo y
d
contenida en la ROC. Al derivar esta parametrización obtenemos dω s = j, por lo que se
suele dar una expresión alternativa para la transformada inversa de Laplace:
Para cualquier σ0 tal que la recta paralela al eje imaginario <{s} = σ0 está contenida en
la ROC, tenemos la siguiente fórmula de inversión para la transformada de Laplace:
Z ∞ Z σ0 +jω
1 1
sı́ntesis: x(t) = X(σ0 +jω) exp((σ0 +jω)t)dω = X(s) exp(st)ds. (28)
2π −∞ 2πj σ0 −j∞
11
Figura 6: Ejemplo de señal continua a tramos. Este tipo de señales no contiene distribu-
ciones delta de Dirac ni tampoco singularidades.
12
• Las señales continuas a tramos, y las combinaciones de éstas con deltas de Dirac,
tienen transformada de Laplace y cumplen las propiedades que estudiaremos.
y por tanto s0 pertenece a Rz . No obstante, Rz podrı́a ser más amplia que Rx ∩Ry . A
modo de ejemplo sencillo, consideremos x(t) = exp(−t)u(t) e y(t) = − exp(−t)u(t), am-
bas con ROC: Rx = Ry = {∀s ∈ C : <{s} > −1}. Si sumamos estas señales obtenemos
z(t) = 0, cuya ROC es el plano complejo completo.
13
por tanto s ∈ Ry si y sólo si s − s0 ∈ Rx . En otras palabras, Ry es la misma banda
que Rx pero desplazada σ0 hacia la derecha si <{s} > 0 (o |σ0 | hacia la izquierda si
<{s} < 0). Aunque esta propiedad es antagónica a la anterior, con transformada de
Laplace no hablaremos de dualidad. La transformada de Laplace transforma señales
de variable real en funciones de variable compleja, y por tanto no hay dualidad posible.
1
Escalado en el tiempo. Si y(t) = x(at), a ∈ R, entonces Y (s) = |a| X( as ) (demostrar).
Respecto a la ROC:
Z ∞ Z ∞ Z ∞
t dt
|y(t)| exp(−σt)dt = |x(at)| exp(−σt)dt = |x(t)| exp −σ , (35)
−∞ −∞ −∞ a |a|
por tanto s ∈ Ry si y sólo si s/a ∈ Rx . En otras palabras, Ry es una versión escalada de
Rx (expandida si |a| > 1 o contraı́da en caso contrario). También puede aparecer reflejada
en el eje horizontal si a < 0.
con Kx y Ky dos constantes positivas finitas, que han de existir puesto que estamos
asumiendo que s ∈ Rx y s ∈ Ry . Podemos concluir entonces que Rz ⊇ Rx ∩ Ry . Al igual
que pasaba con la suma de transformadas de Laplace, la ROC puede ser mayor que la
intersección de las ROC originales; puesto que tenemos el producto de las transformadas
de Laplace, la ROC resultante será mayor que la intersección cuando los ceros de una de
las transformadas cancelen los polos de la otra. Por ejemplo:
14
Ejercicio: Para las señales x(t) = exp(−2t)u(t) + exp(t)u(−t) e y(t) = 2δ(t) + δ 0 (t):
1. Calcular la transformada de Laplace de ambas señales,
R ∞ X(s) e Y (s), junto con las
ROC correspondientes, Rx y Ry . Recodemos que −∞ δ 0 (t)f (t) = −f 0 (0).
2. Representar los polos y los ceros de X(s) e Y (s) junto con sus ROC.
3. Calcular la convolución entre ambas señales, z(t) = (x ∗ y)(t).
4. Calcular la transformada de Laplace de z(t), Z(S), junto con su ROC, Rz , y com-
probar que de hecho Rz es más amplia que Rx ∩ Ry .
d
Diferenciación en tiempo. Si y(t) = dt x(t), y si tanto x(t) como y(t) admiten transfor-
mada de Laplace, entonces Y (s) = sX(s). En este caso el cálculo de la ROC se hace más
sencillo utilizando la propiedad anterior: podemos ver la derivada como una convolución
con la derivada de la delta de Dirac, cuya transformada de Laplace será:
Z ∞
0 0 d
L{δ (t)}(s) = δ (s) exp(−st)dt = − exp(−st) = s, (38)
−∞ dt t=0
siendo la integral convergente ∀s ∈ C. Derivar x(t) será equivalente a convolucionar con
δ 0 (t), y por tanto Y (s) = sX(s) con Ry ⊇ Rx ∩ C = Rx . Ry es al menos Rx , pudiendo ser
incluso más amplia (considerar por ejemplo el escalón de Heaviside y su derivada, la delta
d
de Dirac). Es importante notar que exigimos que tanto x(t) como y(t) = dt x(t) admi-
1
tan transformada de Laplace; por ejemplo, x(t) = (t−1)1/2 u(t) presenta una singularidad
integrable en t = 1, mientras que su derivada, y(t) = −jδ(t) − 21 (t−1)
1
3/2 , presenta una
d
Diferenciación en s. Si y(t) = t·x(t), entonces Y (s) = − ds X(s) (demostrar). La ROC
Ry será idéntica a Rx :
Z ∞ Z ∞
|y(t)| exp(−σt)dt = |x(t)||t| exp(−σt)dt. (39)
−∞ −∞
El carácter de la integral anterior viene dado por la caı́da de las colas del integrando.
Puesto que una exponencial negativa cae más rápido que cualquier potencia tn , y una
exponencial positiva crece más rápido que cualquier potencia tn , el producto del integrando
por |t| no afecta al carácter de la integral. Esta propiedad se podrá generalizar a:
dn
L{tn x(n)}(s) = (−1)n X(s), ROC = Rx . (40)
dsn
Rt
Integración en tiempo. Supongamos que y(t) = −∞ x(τ )dτ . Al igual que la diferen-
ciación en tiempo, la forma más sencilla de deducir esta propiedad es ver la integración
como una convolución como la función de Heaviside, cuya transformada de Laplace es
U (s) = 1s , <{s} > 0. Por tanto, Y (s) = 1s X(s), y Ry ⊇ Rx ∩ {∀s ∈ C : <{s} > 0}.
Comparemos la deducción de la propiedad de integración en este caso con su deducción
en el caso de la transformada de Fourier. Mientras que aquı́ se ha obtenido de forma
trivial a partir de la propiedad de convolución, su obtención para la transformada de
Fourier es mucho más complicada, requiriendo técnicas avanzadas de variable compleja.
Esta consideración puede servir para ilustrar el hecho de que en general la operativa con
transformada de Laplace será más sencilla.
15
4.3 Teoremas del valor inicial y valor final
Además de las ya mencionadas, vamos a ver dos propiedades adicionales de la transformada
de Laplace que no tienen equivalente alguno en la transformada de Fourier. A diferencia
de las vistas hasta ahora, las dos relaciones siguientes no son válidas para todas las
señales de interés. En concreto, sólo serán aplicables cuando:
• La señal es causal, en el sentido de que x(t) = 0, ∀t < 0.
• La señal no contiene distribuciones del tipo delta de Dirac.
• La señal puede tener unaR t singularidad (limt→0+ x(t) = ±∞) en 0, pero ésta ha de
ser integrable, es decir, 0 0 x(t)dt < ∞.
Entonces se cumple:
Una vez determinada la función e transferencia, en caso de existir ésta, podemos estudiar
las propiedades del sistema. Respecto a la memoria, recordemos que un sistema LTI es
sin memoria si y sólo si su respuesta al impulso es una delta de Dirac o combinación de
derivadas de la delta de Dirac en el origen. Por tanto:
16
Propiedad z(t) Z(s) ROC (Rz )
1
Estas propiedades sólo son válidas para señales que:
1. Sean causales: x(t) = 0, ∀t < 0.
2. No contengan distribuciones del tipo delta de Dirac en el origen.
R
3. No contengan singularidades no integrables en el origen: 0 |x(t)|dt < ∞.
17
Un sistema LTI será estable si y sólo si la ROC de su transformada de Laplace contiene
el eje imaginario jω, o bien si la frontera de la ROC es el propio eje jω y éste no contiene
polos de la transformada de Laplace.
Supongamos que la respuesta al impulso h(t) no contiene distribuciones del tipo delta
de Dirac en el origen, y que tampoco contiene singularidades de orden superior (no inte-
4
Recordemos que con transformada Z la causalidad se traduce en que z → ∞ no sea un polo de H(z)
para un h[n] cualquiera. Esta propiedad no se puede extrapolar al caso de la transformada de Laplace
debido a que las respuestas al impulso de tipo delta de Dirac no son funciones (en el caso discreto δ[n] es
una señal más, sin ninguna peculiaridad). Por ejemplo, podemos calcular la función de transferencia para
h(t) = δ 0 (t): H(s) = s, que presenta un polo en s → ∞ a pesar de ser causal.
18
grables) de forma que su transformada de Laplace es convergente para una determinada
ROC. El sistema es causal si y sólo si la ROC es un semiplano que se extiende hacia la
derecha, y si además se cumple: | lims→∞ X(s)| < ∞. Es anticausal si la ROC es un
semiplano hacia la izquierda y se cumple | lims→−∞ X(s)| < ∞.
19
2
3s −7s+5
Figura 7: La ROC de H(s) = (s−1)(s−2) 2 puede ser una de las tres posibles mostradas
En el caso que nos ocupa, la ROC puede ser una de las tres representadas en la Figura. 7:
• Si asumimos que la respuesta al impulso es causal, la ROC ha de ser un semiplano
hacia la derecha, y por tanto será la zona sombreada en azul. Estamos asumiendo
pues que la ROC de las tres fracciones simples se extiende hacia la derecha, y por lo
tanto tenemos tres exponenciales causales:
h(t) = exp(t)u(t) + 2 exp(2t)u(t) + 3t exp(2t)u(t). (50)
20
• Si asumimos que la respuesta al impulso no es causal ni anticausal, la ROC sólo
puede ser la banda sombreada en verde comprendida entre los dos polos. Estamos
asumiendo pues que la ROC de la primera fracción simple se extiende hacia la
derecha (exponencial causal), y las de las otras fracciones simples hacia la izquierda
(exponenciales anticausales):
Como comentario final, notemos que en este caso se ha identificado causalidad con semi-
planos hacia la derecha y anticausalidad con semiplanos hacia la izquierda. Esto sólo es
válido debido a que estamos tratando con funciones de transferencia que son cocientes de
polinomios en s, es decir, tratamos con sistemas descritos por ecuaciones diferenciales. En
un caso general no se puede asociar causalidad con una determinada forma de la ROC.
21