6.4 Dirac Delta and Impulse Response

Rectangular Pulse

Often in applications we study a physical system by putting in a short pulse and then seeing what the system does. The resulting behavior is often called impulse response. Let us see what we mean by a pulse. The simplest kind of a pulse is a simple rectangular pulse defined by

$0 i f t less than a,; M i f a less than or equal to t less than b,; 0 i f b less than or equal to t.$

Notice that

$φ (t) = M (u (t - a) - u (t - b)),$

where

u (t)

is the unit step function (see Figure

1

for a graph).

Graph that is on the x-axis from (0,0) to (0,0.5), goes up to (0,2), then horizontal to (1,2) then down to the x-axis again. — Figure $1$ : Sample square pulse with $a equals 0.5$ , $b equals 1$ , and $M equals 2$ .

Let us take the Laplace transform of a square pulse,

$\begin{matrix} ℒ {φ (t)} = ℒ {M (u (t - a) - u (t - b))} = M \frac{e^{- a s} - e^{- b s}}{s} . \end{matrix}$

For simplicity we let $a equals 0$ and it is convenient to set $M = \frac{1}{b}$ to have

$∫ 0 infinity φ open parenthesis t close parenthesis d t equals 1$

That is, to have the pulse have “unit mass.” For such a pulse we compute

$ℒ {φ (t)} = ℒ {\frac{u (t) - u (t - b)}{b}} = \frac{1 - e^{- b s}}{b s} .$

We generally want $b$ to be very small. That is, we wish to have the pulse be very short and very tall. By letting $b$ go to zero we arrive at the concept of the Dirac delta function.

Function graph showing y = (1/(2*b))*(x/abs(x) - (x-b)/abs(x-b)). Adjustable parameter: Pulse width b (b) = 0.5. Viewing window: x from -2.27 to 7.47, y from -0.51 to 5.51. — The book's own normalized pulse φ(t) = (u(t) − u(t−b))/b — the a = 0, M = 1/b case this section sets up right before it takes the limit — with time t on the horizontal axis. Drag b from 4 down to 0.2 and watch the rectangle go short-and-wide → tall-and-thin: the height 1/b is not chosen, it is forced by the area staying exactly 1, which is precisely the limit that defines δ(t). At b = 0.5 the rectangle is 0.5 wide and 2 tall — the same shape as the book's own sample square pulse (a = 0.5, b = 1, M = 2), slid back to start at the origin. There are no vertical sides at t = 0 and t = b: those are genuine jump discontinuities, not gaps in the drawing.

6.4.2Delta Function

The Dirac delta functionto the power 1 is not exactly a function; it is sometimes called a generalized function. We avoid unnecessary details and simply say that it is an object that does not really make sense unless we integrate it. The motivation is that we would like a “function” $δ (t)$ such that for any continuous function $f (t)$ we have

$\int_{- \infty}^{\infty} δ (t) f (t) d t = f (0)$

The formula should hold if we integrate over any interval that contains 0, not just $(- \infty, \infty)$ . So $δ (t)$ is a “function” with all its “mass” at the single point $t equals 0$ . In other words, for any interval $[ c, d ]$

$1 i f t h e i n t e r v a l [ c, d ] c o n t a i n s 0, i. e. c less than or equal to 0 less than or equal to d,; 0 o t h e r w i s e.$

Unfortunately there is no such function in the classical sense. You could informally think that $δ (t)$ is zero for $t not equal to 0$ and somehow infinite at $t equals 0$ .

A good way to think about $δ (t)$ is as a limit of short pulses whose integral is 1. For example, suppose that we have a square pulse $φ (t)$ as above with $a equals 0$ , $M = \frac{1}{b}$ , that is $φ (t) = \frac{u (t) - u (t - b)}{b}$ .

Compute

$\int_{- \infty}^{\infty} φ (t) f (t) d t = \int_{- \infty}^{\infty} \frac{u (t) - u (t - b)}{b} f (t) d t = \frac{1}{b} \int_{0}^{b} f (t) d t .$

If $f (t)$ is continuous at $t equals 0$ , then for very small $b$ , the function $f (t)$ is approximately equal to $f open parenthesis 0 close parenthesis$ on the interval $[ 0, b close parenthesis$ . We approximate the integral

$\frac{1}{b} \int_{0}^{b} f (t) d t \approx \frac{1}{b} \int_{0}^{b} f (0) d t = f (0) .$

Therefore,

${lim}_{b \to 0} \int_{- \infty}^{\infty} φ (t) f (t) d t = {lim}_{b \to 0} \frac{1}{b} \int_{0}^{b} f (t) d t = f (0)$

Let us therefore accept $δ (t)$ as an object that is possible to integrate. We often want to shift $δ$ to another point, for example $δ (t - a)$ . In that case we have

$\int_{- \infty}^{\infty} δ (t - a) f (t) d t = f (a) .$

Note that $δ (a - t)$ is the same object as $δ (t - a)$ . In other words, the convolution of $δ (t)$ with $f (t)$ is again $f (t)$ ,

$(f * δ) (t) = \int_{0}^{t} δ (t - s) f (s) d s = f (t)$

As we can integrate $δ (t)$ , let us compute its Laplace transform.

$ℒ {δ (t - a)} = \int_{0}^{\infty} e^{- s t} δ (t - a) d t = e^{- a s} .$

In particular,

$ℒ { δ open parenthesis t close parenthesis } equals 1.$

Note

Notice that the Laplace transform of $δ (t - a)$ looks like the Laplace transform of the derivative of the Heaviside function $u (t - a)$ , if we could differentiate the Heaviside function. First notice

$ℒ {δ (t - a)} = \frac{e^{- a s}}{s} .$

To obtain what the Laplace transform of the derivative would be we multiply by $s$ , to obtain $e^{- a s}$ , which is the Laplace transform of $δ (t - a)$ . We see the same thing using integration,

$\int_{0}^{t} δ (s - a) d s = u (t - a)$

So in a certain sense

$\frac{d}{d t} [u (t - a)] = δ (t - a)$

This line of reasoning allows us to talk about derivatives of functions with jump discontinuities. We can think of the derivative of the Heaviside function $u (t - a)$ as being somehow infinite at $a$ , which is precisely our intuitive understanding of the delta function.

Example 1

Compute

$ℒ^{- 1} {\frac{s + 1}{s}} .$

So far we have always looked at proper rational functions in the $s$ variable. That is, the numerator was always of lower degree than the denominator. Not so with $\frac{s + 1}{s}$ . We write,

$ℒ^{- 1} {\frac{s + 1}{s}} = ℒ^{- 1} {1 + \frac{1}{s}} = ℒ^{- 1} {1} + ℒ^{- 1} {\frac{1}{s}} = δ (t) + 1 .$

The resulting object is a generalized function and only makes sense when put underneath an integral.

Impulse Response

As we said before, in the differential equation $L x = f (t)$ , we think of $f (t)$ as input, and $x (t)$ as the output. Often it is important to find the response to an impulse, and then we use the delta function in place of $f (t)$ . The solution to

$L x = δ (t)$

is called the impulse response.

Example 2

Solve (find the impulse response)

$x^{″} + ω_{0}^{2} x = δ (t), x (0) = 0, x^{'} (0) = 0 .$

(6.4.1)

We first apply the Laplace transform to the equation. Denote the transform of $x (t)$ by $X (s)$ .

$s^{2} X (s) + ω_{0}^{2} X (s) = 1, and so X (s) = \frac{1}{s^{2} + ω_{0}^{2}} .$

Taking the inverse Laplace transform we obtain

$x (t) = \frac{\sin (ω_{0} t)}{ω_{0}} .$

Let us notice something about the above example. In Example 6.3.4, we found that when the input was $f (t)$ , then the solution to $L x = f (t)$ was given by

$x (t) = \int_{0}^{t} f (τ) \frac{\sin (ω_{0} (t - τ))}{ω_{0}} d τ .$

Notice that the solution for an arbitrary input is given as convolution with the impulse response. Let us see why. The key is to notice that for functions $x (t)$ and $f (t)$ ,

$(x * f)^{″} (t) = \frac{d^{2}}{d t^{2}} [\int_{0}^{t} f (τ) x (t - τ) d τ] = \int_{0}^{t} f (τ) x^{″} (t - τ) d τ = (x^{″} * f) (t) .$

We simply differentiate twice under the integral,squared the details are left as an exercise. And so if we convolve the entire equation (6.4.1), the left hand side becomes

$(x^{″} + ω_{0}^{2} x) * f = (x^{″} * f) + ω_{0}^{2} (x * f) = (x * f)^{″} + ω_{0}^{2} (x * f) .$

The right hand side becomes

$open parenthesis δ * f close parenthesis open parenthesis t close parenthesis equals f open parenthesis t close parenthesis.$

Therefore $y open parenthesis t close parenthesis equals open parenthesis x * f close parenthesis open parenthesis t close parenthesis$ is the solution to

$y^{″} + ω_{0}^{2} y = f (t) .$

This procedure works in general for other linear equations $L x = f (t)$ . If you determine the impulse response, you also know how to obtain the output $x (t)$ for any input $f (t)$ by simply convolving the impulse response and the input $f (t)$ .

Three-Point Beam Bending

Let us give another quite different example where the delta function turns up: Representing point loads on a steel beam. Consider a beam of length $L$ , resting on two simple supports at the ends. Let $x$ denote the position on the beam, and let $y (x)$ denote the deflection of the beam in the vertical direction. The deflection $y (x)$ satisfies the Euler-Bernoulli equation,cubed

$E I \frac{d^{4} y}{d x^{4}} = F (x),$

where

E

and

I

are constantsto the power 4 and

F (x)

is the force applied per unit length at position

x

. The situation we are interested in is when the force is applied at a single point as in Figure

2

Picture of a bar that is pushed up on the ends and down in the center. — Figure $2$ : Three-point bending.

In this case the equation becomes

$E I \frac{d^{4} y}{d x^{4}} = - F δ (x - a)$

where $x = a$ is the point where the mass is applied. $F$ is the force applied and the minus sign indicates that the force is downward, that is, in the negative $y$ direction. The end points of the beam satisfy the conditions,

$\begin{matrix} y (0) = 0, y^{″} (0) = 0, \\ y (L) = 0, y^{″} (L) = 0 . \end{matrix}$

See Section 5.2, for further information about endpoint conditions applied to beams.

Example 3

Suppose that length of the beam is $2$ , and suppose that $E I equals 1$ for simplicity. Further suppose that the force $F equals 1$ is applied at $x equals 1$ . That is, we have the equation

$\frac{d^{4} y}{d x^{4}} = - δ (x - 1)$

and the endpoint conditions are

$y (0) = 0, y^{″} (0) = 0, y (2) = 0, y^{″} (2) = 0 .$

We could integrate, but using the Laplace transform is even easier. We apply the transform in the $x$ variable rather than the $t$ variable. Let us again denote the transform of $y (x)$ as $Y (s)$ .

$s^{4} Y (s) - s^{3} y (0) - s^{2} y^{'} (0) - s y^{″} (0) - y^{‴} (0) = - e^{- s} .$

We notice that $y open parenthesis 0 close parenthesis equals 0$ and $y double prime open parenthesis 0 close parenthesis equals 0$ . Let us call $C sub 1 equals y prime open parenthesis 0 close parenthesis$ and $C^{2} = y^{‴} (0)$ . We solve for $Y (s)$ ,

$Y (s) = \frac{- e^{- s}}{s^{4}} + \frac{C_{1}}{s^{2}} + \frac{C_{2}}{s^{4}} .$

We take the inverse Laplace transform utilizing the second shifting property Equation (6.2.14) to take the inverse of the first term.

$y (x) = \frac{- (x - 1)^{3}}{6} u (x - 1) + C_{1} x + \frac{C_{2}}{6} x^{3} .$

We still need to apply two of the endpoint conditions. As the conditions are at $x equals 2$ we can simply replace $u (x - 1) = 1$ when taking the derivatives. Therefore,

$0 = y (2) = \frac{- (2 - 1)^{3}}{6} + C_{1} (2) + \frac{C_{2}}{6} 2^{3} = \frac{- 1}{6} + 2 C_{1} + \frac{4}{3 C_{2}},$

and

$0 = y^{″} (2) = \frac{- 3 \cdot 2 \cdot (2 - 1)}{6} + \frac{C_{2}}{6} 3 \cdot 2 \cdot 2 = - 1 + 2 C_{2}$

Hence $C_{2} = \frac{1}{2}$ and solving for $C sub 1$ using the first equation we obtain $C_{1} = \frac{- 1}{4}$ . Our solution for the beam deflection is

$y (x) = \frac{- {(x - 1)}^{3}}{6} u (x - 1) - \frac{x}{4} + \frac{x^{3}}{12} .$

Footnotes

[1] Named after the English physicist and mathematician Paul Adrien Maurice Dirac (1902–1984).

[2] You should really think of the integral going over $(- \infty, \infty)$ rather than over $[ 0, t ]$ and simply assume that $f (t)$ and $x (t)$ are continuous and zero for negative.

[3] Named for the Swiss mathematicians Jacob Bernoulli (1654–1705), Daniel Bernoulli —nephew of Jacob— (1700–1782), and Leonhard Paul Euler (1707–1783).

[4] $E$ is the elastic modulus and $I$ is the second moment of area. Let us not worry about the details and simply think of these as some given constants.

Adapted from Differential Equations for Engineers by Jiří Lebl (Oklahoma State University), hosted on LibreTexts (math.libretexts.org) and licensed under CC BY-SA 4.0. Changes were made. License: CC-BY-SA-4.0.