text:
Dirac delta function
From Wikipedia, the free encyclopedia
Jump to: navigation, search
"Delta function" redirects here. For other uses, see Delta function (disambiguation).
Schematic representation of the Dirac delta function by a line surmounted by an arrow. The height of the arrow is usually meant to specify the value of any multiplicative constant, which will give the area under the function. The other convention is to write the area next to the arrowhead.
The Dirac delta function as the limit (in the sense of distributions) of the sequence of zero-centered normal distributions δ a ( x ) = 1 | a | π e − ( x / a ) 2 {\displaystyle \delta _{a}(x)={\frac {1}{\left|a\right|{\sqrt {\pi }}}}\mathrm {e} ^{-(x/a)^{2}}} as a → 0 {\displaystyle a\rightarrow 0} .
In science and mathematics, the Dirac delta function, or δ function, is a generalized function, or distribution that was historically introduced by the physicist Paul Dirac for modelling the density of an idealized point mass or point charge, as a function that is equal to zero everywhere except for zero and whose integral over the entire real line is equal to one.[1][2][3] As there is no function that has these properties, the computations that were done by the theoretical physicists appeared to mathematicians as nonsense, until the introduction of distributions by Laurent Schwartz, for formalizing and validating mathematically these computations. This amounts, in particular, to define Dirac delta function as a linear functional, which maps every function to its value at zero.[4][5] The Kronecker delta function, which is usually defined on a discrete domain and takes values 0 and 1, is a discrete analog of Dirac delta function.
In engineering and signal processing, the delta function, also known as the unit impulse symbol,[6] may be regarded through its Laplace transform, as coming from the boundary values of a complex analytic function of a complex variable. The formal rules obeyed by this function are part of the operational calculus, a standard tool kit of physics and engineering. In many applications, the Dirac delta is regarded as a kind of limit (a weak limit) of a sequence of functions having a tall spike at the origin. The approximating functions of the sequence are thus "approximate" or "nascent" delta functions.
Contents
[hide]
1 Motivation and overview
2 History
3 Definitions
3.1 As a measure
3.2 As a distribution
3.3 Generalizations
4 Properties
4.1 Scaling and symmetry
4.2 Algebraic properties
4.3 Translation
4.4 Composition with a function
4.5 Properties in n dimensions
5 Fourier transform
6 Distributional derivatives
6.1 Higher dimensions
7 Representations of the delta function
7.1 Approximations to the identity
7.2 Probabilistic considerations
7.3 Semigroups
7.4 Oscillatory integrals
7.5 Plane wave decomposition
7.6 Fourier kernels
7.7 Hilbert space theory
7.7.1 Spaces of holomorphic functions
7.7.2 Resolutions of the identity
7.8 Infinitesimal delta functions
8 Dirac comb
9 Sokhotski–Plemelj theorem
10 Relationship to the Kronecker delta
11 Applications
11.1 Probability theory
11.2 Quantum mechanics
11.3 Structural mechanics
12 See also
13 Notes
14 References
15 External links
Motivation and overview[edit]
The graph of the delta function is usually thought of as following the whole x-axis and the positive y-axis. The Dirac delta is used to model a tall narrow spike function (an impulse), and other similar abstractions such as a point charge, point mass or electron point. For example, to calculate the dynamics of a billiard ball being struck, one can approximate the force of the impact by a delta function. In doing so, one not only simplifies the equations, but one also is able to calculate the motion of the ball by only considering the total impulse of the collision without a detailed model of all of the elastic energy transfer at subatomic levels (for instance).
To be specific, suppose that a billiard ball is at rest. At time t = 0 {\displaystyle t=0} it is struck by another ball, imparting it with a momentum P, in k g m / s {\displaystyle kg\,m/s} . The exchange of momentum is not actually instantaneous, being mediated by elastic processes at the molecular and subatomic level, but for practical purposes it is convenient to consider that energy transfer as effectively instantaneous. The force therefore is P δ ( t ) {\displaystyle P\delta (t)} . (The units of δ ( t ) {\displaystyle \delta (t)} are s − 1 {\displaystyle s^{-1}} .)
To model this situation more rigorously, suppose that the force instead is uniformly distributed over a small time interval Δ t {\displaystyle \Delta t} . That is,
F Δ t ( t ) = { P / Δ t 0 < t < Δ t 0 otherwise . {\displaystyle F_{\Delta t}(t)={\begin{cases}P/\Delta t&0<t<\Delta t\\0&{\text{otherwise}}\end{cases}}.}
Then the momentum at any time t is found by integration:
p ( t ) = ∫ 0 t F Δ t ( τ ) d τ = { P t > Δ t P t / Δ t 0 < t < Δ t 0 otherwise. {\displaystyle p(t)=\int _{0}^{t}F_{\Delta t}(\tau )\,d\tau ={\begin{cases}P&t>\Delta t\\Pt/\Delta t&0<t<\Delta t\\0&{\text{otherwise.}}\end{cases}}}
Now, the model situation of an instantaneous transfer of momentum requires taking the limit as Δ t → 0 {\displaystyle \Delta t\to 0} , giving
p ( t ) = { P t > 0 0 t ≤ 0. {\displaystyle p(t)={\begin{cases}P&t>0\\0&t\leq 0.\end{cases}}}
Here the functions F Δ t {\displaystyle F_{\Delta t}} are thought of as useful approximations to the idea of instantaneous transfer of momentum.
The delta function allows us to construct an idealized limit of these approximations. Unfortunately, the actual limit of the functions (in the sense of ordinary calculus) lim Δ t → 0 F Δ t {\displaystyle \lim _{\Delta t\to 0}F_{\Delta t}} . is zero everywhere but a single point, where it is infinite. To make proper sense of the delta function, we should instead insist that the property
∫ − ∞ ∞ F Δ t ( t ) d t = P , {\displaystyle \int _{-\infty }^{\infty }F_{\Delta t}(t)\,dt=P,}
which holds for all Δ t > 0 {\displaystyle \Delta t>0} , should continue to hold in the limit. So, in the equation F ( t ) = P δ ( t ) = lim Δ t → 0 F Δ t ( t ) {\displaystyle F(t)=P\delta (t)=\lim _{\Delta t\to 0}F_{\Delta t}(t)} , it is understood that the limit is always taken outside the integral.
In applied mathematics, as we have done here, the delta function is often manipulated as a kind of limit (a weak limit) of a sequence of functions, each member of which has a tall spike at the origin: for example, a sequence of Gaussian distributions centered at the origin with variance tending to zero.
Despite its name, the delta function is not truly a function, at least not a usual one with range in real numbers. For example, the objects f(x) = δ(x) and g(x) = 0 are equal everywhere except at x = 0 yet have integrals that are different. According to Lebesgue integration theory, if f and g are functions such that f = g almost everywhere, then f is integrable if and only if g is integrable and the integrals of f and g are identical. A rigorous approach to regarding the Dirac delta function as a mathematical object in its own right requires measure theory or the theory of distributions.
History[edit]
Joseph Fourier presented what is now called the Fourier integral theorem in his treatise Théorie analytique de la chaleur in the form:[7]
f ( x ) = 1 2 π ∫ − ∞ ∞ d α f ( α ) ∫ − ∞ ∞ d p cos  ( p x − p α ) , {\displaystyle f(x)={\frac {1}{2\pi }}\int _{-\infty }^{\infty }\ \ d\alpha f(\alpha )\ \int _{-\infty }^{\infty }dp\ \cos(px-p\alpha )\ ,}
which is tantamount to the introduction of the δ-function in the form:[8]
δ ( x − α ) = 1 2 π ∫ − ∞ ∞ d p cos  ( p x − p α ) . {\displaystyle \delta (x-\alpha )={\frac {1}{2\pi }}\int _{-\infty }^{\infty }dp\ \cos(px-p\alpha )\ .}
Later, Augustin Cauchy expressed the theorem using exponentials:[9][10]
f ( x ) = 1 2 π ∫ − ∞ ∞ e i p x ( ∫ − ∞ ∞ e − i p α f ( α ) d α ) d p . {\displaystyle f(x)={\frac {1}{2\pi }}\int _{-\infty }^{\infty }\ e^{ipx}\left(\int _{-\infty }^{\infty }e^{-ip\alpha }f(\alpha )\ d\alpha \right)\ dp.}
Cauchy pointed out that in some circumstances the order of integration in this result is significant (contrast Fubini's theorem).[11][12]
As justified using the theory of distributions, the Cauchy equation can be rearranged to resemble Fourier's original formulation and expose the δ-function as:
f ( x ) = 1 2 π ∫ − ∞ ∞ e i p x ( ∫ − ∞ ∞ e − i p α f ( α ) d α ) d p = 1 2 π ∫ − ∞ ∞ ( ∫ − ∞ ∞ e i p x e − i p α d p ) f ( α ) d α = ∫ − ∞ ∞ δ ( x − α ) f ( α ) d α , {\displaystyle {\begin{aligned}f(x)&={\frac {1}{2\pi }}\int _{-\infty }^{\infty }e^{ipx}\left(\int _{-\infty }^{\infty }e^{-ip\alpha }f(\alpha )\ d\alpha \right)\ dp\\&={\frac {1}{2\pi }}\int _{-\infty }^{\infty }\left(\int _{-\infty }^{\infty }e^{ipx}e^{-ip\alpha }\ dp\right)f(\alpha )\ d\alpha =\int _{-\infty }^{\infty }\delta (x-\alpha )f(\alpha )\ d\alpha ,\end{aligned}}}
where the δ-function is expressed as:
δ ( x − α ) = 1 2 π ∫ − ∞ ∞ e i p ( x − α ) d p . {\displaystyle \delta (x-\alpha )={\frac {1}{2\pi }}\int _{-\infty }^{\infty }e^{ip(x-\alpha )}\ dp\ .}
A rigorous interpretation of the exponential form and the various limitations upon the function f necessary for its application extended over several centuries. The problems with a classical interpretation are explained as follows:[13]
The greatest drawback of the classical Fourier transformation is a rather narrow class of functions (originals) for which it can be effectively computed. Namely, it is necessary that these functions decrease sufficiently rapidly to zero (in the neighborhood of infinity) in order to ensure the existence of the Fourier integral. For example, the Fourier transform of such simple functions as polynomials does not exist in the classical sense. The extension of the classical Fourier transformation to distributions considerably enlarged the class of functions that could be transformed and this removed many obstacles.
Further developments included generalization of the Fourier integral, "beginning with Plancherel's pathbreaking L2-theory (1910), continuing with Wiener's and Bochner's works (around 1930) and culminating with the amalgamation into L. Schwartz's theory of distributions (1945) ...",[14] and leading to the formal development of the Dirac delta function.
An infinitesimal formula for an infinitely tall, unit impulse delta function (infinitesimal version of Cauchy distribution) explicitly appears in an 1827 text of Augustin Louis Cauchy.[15] Siméon Denis Poisson considered the issue in connection with the study of wave propagation as did Gustav Kirchhoff somewhat later. Kirchhoff and Hermann von Helmholtz also introduced the unit impulse as a limit of Gaussians, which also corresponded to Lord Kelvin's notion of a point heat source. At the end of the 19th century, Oliver Heaviside used formal Fourier series to manipulate the unit impulse.[16] The Dirac delta function as such was introduced as a "convenient notation" by Paul Dirac in his influential 1930 book The Principles of Quantum Mechanics.[17] He called it the "delta function" since he used it as a continuous analogue of the discrete Kronecker delta.
Definitions[edit]
The Dirac delta can be loosely thought of as a function on the real line which is zero everywhere except at the origin, where it is infinite,
δ ( x ) = { + ∞ , x = 0 0 , x ≠ 0 {\displaystyle \delta (x)={\begin{cases}+\infty ,&x=0\\0,&x\neq 0\end{cases}}}
and which is also constrained to satisfy the identity
∫ − ∞ ∞ δ ( x ) d x = 1. {\displaystyle \int _{-\infty }^{\infty }\delta (x)\,dx=1.} [18]
This is merely a heuristic characterization. The Dirac delta is not a function in the traditional sense as no function defined on the real numbers has these properties.[17] The Dirac delta function can be rigorously defined either as a distribution or as a measure.
As a measure[edit]
One way to rigorously capture the notion of the Dirac delta function is to define a measure, which accepts as an argument a subset A of the real line R, and returns δ(A) = 1 if 0 ∈ A, and δ(A) = 0 otherwise.[19] If the delta function is conceptualized as modeling an idealized point mass at 0, then δ(A) represents the mass contained in the set A. One may then define the integral against δ as the integral of a function against this mass distribution. Formally, the Lebesgue integral provides the necessary analytic device. The Lebesgue integral with respect to the measure δ satisfies
∫ − ∞ ∞ f ( x ) δ { d x } = f ( 0 ) {\displaystyle \int _{-\infty }^{\infty }f(x)\,\delta \{dx\}=f(0)}
for all continuous compactly supported functions f. The measure δ is not absolutely continuous with respect to the Lebesgue measure — in fact, it is a singular measure. Consequently, the delta measure has no Radon–Nikodym derivative — no true function for which the property
∫ − ∞ ∞ f ( x ) δ ( x ) d x = f ( 0 ) {\displaystyle \int _{-\infty }^{\infty }f(x)\delta (x)\,dx=f(0)}
holds.[20] As a result, the latter notation is a convenient abuse of notation, and not a standard (Riemann or Lebesgue) integral.
As a probability measure on R, the delta measure is characterized by its cumulative distribution function, which is the unit step function[21]
H ( x ) = { 1 if x ≥ 0 0 if x < 0. {\displaystyle H(x)={\begin{cases}1&{\text{if }}x\geq 0\\0&{\text{if }}x<0.\end{cases}}}
This means that H(x) is the integral of the cumulative indicator function 1(−∞, x] with respect to the measure δ; to wit,
H ( x ) = ∫ R 1 ( − ∞ , x ] ( t ) δ { d t } = δ ( − ∞ , x ] . {\displaystyle H(x)=\int _{\mathbf {R} }\mathbf {1} _{(-\infty ,x]}(t)\,\delta \{dt\}=\delta (-\infty ,x].}
Thus in particular the integral of the delta function against a continuous function can be properly understood as a Riemann–Stieltjes integral:[22]
∫ − ∞ ∞ f ( x ) δ { d x } = ∫ − ∞ ∞ f ( x ) d H ( x ) . {\displaystyle \int _{-\infty }^{\infty }f(x)\delta \{dx\}=\int _{-\infty }^{\infty }f(x)\,dH(x).}
All higher moments of δ are zero. In particular, characteristic function and moment generating function are both equal to one.
As a distribution[edit]
In the theory of distributions, a generalized function is considered not a function in itself but only in relation to how it affects other functions when "integrated" against them. In keeping with this philosophy, to define the delta function properly, it is enough to say what the "integral" of the delta function is against a sufficiently "good" test function [φ]. Test functions are also known as bump functions. If the delta function is already understood as a measure, then the Lebesgue integral of a test function against that measure supplies the necessary integral.
A typical space of test functions consists of all smooth functions on R with compact support that have as many derivatives as required. As a distribution, the Dirac delta is a linear functional on the space of test functions and is defined by[23]
δ [ φ ] = φ ( 0 ) {\displaystyle \delta [\varphi ]=\varphi (0)}
(1)
for every test function φ.
For δ to be properly a distribution, it must be continuous in a suitable topology on the space of test functions. In general, for a linear functional S on the space of test functions to define a distribution, it is necessary and sufficient that, for every positive integer N there is an integer MN and a constant CN such that for every test function φ, one has the inequality[24]
| S [ φ ] | ≤ C N ∑ k = 0 M N sup x ∈ [ − N , N ] | φ ( k ) ( x ) | . {\displaystyle |S[\varphi ]|\leq C_{N}\sum _{k=0}^{M_{N}}\sup _{x\in [-N,N]}|\varphi ^{(k)}(x)|.}
With the δ distribution, one has such an inequality (with CN = 1) with MN = 0 for all N. Thus δ is a distribution of order zero. It is, furthermore, a distribution with compact support (the support being {0}).
The delta distribution can also be defined in a number of equivalent ways. For instance, it is the distributional derivative of the Heaviside step function. This means that, for every test function φ, one has
δ [ φ ] = − ∫ − ∞ ∞ φ ′ ( x ) H ( x ) d x . {\displaystyle \delta [\varphi ]=-\int _{-\infty }^{\infty }\varphi '(x)H(x)\,dx.}
Intuitively, if integration by parts were permitted, then the latter integral should simplify to
∫ − ∞ ∞ φ ( x ) H ′ ( x ) d x = ∫ − ∞ ∞ φ ( x ) δ ( x ) d x , {\displaystyle \int _{-\infty }^{\infty }\varphi (x)H'(x)\,dx=\int _{-\infty }^{\infty }\varphi (x)\delta (x)\,dx,}
and indeed, a form of integration by parts is permitted for the Stieltjes integral, and in that case one does have
− ∫ − ∞ ∞ φ ′ ( x ) H ( x ) d x = ∫ − ∞ ∞ φ ( x ) d H ( x ) . {\displaystyle -\int _{-\infty }^{\infty }\varphi '(x)H(x)\,dx=\int _{-\infty }^{\infty }\varphi (x)\,dH(x).}
In the context of measure theory, the Dirac measure gives rise to a distribution by integration. Conversely, equation (1) defines a Daniell integral on the space of all compactly supported continuous functions φ which, by the Riesz representation theorem, can be represented as the Lebesgue integral of φ with respect to some Radon measure.
Generally, when the term "Dirac delta function" is used, it is in the sense of distributions rather than measures, the Dirac measure being among several terms for the corresponding notion in measure theory. Some sources may also use the term Dirac delta distribution.
Generalizations[edit]
The delta function can be defined in n-dimensional Euclidean space Rn as the measure such that
∫ R n f ( x ) δ { d x } = f ( 0 ) {\displaystyle \int _{\mathbf {R} ^{n}}f(\mathbf {x} )\delta \{d\mathbf {x} \}=f(\mathbf {0} )}
for every compactly supported continuous function f. As a measure, the n-dimensional delta function is the product measure of the 1-dimensional delta functions in each variable separately. Thus, formally, with x = (x1, x2, ..., xn), one has[6]
δ ( x ) = δ ( x 1 ) δ ( x 2 ) ⋯ δ ( x n ) . {\displaystyle \delta (\mathbf {x} )=\delta (x_{1})\delta (x_{2})\cdots \delta (x_{n}).}
(2)
The delta function can also be defined in the sense of distributions exactly as above in the one-dimensional case.[25] However, despite widespread use in engineering contexts, (2) should be manipulated with care, since the product of distributions can only be defined under quite narrow circumstances.[26]
The notion of a Dirac measure makes sense on any set.[19] Thus if X is a set, x0 ∈ X is a marked point, and Σ is any sigma algebra of subsets of X, then the measure defined on sets A ∈ Σ by
δ x 0 ( A ) = { 1 if x 0 ∈ A 0 if x 0 ∉ A {\displaystyle \delta _{x_{0}}(A)={\begin{cases}1&{\text{if }}x_{0}\in A\\0&{\text{if }}x_{0}\notin A\end{cases}}}
is the delta measure or unit mass concentrated at x0.
Another common generalization of the delta function is to a differentiable manifold where most of its properties as a distribution can also be exploited because of the differentiable structure. The delta function on a manifold M centered at the point x0 ∈ M is defined as the following distribution:
δ x 0 [ φ ] = φ ( x 0 ) {\displaystyle \delta _{x_{0}}[\varphi ]=\varphi (x_{0})}
(3)
for all compactly supported smooth real-valued functions φ on M.[27] A common special case of this construction is when M is an open set in the Euclidean space Rn.
On a locally compact Hausdorff space X, the Dirac delta measure concentrated at a point x is the Radon measure associated with the Daniell integral (3) on compactly supported continuous functions φ. At this level of generality, calculus as such is no longer possible, however a variety of techniques from abstract analysis are available. For instance, the mapping x 0 ↦ δ x 0 {\displaystyle x_{0}\mapsto \delta _{x_{0}}} is a continuous embedding of X into the space of finite Radon measures on X, equipped with its vague topology. Moreover, the convex hull of the image of X under this embedding is dense in the space of probability measures on X.[28]
Properties[edit]
Scaling and symmetry[edit]
The delta function satisfies the following scaling property for a non-zero scalar α:[29]
∫ − ∞ ∞ δ ( α x ) d x = ∫ − ∞ ∞ δ ( u ) d u | α | = 1 | α | {\displaystyle \int _{-\infty }^{\infty }\delta (\alpha x)\,dx=\int _{-\infty }^{\infty }\delta (u)\,{\frac {du}{|\alpha |}}={\frac {1}{|\alpha |}}}
and so
δ ( α x ) = δ ( x ) | α | . {\displaystyle \delta (\alpha x)={\frac {\delta (x)}{|\alpha |}}.}
(4)
In particular, the delta function is an even distribution, in the sense that
δ ( − x ) = δ ( x ) {\displaystyle \delta (-x)=\delta (x)}
which is homogeneous of degree −1.
Algebraic properties[edit]
The distributional product of δ with x is equal to zero:
x δ ( x ) = 0. {\displaystyle x\delta (x)=0.}
Conversely, if xf(x) = xg(x), where f and g are distributions, then
f ( x ) = g ( x ) + c δ ( x ) {\displaystyle f(x)=g(x)+c\delta (x)}
for some constant c.[30]
Translation[edit]
The integral of the time-delayed Dirac delta is given by:
∫ − ∞ ∞ f ( t ) δ ( t − T ) d t = f ( T ) . {\displaystyle \int _{-\infty }^{\infty }f(t)\delta (t-T)\,dt=f(T).}
This is sometimes referred to as the sifting property[31] or the sampling property. The delta function is said to "sift out" the value at t = T.
It follows that the effect of convolving a function f(t) with the time-delayed Dirac delta is to time-delay f(t) by the same amount:
( f ( t ) ∗ δ ( t − T ) ) {\displaystyle (f(t)*\delta (t-T))} = d e f ∫ − ∞ ∞ f ( τ ) δ ( t − T − τ ) d τ {\displaystyle \ {\stackrel {\mathrm {def} }{=}}\ \int _{-\infty }^{\infty }f(\tau )\delta (t-T-\tau )\,d\tau }
= ∫ − ∞ ∞ f ( τ ) δ ( τ − ( t − T ) ) d τ {\displaystyle =\int \limits _{-\infty }^{\infty }f(\tau )\delta (\tau -(t-T))\,d\tau } (using (4): δ ( − x ) = δ ( x ) {\displaystyle \delta (-x)=\delta (x)} )
= f ( t − T ) . {\displaystyle =f(t-T).}
This holds under the precise condition that f be a tempered distribution (see the discussion of the Fourier transform below). As a special case, for instance, we have the identity (understood in the distribution sense)
∫ − ∞ ∞ δ ( ξ − x ) δ ( x − η ) d x = δ ( ξ − η ) . {\displaystyle \int _{-\infty }^{\infty }\delta (\xi -x)\delta (x-\eta )\,dx=\delta (\xi -\eta ).}
Composition with a function[edit]
More generally, the delta distribution may be composed with a smooth function g(x) in such a way that the familiar change of variables formula holds, that
∫ R δ ( g ( x ) ) f ( g ( x ) ) | g ′ ( x ) | d x = ∫ g ( R ) δ ( u ) f ( u ) d u {\displaystyle \int _{\mathbf {R} }\delta {\bigl (}g(x){\bigr )}f{\bigl (}g(x){\bigr )}|g'(x)|\,dx=\int _{g(\mathbf {R} )}\delta (u)f(u)\,du}
provided that g is a continuously differentiable function with g′ nowhere zero.[32] That is, there is a unique way to assign meaning to the distribution δ ∘ g {\displaystyle \delta \circ g} so that this identity holds for all compactly supported test functions f. Therefore, the domain must be broken up to exclude the g′ = 0 point. This distribution satisfies δ(g(x)) = 0 if g is nowhere zero, and otherwise if g has a real root at x0, then
δ ( g ( x ) ) = δ ( x − x 0 ) | g ′ ( x 0 ) | . {\displaystyle \delta (g(x))={\frac {\delta (x-x_{0})}{|g'(x_{0})|}}.}
It is natural therefore to define the composition δ(g(x)) for continuously differentiable functions g by
δ ( g ( x ) ) = ∑ i δ ( x − x i ) | g ′ ( x i ) | {\displaystyle \delta (g(x))=\sum _{i}{\frac {\delta (x-x_{i})}{|g'(x_{i})|}}}
where the sum extends over all roots of g(x), which are assumed to be simple.[32] Thus, for example
δ ( x 2 − α 2 ) = 1 2 | α | [ δ ( x + α ) + δ ( x − α ) ] . {\displaystyle \delta \left(x^{2}-\alpha ^{2}\right)={\frac {1}{2|\alpha |}}{\Big [}\delta \left(x+\alpha \right)+\delta \left(x-\alpha \right){\Big ]}.}
In the integral form the generalized scaling property may be written as
∫ − ∞ ∞ f ( x ) δ ( g ( x ) ) d x = ∑ i f ( x i ) | g ′ ( x i ) | . {\displaystyle \int _{-\infty }^{\infty }f(x)\,\delta (g(x))\,dx=\sum _{i}{\frac {f(x_{i})}{|g'(x_{i})|}}.}
Properties in n dimensions[edit]
The delta distribution in an n-dimensional space satisfies the following scaling property instead:
δ ( α x ) = | α | − n δ ( x ) {\displaystyle \delta (\alpha \mathbf {x} )=|\alpha |^{-n}\delta (\mathbf {x} )}
so that δ is a homogeneous distribution of degree −n. Under any reflection or rotation ρ, the delta function is invariant:
δ ( ρ x ) = δ ( x ) . {\displaystyle \delta (\rho \mathbf {x} )=\delta (\mathbf {x} ).}
As in the one-variable case, it is possible to define the composition of δ with a bi-Lipschitz function[33] g: Rn → Rn uniquely so that the identity
∫ R n δ ( g ( x ) ) f ( g ( x ) ) | det g ′ ( x ) | d x = ∫ g ( R n ) δ ( u ) f ( u ) d u {\displaystyle \int _{\mathbf {R} ^{n}}\delta (g(\mathbf {x} ))\,f(g(\mathbf {x} ))\,|\det g'(\mathbf {x} )|\,d\mathbf {x} =\int _{g(\mathbf {R} ^{n})}\delta (\mathbf {u} )f(\mathbf {u} )\,d\mathbf {u} }
for all compactly supported functions f.
Using the coarea formula from geometric measure theory, one can also define the composition of the delta function with a submersion from one Euclidean space to another one of different dimension; the result is a type of current. In the special case of a continuously differentiable function g: Rn → R such that the gradient of g is nowhere zero, the following identity holds[34]
∫ R n f ( x ) δ ( g ( x ) ) d x = ∫ g − 1 ( 0 ) f ( x ) | ∇ g | d σ ( x ) {\displaystyle \int _{\mathbf {R} ^{n}}f(\mathbf {x} )\,\delta (g(\mathbf {x} ))\,d\mathbf {x} =\int _{g^{-1}(0)}{\frac {f(\mathbf {x} )}{|\mathbf {\nabla } g|}}\,d\sigma (\mathbf {x} )}
where the integral on the right is over g−1(0), the (n − 1)-dimensional surface defined by g(x) = 0 with respect to the Minkowski content measure. This is known as a simple layer integral.
More generally, if S is a smooth hypersurface of Rn, then we can associate to S the distribution that integrates any compactly supported smooth function g over S:
δ S [ g ] = ∫ S g ( s ) d σ ( s ) {\displaystyle \delta _{S}[g]=\int _{S}g(\mathbf {s} )\,d\sigma (\mathbf {s} )}
where σ is the hypersurface measure associated to S. This generalization is associated with the potential theory of simple layer potentials on S. If D is a domain in Rn with smooth boundary S, then δS is equal to the normal derivative of the indicator function of D in the distribution sense:
− ∫ R n g ( x ) ∂ 1 D ( x ) ∂ n d x = ∫ S g ( s ) d σ ( s ) , {\displaystyle -\int _{\mathbf {R} ^{n}}g(\mathbf {x} )\,{\frac {\partial 1_{D}(\mathbf {x} )}{\partial n}}\;d\mathbf {x} =\int _{S}\,g(\mathbf {s} )\;d\sigma (\mathbf {s} ),}
where n is the outward normal.[35][36] For a proof, see e.g. the article on the surface delta function.
Fourier transform[edit]
The delta function is a tempered distribution, and therefore it has a well-defined Fourier transform. Formally, one finds[37]
δ ^ ( ξ ) = ∫ − ∞ ∞ e − 2 π i x ξ δ ( x ) d x = 1. {\displaystyle {\hat {\delta }}(\xi )=\int _{-\infty }^{\infty }e^{-2\pi ix\xi }\delta (x)\,dx=1.}
Properly speaking, the Fourier transform of a distribution is defined by imposing self-adjointness of the Fourier transform under the duality pairing ⟨ ⋅ , ⋅ ⟩ {\displaystyle \langle \cdot ,\cdot \rangle } of tempered distributions with Schwartz functions. Thus δ ^ {\displaystyle {\hat {\delta }}} is defined as the unique tempered distribution satisfying
⟨ δ ^ , φ ⟩ = ⟨ δ , φ ^ ⟩ {\displaystyle \langle {\hat {\delta }},\varphi \rangle =\langle \delta ,{\hat {\varphi }}\rangle }
for all Schwartz functions φ. And indeed it follows from this that δ ^ = 1. {\displaystyle {\hat {\delta }}=1.}
As a result of this identity, the convolution of the delta function with any other tempered distribution S is simply S:
S ∗ δ = S . {\displaystyle S*\delta =S.}
That is to say that δ is an identity element for the convolution on tempered distributions, and in fact the space of compactly supported distributions under convolution is an associative algebra with identity the delta function. This property is fundamental in signal processing, as convolution with a tempered distribution is a linear time-invariant system, and applying the linear time-invariant system measures its impulse response. The impulse response can be computed to any desired degree of accuracy by choosing a suitable approximation for δ, and once it is known, it characterizes the system completely. See LTI system theory § Impulse response and convolution.
The inverse Fourier transform of the tempered distribution f(ξ) = 1 is the delta function. Formally, this is expressed
∫ − ∞ ∞ 1 ⋅ e 2 π i x ξ d ξ = δ ( x ) {\displaystyle \int _{-\infty }^{\infty }1\cdot e^{2\pi ix\xi }\,d\xi =\delta (x)}
and more rigorously, it follows since
⟨ 1 , f ∨ ⟩ = f ( 0 ) = ⟨ δ , f ⟩ {\displaystyle \langle 1,f^{\vee }\rangle =f(0)=\langle \delta ,f\rangle }
for all Schwartz functions f.
In these terms, the delta function provides a suggestive statement of the orthogonality property of the Fourier kernel on R. Formally, one has
∫ − ∞ ∞ e i 2 π ξ 1 t [ e i 2 π ξ 2 t ] ∗ d t = ∫ − ∞ ∞ e − i 2 π ( ξ 2 − ξ 1 ) t d t = δ ( ξ 2 − ξ 1 ) . {\displaystyle \int _{-\infty }^{\infty }e^{i2\pi \xi _{1}t}\left[e^{i2\pi \xi _{2}t}\right]^{*}\,dt=\int _{-\infty }^{\infty }e^{-i2\pi (\xi _{2}-\xi _{1})t}\,dt=\delta (\xi _{2}-\xi _{1}).}
This is, of course, shorthand for the assertion that the Fourier transform of the tempered distribution
f ( t ) = e i 2 π ξ 1 t {\displaystyle f(t)=e^{i2\pi \xi _{1}t}}
is
f ^ ( ξ 2 ) = δ ( ξ 1 − ξ 2 ) {\displaystyle {\hat {f}}(\xi _{2})=\delta (\xi _{1}-\xi _{2})}
which again follows by imposing self-adjointness of the Fourier transform.
By analytic continuation of the Fourier transform, the Laplace transform of the delta function is found to be[38]
∫ 0 ∞ δ ( t − a ) e − s t d t = e − s a . {\displaystyle \int _{0}^{\infty }\delta (t-a)e^{-st}\,dt=e^{-sa}.}
Distributional derivatives[edit]
The distributional derivative of the Dirac delta distribution is the distribution δ′ defined on compactly supported smooth test functions φ by[39]
δ ′ [ φ ] = − δ [ φ ′ ] = − φ ′ ( 0 ) . {\displaystyle \delta '[\varphi ]=-\delta [\varphi ']=-\varphi '(0).}
The first equality here is a kind of integration by parts, for if δ were a true function then
∫ − ∞ ∞ δ ′ ( x ) φ ( x ) d x = − ∫ − ∞ ∞ δ ( x ) φ ′ ( x ) d x . {\displaystyle \int _{-\infty }^{\infty }\delta '(x)\varphi (x)\,dx=-\int _{-\infty }^{\infty }\delta (x)\varphi '(x)\,dx.}
The k-th derivative of δ is defined similarly as the distribution given on test functions by
δ ( k ) [ φ ] = ( − 1 ) k φ ( k ) ( 0 ) . {\displaystyle \delta ^{(k)}[\varphi ]=(-1)^{k}\varphi ^{(k)}(0).}
In particular, δ is an infinitely differentiable distribution.
The first derivative of the delta function is the distributional limit of the difference quotients:[40]
δ ′ ( x ) = lim h → 0 δ ( x + h ) − δ ( x ) h . {\displaystyle \delta '(x)=\lim _{h\to 0}{\frac {\delta (x+h)-\delta (x)}{h}}.}
More properly, one has
δ ′ = lim h → 0 1 h ( τ h δ − δ ) {\displaystyle \delta '=\lim _{h\to 0}{\frac {1}{h}}(\tau _{h}\delta -\delta )}
where τh is the translation operator, defined on functions by τhφ(x) = φ(x + h), and on a distribution S by
( τ h S ) [ φ ] = S [ τ − h φ ] . {\displaystyle (\tau _{h}S)[\varphi ]=S[\tau _{-h}\varphi ].}
In the theory of electromagnetism, the first derivative of the delta function represents a point magnetic dipole situated at the origin. Accordingly, it is referred to as a dipole or the doublet function.[41]
The derivative of the delta function satisfies a number of basic properties, including:
d d x δ ( − x ) = d d x δ ( x ) δ ′ ( − x ) = − δ ′ ( x ) x δ ′ ( x ) = − δ ( x ) . {\displaystyle {\begin{aligned}&{\frac {d}{dx}}\delta (-x)={\frac {d}{dx}}\delta (x)\\[8pt]&\delta '(-x)=-\delta '(x)\\[8pt]&x\delta '(x)=-\delta (x).\end{aligned}}} [42]
Furthermore, the convolution of δ′ with a compactly supported smooth function f is
δ ′ ∗ f = δ ∗ f ′ = f ′ , {\displaystyle \delta '*f=\delta *f'=f',}
which follows from the properties of the distributional derivative of a convolution.
Higher dimensions[edit]
More generally, on an open set U in the n-dimensional Euclidean space Rn, the Dirac delta distribution centered at a point a ∈ U is defined by[43]
δ a [ φ ] = φ ( a ) {\displaystyle \delta _{a}[\varphi ]=\varphi (a)}
for all φ ∈ S(U), the space of all smooth compactly supported functions on U. If α = (α1, ..., αn) is any multi-index and ∂α denotes the associated mixed partial derivative operator, then the αth derivative ∂αδa of δa is given by[43]
⟨ ∂ α δ a , φ ⟩ = ( − 1 ) | α | ⟨ δ a , ∂ α φ ⟩ = ( − 1 ) | α | ∂ α φ ( x ) | x = a for all φ ∈ S ( U ) . {\displaystyle \left\langle \partial ^{\alpha }\delta _{a},\varphi \right\rangle =(-1)^{|\alpha |}\left\langle \delta _{a},\partial ^{\alpha }\varphi \right\rangle =\left.(-1)^{|\alpha |}\partial ^{\alpha }\varphi (x)\right|_{x=a}{\mbox{ for all }}\varphi \in S(U).}
That is, the αth derivative of δa is the distribution whose value on any test function φ is the αth derivative of φ at a (with the appropriate positive or negative sign).
The first partial derivatives of the delta function are thought of as double layers along the coordinate planes. More generally, the normal derivative of a simple layer supported on a surface is a double layer supported on that surface, and represents a laminar magnetic monopole. Higher derivatives of the delta function are known in physics as multipoles.
Higher derivatives enter into mathematics naturally as the building blocks for the complete structure of distributions with point support. If S is any distribution on U supported on the set {a} consisting of a single point, then there is an integer m and coefficients cα such that[44]
S = ∑ | α | ≤ m c α ∂ α δ a . {\displaystyle S=\sum _{|\alpha |\leq m}c_{\alpha }\partial ^{\alpha }\delta _{a}.}
Representations of the delta function[edit]
The delta function can be viewed as the limit of a sequence of functions
δ ( x ) = lim ε → 0 + η ε ( x ) , {\displaystyle \delta (x)=\lim _{\varepsilon \to 0^{+}}\eta _{\varepsilon }(x),}
where ηε(x) is sometimes called a nascent delta function. This limit is meant in a weak sense: either that
lim ε → 0 + ∫ − ∞ ∞ η ε ( x ) f ( x ) d x = f ( 0 ) {\displaystyle \lim _{\varepsilon \to 0^{+}}\int _{-\infty }^{\infty }\eta _{\varepsilon }(x)f(x)\,dx=f(0)}
(5)
for all continuous functions f having compact support, or that this limit holds for all smooth functions f with compact support. The difference between these two slightly different modes of weak convergence is often subtle: the former is convergence in the vague topology of measures, and the latter is convergence in the sense of distributions.
Approximations to the identity[edit]
Typically a nascent delta function ηε can be constructed in the following manner. Let η be an absolutely integrable function on R of total integral 1, and define
η ε ( x ) = ε − 1 η ( x ε ) . {\displaystyle \eta _{\varepsilon }(x)=\varepsilon ^{-1}\eta \left({\frac {x}{\varepsilon }}\right).}
In n dimensions, one uses instead the scaling
η ε ( x ) = ε − n η ( x ε ) . {\displaystyle \eta _{\varepsilon }(x)=\varepsilon ^{-n}\eta \left({\frac {x}{\varepsilon }}\right).}
Then a simple change of variables shows that ηε also has integral 1. One may show that (5) holds for all continuous compactly supported functions f,[45] and so ηε converges weakly to δ in the sense of measures.
The ηε constructed in this way are known as an approximation to the identity.[46] This terminology is because the space L1(R) of absolutely integrable functions is closed under the operation of convolution of functions: f ∗ g ∈ L1(R) whenever f and g are in L1(R). However, there is no identity in L1(R) for the convolution product: no element h such that f ∗ h = f for all f. Nevertheless, the sequence ηε does approximate such an identity in the sense that
f ∗ η ε → f a s ε → 0. {\displaystyle f*\eta _{\varepsilon }\to f\quad {\rm {{as\ }\varepsilon \to 0.}}}
This limit holds in the sense of mean convergence (convergence in L1). Further conditions on the ηε, for instance that it be a mollifier associated to a compactly supported function,[47] are needed to ensure pointwise convergence almost everywhere.
If the initial η = η1 is itself smooth and compactly supported then the sequence is called a mollifier. The standard mollifier is obtained by choosing η to be a suitably normalized bump function, for instance
η ( x ) = { e − 1 1 − | x | 2 if | x | < 1 0 if | x | ≥ 1. {\displaystyle \eta (x)={\begin{cases}e^{-{\frac {1}{1-|x|^{2}}}}&{\text{ if }}|x|<1\\0&{\text{ if }}|x|\geq 1.\end{cases}}}
In some situations such as numerical analysis, a piecewise linear approximation to the identity is desirable. This can be obtained by taking η1 to be a hat function. With this choice of η1, one has
η ε ( x ) = ε − 1 max ( 1 − | x ε | , 0 ) {\displaystyle \eta _{\varepsilon }(x)=\varepsilon ^{-1}\max \left(1-\left|{\frac {x}{\varepsilon }}\right|,0\right)}
which are all continuous and compactly supported, although not smooth and so not a mollifier.
Probabilistic considerations[edit]
In the context of probability theory, it is natural to impose the additional condition that the initial η1 in an approximation to the identity should be positive, as such a function then represents a probability distribution. Convolution with a probability distribution is sometimes favorable because it does not result in overshoot or undershoot, as the output is a convex combination of the input values, and thus falls between the maximum and minimum of the input function. Taking η1 to be any probability distribution at all, and letting ηε(x) = η1(x/ε)/ε as above will give rise to an approximation to the identity. In general this converges more rapidly to a delta function if, in addition, η has mean 0 and has small higher moments. For instance, if η1 is the uniform distribution on [−1/2, 1/2], also known as the rectangular function, then:[48]
η ε ( x ) = 1 ε rect ( x ε ) = { 1 ε , − ε 2 < x < ε 2 0 , otherwise . {\displaystyle \eta _{\varepsilon }(x)={\frac {1}{\varepsilon }}\ {\textrm {rect}}\left({\frac {x}{\varepsilon }}\right)={\begin{cases}{\frac {1}{\varepsilon }},&-{\frac {\varepsilon }{2}}<x<{\frac {\varepsilon }{2}}\\0,&{\text{otherwise}}.\end{cases}}}
Another example is with the Wigner semicircle distribution
η ε ( x ) = { 2 π ε 2 ε 2 − x 2 , − ε < x < ε 0 , otherwise {\displaystyle \eta _{\varepsilon }(x)={\begin{cases}{\frac {2}{\pi \varepsilon ^{2}}}{\sqrt {\varepsilon ^{2}-x^{2}}},&-\varepsilon <x<\varepsilon \\0,&{\text{otherwise}}\end{cases}}}
This is continuous and compactly supported, but not a mollifier because it is not smooth.
Semigroups[edit]
Nascent delta functions often arise as convolution semigroups. This amounts to the further constraint that the convolution of ηε with ηδ must satisfy
η ε ∗ η δ = η ε + δ {\displaystyle \eta _{\varepsilon }*\eta _{\delta }=\eta _{\varepsilon +\delta }}
for all ε, δ > 0. Convolution semigroups in L1 that form a nascent delta function are always an approximation to the identity in the above sense, however the semigroup condition is quite a strong restriction.
In practice, semigroups approximating the delta function arise as fundamental solutions or Green's functions to physically motivated elliptic or parabolic partial differential equations. In the context of applied mathematics, semigroups arise as the output of a linear time-invariant system. Abstractly, if A is a linear operator acting on functions of x, then a convolution semigroup arises by solving the initial value problem
{ ∂ ∂ t η ( t , x ) = A η ( t , x ) , t > 0 lim t → 0 + η ( t , x ) = δ ( x ) {\displaystyle {\begin{cases}{\frac {\partial }{\partial t}}\eta (t,x)=A\eta (t,x),\quad t>0\\\displaystyle \lim _{t\to 0^{+}}\eta (t,x)=\delta (x)\end{cases}}}
in which the limit is as usual understood in the weak sense. Setting ηε(x) = η(ε, x) gives the associated nascent delta function.
Some examples of physically important convolution semigroups arising from such a fundamental solution include the following.
The heat kernel
The heat kernel, defined by
η ε ( x ) = 1 2 π ε e − x 2 2 ε {\displaystyle \eta _{\varepsilon }(x)={\frac {1}{\sqrt {2\pi \varepsilon }}}\mathrm {e} ^{-{\frac {x^{2}}{2\varepsilon }}}}
represents the temperature in an infinite wire at time t > 0, if a unit of heat energy is stored at the origin of the wire at time t = 0. This semigroup evolves according to the one-dimensional heat equation:
∂ u ∂ t = 1 2 ∂ 2 u ∂ x 2 . {\displaystyle {\frac {\partial u}{\partial t}}={\frac {1}{2}}{\frac {\partial ^{2}u}{\partial x^{2}}}.}
In probability theory, ηε(x) is a normal distribution of variance ε and mean 0. It represents the probability density at time t = ε of the position of a particle starting at the origin following a standard Brownian motion. In this context, the semigroup condition is then an expression of the Markov property of Brownian motion.
In higher-dimensional Euclidean space Rn, the heat kernel is
η ε = 1 ( 2 π ε ) n / 2 e − x ⋅ x 2 ε , {\displaystyle \eta _{\varepsilon }={\frac {1}{(2\pi \varepsilon )^{n/2}}}\mathrm {e} ^{-{\frac {x\cdot x}{2\varepsilon }}},}
and has the same physical interpretation, mutatis mutandis. It also represents a nascent delta function in the sense that ηε → δ in the distribution sense as ε → 0.
The Poisson kernel
The Poisson kernel
η ε ( x ) = 1 π ε ε 2 + x 2 = ∫ − ∞ ∞ e 2 π i ξ x − | ε ξ | d ξ {\displaystyle \eta _{\varepsilon }(x)={\frac {1}{\pi }}{\frac {\varepsilon }{\varepsilon ^{2}+x^{2}}}=\int _{-\infty }^{\infty }\mathrm {e} ^{2\pi \mathrm {i} \xi x-|\varepsilon \xi |}\;d\xi }
is the fundamental solution of the Laplace equation in the upper half-plane.[49] It represents the electrostatic potential in a semi-infinite plate whose potential along the edge is held at fixed at the delta function. The Poisson kernel is also closely related to the Cauchy distribution. This semigroup evolves according to the equation
∂ u ∂ t = − ( − ∂ 2 ∂ x 2 ) 1 2 u ( t , x ) {\displaystyle {\frac {\partial u}{\partial t}}=-\left(-{\frac {\partial ^{2}}{\partial x^{2}}}\right)^{\frac {1}{2}}u(t,x)}
where the operator is rigorously defined as the Fourier multiplier
F [ ( − ∂ 2 ∂ x 2 ) 1 2 f ] ( ξ ) = | 2 π ξ | F f ( ξ ) . {\displaystyle {\mathcal {F}}\left[\left(-{\frac {\partial ^{2}}{\partial x^{2}}}\right)^{\frac {1}{2}}f\right](\xi )=|2\pi \xi |{\mathcal {F}}f(\xi ).}
Oscillatory integrals[edit]
In areas of physics such as wave propagation and wave mechanics, the equations involved are hyperbolic and so may have more singular solutions. As a result, the nascent delta functions that arise as fundamental solutions of the associated Cauchy problems are generally oscillatory integrals. An example, which comes from a solution of the Euler–Tricomi equation of transonic gas dynamics,[50] is the rescaled Airy function
ε − 1 3 Ai  ( x ε − 1 3 ) . {\displaystyle \varepsilon ^{-{\frac {1}{3}}}\operatorname {Ai} \left(x\varepsilon ^{-{\frac {1}{3}}}\right).}
Although using the Fourier transform, it is easy to see that this generates a semigroup in some sense—it is not absolutely integrable and so cannot define a semigroup in the above strong sense. Many nascent delta functions constructed as oscillatory integrals only converge in the sense of distributions (an example is the Dirichlet kernel below), rather than in the sense of measures.
Another example is the Cauchy problem for the wave equation in R1+1:[51]
c − 2 ∂ 2 u ∂ t 2 − Δ u = 0 u = 0 , ∂ u ∂ t = δ for t = 0. {\displaystyle {\begin{aligned}c^{-2}{\frac {\partial ^{2}u}{\partial t^{2}}}-\Delta u&=0\\u=0,\quad {\frac {\partial u}{\partial t}}=\delta &\qquad {\text{for }}t=0.\end{aligned}}}
The solution u represents the displacement from equilibrium of an infinite elastic string, with an initial disturbance at the origin.
Other approximations to the identity of this kind include the sinc function (used widely in electronics and telecommunications)
η ε ( x ) = 1 π x sin  ( x ε ) = 1 2 π ∫ − 1 ε 1 ε cos  ( k x ) d k {\displaystyle \eta _{\varepsilon }(x)={\frac {1}{\pi x}}\sin \left({\frac {x}{\varepsilon }}\right)={\frac {1}{2\pi }}\int _{-{\frac {1}{\varepsilon }}}^{\frac {1}{\varepsilon }}\cos(kx)\;dk}
and the Bessel function
η ε ( x ) = 1 ε J 1 ε ( x + 1 ε ) . {\displaystyle \eta _{\varepsilon }(x)={\frac {1}{\varepsilon }}J_{\frac {1}{\varepsilon }}\left({\frac {x+1}{\varepsilon }}\right).}
Plane wave decomposition[edit]
One approach to the study of a linear partial differential equation
L [ u ] = f , {\displaystyle L[u]=f,}
where L is a differential operator on Rn, is to seek first a fundamental solution, which is a solution of the equation
L [ u ] = δ . {\displaystyle L[u]=\delta .}
When L is particularly simple, this problem can often be resolved using the Fourier transform directly (as in the case of the Poisson kernel and heat kernel already mentioned). For more complicated operators, it is sometimes easier first to consider an equation of the form
L [ u ] = h {\displaystyle L[u]=h}
where h is a plane wave function, meaning that it has the form
h = h ( x ⋅ ξ ) {\displaystyle h=h(x\cdot \xi )}
for some vector ξ. Such an equation can be resolved (if the coefficients of L are analytic functions) by the Cauchy–Kovalevskaya theorem or (if the coefficients of L are constant) by quadrature. So, if the delta function can be decomposed into plane waves, then one can in principle solve linear partial differential equations.
Such a decomposition of the delta function into plane waves was part of a general technique first introduced essentially by Johann Radon, and then developed in this form by Fritz John (1955).[52] Choose k so that n + k is an even integer, and for a real number s, put
g ( s ) = Re  [ − s k log  ( − i s ) k ! ( 2 π i ) n ] = { | s | k 4 k ! ( 2 π i ) n − 1 n odd − | s | k log  | s | k ! ( 2 π i ) n n even. {\displaystyle g(s)=\operatorname {Re} \left[{\frac {-s^{k}\log(-is)}{k!(2\pi i)^{n}}}\right]={\begin{cases}{\frac {|s|^{k}}{4k!(2\pi i)^{n-1}}}&n{\text{ odd}}\\&\\-{\frac {|s|^{k}\log |s|}{k!(2\pi i)^{n}}}&n{\text{ even.}}\end{cases}}}
Then δ is obtained by applying a power of the Laplacian to the integral with respect to the unit sphere measure dω of g(x · ξ) for ξ in the unit sphere Sn−1:
δ ( x ) = Δ x n + k 2 ∫ S n − 1 g ( x ⋅ ξ ) d ω ξ . {\displaystyle \delta (x)=\Delta _{x}^{\frac {n+k}{2}}\int _{S^{n-1}}g(x\cdot \xi )\,d\omega _{\xi }.}
The Laplacian here is interpreted as a weak derivative, so that this equation is taken to mean that, for any test function φ,
φ ( x ) = ∫ R n φ ( y ) d y Δ x n + k 2 ∫ S n − 1 g ( ( x − y ) ⋅ ξ ) d ω ξ . {\displaystyle \varphi (x)=\int _{\mathbf {R} ^{n}}\varphi (y)\,dy\,\Delta _{x}^{\frac {n+k}{2}}\int _{S^{n-1}}g((x-y)\cdot \xi )\,d\omega _{\xi }.}
The result follows from the formula for the Newtonian potential (the fundamental solution of Poisson's equation). This is essentially a form of the inversion formula for the Radon transform, because it recovers the value of φ(x) from its integrals over hyperplanes. For instance, if n is odd and k = 1, then the integral on the right hand side is
c n Δ x n + 1 2 ∫ ∫ S n − 1 φ ( y ) | ( y − x ) ⋅ ξ | d ω ξ d y = c n Δ x n + 1 2 ∫ S n − 1 d ω ξ ∫ − ∞ ∞ | p | R φ ( ξ , p + x ⋅ ξ ) d p {\displaystyle c_{n}\Delta _{x}^{\frac {n+1}{2}}\int \int _{S^{n-1}}\varphi (y)|(y-x)\cdot \xi |\,d\omega _{\xi }\,dy=c_{n}\Delta _{x}^{\frac {n+1}{2}}\int _{S^{n-1}}\,d\omega _{\xi }\int _{-\infty }^{\infty }|p|R\varphi (\xi ,p+x\cdot \xi )\,dp}
where Rφ(ξ, p) is the Radon transform of φ:
R φ ( ξ , p ) = ∫ x ⋅ ξ = p f ( x ) d n − 1 x . {\displaystyle R\varphi (\xi ,p)=\int _{x\cdot \xi =p}f(x)\,d^{n-1}x.}
An alternative equivalent expression of the plane wave decomposition, from Gel'fand & Shilov (1966–1968, I, §3.10), is
δ ( x ) = ( n − 1 ) ! ( 2 π i ) n ∫ S n − 1 ( x ⋅ ξ ) − n d ω ξ {\displaystyle \delta (x)={\frac {(n-1)!}{(2\pi i)^{n}}}\int _{S^{n-1}}(x\cdot \xi )^{-n}\,d\omega _{\xi }}
for n even, and
δ ( x ) = 1 2 ( 2 π i ) n − 1 ∫ S n − 1 δ ( n − 1 ) ( x ⋅ ξ ) d ω ξ {\displaystyle \delta (x)={\frac {1}{2(2\pi i)^{n-1}}}\int _{S^{n-1}}\delta ^{(n-1)}(x\cdot \xi )\,d\omega _{\xi }}
for n odd.
Fourier kernels[edit]
See also: Convergence of Fourier series
In the study of Fourier series, a major question consists of determining whether and in what sense the Fourier series associated with a periodic function converges to the function. The nth partial sum of the Fourier series of a function f of period 2π is defined by convolution (on the interval [−π,π]) with the Dirichlet kernel:
D N ( x ) = ∑ n = − N N e i n x = sin  ( ( N + 1 2 ) x ) sin  ( x / 2 ) . {\displaystyle D_{N}(x)=\sum _{n=-N}^{N}e^{inx}={\frac {\sin \left((N+{\tfrac {1}{2}})x\right)}{\sin(x/2)}}.}
Thus,
s N ( f ) ( x ) = D N ∗ f ( x ) = ∑ n = − N N a n e i n x {\displaystyle s_{N}(f)(x)=D_{N}*f(x)=\sum _{n=-N}^{N}a_{n}e^{inx}}
where
a n = 1 2 π ∫ − π π f ( y ) e − i n y d y . {\displaystyle a_{n}={\frac {1}{2\pi }}\int _{-\pi }^{\pi }f(y)e^{-iny}\,dy.}
A fundamental result of elementary Fourier series states that the Dirichlet kernel tends to the a multiple of the delta function as N → ∞. This is interpreted in the distribution sense, that
s N ( f ) ( 0 ) = ∫ R D N ( x ) f ( x ) d x → 2 π f ( 0 ) {\displaystyle s_{N}(f)(0)=\int _{\mathbf {R} }D_{N}(x)f(x)\,dx\to 2\pi f(0)}
for every compactly supported smooth function f. Thus, formally one has
δ ( x ) = 1 2 π ∑ n = − ∞ ∞ e i n x {\displaystyle \delta (x)={\frac {1}{2\pi }}\sum _{n=-\infty }^{\infty }e^{inx}}
on the interval [−π,π].
In spite of this, the result does not hold for all compactly supported continuous functions: that is DN does not converge weakly in the sense of measures. The lack of convergence of the Fourier series has led to the introduction of a variety of summability methods in order to produce convergence. The method of Cesàro summation leads to the Fejér kernel[53]
F N ( x ) = 1 N ∑ n = 0 N − 1 D n ( x ) = 1 N ( sin  N x 2 sin  x 2 ) 2 . {\displaystyle F_{N}(x)={\frac {1}{N}}\sum _{n=0}^{N-1}D_{n}(x)={\frac {1}{N}}\left({\frac {\sin {\frac {Nx}{2}}}{\sin {\frac {x}{2}}}}\right)^{2}.}
The Fejér kernels tend to the delta function in a stronger sense that[54]
∫ R F N ( x ) f ( x ) d x → 2 π f ( 0 ) {\displaystyle \int _{\mathbf {R} }F_{N}(x)f(x)\,dx\to 2\pi f(0)}
for every compactly supported continuous function f. The implication is that the Fourier series of any continuous function is Cesàro summable to the value of the function at every point.
Hilbert space theory[edit]
The Dirac delta distribution is a densely defined unbounded linear functional on the Hilbert space L2 of square integrable functions. Indeed, smooth compactly support functions are dense in L2, and the action of the delta distribution on such functions is well-defined. In many applications, it is possible to identify subspaces of L2 and to give a stronger topology on which the delta function defines a bounded linear functional.
Sobolev spaces
The Sobolev embedding theorem for Sobolev spaces on the real line R implies that any square-integrable function f such that
∥ f ∥ H 1 2 = ∫ − ∞ ∞ | f ^ ( ξ ) | 2 ( 1 + | ξ | 2 ) d ξ < ∞ {\displaystyle \|f\|_{H^{1}}^{2}=\int _{-\infty }^{\infty }|{\hat {f}}(\xi )|^{2}(1+|\xi |^{2})\,d\xi <\infty }
is automatically continuous, and satisfies in particular
δ [ f ] = | f ( 0 ) | < C ∥ f ∥ H 1 . {\displaystyle \delta [f]=|f(0)|<C\|f\|_{H^{1}}.}
Thus δ is a bounded linear functional on the Sobolev space H1. Equivalently δ is an element of the continuous dual space H−1 of H1. More generally, in n dimensions, one has δ ∈ H−s(Rn) provided s > n / 2.
Spaces of holomorphic functions[edit]
In complex analysis, the delta function enters via Cauchy's integral formula which asserts that if D is a domain in the complex plane with smooth boundary, then
f ( z ) = 1 2 π i ∮ ∂ D f ( ζ ) d ζ ζ − z , z ∈ D {\displaystyle f(z)={\frac {1}{2\pi i}}\oint _{\partial D}{\frac {f(\zeta )\,d\zeta }{\zeta -z}},\quad z\in D}
for all holomorphic functions f in D that are continuous on the closure of D. As a result, the delta function δz is represented on this class of holomorphic functions by the Cauchy integral:
δ z [ f ] = f ( z ) = 1 2 π i ∮ ∂ D f ( ζ ) d ζ ζ − z . {\displaystyle \delta _{z}[f]=f(z)={\frac {1}{2\pi i}}\oint _{\partial D}{\frac {f(\zeta )\,d\zeta }{\zeta -z}}.}
More generally, let H2(∂D) be the Hardy space consisting of the closure in L2(∂D) of all holomorphic functions in D continuous up to the boundary of D. Then functions in H2(∂D) uniquely extend to holomorphic functions in D, and the Cauchy integral formula continues to hold. In particular for z ∈ D, the delta function δz is a continuous linear functional on H2(∂D). This is a special case of the situation in several complex variables in which, for smooth domains D, the Szegő kernel plays the role of the Cauchy integral.
Resolutions of the identity[edit]
Given a complete orthonormal basis set of functions {φn} in a separable Hilbert space, for example, the normalized eigenvectors of a compact self-adjoint operator, any vector f can be expressed as:
f = ∑ n = 1 ∞ α n φ n . {\displaystyle f=\sum _{n=1}^{\infty }\alpha _{n}\varphi _{n}.}
The coefficients {αn} are found as:
α n = ⟨ φ n , f ⟩ , {\displaystyle \alpha _{n}=\langle \varphi _{n},f\rangle ,}
which may be represented by the notation:
α n = φ n † f , {\displaystyle \alpha _{n}=\varphi _{n}^{\dagger }f,}
a form of the bra–ket notation of Dirac.[55] Adopting this notation, the expansion of f takes the dyadic form:[56]
f = ∑ n = 1 ∞ φ n ( φ n † f ) . {\displaystyle f=\sum _{n=1}^{\infty }\varphi _{n}\left(\varphi _{n}^{\dagger }f\right).}
Letting I denote the identity operator on the Hilbert space, the expression
I = ∑ n = 1 ∞ φ n φ n † , {\displaystyle I=\sum _{n=1}^{\infty }\varphi _{n}\varphi _{n}^{\dagger },}
is called a resolution of the identity. When the Hilbert space is the space L2(D) of square-integrable functions on a domain D, the quantity:
φ n φ n † , {\displaystyle \varphi _{n}\varphi _{n}^{\dagger },}
is an integral operator, and the expression for f can be rewritten as:
f ( x ) = ∑ n = 1 ∞ ∫ D ( φ n ( x ) φ n ∗ ( ξ ) ) f ( ξ ) d ξ . {\displaystyle f(x)=\sum _{n=1}^{\infty }\int _{D}\,\left(\varphi _{n}(x)\varphi _{n}^{*}(\xi )\right)f(\xi )\,d\xi .}
The right-hand side converges to f in the L2 sense. It need not hold in a pointwise sense, even when f is a continuous function. Nevertheless, it is common to abuse notation and write
f ( x ) = ∫ δ ( x − ξ ) f ( ξ ) d ξ , {\displaystyle f(x)=\int \,\delta (x-\xi )f(\xi )\,d\xi ,}
resulting in the representation of the delta function:[57]
δ ( x − ξ ) = ∑ n = 1 ∞ φ n ( x ) φ n ∗ ( ξ ) . {\displaystyle \delta (x-\xi )=\sum _{n=1}^{\infty }\varphi _{n}(x)\varphi _{n}^{*}(\xi ).}
With a suitable rigged Hilbert space (Φ, L2(D), Φ*) where Φ ⊂ L2(D) contains all compactly supported smooth functions, this summation may converge in Φ*, depending on the properties of the basis φn. In most cases of practical interest, the orthonormal basis comes from an integral or differential operator, in which case the series converges in the distribution sense.[58]
Infinitesimal delta functions[edit]
Cauchy used an infinitesimal α to write down a unit impulse, infinitely tall and narrow Dirac-type delta function δα satisfying ∫ F ( x ) δ α ( x ) = F ( 0 ) {\displaystyle \int F(x)\delta _{\alpha }(x)=F(0)} in a number of articles in 1827.[59] Cauchy defined an infinitesimal in Cours d'Analyse (1827) in terms of a sequence tending to zero. Namely, such a null sequence becomes an infinitesimal in Cauchy's and Lazare Carnot's terminology.
Non-standard analysis allows one to rigorously treat infinitesimals. The article by Yamashita (2007) contains a bibliography on modern Dirac delta functions in the context of an infinitesimal-enriched continuum provided by the hyperreals. Here the Dirac delta can be given by an actual function, having the property that for every real function F one has ∫ F ( x ) δ α ( x ) = F ( 0 ) {\displaystyle \int F(x)\delta _{\alpha }(x)=F(0)} as anticipated by Fourier and Cauchy.
Dirac comb[edit]
Main article: Dirac comb
A Dirac comb is an infinite series of Dirac delta functions spaced at intervals of T
A so-called uniform "pulse train" of Dirac delta measures, which is known as a Dirac comb, or as the Shah distribution, creates a sampling function, often used in digital signal processing (DSP) and discrete time signal analysis. The Dirac comb is given as the infinite sum, whose limit is understood in the distribution sense,
Δ ( x ) = ∑ n = − ∞ ∞ δ ( x − n ) , {\displaystyle \Delta (x)=\sum _{n=-\infty }^{\infty }\delta (x-n),}
which is a sequence of point masses at each of the integers.
Up to an overall normalizing constant, the Dirac comb is equal to its own Fourier transform. This is significant because if f is any Schwartz function, then the periodization of f is given by the convolution
( f ∗ Δ ) ( x ) = ∑ n = − ∞ ∞ f ( x − n ) . {\displaystyle (f*\Delta )(x)=\sum _{n=-\infty }^{\infty }f(x-n).}
In particular,
( f ∗ Δ ) ∧ = f ^ Δ ^ = f ^ Δ {\displaystyle (f*\Delta )^{\wedge }={\hat {f}}{\widehat {\Delta }}={\hat {f}}\Delta }
is precisely the Poisson summation formula.[60]
Sokhotski–Plemelj theorem[edit]
The Sokhotski–Plemelj theorem, important in quantum mechanics, relates the delta function to the distribution p.v.1/x, the Cauchy principal value of the function 1/x, defined by
⟨ p . v .  1 x , φ ⟩ = lim ε → 0 + ∫ | x | > ε φ ( x ) x d x . {\displaystyle \left\langle \operatorname {p.v.} {\frac {1}{x}},\varphi \right\rangle =\lim _{\varepsilon \to 0^{+}}\int _{|x|>\varepsilon }{\frac {\varphi (x)}{x}}\,dx.}
Sokhotsky's formula states that[61]
lim ε → 0 + 1 x ± i ε = p . v .  1 x ∓ i π δ ( x ) , {\displaystyle \lim _{\varepsilon \to 0^{+}}{\frac {1}{x\pm i\varepsilon }}=\operatorname {p.v.} {\frac {1}{x}}\mp i\pi \delta (x),}
Here the limit is understood in the distribution sense, that for all compactly supported smooth functions f,
lim ε → 0 + ∫ − ∞ ∞ f ( x ) x ± i ε d x = ∓ i π f ( 0 ) + lim ε → 0 + ∫ | x | > ε f ( x ) x d x . {\displaystyle \lim _{\varepsilon \to 0^{+}}\int _{-\infty }^{\infty }{\frac {f(x)}{x\pm i\varepsilon }}\,dx=\mp i\pi f(0)+\lim _{\varepsilon \to 0^{+}}\int _{|x|>\varepsilon }{\frac {f(x)}{x}}\,dx.}
Relationship to the Kronecker delta[edit]
The Kronecker delta δij is the quantity defined by
δ i j = { 1 i = j 0 i ≠ j {\displaystyle \delta _{ij}={\begin{cases}1&i=j\\0&i\not =j\end{cases}}}
for all integers i, j. This function then satisfies the following analog of the sifting property: if ( a i ) i ∈ Z {\displaystyle (a_{i})_{i\in \mathbf {Z} }} is any doubly infinite sequence, then
∑ i = − ∞ ∞ a i δ i k = a k . {\displaystyle \sum _{i=-\infty }^{\infty }a_{i}\delta _{ik}=a_{k}.}
Similarly, for any real or complex valued continuous function f on R, the Dirac delta satisfies the sifting property
∫ − ∞ ∞ f ( x ) δ ( x − x 0 ) d x = f ( x 0 ) . {\displaystyle \int _{-\infty }^{\infty }f(x)\delta (x-x_{0})\,dx=f(x_{0}).}
This exhibits the Kronecker delta function as a discrete analog of the Dirac delta function.[62]
Applications[edit]
Probability theory[edit]
In probability theory and statistics, the Dirac delta function is often used to represent a discrete distribution, or a partially discrete, partially continuous distribution, using a probability density function (which is normally used to represent fully continuous distributions). For example, the probability density function f(x) of a discrete distribution consisting of points x = {x1, ..., xn}, with corresponding probabilities p1, ..., pn, can be written as
f ( x ) = ∑ i = 1 n p i δ ( x − x i ) . {\displaystyle f(x)=\sum _{i=1}^{n}p_{i}\delta (x-x_{i}).}
As another example, consider a distribution which 6/10 of the time returns a standard normal distribution, and 4/10 of the time returns exactly the value 3.5 (i.e. a partly continuous, partly discrete mixture distribution). The density function of this distribution can be written as
f ( x ) = 0.6 1 2 π e − x 2 2 + 0.4 δ ( x − 3.5 ) . {\displaystyle f(x)=0.6\,{\frac {1}{\sqrt {2\pi }}}e^{-{\frac {x^{2}}{2}}}+0.4\,\delta (x-3.5).}
The delta function is also used in a completely different way to represent the local time of a diffusion process (like Brownian motion). The local time of a stochastic process B(t) is given by
ℓ ( x , t ) = ∫ 0 t δ ( x − B ( s ) ) d s {\displaystyle \ell (x,t)=\int _{0}^{t}\delta (x-B(s))\,ds}
and represents the amount of time that the process spends at the point x in the range of the process. More precisely, in one dimension this integral can be written
ℓ ( x , t ) = lim ε → 0 + 1 2 ε ∫ 0 t 1 [ x − ε , x + ε ] ( B ( s ) ) d s {\displaystyle \ell (x,t)=\lim _{\varepsilon \to 0^{+}}{\frac {1}{2\varepsilon }}\int _{0}^{t}\mathbf {1} _{[x-\varepsilon ,x+\varepsilon ]}(B(s))\,ds}
where 1[x−ε, x+ε] is the indicator function of the interval [x−ε, x+ε].
Quantum mechanics[edit]
The delta function is expedient in quantum mechanics. The wave function of a particle gives the probability amplitude of finding a particle within a given region of space. Wave functions are assumed to be elements of the Hilbert space L2 of square-integrable functions, and the total probability of finding a particle within a given interval is the integral of the magnitude of the wave function squared over the interval. A set {φn} of wave functions is orthonormal if they are normalized by
⟨ φ n | φ m ⟩ = δ n m {\displaystyle \langle \varphi _{n}|\varphi _{m}\rangle =\delta _{nm}}
where δ here refers to the Kronecker delta. A set of orthonormal wave functions is complete in the space of square-integrable functions if any wave function ψ can be expressed as a combination of the φn:
ψ = ∑ c n φ n , {\displaystyle \psi =\sum c_{n}\varphi _{n},}
with c n = ⟨ φ n | ψ ⟩ {\displaystyle c_{n}=\langle \varphi _{n}|\psi \rangle } . Complete orthonormal systems of wave functions appear naturally as the eigenfunctions of the Hamiltonian (of a bound system) in quantum mechanics that measures the energy levels, which are called the eigenvalues. The set of eigenvalues, in this case, is known as the spectrum of the Hamiltonian. In bra–ket notation, as above, this equality implies the resolution of the identity:
I = ∑ | φ n ⟩ ⟨ φ n | . {\displaystyle I=\sum |\varphi _{n}\rangle \langle \varphi _{n}|.}
Here the eigenvalues are assumed to be discrete, but the set of eigenvalues of an observable may be continuous rather than discrete. An example is the position observable, Qψ(x) = xψ(x). The spectrum of the position (in one dimension) is the entire real line, and is called a continuous spectrum. However, unlike the Hamiltonian, the position operator lacks proper eigenfunctions. The conventional way to overcome this shortcoming is to widen the class of available functions by allowing distributions as well: that is, to replace the Hilbert space of quantum mechanics by an appropriate rigged Hilbert space.[63] In this context, the position operator has a complete set of eigen-distributions, labeled by the points y of the real line, given by
φ y ( x ) = δ ( x − y ) . {\displaystyle \varphi _{y}(x)=\delta (x-y).}
The eigenfunctions of position are denoted by φ y = | y ⟩ {\displaystyle \varphi _{y}=|y\rangle } in Dirac notation, and are known as position eigenstates.
Similar considerations apply to the eigenstates of the momentum operator, or indeed any other self-adjoint unbounded operator P on the Hilbert space, provided the spectrum of P is continuous and there are no degenerate eigenvalues. In that case, there is a set Ω of real numbers (the spectrum), and a collection φy of distributions indexed by the elements of Ω, such that
P φ y = y φ y . {\displaystyle P\varphi _{y}=y\varphi _{y}.}
That is, φy are the eigenvectors of P. If the eigenvectors are normalized so that
⟨ φ y , φ y ′ ⟩ = δ ( y − y ′ ) {\displaystyle \langle \varphi _{y},\varphi _{y'}\rangle =\delta (y-y')}
in the distribution sense, then for any test function ψ,
ψ ( x ) = ∫ Ω c ( y ) φ y ( x ) d y {\displaystyle \psi (x)=\int _{\Omega }c(y)\varphi _{y}(x)\,dy}
where
c ( y ) = ⟨ ψ , φ y ⟩ . {\displaystyle c(y)=\langle \psi ,\varphi _{y}\rangle .}
That is, as in the discrete case, there is a resolution of the identity
I = ∫ Ω | φ y ⟩ ⟨ φ y | d y {\displaystyle I=\int _{\Omega }|\varphi _{y}\rangle \,\langle \varphi _{y}|\,dy}
where the operator-valued integral is again understood in the weak sense. If the spectrum of P has both continuous and discrete parts, then the resolution of the identity involves a summation over the discrete spectrum and an integral over the continuous spectrum.
The delta function also has many more specialized applications in quantum mechanics, such as the delta potential models for a single and double potential well.
Structural mechanics[edit]
The delta function can be used in structural mechanics to describe transient loads or point loads acting on structures. The governing equation of a simple mass–spring system excited by a sudden force impulse I at time t = 0 can be written
m d 2 ξ d t 2 + k ξ = I δ ( t ) , {\displaystyle m{\frac {\mathrm {d} ^{2}\xi }{\mathrm {d} t^{2}}}+k\xi =I\delta (t),}
where m is the mass, ξ the deflection and k the spring constant.
As another example, the equation governing the static deflection of a slender beam is, according to Euler–Bernoulli theory,
E I d 4 w d x 4 = q ( x ) , {\displaystyle EI{\frac {\mathrm {d} ^{4}w}{\mathrm {d} x^{4}}}=q(x),}
where EI is the bending stiffness of the beam, w the deflection, x the spatial coordinate and q(x) the load distribution. If a beam is loaded by a point force F at x = x0, the load distribution is written
q ( x ) = F δ ( x − x 0 ) . {\displaystyle q(x)=F\delta (x-x_{0}).}
As integration of the delta function results in the Heaviside step function, it follows that the static deflection of a slender beam subject to multiple point loads is described by a set of piecewise polynomials.
Also a point moment acting on a beam can be described by delta functions. Consider two opposing point forces F at a distance d apart. They then produce a moment M = Fd acting on the beam. Now, let the distance d approach the limit zero, while M is kept constant. The load distribution, assuming a clockwise moment acting at x = 0, is written
q ( x ) = lim d → 0 ( F δ ( x ) − F δ ( x − d ) ) = lim d → 0 ( M d δ ( x ) − M d δ ( x − d ) ) = M lim d → 0 δ ( x ) − δ ( x − d ) d = M δ ′ ( x ) . {\displaystyle {\begin{aligned}q(x)&=\lim _{d\to 0}{\Big (}F\delta (x)-F\delta (x-d){\Big )}\\&=\lim _{d\to 0}\left({\frac {M}{d}}\delta (x)-{\frac {M}{d}}\delta (x-d)\right)\\&=M\lim _{d\to 0}{\frac {\delta (x)-\delta (x-d)}{d}}\\&=M\delta '(x).\end{aligned}}}
Point moments can thus be represented by the derivative of the delta function. Integration of the beam equation again results in piecewise polynomial deflection.
See also[edit]
Atom (measure theory)
Delta potential
Dirac measure
Fundamental solution
Green's function
Laplacian of the indicator
Notes[edit]
Jump up ^ Arfken & Weber 2000, p. 84
Jump up ^ Dirac 1958, §15 The δ function, p. 58
Jump up ^ Gel'fand & Shilov 1968, Volume I, §1.1
Jump up ^ Gel'fand & Shilov 1968, Volume I, §1.3
Jump up ^ Schwartz 1950, p. 3
^ Jump up to: a b Bracewell 1986, Chapter 5
Jump up ^ JB Fourier (1822). The Analytical Theory of Heat (English translation by Alexander Freeman, 1878 ed.). The University Press. p. 408. , cf p 449 and pp 546–551. The original French text can be found here.
Jump up ^ Hikosaburo Komatsu (2002). "Fourier's hyperfunctions and Heaviside's pseudodifferential operators". In Takahiro Kawai; Keiko Fujita. Microlocal Analysis and Complex Fourier Analysis. World Scientific. p. 200. ISBN 981-238-161-9.
Jump up ^ Tyn Myint-U.; Lokenath Debnath (2007). Linear Partial Differential Equations for Scientists And Engineers (4th ed.). Springer. p. 4. ISBN 0-8176-4393-1.
Jump up ^ Lokenath Debnath; Dambaru Bhatta (2007). Integral Transforms And Their Applications (2nd ed.). CRC Press. p. 2. ISBN 1-58488-575-0.
Jump up ^ Ivor Grattan-Guinness (2009). Convolutions in French Mathematics, 1800–1840: From the Calculus and Mechanics to Mathematical Analysis and Mathematical Physics, Volume 2. Birkhäuser. p. 653. ISBN 3-7643-2238-1.
Jump up ^ See, for example, Des intégrales doubles qui se présentent sous une forme indéterminèe
Jump up ^ Dragiša Mitrović; Darko Žubrinić (1998). Fundamentals of Applied Functional Analysis: Distributions, Sobolev Spaces. CRC Press. p. 62. ISBN 0-582-24694-6.
Jump up ^ Manfred Kracht; Erwin Kreyszig (1989). "On singular integral operators and generalizations". In Themistocles M. Rassias. Topics in Mathematical Analysis: A Volume Dedicated to the Memory of A.L. Cauchy. World Scientific. p. 553. ISBN 9971-5-0666-1.
Jump up ^ Laugwitz 1989, p. 230
Jump up ^ A more complete historical account can be found in van der Pol & Bremmer 1987, §V.4.
^ Jump up to: a b Dirac 1958, §15
Jump up ^ Gel'fand & Shilov 1968, Volume I, §1.1, p. 1
^ Jump up to: a b Rudin 1966, §1.20
Jump up ^ Hewitt & Stromberg 1963, §19.61
Jump up ^ Driggers 2003, p. 2321. See also Bracewell 1986, Chapter 5 for a different interpretation. Other conventions for the assigning the value of the Heaviside function at zero exist, and some of these are not consistent with what follows.
Jump up ^ Hewitt & Stromberg 1965, §9.19
Jump up ^ Strichartz 1994, §2.2
Jump up ^ Hörmander 1983, Theorem 2.1.5
Jump up ^ Hörmander 1983, §3.1
Jump up ^ Strichartz 1994, §2.3; Hörmander 1983, §8.2
Jump up ^ Dieudonné 1972, §17.3.3
Jump up ^ Federer 1969, §2.5.19
Jump up ^ Strichartz 1994, Problem 2.6.2
Jump up ^ Vladimirov 1971, Chapter 2, Example 3(d)
Jump up ^ Weisstein, Eric W. "Sifting Property". MathWorld.
^ Jump up to: a b Gel'fand & Shilov 1966–1968, Vol. 1, §II.2.5
Jump up ^ Further refinement is possible, namely to submersions, although these require a more involved change of variables formula.
Jump up ^ Hörmander 1983, §6.1
Jump up ^ Lange 2012, pp.29–30
Jump up ^ Gelfand & Shilov, p. 212
Jump up ^ In some conventions for the Fourier transform.
Jump up ^ Bracewell 1986
Jump up ^ Gel'fand & Shilov 1966, p. 26
Jump up ^ Gel'fand & Shilov 1966, §2.1
Jump up ^ Weisstein, Eric W. "Doublet Function". MathWorld.
Jump up ^ The property follows by applying a test function and integration by parts.
^ Jump up to: a b Hörmander 1983, p. 56
Jump up ^ Hörmander 1983, p. 56; Rudin 1991, Theorem 6.25
Jump up ^ Stein & Weiss, Theorem 1.18
Jump up ^ Rudin 1991, §II.6.31
Jump up ^ More generally, one only needs η = η1 to have an integrable radially symmetric decreasing rearrangement.
Jump up ^ Saichev & Woyczyński 1997, §1.1 The "delta function" as viewed by a physicist and an engineer, p. 3
Jump up ^ Stein & Weiss 1971, §I.1
Jump up ^ Vallée & Soares 2004, §7.2
Jump up ^ Hörmander 1983, §7.8
Jump up ^ See also Courant & Hilbert 1962, §14.
Jump up ^ Lang 1997, p. 312
Jump up ^ In the terminology of Lang (1997), the Fejér kernel is a Dirac sequence, whereas the Dirichlet kernel is not.
Jump up ^ The development of this section in bra–ket notation is found in (Levin 2002, Coordinate-space wave functions and completeness, pp.=109ff)
Jump up ^ Davis & Thomson 2000, Perfect operators, p.344
Jump up ^ Davis & Thomson 2000, Equation 8.9.11, p. 344
Jump up ^ de la Madrid, Bohm & Gadella 2002
Jump up ^ See Laugwitz (1989).
Jump up ^ Córdoba 1988; Hörmander 1983, §7.2
Jump up ^ Vladimirov 1971, §5.7
Jump up ^ Hartmann 1997, pp. 154–155
Jump up ^ Isham 1995, §6.2
References[edit]
Aratyn, Henrik; Rasinariu, Constantin (2006), A short course in mathematical methods with Maple, World Scientific, ISBN 981-256-461-6 .
Arfken, G. B.; Weber, H. J. (2000), Mathematical Methods for Physicists (5th ed.), Boston, Massachusetts: Academic Press, ISBN 978-0-12-059825-0 .
Bracewell, R. (1986), The Fourier Transform and Its Applications (2nd ed.), McGraw-Hill .
Córdoba, A., "La formule sommatoire de Poisson", C.R. Acad. Sci. Paris, Series I, 306: 373–376 .
Courant, Richard; Hilbert, David (1962), Methods of Mathematical Physics, Volume II, Wiley-Interscience .
Davis, Howard Ted; Thomson, Kendall T (2000), Linear algebra and linear operators in engineering with applications in Mathematica, Academic Press, ISBN 0-12-206349-X
Dieudonné, Jean (1976), Treatise on analysis. Vol. II, New York: Academic Press [Harcourt Brace Jovanovich Publishers], ISBN 978-0-12-215502-4, MR 0530406 .
Dieudonné, Jean (1972), Treatise on analysis. Vol. III, Boston, Massachusetts: Academic Press, MR 0350769
Dirac, Paul (1958), The Principles of Quantum Mechanics (4th ed.), Oxford at the Clarendon Press, ISBN 978-0-19-852011-5 .
Driggers, Ronald G. (2003), Encyclopedia of Optical Engineering, CRC Press, ISBN 978-0-8247-0940-2 .
Duistermaat, Hans; Kolk (2010), Distributions: Theory and applications, Springer .
Federer, Herbert (1969), Geometric measure theory, Die Grundlehren der mathematischen Wissenschaften, 153, New York: Springer-Verlag, pp. xiv+676, ISBN 978-3-540-60656-7, MR 0257325 .
Gannon, Terry (2008), "Vertex operator algebras", Princeton Companion to Mathematics, Princeton University Press .
Gel'fand, I. M.; Shilov, G. E. (1966–1968), Generalized functions, 1–5, Academic Press .
Hartman, William M. (1997), Signals, sound, and sensation, Springer, ISBN 978-1-56396-283-7 .
Hewitt, E; Stromberg, K (1963), Real and abstract analysis, Springer-Verlag .
Hörmander, L. (1983), The analysis of linear partial differential operators I, Grundl. Math. Wissenschaft., 256, Springer, ISBN 3-540-12104-8, MR 0717035 .
Isham, C. J. (1995), Lectures on quantum theory: mathematical and structural foundations, Imperial College Press, ISBN 978-81-7764-190-5 .
John, Fritz (1955), Plane waves and spherical means applied to partial differential equations, Interscience Publishers, New York-London, MR 0075429 .
Lang, Serge (1997), Undergraduate analysis, Undergraduate Texts in Mathematics (2nd ed.), Berlin, New York: Springer-Verlag, ISBN 978-0-387-94841-6, MR 1476913 .
Lange, Rutger-Jan (2012), "Potential theory, path integrals and the Laplacian of the indicator", Journal of High Energy Physics, Springer, 2012 (11): 29–30, Bibcode:2012JHEP...11..032L, arXiv:1302.0864  , doi:10.1007/JHEP11(2012)032 .
Laugwitz, D. (1989), "Definite values of infinite sums: aspects of the foundations of infinitesimal analysis around 1820", Arch. Hist. Exact Sci., 39 (3): 195–245, doi:10.1007/BF00329867 .
Levin, Frank S. (2002), "Coordinate-space wave functions and completeness", An introduction to quantum theory, Cambridge University Press, pp. 109ff, ISBN 0-521-59841-9
Li, Y. T.; Wong, R. (2008), "Integral and series representations of the Dirac delta function", Commun. Pure Appl. Anal., 7 (2): 229–247, MR 2373214, doi:10.3934/cpaa.2008.7.229 .
de la Madrid, R.; Bohm, A.; Gadella, M. (2002), "Rigged Hilbert Space Treatment of Continuous Spectrum", Fortschr. Phys., 50 (2): 185–216, Bibcode:2002ForPh..50..185D, arXiv:quant-ph/0109154  , doi:10.1002/1521-3978(200203)50:2<185::AID-PROP185>3.0.CO;2-S .
McMahon, D. (2005-11-22), "An Introduction to State Space", Quantum Mechanics Demystified, A Self-Teaching Guide, Demystified Series, New York: McGraw-Hill, p. 108, ISBN 0-07-145546-9, doi:10.1036/0071455469, retrieved 2008-03-17 .
van der Pol, Balth.; Bremmer, H. (1987), Operational calculus (3rd ed.), New York: Chelsea Publishing Co., ISBN 978-0-8284-0327-6, MR 904873 .
Rudin, W. (1991), Functional Analysis (2nd ed.), McGraw-Hill, ISBN 0-07-054236-8 .
Soares, Manuel; Vallée, Olivier (2004), Airy functions and applications to physics, London: Imperial College Press .
Saichev, A I; Woyczyński, Wojbor Andrzej (1997), "Chapter1: Basic definitions and operations", Distributions in the Physical and Engineering Sciences: Distributional and fractal calculus, integral transforms, and wavelets, Birkhäuser, ISBN 0-8176-3924-1
Schwartz, L. (1950), Théorie des distributions, 1, Hermann .
Schwartz, L. (1951), Théorie des distributions, 2, Hermann .
Stein, Elias; Weiss, Guido (1971), Introduction to Fourier Analysis on Euclidean Spaces, Princeton University Press, ISBN 0-691-08078-X .
Strichartz, R. (1994), A Guide to Distribution Theory and Fourier Transforms, CRC Press, ISBN 0-8493-8273-4 .
Vladimirov, V. S. (1971), Equations of mathematical physics, Marcel Dekker, ISBN 0-8247-1713-9 .
Weisstein, Eric W. "Delta Function". MathWorld.
Yamashita, H. (2006), "Pointwise analysis of scalar fields: A nonstandard approach", Journal of Mathematical Physics, 47 (9): 092301, Bibcode:2006JMP....47i2301Y, doi:10.1063/1.2339017
Yamashita, H. (2007), "Comment on "Pointwise analysis of scalar fields: A nonstandard approach" [J. Math. Phys. 47, 092301 (2006)]", Journal of Mathematical Physics, 48 (8): 084101, Bibcode:2007JMP....48h4101Y, doi:10.1063/1.2771422
External links[edit]
Hazewinkel, Michiel, ed. (2001) [1994], "Delta-function", Encyclopedia of Mathematics, Springer Science+Business Media B.V. / Kluwer Academic Publishers, ISBN 978-1-55608-010-4
KhanAcademy.org video lesson
The Dirac Delta function, a tutorial on the Dirac delta function.
Video Lectures – Lecture 23, a lecture by Arthur Mattuck.
Dirac Delta Function on PlanetMath
The Dirac delta measure is a hyperfunction
We show the existence of a unique solution and analyze a finite element approximation when the source term is a Dirac delta measure
Non-Lebesgue measures on R. Lebesgue-Stieltjes measure, Dirac delta measure.
[show]
v
t
e
Probability distributions
List
Discrete univariate
with finite support
Benford
Bernoulli
beta-binomial
binomial
categorical
hypergeometric
Poisson binomial
Rademacher
discrete uniform
Zipf
Zipf–Mandelbrot
Discrete univariate
with infinite support
beta negative binomial
Borel
Conway–Maxwell–Poisson
discrete phase-type
Delaporte
extended negative binomial
Gauss–Kuzmin
geometric
logarithmic
negative binomial
parabolic fractal
Poisson
Skellam
Yule–Simon
zeta
Continuous univariate
supported on a bounded interval
arcsine
ARGUS
Balding–Nichols
Bates
beta
beta rectangular
Irwin–Hall
Kumaraswamy
logit-normal
noncentral beta
raised cosine
reciprocal
triangular
U-quadratic
uniform
Wigner semicircle
Continuous univariate
supported on a semi-infinite interval
Benini
Benktander 1st kind
Benktander 2nd kind
beta prime
Burr
chi-squared
chi
Dagum
Davis
exponential-logarithmic
Erlang
exponential
F
folded normal
Flory–Schulz
Fréchet
gamma
gamma/Gompertz
generalized inverse Gaussian
Gompertz
half-logistic
half-normal
Hotelling's T-squared
hyper-Erlang
hyperexponential
hypoexponential
inverse chi-squared
scaled inverse chi-squared
inverse Gaussian
inverse gamma
Kolmogorov
Lévy
log-Cauchy
log-Laplace
log-logistic
log-normal
Lomax
matrix-exponential
Maxwell–Boltzmann
Maxwell–Jüttner
Mittag-Leffler
Nakagami
noncentral chi-squared
Pareto
phase-type
poly-Weibull
Rayleigh
relativistic Breit–Wigner
Rice
shifted Gompertz
truncated normal
type-2 Gumbel
Weibull
Discrete Weibull
Wilks's lambda
Continuous univariate
supported on the whole real line
Cauchy
exponential power
Fisher's z
Gaussian q
generalized normal
generalized hyperbolic
geometric stable
Gumbel
Holtsmark
hyperbolic secant
Johnson's SU
Landau
Laplace
asymmetric Laplace
logistic
noncentral t
normal (Gaussian)
normal-inverse Gaussian
skew normal
slash
stable
Student's t
type-1 Gumbel
Tracy–Widom
variance-gamma
Voigt
Continuous univariate
with support whose type varies
generalized extreme value
generalized Pareto
Marchenko–Pastur
q-exponential
q-Gaussian
q-Weibull
shifted log-logistic
Tukey lambda
Mixed continuous-discrete univariate
rectified Gaussian
Multivariate (joint)
Discrete
Ewens
multinomial
Dirichlet-multinomial
negative multinomial
Continuous
Dirichlet
generalized Dirichlet
multivariate Laplace
multivariate normal
multivariate stable
multivariate t
normal-inverse-gamma
normal-gamma
Matrix-valued
inverse matrix gamma
inverse-Wishart
matrix normal
matrix t
matrix gamma
normal-inverse-Wishart
normal-Wishart
Wishart
Directional
Univariate (circular) directional
Circular uniform
univariate von Mises
wrapped normal
wrapped Cauchy
wrapped exponential
wrapped asymmetric Laplace
wrapped Lévy
Bivariate (spherical)
Kent
Bivariate (toroidal)
bivariate von Mises
Multivariate
von Mises–Fisher
Bingham
Degenerate and singular
Degenerate
Dirac delta function
Singular
Cantor
Families
Circular
compound Poisson
elliptical
exponential
natural exponential
location–scale
maximum entropy
mixture
Pearson
Tweedie
wrapped
[show]
v
t
e
Infinitesimals
History
Adequality
Leibniz's notation
Integral symbol
Criticism of non-standard analysis
The Analyst
The Method of Mechanical Theorems
Cavalieri's principle
Method of indivisibles
Related branches
Non-standard analysis
Non-standard calculus
Internal set theory
Synthetic differential geometry
Constructive non-standard analysis
Infinitesimal strain theory (physics)
Formalizations
Differentials
Hyperreal numbers
Dual numbers
Surreal numbers
Individual concepts
Standard part function
Transfer principle
Hyperinteger
Increment theorem
Monad
Internal set
Levi-Civita field
Hyperfinite set
Law of Continuity
Overspill
Microcontinuity
Transcendental law of homogeneity
Scientists
Gottfried Wilhelm Leibniz
Abraham Robinson
Pierre de Fermat
Augustin-Louis Cauchy
Leonhard Euler
Textbooks
Analyse des Infiniment Petits
Elementary Calculus
Cours d'Analyse
Retrieved from "https://en.wikipedia.org/w/index.php?title=Dirac_delta_function&oldid=798167826"
Categories:
Mathematics of infinitesimals
Fourier analysis
Generalized functions
Measure theory
Digital signal processing
Paul Dirac
Hidden categories:
Good articles
Navigation menu
Personal tools
Not logged in
Talk
Contributions
Create account
Log in
Namespaces
Article
Talk
Variants
Views
Read
Edit
View history
More
Search
Navigation
Main page
Contents
Featured content
Current events
Random article
Donate to Wikipedia
Wikipedia store
Interaction
Help
About Wikipedia
Community portal
Recent changes
Contact page
Tools
What links here
Related changes
Upload file
Special pages
Permanent link
Page information
Wikidata item
Cite this page
Print/export
Create a book
Download as PDF
Printable version
In other projects
Wikimedia Commons
Languages
العربية
Беларуская
Català
Čeština
Dansk
Deutsch
Ελληνικά
Español
Esperanto
فارسی
Français
한국어
हिन्दी
Íslenska
Italiano
עברית
ქართული
Latviešu
Magyar
Nederlands
日本語
Norsk
Polski
Português
Română
Русский
සිංහල
Simple English
Slovenščina
Српски / srpski
Suomi
Svenska
Татарча/tatarça
ไทย
Türkçe
Українська
Tiếng Việt
中文
Edit links
This page was last edited on 31 August 2017, at 10:45.
Text is available under the Creative Commons Attribution-ShareAlike License; additional terms may apply. By using this site, you agree to the Terms of Use and Privacy Policy. Wikipedia® is a registered trademark of the Wikimedia Foundation, Inc., a non-profit organization.
Privacy policy
About Wikipedia
Disclaimers
Contact Wikipedia
Developers
Cookie statement
Mobile view

values: mw-page-base noprint

values: mw-head-base noprint

attrs: id class style

values: This is a good article. Click here for more information.

attrs: class id

values: firstHeading firstHeading en

values: siteSub noprint

values: jump-to-nav mw-jump

attrs: href class title

values:
/wiki/Delta_function_(disambiguation)
mw-disambig
Delta function (disambiguation)

attrs: role class

values:
note
hatnote navigation-not-searchable

attrs: class style

attrs: alt src width height class data-file-width data-file-height

value: mi

values: vertical-align: -2.838ex; width:22.845ex; height:6.176ex;

values: vertical-align: -0.338ex; width:6.038ex; height:2.176ex;

values:
/wiki/Point_charge
mw-redirect
Point charge

values:
/wiki/Linear_functional
mw-redirect
Linear functional

values:
/wiki/Complex_analytic
mw-redirect
Complex analytic

values: vertical-align: -0.338ex; width:5.132ex; height:2.176ex;

values: vertical-align: -0.838ex; width:7.06ex; height:2.843ex;

values: vertical-align: -0.838ex; width:5.495ex; height:2.843ex;

values: vertical-align: -0.838ex; width:3.739ex; height:2.843ex;

values: vertical-align: -0.338ex; width:3.448ex; height:2.676ex;

values: vertical-align: -0.338ex; width:2.796ex; height:2.343ex;

attrs: columnalign rowspacing columnspacing displaystyle

values: vertical-align: -2.505ex; width:31.751ex; height:6.176ex;

values: vertical-align: -3.671ex; margin-left: -0.079ex; width:46.198ex; height:8.509ex;

values: vertical-align: -0.338ex; width:7.594ex; height:2.343ex;

values: vertical-align: -2.505ex; margin-left: -0.079ex; width:19.421ex; height:6.176ex;

values: vertical-align: -0.671ex; width:3.715ex; height:2.509ex;

values: vertical-align: -2.171ex; width:8.559ex; height:4.009ex;

values: {\displaystyle \lim _{\Delta t\to 0}F_{\Delta t}}

values: vertical-align: -2.338ex; width:19.239ex; height:5.843ex;

values: {\displaystyle \int _{-\infty }^{\infty }F_{\Delta t}(t)\,dt=P,}

values: vertical-align: -0.338ex; width:7.078ex; height:2.343ex;

values: vertical-align: -2.171ex; width:27.384ex; height:4.176ex;

values: {\displaystyle F(t)=P\delta (t)=\lim _{\Delta t\to 0}F_{\Delta t}(t)}

values:
/wiki/Gaussian_distribution
mw-redirect
Gaussian distribution

values:
/wiki/Measure_theory
mw-redirect
Measure theory

values:
/wiki/Fourier_integral_theorem
mw-redirect
Fourier integral theorem

values: vertical-align: -2.338ex; width:48.129ex; height:5.843ex;

values: vertical-align: -2.338ex; width:37.918ex; height:5.843ex;

values:
/wiki/Augustin_Cauchy
mw-redirect
Augustin Cauchy

values: vertical-align: -2.505ex; width:46.151ex; height:6.176ex;

values: right left right left right left right left right left right left

values: 0em 2em 0em 2em 0em 2em 0em 2em 0em 2em 0em

values: vertical-align: -5.671ex; width:70.403ex; height:12.509ex;

values: vertical-align: -2.338ex; width:31.832ex; height:5.843ex;

values:
/wiki/Rapidly_decreasing
mw-redirect
Rapidly decreasing

values:
/wiki/Augustin_Louis_Cauchy
mw-redirect
Augustin Louis Cauchy

values:
/wiki/Lord_Kelvin
mw-redirect
Lord Kelvin

values: vertical-align: -2.505ex; width:22.598ex; height:6.176ex;

values: \delta (x)={\begin{cases}+\infty ,&x=0\\0,&x\neq 0\end{cases}}

values: vertical-align: -2.338ex; width:16.981ex; height:5.843ex;

values: \int _{-\infty }^{\infty }\delta (x)\,dx=1.

values:
/wiki/Lebesgue_integral
mw-redirect
Lebesgue integral

values: vertical-align: -2.338ex; width:23.078ex; height:5.843ex;

values: \int _{-\infty }^{\infty }f(x)\,\delta \{dx\}=f(0)

values:
/wiki/Absolutely_continuous
mw-redirect
Absolutely continuous

values:
/wiki/Radon%E2%80%93Nikodym_derivative
mw-redirect
Radon–Nikodym derivative

values: vertical-align: -2.338ex; width:23.902ex; height:5.843ex;

values: \int _{-\infty }^{\infty }f(x)\delta (x)\,dx=f(0)

values:
/wiki/Unit_step_function
mw-redirect
Unit step function

values: vertical-align: -2.505ex; width:22.593ex; height:6.176ex;

values: H(x)={\begin{cases}1&{\text{if }}x\geq 0\\0&{\text{if }}x<0.\end{cases}}

values: vertical-align: -2.338ex; width:40.223ex; height:5.843ex;

values: vertical-align: -2.338ex; width:35.232ex; height:5.843ex;

values: Characteristic function (probability theory)

values:
/wiki/Moment_generating_function
mw-redirect
Moment generating function

values:
/wiki/Smooth_function
mw-redirect
Smooth function

values:
/wiki/Compact_support
mw-redirect
Compact support

values: vertical-align: -0.838ex; width:11.547ex; height:2.843ex;

values: {\displaystyle \delta [\varphi ]=\varphi (0)}

values: border-left:none; border-top:0px none #e5e5e5; border-right:none; border-bottom:none; padding:0.08em;

values: border-collapse:collapse; background:none; margin:0; border:none; width:99%;

values: vertical-align:middle; width:99%; border:none; padding:0.08em;

values:
math_1
reference nourlexpansion
font-weight:bold; font-style:italic;

values: border-collapse:collapse; background:none; margin:0; border:none;

values: vertical-align: -3.171ex; width:33.404ex; height:7.676ex;

values:
/wiki/Distributional_derivative
mw-redirect
Distributional derivative

values: vertical-align: -2.338ex; width:28.326ex; height:5.843ex;

values: \delta [\varphi ]=-\int _{-\infty }^{\infty }\varphi '(x)H(x)\,dx.

values: vertical-align: -2.338ex; width:38.999ex; height:5.843ex;

values: vertical-align: -2.338ex; width:40.839ex; height:5.843ex;

values: vertical-align: -2.338ex; width:22.486ex; height:5.843ex;

values: \int _{\mathbf {R} ^{n}}f(\mathbf {x} )\delta \{d\mathbf {x} \}=f(\mathbf {0} )

values: vertical-align: -0.838ex; width:27.622ex; height:2.843ex;

values:
math_2
reference nourlexpansion
font-weight:bold; font-style:italic;

values:
/wiki/Sigma_algebra
mw-redirect
Sigma algebra

values: vertical-align: -2.505ex; width:24.72ex; height:6.176ex;

values:
/wiki/Differentiable_structure
mw-redirect
Differentiable structure

values: vertical-align: -1.005ex; width:14.777ex; height:3.009ex;

values: \delta _{x_{0}}[\varphi ]=\varphi (x_{0})

values:
math_3
reference nourlexpansion
font-weight:bold; font-style:italic;

values:
/wiki/Locally_compact_Hausdorff_space
mw-redirect
Locally compact Hausdorff space

values: vertical-align: -1.005ex; width:9.087ex; height:3.009ex;

values: vertical-align: -2.671ex; width:36.531ex; height:6.176ex;

values: vertical-align: -2.671ex; width:14.559ex; height:6.509ex;

values: \delta (\alpha x)={\frac {\delta (x)}{|\alpha |}}.

values:
math_4
reference nourlexpansion
font-weight:bold; font-style:italic;

values:
/wiki/Even_function
mw-redirect
Even function

values: vertical-align: -0.838ex; width:13.386ex; height:2.843ex;

values: vertical-align: -0.838ex; width:10.509ex; height:2.843ex;

values: vertical-align: -0.838ex; width:19.963ex; height:2.843ex;

values: vertical-align: -2.338ex; width:28.061ex; height:5.843ex;

values: \int _{-\infty }^{\infty }f(t)\delta (t-T)\,dt=f(T).

values: vertical-align: -0.838ex; width:16.242ex; height:2.843ex;

values: vertical-align: -2.338ex; width:28.143ex; height:5.843ex;

values: vertical-align: -4.005ex; width:27.018ex; height:8.843ex;

values: =\int \limits _{-\infty }^{\infty }f(\tau )\delta (\tau -(t-T))\,d\tau

values: vertical-align: -0.838ex; width:11.588ex; height:2.843ex;

values: vertical-align: -2.338ex; width:35.92ex; height:5.843ex;

values: Edit section: Composition with a function

values: vertical-align: -2.838ex; width:48.037ex; height:6.343ex;

values:
/wiki/Continuously_differentiable
mw-redirect
Continuously differentiable

values: vertical-align: -0.671ex; width:4.391ex; height:2.676ex;

values:
/wiki/Root_of_a_function
mw-redirect
Root of a function

values: vertical-align: -2.671ex; width:21.27ex; height:6.509ex;

values: vertical-align: -3.005ex; width:24.111ex; height:6.843ex;

values:
/wiki/Simple_root
mw-redirect
Simple root

values: vertical-align: -2.671ex; width:43.382ex; height:6.009ex;

values: vertical-align: -3.005ex; width:35.313ex; height:6.843ex;

values: Edit section: Properties in n dimensions

values: vertical-align: -0.838ex; width:18.553ex; height:3.009ex;

values: \delta (\alpha \mathbf {x} )=|\alpha |^{-n}\delta (\mathbf {x} )

values: vertical-align: -0.838ex; width:13.6ex; height:2.843ex;

values: \delta (\rho \mathbf {x} )=\delta (\mathbf {x} ).

values:
/wiki/Lipschitz_function
mw-redirect
Lipschitz function

values: vertical-align: -2.838ex; width:54.939ex; height:6.343ex;

values: vertical-align: -2.838ex; width:40.964ex; height:6.676ex;

values: vertical-align: -2.338ex; width:20.755ex; height:5.843ex;

values: \delta _{S}[g]=\int _{S}g(\mathbf {s} )\,d\sigma (\mathbf {s} )

values:
/wiki/Simple_layer_potential
mw-redirect
Simple layer potential

values:
/wiki/Normal_derivative
mw-redirect
Normal derivative

values: vertical-align: -2.338ex; width:40.044ex; height:6.176ex;

values: vertical-align: -2.338ex; width:30.922ex; height:5.843ex;

values: vertical-align: -0.838ex; width:4.189ex; height:2.843ex;

values:
/wiki/Schwartz_functions
mw-redirect
Schwartz functions

values: vertical-align: -0.338ex; width:1.313ex; height:3.009ex;

values: vertical-align: -0.838ex; width:14.306ex; height:3.509ex;

values: vertical-align: -0.338ex; width:6.252ex; height:3.009ex;

values: vertical-align: -0.338ex; width:10.05ex; height:2.343ex;

values:
/wiki/Linear_time-invariant_system
mw-redirect
Linear time-invariant system

values: vertical-align: -2.338ex; width:23.077ex; height:5.843ex;

values: \int _{-\infty }^{\infty }1\cdot e^{2\pi ix\xi }\,d\xi =\delta (x)

values: vertical-align: -0.838ex; width:22.447ex; height:3.009ex;

values: \langle 1,f^{\vee }\rangle =f(0)=\langle \delta ,f\rangle

values: vertical-align: -2.338ex; width:57.431ex; height:5.843ex;

values: vertical-align: -0.838ex; width:12.924ex; height:3.176ex;

values: vertical-align: -0.838ex; width:18.66ex; height:3.509ex;

values: vertical-align: -2.338ex; width:27.064ex; height:5.843ex;

values: Edit section: Distributional derivatives

values: vertical-align: -0.838ex; width:24.936ex; height:3.009ex;

values: \delta '[\varphi ]=-\delta [\varphi ']=-\varphi '(0).

values: vertical-align: -2.338ex; width:40.847ex; height:5.843ex;

values: vertical-align: -0.838ex; width:22.908ex; height:3.343ex;

values: vertical-align: -2.171ex; width:28.86ex; height:6.009ex;

values: \delta '(x)=\lim _{h\to 0}{\frac {\delta (x+h)-\delta (x)}{h}}.

values: vertical-align: -2.171ex; width:19.884ex; height:5.509ex;

values: \delta '=\lim _{h\to 0}{\frac {1}{h}}(\tau _{h}\delta -\delta )

values: vertical-align: -0.838ex; width:20.018ex; height:2.843ex;

values: (\tau _{h}S)[\varphi ]=S[\tau _{-h}\varphi ].

values: vertical-align: -7.171ex; width:20.943ex; height:15.509ex;

values: vertical-align: -0.671ex; width:19.436ex; height:2.843ex;

values: vertical-align: -0.838ex; width:12.707ex; height:2.843ex;

values: vertical-align: -1.338ex; width:67.661ex; height:3.843ex;

values: vertical-align: -3.671ex; width:17.494ex; height:6.176ex;

values: S=\sum _{|\alpha |\leq m}c_{\alpha }\partial ^{\alpha }\delta _{a}.

values: Edit section: Representations of the delta function

values: vertical-align: -2.338ex; width:18.186ex; height:4.343ex;

values: {\displaystyle \delta (x)=\lim _{\varepsilon \to 0^{+}}\eta _{\varepsilon }(x),}

values: vertical-align: -2.338ex; width:29.864ex; height:5.843ex;

values:
math_5
reference nourlexpansion
font-weight:bold; font-style:italic;

values: Edit section: Approximations to the identity

values: vertical-align: -2.005ex; width:19.477ex; height:5.009ex;

values: vertical-align: -2.005ex; width:19.641ex; height:5.009ex;

values: vertical-align: -0.838ex; width:22.131ex; height:2.676ex;

values: f*\eta _{\varepsilon }\to f\quad {\rm {{as\ }\varepsilon \to 0.}}

values:
/wiki/Mean_convergence
mw-redirect
Mean convergence

values: vertical-align: -3.089ex; margin-bottom: -0.249ex; width:30.048ex; height:7.843ex;

values:
/wiki/Hat_function
mw-redirect
Hat function

values: vertical-align: -2.005ex; width:29.939ex; height:5.009ex;

values: Edit section: Probabilistic considerations

values: vertical-align: -2.671ex; width:43.29ex; height:6.509ex;

values: vertical-align: -3.171ex; width:38.381ex; height:7.509ex;

values: vertical-align: -0.838ex; width:13.839ex; height:2.176ex;

values: \eta _{\varepsilon }*\eta _{\delta }=\eta _{\varepsilon +\delta }

values:
/wiki/Partial_differential_equations
mw-redirect
Partial differential equations

values: vertical-align: -3.505ex; width:31.055ex; height:8.176ex;

values: vertical-align: -2.838ex; width:19.816ex; height:6.676ex;

values: vertical-align: -2.171ex; width:13.872ex; height:6.009ex;

values: vertical-align: -2.838ex; width:20.411ex; height:6.176ex;

values: vertical-align: -2.338ex; width:38.68ex; height:5.843ex;

values:
/wiki/Laplace_equation
mw-redirect
Laplace equation

values:
/wiki/Electrostatic_potential
mw-redirect
Electrostatic potential

values: vertical-align: -2.505ex; width:26.406ex; height:7.509ex;

values:
/wiki/Fourier_multiplier
mw-redirect
Fourier multiplier

values: vertical-align: -3.838ex; width:36.295ex; height:8.843ex;

values: Hyperbolic partial differential equations

values:
/wiki/Gas_dynamics
mw-redirect
Gas dynamics

values: vertical-align: -2.505ex; width:16.481ex; height:6.176ex;

values: vertical-align: -5.338ex; width:31.246ex; height:11.843ex;

values: vertical-align: -3.005ex; width:42.409ex; height:7.509ex;

values: vertical-align: -2.505ex; width:24.576ex; height:6.176ex;

values: vertical-align: -0.838ex; width:9.303ex; height:2.843ex;

values: vertical-align: -0.838ex; width:9.073ex; height:2.843ex;

values: vertical-align: -0.838ex; width:8.706ex; height:2.843ex;

values: vertical-align: -0.838ex; width:11.708ex; height:2.843ex;

values: vertical-align: -6.171ex; width:51.829ex; height:13.509ex;

values:
/wiki/Sphere_measure
mw-redirect
Sphere measure

values: vertical-align: -2.338ex; width:30.708ex; height:6.343ex;

values: vertical-align: -2.338ex; width:49.037ex; height:6.343ex;

values: vertical-align: -2.338ex; width:85.894ex; height:6.176ex;

values: vertical-align: -2.671ex; width:29.113ex; height:6.176ex;

values: R\varphi (\xi ,p)=\int _{x\cdot \xi =p}f(x)\,d^{n-1}x.

values: vertical-align: -2.671ex; width:34.235ex; height:6.509ex;

values: vertical-align: -2.671ex; width:39.566ex; height:6.176ex;

values: vertical-align: -3.171ex; width:39.221ex; height:8.843ex;

values: vertical-align: -3.171ex; width:37.333ex; height:7.509ex;

values: vertical-align: -2.338ex; width:26.198ex; height:5.843ex;

values: vertical-align: -2.338ex; width:40.026ex; height:5.843ex;

values: vertical-align: -3.005ex; width:20.522ex; height:7.009ex;

values: \delta (x)={\frac {1}{2\pi }}\sum _{n=-\infty }^{\infty }e^{inx}

values:
/wiki/Summability_methods
mw-redirect
Summability methods

values: vertical-align: -3.171ex; width:41.881ex; height:8.343ex;

values: vertical-align: -2.338ex; width:27.565ex; height:5.843ex;

values: \int _{\mathbf {R} }F_{N}(x)f(x)\,dx\to 2\pi f(0)

values:
/wiki/Densely_defined
mw-redirect
Densely defined

values:
/wiki/Square_integrable_function
mw-redirect
Square integrable function

values:
/wiki/Bounded_linear_functional
mw-redirect
Bounded linear functional

values:
/wiki/Sobolev_embedding_theorem
mw-redirect
Sobolev embedding theorem

values: vertical-align: -2.338ex; width:38.491ex; height:5.843ex;

values: vertical-align: -1.005ex; width:24.121ex; height:3.009ex;

values:
/wiki/Continuous_dual_space
mw-redirect
Continuous dual space

values: Edit section: Spaces of holomorphic functions

values: vertical-align: -2.505ex; width:33.139ex; height:6.343ex;

values: vertical-align: -2.505ex; width:32.31ex; height:6.343ex;

values: Edit section: Resolutions of the identity

values: vertical-align: -3.005ex; width:14.288ex; height:7.009ex;

values: f=\sum _{n=1}^{\infty }\alpha _{n}\varphi _{n}.

values: vertical-align: -0.838ex; width:13.41ex; height:2.843ex;

values: \alpha _{n}=\langle \varphi _{n},f\rangle ,

values: vertical-align: -0.838ex; width:10.536ex; height:3.509ex;

values:
/wiki/Dyadic_tensor
mw-redirect
Dyadic tensor

values: vertical-align: -3.005ex; width:19.18ex; height:7.009ex;

values: f=\sum _{n=1}^{\infty }\varphi _{n}\left(\varphi _{n}^{\dagger }f\right).

values:
/wiki/Identity_operator
mw-redirect
Identity operator

values: vertical-align: -3.005ex; width:14.214ex; height:7.009ex;

values: I=\sum _{n=1}^{\infty }\varphi _{n}\varphi _{n}^{\dagger },

values: vertical-align: -0.838ex; width:6.17ex; height:3.509ex;

values: vertical-align: -3.005ex; width:36.241ex; height:7.009ex;

values: vertical-align: -2.338ex; width:26.138ex; height:5.843ex;

values: f(x)=\int \,\delta (x-\xi )f(\xi )\,d\xi ,

values: vertical-align: -3.005ex; width:27.193ex; height:7.009ex;

values: \delta (x-\xi )=\sum _{n=1}^{\infty }\varphi _{n}(x)\varphi _{n}^{*}(\xi ).

values: Edit section: Infinitesimal delta functions

values: vertical-align: -2.338ex; width:20.881ex; height:5.843ex;

values:
/wiki/Infinite_sum
mw-redirect
Infinite sum

values: vertical-align: -3.005ex; width:22.971ex; height:7.009ex;

values: \Delta (x)=\sum _{n=-\infty }^{\infty }\delta (x-n),

values: vertical-align: -3.005ex; width:28.526ex; height:7.009ex;

values: (f*\Delta )(x)=\sum _{n=-\infty }^{\infty }f(x-n).

values: vertical-align: -0.838ex; width:22.524ex; height:3.509ex;

values: (f*\Delta )^{\wedge }={\hat {f}}{\widehat {\Delta }}={\hat {f}}\Delta

values: vertical-align: -2.838ex; width:35.853ex; height:6.676ex;

values: vertical-align: -2.338ex; width:31.591ex; height:5.676ex;

values: vertical-align: -2.838ex; width:53.435ex; height:6.676ex;

values: Edit section: Relationship to the Kronecker delta

values: vertical-align: -2.505ex; width:16.535ex; height:6.176ex;

values: vertical-align: -0.838ex; width:6.951ex; height:2.843ex;

values: vertical-align: -3.005ex; width:16.048ex; height:7.009ex;

values: \sum _{i=-\infty }^{\infty }a_{i}\delta _{ik}=a_{k}.

values: vertical-align: -2.338ex; width:31.041ex; height:5.843ex;

values:
/wiki/Discrete_distribution
mw-redirect
Discrete distribution

values:
/wiki/Continuous_distribution
mw-redirect
Continuous distribution

values: vertical-align: -3.005ex; width:23.193ex; height:7.009ex;

values: vertical-align: -2.838ex; width:38.262ex; height:6.509ex;

values: vertical-align: -2.338ex; width:27.349ex; height:6.176ex;

values: vertical-align: -2.338ex; width:39.426ex; height:6.176ex;

values: vertical-align: -0.838ex; width:15.285ex; height:2.843ex;

values: \langle \varphi _{n}|\varphi _{m}\rangle =\delta _{nm}

values: vertical-align: -1.338ex; width:14.042ex; height:3.843ex;

values: vertical-align: -0.838ex; width:12.12ex; height:2.843ex;

values: vertical-align: -1.338ex; width:17.359ex; height:3.843ex;

values: I=\sum |\varphi _{n}\rangle \langle \varphi _{n}|.

values: vertical-align: -1.005ex; width:17.77ex; height:3.009ex;

values: {\displaystyle \varphi _{y}(x)=\delta (x-y).}

values: vertical-align: -1.005ex; width:8.435ex; height:3.009ex;

values: vertical-align: -1.005ex; width:11.863ex; height:2.843ex;

values: {\displaystyle P\varphi _{y}=y\varphi _{y}.}

values: vertical-align: -1.005ex; width:20.468ex; height:3.176ex;

values: \langle \varphi _{y},\varphi _{y'}\rangle =\delta (y-y')

values: vertical-align: -2.338ex; width:23.47ex; height:5.843ex;

values: \psi (x)=\int _{\Omega }c(y)\varphi _{y}(x)\,dy

values: vertical-align: -1.005ex; width:14.765ex; height:3.009ex;

values: c(y)=\langle \psi ,\varphi _{y}\rangle .

values: vertical-align: -2.338ex; width:18.893ex; height:5.843ex;

values: I=\int _{\Omega }|\varphi _{y}\rangle \,\langle \varphi _{y}|\,dy

values: vertical-align: -2.171ex; width:20.093ex; height:6.176ex;

values: m{\frac {\mathrm {d} ^{2}\xi }{\mathrm {d} t^{2}}}+k\xi =I\delta (t),

values:
/wiki/Spring_constant
mw-redirect
Spring constant

values: vertical-align: -2.171ex; width:15.86ex; height:6.176ex;

values: {\displaystyle EI{\frac {\mathrm {d} ^{4}w}{\mathrm {d} x^{4}}}=q(x),}

values: vertical-align: -0.838ex; width:19.25ex; height:2.843ex;

values: vertical-align: -9.509ex; margin-bottom: -0.329ex; width:37.251ex; height:20.843ex;

values:
nofollow
external text
//arxiv.org/abs/1302.0864

values:
nofollow
external text
//doi.org/10.1007%2FJHEP11%282012%29032

values:
nofollow
external text
//doi.org/10.1007%2FBF00329867

values:
nofollow
external text
//doi.org/10.3934%2Fcpaa.2008.7.229

values:
nofollow
external text
//arxiv.org/abs/quant-ph/0109154

values:
nofollow
external text
//doi.org/10.1036%2F0071455469

values:
/wiki/Elias_Stein
mw-redirect
Elias Stein

values:
nofollow
external text
//doi.org/10.1063%2F1.2339017

values:
nofollow
external text
//doi.org/10.1063%2F1.2771422

values: Generalized inverse Gaussian distribution

values: Maximum entropy probability distribution

values: nowraplinks collapsible autocollapse navbox-inner

values:
/wiki/Method_of_indivisibles
mw-redirect
Method of indivisibles

values:
navbox-image
6
width:1px;padding:0px 0px 0px 2px

values: Analyse des Infiniment Petits pour l'Intelligence des Lignes Courbes

values: Elementary Calculus: An Infinitesimal Approach

values: navigation navbox Infinitesimals padding:3px

values: Discussion about edits from this IP address [ctrl-alt-n]

values: A list of edits made from this IP address [ctrl-alt-y]

values: You are encouraged to create an account and log in; however, it is not mandatory

values: You're encouraged to log in; however, it's not mandatory. [ctrl-alt-o]

values:
/wiki/Dirac_delta_function
View the content page [ctrl-alt-c]
c

values: Discussion about the content page [ctrl-alt-t]

values: p-namespaces navigation vectorTabs p-namespaces-label

values:
p-variants
navigation
vectorMenu emptyPortlet
p-variants-label

values: Past revisions of this page [ctrl-alt-h]

values: p-views navigation vectorTabs p-views-label

values:
p-cactions
navigation
vectorMenu emptyPortlet
p-cactions-label

values:
mw-wiki-logo
/wiki/Main_Page
Visit the main page

values:
/wiki/Main_Page
Visit the main page [ctrl-alt-z]
z

values: Featured content – the best of Wikipedia

values: Find background information on current events

values: دالة ديراك – Arabic

values: Дэльта-функцыя – Belarusian

values: Κρουστική συνάρτηση – Greek

values: تابع دلتای دیراک – Persian

values: 디랙 델타 함수 – Korean

values: डिरैक डेल्टा फलन – Hindi

values: פונקציית דלתא של דיראק – Hebrew

values: დირაკის დელტა ფუნქცია – Georgian

values: ディラックのデルタ関数 – Japanese

values: Дельта-функция – Russian

values: ඩිරැක් ඩෙල්ටා ශ්‍රිතය – Sinhala

values: Дельта-функция – Tatar

values: ดิแรกเดลตาฟังก์ชัน – Thai

values: Дельта-функція Дірака – Ukrainian

values: 狄拉克δ函数 – Chinese

values: interlanguage-link interwiki-zh badge-Q17437798 badge-goodarticle

values: /static/images/wikimedia-button-1.5x.png 1.5x, /static/images/wikimedia-button-2x.png 2x

values: /static/images/poweredby_mediawiki_132x47.png 1.5x, /static/images/poweredby_mediawiki_176x62.png 2x

values:
suggestions
display: none; font-size: 13px;

values: mediawiki ltr sitedir-ltr mw-hide-empty-elt ns-0 ns-subject page-Dirac_delta_function rootpage-Dirac_delta_function skin-vector action-view