The Laplace transform of f ( t ) f(t) f ( t ) (defined for t ≥ 0 t \geq 0 t ≥ 0 ) is
L { f ( t ) } = F ( s ) = ∫ 0 ∞ e − s t f ( t ) d t \mathcal{L}\{f(t)\} = F(s) = \int_0^{\infty} e^{-st} f(t)\, dt L { f ( t )} = F ( s ) = ∫ 0 ∞ e − s t f ( t ) d t
The transform exists when f f f is piecewise continuous on [ 0 , ∞ ) [0, \infty) [ 0 , ∞ ) and of exponential order : ∣ f ( t ) ∣ ≤ M e a t |f(t)| \leq Me^{at} ∣ f ( t ) ∣ ≤ M e a t for some M , a > 0 M, a > 0 M , a > 0 .
Theorem 5.1 (Linearity). L { a f + b g } = a L { f } + b L { g } \mathcal{L}\{af + bg\} = a\mathcal{L}\{f\} + b\mathcal{L}\{g\} L { a f + b g } = a L { f } + b L { g } .
Theorem 5.2 (First Shifting). L { e a t f ( t ) } = F ( s − a ) \mathcal{L}\{e^{at}f(t)\} = F(s - a) L { e a t f ( t )} = F ( s − a ) .
Theorem 5.3 (Second Shifting). L { u c ( t ) f ( t − c ) } = e − c s F ( s ) \mathcal{L}\{u_c(t)f(t - c)\} = e^{-cs}F(s) L { u c ( t ) f ( t − c )} = e − cs F ( s ) Where u c ( t ) u_c(t) u c ( t ) is The unit step function.
Theorem 5.4 (Derivative). L { f " ( t ) } = s F ( s ) − f ( 0 ) \mathcal{L}\{f"(t)\} = sF(s) - f(0) L { f " ( t )} = s F ( s ) − f ( 0 ) .
Theorem 5.5 (n-th Derivative). L { f ( n ) ( t ) } = s n F ( s ) − s n − 1 f ( 0 ) − ⋯ − f ( n − 1 ) ( 0 ) \mathcal{L}\{f^{(n)}(t)\} = s^n F(s) - s^{n-1}f(0) - \cdots - f^{(n-1)}(0) L { f ( n ) ( t )} = s n F ( s ) − s n − 1 f ( 0 ) − ⋯ − f ( n − 1 ) ( 0 ) .
Theorem 5.6 (Integration). L { ∫ 0 t f ( τ ) d τ } = F ( s ) s \mathcal{L}\left\{\int_0^t f(\tau)\, d\tau\right\} = \frac{F(s)}{s} L { ∫ 0 t f ( τ ) d τ } = s F ( s ) .
Theorem 5.7 (Convolution). L { f ∗ g } = F ( s ) G ( s ) \mathcal{L}\{f * g\} = F(s)G(s) L { f ∗ g } = F ( s ) G ( s ) Where ( f ∗ g ) ( t ) = ∫ 0 t f ( τ ) g ( t − τ ) d τ (f * g)(t) = \int_0^t f(\tau)g(t - \tau)\, d\tau ( f ∗ g ) ( t ) = ∫ 0 t f ( τ ) g ( t − τ ) d τ .
Proof of Linearity. L { a f + b g } = ∫ 0 ∞ e − s t [ a f ( t ) + b g ( t ) ] d t = a ∫ 0 ∞ e − s t f ( t ) d t + b ∫ 0 ∞ e − s t g ( t ) d t = a F ( s ) + b G ( s ) \mathcal{L}\{af + bg\} = \int_0^{\infty} e^{-st}[af(t) + bg(t)]\, dt = a\int_0^{\infty} e^{-st}f(t)\, dt + b\int_0^{\infty} e^{-st}g(t)\, dt = aF(s) + bG(s) L { a f + b g } = ∫ 0 ∞ e − s t [ a f ( t ) + b g ( t )] d t = a ∫ 0 ∞ e − s t f ( t ) d t + b ∫ 0 ∞ e − s t g ( t ) d t = a F ( s ) + b G ( s ) . ■ \blacksquare ■
Proof of First Shifting Theorem. L { e a t f ( t ) } = ∫ 0 ∞ e − s t e a t f ( t ) d t = ∫ 0 ∞ e − ( s − a ) t f ( t ) d t = F ( s − a ) \mathcal{L}\{e^{at}f(t)\} = \int_0^{\infty} e^{-st} e^{at} f(t)\, dt = \int_0^{\infty} e^{-(s-a)t} f(t)\, dt = F(s - a) L { e a t f ( t )} = ∫ 0 ∞ e − s t e a t f ( t ) d t = ∫ 0 ∞ e − ( s − a ) t f ( t ) d t = F ( s − a ) . ■ \blacksquare ■
Proof of Derivative Property. Integrating by parts with u = e − s t u = e^{-st} u = e − s t , d v = f ′ ( t ) d t dv = f'(t)\, dt d v = f ′ ( t ) d t :
L { f ′ ( t ) } = [ e − s t f ( t ) ] 0 ∞ + s ∫ 0 ∞ e − s t f ( t ) d t = − f ( 0 ) + s F ( s ) = s F ( s ) − f ( 0 ) \mathcal{L}\{f'(t)\} = \left[e^{-st}f(t)\right]_0^{\infty} + s\int_0^{\infty} e^{-st}f(t)\, dt = -f(0) + sF(s) = sF(s) - f(0) L { f ′ ( t )} = [ e − s t f ( t ) ] 0 ∞ + s ∫ 0 ∞ e − s t f ( t ) d t = − f ( 0 ) + s F ( s ) = s F ( s ) − f ( 0 ) .
The boundary term vanishes at ∞ \infty ∞ because f f f is of exponential order. ■ \blacksquare ■
Proof of Integration Property. Let g ( t ) = ∫ 0 t f ( τ ) d τ g(t) = \int_0^t f(\tau)\, d\tau g ( t ) = ∫ 0 t f ( τ ) d τ . Then g ′ ( t ) = f ( t ) g'(t) = f(t) g ′ ( t ) = f ( t ) and g ( 0 ) = 0 g(0) = 0 g ( 0 ) = 0 . L { f ( t ) } = L { g ′ ( t ) } = s L { g ( t ) } − g ( 0 ) = s L { g ( t ) } \mathcal{L}\{f(t)\} = \mathcal{L}\{g'(t)\} = s\mathcal{L}\{g(t)\} - g(0) = s\mathcal{L}\{g(t)\} L { f ( t )} = L { g ′ ( t )} = s L { g ( t )} − g ( 0 ) = s L { g ( t )} . Therefore L { g ( t ) } = F ( s ) / s \mathcal{L}\{g(t)\} = F(s)/s L { g ( t )} = F ( s ) / s . ■ \blacksquare ■
f ( t ) f(t) f ( t ) F ( s ) = L { f ( t ) } F(s) = \mathcal{L}\{f(t)\} F ( s ) = L { f ( t )} 1 1 1 1 / s 1/s 1/ s t n t^n t n n ! / s n + 1 n!/s^{n+1} n ! / s n + 1 e a t e^{at} e a t 1 / ( s − a ) 1/(s - a) 1/ ( s − a ) t n e a t t^n e^{at} t n e a t n ! / ( s − a ) n + 1 n!/(s - a)^{n+1} n ! / ( s − a ) n + 1 sin ( b t ) \sin(bt) sin ( b t ) b / ( s 2 + b 2 ) b/(s^2 + b^2) b / ( s 2 + b 2 ) cos ( b t ) \cos(bt) cos ( b t ) s / ( s 2 + b 2 ) s/(s^2 + b^2) s / ( s 2 + b 2 ) e a t sin ( b t ) e^{at}\sin(bt) e a t sin ( b t ) b / ( ( s − a ) 2 + b 2 ) b/((s-a)^2 + b^2) b / (( s − a ) 2 + b 2 ) e a t cos ( b t ) e^{at}\cos(bt) e a t cos ( b t ) ( s − a ) / ( ( s − a ) 2 + b 2 ) (s-a)/((s-a)^2 + b^2) ( s − a ) / (( s − a ) 2 + b 2 ) u c ( t ) u_c(t) u c ( t ) e − c s / s e^{-cs}/s e − cs / s δ ( t − c ) \delta(t - c) δ ( t − c ) e − c s e^{-cs} e − cs
Procedure:
Take the Laplace transform of both sides of the ODE. Solve the resulting algebraic equation for Y ( s ) = L { y ( t ) } Y(s) = \mathcal{L}\{y(t)\} Y ( s ) = L { y ( t )} . Find the inverse Laplace transform to obtain y ( t ) y(t) y ( t ) . Problem. Solve y ′ ′ + 4 y = sin ( 2 t ) y'' + 4y = \sin(2t) y ′′ + 4 y = sin ( 2 t ) , y ( 0 ) = 0 y(0) = 0 y ( 0 ) = 0 , y ′ ( 0 ) = 0 y'(0) = 0 y ′ ( 0 ) = 0 .
Solution. Taking Laplace transforms:
s 2 Y + 4 Y = 2 s 2 + 4 s^2 Y + 4Y = \frac{2}{s^2 + 4} s 2 Y + 4 Y = s 2 + 4 2
( s 2 + 4 ) Y = 2 s 2 + 4 (s^2 + 4)Y = \frac{2}{s^2 + 4} ( s 2 + 4 ) Y = s 2 + 4 2
Y = 2 ( s 2 + 4 ) 2 Y = \frac{2}{(s^2 + 4)^2} Y = ( s 2 + 4 ) 2 2
To invert, use the convolution theorem: L − 1 { 1 s 2 + 4 } = 1 2 sin ( 2 t ) \mathcal{L}^{-1}\left\{\frac{1}{s^2 + 4}\right\} = \frac{1}{2}\sin(2t) L − 1 { s 2 + 4 1 } = 2 1 sin ( 2 t ) .
y ( t ) = 1 2 sin ( 2 t ) ∗ 1 2 sin ( 2 t ) = 1 4 ∫ 0 t sin ( 2 τ ) sin ( 2 ( t − τ ) ) d τ y(t) = \frac{1}{2}\sin(2t) * \frac{1}{2}\sin(2t) = \frac{1}{4}\int_0^t \sin(2\tau)\sin(2(t - \tau))\, d\tau y ( t ) = 2 1 sin ( 2 t ) ∗ 2 1 sin ( 2 t ) = 4 1 ∫ 0 t sin ( 2 τ ) sin ( 2 ( t − τ )) d τ
Using sin A sin B = 1 2 ( cos ( A − B ) − cos ( A + B ) ) \sin A \sin B = \frac{1}{2}(\cos(A - B) - \cos(A + B)) sin A sin B = 2 1 ( cos ( A − B ) − cos ( A + B )) :
y ( t ) = 1 8 ∫ 0 t [ cos ( 4 τ − 2 t ) − cos ( 2 t ) ] d τ = 1 8 [ sin ( 4 τ − 2 t ) 4 ] 0 t − t 8 cos ( 2 t ) y(t) = \frac{1}{8}\int_0^t [\cos(4\tau - 2t) - \cos(2t)]\, d\tau = \frac{1}{8}\left[\frac{\sin(4\tau - 2t)}{4}\right]_0^t - \frac{t}{8}\cos(2t) y ( t ) = 8 1 ∫ 0 t [ cos ( 4 τ − 2 t ) − cos ( 2 t )] d τ = 8 1 [ 4 s i n ( 4 τ − 2 t ) ] 0 t − 8 t cos ( 2 t )
= 1 32 [ sin ( 2 t ) − sin ( − 2 t ) ] − t 8 cos ( 2 t ) = sin ( 2 t ) 16 − t cos ( 2 t ) 8 = \frac{1}{32}[\sin(2t) - \sin(-2t)] - \frac{t}{8}\cos(2t) = \frac{\sin(2t)}{16} - \frac{t\cos(2t)}{8} = 32 1 [ sin ( 2 t ) − sin ( − 2 t )] − 8 t cos ( 2 t ) = 16 s i n ( 2 t ) − 8 t c o s ( 2 t )
y ( t ) = sin ( 2 t ) − 2 t cos ( 2 t ) 16 y(t) = \frac{\sin(2t) - 2t\cos(2t)}{16} y ( t ) = 16 s i n ( 2 t ) − 2 t c o s ( 2 t ) . ■ \blacksquare ■
Problem. Find \mathcal{L}^{-1}\left\{\frac{1}{s(s^2 + 4)} Ight\} .
Solution Solution. Partial fractions: 1 s ( s 2 + 4 ) = A s + B s + C s 2 + 4 \frac{1}{s(s^2 + 4)} = \frac{A}{s} + \frac{Bs + C}{s^2 + 4} s ( s 2 + 4 ) 1 = s A + s 2 + 4 B s + C .
1 = A ( s 2 + 4 ) + ( B s + C ) s = ( A + B ) s 2 + C s + 4 A 1 = A(s^2 + 4) + (Bs + C)s = (A + B)s^2 + Cs + 4A 1 = A ( s 2 + 4 ) + ( B s + C ) s = ( A + B ) s 2 + C s + 4 A .
s 2 s^2 s 2 : A + B = 0 A + B = 0 A + B = 0 , s s s : C = 0 C = 0 C = 0 Constant: 4 A = 1 ⟹ A = 1 / 4 4A = 1 \implies A = 1/4 4 A = 1 ⟹ A = 1/4 , B = − 1 / 4 B = -1/4 B = − 1/4 .
1 s ( s 2 + 4 ) = 1 / 4 s − s / 4 s 2 + 4 \frac{1}{s(s^2 + 4)} = \frac{1/4}{s} - \frac{s/4}{s^2 + 4} s ( s 2 + 4 ) 1 = s 1/4 − s 2 + 4 s /4 .
f ( t ) = 1 4 − 1 4 cos 2 t f(t) = \frac{1}{4} - \frac{1}{4}\cos 2t f ( t ) = 4 1 − 4 1 cos 2 t . ■ \blacksquare ■
### 5.7 Convolution Theorem
Theorem 5.8 (Convolution). If F ( s ) = L { f ( t ) } F(s) = \mathcal{L}\{f(t)\} F ( s ) = L { f ( t )} and G ( s ) = L { g ( t ) } G(s) = \mathcal{L}\{g(t)\} G ( s ) = L { g ( t )} Then
L { f ∗ g } = F ( s ) G ( s ) \mathcal{L}\{f * g\} = F(s)G(s) L { f ∗ g } = F ( s ) G ( s )
Where ( f ∗ g ) ( t ) = ∫ 0 t f ( τ ) g ( t − τ ) d τ (f * g)(t) = \int_0^t f(\tau)g(t - \tau)\, d\tau ( f ∗ g ) ( t ) = ∫ 0 t f ( τ ) g ( t − τ ) d τ .
Proof. F ( s ) G ( s ) = ∫ 0 ∞ e − s τ f ( τ ) d τ ⋅ ∫ 0 ∞ e − s u g ( u ) d u F(s)G(s) = \int_0^{\infty} e^{-s\tau} f(\tau)\, d\tau \cdot \int_0^{\infty} e^{-su} g(u)\, du F ( s ) G ( s ) = ∫ 0 ∞ e − s τ f ( τ ) d τ ⋅ ∫ 0 ∞ e − s u g ( u ) d u
= ∫ 0 ∞ ∫ 0 ∞ e − s ( τ + u ) f ( τ ) g ( u ) d u d τ = \int_0^{\infty} \int_0^{\infty} e^{-s(\tau+u)} f(\tau) g(u)\, du\, d\tau = ∫ 0 ∞ ∫ 0 ∞ e − s ( τ + u ) f ( τ ) g ( u ) d u d τ
Substitute t = τ + u t = \tau + u t = τ + u , u = t − τ u = t - \tau u = t − τ . For fixed τ \tau τ , u u u ranges from 0 0 0 to ∞ \infty ∞ So t t t ranges From τ \tau τ to ∞ \infty ∞ .
= ∫ 0 ∞ f ( τ ) [ ∫ τ ∞ e − s t g ( t − τ ) d t ] d τ = \int_0^{\infty} f(\tau)\left[\int_{\tau}^{\infty} e^{-st} g(t - \tau)\, dt\right] d\tau = ∫ 0 ∞ f ( τ ) [ ∫ τ ∞ e − s t g ( t − τ ) d t ] d τ
Swap the order of integration (Fubini):
= ∫ 0 ∞ e − s t [ ∫ 0 t f ( τ ) g ( t − τ ) d τ ] d t = L { ( f ∗ g ) ( t ) } = \int_0^{\infty} e^{-st}\left[\int_0^t f(\tau) g(t - \tau)\, d\tau\right] dt = \mathcal{L}\{(f * g)(t)\} = ∫ 0 ∞ e − s t [ ∫ 0 t f ( τ ) g ( t − τ ) d τ ] d t = L {( f ∗ g ) ( t )} . ■ \blacksquare ■
Problem. Find L − 1 { 3 s + 1 s 2 + 4 s + 13 } \mathcal{L}^{-1}\left\{\frac{3s + 1}{s^2 + 4s + 13}\right\} L − 1 { s 2 + 4 s + 13 3 s + 1 } .
Solution Solution. Complete the square: s 2 + 4 s + 13 = ( s + 2 ) 2 + 9 s^2 + 4s + 13 = (s + 2)^2 + 9 s 2 + 4 s + 13 = ( s + 2 ) 2 + 9 .
3 s + 1 s 2 + 4 s + 13 = 3 ( s + 2 ) − 5 ( s + 2 ) 2 + 9 = 3 ⋅ s + 2 ( s + 2 ) 2 + 9 − 5 3 ⋅ 3 ( s + 2 ) 2 + 9 \frac{3s + 1}{s^2 + 4s + 13} = \frac{3(s + 2) - 5}{(s + 2)^2 + 9} = 3 \cdot \frac{s + 2}{(s + 2)^2 + 9} - \frac{5}{3} \cdot \frac{3}{(s + 2)^2 + 9} s 2 + 4 s + 13 3 s + 1 = ( s + 2 ) 2 + 9 3 ( s + 2 ) − 5 = 3 ⋅ ( s + 2 ) 2 + 9 s + 2 − 3 5 ⋅ ( s + 2 ) 2 + 9 3
Using the shifting theorem:
L − 1 { s + 2 ( s + 2 ) 2 + 9 } = e − 2 t cos ( 3 t ) \mathcal{L}^{-1}\left\{\frac{s + 2}{(s + 2)^2 + 9}\right\} = e^{-2t}\cos(3t) L − 1 { ( s + 2 ) 2 + 9 s + 2 } = e − 2 t cos ( 3 t )
L − 1 { 3 ( s + 2 ) 2 + 9 } = e − 2 t sin ( 3 t ) \mathcal{L}^{-1}\left\{\frac{3}{(s + 2)^2 + 9}\right\} = e^{-2t}\sin(3t) L − 1 { ( s + 2 ) 2 + 9 3 } = e − 2 t sin ( 3 t )
f ( t ) = 3 e − 2 t cos ( 3 t ) − 5 3 e − 2 t sin ( 3 t ) = e − 2 t ( 3 cos ( 3 t ) − 5 3 sin ( 3 t ) ) f(t) = 3e^{-2t}\cos(3t) - \frac{5}{3}e^{-2t}\sin(3t) = e^{-2t}\left(3\cos(3t) - \frac{5}{3}\sin(3t)\right) f ( t ) = 3 e − 2 t cos ( 3 t ) − 3 5 e − 2 t sin ( 3 t ) = e − 2 t ( 3 cos ( 3 t ) − 3 5 sin ( 3 t ) ) . ■ \blacksquare ■
Problem. Solve y ′ ′ + 3 y ′ + 2 y = u 2 ( t ) y'' + 3y' + 2y = u_2(t) y ′′ + 3 y ′ + 2 y = u 2 ( t ) , y ( 0 ) = 0 y(0) = 0 y ( 0 ) = 0 , y ′ ( 0 ) = 0 y'(0) = 0 y ′ ( 0 ) = 0 .
Solution Solution. Take Laplace transforms:
s 2 Y + 3 s Y + 2 Y = e − 2 s s s^2 Y + 3sY + 2Y = \frac{e^{-2s}}{s} s 2 Y + 3 s Y + 2 Y = s e − 2 s
( s 2 + 3 s + 2 ) Y = e − 2 s s (s^2 + 3s + 2)Y = \frac{e^{-2s}}{s} ( s 2 + 3 s + 2 ) Y = s e − 2 s
Y = e − 2 s s ( s + 1 ) ( s + 2 ) Y = \frac{e^{-2s}}{s(s+1)(s+2)} Y = s ( s + 1 ) ( s + 2 ) e − 2 s
First find the partial fraction decomposition of 1 s ( s + 1 ) ( s + 2 ) \frac{1}{s(s+1)(s+2)} s ( s + 1 ) ( s + 2 ) 1 :
1 s ( s + 1 ) ( s + 2 ) = A s + B s + 1 + C s + 2 \frac{1}{s(s+1)(s+2)} = \frac{A}{s} + \frac{B}{s+1} + \frac{C}{s+2} s ( s + 1 ) ( s + 2 ) 1 = s A + s + 1 B + s + 2 C
1 = A ( s + 1 ) ( s + 2 ) + B s ( s + 2 ) + C s ( s + 1 ) 1 = A(s+1)(s+2) + Bs(s+2) + Cs(s+1) 1 = A ( s + 1 ) ( s + 2 ) + B s ( s + 2 ) + C s ( s + 1 )
s = 0 s = 0 s = 0 : 1 = 2 A ⟹ A = 1 / 2 1 = 2A \implies A = 1/2 1 = 2 A ⟹ A = 1/2 s = − 1 s = -1 s = − 1 : 1 = − B ⟹ B = − 1 1 = -B \implies B = -1 1 = − B ⟹ B = − 1 s = − 2 s = -2 s = − 2 : 1 = 2 C ⟹ C = 1 / 2 1 = 2C \implies C = 1/2 1 = 2 C ⟹ C = 1/2
So L − 1 { 1 s ( s + 1 ) ( s + 2 ) } = 1 2 − e − t + 1 2 e − 2 t \mathcal{L}^{-1}\left\{\frac{1}{s(s+1)(s+2)}\right\} = \frac{1}{2} - e^{-t} + \frac{1}{2}e^{-2t} L − 1 { s ( s + 1 ) ( s + 2 ) 1 } = 2 1 − e − t + 2 1 e − 2 t .
By the second shifting theorem:
y ( t ) = u 2 ( t ) [ 1 2 − e − ( t − 2 ) + 1 2 e − 2 ( t − 2 ) ] y(t) = u_2(t)\left[\frac{1}{2} - e^{-(t-2)} + \frac{1}{2}e^{-2(t-2)}\right] y ( t ) = u 2 ( t ) [ 2 1 − e − ( t − 2 ) + 2 1 e − 2 ( t − 2 ) ] for t ≥ 0 t \geq 0 t ≥ 0 . ■ \blacksquare ■
The Heaviside (unit step) function is defined as
u c ( t ) = { 0 t < c 1 t ≥ c u_c(t) = \begin{cases} 0 & t \lt c \\ 1 & t \geq c \end{cases} u c ( t ) = { 0 1 t < c t ≥ c
It models a sudden switch being turned on at time t = c t = c t = c . The second shifting theorem states:
L { u c ( t ) f ( t − c ) } = e − c s F ( s ) \mathcal{L}\{u_c(t)f(t - c)\} = e^{-cs}F(s) L { u c ( t ) f ( t − c )} = e − cs F ( s )
Conversely, if Y ( s ) = e − c s G ( s ) Y(s) = e^{-cs}G(s) Y ( s ) = e − cs G ( s ) Then y ( t ) = u c ( t ) ⋅ g ( t − c ) y(t) = u_c(t) \cdot g(t - c) y ( t ) = u c ( t ) ⋅ g ( t − c ) where g = L − 1 { G } g = \mathcal{L}^{-1}\{G\} g = L − 1 { G } .
The Dirac delta function δ ( t − c ) \delta(t - c) δ ( t − c ) satisfies ∫ 0 ∞ δ ( t − c ) f ( t ) d t = f ( c ) \int_0^{\infty} \delta(t - c)f(t)\, dt = f(c) ∫ 0 ∞ δ ( t − c ) f ( t ) d t = f ( c ) for continuous f f f . It models an instantaneous Impulse.
L { δ ( t − c ) } = e − c s \mathcal{L}\{\delta(t - c)\} = e^{-cs} L { δ ( t − c )} = e − cs .
Properties:
δ ( t − c ) = 0 \delta(t - c) = 0 δ ( t − c ) = 0 for t ≠ c t \neq c t = c ∫ − ∞ ∞ δ ( t − c ) d t = 1 \int_{-\infty}^{\infty} \delta(t - c)\, dt = 1 ∫ − ∞ ∞ δ ( t − c ) d t = 1 ∫ − ∞ ∞ δ ( t − c ) f ( t ) d t = f ( c ) \int_{-\infty}^{\infty} \delta(t - c)f(t)\, dt = f(c) ∫ − ∞ ∞ δ ( t − c ) f ( t ) d t = f ( c ) (sifting property)Problem. An undamped spring-mass system (m = 1 m = 1 m = 1 , k = 4 k = 4 k = 4 ) is struck by an impulse at t = 0 t = 0 t = 0 : y ′ ′ + 4 y = δ ( t ) y'' + 4y = \delta(t) y ′′ + 4 y = δ ( t ) , y ( 0 ) = 0 y(0) = 0 y ( 0 ) = 0 , y ′ ( 0 ) = 0 y'(0) = 0 y ′ ( 0 ) = 0 . Find y ( t ) y(t) y ( t ) .
Solution Solution. Taking Laplace transforms:
s 2 Y + 4 Y = 1 s^2 Y + 4Y = 1 s 2 Y + 4 Y = 1
Y = 1 s 2 + 4 Y = \frac{1}{s^2 + 4} Y = s 2 + 4 1
y ( t ) = 1 2 sin ( 2 t ) y(t) = \frac{1}{2}\sin(2t) y ( t ) = 2 1 sin ( 2 t ) .
This is the impulse response (Green’s function) for the operator D 2 + 4 D^2 + 4 D 2 + 4 . ■ \blacksquare ■
Problem. Solve y ′ ′ − 6 y ′ + 9 y = 6 e 3 t y'' - 6y' + 9y = 6e^{3t} y ′′ − 6 y ′ + 9 y = 6 e 3 t , y ( 0 ) = 0 y(0) = 0 y ( 0 ) = 0 , y ′ ( 0 ) = 0 y'(0) = 0 y ′ ( 0 ) = 0 .
Solution Solution. Taking Laplace transforms:
s 2 Y − 6 s Y + 9 Y = 6 s − 3 s^2 Y - 6sY + 9Y = \frac{6}{s - 3} s 2 Y − 6 s Y + 9 Y = s − 3 6
( s − 3 ) 2 Y = 6 s − 3 (s - 3)^2 Y = \frac{6}{s - 3} ( s − 3 ) 2 Y = s − 3 6
Y = 6 ( s − 3 ) 3 Y = \frac{6}{(s - 3)^3} Y = ( s − 3 ) 3 6
Using the table: L { t n e a t } = n ! ( s − a ) n + 1 \mathcal{L}\{t^n e^{at}\} = \frac{n!}{(s-a)^{n+1}} L { t n e a t } = ( s − a ) n + 1 n ! .
y ( t ) = 6 ⋅ t 2 e 3 t 2 ! = 3 t 2 e 3 t y(t) = 6 \cdot \frac{t^2 e^{3t}}{2!} = 3t^2 e^{3t} y ( t ) = 6 ⋅ 2 ! t 2 e 3 t = 3 t 2 e 3 t . ■ \blacksquare ■
The Laplace transform converts differential equations into algebraic equations. Instead of solving a second-order ODE (which requires finding homogeneous and particular solutions, applying initial conditions, etc.), you transform the ODE into a polynomial equation for Y ( s ) Y(s) Y ( s ) , solve for Y ( s ) Y(s) Y ( s ) , and then invert.
The key insight: The derivative property L { f ′ ( t ) } = s F ( s ) − f ( 0 ) \mathcal{L}\{f'(t)\} = sF(s) - f(0) L { f ′ ( t )} = s F ( s ) − f ( 0 ) converts differentiation into multiplication by s s s . This means a differential equation like y ′ ′ + 3 y ′ + 2 y = e t y'' + 3y' + 2y = e^t y ′′ + 3 y ′ + 2 y = e t becomes ( s 2 + 3 s + 2 ) Y ( s ) = 1 s − 1 + initial condition terms (s^2 + 3s + 2)Y(s) = \frac{1}{s-1} + \text{initial condition terms} ( s 2 + 3 s + 2 ) Y ( s ) = s − 1 1 + initial condition terms , which is just algebra.
Why the convolution theorem matters: The convolution ( f ∗ g ) ( t ) = ∫ 0 t f ( τ ) g ( t − τ ) d τ (f * g)(t) = \int_0^t f(\tau)g(t - \tau)\,d\tau ( f ∗ g ) ( t ) = ∫ 0 t f ( τ ) g ( t − τ ) d τ is hard to compute directly. But L { f ∗ g } = F ( s ) G ( s ) \mathcal{L}\{f * g\} = F(s)G(s) L { f ∗ g } = F ( s ) G ( s ) means that convolution in the time domain becomes multiplication in the frequency domain. This is the same principle behind Fourier transforms and signal processing.
Physical interpretation: The Laplace variable s s s can be interpreted as a complex frequency. F ( s ) F(s) F ( s ) decomposes the signal f ( t ) f(t) f ( t ) into exponential components e s t e^{st} e s t . The poles of F ( s ) F(s) F ( s ) (values of s s s where F ( s ) F(s) F ( s ) blows up) determine the natural modes of the system, while the zeros determine the frequencies that are suppressed.
Problem. Solve the system: x ′ = 3 x − 2 y y ′ = 2 x − 2 y \begin{aligned} x' &= 3x - 2y \\ y' &= 2x - 2y \end{aligned} x ′ y ′ = 3 x − 2 y = 2 x − 2 y with x ( 0 ) = 1 x(0) = 1 x ( 0 ) = 1 , y ( 0 ) = 0 y(0) = 0 y ( 0 ) = 0 .
Solution Take Laplace transforms: s X − 1 = 3 X − 2 Y s Y = 2 X − 2 Y \begin{aligned} sX - 1 &= 3X - 2Y \\ sY &= 2X - 2Y \end{aligned} s X − 1 s Y = 3 X − 2 Y = 2 X − 2 Y
Rearranging: ( s − 3 ) X + 2 Y = 1 − 2 X + ( s + 2 ) Y = 0 \begin{aligned} (s - 3)X + 2Y &= 1 \\ -2X + (s + 2)Y &= 0 \end{aligned} ( s − 3 ) X + 2 Y − 2 X + ( s + 2 ) Y = 1 = 0
From the second equation: X = ( s + 2 ) Y 2 X = \frac{(s+2)Y}{2} X = 2 ( s + 2 ) Y .
Substituting into the first: ( s − 3 ) ⋅ ( s + 2 ) Y 2 + 2 Y = 1 (s - 3) \cdot \frac{(s+2)Y}{2} + 2Y = 1 ( s − 3 ) ⋅ 2 ( s + 2 ) Y + 2 Y = 1
( s − 3 ) ( s + 2 ) Y + 4 Y 2 = 1 \frac{(s-3)(s+2)Y + 4Y}{2} = 1 2 ( s − 3 ) ( s + 2 ) Y + 4 Y = 1
( s 2 − s − 6 + 4 ) Y = 2 (s^2 - s - 6 + 4)Y = 2 ( s 2 − s − 6 + 4 ) Y = 2
( s 2 − s − 2 ) Y = 2 (s^2 - s - 2)Y = 2 ( s 2 − s − 2 ) Y = 2
Y = 2 ( s − 2 ) ( s + 1 ) = 2 / 3 s − 2 − 2 / 3 s + 1 Y = \frac{2}{(s-2)(s+1)} = \frac{2/3}{s-2} - \frac{2/3}{s+1} Y = ( s − 2 ) ( s + 1 ) 2 = s − 2 2/3 − s + 1 2/3
y ( t ) = 2 3 e 2 t − 2 3 e − t y(t) = \frac{2}{3}e^{2t} - \frac{2}{3}e^{-t} y ( t ) = 3 2 e 2 t − 3 2 e − t
For X X X : X = ( s + 2 ) 2 ⋅ 2 ( s − 2 ) ( s + 1 ) = s + 2 ( s − 2 ) ( s + 1 ) = 4 / 3 s − 2 − 1 / 3 s + 1 X = \frac{(s+2)}{2} \cdot \frac{2}{(s-2)(s+1)} = \frac{s+2}{(s-2)(s+1)} = \frac{4/3}{s-2} - \frac{1/3}{s+1} X = 2 ( s + 2 ) ⋅ ( s − 2 ) ( s + 1 ) 2 = ( s − 2 ) ( s + 1 ) s + 2 = s − 2 4/3 − s + 1 1/3
x ( t ) = 4 3 e 2 t − 1 3 e − t x(t) = \frac{4}{3}e^{2t} - \frac{1}{3}e^{-t} x ( t ) = 3 4 e 2 t − 3 1 e − t
Verification: x ( 0 ) = 4 / 3 − 1 / 3 = 1 x(0) = 4/3 - 1/3 = 1 x ( 0 ) = 4/3 − 1/3 = 1 ✓ \checkmark ✓ , y ( 0 ) = 2 / 3 − 2 / 3 = 0 y(0) = 2/3 - 2/3 = 0 y ( 0 ) = 2/3 − 2/3 = 0 ✓ \checkmark ✓ .
x ′ = 8 3 e 2 t + 1 3 e − t x' = \frac{8}{3}e^{2t} + \frac{1}{3}e^{-t} x ′ = 3 8 e 2 t + 3 1 e − t , 3 x − 2 y = 4 e 2 t − e − t − 4 3 e 2 t + 4 3 e − t = 8 3 e 2 t + 1 3 e − t 3x - 2y = 4e^{2t} - e^{-t} - \frac{4}{3}e^{2t} + \frac{4}{3}e^{-t} = \frac{8}{3}e^{2t} + \frac{1}{3}e^{-t} 3 x − 2 y = 4 e 2 t − e − t − 3 4 e 2 t + 3 4 e − t = 3 8 e 2 t + 3 1 e − t ✓ \checkmark ✓ . ■ \blacksquare ■
The Laplace transform converts differential equations into algebraic equations. By transforming a function f ( t ) f(t) f ( t ) into F ( s ) = ∫ 0 ∞ e − s t f ( t ) d t F(s) = \int_0^\infty e^{-st}f(t)\,dt F ( s ) = ∫ 0 ∞ e − s t f ( t ) d t , differentiation in the time domain becomes multiplication by s s s in the frequency domain. This means a linear ODE with constant coefficients becomes a polynomial equation in s s s , which can be solved algebraically and then inverted to recover the solution.
The Laplace transform is particularly powerful for several reasons. It handles discontinuous forcing functions (step functions, impulses) logically through the second shifting theorem and the convolution theorem. It converts initial value problems into algebraic problems without needing to find the homogeneous solution first. The convolution theorem says that convolution in time becomes multiplication in the frequency domain, which is why the Laplace transform is the tool of choice for systems theory, control engineering, and circuit analysis. The transform exists for functions of exponential order, which covers essentially all physically reasonable inputs.
The Picard-Lindelöf theorem can be proved constructively via Picard iteration . For the IVP y ′ = f ( x , y ) y' = f(x, y) y ′ = f ( x , y ) , y ( x 0 ) = y 0 y(x_0) = y_0 y ( x 0 ) = y 0 Define the sequence
ϕ 0 ( x ) = y 0 , ϕ n + 1 ( x ) = y 0 + ∫ x 0 x f ( t , ϕ n ( t ) ) d t \phi_0(x) = y_0, \quad \phi_{n+1}(x) = y_0 + \int_{x_0}^x f(t, \phi_n(t))\, dt ϕ 0 ( x ) = y 0 , ϕ n + 1 ( x ) = y 0 + ∫ x 0 x f ( t , ϕ n ( t )) d t
If f f f and ∂ f / ∂ y \partial f/\partial y ∂ f / ∂ y are continuous, one shows by induction that ( ϕ n ) (\phi_n) ( ϕ n ) is uniformly Cauchy on some interval [ x 0 − h , x 0 + h ] [x_0 - h, x_0 + h] [ x 0 − h , x 0 + h ] Hence converges uniformly to a function ϕ \phi ϕ . Passing to the limit in the integral equation shows ϕ \phi ϕ satisfies the ODE. Uniqueness follows From the Gronwall inequality applied to the difference of two solutions.
A[5_Laplace Transforms] --> B[Key Concepts]
A --> D[Practical Applications]
B --> E[Fundamental definitions]
D --> G[Real-world usage]
Second-Order Linear ODEs : The Laplace transform converts the characteristic equation approach into algebraic manipulation in the frequency domain.
Fourier Series : Both Fourier and Laplace transforms decompose signals into frequency components, but Laplace handles transient and aperiodic signals.
Complex Numbers Review : The Laplace variable s s s is complex, and poles and zeros in the complex plane determine system behaviour.
Systems of ODEs : The Laplace transform extends logically to systems of ODEs, converting coupled differential equations into algebraic systems.
Classical Mechanics
Electromagnetism