Transformada inversa de LaplaceLas transformadas de Laplace de las funciones que hemos estudiado en esta página se resumen en la tabla siguiente:

f(t) F(s)=∫0∞e−stf(t)dt

c1f1(t)+c2f2(t) c1F1(s)+c2F2(s)exp(a·t) 1s−a

cos(ωt) ss2+ω2

sin(ωt) ωs2+ω2

tn n!sn+1

exp(at)·f(t)

exp(at)·cos(ωt)F(s-a)

s−a(s−a)2+ω2

u(t-a) exp(-as)/su(t-a)·f(t-a) exp(-as)·F(s)δ(t-a) exp(-as)f'(t) (derivada primera) s·F(s)-f(0)f''(t) (derivada segunda) s2·F(s)-s·f(0)-s·f'(0)g(t)=∫0tf(τ)dτ (integral) F(s)/sf(t)=f(t+p), (función periódica)

11−e−sp∫0pe−stf(t)dt

f(at) 1aF(sa)

tnf(t) (−1)ndndsnF(s)

Como vimos en la página Fracciones polinómicas, la función residue nos permite descomponer una fracción polinómica en suma de fracciones más simples. Modificaremos cada fracción para buscar en la tabla la función f(t), es decir su correspondiente transformada inversa de Laplace, comprobaremos el resultado hecho a mano con la llamada a la funciónilaplace.

1.-Raíces reales distintas

Expresamos la primera fracción en términos de la variable s en vez de x, y sustituímos los decimales periódicos por las fracciones equivalentes.

s−2s3−s2−6s=1151s−3−251s+2+131s

Como vemos en la tabla, la transformada inversa de Laplace se escribe

f(t)=exp(3t)/15-2·exp(-2t)/5+1/3

Obtenemos una expresión similar empleando la función MATLAB ilaplace

>> syms s;>> fs=(s-2)/(s^3-s^2-6*s);>> ft=ilaplace(f)ft =exp(3*t)/15 - 2/(5*exp(2*t)) + 1/3

2.-Raíces complejas distintas

s2−2s+1s3+3s2+4s+2=−3s+7(s+1)2+1+4s+1=−3s+1(s+1)2+12−4(s+1)2+12+4s+1

Como vemos en la tabla, la transformada inversa de Laplace se escribe

f(t)=-3·exp(-t)·cos(t)-4·exp(-t)sin(t)+4·exp(-t)=(4-3·cos(t)-4·sin(t))·exp(-t)

No es necesario sumar las dos fracciones con las raíces conjugadas para convertirla en una única fracción racional

s2−2s+1s3+3s2+4s+2=−1.5+2is−(−1+i)+−1.5−2is−(−1−i)+4s+1

si creamos una nueva entrada en la tabla de las transformadas inversas de Laplace. Sea

F(s)=c+dis−(a+bi)+c−dis−(a−bi)

Como vemos en la tabla, la transformada inversa de Laplace se escribe

f(t)=(c+di)e(a+bi)t+(c−di)e(a−bi)t=eat((c+di)eibt+(c−di)e−ibt)=eat(2ceibt+e−ibt2+2ideibt−e−ibt2)=2eat(c⋅cos(bt)−d⋅sin(bt))Con c=-1.5, d=2, a=-1, b=1, obtenemos

f(t)=(2·exp(-t))(-1.5·cos(t)-2·sin(t))+4·exp(-t))=(4-3·cos(t)-4·sin(t))·exp(-t)

>> syms s;>> fs=(s^2-2*s+1)/(s^3+3*s^2+4*s+2);>> ft=ilaplace(fs)ft =4/exp(t) - (3*(cos(t) + (4*sin(t))/3))/exp(t)>> simplify(ft)ans =-(3*cos(t) + 4*sin(t) - 4)/exp(t)

3.-Raíces repetidas

5s−1s3−3s−2=1s−2+−1s+1+2(s+1)2

Como vemos en la tabla, la transformada inversa de Laplace se escribe

f(t)=exp(2t)-exp(-t)+2t·exp(-t)

>> syms s;>> fs=(5*s-1)/(s^3-3*s-2);>> ft=ilaplace(fs)ft =exp(2*t) - 1/exp(t) + (2*t)/exp(t)

4.-Fracciones impropias

s2+2s+3s2+s−6=1−651s+3+1151s−2

Como vemos en la tabla, la transformada inversa de Laplace se escribe

f(t)=δ(t)+-6·exp(-3t)/5+11·exp(2t)/5

>> syms s;>> fs=(s^2+2*s+3)/(s^2+s-6);>> ft=ilaplace(fs)ft =(11*exp(2*t))/5 - 6/(5*exp(3*t)) + dirac(t)

Ejercicios

Descomponer una fracción en suma de fracciones más simples mediante residue. Determinar las expresiones de las transformadas inversas de Laplace a partir de la tabla y comprobar los resultados con la función ilaplace.3s2+2s+5s3+12s2+44s+48s+3s3+5s2+12s+8s2+3s+1s5+7s4+19s3+25s2+16s+42s4+s3−2ss4+7s3+18s2+20s+8

Ecuaciones diferenciales lineales

1.-Vamos a utilizar la transformada de Laplace para calcular la solución de la ecuación diferencial con las condiciones iniciales especificadas.

d2xdt2+4dxdt+3x=0 t=0{x=3dxdt=1(s2F(s)−sx(0)−x'(0))+4(sF(s)−x(0))+3F(s)=0(s2F(s)−3s−1)+4(sF(s)−3)+3F(s)=0s2F(s)+4sF(s)−3s−13+3F(s)=0F(s)=3s+13s2+4s+3

Transformamos la fracción racional F(s) en suma de fracciones más simples, factorizando el denominador y aplicamos la transformada inversa de Laplace

3s+13s2+4s+3=As+1+Bs+3 {A=5B=−23s+13s2+4s+3=5s+1−2s+3x(t)=5e−t−2e−3t>> syms s;>> Fs=(3*s+13)/(s^2+4*s+3);>> x=ilaplace(Fs)x =5/exp(t) - 2/exp(3*t)

2.-Resolver la ecuación diferencial con las condiciones iniciales especificadas

d2xdt2+4x=g(t) {t 0≤t≤ππ t>πt=0 {x=0dxdt=0

g(t)=(u(t-0)-u(t-π))·t+u(t-π)·π=t-(t-π)·u(t-π)

>> g=@(t) t+(pi-t).*heaviside(t-pi);>> t=linspace(0,10,100);>> plot(t,g(t))

La transformada de Laplace es

(s2F(s)−s⋅x(0)−x'(0))+4F(s)=1s2−e−πss2s2F(s)+4F(s)=1s2(1−e−πs)F(s)=1−e−πss2(s2+4)

Expresamos la fracción en suma de fracciones más simples y calculamos la transformada inversa de Laplace

1s2(s2+4)=As+Bs2+4+Cs+Ds2 ⎧⎩⎨⎪⎪⎪⎪⎪⎪⎪⎪A=0B=1/4C=0D=−1/4F(s)=1−e−πss2(s2+4)=(1−e−πs)(14s2−141s2+4)=14s2−141s2+4−e−πs4s2+14e−πss2+4=141s2−182s2+22−14e−πss2+e−πs82s2

+22f(t)=14t−18sin(2t)−14u(t−π)⋅(t−π)⋅+18u(t−π)⋅sin(2(t−π))>> syms s;>> Fs=(1-exp(-sym('pi')*s))/(s^2*(s^2+4));>> x=ilaplace(Fs)

x =t/4 - sin(2*t)/8 + heaviside(t - pi)*(pi/4 - t/4 + sin(2*t)/8)

Sistema de dos ecuaciones diferenciales

Sea el sistema de dos ecuaciones diferenciales lineales de primer orden con coeficientes constantes que estudiamos en la página Sistema no homogéneo de ecuaciones diferenciales.

{dxdt+2x+4y=1+4tdydt+x−y=32t2t=0 {x=0y=1

Aplicamos la transformada de Laplace a cada una de las ecuaciones y tenemos en cuenta las condiciones iniciales, es decir, el valor inicial x(0) e y(0) en el instante t=0

{sX(s)−x(0)+2X(s)+4Y(s)=1s+4s2sY(s)−y(0)+X(s)−Y(s)=322s3{(s+2)X(s)+4Y(s)=1s+4s2X(s)+(s−1)Y(s)=3s3+1X(s)=−3s3−3s2+4s+12s3(s−2)(s+3)Y(s)=s4+2s3−s2−s+6s3(s−2)(s+3)

Resolvemos el sistema de dos ecuaciones con dos incógnitas, despejando X(s) e Y(s) mediante el operador división por la izquierda. Calculamos la transformada inversa de Laplace mediante la función ilaplace

>> syms s;>> A=[(s+2),4;1,(s-1)];>> b=[1/s+4/s^2;1+3/s^3];>> X=A\bX = -(3*s^3 - 3*s^2 + 4*s + 12)/(s^3*(s^2 + s - 6)) (s^4 + 2*s^3 - s^2 - s + 6)/(s^3*(s^2 + s - 6)) >> ilaplace(X)ans = t - (4*exp(2*t))/5 + 4/(5*exp(3*t)) + t^2 (4*exp(2*t))/5 + 1/(5*exp(3*t)) - t^2/2

Obtenemos el mismo resultado que con la función dsolve.