ODE: Impulse Functions and Delta Forcing

Engineering systems often experience forces that act over an extremely short duration—a hammer striking a mass, a voltage spike in a circuit, or a sudden impact in a structure. Modeling these instantaneous events with ordinary functions is impossible, as it would require an infinite force over an infinitesimal time. This conceptual and mathematical challenge is resolved by the Dirac delta function, a generalized function that provides a powerful framework for modeling impulsive forcing in ordinary differential equations (ODEs). Mastering this tool allows you to analyze systems subject to shocks, understand signal sampling, and solve a broader class of initial value problems efficiently.

The Dirac Delta Function: A Generalized Concept

The Dirac delta function, denoted $\delta (t)$, is not a function in the classical sense but a distribution or generalized function. It is defined by two fundamental properties. First, it is zero everywhere except at the origin: $\delta (t)=0$ for $t\ne 0$. Second, its integral over the entire real line is one: ${\int }_{-\infty }^{\infty }\delta (t)dt=1$. You can think of it as the limit of a sequence of pulses (like tall, narrow rectangles or Gaussian curves) as their width approaches zero while their area remains fixed at one.

The most critical operational property for solving ODEs is the sifting property. For a continuous function $f(t)$, the integral ${\int }_{-\infty }^{\infty }f(t)\delta (t-a)dt=f(a)$. The delta function $\delta (t-a)$ "sifts out" the value of $f$ at the point $t=a$. This property is the workhorse for handling delta forcing in integrals and transforms. Physically, if $\delta (t-a)$ represents an impulse at time $t=a$, then its integral, ${\int }_{-\infty }^t\delta (\tau -a)d\tau$, is the unit step function $u(t-a)$, which jumps from 0 to 1 at $t=a$.

Laplace Transform of the Delta Function

The Laplace transform is the most effective technique for solving linear ODEs with impulsive forcing. Applying the sifting property directly yields the transform. For a delta function shifted to time $t=a\ge 0$: $$\mathcal{L}\{\delta (t-a)\}={\int }_0^{\infty }{e}^{-st}\delta (t-a)dt=e^{-as}.$$ The key case is for an impulse at $t=0^{+}$: $$\mathcal{L}\{\delta (t)\}=1.$$ This simple result is profound: in the Laplace domain, an instantaneous unit impulse is represented by the constant 1. This transforms a differential equation with a delta forcing term into an algebraic equation that is straightforward to solve.

Solving ODEs with Impulsive Forcing

Consider a standard second-order system, such as a spring-mass-damper, struck by a hammer at $t=a$. The ODE is: $$m{y}^{\prime \prime }+c{y}^{\prime }+ky=\delta (t-a),$$ with given initial conditions, often rest conditions $y(0)=0,y^{\prime }(0)=0$. The solution workflow using Laplace transforms is methodical:

1. Take the Laplace transform of the entire equation:

$$m[s^{2}Y(s)-sy(0)-y^{\prime }(0)]+c[sY(s)-y(0)]+kY(s)=e^{-as}.$$

1. Apply the initial conditions and solve for $Y(s)$:

$$Y(s)=\frac{e^{-as}}{m{s}^2+cs+k}.$$

1. The term $e^{-as}$ in the $s$-domain corresponds to a time shift. Let $G(s)=1/(m{s}^2+cs+k)$ be the transfer function. Its inverse transform, $g(t)={\mathcal{L}}^{-1}\{G(s)\}$, is the impulse response of the system—the solution when the impulse occurs at $t=0$.
2. By the time-shift property, the final solution is:

$$y(t)=g(t-a)u(t-a).$$ The solution is zero up until $t=a$, and for $t>a$, it is precisely the system's natural impulse response, initiating from the moment of impact. The impulse instantly imparts energy to the system, which then evolves according to its inherent dynamics (damped oscillation, pure oscillation, or exponential decay).

Physical Interpretation and Momentum Transfer

The physical interpretation of the Dirac delta is best understood through mechanics. Newton's second law states $F=ma=m\frac{dv}{dt}$, or $Fdt=mdv$. The integral of force over time is the impulse, which equals the change in momentum. A force $F(t)$ modeled as $\delta (t-a)$ has an impulse of 1: $$\text{Impulse}=\int F(t)dt=\int \delta (t-a)dt=1.$$ Applying this to $m\frac{dv}{dt}=\delta (t-a)$ and integrating over an infinitesimal interval around $t=a$ shows the velocity jumps discontinuously: $\Delta v=v(a^{+})-v(a^{-})=1/m$. The position $y(t)$, however, remains continuous. The delta function thus models an instantaneous transfer of momentum (or an equivalent sudden change in velocity) without a finite displacement. This is the essence of an ideal impact.

Applications: Impact Problems and Sampling Theory

The two primary application domains are mechanical impact and signal processing. In impact problems, such as a particle striking a barrier or a hammer hitting a nail, the delta function idealizes the contact force. The analysis provides the resulting motion or stress waves. For example, determining the subsequent vibrations of a beam after being tapped is a direct application.

In sampling theory, foundational to digital signal processing, the delta function is used to represent the ideal sampling of a continuous signal $f(t)$. The sampled signal is a train of impulses: $f_s(t)={\sum }_{n=-\infty }^{\infty }f(nT)\delta (t-nT)$, where $T$ is the sampling interval. Each delta function sifts out the signal's value at the sampling instant $t=nT$. This model allows engineers to analyze the frequency content of sampled signals and derive critical results like the Nyquist-Shannon sampling theorem.

Common Pitfalls

1. Treating $\delta (t)$ as a pointwise function: A common error is evaluating $\delta (0)$ as "infinity." The delta function is only meaningful under an integral sign via the sifting property. Remember, it's a distribution, not a real-valued function.
2. Misapplying initial conditions with impulses at $t=0$: When an impulse $\delta (t)$ occurs at the initial time, the standard initial conditions $y(0),y^{\prime }(0)$ refer to the instant just before the impulse, $t=0^{-}$. The impulse causes an instantaneous change. The correct method is to take the Laplace transform using the initial conditions at $t=0^{-}$, which are typically given. The transform automatically incorporates the jump.
3. Forgetting the time shift in the final solution: After finding $Y(s)=e^{-as}G(s)$, the inverse transform is $g(t-a)u(t-a)$, not simply $g(t)$. The system is quiescent until the impulse occurs at $t=a$.
4. Ignoring system linearity: The superposition principle holds. The solution for multiple impulses, e.g., $\delta (t-a)+2\delta (t-b)$, is the sum of the individual time-shifted impulse responses: $y(t)=g(t-a)u(t-a)+2g(t-b)u(t-b)$.

Summary

• The Dirac delta function $\delta (t-a)$ is a mathematical idealization of a unit impulse occurring instantaneously at time $t=a$, defined by its sifting property: $\int f(t)\delta (t-a)dt=f(a)$.
• Its Laplace transform is $\mathcal{L}\{\delta (t-a)\}=e^{-as}$, simplifying the solution of linear ODEs with impulsive forcing to an algebraic process in the $s$-domain.
• Solving an ODE like $m{y}^{\prime \prime }+c{y}^{\prime }+ky=\delta (t-a)$ yields a solution $y(t)=g(t-a)u(t-a)$, where $g(t)$ is the system's impulse response, capturing its natural dynamics triggered at $t=a$.
• Physically, a delta force causes a discontinuous jump in velocity (a momentum change of 1 unit) while leaving position continuous, perfectly modeling idealized instantaneous impacts.
• Key applications span mechanical impact analysis (e.g., vibration initiation) and signal sampling theory, where a train of delta functions is used to mathematically represent the conversion of a continuous signal into a discrete sequence.