La PizarraMatemática.

Matemáticas para ingeniería

Transformada de Laplace y resolución de ecuaciones

Una ecuación diferencial mezcla una función desconocida con sus derivadas. Laplace reúne esos términos en una ecuación algebraica y conserva los datos iniciales. Antes de usarla como herramienta veremos qué integral calcula, para qué valores existe y de dónde salen los términos que recuerdan cómo empezó el problema.

1. Laplace unilateral mira desde el inicio y amortigua antes de sumar

Definimos F(s)=∫₀∞f(t)exp(−st)dt para funciones ordinarias de t≥0.

Es la versión unilateral: no integra tiempos negativos. Los datos iniciales serán por la derecha, f(0+). Los impulsos en el origen necesitan otra definición y no se añaden silenciosamente a esta convención.

La parte real del parámetro decide la amortiguación
Amortiguar antes de integrarf(t)=exp(t)σ=3: f(t)exp(−3t)=exp(−2t)101Rojo: crece sin límite.Azul: queda área finita.F(3)=∫₀∞exp(−2t)dt=1/2
Para f(t)=exp(t), multiplicar por exp(−σt) con σ=3 produce exp(−2t). La señal original crece, pero la integral amortiguada vale un medio. Las curvas usan los mismos ejes y la misma escala; el gráfico muestra t entre cero y uno y medio, mientras la integral se extiende hasta infinito.
Escribimos s=σ+iω.

El factor exp(−σt) cambia el tamaño y exp(−iωt) oscila. La parte real puede amortiguar el crecimiento lo suficiente para obtener área finita. s es un parámetro de una integral, no una sustitución de la letra t. La transformada es una función de ese parámetro junto con su región de convergencia.

Una señal creciente puede transformarse si se amortigua lo suficiente

f(t)=exp(t), evaluar F(3).

  1. Combinar exponentes

    f(t)exp(−3t)=exp(−2t).
  2. Integral finita

    ∫₀∞exp(−2t)dt=[−exp(−2t)/2]₀∞=1/2.
  3. Generalizar

    Con s complejo, la misma integración da F(s)=1/(s−1) cuando Re(s)>1, donde el término del extremo infinito sí desaparece.

Resultado

La señal original crece, pero el integrando amortiguado decae.

Lo comprobamos

Con s=1 el integrando sería uno y el área divergiría: ese punto no pertenece a la región.

La parte imaginaria no cambia la condición de convergencia absoluta

f(t)=1 y s=2+i.

  1. Módulo

    |exp(−(2+i)t)|=exp(−2t), integrable en [0,∞).
  2. Integrar

    F(2+i)=1/(2+i).
  3. Separar componentes

    Multiplicamos numerador y denominador por 2−i: F(2+i)=(2−i)/5.

    La transformada puede ser compleja aunque f sea real.

Resultado

La convergencia absoluta se controla con σ; la oscilación introduce fase.

Lo comprobamos

No se reemplaza el número complejo por su módulo al calcular el valor de la integral.

2. La región de convergencia forma parte del resultado

Una condición suficiente es que f sea continua por tramos en intervalos finitos, integrable cerca de cero y de orden exponencial: para tiempos grandes, |f(t)|≤Mexp(at). Entonces la integral converge absolutamente si Re(s)>a, al comparar con Mexp(−(Re(s)−a)t). Esa cota da un semiplano suficiente, que puede ser menor que la región real de la función.

La fracción obtenida por álgebra puede tener valores fuera de la región donde representa la integral. Una continuación analítica no es la misma afirmación que convergencia de la integral original. Tampoco toda función continua tiene orden exponencial: hay crecimientos que ninguna amortiguación exp(−σt), con σ finito, consigue controlar.

Una cota cómoda puede ser poco ajustada

f(t)=exp(−2t).

  1. Cota sencilla

    Para t≥0, |f(t)|≤1.
    El criterio con a=0 garantiza convergencia para Re(s)>0.
  2. Integral exacta

    F(s)=∫₀∞exp(−(s+2)t)dt=1/(s+2) cuando Re(s)>−2.
  3. Un punto adicional

    s=−1 está fuera de la garantía inicial Re(s)>0, pero la integral vale uno porque el integrando es exp(−t). La región verdadera era más amplia.

Resultado

No se confunde una condición suficiente con una caracterización exacta.

Lo comprobamos

En Re(s)=−2 falla la convergencia absoluta de esta exponencial.

Que una fracción tenga valor no significa que la integral converja allí

f(t)=sin t tiene F(s)=1/(s²+1) para Re(s)>0.

  1. Sustituir formalmente

    La fracción evaluada en s=0 vale uno.
  2. Revisar la integral original

    ∫₀ᴿsin t dt=1−cos R.
    No tiene límite cuando R→∞.
  3. Conclusión

    Uno es el valor de la expresión continuada, pero no el valor de una integral impropia convergente en s=0. De modo parecido, exp(t²) no tiene transformada absoluta para ningún σ finito: t²−σt acaba siendo al menos t²/2.

Resultado

La región impide presentar un cálculo algebraico como una integral que no existe.

Lo comprobamos

No estamos usando sumación regularizada de la integral oscilante; sería otra noción que habría que declarar.

3. Las primeras fórmulas salen de integrar y conservar los factores

L{1}=1/s para Re(s)>0.

Si I_n=∫₀∞tⁿexp(−st)dt, integrar por partes da I_n=(n/s)I_(n−1) para n≥1 en ese mismo semiplano; los términos de frontera se anulan.

Desde I₀=1/s obtenemos L{tⁿ}=n!/s⁽ⁿ⁺¹⁾.

El factorial viene de repetir las multiplicaciones n, n−1,…,1.

La linealidad permite combinar funciones en una región común. Además L{exp(at)f(t)}(s)=F(s−a), porque los exponentes se suman. Para a real desplaza el semiplano horizontalmente. Por Euler, L{cos(bt)}=s/(s²+b²) y L{sin(bt)}=b/(s²+b²) para b real y Re(s)>0; la fórmula debe conservar el b del numerador del seno.

Deducir una potencia sin memorizar una tabla entera

Calcular L{t²}.

  1. Primer paso

    I₁=I₀/s=1/s².
  2. Segundo paso

    I₂=(2/s)I₁=2/s³.
  3. Comprobar un valor

    En s=2, ∫₀∞t²exp(−2t)dt=2/8=1/4.

    El dos procede de derivar t², no de una normalización de Fourier.

Resultado

L{t²}=2/s³, Re(s)>0.

Lo comprobamos

Para n=0 la función constante uno sigue teniendo 1/s; no se aplica una recurrencia con índice negativo.

Combinar desplazamiento exponencial y linealidad

f(t)=t exp(−3t)+2sin(2t).

  1. Primer término

    L{t}=1/s²; sustituir s por s+3 da 1/(s+3)², válido para Re(s)>−3.
  2. Segundo término

    L{sin(2t)}=2/(s²+4).

    Multiplicar la señal por dos da 4/(s²+4).

  3. Región común

    F(s)=1/(s+3)²+4/(s²+4) para Re(s)>0.

    La parte sinusoidal obliga a mantener ese semiplano para la integral absoluta conjunta.

Resultado

Se transforma cada término con sus propios factores y se reúne el marco común.

Lo comprobamos

El dos de la altura y el dos de la frecuencia son datos diferentes, aunque aquí coincidan.

4. La derivada conserva el valor inicial como un término de frontera

Con regularidad suficiente, como f localmente absolutamente continua y crecimiento adecuado de f y f′, integrar por partes da L{f′}=[f exp(−st)]₀∞+sF=sF−f(0+), en un semiplano donde la frontera infinita se anula. El signo menos viene del extremo inferior. Omitirlo normalmente cambia el problema inicial.

Repitiendo bajo las hipótesis correspondientes, L{f″}=s²F−sf(0+)−f′(0+).

Cada dato debe tener su lugar. Una discontinuidad requiere tratar los tramos y, si corresponde, impulsos distribucionales: no se aplica la fórmula ordinaria a una señal con salto fingiendo que es derivable en todo punto.

Ver la cancelación exacta del término inicial

f(t)=t+2, con F=1/s²+2/s.

  1. Derivada temporal

    f′=1, cuya transformada es 1/s.
  2. Fórmula con dato

    sF−f(0+)=s(1/s²+2/s)−2.
  3. Simplificar

    1/s+2−2=1/s.

    Sin restar dos quedaría una constante extra que no pertenece a la transformada de la derivada.

Resultado

El dato inicial sale de la integral; no es una corrección inventada.

Lo comprobamos

El cálculo se hace donde Re(s)>0 y la frontera amortiguada se anula.

La segunda derivada recuerda dos datos

f(t)=cos(2t); f(0)=1, f′(0)=0 y F=s/(s²+4).

  1. Transformar la derivada segunda

    L{f″}=s²F−s·1−0.
  2. Sustituir F

    s³/(s²+4)−s=−4s/(s²+4).
  3. Comparar en tiempo

    f″=−4cos(2t), cuya transformada es −4F.

    Ambas cuentas coinciden.

Resultado

Cada dato conserva su signo y su potencia de s.

Lo comprobamos

El término de f(0) es s, no uno; corresponde a la segunda derivada.

5. Resolvemos una ecuación completa, no solo una fracción

El orden es: escribir ecuación y datos; transformar cada término; reunir los factores de Y(s); despejar; invertir; comprobar en tiempo. Una tabla ayuda a reconocer formas conocidas, pero no reemplaza ese recorrido. Las fracciones simples se eligen para convertir el resultado en una suma de transformadas ya comprendidas.

La unicidad de Laplace requiere una clase adecuada, como funciones continuas por tramos y de orden exponencial; las integrales no distinguen cambios aislados. Para una solución continua de una ecuación ordinaria con datos fijados, la comprobación temporal elimina esa ambigüedad puntual. El dominio t≥0 y el semiplano deben conservarse durante la resolución.

Entrada constante con estado inicial distinto de cero

y′+3y=6, y(0)=1.

  1. Transformar

    sY−1+3Y=6/s, luego (s+3)Y=1+6/s.
  2. Separar

    Y=(s+6)/[s(s+3)]=2/s−1/(s+3).
    Multiplicar por el denominador comprueba 2(s+3)−s=s+6.
  3. Invertir

    y(t)=2−exp(−3t).
    En cero vale uno; y′=3exp(−3t) y y′+3y=6.
Resultado

La solución combina un nivel final y un transitorio determinado por el dato inicial.

Lo comprobamos

La integral de Y corresponde a Re(s)>0 por la constante final dos.

Dos datos iniciales fijan las dos componentes de una oscilación

y″+4y=0, y(0)=1, y′(0)=2.

  1. Transformar

    s²Y−s−2+4Y=0, así Y=(s+2)/(s²+4).
  2. Reconocer términos

    Y=s/(s²+4)+2/(s²+4); las inversas son cos(2t) y sin(2t).
  3. Verificar

    y=cos(2t)+sin(2t).

    Vale uno en cero, su derivada vale dos allí y su segunda derivada es −4y.

Resultado

Los dos datos aparecen antes de despejar y ambos se verifican al final.

Lo comprobamos

Si se omitiera el término −2, desaparecería el seno y fallaría el dato de la derivada.

6. Un retraso necesita el escalón que activa la señal

Sea u(t−T) el escalón que vale cero antes de T y uno después; su valor aislado en T no afecta a la transformada.

Para T≥0, L{u(t−T)f(t−T)}=exp(−sT)F(s).
Se demuestra integrando desde T y poniendo v=t−T.

El escalón importa: no se evalúa f antes de su comienzo.

Retrasar la entrada no equivale a saltar instantáneamente al valor final
La entrada empieza en t=1Entrada: u(t−1)1y′+2y=u(t−1), y(0)=01/214La salida sube sin saltar.
El escalón entra en t=1. La solución es cero antes y [1−exp(−2(t−1))]/2 después. Es continua en uno y se acerca al nivel un medio. Se han dibujado la entrada y la salida con la misma escala vertical por unidad; la línea discontinua inferior marca el límite de la salida.

La convolución causal (f*g)(t)=∫₀ᵗf(τ)g(t−τ)dτ se transforma en FG cuando el amortiguamiento permite las integrales y Fubini, por ejemplo para funciones de orden exponencial en un semiplano común suficiente. En una ecuación con dato inicial, la respuesta a la entrada se añade a la contribución del estado inicial; el producto FG solo no borra esa información.

La respuesta empieza cuando llega la entrada

y′+2y=u(t−1), y(0)=0.

  1. Transformar

    (s+2)Y=exp(−s)/s; Y=exp(−s)/[s(s+2)].
  2. Resolver sin retraso

    1/[s(s+2)]=(1/2)[1/s−1/(s+2)], cuya inversa es [1−exp(−2t)]/2.
  3. Activar y desplazar

    y=0 antes de uno; después y=[1−exp(−2(t−1))]/2.

    Es continua en t=1 y su derivada cambia para satisfacer la entrada por tramos.

Resultado

La salida no salta al valor final: empieza en cero y se aproxima a un medio.

Lo comprobamos

La ecuación se satisface en los intervalos abiertos de cada tramo, con continuidad en el instante de activación.

El producto de transformadas se puede comprobar integrando

f(t)=1, g(t)=exp(−2t), para t≥0.

  1. Convolución

    (f*g)(t)=∫₀ᵗexp(−2(t−τ))dτ.
  2. Calcular

    Con v=t−τ obtenemos ∫₀ᵗexp(−2v)dv=[1−exp(−2t)]/2.
  3. Transformadas

    F=1/s y G=1/(s+2); su producto es 1/[s(s+2)], igual al de la respuesta recién calculada.
Resultado

La operación algebraica tiene una interpretación como suma de respuestas pasadas.

Lo comprobamos

El límite superior t refleja causalidad; no se integra sobre tiempos futuros.

Una entrada que se apaga requiere restar otra respuesta retrasada

y′+2y=u(t−1)−u(t−3), y(0)=0.

  1. Superponer

    Sea r(v)=[1−exp(−2v)]/2 para v≥0.
    Entonces y=u(t−1)r(t−1)−u(t−3)r(t−3).
  2. Entre uno y tres

    Solo está activa la primera respuesta. Después de tres, ambas contribuyen y y=[exp(−2(t−3))−exp(−2(t−1))]/2.

  3. Comprobar apagado

    Para t>3 se cumple y′+2y=0 y la salida decae a cero.

    En tres es continua: no se reinicia artificialmente el valor de y.

Resultado

Apagar la entrada deja un estado que todavía debe decaer.

Lo comprobamos

El signo menos corresponde al escalón que retira la entrada, no a borrar toda la historia previa.

7. Los polos se interpretan con región, versión y teorema

La transformada bilateral integra toda ℝ y puede asociar una misma fracción a señales de lados temporales distintos. La unilateral de esta lección integra desde cero y conserva datos iniciales en las derivadas. No mezclamos una región bilateral con fórmulas unilaterales sin indicar el cambio.

Para una señal continua por la derecha en cero y de crecimiento adecuado, el teorema del valor inicial permite lim real s→+∞ sF(s)=f(0+). El valor final tiene condiciones distintas. En el caso racional propio ordinario, un criterio suficiente es que todos los polos de sF estén estrictamente en el semiplano izquierdo; entonces lim_t→∞f(t)=lim_s→0sF(s).

Polos imaginarios o derechos invalidan ese criterio. La estabilidad interna de un modelo de estados necesita además estudiar sus modos, incluidos los que una función externa podría ocultar.

La misma fracción bilateral representa dos soportes distintos

F(s)=1/(s+2).

  1. Señal derecha

    f(t)=exp(−2t) para t≥0 y cero antes tiene región Re(s)>−2.

  2. Señal izquierda

    g(t)=−exp(−2t) para t<0 y cero después también integra a 1/(s+2), ahora en Re(s)<−2.

  3. Diferencia unilateral

    La transformada unilateral de esa g ordinaria sería cero, porque g es cero para t≥0. No se traslada la fracción bilateral ignorando el soporte.

Resultado

La fórmula y su región identifican juntas la representación bilateral.

Lo comprobamos

No se toma un punto a ambos lados de la frontera como si fuera la misma integral convergente.

Un límite algebraico puede acertar una señal y fallar otra

Comparar y=2−exp(−3t) con f=cos t.

  1. Caso convergente

    Y=2/s−1/(s+3). sY solo tiene polo en −3 y lim_s→0sY=2: coincide con el valor final temporal.

  2. Caso oscilante

    F=s/(s²+1), de modo que sF→0 al hacer s→0.

    Pero cos t no tiene límite final.

  3. Detectar la hipótesis ausente

    sF conserva polos en ±i, fuera del semiplano estrictamente izquierdo. No estaba permitido usar el criterio del valor final.

Resultado

El teorema exige revisar sus polos antes de calcular el límite.

Lo comprobamos

En cambio, lim_s→+∞sF=1 sí recupera el valor inicial de cos t bajo las hipótesis iniciales.

8. La ecuación y todos los datos se comprueban en el tiempo

Después de invertir, derivamos la función candidata y la sustituimos en la ecuación. Revisamos cada dato inicial y, para entradas por tramos, cada intervalo y sus uniones. Una respuesta puede satisfacer la ecuación homogénea y fallar el dato inicial; otra puede coincidir al inicio y fallar la ecuación.

Son controles independientes.

También es útil contrastar con un método de EDO, como factor integrante, o volver a transformar la respuesta. Se mantienen región y clase de funciones. Laplace organiza la información del problema completo; no reemplaza una solución por una fracción sin interpretación temporal.

Dos exponentes pueden resolver la ecuación, pero no el mismo dato

y′+3y=0, y(0)=4. Candidatas exp(−3t) y 4exp(−3t).

  1. Ecuación

    Ambas cumplen y′+3y=0 al derivar y cancelar términos.
  2. Dato inicial

    La primera vale uno en cero y falla. La segunda vale cuatro y pasa.

  3. Laplace

    (s+3)Y=4, por tanto Y=4/(s+3).

    La constante del numerador conservaba precisamente el dato cuatro.

Resultado

Satisfacer la ecuación no fija por sí solo cuál solución inicial queremos.

Lo comprobamos

El factor integrante también da y exp(3t)=4, una comprobación independiente.

Comprobar una solución con entrada y estado inicial

y=2−exp(−3t), propuesta para y′+3y=6, y(0)=1.

  1. Dato

    y(0)=2−1=1.
  2. Ecuación

    y′=3exp(−3t); y′+3y=3exp(−3t)+6−3exp(−3t)=6.
  3. Largo plazo y transformada

    y→2 y Y=2/s−1/(s+3) en Re(s)>0.

    El estado inicial, la ecuación y el valor final encajan sin omitir condiciones.

Resultado

La candidata resuelve el problema completo para t≥0.

Lo comprobamos

No se exige que Y tenga región Re(s)>−3: la constante dos impone Re(s)>0 para la integral unilateral absoluta.

Ahora lo resolvemos de principio a fin.

EJEMPLO RESUELTO 1

Decaimiento con dato inicial

Resuelve y′+2y=0, y(0)=3 mediante Laplace.

  1. Datos

    Ecuación y′+2y=0, con y(0)=3 y t≥0.
  2. Transformar derivada

    L{y′}=sY−3.
  3. Reunir

    sY−3+2Y=0, así (s+2)Y=3.
  4. Despejar

    Y=3/(s+2), con región Re(s)>−2 para la solución resultante.
  5. Invertir

    y=3exp(−2t).
  6. Comprobar

    y(0)=3 y y′+2y=−6exp(−2t)+6exp(−2t)=0.
Resultado

y(t)=3e^(−2t), t≥0.

Lo comprobamos

y(0)=3 y y′+2y=−6e^(−2t)+6e^(−2t)=0.

EJEMPLO RESUELTO 2

Entrada constante

Resuelve y′+y=1, y(0)=0.

  1. Datos

    Ecuación y′+y=1, con y(0)=0.
  2. Transformar

    sY−0+Y=1/s.
  3. Despejar

    Y=1/[s(s+1)].
  4. Separar fracciones

    Y=1/s−1/(s+1).
    Al reunirlas, el numerador es (s+1)−s=1.
  5. Invertir

    y=1−exp(−t), con transformada en Re(s)>0.
  6. Verificar

    y(0)=0 y y′+y=exp(−t)+1−exp(−t)=1.
Resultado

y(t)=1−e^(−t).

Lo comprobamos

y(0)=0 y y′+y=e^(−t)+1−e^(−t)=1.

Los símbolos que hemos usado

s=σ+iω
Variable compleja: parte real de amortiguamiento y parte imaginaria de oscilación.
F(s)
Transformada unilateral de f desde t=0.
Orden exponencial
Cota de crecimiento por M exp(at) para tiempos suficientemente grandes.
f(0+)
Límite por la derecha en el instante inicial.
Región de convergencia
Valores de s donde la integral definida converge en el marco considerado.
Inversión
Recuperar una función temporal a partir de F con sus condiciones.
Escalón retardado
Función que activa una señal a partir de un tiempo fijado.

Ahora te toca a ti.

Tienes cinco ejercicios. Inténtalos; las pistas te ayudan si te atascas.

Practicar este tema →

Para ampliar