Skip to content

The Laplace transform

Divide a signal by a test exponential, find which tests make the area exist, and read the Fourier transform off one line of the resulting plane.

Before this6.3 · 8.3 · 5 more
Chapter 9 · Lesson 1 of 3

First, the picture

Here a growing signal, e2tu(t)e^{2t}u(t), is divided by a test exponential that grows faster and faster, and the area under the ratio is shaded. Before you press play, decide: how fast must the test grow for the area to stop growing? Watch the readout “area”.

Test exponential

The signal e^(2t), divided by a test exponential e^(σt).

Divided by a flat test: the area never stops growing.

0.00 / 12.00 s
Describe this picture

Two plots: e2te^{2t} and the test exponential eσte^{\sigma t} on top, and their ratio e2t/eσte^{2t}/e^{\sigma t} below, with the area under the ratio shaded. For the first ten seconds σ\sigma rises steadily from 0 to 4, and the readout “area” follows: “grows without limit” while the test is slower than the signal, “grows by 1 every second, no limit” at σ=2\sigma=2, and then a number, “1.000” at σ=3\sigma=3 and “0.500” at σ=4\sigma=4. At 10 s the axis “decay rate σ” fades in under the plot, with a mark at +2+2 saying ”e2te^{2t} grows at +2” and the region σ>2\sigma>2 shaded and labelled “tests that grow faster: the area exists”; the clip holds that picture until 12 s. After the clip, σ can be dragged along that axis.

Divide by a faster exponential

Some signals grow, such as e2tu(t)e^{2t}u(t). The Fourier integral of From series to transform (8.1) cannot be added up for them, because the area under the winding never settles. The spikes of Transforms of common signals (8.3) rescue signals that hold steady, but not signals that grow.

The Laplace transform repairs this in one step. Divide the signal by a test exponential first, and keep track of which tests make the area exist. The border between the tests that work and the tests that fail turns out to be the signal’s own growth rate.

I need one meaning for the letter σ\sigma, and I will keep it for the whole page. It is the one it had in Complex exponentials & phasors (3.4) and First- and second-order systems (6.3): the rate of an exponential eσte^{\sigma t}. A positive σ\sigma grows, a negative σ\sigma dies, and the signal e2te^{2t} grows at rate +2+2.

In the picture at the top of the page I divided e2tu(t)e^{2t}u(t) by a test exponential eσte^{\sigma t}, with σ\sigma rising from 0 to 4. Divided by a flat test, σ=0\sigma=0, the area never stops growing. A test that grows slower than the signal leaves a ratio that still grows, and so does the area. At σ=2\sigma=2 the ratio is flat at 1, and the area grows by 1 every second, with no limit.

A test that grows faster makes the ratio die out, and then the area is a number: 1.000 at σ=3\sigma=3 and 0.500 at σ=4\sigma=4. The test must grow faster than the signal. The border is the signal’s own rate, +2, marked on the axis “decay rate σ” of 6.3 at the end of the clip. After the clip, drag σ along that axis and watch the readout.

The numbers agree with a one-line integral. The ratio is e(2−σ)te^{(2-\sigma)t}, and for σ>2\sigma>2 the exponent is negative, so

∫0∞e(2−σ)t dt=1σ−2.\int_0^\infty e^{(2-\sigma)t}\,dt=\frac{1}{\sigma-2}.

At σ=3\sigma=3 this is 11, and at σ=4\sigma=4 it is 0.50.5; at σ=2.5\sigma=2.5 it would be 22. For σ≤2\sigma\le2 the exponent is zero or positive and the integral has no limit.

The definition

Now I let the test spin as well. A test exponential este^{st} has the two numbers of 3.4: s=σ+jωs=\sigma+j\omega, with σ\sigma the rate at which it grows or dies and ω\omega the rate at which it spins. Dividing by este^{st} is multiplying by e−ste^{-st}. The Laplace transform adds up the result:

X(s)=∫−∞∞x(t) e−st dt,s=σ+jω.X(s)=\int_{-\infty}^{\infty}x(t)\,e^{-st}\,dt,\qquad s=\sigma+j\omega .

For each test exponential este^{st}, X(s)X(s) says how much of xx matches it. At a fixed σ\sigma the integral is 8.1’s Fourier transform of the signal x(t)/eσtx(t)/e^{\sigma t}, evaluated at ω\omega, because e−st=e−σte−jωte^{-st}=e^{-\sigma t}e^{-j\omega t}. So the new transform is the old one, applied after dividing by a test exponential.

The set of values of σ\sigma where the integral has a value is the region of convergence, written ROC. For our signal it is Re s>2\mathrm{Re}\,s>2, and there X(s)=1s−2X(s)=\dfrac{1}{s-2}, the number we computed above with σ=Re s\sigma=\mathrm{Re}\,s. Here Re s\mathrm{Re}\,s is the real part of ss.

Slicing the plane gives the Fourier transform

Take a signal that does not grow: x(t)=e−tu(t)x(t)=e^{-t}u(t), the decaying exponential of 8.1. Its transform is

X(s)=∫0∞e−te−st dt=1s+1,Re s>−1,X(s)=\int_0^\infty e^{-t}e^{-st}\,dt=\frac{1}{s+1},\qquad \mathrm{Re}\,s>-1,

because the exponent is −(s+1)t-(s+1)t and it dies out only when Re (s+1)>0\mathrm{Re}\,(s+1)>0. The formula blows up at s=−1s=-1. A value of ss where X(s)X(s) blows up is a pole of the signal’s X(s)X(s). It sits at the rate of one of the signal’s own exponentials, as the poles of First- and second-order systems (6.3) sat at the rates of a system’s free motion.

I draw XX on the same plane of rates as 6.3, decay rate σ\sigma across and spin rate ω\omega up. At every point I shade the length ∣X(s)∣\lvert X(s)\rvert, as in Complex numbers for signals (3.3): more colour means bigger. I shade only where the integral has a value, which is the ROC, and leave the rest grey. The edge of the ROC is a line at −1-1, and the colour is strongest toward the pole on it.

The next instrument adds a vertical line at one σ\sigma. The values of ∣X∣\lvert X\rvert along that line, plotted against ω\omega, are the profile. They are the Fourier transform of x(t)e−σtx(t)e^{-\sigma t}, since a fixed σ\sigma is the situation of the last section. Watch the profile as the line moves from σ=3\sigma=3 to σ=0\sigma=0.

Plane and slice

X(s) = 1/(s + 1) for e^(−t)u(t), shaded only where the integral has a value.

A slice far from the bright spot: a low profile.

0.00 / 9.50 s
Describe this picture

Two panels. The plane has the axes “decay rate σ” and “spin rate ω”, shaded by ∣X(s)∣\lvert X(s)\rvert with a colour bar labelled “|X(s)|: more colour is bigger” from 0 through 1 to “2 or more”; the part without a value is grey and labelled “no value here”, and the pole is marked “pole −1”. The profile panel, “|X| along the line”, has its spin-rate axis running the same way as the plane’s, with ticks at ω=±1,±2\omega=\pm1,\pm2 level with the plane’s; on narrow screens the panels stack, and the profile is drawn with ω\omega horizontal, joined to the plane by faint lines at the same ticks. The line moves from σ=3\sigma=3 to σ=0\sigma=0 in eight seconds. The readouts “peak” and “at ω = 1” follow it: 0.250 and 0.243 at σ=3\sigma=3, 0.400 and 0.371 at σ=1.5\sigma=1.5, 1.000 and 0.707 at σ=0\sigma=0, where the caption gives 0.7071. When the clip ends, the line can be dragged.

Far from the bright spot, at σ=3\sigma=3, the profile is low: its peak is 0.2500.250 at ω=0\omega=0, and it is 0.2430.243 at ω=1\omega=1. Closer to the pole the profile rises, to a peak of 0.4000.400 at σ=1.5\sigma=1.5, with 0.3710.371 at ω=1\omega=1. At the end the line sits on the spin-rate axis, σ=0\sigma=0. The peak is 1.0001.000, and the value at ω=1\omega=1 is 0.70710.7071.

On the spin-rate axis the slice is 8.1’s curve, with ∣X(j1)∣=0.7071\lvert X(j1)\rvert=0.7071 as in 8.1. So the Fourier transform is the slice σ=0\sigma=0 of the Laplace transform, which is X(jω)=X(s)∣s=jωX(j\omega)=X(s)\big\rvert_{s=j\omega}. This holds only when the line σ=0\sigma=0 is inside the ROC. Here it is, because 0>−10>-1. When the clip ends, drag the line yourself.

The numbers come from the formula. On a vertical line at σ\sigma,

∣X(σ+jω)∣=1(σ+1)2+ω2.\lvert X(\sigma+j\omega)\rvert=\frac{1}{\sqrt{(\sigma+1)^2+\omega^2}} .

At σ=3\sigma=3 the peak is 1/4=0.251/4=0.25, and at ω=1\omega=1 it is 1/17=0.24251/\sqrt{17}=0.2425. At σ=1.5\sigma=1.5 the peak is 1/2.5=0.41/2.5=0.4, and at ω=1\omega=1 it is 1/7.25=0.37141/\sqrt{7.25}=0.3714. At σ=0\sigma=0 the value at ω=1\omega=1 is 1/2=0.70711/\sqrt2=0.7071.

One formula, two signals: the region decides

Is the formula enough to name the signal? Take the left-sided signal −e−tu(−t)-e^{-t}u(-t), which is −e−t-e^{-t} for t<0t<0 and zero afterwards. It lives before time 0, and going back in time it grows without limit: at t=−3t=-3 it equals −e3=−20.09-e^{3}=-20.09. To tame it, the test exponential eσte^{\sigma t} must grow even faster going back, which needs σ<−1\sigma<-1.

Its transform is

X(s)=−∫−∞0e−(s+1)t dt=1s+1,Re s<−1.X(s)=-\int_{-\infty}^{0}e^{-(s+1)t}\,dt=\frac{1}{s+1},\qquad \mathrm{Re}\,s<-1 .

The integral has a value only if e−(s+1)te^{-(s+1)t} dies out as t→−∞t\to-\infty, which needs Re (s+1)<0\mathrm{Re}\,(s+1)<0. The formula is the same 1/(s+1)1/(s+1) as for e−tu(t)e^{-t}u(t). The region is on the other side of the pole.

The next instrument puts the two signals side by side. For each, a line sweeps σ\sigma from −3-3 to 33, and a light beside it says whether the integral adds up. Watch which side of −1-1 each region falls on.

Plane and slice: two signals

e^(−t)u(t) and −e^(−t)u(−t), the same formula 1/(s + 1).

e^(−t)u(t): it starts at 0 and dies out.

0.00 / 10.00 s
Describe this picture

Each signal has a time plot and a small copy of the plane; the picture has no control. A vertical line sweeps σ\sigma from −3-3 to 33 in three seconds, and a light beside it says “adds up” or “no value”. Where the integral adds up, the plane gets a plain flat fill, with no ∣X∣\lvert X\rvert map this time, and the rest stays grey. The first signal’s line sweeps from 0.5 s to 3.5 s, and its plane is then shaded from −1-1 onward. At 5.1 s the second signal fades in with its own plane, and its line sweeps from 5.4 s to 8.4 s; the first plane stays on screen. From 8.4 s both planes are labelled “1/(s+1)”, the second shaded up to −1-1, and the clip holds.

The first signal, e−tu(t)e^{-t}u(t), starts at 0 and dies out, and it adds up for decay rates to the right of −1-1. The second, −e−tu(−t)-e^{-t}u(-t), lives before 0 and grows going back, so the test must grow even faster going back: it adds up only to the left of −1-1. Same formula, opposite regions. The formula alone does not tell you the signal; the region does.

Only the first region contains the spin-rate axis, so only the first signal has a Fourier transform. In both pictures the region stops at the pole, because the integral cannot have a value where the formula blows up.

A signal that extends in both directions needs both conditions at once. For e−∣t∣e^{-\lvert t\rvert} the region is a strip, −1<Re s<1-1<\mathrm{Re}\,s<1; the worked example below computes it.

There is one rule for the general case. A right-sided signal is zero before some time, and a left-sided signal is zero after some time. For a right-sided signal that is a sum of exponentials, the test must outgrow every exponential in it, so its ROC is everything to the right of its rightmost pole. A left-sided signal has the mirror rule. The ROC never contains a pole.

I call σ1\sigma_1 the rate of the pole that forms the border: the rightmost pole for a right-sided signal, the leftmost for a left-sided one. In the first example σ1=2\sigma_1=2, and in the last two it is −1-1.

Where the clock starts

For a signal that starts at time zero, as in Kinds of signals (1.2), the whole integral can start at zero. But what if the signal has an impulse exactly at t=0t=0? Write the lower limit as 0−0^-, a time just before zero, and the impulse is counted. This is the unilateral transform:

X(s)=∫0−∞x(t) e−st dt.X(s)=\int_{0^-}^{\infty}x(t)\,e^{-st}\,dt .

An impulse is an arrow read by its area, as in 8.3, and the arrow at zero has area 1. A hammer blow at the instant the stopwatch starts must count, and the lower limit 0−0^- makes it count. The next lesson, Properties and the inverse Laplace transform (9.2), uses this to put starting values into the algebra.

The next instrument takes the signal δ(t)+e−tu(t)\delta(t)+e^{-t}u(t) at s=1s=1, so the integrand is x(t)e−tx(t)e^{-t}, and adds it up from a lower edge that slides from −1-1 to +1+1. Press play and watch the readout as the edge passes 0.

Where the clock starts

δ(t) + e^(−t)u(t), added up from a lower edge at s = 1.

Starting before 0: the impulse (1) plus the rest (0.5) is 1.5.

0.00 / 9.50 s
Describe this picture

The integrand of δ(t)+e−tu(t)\delta(t)+e^{-t}u(t) at s=1s=1, with the impulse as an arrow labelled “1”, shaded from a lower edge onward. The edge slides from −1-1 to +1+1 in eight seconds, and the readout “integral from the edge” reads 1.500 at −1-1. With the edge at 0 it gives way to two readouts, “edge just before 0 (0⁻): 1.5” and “just after (0⁺): 0.5”, which stay in quieter type after the edge passes 0. At +1+1 it reads 0.068. The clip then holds, and the edge can be dragged.

Starting before 0, the impulse (1) plus the rest (0.5) is 1.5. As the edge passes 0 the integral drops from 1.5 to 0.5: the drop is the impulse’s area, 1. Start later, at +1+1, and you lose what came before: 0.068 is left. When the clip holds, drag the edge yourself. The three numbers agree with the integrals. The impulse gives 11, and the exponential part gives ∫0∞e−2t dt=0.5\int_0^\infty e^{-2t}\,dt=0.5. From t=1t=1 the exponential part gives ∫1∞e−2t dt=12e−2=0.0677\int_1^\infty e^{-2t}\,dt=\tfrac12e^{-2}=0.0677.

The maths behind it · eigenvectors of LTI systems

A growing or dying spiral este^{st} goes through any linear time-invariant system and comes out as the same spiral times one number, just as ejωte^{j\omega t} did in Fourier series and LTI systems (7.4). A vector that a matrix only stretches is an eigenvector. These spirals are the eigenvectors of every such system, and X(s)X(s) measures how much of each a signal holds. The ROC is where that measurement has a value.

The maths behind it · moment-generating functions

Take a waiting time with density p(t)p(t) on t≥0t\ge0, a curve of area 1 that says how likely each wait is. Then ∫0∞p(t)e−stdt\int_0^\infty p(t)e^{-st}dt is the average of e−s⋅(waiting time)e^{-s\cdot(\text{waiting time})} over many waits. Statisticians use the same average with the sign of ss flipped and call it the moment-generating function. For Re s≥0\mathrm{Re}\,s\ge0 it always has a value; how far it reaches toward negative Re s\mathrm{Re}\,s depends on how fast the density’s tail dies, which is its ROC.

Worked example

  1. Growing exponential. e2tu(t)→1s−2e^{2t}u(t)\to\dfrac{1}{s-2}, Re s>2\mathrm{Re}\,s>2. At s=2.5s=2.5 it is 22, at s=3s=3 it is 11, and at s=4s=4 it is 0.50.5.
  2. Step. u(t)→∫0∞e−st dt=1su(t)\to\displaystyle\int_0^\infty e^{-st}\,dt=\dfrac1s, Re s>0\mathrm{Re}\,s>0. At s=2s=2 it is 0.50.5.
  3. Decaying exponential. e−tu(t)→1s+1e^{-t}u(t)\to\dfrac{1}{s+1}, Re s>−1\mathrm{Re}\,s>-1. On s=jωs=j\omega it is 1/(1+jω)1/(1+j\omega), so ∣X(j1)∣=0.7071\lvert X(j1)\rvert=0.7071. The slices are those of the second instrument: 0.250.25 and 0.24250.2425 at σ=3\sigma=3, 0.40.4 and 0.37140.3714 at σ=1.5\sigma=1.5.
  4. Two-sided signal. For e−∣t∣e^{-\lvert t\rvert} split the integral at 0. The part t>0t>0 is e−tu(t)e^{-t}u(t), giving 11+s\dfrac{1}{1+s} for Re s>−1\mathrm{Re}\,s>-1. The part t<0t<0 is ete^{t} for t<0t<0, giving 11−s\dfrac{1}{1-s} for Re s<1\mathrm{Re}\,s<1. Both must hold at once, so X(s)=11−s+11+s=21−s2,−1<Re s<1.X(s)=\frac{1}{1-s}+\frac{1}{1+s}=\frac{2}{1-s^2},\qquad -1<\mathrm{Re}\,s<1 . At s=0s=0 this is 22, the area under e−∣t∣e^{-\lvert t\rvert}. At s=0.5s=0.5 it is 2.6672.667.
  5. Left-sided signal. −e−tu(−t)→1s+1-e^{-t}u(-t)\to\dfrac{1}{s+1}, Re s<−1\mathrm{Re}\,s<-1. At s=−2s=-2 it is −1-1, and at t=−3t=-3 the signal is −e3=−20.09-e^{3}=-20.09.
  6. Unilateral and impulses. δ(t)+e−tu(t)→1+1s+1\delta(t)+e^{-t}u(t)\to1+\dfrac{1}{s+1}. At s=1s=1 this is 1.51.5. From 0+0^+ it is 0.50.5, and from t=1t=1 it is 12e−2=0.0677\tfrac12e^{-2}=0.0677.
  7. Cosine and sine from two arrows. Write cos⁡ω1t\cos\omega_1t as half the sum of two arrows, 12(ejω1t+e−jω1t)\tfrac12\left(e^{j\omega_1t}+e^{-j\omega_1t}\right). Each has a transform of the form e−atu(t)→1s+ae^{-at}u(t)\to\dfrac{1}{s+a}, with a=∓jω1a=\mp j\omega_1, valid for Re s>0\mathrm{Re}\,s>0. Here ω1\omega_1 is the frequency of one tone, as in 8.3. Then cos⁡ω1t u(t)→12(1s−jω1+1s+jω1)=ss2+ω12.\cos\omega_1t\,u(t)\to\tfrac12\left(\frac{1}{s-j\omega_1}+\frac{1}{s+j\omega_1}\right)=\frac{s}{s^2+\omega_1^2}. For the sine, subtract and divide by 2j2j: sin⁡ω1t u(t)→12j(1s−jω1−1s+jω1)=ω1s2+ω12.\sin\omega_1t\,u(t)\to\frac{1}{2j}\left(\frac{1}{s-j\omega_1}-\frac{1}{s+j\omega_1}\right)=\frac{\omega_1}{s^2+\omega_1^2}. At ω1=1\omega_1=1 and s=1s=1 both are 0.50.5.

Where you’ll meet this

The unilateral transform is the one engineers use to solve circuits that start from a charged state, and 9.2 does exactly that with the capacitor of Differential equations and analog systems (6.2). The poles you saw as marks on a plane are the system’s own rates of motion, and Poles, zeros and the s-plane (9.3) reads a system’s behaviour off their positions.

The same idea, with sums in place of integrals, becomes the z-transform for sampled signals in a later chapter. There the ROC will matter for exactly the reason it did here: the formula alone does not name the signal.

Reference card

QuantityFormulaNotes
Bilateral Laplace transformX(s)=∫−∞∞x(t)e−stdtX(s)=\displaystyle\int_{-\infty}^{\infty}x(t)e^{-st}dt, s=σ+jωs=\sigma+j\omegathe Fourier transform of x(t)e−σtx(t)e^{-\sigma t}; needs an ROC; σ\sigma is always the rate of eσte^{\sigma t}
UnilateralX(s)=∫0−∞x(t)e−stdtX(s)=\displaystyle\int_{0^-}^{\infty}x(t)e^{-st}dtcounts an impulse at 0; starting values (9.2)
Relation to FourierX(jω)=X(s)∣s=jωX(j\omega)=X(s)\big\rvert_{s=j\omega}only if the ROC contains the spin-rate axis
ROC shapesright-sided: Re s>σ1\mathrm{Re}\,s>\sigma_1; left-sided: Re s<σ1\mathrm{Re}\,s<\sigma_1; two-sided: a stripσ1\sigma_1: the rate of the border pole (rightmost for right-sided, leftmost for left-sided); never contains a pole
δ(t)\delta(t)11all ss
u(t)u(t)1/s1/sRe s>0\mathrm{Re}\,s>0
e−atu(t)e^{-at}u(t)1s+a\dfrac{1}{s+a}Re s>−a\mathrm{Re}\,s>-a
−e−atu(−t)-e^{-at}u(-t)1s+a\dfrac{1}{s+a}Re s<−a\mathrm{Re}\,s<-a
cos⁡ω1t u(t)\cos\omega_1t\,u(t)ss2+ω12\dfrac{s}{s^2+\omega_1^2}Re s>0\mathrm{Re}\,s>0; worked example 7
sin⁡ω1t u(t)\sin\omega_1t\,u(t)ω1s2+ω12\dfrac{\omega_1}{s^2+\omega_1^2}Re s>0\mathrm{Re}\,s>0; worked example 7

End of lesson 9.1

Where to go next.

Phasorium
LibraryEvery lesson, in order

Parts

About Phasorium
Look