Q: How do you define the derivatives of the Heaviside, Sign, Absolute Value, and Delta functions? How do they relate to one another?

Physicist: These are four standard reference functions.  In the same way that there are named mathematical constants, like π or e, there are named mathematical functions.  These are among the more famous (after the spotlight hogging trig functions).

The Sign, Delta, Absolute Value, and Heaviside functions.

The Sign, Delta, Absolute Value, and Heaviside functions.  The graphs on top are the slope of the graphs on the bottom (and slope=derivative).

The absolute value function flips the sign of negative numbers, and leaves positive numbers alone.  The sign function is 1 for positive numbers and -1 for negative numbers.  The Heaviside function is very similar; 1 for positive numbers and 0 for negative numbers.  By the way, the Heaviside function, rather than being named after its shape, is named after Oliver Heaviside, who was awesome.

The delta function is a whole other thing.  The delta function is zero everywhere other than at x=0 and at x=0 it’s infinite but there’s “one unit of area” under that spike.  Technically the delta function isn’t a function because it can’t be defined at zero.  The “Dirac delta function” is used a lot in physics (Dirac was a physicist) to do things like describe the location of the charge of single particles.  An electron has one unit of charge, but it’s smaller than basically anything, so describing it as a unit of charge located in exactly one point usually works fine (and if it doesn’t, don’t use a delta function).  This turns out to be a lot easier than modeling a particle as… just about anything else.

The derivative of a function is the slope of that function.  So, the derivative of |x| is 1 for positive numbers (45° up), and -1 for negative numbers (45° down).  But that’s the sign function!  Notice that at x=0, |x| has a kink and the slope can’t be defined (hence the open circles in the graph of sgn(x)).

The derivative of the Heaviside function is clearly zero for x≠0 (it’s completely level), but weird stuff happens at x=0.  There, if you were to insist that somehow the slope exists, you would find that no finite number does the job (vertical lines are “infinitely steep”).  But that sounds a bit like the delta function; zero everywhere, except for an infinite spike at x=0.

It is possible (even useful!) to define the delta function as δ(x) = H'(x).  Using that, you find that sgn'(x) = 2δ(x), simply because the jump is twice the size.  However, how you define derivatives for discontinuous functions is a whole thing, so that’ll be left in the answer gravy.

Answer Gravy: The Dirac delta function really got under the skin of a lot of mathematicians.  Many of them flatly refuse to even call it a “function” (since technically it doesn’t meet all the requirements).  Math folk are a skittish bunch, and when a bunch of handsome/beautiful pioneers (physicists) are using a function that isn’t definable, mathematicians can’t help but be helpful.  When they bother, physicists usually define the delta function as the limit of a series of progressively thinner and taller functions (usually Gaussians).

One of the simplest ways to construct the delta function is the series of functions shown here (and infinitely many more).

One of the simplest ways to construct the delta function is the series of functions fn(x) = n, 0 ≤ x ≤ 1/n (the first four of which are shown).  The area under each of these is 1, and most/many of the important properties of delta functions can be derived by looking at the limit as n→∞.

Mathematicians take a different tack.  For those brave cognoscenti, the delta function isn’t a function at all; instead it’s a “distribution”, which is a member of the dual of function space, and it’s used to define a “bounded linear functional”.

So that’s one issue cleared up.

A “functional” takes an entire function as input, and spits out a single number as output.  When you require that the functional is linear (and why not?), you’ll find that the only real option is for the functional to take the form F(g) = \int f(x)g(x)\,dx.  This is because of the natural linearity of the integral:

\begin{array}{ll}F(g + h) \\= \int f(x)\left(g(x)+h(x)\right)\,dx \\= \int f(x)g(x)\,dx + \int f(x)h(x)\,dx \\= F(g) + F(h)\end{array}

In F(g) = \int f(x)g(x)\,dx, F is the functional, f(x) is the distribution corresponding to that functional, and g(x) is the function being acted upon.  The delta function is the distribution corresponding to the functional which simply returns the value at zero.  That is, \int \delta(x)g(x)\,dx = g(0).  So finally, what in the crap does “returning the value at zero” have to do with the derivative of the Heaviside function?  As it happens: buckets!

Assume here that A<0<B,

\begin{array}{ll}    \int_A^B \delta(x)g(x)\,dx \\[2mm]    = H(x)g(x)\big|_A^B - \int_A^B H(x)g^\prime(x)\,dx & \textrm{(integration by parts)} \\[2mm]    = H(B)g(B) - H(A)g(A) - \int_A^B H(x)g^\prime(x)\,dx \\[2mm]    = H(B)g(B) - \int_0^B H(x)g^\prime(x)\,dx & (H(x)=0, x<0) \\[2mm]    = g(B) - \int_0^B g^\prime(x)\,dx & (H(x)=1, x>0) \\[2mm]    = g(B) - g(x)\big|_0^B & \textrm{(fundamental theorem of calculus)} \\[2mm]    = g(B) - \left[g(B) - g(0) \right] \\[2mm]    =g(0)    \end{array}

Running through the same process again, you’ll find that this is a halfway decent way of going a step further and defining the derivative of the delta function, δ'(x).

\begin{array}{ll}    \int_A^B \delta^\prime(x)g(x)\,dx \\[2mm]    = \delta(x)g(x)\big|_A^B - \int_A^B \delta(x)g^\prime(x)\,dx \\[2mm]    = - \int_A^B \delta(x)g^\prime(x)\,dx \\[2mm]    = -g^\prime(0)    \end{array}

δ'(x) also isn’t a function, but is instead another profoundly abstract distribution.  And yes: this can be done ad nauseam (or at least ad queasyam) to create distributions that grab higher and higher derivatives of input functions.

This entry was posted in -- By the Physicist, Equations, Math. Bookmark the permalink.

3 Responses to Q: How do you define the derivatives of the Heaviside, Sign, Absolute Value, and Delta functions? How do they relate to one another?

  1. Flavian Popa says:

    Thanks for posting. One question, regarding the derivative of Dirac function, we actually cannot explicitly define it in every point of the interval [A;B], we just define it based on the integral relation involving δ’(x) and g(x), right? I mean this is what makes distributions so “abstract” mathematical concepts..if I am not mistaking
    Anyway great thanks again for amazing Maths articles posting.

  2. The Physicist The Physicist says:

    That’s exactly how it’s so damn abstract.
    If you describe δ(x) as the limit of a string of functions, fn(x), then you can try to describe δ’(x) as the limit of the string of functions fn‘(x). When you do this you find that δ’(x) is an infinite positive spike immediately followed by, and at the same location as, an infinite negative spike. Somehow this seems even worse than the difficulties with the original delta function.
    Easiest to stick to the (more useful/consistent) abstraction.

  3. Andres says:

    You can find the derivative of absolute value of x by using some kind of substitution. Notice that |x| = sqrt(x^2). Differentiating the second one and using chain rule we have:

    2x/2sqrt(x^2)=x/sqrt(x^2). As |x| = sqrt(x^2).


    Notice is equal to 1 for all x>0, -1 for all x<0 and that is not defined at x=0. (Not differentiable at x=0).

Leave a Reply

Your email address will not be published. Required fields are marked *