Historical Context & Motivation
In the real world, forces and inputs don't always change smoothly. Imagine flipping a light switch — the voltage jumps instantly from zero to full power. Or think about a hammer striking a nail — the force is enormous for a tiny fraction of a second, then disappears. These kinds of sudden, discontinuous inputs show up everywhere in engineering, physics, and biology. Standard methods for solving ordinary differential equations (ODEs) struggle with these abrupt changes, because those techniques typically assume the input function is smooth and well-behaved.
Mathematicians and engineers developed special tools to handle these situations. The Heaviside step function models sudden on/off switches, while the Dirac delta function models instantaneous impulses. Combined with Laplace transforms, these tools convert messy, piecewise differential equations into straightforward algebra problems. The timeline below shows how these ideas came together.
The central question this lesson addresses is: How do we solve a differential equation when the input suddenly switches on, switches off, or delivers an instantaneous shock? The answer lies in translating these tricky inputs into the Laplace domain, where they become simple algebraic expressions.
Core Principles & Definitions
Before diving into calculations, you need to understand four foundational ideas that make everything else in this lesson work. These concepts build on one another: the step function models switches, the delta function models impulses, Laplace transforms convert them to algebra, and the shifting theorem handles delayed inputs.
The Heaviside Step Function u(t − a)
The Dirac Delta Function δ(t − a)
The Laplace Transform
The Second Shifting Theorem
Visualizing Discontinuous & Impulsive Inputs
To truly understand these functions, you need to see what they look like. The diagram below shows the Heaviside step function u(t − 2), which stays at 0 until t = 2 and then jumps to 1, alongside the Dirac delta function δ(t − 2), which is an infinitely tall, infinitely thin spike at t = 2. Both are plotted on the same time axis so you can compare their behavior.
Notice how different these two functions are. The step function represents a permanent change — once the switch flips on, it stays on. The delta function represents a momentary jolt — all the energy is concentrated at a single instant. In practice, you can build any piecewise-defined input by combining step functions. For example, a force that's on from t = 2 to t = 5 can be written as u(t − 2) − u(t − 5). This is why these building blocks are so powerful.
Mathematical Framework
Now let's formalize the mathematics. The strategy for solving ODEs with discontinuous or impulsive inputs follows a consistent three-step process: (1) take the Laplace transform of both sides, (2) solve the resulting algebraic equation for Y(s), and (3) apply the inverse Laplace transform to get y(t). The key formulas you need are listed below.
Building Piecewise Functions with Step Functions
One of the most practical skills in this topic is learning to rewrite a piecewise-defined function using step functions. In standard math notation, you might define a function with cases — for example, g(t) = 0 for t < 1, g(t) = 5 for 1 ≤ t < 4, and g(t) = 0 for t ≥ 4. Using step functions, this becomes g(t) = 5[u(t − 1) − u(t − 4)]. The step u(t − 1) switches the 5 on at t = 1, and subtracting u(t − 4) switches it back off at t = 4.
This technique generalizes naturally. Any function that changes its formula at specific times can be rewritten using step functions. For instance, suppose a force ramps up linearly starting at t = 3. You'd write f(t) = (t − 3) · u(t − 3). The factor (t − 3) gives the ramp shape, and u(t − 3) ensures it only starts at t = 3. This form is exactly what the second shifting theorem is designed to handle.
| Piecewise Input | Step-Function Form | Laplace Transform |
|---|---|---|
| Constant c from t = a onward | c · u(t − a) | c · e^(−as) / s |
| Constant c from t = a to t = b | c · [u(t − a) − u(t − b)] | (c/s) · [e^(−as) − e^(−bs)] |
| Ramp (t − a) starting at t = a | (t − a) · u(t − a) | e^(−as) / s² |
| Impulse at t = a | δ(t − a) | e^(−as) |
Worked Example — ODE with a Step-Function Input
Let's solve a complete problem from start to finish. Consider a spring-mass system where a constant force of 5 is suddenly applied at t = 2. The governing equation is:
Strengths & Limitations of the Laplace Transform Method
The Laplace transform approach is not the only way to solve ODEs, so it's important to understand when it shines and when other methods might be preferred. The table below compares the Laplace method with one alternative you may encounter in a differential equations course: the method of undetermined coefficients, which is used for smooth (non-discontinuous) inputs.
| Feature | Laplace Transforms | Undetermined Coefficients |
|---|---|---|
| Discontinuous inputs | Handles naturally via step/delta functions | Cannot handle directly; must split into intervals |
| Impulsive inputs (delta) | Simple: ℒ{δ(t − a)} = e^(−as) | Not applicable |
| Initial conditions | Built into the transform automatically | Applied separately after finding general solution |
| Smooth, polynomial/exponential inputs | Works but may be more work than needed | Usually faster and more direct |
| Algebraic complexity | Partial fraction decomposition can get tedious | System of equations for coefficients |
Connections to Transfer Functions & Advanced Theory
The techniques you've learned in this lesson are the gateway to deeper topics in engineering and applied mathematics. When you take the Laplace transform of an ODE and solve for Y(s), the ratio Y(s)/G(s) — where G(s) is the Laplace transform of the input — is called the transfer function H(s). The transfer function captures everything about how the system responds to inputs, independent of what the actual input is — it is a property of the system itself, not of any particular input. This idea is foundational in control systems engineering, signal processing, and circuit analysis, and is a major topic in upper-division engineering courses.
| This Lesson | Advanced Extension |
|---|---|
| Step function u(t − a) | Step response of systems; characterizing system behavior |
| Delta function δ(t − a) | Impulse response h(t); convolution integral y(t) = h(t) * g(t) |
| Solving Y(s) algebraically | Transfer functions H(s) = Y(s)/G(s); Bode plots; stability analysis |
| Second shifting theorem | Delay differential equations; time-delay systems in control |
If this topic interests you, the natural next step is studying the convolution theorem, which states that ℒ{f * g} = F(s)·G(s). Here, f * g denotes the convolution integral — a way of combining two functions that measures how much one function overlaps with a reversed and shifted copy of the other. This powerful result means that the output of any linear system can be computed by convolving the input with the system's impulse response — the response to δ(t). You've already taken the first step toward understanding this by working with δ(t) in this lesson.
Practice Problems
Lesson Summary
This lesson introduced how to solve ordinary differential equations with discontinuous or impulsive inputs using Laplace transforms — a core technique from undergraduate differential equations and engineering mathematics. We saw that the Heaviside step function u(t − a) models sudden switches that turn on at time a, while the Dirac delta function δ(t − a) models instantaneous impulses. Their Laplace transforms — e^(−as)/s and e^(−as), respectively — are elegantly simple. By combining these functions, any piecewise input can be expressed in a form suitable for the Laplace approach.
The solution process follows three steps: take the Laplace transform of both sides (converting derivatives to algebra), solve for Y(s) using algebra and partial fraction decomposition, then apply the inverse Laplace transform with the second shifting theorem to recover y(t). This method automatically incorporates initial conditions and handles jumps and spikes that would be extremely difficult to manage with other techniques. Looking ahead, these ideas lead to transfer functions and convolution — cornerstones of modern engineering analysis typically studied in upper-division college courses.