5.1 Definition
The Laplace transform of f(t) (defined for t≥0) is
L{f(t)}=F(s)=∫0∞e−stf(t)dt
The transform exists when f is piecewise continuous on [0,∞) and of exponential order: ∣f(t)∣≤Meat for some M,a>0.
5.2 Basic Properties
Theorem 5.1 (Linearity). L{af+bg}=aL{f}+bL{g}.
Theorem 5.2 (First Shifting). L{eatf(t)}=F(s−a).
Theorem 5.3 (Second Shifting). L{uc(t)f(t−c)}=e−csF(s)Where uc(t) is The unit step function.
Theorem 5.4 (Derivative). L{f"(t)}=sF(s)−f(0).
Theorem 5.5 (n-th Derivative). L{f(n)(t)}=snF(s)−sn−1f(0)−⋯−f(n−1)(0).
Theorem 5.6 (Integration). L{∫0tf(τ)dτ}=sF(s).
Theorem 5.7 (Convolution). L{f∗g}=F(s)G(s)Where (f∗g)(t)=∫0tf(τ)g(t−τ)dτ.
5.3 Proofs of Key Properties
Proof of Linearity. L{af+bg}=∫0∞e−st[af(t)+bg(t)]dt=a∫0∞e−stf(t)dt+b∫0∞e−stg(t)dt=aF(s)+bG(s). ■
Proof of First Shifting Theorem. L{eatf(t)}=∫0∞e−steatf(t)dt=∫0∞e−(s−a)tf(t)dt=F(s−a). ■
Proof of Derivative Property. Integrating by parts with u=e−st, dv=f′(t)dt:
L{f′(t)}=[e−stf(t)]0∞+s∫0∞e−stf(t)dt=−f(0)+sF(s)=sF(s)−f(0).
The boundary term vanishes at ∞ because f is of exponential order. ■
Proof of Integration Property. Let g(t)=∫0tf(τ)dτ. Then g′(t)=f(t) and g(0)=0. L{f(t)}=L{g′(t)}=sL{g(t)}−g(0)=sL{g(t)}. Therefore L{g(t)}=F(s)/s. ■
| f(t) | F(s)=L{f(t)} |
|---|
| 1 | 1/s |
| tn | n!/sn+1 |
| eat | 1/(s−a) |
| tneat | n!/(s−a)n+1 |
| sin(bt) | b/(s2+b2) |
| cos(bt) | s/(s2+b2) |
| eatsin(bt) | b/((s−a)2+b2) |
| eatcos(bt) | (s−a)/((s−a)2+b2) |
| uc(t) | e−cs/s |
| δ(t−c) | e−cs |
Procedure:
- Take the Laplace transform of both sides of the ODE.
- Solve the resulting algebraic equation for Y(s)=L{y(t)}.
- Find the inverse Laplace transform to obtain y(t).
5.6 Worked Example: Solving an IVP
Problem. Solve y′′+4y=sin(2t), y(0)=0, y′(0)=0.
Solution. Taking Laplace transforms:
s2Y+4Y=s2+42
(s2+4)Y=s2+42
Y=(s2+4)22
To invert, use the convolution theorem: L−1{s2+41}=21sin(2t).
y(t)=21sin(2t)∗21sin(2t)=41∫0tsin(2τ)sin(2(t−τ))dτ
Using sinAsinB=21(cos(A−B)−cos(A+B)):
y(t)=81∫0t[cos(4τ−2t)−cos(2t)]dτ=81[4sin(4τ−2t)]0t−8tcos(2t)
=321[sin(2t)−sin(−2t)]−8tcos(2t)=16sin(2t)−8tcos(2t)
y(t)=16sin(2t)−2tcos(2t). ■
5.6b Worked Example: Partial Fractions for Inverse Laplace
Problem. Find \mathcal{L}^{-1}\left\{\frac{1}{s(s^2 + 4)} Ight\}.
Solution
Solution. Partial fractions: s(s2+4)1=sA+s2+4Bs+C.
1=A(s2+4)+(Bs+C)s=(A+B)s2+Cs+4A.
s2: A+B=0, s: C=0Constant: 4A=1⟹A=1/4, B=−1/4.
s(s2+4)1=s1/4−s2+4s/4.
f(t)=41−41cos2t. ■
### 5.7 Convolution Theorem
Theorem 5.8 (Convolution). If F(s)=L{f(t)} and G(s)=L{g(t)} Then
L{f∗g}=F(s)G(s)
Where (f∗g)(t)=∫0tf(τ)g(t−τ)dτ.
Proof. F(s)G(s)=∫0∞e−sτf(τ)dτ⋅∫0∞e−sug(u)du
=∫0∞∫0∞e−s(τ+u)f(τ)g(u)dudτ
Substitute t=τ+u, u=t−τ. For fixed τ, u ranges from 0 to ∞So t ranges From τ to ∞.
=∫0∞f(τ)[∫τ∞e−stg(t−τ)dt]dτ
Swap the order of integration (Fubini):
=∫0∞e−st[∫0tf(τ)g(t−τ)dτ]dt=L{(f∗g)(t)}. ■
Problem. Find L−1{s2+4s+133s+1}.
Solution
Solution. Complete the square: s2+4s+13=(s+2)2+9.
s2+4s+133s+1=(s+2)2+93(s+2)−5=3⋅(s+2)2+9s+2−35⋅(s+2)2+93
Using the shifting theorem:
L−1{(s+2)2+9s+2}=e−2tcos(3t)
L−1{(s+2)2+93}=e−2tsin(3t)
f(t)=3e−2tcos(3t)−35e−2tsin(3t)=e−2t(3cos(3t)−35sin(3t)). ■
5.9 Worked Example: IVP with Discontinuous Forcing
Problem. Solve y′′+3y′+2y=u2(t), y(0)=0, y′(0)=0.
Solution
Solution. Take Laplace transforms:
s2Y+3sY+2Y=se−2s
(s2+3s+2)Y=se−2s
Y=s(s+1)(s+2)e−2s
First find the partial fraction decomposition of s(s+1)(s+2)1:
s(s+1)(s+2)1=sA+s+1B+s+2C
1=A(s+1)(s+2)+Bs(s+2)+Cs(s+1)
s=0: 1=2A⟹A=1/2 s=−1: 1=−B⟹B=−1 s=−2: 1=2C⟹C=1/2
So L−1{s(s+1)(s+2)1}=21−e−t+21e−2t.
By the second shifting theorem:
y(t)=u2(t)[21−e−(t−2)+21e−2(t−2)] for t≥0. ■
5.10 The Heaviside Function
The Heaviside (unit step) function is defined as
uc(t)={01t<ct≥c
It models a sudden switch being turned on at time t=c. The second shifting theorem states:
L{uc(t)f(t−c)}=e−csF(s)
Conversely, if Y(s)=e−csG(s)Then y(t)=uc(t)⋅g(t−c) where g=L−1{G}.
5.11 The Dirac Delta Function
The Dirac delta function δ(t−c) satisfies ∫0∞δ(t−c)f(t)dt=f(c) for continuous f. It models an instantaneous Impulse.
L{δ(t−c)}=e−cs.
Properties:
- δ(t−c)=0 for t=c
- ∫−∞∞δ(t−c)dt=1
- ∫−∞∞δ(t−c)f(t)dt=f(c) (sifting property)
5.12 Worked Example: Impulse Response
Problem. An undamped spring-mass system (m=1, k=4) is struck by an impulse at t=0: y′′+4y=δ(t), y(0)=0, y′(0)=0. Find y(t).
Solution
Solution. Taking Laplace transforms:
s2Y+4Y=1
Y=s2+41
y(t)=21sin(2t).
This is the impulse response (Green’s function) for the operator D2+4. ■
5.13 Worked Example: IVP with Laplace
Problem. Solve y′′−6y′+9y=6e3t, y(0)=0, y′(0)=0.
Solution
Solution. Taking Laplace transforms:
s2Y−6sY+9Y=s−36
(s−3)2Y=s−36
Y=(s−3)36
Using the table: L{tneat}=(s−a)n+1n!.
y(t)=6⋅2!t2e3t=3t2e3t. ■
:::caution Common Pitfall The Laplace transform of y′(t) is sY(s)−y(0)Not sY(s). The Initial conditions are built into the transform. Forgetting them leads to incorrect solutions. :::
:::caution Common Pitfall When using the second shifting theorem, the time-shifted function must Be written as uc(t)f(t−c)Not uc(t)f(t). The function f must be shifted by the same Amount as the step. :::
5.15 Proof Sketch: Picard Iteration
The Picard-Lindelöf theorem can be proved constructively via Picard iteration. For the IVP y′=f(x,y), y(x0)=y0Define the sequence
ϕ0(x)=y0,ϕn+1(x)=y0+∫x0xf(t,ϕn(t))dt
If f and ∂f/∂y are continuous, one shows by induction that (ϕn) is uniformly Cauchy on some interval [x0−h,x0+h]Hence converges uniformly to a function ϕ. Passing to the limit in the integral equation shows ϕ satisfies the ODE. Uniqueness follows From the Gronwall inequality applied to the difference of two solutions.