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.
Si necesitas repasar:
1. Laplace unilateral mira desde el inicio y amortigua antes de sumar
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.
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).
Combinar exponentes
f(t)exp(−3t)=exp(−2t).Integral finita
∫₀∞exp(−2t)dt=[−exp(−2t)/2]₀∞=1/2.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.
La señal original crece, pero el integrando amortiguado decae.
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.
Módulo
|exp(−(2+i)t)|=exp(−2t), integrable en [0,∞).Integrar
F(2+i)=1/(2+i).Separar componentes
Multiplicamos numerador y denominador por 2−i: F(2+i)=(2−i)/5.La transformada puede ser compleja aunque f sea real.
La convergencia absoluta se controla con σ; la oscilación introduce fase.
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).
Cota sencilla
Para t≥0, |f(t)|≤1.El criterio con a=0 garantiza convergencia para Re(s)>0.Integral exacta
F(s)=∫₀∞exp(−(s+2)t)dt=1/(s+2) cuando Re(s)>−2.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.
No se confunde una condición suficiente con una caracterización exacta.
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.
Sustituir formalmente
La fracción evaluada en s=0 vale uno.Revisar la integral original
∫₀ᴿsin t dt=1−cos R.No tiene límite cuando R→∞.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.
La región impide presentar un cálculo algebraico como una integral que no existe.
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
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.
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²}.
Primer paso
I₁=I₀/s=1/s².Segundo paso
I₂=(2/s)I₁=2/s³.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.
L{t²}=2/s³, Re(s)>0.
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).
Primer término
L{t}=1/s²; sustituir s por s+3 da 1/(s+3)², válido para Re(s)>−3.Segundo término
L{sin(2t)}=2/(s²+4).Multiplicar la señal por dos da 4/(s²+4).
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.
Se transforma cada término con sus propios factores y se reúne el marco común.
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.
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.
Derivada temporal
f′=1, cuya transformada es 1/s.Fórmula con dato
sF−f(0+)=s(1/s²+2/s)−2.Simplificar
1/s+2−2=1/s.Sin restar dos quedaría una constante extra que no pertenece a la transformada de la derivada.
El dato inicial sale de la integral; no es una corrección inventada.
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).
Transformar la derivada segunda
L{f″}=s²F−s·1−0.Sustituir F
s³/(s²+4)−s=−4s/(s²+4).Comparar en tiempo
f″=−4cos(2t), cuya transformada es −4F.Ambas cuentas coinciden.
Cada dato conserva su signo y su potencia de s.
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.
Transformar
sY−1+3Y=6/s, luego (s+3)Y=1+6/s.Separar
Y=(s+6)/[s(s+3)]=2/s−1/(s+3).Multiplicar por el denominador comprueba 2(s+3)−s=s+6.Invertir
y(t)=2−exp(−3t).En cero vale uno; y′=3exp(−3t) y y′+3y=6.
La solución combina un nivel final y un transitorio determinado por el dato inicial.
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.
Transformar
s²Y−s−2+4Y=0, así Y=(s+2)/(s²+4).Reconocer términos
Y=s/(s²+4)+2/(s²+4); las inversas son cos(2t) y sin(2t).Verificar
y=cos(2t)+sin(2t).Vale uno en cero, su derivada vale dos allí y su segunda derivada es −4y.
Los dos datos aparecen antes de despejar y ambos se verifican al final.
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.
El escalón importa: no se evalúa f antes de su comienzo.
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.
Transformar
(s+2)Y=exp(−s)/s; Y=exp(−s)/[s(s+2)].Resolver sin retraso
1/[s(s+2)]=(1/2)[1/s−1/(s+2)], cuya inversa es [1−exp(−2t)]/2.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.
La salida no salta al valor final: empieza en cero y se aproxima a un medio.
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.
Convolución
(f*g)(t)=∫₀ᵗexp(−2(t−τ))dτ.Calcular
Con v=t−τ obtenemos ∫₀ᵗexp(−2v)dv=[1−exp(−2t)]/2.Transformadas
F=1/s y G=1/(s+2); su producto es 1/[s(s+2)], igual al de la respuesta recién calculada.
La operación algebraica tiene una interpretación como suma de respuestas pasadas.
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.
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).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.
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.
Apagar la entrada deja un estado que todavía debe decaer.
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).
Señal derecha
f(t)=exp(−2t) para t≥0 y cero antes tiene región Re(s)>−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.
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.
La fórmula y su región identifican juntas la representación bilateral.
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.
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.
Caso oscilante
F=s/(s²+1), de modo que sF→0 al hacer s→0.Pero cos t no tiene límite final.
Detectar la hipótesis ausente
sF conserva polos en ±i, fuera del semiplano estrictamente izquierdo. No estaba permitido usar el criterio del valor final.
El teorema exige revisar sus polos antes de calcular el límite.
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).
Ecuación
Ambas cumplen y′+3y=0 al derivar y cancelar términos.Dato inicial
La primera vale uno en cero y falla. La segunda vale cuatro y pasa.
Laplace
(s+3)Y=4, por tanto Y=4/(s+3).La constante del numerador conservaba precisamente el dato cuatro.
Satisfacer la ecuación no fija por sí solo cuál solución inicial queremos.
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.
Dato
y(0)=2−1=1.Ecuación
y′=3exp(−3t); y′+3y=3exp(−3t)+6−3exp(−3t)=6.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.
La candidata resuelve el problema completo para t≥0.
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.
Datos
Ecuación y′+2y=0, con y(0)=3 y t≥0.Transformar derivada
L{y′}=sY−3.Reunir
sY−3+2Y=0, así (s+2)Y=3.Despejar
Y=3/(s+2), con región Re(s)>−2 para la solución resultante.Invertir
y=3exp(−2t).Comprobar
y(0)=3 y y′+2y=−6exp(−2t)+6exp(−2t)=0.
y(t)=3e^(−2t), t≥0.
y(0)=3 y y′+2y=−6e^(−2t)+6e^(−2t)=0.
EJEMPLO RESUELTO 2
Entrada constante
Resuelve y′+y=1, y(0)=0.
Datos
Ecuación y′+y=1, con y(0)=0.Transformar
sY−0+Y=1/s.Despejar
Y=1/[s(s+1)].Separar fracciones
Y=1/s−1/(s+1).Al reunirlas, el numerador es (s+1)−s=1.Invertir
y=1−exp(−t), con transformada en Re(s)>0.Verificar
y(0)=0 y y′+y=exp(−t)+1−exp(−t)=1.
y(t)=1−e^(−t).
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 →