Heaviside step function

${\displaystyle {\tfrac {1}{2}}+{\tfrac {1}{2}}\tanh(kx)={\frac {1}{1+e^{-2kx}}}}$
approaches the step function as k → ∞.

For a smooth approximation to the step function, one can use the logistic function

$${\displaystyle H(x)\approx {\tfrac {1}{2}}+{\tfrac {1}{2}}\tanh kx={\frac {1}{1+e^{-2kx}}},}$$

where a larger k corresponds to a sharper transition at x = 0. If we take H(0) = 1/2, equality holds in the limit:

$${\displaystyle H(x)=\lim _{k\to \infty }{\tfrac {1}{2}}(1+\tanh kx)=\lim _{k\to \infty }{\frac {1}{1+e^{-2kx}}}.}$$

There are many other smooth, analytic approximations to the step function.[2] Among the possibilities are:

$${\displaystyle {\begin{aligned}H(x)&=\lim _{k\to \infty }\left({\tfrac {1}{2}}+{\tfrac {1}{\pi }}\arctan kx\right)\\H(x)&=\lim _{k\to \infty }\left({\tfrac {1}{2}}+{\tfrac {1}{2}}\operatorname {erf} kx\right)\end{aligned}}}$$

These limits hold pointwise and in the sense of distributions. In general, however, pointwise convergence need not imply distributional convergence, and vice versa distributional convergence need not imply pointwise convergence. (However, if all members of a pointwise convergent sequence of functions are uniformly bounded by some "nice" function, then convergence holds in the sense of distributions too.)

In general, any cumulative distribution function of a continuous probability distribution that is peaked around zero and has a parameter that controls for variance can serve as an approximation, in the limit as the variance approaches zero. For example, all three of the above approximations are cumulative distribution functions of common probability distributions: the logistic, Cauchy and normal distributions, respectively.

Often an integral representation of the Heaviside step function is useful:

$${\displaystyle {\begin{aligned}H(x)&=\lim _{\varepsilon \to 0^{+}}-{\frac {1}{2\pi i}}\int _{-\infty }^{\infty }{\frac {1}{\tau +i\varepsilon }}e^{-ix\tau }d\tau \\&=\lim _{\varepsilon \to 0^{+}}{\frac {1}{2\pi i}}\int _{-\infty }^{\infty }{\frac {1}{\tau -i\varepsilon }}e^{ix\tau }d\tau .\end{aligned}}}$$

where the second representation is easy to deduce from the first, given that the step function is real and thus is its own complex conjugate.

Since H is usually used in integration, and the value of a function at a single point does not affect its integral, it rarely matters what particular value is chosen of H(0). Indeed when H is considered as a distribution or an element of L^∞ (see L^p space) it does not even make sense to talk of a value at zero, since such objects are only defined almost everywhere. If using some analytic approximation (as in the examples above) then often whatever happens to be the relevant limit at zero is used.

There exist various reasons for choosing a particular value.

• H(0) = 1/2 is often used since the graph then has rotational symmetry; put another way, H − 1/2 is then an odd function. In this case the following relation with the sign function holds for all x:
  $${\displaystyle H(x)={\tfrac {1}{2}}(1+\operatorname {sgn} x).}$$

  
• H(0) = 1 is used when H needs to be right-continuous. For instance cumulative distribution functions are usually taken to be right continuous, as are functions integrated against in Lebesgue–Stieltjes integration. In this case H is the indicator function of a closed semi-infinite interval:
  $${\displaystyle H(x)=\mathbf {1} _{[0,\infty )}(x).}$$

  
  The corresponding probability distribution is the degenerate distribution.
• H(0) = 0 is used when H needs to be left-continuous. In this case H is an indicator function of an open semi-infinite interval:
  $${\displaystyle H(x)=\mathbf {1} _{(0,\infty )}(x).}$$

  
• In functional-analysis contexts from optimization and game theory, it is often useful to define the Heaviside function as a set-valued function to preserve the continuity of the limiting functions and ensure the existence of certain solutions. In these cases, the Heaviside function returns a whole interval of possible solutions, H(0) = [0,1].

An alternative form of the unit step, defined instead as a function ${\displaystyle H:\mathbb {Z} \rightarrow \mathbb {R} }$ (that is, taking in a discrete variable n), is:

$${\displaystyle H[n]={\begin{cases}0,&n<0,\\1,&n\geq 0,\end{cases}}}$$

or using the half-maximum convention:[3]

$${\displaystyle H[n]={\begin{cases}0,&n<0,\\{\tfrac {1}{2}},&n=0,\\1,&n>0,\end{cases}}}$$

where n is an integer. If n is an integer, then n < 0 must imply that n ≤ −1, while n > 0 must imply that the function attains unity at n = 1. Therefore the "step function" exhibits ramp-like behavior over the domain of [−1, 1], and cannot authentically be a step function, using the half-maximum convention.

Unlike the continuous case, the definition of H[0] is significant.

The discrete-time unit impulse is the first difference of the discrete-time step

$${\displaystyle \delta [n]=H[n]-H[n-1].}$$

This function is the cumulative summation of the Kronecker delta:

$${\displaystyle H[n]=\sum _{k=-\infty }^{n}\delta [k]}$$

where

$${\displaystyle \delta [k]=\delta _{k,0}}$$

is the discrete unit impulse function.

The ramp function is an antiderivative of the Heaviside step function:

$${\displaystyle \int _{-\infty }^{x}H(\xi )\,d\xi =xH(x)=\max\{0,x\}\,.}$$

The distributional derivative of the Heaviside step function is the Dirac delta function:

$${\displaystyle {\frac {dH(x)}{dx}}=\delta (x)\,.}$$

The Fourier transform of the Heaviside step function is a distribution. Using one choice of constants for the definition of the Fourier transform we have

$${\displaystyle {\hat {H}}(s)=\lim _{N\to \infty }\int _{-N}^{N}e^{-2\pi ixs}H(x)\,dx={\frac {1}{2}}\left(\delta (s)-{\frac {i}{\pi }}\operatorname {p.v.} {\frac {1}{s}}\right).}$$

Here p.v.1/s is the distribution that takes a test function φ to the Cauchy principal value of ${\displaystyle \textstyle \int _{-\infty }^{\infty }{\frac {\varphi (s)}{s}}\,ds}$. The limit appearing in the integral is also taken in the sense of (tempered) distributions.

The Laplace transform of the Heaviside step function is a meromorphic function. Using the unilateral Laplace transform we have:

$${\displaystyle {\begin{aligned}{\hat {H}}(s)&=\lim _{N\to \infty }\int _{0}^{N}e^{-sx}H(x)\,dx\\&=\lim _{N\to \infty }\int _{0}^{N}e^{-sx}\,dx\\&={\frac {1}{s}}\end{aligned}}}$$

When the bilateral transform is used, the integral can be split in two parts and the result will be the same.

The Heaviside step function can be represented as a hyperfunction as

$${\displaystyle H(x)=\left(1-{\frac {1}{2\pi i}}\log z,\ -{\frac {1}{2\pi i}}\log z\right).}$$

where log z is the principal value of the complex logarithm of z.

It can also be expressed for x ≠ 0 in terms of the absolute value function as

$${\displaystyle H(x)={\frac {x+|x|}{2x}}\,.}$$

Wikimedia Commons has media related to Heaviside function.

• Digital Library of Mathematical Functions, NIST, .
• Berg, Ernst Julius (1936). "Unit function". Heaviside's Operational Calculus, as applied to Engineering and Physics. McGraw-Hill Education. p. 5.
• Calvert, James B. (2002). "Heaviside, Laplace, and the Inversion Integral". University of Denver.
• Davies, Brian (2002). "Heaviside step function". Integral Transforms and their Applications (3rd ed.). Springer. p. 28.
• Duff, George F. D.; Naylor, D. (1966). "Heaviside unit function". Differential Equations of Applied Mathematics. John Wiley & Sons. p. 42.