The Laplace Transform

Transforming differential equations into algebra

Everything we have done so far in this course has involved solving differential equations directly in the time domain. We manipulated derivatives, found integrating factors, guessed solutions, and carefully tracked initial conditions at the end.

The Laplace transform offers a radically different approach. Instead of working with functions of time and their derivatives, we transform the entire problem into a new domain where differential equations become algebraic equations. Derivatives become multiplication. Initial conditions appear automatically. The calculus disappears.

This is not a trick. It is a change of perspective, and perspectives matter. Engineers do not just use the Laplace transform because it simplifies calculations; they think in terms of it. The "s-domain" is where transfer functions live, where system stability becomes visible, and where the frequency response of circuits and mechanical systems reveals itself at a glance.

The Definition

The Laplace transform of a function 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

This is an improper integral, meaning the upper limit is infinity. To make sense of it:

0estf(t)dt=limT0Testf(t)dt\int_0^\infty e^{-st} f(t) \, dt = \lim_{T \to \infty} \int_0^T e^{-st} f(t) \, dt

The integral converges if this limit exists and is finite. It diverges if the limit is infinite or does not exist. The exponential este^{-st} acts as a weighting factor that decays as tt increases, and this decay helps the integral converge for many functions.

Interactive: The Defining Integral

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

For f(t)=etf(t) = e^{-t}, the integral converges to:

F(1.0)=11.0+1=0.500F(1.0) = \frac{1}{1.0 + 1} = 0.500

The shaded area represents the integral. As s increases, the integrand decays faster, making the integral smaller.

The integral converges when the exponential decay este^{-st} outpaces the growth of f(t)f(t). For a function that grows no faster than ecte^{ct} for some constant cc, the transform exists for all s>cs > c. This threshold cc is called the abscissa of convergence.

For example:

  • f(t)=1f(t) = 1 grows slower than any exponential, so the transform exists for s>0s > 0
  • f(t)=e3tf(t) = e^{3t} requires s>3s > 3 for convergence
  • f(t)=et2f(t) = e^{t^2} grows faster than any exponential, so the transform does not exist

Think of the transform as extracting information about the original function. The weighting este^{-st} suppresses the function's behavior at large tt while preserving information about its early behavior and overall character. Larger values of ss weight early-time behavior more heavily, since este^{-st} decays faster.

Computing Basic Transforms

Rather than computing every transform from the integral, we build a table of common transforms that we can use repeatedly. Let us derive a few to see how the process works.

Transform of 1:

L{1}=0est1dt=[1sest]0=0(1s)=1s\mathcal{L}\{1\} = \int_0^\infty e^{-st} \cdot 1 \, dt = \left[ -\frac{1}{s} e^{-st} \right]_0^\infty = 0 - \left( -\frac{1}{s} \right) = \frac{1}{s}

This holds for s>0s > 0, where the exponential decays to zero as tt \to \infty.

Transform of eate^{at}:

L{eat}=0esteatdt=0e(sa)tdt=1sa\mathcal{L}\{e^{at}\} = \int_0^\infty e^{-st} e^{at} \, dt = \int_0^\infty e^{-(s-a)t} \, dt = \frac{1}{s-a}

This holds for s>as > a. Notice the elegant pattern: the exponential function eate^{at} has a transform with a pole at s=as = a. The parameter in the time domain becomes a shift in the s-domain.

Transform of sin(bt)\sin(bt):

This requires integration by parts twice, but the result is:

L{sin(bt)}=bs2+b2\mathcal{L}\{\sin(bt)\} = \frac{b}{s^2 + b^2}

The oscillation frequency bb appears in both numerator and denominator. The denominator s2+b2s^2 + b^2 has roots at s=±bis = \pm bi, pure imaginary numbers reflecting the purely oscillatory nature of sine.

Interactive: See Transforms Side by Side

Time Domain

f(t)=1f(t) = 1

s-Domain (Laplace)

F(s)=1sF(s) = \frac{1}{s}

The blue curve in time domain transforms to the red curve in s-domain.

Select different functions to see their time-domain behavior and corresponding s-domain transforms. Notice how the character of each function, whether it grows, decays, or oscillates, is reflected in the shape of its transform.

The Transform Table

Working with Laplace transforms, you will refer to a table constantly. Here are the essential transforms:

Common Laplace Transform Pairs

f(t)f(t)F(s)=L{f(t)}F(s) = \mathcal{L}\{f(t)\}Region
111s\frac{1}{s}s > 0
tt1s2\frac{1}{s^2}s > 0
tnt^nn!sn+1\frac{n!}{s^{n+1}}s > 0
eate^{at}1sa\frac{1}{s-a}s > a
sin(bt)\sin(bt)bs2+b2\frac{b}{s^2+b^2}s > 0
cos(bt)\cos(bt)ss2+b2\frac{s}{s^2+b^2}s > 0
eatsin(bt)e^{at}\sin(bt)b(sa)2+b2\frac{b}{(s-a)^2+b^2}s > a
eatcos(bt)e^{at}\cos(bt)sa(sa)2+b2\frac{s-a}{(s-a)^2+b^2}s > a

Click any row to learn more about that transform.

These transforms form the building blocks. Memorizing them, or at least becoming deeply familiar with them, is essential. The integral definition is important for understanding, but practical calculations rely on the table.

Linearity

The Laplace transform is linear. If ff and gg have transforms FF and GG, and aa and bb are constants, then:

L{af(t)+bg(t)}=aF(s)+bG(s)\mathcal{L}\{af(t) + bg(t)\} = aF(s) + bG(s)

This property is crucial. It allows us to transform complicated functions by breaking them into simpler pieces. The transform of a sum is the sum of the transforms.

Interactive: Linearity in Action

L{af(t)+bg(t)}=aL{f}+bL{g}\mathcal{L}\{af(t) + bg(t)\} = a\mathcal{L}\{f\} + b\mathcal{L}\{g\}

Time Domain

f(t)=t, g(t)=e^(-t), af+bg

s-Domain

2.01s2+1.01s+12.0 \cdot \frac{1}{s^2} + 1.0 \cdot \frac{1}{s+1}

Linearity is essential: you can transform complicated functions by breaking them into simpler pieces, transforming each, and combining the results.

Adjust the coefficients to see how linear combinations in the time domain correspond to linear combinations in the s-domain. The structure is preserved.

The Key Property: Transforming Derivatives

Here is where the magic happens. Consider the transform of a derivative:

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

Using integration by parts with u=estu = e^{-st} and dv=f(t)dtdv = f'(t) \, dt:

L{f(t)}=[estf(t)]0+s0estf(t)dt\mathcal{L}\{f'(t)\} = \left[ e^{-st} f(t) \right]_0^\infty + s \int_0^\infty e^{-st} f(t) \, dt

The first term evaluates to f(0)-f(0) (assuming the limit at infinity vanishes). The second term is sF(s)sF(s). Therefore:

L{f(t)}=sF(s)f(0)\mathcal{L}\{f'(t)\} = sF(s) - f(0)

Differentiation in time becomes multiplication by ss in the transform domain, with a correction for the initial value. This is the fundamental property that makes the Laplace transform so powerful for differential equations.

For the second derivative, apply the rule twice:

L{f(t)}=sL{f(t)}f(0)=s(sF(s)f(0))f(0)=s2F(s)sf(0)f(0)\mathcal{L}\{f''(t)\} = s\mathcal{L}\{f'(t)\} - f'(0) = s(sF(s) - f(0)) - f'(0) = s^2 F(s) - sf(0) - f'(0)

Interactive: Differentiation Becomes Multiplication

Time Domain

f(t)

s-Domain

F(s), sF(s)
L{f(t)}=sF(s)f(0)\mathcal{L}\{f'(t)\} = sF(s) - f(0)

Differentiation in time becomes multiplication by s in the transform domain, minus the initial value. This is why the Laplace transform converts differential equations into algebraic equations.

Watch how differentiation in the time domain corresponds to multiplication by ss and subtraction of initial values in the s-domain. The derivative property is what transforms differential equations into algebraic equations.

From Differential Equation to Algebra

Consider the initial value problem:

y+3y+2y=0,y(0)=1,y(0)=0y'' + 3y' + 2y = 0, \quad y(0) = 1, \quad y'(0) = 0

Let Y(s)=L{y(t)}Y(s) = \mathcal{L}\{y(t)\}. Transform each term:

  • L{y}=s2Ysy(0)y(0)=s2Ys\mathcal{L}\{y''\} = s^2 Y - sy(0) - y'(0) = s^2 Y - s
  • L{y}=sYy(0)=sY1\mathcal{L}\{y'\} = sY - y(0) = sY - 1
  • L{y}=Y\mathcal{L}\{y\} = Y

The transformed equation is:

(s2Ys)+3(sY1)+2Y=0(s^2 Y - s) + 3(sY - 1) + 2Y = 0

Collect terms:

(s2+3s+2)Y=s+3(s^2 + 3s + 2)Y = s + 3

Solve for YY:

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

No calculus involved. The differential equation has become an algebraic equation, and we solved it by factoring and rearranging. The initial conditions appeared naturally, not as an afterthought.

Step-by-Step: DE to Algebraic Equation

Step 1: Original Differential Equation

Time domain:

y+3y+2y=0y'' + 3y' + 2y = 0

We start with a second-order linear ODE with initial conditions y(0) = 1, y'(0) = 0.

The Laplace transform converts a differential equation (calculus) into an algebraic equation (algebra). Initial conditions are built into the process, not added afterward.

The remaining step is to find y(t)y(t) from Y(s)Y(s), the inverse transform, which is the subject of the next chapter.

The Variable ss: Complex Frequency

What is ss, really? While we often work with ss as if it were a positive real number, ss is properly a complex variable: s=σ+iωs = \sigma + i\omega.

The real part σ\sigma relates to exponential growth or decay. The imaginary part ω\omega relates to oscillation frequency. In this light, the Laplace transform decomposes a signal into components of the form est=eσteiωte^{st} = e^{\sigma t} e^{i\omega t}, each with a specific growth rate and frequency.

For our purposes, you can think of ss as a parameter that must be "large enough" for the transform integral to converge. When solving differential equations, we manipulate Y(s)Y(s) algebraically and only worry about the precise domain of convergence when it matters for the inverse transform.

Engineers call the s-domain the frequency domain, though this is slightly imprecise since ss encodes both frequency (through its imaginary part) and damping (through its real part). The Laplace transform generalizes the Fourier transform to handle transient, non-periodic signals and systems with initial conditions.

Why This Approach Works

The Laplace transform works because it converts operations that are hard (differentiation, solving ODEs) into operations that are easy (multiplication, algebra).

The price we pay is twofold: first, we need to learn the transform pairs and properties; second, we need to convert back from Y(s)Y(s) to y(t)y(t). But this price is well worth paying. Once you are fluent with the method, you can solve linear constant-coefficient ODEs, even complicated ones, almost mechanically.

Moreover, the s-domain representation reveals things that the time domain hides. The poles of Y(s)Y(s), where the denominator is zero, tell you about the natural frequencies of the system. The behavior of Y(s)Y(s) for large ss tells you about high-frequency response. Transfer functions, the ratio of output to input in the s-domain, become the central objects of study in control theory and signal processing.

What Comes Next

We have seen how to transform a differential equation into an algebraic equation for Y(s)Y(s). But how do we get back to y(t)y(t)? That is the topic of the next chapter: the inverse Laplace transform. There, we will learn partial fraction decomposition, use our transform table in reverse, and complete the solution process.

The inverse transform is where the algebra we have created here becomes a solution we can use.

Key Takeaways

  • The Laplace transform converts f(t)f(t) into F(s)F(s) via L{f}=0estf(t)dt\mathcal{L}\{f\} = \int_0^\infty e^{-st} f(t) \, dt
  • Key transforms: L{1}=1/s\mathcal{L}\{1\} = 1/s, L{tn}=n!/sn+1\mathcal{L}\{t^n\} = n!/s^{n+1}, L{eat}=1/(sa)\mathcal{L}\{e^{at}\} = 1/(s-a), L{sinbt}=b/(s2+b2)\mathcal{L}\{\sin bt\} = b/(s^2+b^2), L{cosbt}=s/(s2+b2)\mathcal{L}\{\cos bt\} = s/(s^2+b^2)
  • The derivative property L{f}=sF(s)f(0)\mathcal{L}\{f'\} = sF(s) - f(0) is what makes the transform useful for DEs
  • Differential equations become algebraic equations in Y(s)Y(s), with initial conditions built in
  • The transform is linear: L{af+bg}=aF+bG\mathcal{L}\{af + bg\} = aF + bG
  • The s-domain provides insight into system behavior: poles, frequencies, stability