Skip to content

Laplace Transforms | Mathematics

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.

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.

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

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}

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).

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

Section titled “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 \infty So 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

Section titled “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

Section titled “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

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\}.

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)

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

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

Section titled “5.14 Common Pitfalls for Laplace Transforms”

5.15 Intuition: Why Are Laplace Transforms Useful?

Section titled “5.15 Intuition: Why Are Laplace Transforms Useful?”

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), solve for Y(s)Y(s), and then invert.

The key insight: The derivative property L{f(t)}=sF(s)f(0)\mathcal{L}\{f'(t)\} = sF(s) - f(0) converts differentiation into multiplication by ss. This means a differential equation like y+3y+2y=ety'' + 3y' + 2y = e^t becomes (s2+3s+2)Y(s)=1s1+initial condition terms(s^2 + 3s + 2)Y(s) = \frac{1}{s-1} + \text{initial condition terms}, which is just algebra.

Why the convolution theorem matters: The convolution (fg)(t)=0tf(τ)g(tτ)dτ(f * g)(t) = \int_0^t f(\tau)g(t - \tau)\,d\tau is hard to compute directly. But L{fg}=F(s)G(s)\mathcal{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 ss can be interpreted as a complex frequency. F(s)F(s) decomposes the signal f(t)f(t) into exponential components este^{st}. The poles of F(s)F(s) (values of ss where F(s)F(s) blows up) determine the natural modes of the system, while the zeros determine the frequencies that are suppressed.

5.16 Worked Example: System of ODEs via Laplace

Section titled “5.16 Worked Example: System of ODEs via Laplace”

Problem. Solve the system: x=3x2yy=2x2y\begin{aligned} x' &= 3x - 2y \\ y' &= 2x - 2y \end{aligned} with x(0)=1x(0) = 1, y(0)=0y(0) = 0.

Solution

Take Laplace transforms: sX1=3X2YsY=2X2Y\begin{aligned} sX - 1 &= 3X - 2Y \\ sY &= 2X - 2Y \end{aligned}

Rearranging: (s3)X+2Y=12X+(s+2)Y=0\begin{aligned} (s - 3)X + 2Y &= 1 \\ -2X + (s + 2)Y &= 0 \end{aligned}

From the second equation: X=(s+2)Y2X = \frac{(s+2)Y}{2}.

Substituting into the first: (s3)(s+2)Y2+2Y=1(s - 3) \cdot \frac{(s+2)Y}{2} + 2Y = 1

(s3)(s+2)Y+4Y2=1\frac{(s-3)(s+2)Y + 4Y}{2} = 1

(s2s6+4)Y=2(s^2 - s - 6 + 4)Y = 2

(s2s2)Y=2(s^2 - s - 2)Y = 2

Y=2(s2)(s+1)=2/3s22/3s+1Y = \frac{2}{(s-2)(s+1)} = \frac{2/3}{s-2} - \frac{2/3}{s+1}

y(t)=23e2t23ety(t) = \frac{2}{3}e^{2t} - \frac{2}{3}e^{-t}

For XX: X=(s+2)22(s2)(s+1)=s+2(s2)(s+1)=4/3s21/3s+1X = \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(t)=43e2t13etx(t) = \frac{4}{3}e^{2t} - \frac{1}{3}e^{-t}

Verification: x(0)=4/31/3=1x(0) = 4/3 - 1/3 = 1 \checkmark, y(0)=2/32/3=0y(0) = 2/3 - 2/3 = 0 \checkmark.

x=83e2t+13etx' = \frac{8}{3}e^{2t} + \frac{1}{3}e^{-t}, 3x2y=4e2tet43e2t+43et=83e2t+13et3x - 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} \checkmark. \blacksquare

5.17 Intuition: Why Are Laplace Transforms Useful?

Section titled “5.17 Intuition: Why Are Laplace Transforms Useful?”

The Laplace transform converts differential equations into algebraic equations. By transforming a function f(t)f(t) into F(s)=0estf(t)dtF(s) = \int_0^\infty e^{-st}f(t)\,dt, differentiation in the time domain becomes multiplication by ss in the frequency domain. This means a linear ODE with constant coefficients becomes a polynomial equation in ss, 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(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.

flowchart TD
A[5_Laplace Transforms] --> B[Key Concepts]
A --> C[Core Principles]
A --> D[Practical Applications]
B --> E[Fundamental definitions]
C --> F[Design patterns]
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 ss 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