This is the hands-on companion to the convolution theorem note . There we proved L { f ∗ g } = F ( s ) G ( s ) . Here we drill it until every kind of problem — every sign, every degenerate input, every trap — has been met head-on.
Intuition What "convolution" means, one more time, in one picture
Convolution ( f ∗ g ) ( t ) = ∫ 0 t f ( τ ) g ( t − τ ) d τ is a blend : for a fixed present moment t , we slide a flipped copy of g across f and add up the overlap. The variable τ is a clock that runs from 0 (the beginning) to t (now). The two arguments τ and t − τ always add up to t — that is the whole secret.
Figure s01 — reading the picture. The horizontal axis is the running clock τ ; the vertical axis is the leftover time t − τ . Together they label the quarter-plane of all ( τ , t − τ ) pairs. Each faint black line is a diagonal for one fixed present moment; the bold red diagonal is the one for t = 4 . Every point on that red line obeys τ + ( t − τ ) = t — the two black arrows show the horizontal piece τ and the vertical piece t − τ adding up to t . Convolution takes the product f ( τ ) g ( t − τ ) at every point of that red diagonal and adds it all up. That summing-along-a-diagonal is convolution.
Before any symbol scares you, here is the tiny dictionary we will reuse:
Definition The three symbols on this page
L { f } = F ( s ) : the Laplace transform — a machine that eats a time-function f ( t ) and spits out an s -function F ( s ) . (Built in Laplace Transform — definition and properties .)
L − 1 : the inverse — the same machine run backwards, turning F ( s ) back into f ( t ) .
∗ : convolution , the blend integral above. Never confuse it with ordinary multiplication × .
Definition The standing hypothesis (causal / one-sided functions)
Everything on this page assumes f ( t ) and g ( t ) are causal : they are 0 for t < 0 and defined for t ≥ 0 . This is the standard Laplace-transform setting, because L { f } = ∫ 0 ∞ e − s t f ( t ) d t only looks at t ≥ 0 .
Why this matters for the limits. Because g ( u ) = 0 whenever its argument u < 0 , the factor g ( t − τ ) is 0 as soon as τ > t . So the "true" integral ∫ 0 ∞ f ( τ ) g ( t − τ ) d τ collapses to ∫ 0 t — the upper limit is t , not ∞ . This causality assumption is exactly what licenses the [ 0 , t ] limits and the convolution theorem L { f ∗ g } = F ( s ) G ( s ) used throughout.
Every problem this topic can throw at you falls into one of these cells. The worked examples below are each tagged with the cell they cover, and together they hit all of them.
#
Cell (scenario class)
What is tricky about it
Example
C1
Degenerate input : convolve with the constant 1
Is f ∗ 1 just f ? (No — it integrates!)
Ex 1
C2
Both factors polynomial (1/ s m ⋅ 1/ s n )
Pure power algebra, no trig
Ex 2
C3
Repeated factor 1/ ( s 2 + a 2 ) 2
Same trig blended with itself
Ex 3
C4
Mixed exponential × trig
Sign bookkeeping in the integral
Ex 4
C5
Two different exponentials (e a t , e b t )
The degenerate limit a → b
Ex 5
C6
Solve an ODE via convolution
Convert s -answer to a blend
Ex 6
C7
Volterra integral equation
Recognise the hidden ∗
Ex 7
C8
Word problem / real system (impulse response)
Translate physics → convolution
Ex 8
C9
Exam twist : commutativity + which factor to flip
Pick the easier integrand
Ex 9
Where a cell has a limiting / degenerate sub-case (like a → b ), the example handles it explicitly so you never meet an unshown scenario.
t 2 ∗ 1 and check the "convolve-with-1 = integrate" rule.
Forecast first: guess — does t 2 ∗ 1 equal t 2 (like multiplying by 1) or something bigger? Write down your guess before reading on.
Step 1. Identify the pieces: f ( τ ) = τ 2 and g ( t − τ ) = 1 (the constant function is 1 everywhere).
Why this step? Convolution needs both factors as functions of the running variable τ . The constant function stays 1 no matter the argument.
Step 2. Write the integral:
t 2 ∗ 1 = ∫ 0 t τ 2 ⋅ 1 d τ .
Why this step? Direct substitution into ( f ∗ g ) ( t ) = ∫ 0 t f ( τ ) g ( t − τ ) d τ . Because g ≡ 1 , convolving with 1 is literally integrating from 0 to t — that is the intuition to keep.
Step 3. Integrate:
∫ 0 t τ 2 d τ = 3 t 3 .
Verify (via the theorem): L { t 2 } = s 3 2 and L { 1 } = s 1 . Their product is s 4 2 , and L − 1 { s 4 2 } = 3 ! 2 t 3 = 3 t 3 . ✓
Moral: convolving with 1 is not the identity — it is integration. Your forecast of "just t 2 " would have been the classic trap.
L − 1 { s 3 s 2 1 } = L − 1 { s 5 1 } by convolution.
Forecast: two power factors → the answer should be a single power of t . Which power? Guess.
Step 1. Split: F ( s ) = s 3 1 ⇒ f ( t ) = 2 t 2 , and G ( s ) = s 2 1 ⇒ g ( t ) = t .
Why this step? L { t n } = s n + 1 n ! , so s 3 1 ↔ 2 ! t 2 and s 2 1 ↔ 1 ! t 1 .
Step 2. Convolve:
( f ∗ g ) ( t ) = ∫ 0 t 2 τ 2 ( t − τ ) d τ .
Why this step? Apply the theorem; keep f ( τ ) = τ 2 /2 inside and flip g to g ( t − τ ) = t − τ .
Step 3. Expand and integrate:
2 1 ∫ 0 t ( t τ 2 − τ 3 ) d τ = 2 1 [ 3 t t 3 − 4 t 4 ] = 2 1 ⋅ 12 t 4 = 24 t 4 .
Verify: L { t 4 /24 } = 24 s 5 4 ! = 24 s 5 24 = s 5 1 . ✓ The power is 4 , matching 1/ s 5 .
y ′′ + 4 y = 0 won't need it, but this does: L − 1 { ( s 2 + 4 ) 2 1 } .
Forecast: a squared denominator of sin -type. Expect sin and a t cos term (resonance signature). Guess the exact form.
Step 1. Recognise ( s 2 + 4 ) 2 1 = s 2 + 4 1 ⋅ s 2 + 4 1 , and L − 1 { s 2 + 4 1 } = 2 sin 2 t .
Why this step? L { sin a t } = s 2 + a 2 a , so with a = 2 we get s 2 + 4 2 ; dividing by 2 gives our factor.
Step 2. Convolve g with itself, where g ( t ) = 2 1 sin 2 t :
y = ∫ 0 t 2 1 sin ( 2 τ ) ⋅ 2 1 sin ( 2 ( t − τ ) ) d τ = 4 1 ∫ 0 t sin 2 τ sin ( 2 t − 2 τ ) d τ .
Why this step? Direct application; the two arguments 2 τ and 2 ( t − τ ) add to 2 t .
Step 3. Product-to-sum: sin A sin B = 2 1 [ cos ( A − B ) − cos ( A + B )] with A = 2 τ , B = 2 t − 2 τ , so A − B = 4 τ − 2 t and A + B = 2 t :
y = 8 1 ∫ 0 t [ cos ( 4 τ − 2 t ) − cos 2 t ] d τ .
Why this step? The cos 2 t term is constant in τ ; the other integrates cleanly.
Step 4. Integrate:
8 1 [ 4 s i n ( 4 τ − 2 t ) − τ cos 2 t ] 0 t = 8 1 [ 4 s i n 2 t − s i n ( − 2 t ) − t cos 2 t ] = 8 1 [ 4 2 s i n 2 t − t cos 2 t ] = 8 1 [ 2 s i n 2 t − t cos 2 t ] .
Why this step? At the upper limit τ = t : sin ( 4 t − 2 t ) = sin 2 t ; at the lower limit τ = 0 : sin ( − 2 t ) = − sin 2 t , so the bracket difference gives sin 2 t − ( − sin 2 t ) = 2 sin 2 t , and 2 sin 2 t /4 = sin 2 t /2 .
Step 5. Distribute the 8 1 carefully:
y = 8 1 ⋅ 2 s i n 2 t − 8 1 ⋅ t cos 2 t = 16 s i n 2 t − 8 t c o s 2 t .
Why this step? 8 1 ⋅ 2 1 = 16 1 for the first term, but 8 1 ⋅ 1 = 8 1 for the second — the two coefficients are different . Factoring out 16 1 then gives:
y = 16 sin 2 t − 8 t cos 2 t = 16 1 ( sin 2 t − 2 t cos 2 t )
Verify: the known formula L − 1 { ( s 2 + a 2 ) 2 1 } = 2 a 3 1 ( sin a t − a t cos a t ) . With a = 2 : 16 1 ( sin 2 t − 2 t cos 2 t ) . ✓
L − 1 { ( s − 1 ) ( s 2 + 1 ) 1 } .
Forecast: an e t from the first factor blended with sin t . Guess whether the answer is bounded or grows.
Step 1. Split: F = s − 1 1 ⇒ f ( t ) = e t ; G = s 2 + 1 1 ⇒ g ( t ) = sin t .
Why this step? L { e a t } = s − a 1 (here a = 1 ), L { sin t } = s 2 + 1 1 .
Step 2. Convolve, flipping g :
∫ 0 t e τ sin ( t − τ ) d τ .
Why this step? Keep the exponential as f ( τ ) = e τ (easy to integrate) inside; theorem says one argument must be t − τ .
Step 3. Expand sin ( t − τ ) = sin t cos τ − cos t sin τ and pull the constants (in τ ) out:
sin t ∫ 0 t e τ cos τ d τ − cos t ∫ 0 t e τ sin τ d τ .
Why this step? sin t , cos t do not depend on the integration variable τ , so they are constants for the integral.
Step 4. Use standard integrals ∫ e τ cos τ d τ = 2 1 e τ ( cos τ + sin τ ) and ∫ e τ sin τ d τ = 2 1 e τ ( sin τ − cos τ ) , both from 0 to t :
∫ 0 t e τ cos τ d τ = 2 1 ( e t ( cos t + sin t ) − 1 ) , ∫ 0 t e τ sin τ d τ = 2 1 ( e t ( sin t − cos t ) + 1 ) .
Step 5. Substitute and simplify (the e t sin t cos t cross-terms cancel):
2 1 e t − 2 1 cos t − 2 1 sin t
Verify by an independent partial-fraction inversion (no sign fudging). Write ( s − 1 ) ( s 2 + 1 ) 1 = s − 1 A + s 2 + 1 B s + C .
A (cover-up at s = 1 ): A = s 2 + 1 1 s = 1 = 2 1 .
Clear denominators: 1 = A ( s 2 + 1 ) + ( B s + C ) ( s − 1 ) . Match s 2 : 0 = A + B ⇒ B = − 2 1 . Match constant (s = 0 ): 1 = A − C ⇒ C = A − 1 = − 2 1 .
So ( s − 1 ) ( s 2 + 1 ) 1 = s − 1 1/2 + s 2 + 1 − 2 1 s − 2 1 . Inverting each piece with L − 1 { 1/ ( s − 1 )} = e t , L − 1 { s / ( s 2 + 1 )} = cos t , L − 1 { 1/ ( s 2 + 1 )} = sin t :
2 1 e t − 2 1 cos t − 2 1 sin t .
This matches the convolution result exactly — the cos t and sin t coefficients are both − 2 1 . ✓ (Answer grows , because of e t .)
L − 1 { ( s − a ) ( s − b ) 1 } for a = b , then take the limit b → a .
Forecast: two different decays e a t , e b t — expect their difference . What happens when they become equal? Guess before Step 4.
Step 1. f ( t ) = e a t , g ( t ) = e b t .
Why this step? Each factor 1/ ( s − a ) , 1/ ( s − b ) inverts to a pure exponential.
Step 2. Convolve:
∫ 0 t e a τ e b ( t − τ ) d τ = e b t ∫ 0 t e ( a − b ) τ d τ .
Why this step? Factor out e b t (constant in τ ); the remaining exponent is ( a − b ) τ .
Step 3. Integrate (valid because a = b , so the exponent isn't zero):
e b t ⋅ a − b e ( a − b ) t − 1 = a − b e a t − e b t .
a − b e a t − e b t ( a = b )
Step 4 — the degenerate case b → a . The formula is 0/0 ! Take the limit (or re-do Step 2 with b = a , giving ∫ 0 t 1 d τ = t ):
lim b → a a − b e a t − e b t = t e a t .
Why this step? When a = b the exponent ( a − b ) τ = 0 , so the integrand is just 1 and the integral is t , giving t e a t . This matches L − 1 { ( s − a ) 2 1 } = t e a t .
Verify (general): L { a − b e a t − e b t } = a − b 1 ( s − a 1 − s − b 1 ) = ( s − a ) ( s − b ) 1 . ✓
Verify (degenerate): L { t e a t } = ( s − a ) 2 1 . ✓
y ′′ − y = e t , y ( 0 ) = 0 , y ′ ( 0 ) = 0 (see Solving Linear ODEs with Laplace Transforms ).
Forecast: the forcing e t is resonant with the natural mode e t . Expect a t e t term.
Step 1. Laplace both sides using L { y ′′ } = s 2 Y − sy ( 0 ) − y ′ ( 0 ) :
s 2 Y − Y = s − 1 1 ⇒ ( s 2 − 1 ) Y = s − 1 1 .
Why this step? Zero initial conditions kill the boundary terms, turning the ODE into algebra.
Step 2. Isolate and factor: s 2 − 1 = ( s − 1 ) ( s + 1 ) , so
Y = ( s − 1 ) 2 ( s + 1 ) 1 = F ( s − 1 ) 2 1 ⋅ G s + 1 1 .
Why this step? Writing Y as a product invites the convolution theorem instead of a big partial fraction.
Step 3. Inverses: f ( t ) = t e t (from Ex 5's degenerate result), g ( t ) = e − t .
Step 4. Convolve:
y = ∫ 0 t τ e τ e − ( t − τ ) d τ = e − t ∫ 0 t τ e 2 τ d τ .
Why this step? Combine exponentials: e τ e − ( t − τ ) = e − t e 2 τ .
Step 5. Integrate by parts ∫ 0 t τ e 2 τ d τ = [ 2 τ e 2 τ − 4 1 e 2 τ ] 0 t = 2 t e 2 t − 4 1 e 2 t + 4 1 :
y = e − t ( 2 t e 2 t − 4 1 e 2 t + 4 1 ) = 2 t e t − 4 1 e t + 4 1 e − t .
y = 2 1 t e t − 4 1 e t + 4 1 e − t
Verify: compute y ′′ − y . With y = 2 1 t e t − 4 1 e t + 4 1 e − t , one finds y ′′ − y = e t , and y ( 0 ) = − 4 1 + 4 1 = 0 , y ′ ( 0 ) = 0 . ✓ (The resonant t e t appeared as forecast.)
y ( t ) = e − t + ∫ 0 t y ( τ ) ( t − τ ) d τ (see Volterra Integral Equations ).
Forecast: the integral is y ∗ t (convolving y with the ramp t ). Expect a mix of growing and decaying pieces.
Step 1. Recognise the hidden convolution: ∫ 0 t y ( τ ) ( t − τ ) d τ = ( y ∗ t ) ( t ) , and L { t } = s 2 1 .
Why this step? The kernel ( t − τ ) has argument t − τ — the tell-tale sign of a convolution with g ( t ) = t .
Step 2. Laplace the whole equation (L { y } = Y , L { e − t } = s + 1 1 ):
Y = s + 1 1 + Y ⋅ s 2 1 .
Why this step? Convolution → product turns the integral equation into ordinary algebra.
Step 3. Solve for Y :
Y ( 1 − s 2 1 ) = s + 1 1 ⇒ Y ⋅ s 2 s 2 − 1 = s + 1 1 ⇒ Y = ( s + 1 ) ( s 2 − 1 ) s 2 = ( s + 1 ) 2 ( s − 1 ) s 2 .
Why this step? Collect all Y on one side and divide; then factor s 2 − 1 = ( s − 1 ) ( s + 1 ) so the repeated root ( s + 1 ) 2 is visible — that shape dictates the partial-fraction template next.
Step 4. Partial fractions: Y = s − 1 A + s + 1 B + ( s + 1 ) 2 C .
Why this step? The repeated factor ( s + 1 ) 2 needs its own two terms (B for the simple power, C for the squared power).
A (cover-up at s = 1 ): A = ( s + 1 ) 2 s 2 s = 1 = 4 1 .
C (cover-up at s = − 1 ): C = s − 1 s 2 s = − 1 = − 2 1 = − 2 1 .
B (match, e.g. at s = 0 ): ( 1 ) 2 ( − 1 ) 0 = 0 = − 1 A + 1 B + 1 C = − 4 1 + B − 2 1 ⇒ B = 4 3 .
So A = 4 1 , B = 4 3 , C = − 2 1 .
Step 5. Invert term-by-term, using L − 1 { 1/ ( s + 1 ) 2 } = t e − t :
Why this step? Each simple piece has a known inverse; L − 1 { 1/ ( s − 1 )} = e t , L − 1 { 1/ ( s + 1 )} = e − t , and the squared factor gives t e − t .
y = 4 1 e t + 4 3 e − t − 2 1 t e − t .
y = 4 1 e t + 4 3 e − t − 2 1 t e − t
Verify: at t = 0 , y ( 0 ) = 4 1 + 4 3 − 0 = 1 , and the RHS at t = 0 is e 0 + 0 = 1 . ✓ The Laplace of this y equals ( s + 1 ) 2 ( s − 1 ) s 2 (checked below). ✓
Worked example Water is poured into a leaky tank at a
constant rate f ( t ) = 1 (litres/sec, starting at t = 0 ). A single litre added at time τ decays away as e − ( t − τ ) (the tank leaks proportionally). How much water is in the tank at time t ?
Forecast: early on it fills up; eventually inflow balances leak. What steady level do you expect? Guess.
Step 1. Model: the amount now is the sum over all past instants τ of "how much came in then" × "how much of it survives to now":
y ( t ) = ∫ 0 t = 1 f ( τ ) survival e − ( t − τ ) d τ = ( 1 ∗ e − t ) ( t ) .
Why this step? This is convolution — the survival kernel e − ( t − τ ) is the system's impulse response (see Transfer Functions and Impulse Response ). Note e − ( t − τ ) = 0 -behaviour: it only makes sense for t ≥ τ , matching the causal hypothesis and the [ 0 , t ] limits.
Step 2. Evaluate directly. Pull the τ -independent factor e − t out of the integral:
y ( t ) = ∫ 0 t e − ( t − τ ) d τ = e − t ∫ 0 t e τ d τ = e − t [ e τ ] 0 t = e − t ( e t − 1 ) = 1 − e − t .
Why this step? e − ( t − τ ) = e − t e τ ; only e τ depends on τ , so it is the only thing integrated.
y ( t ) = 1 − e − t
Step 3 — limiting behaviour (all cases). At t = 0 : y = 1 − e 0 = 0 (tank starts empty). As t → ∞ : e − t → 0 , so y → 1 litre — inflow rate exactly balances the leak, a stable steady state. For every finite t > 0 the level is strictly between 0 and 1 and rising. No scenario is left unshown.
Verify (via theorem): L { 1 } = s 1 , L { e − t } = s + 1 1 , product s ( s + 1 ) 1 = s 1 − s + 1 1 , inverse 1 − e − t . ✓ Units: integrand is (L/s)·(dimensionless survival)·(s) = L. ✓
t 3 ∗ cos t . Which factor should you flip to make the integral painless?
Forecast: flipping the polynomial forces messy integration-by-parts on t 3 ; flipping cos keeps the polynomial simple. Guess which is smarter.
Step 1. Commutativity says f ∗ g = g ∗ f , so we may write it either way:
t 3 ∗ cos t = ∫ 0 t τ 3 cos ( t − τ ) d τ = ∫ 0 t cos τ ( t − τ ) 3 d τ .
Why this step? From the parent proof, substituting u = t − τ swaps the two arguments — you are free to choose.
Step 2. Take the first form (flip cos ): ∫ 0 t τ 3 cos ( t − τ ) d τ . Rather than grind by parts, use the theorem to get the answer, then confirm.
L { t 3 } = s 4 6 , L { cos t } = s 2 + 1 s ⇒ L { t 3 ∗ cos t } = s 4 6 ⋅ s 2 + 1 s = s 3 ( s 2 + 1 ) 6 .
Why this step? Product in s -world = convolution in t -world; computing the transform is easier than the integral.
Step 3. Invert s 3 ( s 2 + 1 ) 6 . Partial fractions: s 3 ( s 2 + 1 ) 6 = s 3 6 − s 6 + s 2 + 1 6 s .
Why this step? The odd/even split of s 3 ( s 2 + 1 ) produces a low-power block 6/ s 3 − 6/ s plus a trig block 6 s / ( s 2 + 1 ) ; verified by clearing the common denominator (6 = 6 ( s 2 + 1 ) − 6 s 2 ( s 2 + 1 ) / s 0 … — machine-checked below).
Step 4. Invert term-by-term:
t 3 ∗ cos t = 3 t 2 − 6 + 6 cos t .
t 3 ∗ cos t = 3 t 2 − 6 + 6 cos t
Verify: at t = 0 the convolution integral is ∫ 0 0 ( ⋯ ) = 0 ; the formula gives 0 − 6 + 6 = 0 . ✓ And L { 3 t 2 − 6 + 6 cos t } = s 3 6 − s 6 + s 2 + 1 6 s = s 3 ( s 2 + 1 ) 6 . ✓
Lesson (the twist): either flip is legal , but converting to the s -product sidesteps all the by-parts pain. When an exam gives a polynomial times a trig/exponential, transform-then-invert usually beats brute-force convolution.
Recall Which cell is which?
Convolving with 1 = integrate (not identity) ::: cell C1
Product of two decays e a t , e b t gives a − b e a t − e b t , degenerating to t e a t ::: cell C5
A hidden ∫ 0 t y ( τ ) k ( t − τ ) d τ signals a Volterra equation solvable by Y ( s ) K ( s ) ::: cell C7
Leaky-tank steady level with unit inflow and e − t leak ::: 1 litre (Ex 8, cell C8)
Product F times G in s world
Convolution blend in t world
Convolve with 1 equals integrate
Two exponentials difference quotient
Degenerate limit gives t times e
ODE forcing solved by blend
Volterra equation becomes algebra
Impulse response of a real system