Skip to content

Laplace Transforms

5.1 Definition

The Laplace transform of f(t)f(t) (defined for t0t \geq 0) is

L{f(t)}=F(s)=0estf(t)dt\mathcal{L}\{f(t)\} = F(s) = \int_0^{\infty} e^{-st} f(t)\, dt

The transform exists when ff is piecewise continuous on [0,)[0, \infty) and of exponential order: f(t)Meat|f(t)| \leq Me^{at} for some M,a>0M, a > 0.

5.2 Basic Properties

Theorem 5.1 (Linearity). L{af+bg}=aL{f}+bL{g}\mathcal{L}\{af + bg\} = a\mathcal{L}\{f\} + b\mathcal{L}\{g\}.

Theorem 5.2 (First Shifting). L{eatf(t)}=F(sa)\mathcal{L}\{e^{at}f(t)\} = F(s - a).

Theorem 5.3 (Second Shifting). L{uc(t)f(tc)}=ecsF(s)\mathcal{L}\{u_c(t)f(t - c)\} = e^{-cs}F(s)Where uc(t)u_c(t) is The unit step function.

Theorem 5.4 (Derivative). L{f"(t)}=sF(s)f(0)\mathcal{L}\{f"(t)\} = sF(s) - f(0).

Theorem 5.5 (n-th Derivative). L{f(n)(t)}=snF(s)sn1f(0)f(n1)(0)\mathcal{L}\{f^{(n)}(t)\} = s^n F(s) - s^{n-1}f(0) - \cdots - f^{(n-1)}(0).

Theorem 5.6 (Integration). L{0tf(τ)dτ}=F(s)s\mathcal{L}\left\{\int_0^t f(\tau)\, d\tau\right\} = \frac{F(s)}{s}.

Theorem 5.7 (Convolution). L{fg}=F(s)G(s)\mathcal{L}\{f * g\} = F(s)G(s)Where (fg)(t)=0tf(τ)g(tτ)dτ(f * g)(t) = \int_0^t f(\tau)g(t - \tau)\, d\tau.

5.3 Proofs of Key Properties

Proof of Linearity. L{af+bg}=0est[af(t)+bg(t)]dt=a0estf(t)dt+b0estg(t)dt=aF(s)+bG(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). \blacksquare

Proof of First Shifting Theorem. L{eatf(t)}=0esteatf(t)dt=0e(sa)tf(t)dt=F(sa)\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). \blacksquare

Proof of Derivative Property. Integrating by parts with u=estu = e^{-st}, dv=f(t)dtdv = f'(t)\, dt:

L{f(t)}=[estf(t)]0+s0estf(t)dt=f(0)+sF(s)=sF(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).

The boundary term vanishes at \infty because ff is of exponential order. \blacksquare

Proof of Integration Property. Let g(t)=0tf(τ)dτg(t) = \int_0^t f(\tau)\, d\tau. Then g(t)=f(t)g'(t) = f(t) and g(0)=0g(0) = 0. L{f(t)}=L{g(t)}=sL{g(t)}g(0)=sL{g(t)}\mathcal{L}\{f(t)\} = \mathcal{L}\{g'(t)\} = s\mathcal{L}\{g(t)\} - g(0) = s\mathcal{L}\{g(t)\}. Therefore L{g(t)}=F(s)/s\mathcal{L}\{g(t)\} = F(s)/s. \blacksquare

5.4 Common Laplace Transforms

f(t)f(t)F(s)=L{f(t)}F(s) = \mathcal{L}\{f(t)\}
111/s1/s
tnt^nn!/sn+1n!/s^{n+1}
eate^{at}1/(sa)1/(s - a)
tneatt^n e^{at}n!/(sa)n+1n!/(s - a)^{n+1}
sin(bt)\sin(bt)b/(s2+b2)b/(s^2 + b^2)
cos(bt)\cos(bt)s/(s2+b2)s/(s^2 + b^2)
eatsin(bt)e^{at}\sin(bt)b/((sa)2+b2)b/((s-a)^2 + b^2)
eatcos(bt)e^{at}\cos(bt)(sa)/((sa)2+b2)(s-a)/((s-a)^2 + b^2)
uc(t)u_c(t)ecs/se^{-cs}/s
δ(tc)\delta(t - c)ecse^{-cs}

5.5 Solving IVPs with Laplace Transforms

Procedure:

  1. Take the Laplace transform of both sides of the ODE.
  2. Solve the resulting algebraic equation for Y(s)=L{y(t)}Y(s) = \mathcal{L}\{y(t)\}.
  3. Find the inverse Laplace transform to obtain y(t)y(t).

5.6 Worked Example: Solving an IVP

Problem. Solve y+4y=sin(2t)y'' + 4y = \sin(2t), y(0)=0y(0) = 0, y(0)=0y'(0) = 0.

Solution. Taking Laplace transforms:

s2Y+4Y=2s2+4s^2 Y + 4Y = \frac{2}{s^2 + 4}

(s2+4)Y=2s2+4(s^2 + 4)Y = \frac{2}{s^2 + 4}

Y=2(s2+4)2Y = \frac{2}{(s^2 + 4)^2}

To invert, use the convolution theorem: L1{1s2+4}=12sin(2t)\mathcal{L}^{-1}\left\{\frac{1}{s^2 + 4}\right\} = \frac{1}{2}\sin(2t).

y(t)=12sin(2t)12sin(2t)=140tsin(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

Using sinAsinB=12(cos(AB)cos(A+B))\sin A \sin B = \frac{1}{2}(\cos(A - B) - \cos(A + B)):

y(t)=180t[cos(4τ2t)cos(2t)]dτ=18[sin(4τ2t)4]0tt8cos(2t)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)

=132[sin(2t)sin(2t)]t8cos(2t)=sin(2t)16tcos(2t)8= \frac{1}{32}[\sin(2t) - \sin(-2t)] - \frac{t}{8}\cos(2t) = \frac{\sin(2t)}{16} - \frac{t\cos(2t)}{8}

y(t)=sin(2t)2tcos(2t)16y(t) = \frac{\sin(2t) - 2t\cos(2t)}{16}. \blacksquare

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: 1s(s2+4)=As+Bs+Cs2+4\frac{1}{s(s^2 + 4)} = \frac{A}{s} + \frac{Bs + C}{s^2 + 4}.

1=A(s2+4)+(Bs+C)s=(A+B)s2+Cs+4A1 = A(s^2 + 4) + (Bs + C)s = (A + B)s^2 + Cs + 4A.

s2s^2: A+B=0A + B = 0, ss: C=0C = 0Constant: 4A=1    A=1/44A = 1 \implies A = 1/4, B=1/4B = -1/4.

1s(s2+4)=1/4ss/4s2+4\frac{1}{s(s^2 + 4)} = \frac{1/4}{s} - \frac{s/4}{s^2 + 4}.

f(t)=1414cos2tf(t) = \frac{1}{4} - \frac{1}{4}\cos 2t. \blacksquare

### 5.7 Convolution Theorem

Theorem 5.8 (Convolution). If F(s)=L{f(t)}F(s) = \mathcal{L}\{f(t)\} and G(s)=L{g(t)}G(s) = \mathcal{L}\{g(t)\} Then

L{fg}=F(s)G(s)\mathcal{L}\{f * g\} = F(s)G(s)

Where (fg)(t)=0tf(τ)g(tτ)dτ(f * g)(t) = \int_0^t f(\tau)g(t - \tau)\, d\tau.

Proof. F(s)G(s)=0esτf(τ)dτ0esug(u)duF(s)G(s) = \int_0^{\infty} e^{-s\tau} f(\tau)\, d\tau \cdot \int_0^{\infty} e^{-su} g(u)\, du

=00es(τ+u)f(τ)g(u)dudτ= \int_0^{\infty} \int_0^{\infty} e^{-s(\tau+u)} f(\tau) g(u)\, du\, d\tau

Substitute t=τ+ut = \tau + u, u=tτu = t - \tau. For fixed τ\tau, uu ranges from 00 to \inftySo tt ranges From τ\tau to \infty.

=0f(τ)[τestg(tτ)dt]dτ= \int_0^{\infty} f(\tau)\left[\int_{\tau}^{\infty} e^{-st} g(t - \tau)\, dt\right] d\tau

Swap the order of integration (Fubini):

=0est[0tf(τ)g(tτ)dτ]dt=L{(fg)(t)}= \int_0^{\infty} e^{-st}\left[\int_0^t f(\tau) g(t - \tau)\, d\tau\right] dt = \mathcal{L}\{(f * g)(t)\}. \blacksquare

5.8 Worked Example: Inverse Laplace Transform

Problem. Find L1{3s+1s2+4s+13}\mathcal{L}^{-1}\left\{\frac{3s + 1}{s^2 + 4s + 13}\right\}.

Solution

Solution. Complete the square: s2+4s+13=(s+2)2+9s^2 + 4s + 13 = (s + 2)^2 + 9.

3s+1s2+4s+13=3(s+2)5(s+2)2+9=3s+2(s+2)2+9533(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}

Using the shifting theorem:

L1{s+2(s+2)2+9}=e2tcos(3t)\mathcal{L}^{-1}\left\{\frac{s + 2}{(s + 2)^2 + 9}\right\} = e^{-2t}\cos(3t)

L1{3(s+2)2+9}=e2tsin(3t)\mathcal{L}^{-1}\left\{\frac{3}{(s + 2)^2 + 9}\right\} = e^{-2t}\sin(3t)

f(t)=3e2tcos(3t)53e2tsin(3t)=e2t(3cos(3t)53sin(3t))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). \blacksquare

5.9 Worked Example: IVP with Discontinuous Forcing

Problem. Solve y+3y+2y=u2(t)y'' + 3y' + 2y = u_2(t), y(0)=0y(0) = 0, y(0)=0y'(0) = 0.

Solution

Solution. Take Laplace transforms:

s2Y+3sY+2Y=e2sss^2 Y + 3sY + 2Y = \frac{e^{-2s}}{s}

(s2+3s+2)Y=e2ss(s^2 + 3s + 2)Y = \frac{e^{-2s}}{s}

Y=e2ss(s+1)(s+2)Y = \frac{e^{-2s}}{s(s+1)(s+2)}

First find the partial fraction decomposition of 1s(s+1)(s+2)\frac{1}{s(s+1)(s+2)}:

1s(s+1)(s+2)=As+Bs+1+Cs+2\frac{1}{s(s+1)(s+2)} = \frac{A}{s} + \frac{B}{s+1} + \frac{C}{s+2}

1=A(s+1)(s+2)+Bs(s+2)+Cs(s+1)1 = A(s+1)(s+2) + Bs(s+2) + Cs(s+1)

s=0s = 0: 1=2A    A=1/21 = 2A \implies A = 1/2 s=1s = -1: 1=B    B=11 = -B \implies B = -1 s=2s = -2: 1=2C    C=1/21 = 2C \implies C = 1/2

So L1{1s(s+1)(s+2)}=12et+12e2t\mathcal{L}^{-1}\left\{\frac{1}{s(s+1)(s+2)}\right\} = \frac{1}{2} - e^{-t} + \frac{1}{2}e^{-2t}.

By the second shifting theorem:

y(t)=u2(t)[12e(t2)+12e2(t2)]y(t) = u_2(t)\left[\frac{1}{2} - e^{-(t-2)} + \frac{1}{2}e^{-2(t-2)}\right] for t0t \geq 0. \blacksquare

5.10 The Heaviside Function

The Heaviside (unit step) function is defined as

uc(t)={0t<c1tcu_c(t) = \begin{cases} 0 & t \lt c \\ 1 & t \geq c \end{cases}

It models a sudden switch being turned on at time t=ct = c. The second shifting theorem states:

L{uc(t)f(tc)}=ecsF(s)\mathcal{L}\{u_c(t)f(t - c)\} = e^{-cs}F(s)

Conversely, if Y(s)=ecsG(s)Y(s) = e^{-cs}G(s)Then y(t)=uc(t)g(tc)y(t) = u_c(t) \cdot g(t - c) where g=L1{G}g = \mathcal{L}^{-1}\{G\}.

5.11 The Dirac Delta Function

The Dirac delta function δ(tc)\delta(t - c) satisfies 0δ(tc)f(t)dt=f(c)\int_0^{\infty} \delta(t - c)f(t)\, dt = f(c) for continuous ff. It models an instantaneous Impulse.

L{δ(tc)}=ecs\mathcal{L}\{\delta(t - c)\} = e^{-cs}.

Properties:

  • δ(tc)=0\delta(t - c) = 0 for tct \neq c
  • δ(tc)dt=1\int_{-\infty}^{\infty} \delta(t - c)\, dt = 1
  • δ(tc)f(t)dt=f(c)\int_{-\infty}^{\infty} \delta(t - c)f(t)\, dt = f(c) (sifting property)

5.12 Worked Example: Impulse Response

Problem. An undamped spring-mass system (m=1m = 1, k=4k = 4) is struck by an impulse at t=0t = 0: y+4y=δ(t)y'' + 4y = \delta(t), y(0)=0y(0) = 0, y(0)=0y'(0) = 0. Find y(t)y(t).

Solution

Solution. Taking Laplace transforms:

s2Y+4Y=1s^2 Y + 4Y = 1

Y=1s2+4Y = \frac{1}{s^2 + 4}

y(t)=12sin(2t)y(t) = \frac{1}{2}\sin(2t).

This is the impulse response (Green’s function) for the operator D2+4D^2 + 4. \blacksquare

5.13 Worked Example: IVP with Laplace

Problem. Solve y6y+9y=6e3ty'' - 6y' + 9y = 6e^{3t}, y(0)=0y(0) = 0, y(0)=0y'(0) = 0.

Solution

Solution. Taking Laplace transforms:

s2Y6sY+9Y=6s3s^2 Y - 6sY + 9Y = \frac{6}{s - 3}

(s3)2Y=6s3(s - 3)^2 Y = \frac{6}{s - 3}

Y=6(s3)3Y = \frac{6}{(s - 3)^3}

Using the table: L{tneat}=n!(sa)n+1\mathcal{L}\{t^n e^{at}\} = \frac{n!}{(s-a)^{n+1}}.

y(t)=6t2e3t2!=3t2e3ty(t) = 6 \cdot \frac{t^2 e^{3t}}{2!} = 3t^2 e^{3t}. \blacksquare

5.14 Common Pitfalls for Laplace Transforms

:::caution Common Pitfall The Laplace transform of y(t)y'(t) is sY(s)y(0)sY(s) - y(0)Not sY(s)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(tc)u_c(t)f(t - c)Not uc(t)f(t)u_c(t)f(t). The function ff 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' = f(x, y), y(x0)=y0y(x_0) = y_0Define the sequence

ϕ0(x)=y0,ϕn+1(x)=y0+x0xf(t,ϕn(t))dt\phi_0(x) = y_0, \quad \phi_{n+1}(x) = y_0 + \int_{x_0}^x f(t, \phi_n(t))\, dt

If ff and f/y\partial f/\partial y are continuous, one shows by induction that (ϕn)(\phi_n) is uniformly Cauchy on some interval [x0h,x0+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.