TRANSFORMADA DE LAPLACE

TRANSFORMADA DE LAPLACE


Apunte
EDO Laplace Inversa Derivada Ejercicios

Estos apuntes fueron basados en la clase del Prof. Arthur Mattuck (1930-2021), especialista en geometría algebraica en el Massachusetts Institute of Technology. [1]

Introducción

Existen muchos misterios en torno a la transformada de Laplace que la gente se pregunta. Es bueno hacernos ideas de dónde vienen las cosas.

¿De dónde viene?

Una buena forma de comenzar a entenderla es con series de potencias:

0anxn\sum_{0}^\infty a_n x^n

Para identificar la suma con los coeficientes:

0anxn=A(x)\sum_{0}^\infty a_n x^n = A(x)

Usando notación de computadora, se piensa como la función de la variable discreta. En otras palabras, es una función que asigna a n=0,1,2,3,n=0,1,2,3,\ldots:

0a(n)xn=A(x)\sum_{0}^\infty a(n) x^n = A(x)

Esto sigue significando que el número real nn sigue asociado.

Tomamos la función discreta que nos da la secuencia de coeficientes de la serie de potencias y la asociamos con la suma de la serie de potencias:

a(n)A(x)a(n) {\color{magenta}\rightsquigarrow} A(x)

Suponiendo que se tiene la función constante 1:

1A(x)1 {\color{magenta}\rightsquigarrow} A(x)

Al reemplazar en la serie, se tiene algo así:

0a(n)xn=10+11+12+13+\sum_{0}^\infty a(n) x^n = 1^0 + 1^1 + 1^2 + 1^3 + \ldots
¿Sabes cuál es el resultado? 111x1 {\color{magenta}\rightsquigarrow} \frac{1}{1-x}

Pero esto no es cierto para todo xx. Solo es cierto cuando xx es tal que esa serie converge. Eso solo es cierto cuando xx se encuentra entre -1 y 1:

111x:x<11 {\color{magenta}\rightsquigarrow} \frac{1}{1-x} \quad : \quad |x| < 1

¿A qué converge? Si xx es mayor que uno, la respuesta es que no converge.

Ahora, para este factorial:

1n!ex\frac{1}{n!} {\color{magenta}\rightsquigarrow} e^x

De este peculiar punto de vista, podemos pensar en una potencia como la suma de la operación en la suma de una serie de potencias, tomando una función discreta definida para enteros positivos o enteros no negativos y haciendo este divertido proceso. Y de ahí surge una función continua de algún tipo.

Observe que lo que entra es la variable nn, pero lo que sale es la variable xx.

BUENO, ESO ES PERFECTAMENTE NATURAL.

Así es como se establece una serie de potencias. La pregunta que hacemos es: esta es una situación discreta. Supongamos que hiciéramos la suma continua en lugar de discreta. Entonces quiero el análogo continuo de lo que hice allí.

¿Qué es el análogo continuo?

Reemplazaremos n=0,1,2,n=0,1,2,\ldots, esto es una variable discreta, por una continua. Lo reemplazamos por una variable continua tt que va de 0 a infinito:

0t<0 \leq t < \infty

Claramente se puede sumar de la forma habitual sobre todos los números reales, pero algo cambia. EL PROCESO PARA REEMPLAZAR LA SUMA SOBRE TODOS LOS NÚMEROS REALES ES LA INTEGRACIÓN:

  0a(t)xtdt=A(x)  \boxed{\;\int_{0}^\infty a(t) x^t dt = A(x)\;}

Al igual que antes, terminará como una función de xx. El problema es que cuando uno quiere integrar o derivar, no quiere tener como base de una exponencial a xx.

Entonces podemos transformar:

x=elnxxt=(elnx)tx = e^{\ln{x}} \\ x^t = (e^{\ln{x}})^t

Si reemplazáramos la integral por 1, no convergería. Por lo tanto, para que converja xx debe ser menor que 1. Esto para que mientras más crezca en tt, se volverá más pequeño.

Recordar que esta es una integral impropia, eso necesita un tratamiento. Si tuviéramos un negativo en xx, deberíamos considerar tratar con imaginarios:

0a(t)xtdt=A(x)\int_{0}^\infty a(t) x^t dt = A(x)

Entonces:

0<x<10 < x < 1

Además, deberíamos considerar los valores en:

x=elnxxt=(elnx)tx = e^{\ln{x}} \\ x^t = (e^{\ln{x}})^t

Para esos valores de xx, al reemplazar siempre tendremos valores negativos:

0<x<1lnx<00 < x < 1 \quad \Rightarrow \quad \ln{x} < 0

Entonces, deberíamos hacerla negativa:

s=lnx-s = \ln{x}

Reemplazando y considerando que los cambios fueron hechos por estética:

0f(t)(es)tdt=F(s)\int_{0}^\infty f(t) (e^{-s})^t dt = F(s) 0f(t)(es)tdt=F(s)\boxed{\int_{0}^\infty f(t) (e^{-s})^t dt = F(s)}

Esto no es más que el análogo continuo de una suma de series de potencias. Considerando que teníamos la versión discreta:

0a(t)xtdt=A(x)\boxed{\int_{0}^\infty a(t) x^t dt = A(x)}

Comienza por una función definida para valores positivos de tt y lo convierte en una función de ss. A esto se le llama:

TRANSFORMADA DE LAPLACE.

Uno de los principales problemas es que no se entiende bien porque terminamos con una función de ss si partimos con una función de tt.

Entonces es importante saber diferenciar entre una transformada y un operador:

f(t)  TRANSFORMADA  F(s)\Large f(t) \longrightarrow \boxed{\color{green}~~\mathbf{TRANSFORMADA}~~} \longrightarrow \color{green}F(s) f(t)         OPERADOR         g(t)\Large f(t) \longrightarrow \boxed{\color{magenta}~~~~~~~~~\mathbf{OPERADOR}~~~~~~~~~} \longrightarrow \color{magenta}g(t)

Esto es una transformación lineal.

Notación

LL{\Large \mathcal{L}} \qquad \lor \qquad {\Large\mathscr{L}}

Dos formas de expresarlo:

L(f(t))=F(s)f(t)F(s)\large \begin{align} \mathscr{L}(\,f(t)\,) =& F(s) \\ f(t) {\color{magenta}\LARGE\rightsquigarrow} & F(s) \end{align}

Transformación lineal

L(f+g)=F(s)+G(s)=L(f)+L(g)L(cf)=cL(f)\mathcal{L}(f+g) = F(s) + G(s) = \mathcal{L}(f) + \mathcal{L}(g) \\ \mathcal{L}(cf) = c \,\mathcal{L}(f)

En otras palabras, se está sumando ambas integrales. La integral en sí misma es un operador lineal.

001?0 \rightsquigarrow 0 \\ 1 \rightsquigarrow ? 0estdt=limR0Restdt=limR0Restdt=limRests0R=limResRses(0)s=limResR01s=1s\begin{split} \int_{0}^\infty e^{-st} dt &= \displaystyle\lim_{R \to \infty} \int_{0}^R e^{-st} dt \\ &= \displaystyle\lim_{R \to \infty} \int_{0}^R e^{-st} dt \\ &= \displaystyle\lim_{R \to \infty} \cfrac{e^{-st}}{-s} \Big|_{0}^R \\ &= \displaystyle\lim_{R \to \infty} \cfrac{e^{-sR}}{-s} \,-\, \cfrac{e^{-s(0)}}{-s} \\ &= \displaystyle\lim_{R \to \infty} \cfrac{{\color{red}\cancel{{\color{gray}e^{-sR}}}^{0}}-1}{-s} \\ &= \frac{1}{s} \end{split} 11s1 \rightsquigarrow \frac{1}{s}

¿Está bien?

La única forma en la que tiende a 0 es cuando ss es un número positivo, por lo tanto s>0s>0:

11s:s>01 \rightsquigarrow \frac{1}{s} \quad : \quad s>0 eat?e^{at} \rightsquigarrow ? \\ eatf(t)?e^{at} f(t) \rightsquigarrow ? 0eatf(t)estdt=0f(t)e(sa)tdt=F(sa)\begin{split} \int_{0}^\infty e^{-at} f(t) \, e^{-st}\, dt &= \int_{0}^\infty f(t) \, e^{-(s-a)t} \, dt \\ &= F(s-a) \end{split} eatf(t)F(sa):s>af(s),s>0e^{at} f(t) \rightsquigarrow F(s-a) \quad : \quad s > a \\ f(s), s > 0 0eatf(t)estdt=0f(t)e(sa)tdt=F(sa):s>a    sa>0\begin{split} \int_{0}^\infty e^{-at} f(t) \, e^{-st}\, dt &= \int_{0}^\infty f(t) \, e ^{-(s-a)t} \, dt \\ &= F(s-a) \quad : \quad s > a ~~\lor~~ s - a > 0 \end{split} eatf(t)F(sa):s>aAsumiendo F(s)  :  s>0e^{at} f(t) \rightsquigarrow F(s-a) \quad : \quad s > a \\ \text{Asumiendo } F(s) \;:\; s > 0

Esto se puede llamar como fórmula de cambio de exponencial. Lo que puedes hacer con uno, lo puedes hacer con el otro. Teníamos que:

eatf(t)?e^{at} f(t) \rightsquigarrow ?

Entonces:

11s:s>01 \rightsquigarrow \frac{1}{s} \quad : \quad s > 0 eatf(t)F(sa)e^{at} f(t) \rightsquigarrow F(s-a) \quad eatf(t)1sa:s>ae^{at} f(t) \rightsquigarrow \frac{1}{s-a} \quad : \quad s > a

Para seno y coseno, la fórmula funciona cuando aa es un número complejo:

e(α+βi)t1s(α+βi):s>ae^{(\alpha+\beta i)t} \rightsquigarrow \frac{1}{s-(\alpha+\beta i)} \quad : \quad s > a sinat?cosat?\sin{at} \rightsquigarrow ? \\ \cos{at} \rightsquigarrow ?

Podemos usar la fórmula de Euler al revés:

cosat=eiat+eiat2\cos{at} = \cfrac{e^{iat} + e^{-iat}}{2}

Simplemente es expresarla al revés, es decir, expresar cosenos y senos en términos de exponenciales complejos:

L(cosat)=12(1sia+1s+ia)\mathscr{L}(\cos{at}) = \frac{1}{2}\left( \frac{1}{s-ia} + \frac{1}{s+ia} \right)

El término a la izquierda es REAL, el de la derecha parece COMPLEJO, pero no lo es. Automáticamente podemos saber que no lo es. Se pueden combinar los términos y calcular. Existen dos métodos para ver que algo es real. Puedes calcularlo y ver que su parte imaginaria es 0, que es el truco que usamos anteriormente. La otra es sin cálculos:

Si cambias ii por i-i y obtienes lo mismo, debe ser real, es decir:

12(1sia+1s+ia)12(1s(i)a+1s+(i)a)\frac{1}{2}\left( \frac{1}{s-ia} + \frac{1}{s+ia} \right) \quad \Rightarrow \quad \frac{1}{2}\left( \frac{1}{s-(-i)a} + \frac{1}{s+(-i)a} \right)

Si calculamos:

L(cosat)=12(1sia+1s+ia)=122ss2+a2\begin{split} \mathscr{L}(\cos{at}) &= \frac{1}{2}\left( \frac{1}{s-ia} + \frac{1}{s+ia} \right) \\ &= \frac{1}{2} \, \frac{2s}{s^2 + a^2} \end{split}

Entonces:

cosatss2+a2:s>0\cos{at} \rightsquigarrow \frac{s}{s^2 + a^2} \quad : \quad s > 0 \\ sinatas2+a2:s>0\sin{at} \rightsquigarrow \frac{a}{s^2 + a^2} \quad : \quad s > 0 \\

Lo difícil de todo esto es recuperar la inversa. Para la inversa de la transformada se podría pensar que existen tablas, pero si las tuviéramos, serían muy extensas y no sería muy útil.

Si tenemos:

1s(s+3)   L1   ?\frac{1}{s(s+3)}~~~\overbrace{\rightsquigarrow}^{\mathscr{L}^{-1}}~~~ ?

Para calcular esta inversa debemos llevarla a algo conocido que esté en las tablas. Para ello se aplican descomposiciones de fracciones parciales.

Primero, cubro el término 1/s1/s:

1s(s+3)=1/3s1s+3\frac{1}{s(s+3)} = \frac{1/3}{s} \cdot \frac{1}{s+3}

Ahora se cubre s+3s+3, y pongo s=3s = -3 porque eso se debe hacer:

1s(s+3)=1/3s1/3s+3\frac{1}{s(s+3)} = \frac{1/3}{s} \cdot \frac{-1/3}{s+3}

Conclusión

La transformada de Laplace es una herramienta poderosa en matemáticas y física que permite convertir funciones del dominio del tiempo en funciones del dominio de la frecuencia, facilitando así la resolución de ecuaciones diferenciales y sistemas dinámicos. Entender su origen y propiedades, como la linealidad y las transformaciones de funciones exponenciales, es esencial para su aplicación efectiva. La descomposición en fracciones parciales es fundamental para encontrar las inversas de estas transformadas, lo que subraya la importancia de esta técnica en el análisis matemático y la ingeniería.


Referencia: