在科学和数学 中,狄拉克δ 函数 或简称δ 函数 (译名德尔塔函数 、得耳他函数 )是在实数线上定义的一个广义函数 或分布 。它在除零以外的点上都等于零,且其在整个定义域上的积分 等于1。[ 2] δ 函数有时可看作是在原点处无限高、无限细,但是总面积为1的一个尖峰,在物理上代表了理想化的质点 或点电荷 的密度。
狄拉克δ 函数示意图。直线上箭头的高度一般用于指定δ 函数前任何乘法常数的值,亦即等于函数下方的面积。另一种惯例是把面积值写在箭头的旁边。
狄拉克δ 函数是以零为中心的正态分布
δ
a
(
x
)
=
1
a
π
e
−
x
2
/
a
2
{\displaystyle \delta _{a}(x)={\frac {1}{a{\sqrt {\pi }}}}\mathrm {e} ^{-x^{2}/a^{2}}}
随
a
→
0
{\displaystyle a\rightarrow 0}
的(分布 意义上的)极限 。
从纯数学的观点来看,狄拉克δ 函数并非严格意义上的函数 ,因为任何在扩展实数线 上定义的函数,如果在一个点以外的地方都等于零,其总积分必须为零。δ 函数只有在出现在积分以内的时候才有实质的意义。根据这一点,δ 函数一般可以当做普通函数一样使用。
狄拉克δ 函数得名自物理学家 保罗·狄拉克 ,其形式上所遵守的规则属于运算微积分 的一部分,是物理学和工程学的标准工具。包括δ 函数在内的运算微积分方法,在20世纪初受到数学家的质疑,直到1950年代洛朗·施瓦茨 才发展出一套令人满意的严谨理论。严谨地来说,δ 函数必须定义为一个分布,对应于支撑集 为原点的概率测度 。
在许多应用中,均将δ 视为由在原点处有尖峰的函数所组成的序列 的极限(弱极限 ),而序列中的函数则可作为对δ 函数的近似。在信号处理 上,δ 函数常称为单位脉冲符号 或单位脉冲函数 。克罗内克δ 函数 是对应于狄拉克δ 函数的离散函数,其定义域为离散集,值域可以是0或者1。
δ 函数的图形 通常可以视为整条x 轴和正y 轴。虽然称为函数,但δ 函数并非真正的函数,至少它的值域不在实数 以内。例如,f (x ) = δ (x ) 和g (x ) = 0 这两个数学对象 除了在x = 0 以外都有相同的值,但其积分却不相同。根据勒贝格积分 理论,若f 和g 为函数,使得f = g 几乎处处 成立,则f 可积当且仅当 g 可积且f 和g 的积分相同。要严谨处理δ 函数,须用到测度论 或分布。
δ 函数可以代表一个既高又窄的尖峰函数(脉冲),用以描述点电荷 和质点 等抽象化 的概念。举例来说,要描述球杆击球的动力学 问题,可以用δ 函数描述击球那一刻的力 。不但各种方程会因此简化,而且只需球杆传递的总冲量就能算出球击出后的运动,而不须考虑球杆向球传递能量的复杂具体情况。
在应用数学中,δ 函数往往能看作是某函数序列的极限(弱极限),该序列中的每一项都在原点处有一个尖峰,例如以零为中心、方差 趋向零的高斯分布 序列。
约瑟夫·傅里叶 在他的《热分析理论》(法语:Théorie analytique de la chaleur )中呈现了以下的方程,今天称为傅里叶积分定理:
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 )\ ,}
这相当于以这种方式引入了δ 函数:
δ
(
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 )\ .}
之后,奥古斯丁·路易·柯西 用指数函数表达了这一定理:
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.}
柯西指出,在某些情况下,积分的计算顺序会影响计算结果。[ a]
分布理论 允许重新排列柯西方程,使它更接近以上傅里叶的方程:
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}}}
其中δ 函数可表达为:
δ
(
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\ .}
在接下来的几个世纪,数学家才逐渐理解这一指数形式的严谨含义,以及方程中的函数f 所需的条件。用旧有的数学观念来理解,会有以下的问题:
经典傅里叶变换的最大缺点在于,能够有效计算的只有很狭窄的一类函数。这些函数必须(在无限的邻域 内)足够快地降至零,才能保证傅里叶积分值的存在。例如,连多项式 这种如此简单的函数,也不存在经典意义上的傅里叶变换。经典傅里叶变换扩展至分布,大大增加了能够进行变换的函数类型,移除了诸多障碍。
之后,米歇尔·普朗歇尔 (Michel Plancherel)开创性的L 2 理论(1910年)、诺伯特·维纳 和萨洛蒙·博赫纳 (Salomon Bochner)的贡献(1930年前后)以及最后洛朗·施瓦茨 归纳这一切的分布理论(1945年)进一步推广了傅里叶积分,并建立了狄拉克δ 函数的严格定义。
1827年,柯西首次明确写下一个无限高的单位脉冲函数(柯西分布 的无限小 版本)。西莫恩·德尼·帕松 和古斯塔夫·基尔霍夫 之后在研究波传播的时候,考虑过这一函数。基尔霍夫与赫尔曼·冯·亥姆霍兹 将单位脉冲描述为高斯分布的极限,这也符合开尔文勋爵 对点热源的描述。19世纪末,奥利弗·亥维赛 利用形式上的傅里叶级数 对单位脉冲进行操作。[ b] 1930年,保罗·狄拉克 在影响深远的《量子力学原理》中引入了δ 函数作为一种“方便的记号”,故此该函数今天以他命名。
笼统地来说,δ 函数是在实数线上的一个函数,在原点上无限,在所有其他点上为零,
δ
(
x
)
=
{
+
∞
,
x
=
0
0
,
x
≠
0
{\displaystyle \delta (x)={\begin{cases}+\infty ,&x=0\\0,&x\neq 0\end{cases}}}
并同时满足以下条件
∫
−
∞
∞
δ
(
x
)
d
x
=
1.
{\displaystyle \int _{-\infty }^{\infty }\delta (x)\,dx=1.}
[ 19]
这只是一个概略的表述:δ 函数并不是一个严格意义上的函数,没有任何定义在实数线上的函数能满足以上的条件。更严谨地来说,δ 函数可以定义为分布 或测度 。
测度 是其中一种严谨定义δ 函数的方法。作为一个测度,δ 函数取一个实线R 的子集A ,当0 ∈ A 时输出δ (A ) = 1 ,否则δ (A ) = 0 。如果把δ 函数想象成位于0的一个理想化的质点,则δ (A )代表集合A 所包含的质量。一个函数相对于δ 的积分便可以定义为相对于这个测度的勒贝格积分 。对于所有连续紧支撑函数f ,这一积分满足:
∫
−
∞
∞
f
(
x
)
δ
{
d
x
}
=
f
(
0
)
.
{\displaystyle \int _{-\infty }^{\infty }f(x)\,\delta \{dx\}=f(0).}
测度δ 相对于勒贝格测度 不绝对连续 ,它其实是一个奇异测度 。因此,它并不具有拉东-尼科迪姆导数 ,也就是不存在满足以下条件的函数δ :
∫
−
∞
∞
f
(
x
)
δ
(
x
)
d
x
=
f
(
0
)
.
{\displaystyle \int _{-\infty }^{\infty }f(x)\delta (x)\,dx=f(0).}
虽然这种写法仍非常常见,但是它实际上只是一种方便的记号,而不是任何有良好定义的(黎曼 或勒贝格 )积分。
作为R 上的概率测度 ,狄拉克测度可以通过它的累积分布函数 ──单位阶跃函数 ──来定义:[ c]
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}}}
换句话说,H (x )是积累指示函数 1 (−∞, x ] 相对于测度δ 的积分:
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].}
δ 函数相对于一个连续函数的积分可以通过黎曼-斯蒂尔杰斯积分 严格定义:
∫
−
∞
∞
f
(
x
)
δ
{
d
x
}
=
∫
−
∞
∞
f
(
x
)
d
H
(
x
)
.
{\displaystyle \int _{-\infty }^{\infty }f(x)\delta \{dx\}=\int _{-\infty }^{\infty }f(x)\,dH(x).}
δ 函数的所有更高矩 都是零。其特征函数 和矩母函数 都等于1。
在分布理论 中,一个广义函数并不像普通函数一样直接定义,而是在它相对其他函数积分的时候,以它如何影响这一积分来定义。沿着这条思路,只须定义δ 函数相对某个足够“良好”的测试函数的“积分”就足够了。如果δ 函数已经定义为测度,则这种积分可以是测试函数相对于这δ 测度的勒贝格积分。
测试函数空间一般可包括所有
R
{\displaystyle \mathbf {R} }
上的紧支撑 光滑函数 。作为一个分布,δ 函数是在测试函数空间上的线性泛函 ,定义为
∀
φ
,
δ
[
φ
]
=
φ
(
0
)
.
{\displaystyle \forall \varphi ,\delta [\varphi ]=\varphi (0).}
1
若要使δ 成为一个正式的分布,它必须要在测试函数空间上相对某个合适拓扑 为连续的。在测试函数空间上的线性泛函
S
{\displaystyle S}
要能够良好定义一个分布,其必要和充分条件是,对于每个正整数
N
{\displaystyle N}
,有整数
M
N
{\displaystyle M_{N}}
和常数
C
N
{\displaystyle C_{N}}
,使得对每个测试函数
φ
{\displaystyle \varphi }
,以下不等式都成立:[ 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)|.}
当
S
{\displaystyle S}
就是δ 分布时,对所有的
N
{\displaystyle N}
取
C
N
=
1
,
M
N
=
0
{\displaystyle C_{N}=1,M_{N}=0}
,就能满足这条不等式。因此,δ 是级数为零的分布。它也是一个紧支撑分布,其支撑集 是
{
0
}
{\displaystyle \{0\}}
。
δ分布有几种等价的定义。例如,它是亥维赛阶跃函数 的分布导数,也就是说,对于任何测试函数
φ
{\displaystyle \varphi }
,
δ
[
φ
]
=
−
∫
−
∞
∞
φ
′
(
x
)
H
(
x
)
d
x
.
{\displaystyle \delta [\varphi ]=-\int _{-\infty }^{\infty }\varphi '(x)H(x)\,\mathrm {d} x.}
直观而言,如果允许分部积分法 ,以上的积分就会简化为
∫
−
∞
∞
φ
(
x
)
H
′
(
x
)
d
x
=
∫
−
∞
∞
φ
(
x
)
δ
(
x
)
d
x
.
{\displaystyle \int _{-\infty }^{\infty }\varphi (x)H'(x)\,\mathrm {d} x=\int _{-\infty }^{\infty }\varphi (x)\delta (x)\,\mathrm {d} x.}
而利用黎曼-斯蒂尔杰斯积分 的分部积分法可以得到
−
∫
−
∞
∞
φ
′
(
x
)
H
(
x
)
d
x
=
∫
−
∞
∞
φ
(
x
)
d
H
(
x
)
.
{\displaystyle -\int _{-\infty }^{\infty }\varphi '(x)H(x)\,\mathrm {d} x=\int _{-\infty }^{\infty }\varphi (x)\,\mathrm {d} H(x).}
在测度论中,狄拉克测度通过积分产生分布。相反,公式(1 )在所有紧支撑连续函数
φ
{\displaystyle \varphi }
空间上定义了一个Daniell积分 ,且根据里斯表示定理 ,该积分可以表示为
φ
{\displaystyle \varphi }
相对于某拉东测度 的勒贝格积分 。
当讲到“狄拉克δ 函数”时,一般指的是分布,而不是测度。因此一些文献也会称之为“狄拉克δ 分布”。测度论中相对应的概念则称为狄拉克测度 。
在n 维欧几里得空间 R n 中,狄拉克δ 函数可以定义为一个测度,使得对于所有紧支撑连续函数f ,满足
∫
R
n
f
(
x
)
δ
{
d
x
}
=
f
(
0
)
.
{\displaystyle \int _{\mathbf {R} ^{n}}f(\mathbf {x} )\delta \{d\mathbf {x} \}=f(\mathbf {0} ).}
作为一个测度,n 维δ 函数是每个独立变量的1维δ 函数的积测度 。也就是说,若x = (x 1 , x 2 , ..., x n ) ,则
δ
(
x
)
=
δ
(
x
1
)
δ
(
x
2
)
⋯
δ
(
x
n
)
.
{\displaystyle \delta (\mathbf {x} )=\delta (x_{1})\delta (x_{2})\cdots \delta (x_{n}).}
2
上文的1维δ 分布也存在n 维推广。[ 25] 尽管在物理学和工程学中应用广泛,公式(2 )还是必须小心使用,因为多个分布的积只有在较狭窄的条件下才有良好的定义。[ 27]
狄拉克测度这个概念可以定义在任何集合上。设X 是集合,x 0 ∈ X ,Σ为X 子集上的任何σ-代数 ,则对每个集合A ∈ Σ 可以定义测度:
δ
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}}}
这就是单位质量集中在x 0 处的狄拉克测度。
δ 函数也可以推广至微分流形 ,由于具有微分结构 ,因此能保留它作为分布的一些性质。流形M 上以x 0 ∈ M 为中心的δ 函数可定义为以下分布:
δ
x
0
[
φ
]
=
φ
(
x
0
)
{\displaystyle \delta _{x_{0}}[\varphi ]=\varphi (x_{0})}
3
对于所有M 上的紧支撑光滑实数值函数φ 。一个常见的特殊情况是,M 是欧几里得空间R n 中的一个开集 。
在局部紧 豪斯多夫空间 X 中,集中在点x 的狄拉克测度是对应于对紧支撑连续函数φ 的Daniell积分(3 )的拉东测度 。推广到这一层次,已经无法进行普通的微积分,不过仍然可以使用抽象分析中的许多工具。例如,映射
x
0
↦
δ
x
0
{\displaystyle x_{0}\mapsto \delta _{x_{0}}}
是把X 嵌入 到包含所有在X 上的有限拉东测度的空间(具有淡拓扑)的一个连续函数。而且,X 在这一嵌入下的值域的凸包 ,在在X 上的概率测度空间中是一个稠密集 。
对非零标量
α
{\displaystyle \alpha }
,δ 函数有以下缩放性质:
∫
−
∞
∞
δ
(
α
x
)
d
x
=
∫
−
∞
∞
δ
(
u
)
d
u
|
α
|
=
1
|
α
|
{\displaystyle \int _{-\infty }^{\infty }\delta (\alpha x)\,\mathrm {d} x=\int _{-\infty }^{\infty }\delta (u)\,{\frac {\mathrm {d} u}{|\alpha |}}={\frac {1}{|\alpha |}}}
所以
δ
(
α
x
)
=
δ
(
x
)
|
α
|
.
{\displaystyle \delta (\alpha x)={\frac {\delta (x)}{|\alpha |}}.}
4
δ 函数是一个偶分布 ,也就是说
δ
(
−
x
)
=
δ
(
x
)
.
{\displaystyle \delta (-x)=\delta (x).}
因此δ 函数属于−1阶齐次函数 。
δ 分布在
n
{\displaystyle n}
维空间中的缩放性质如下:
δ
(
α
x
)
=
|
α
|
−
n
δ
(
x
)
{\displaystyle \delta (\alpha \mathbf {x} )=|\alpha |^{-n}\delta (\mathbf {x} )}
从而δ 是一个
n
{\displaystyle n}
阶齐次分布 。δ 分布在任何镜射 和旋转
ρ
{\displaystyle \rho }
下不变:
δ
(
ρ
x
)
=
δ
(
x
)
.
{\displaystyle \delta (\rho \mathbf {x} )=\delta (\mathbf {x} ).}
和单变量时一样,可以唯一地定义δ 与双利普希茨函数
g
:
R
n
→
R
n
{\displaystyle g:\mathbf {R} ^{n}\to \mathbf {R} ^{n}}
的复合,使得恒等式
∫
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} )|\,\mathrm {d} \mathbf {x} =\int _{g(\mathbf {R} ^{n})}\delta (\mathbf {u} )f(\mathbf {u} )\,\mathrm {d} \mathbf {u} }
对于所有紧支撑函数
f
{\displaystyle f}
成立。
利用几何测度论中的余面积公式,可以定义δ 函数与从一个欧几里得空间到另一个不同维度的空间的浸没的复合,所产生的结果是一种流 。在连续可微函数
g
:
R
n
→
R
n
{\displaystyle g:\mathbf {R} ^{n}\to \mathbf {R} ^{n}}
满足
g
{\displaystyle g}
的梯度 处处非零的特殊情况下,以下恒等式成立:[ 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} ))\,\mathrm {d} \mathbf {x} =\int _{g^{-1}(0)}{\frac {f(\mathbf {x} )}{|\mathbf {\nabla } g|}}\,\mathrm {d} \sigma (\mathbf {x} )}
其中右边的积分范围是
g
−
1
(
0
)
{\displaystyle g^{-1}(0)}
,即
g
(
x
)
=
0
{\displaystyle g(\mathbf {x} )=0}
所定义的一个
(
n
−
1
)
{\displaystyle (n-1)}
维曲面(依闵可夫斯基容度)。这叫做单层(simple layer)积分。
更一般地来说,若
S
{\displaystyle S}
是
R
n
{\displaystyle \mathbf {R} ^{n}}
中的光滑超曲面 ,则可以把
S
{\displaystyle S}
联系到在
S
{\displaystyle S}
上对任何紧支撑光滑函数
g
{\displaystyle g}
积分的分布:
δ
S
[
g
]
=
∫
S
g
(
s
)
d
σ
(
s
)
{\displaystyle \delta _{S}[g]=\int _{S}g(\mathbf {s} )\,\mathrm {d} \sigma (\mathbf {s} )}
其中
σ
{\displaystyle \sigma }
是联系到
S
{\displaystyle S}
的超曲面测度。这种推广与
S
{\displaystyle S}
上的单层位势的位势论 相关。设
D
{\displaystyle D}
是
R
n
{\displaystyle \mathbf {R} ^{n}}
中具有光滑边缘
S
{\displaystyle S}
的区域 ,则
δ
S
{\displaystyle \delta _{S}}
等于
D
{\displaystyle D}
的指示函数 的(分布意义上的)法向导数 :
−
∫
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}}\;\mathrm {d} \mathbf {x} =\int _{S}\,g(\mathbf {s} )\;\mathrm {d} \sigma (\mathbf {s} ),}
其中
n
{\displaystyle n}
是向外法线 。[ 36]
δ 函数属于缓增函数,所以拥有良好定义的傅里叶变换 。正式地说,(在一些傅里叶变换惯例下)有
δ
^
(
ξ
)
=
∫
−
∞
∞
e
−
2
π
i
x
ξ
δ
(
x
)
d
x
=
1.
{\displaystyle {\hat {\delta }}(\xi )=\int _{-\infty }^{\infty }e^{-2\pi ix\xi }\delta (x)\,dx=1.}
一个分布的傅里叶变换的定义是,在缓增分布与速降函数 的对偶配对
⟨
⋅
,
⋅
⟩
{\displaystyle \langle \cdot ,\cdot \rangle }
下,要求傅里叶变换是自伴 的。从而,
δ
^
{\displaystyle {\hat {\delta }}}
定义为满足以下条件的唯一缓增分布:
⟨
δ
^
,
φ
⟩
=
⟨
δ
,
φ
^
⟩
{\displaystyle \langle {\hat {\delta }},\varphi \rangle =\langle \delta ,{\hat {\varphi }}\rangle }
对于一切速降函数φ 。从此可推导,的确
δ
^
=
1
{\displaystyle {\hat {\delta }}=1}
。
这条恒等式意味着,δ 函数与任何其他缓增分布S 的卷积 即等于S :
S
∗
δ
=
S
.
{\displaystyle S*\delta =S.}
这意味着,δ 是缓增分布上的卷积的单位元 。而且,在卷积下的紧支撑分布空间是一个以δ 函数为单位元的结合代数 。这在信号处理 应用中尤其重要,因为与缓增分布的卷积属于线性时不变系统 ,而基于δ 函数的线性时不变系统可以测量该缓增分布的冲激响应 。只要对δ 作适当的近似,就可以以任意要求的程度算出冲激响应。一旦知道冲激响应,就能完全描述整个系统的特征。详见线性时不变系统理论:冲激响应和卷积 。
缓增函数f (ξ ) = 1的反傅里叶变换等于δ 函数。更正式地表达,
∫
−
∞
∞
1
⋅
e
2
π
i
x
ξ
d
ξ
=
δ
(
x
)
.
{\displaystyle \int _{-\infty }^{\infty }1\cdot e^{2\pi ix\xi }\,d\xi =\delta (x).}
更严谨地来说,有
⟨
1
,
f
∨
⟩
=
f
(
0
)
=
⟨
δ
,
f
⟩
{\displaystyle \langle 1,f^{\vee }\rangle =f(0)=\langle \delta ,f\rangle }
对于一切速降函数f 。
这样,δ 函数暗示著在R 上的傅里叶核的正交性,即:
∫
−
∞
∞
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}).}
换句话说,缓增分布
f
(
t
)
=
e
i
2
π
ξ
1
t
{\displaystyle f(t)=e^{i2\pi \xi _{1}t}}
的傅里叶变换是
f
^
(
ξ
2
)
=
δ
(
ξ
1
−
ξ
2
)
.
{\displaystyle {\hat {f}}(\xi _{2})=\delta (\xi _{1}-\xi _{2}).}
这同样可以通过对傅里叶变换要求自伴性而得出。
利用傅里叶变换的解析延拓 ,可以得出δ 函数的拉普拉斯变换 :
∫
0
∞
δ
(
t
−
a
)
e
−
s
t
d
t
=
e
−
s
a
.
{\displaystyle \int _{0}^{\infty }\delta (t-a)e^{-st}\,dt=e^{-sa}.}
狄拉克δ 分布的分布导数是一个分布δ ′,它对于所有紧支撑光滑测试函数φ 定义为[ 38]
δ
′
[
φ
]
=
−
δ
[
φ
′
]
=
−
φ
′
(
0
)
.
{\displaystyle \delta '[\varphi ]=-\delta [\varphi ']=-\varphi '(0).}
此处第一个等号类似于分部积分,因为若δ 是个真正的函数,则
∫
−
∞
∞
δ
′
(
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.}
δ 的k 阶导数的定义大同小异,对任何测试函数φ ,
δ
(
k
)
[
φ
]
=
(
−
1
)
k
φ
(
k
)
(
0
)
.
{\displaystyle \delta ^{(k)}[\varphi ]=(-1)^{k}\varphi ^{(k)}(0).}
从而δ 是个无限可微分布。
δ 函数的一阶微分是差商
δ
′
(
x
)
=
lim
h
→
0
δ
(
x
+
h
)
−
δ
(
x
)
h
{\displaystyle \delta '(x)=\lim _{h\to 0}{\frac {\delta (x+h)-\delta (x)}{h}}}
的分布极限。[ 39] 更准确地说,有
δ
′
=
lim
h
→
0
1
h
(
τ
h
δ
−
δ
)
{\displaystyle \delta '=\lim _{h\to 0}{\frac {1}{h}}(\tau _{h}\delta -\delta )}
其中τh 是平移算子,对于函数的定义是τ h φ (x ) = φ (x + h ) ,而对于分布S 的定义是
(
τ
h
S
)
[
φ
]
=
S
[
τ
−
h
φ
]
.
{\displaystyle (\tau _{h}S)[\varphi ]=S[\tau _{-h}\varphi ].}
在电磁学 中,δ 函数的一阶导数代表一个位于原点的点磁偶极 ,因此也称为“偶极”或偶函数 。
δ 函数的导数满足一些基本性质,包括:
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}}}
这些性质都可以通过对测试函数积分,并运用分部积分法推导而出。
另外,δ ′与紧支撑光滑测试函数f 的卷积为
δ
′
∗
f
=
δ
∗
f
′
=
f
′
,
{\displaystyle \delta '*f=\delta *f'=f',}
这来自卷积的分布导数的性质。
更一般地说,在n 维欧几里得空间 R n 中的开集 U 上,以点a ∈ U 为中心的狄拉克δ 分布定义为[ 41]
δ
a
[
φ
]
=
φ
(
a
)
{\displaystyle \delta _{a}[\varphi ]=\varphi (a)}
对于一切φ ∈ S (U ) ,其中S (U )是所有U 上的紧支撑光滑函数的空间。若α = (α 1 , ..., α n ) 是任何多重指标 ,而∂α 表示相关的混合偏导数 算子,则δ a 的α 阶导数∂α δ a 是[ 41]
⟨
∂
α
δ
a
,
φ
⟩
=
(
−
1
)
|
α
|
⟨
δ
a
,
∂
α
φ
⟩
=
(
−
1
)
|
α
|
∂
α
φ
(
x
)
|
x
=
a
{\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}}
对于一切φ ∈ S (U ) 。也就是说,δ a 的α 阶导数是个分布,它在任何测试函数φ 的值等于φ 在点a 的α 阶导数(加上合适的正负号)。
δ 函数的一阶偏导数可以视为沿着坐标平面的双层 。更一般地来讲,支撑在一个曲面上的单层的法向导数 是在该曲面上的双层,并表示一个层磁单极。δ 函数的更高阶导数,在物理学里称为多极 。
高阶导数很自然地能够建构具有单元素支撑集的分布的完整结构。若S 是任何在U 上、支撑集为一个点{a }的分布,那么存在整数m 和一组系数c α ,使得[ 41]
S
=
∑
|
α
|
≤
m
c
α
∂
α
δ
a
.
{\displaystyle S=\sum _{|\alpha |\leq m}c_{\alpha }\partial ^{\alpha }\delta _{a}.}
δ 函数可以视为一个函数序列的极限:
δ
(
x
)
=
lim
ε
→
0
+
η
ε
(
x
)
,
{\displaystyle \delta (x)=\lim _{\varepsilon \to 0^{+}}\eta _{\varepsilon }(x),\,}
其中ηε (x )有时称为初生δ 函数 。这一极限是个弱极限:对于一切紧支撑连续函数 f ,有
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
或者,这个极限对于一切紧支撑光滑函数 f 都存在。这两种不同的弱收敛模式往往有十分微妙的差异,前者是依测度的淡拓扑收敛,而后者则是分布 的收敛。
通常一个初生δ 函数ηε 可以如下建构。设η 是一个R 上的总积分为1的绝对可积函数,并定义
η
ε
(
x
)
=
ε
−
1
η
(
x
ε
)
.
{\displaystyle \eta _{\varepsilon }(x)=\varepsilon ^{-1}\eta \left({\frac {x}{\varepsilon }}\right).}
在n 维当中,改用以下缩放
η
ε
(
x
)
=
ε
−
n
η
(
x
ε
)
.
{\displaystyle \eta _{\varepsilon }(x)=\varepsilon ^{-n}\eta \left({\frac {x}{\varepsilon }}\right).}
在简单的变量更换之后,可见ηε 的积分同样等于1。不难证明,(5 )对于一切连续紧支撑函数f 都成立,从而ηε 作为一个测度向δ 弱收敛。
这样建构的ηε 叫做对单位元的近似 ,因为包含所有绝对可积函数的空间L 1 (R )在函数卷积 这一作用下闭合:f ∗ g ∈ L 1 (R ) ,当f 和g 都属于L 1 (R )。然而,L 1 (R )在卷积下并没有单位元:没有任何元素h 使得f ∗ h = f 对于所有f 都成立。但序列ηε 仍然能够近似这种单位元,就是说
f
∗
η
ε
→
f
a
s
ε
→
0.
{\displaystyle f*\eta _{\varepsilon }\to f\quad {\rm {{as\ }\varepsilon \to 0.}}}
在平均收敛 (即L 1 中的收敛)下,此极限存在。要确保几乎处处 点收敛,还需要对ηε 加上更多的前提,比如它必须是对应于一个紧支撑函数的柔化函数 。
如果最初的η = η 1 本身已经是光滑的而且具有紧支撑,那么整个序列就叫做柔化序列。标准柔化序列可以通过选择某个适当归一化的脉冲函数 η 来定义,例如
η
(
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}}}
在数值分析 等的一些情况下,以分段线性函数 对单位元进行近似会更加有用。可以定义η 1 是一个三角形函数 ,然后有
η
ε
(
x
)
=
ε
−
1
max
(
1
−
|
x
ε
|
,
0
)
{\displaystyle \eta _{\varepsilon }(x)=\varepsilon ^{-1}\max \left(1-\left|{\frac {x}{\varepsilon }}\right|,0\right)}
它们全部都是连续的且具有紧支撑,但不是光滑的,从而也不是柔化函数。
在概率论 中,很自然地会加入一项额外条件,亦即要求对单位元近似的初始η 1 是正的,因为这样它就会代表一个概率分布 。有时需要和概率分布 来卷积的原因是,输出值是输入值的凸组合 ,因此处于最高与最低输入值之间,从而能够避免过冲 或下冲。取η 1 为任何随意的概率分布,并如上设ηε (x ) = η 1 (x /ε )/ε ,就可以取得对单位元的近似。此外,如果η 的平均值为0,且更高矩 较小,那么序列就会向δ 函数收敛得更快。比如,设η 1 是[−1/2, 1/2] 上的均匀分布 ,即矩形函数 ,则:[ 45]
η
ε
(
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}}}
又以维格纳半圆分布 举例,
η
ε
(
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}}}
这是连续的,且具有紧支撑,但因为不是光滑的,所以也不是柔化函数。
初生δ 函数往往以卷积半群 的身份出现。这相等于另加一个条件──ηε 与ηδ 的卷积必须满足
η
ε
∗
η
δ
=
η
ε
+
δ
{\displaystyle \eta _{\varepsilon }*\eta _{\delta }=\eta _{\varepsilon +\delta }}
对于所有ε 和δ > 0 。L 1 中形成初生δ 函数的卷积半群一定是对单位元的近似(用上式的意思),但半群设下了限制性颇强的条件。
在实践当中,对单位元的半群近似出现在物理学所启发的椭圆型 和抛物型偏微分方程 中,作为基本解或格林函数 。在应用数学 中,半群是线性时不变系统 的输出。抽象地来说,若A 是一个作用在x 的函数上的线性算子,则在对初值问题 求解时会出现卷积半群:
{
∂
∂
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}}}
其中极限同样是弱极限。设ηε (x ) = η (ε , x ) ,得出相关的初生δ 函数。
下面将列出若干有物理意义的、在此类基本解中所出现的卷积半群。
热核
热核 ,定义是
η
ε
(
x
)
=
1
2
π
ε
e
−
x
2
2
ε
{\displaystyle \eta _{\varepsilon }(x)={\frac {1}{\sqrt {2\pi \varepsilon }}}\mathrm {e} ^{-{\frac {x^{2}}{2\varepsilon }}}}
代表的是一条无限长丝在t > 0时的温度,如果在t = 0时在原点处储藏了一个单位的热能。此半群依照一维热传导方程 演变:
∂
u
∂
t
=
1
2
∂
2
u
∂
x
2
.
{\displaystyle {\frac {\partial u}{\partial t}}={\frac {1}{2}}{\frac {\partial ^{2}u}{\partial x^{2}}}.}
在概率论中,ηε (x )是一个平均值为0、方差 为ε 的正态分布 。它代表了一个粒子从原点开始,作标准布朗运动 ,在t = ε 这一刻的概率密度函数 。在这种情况下,半群条件也就体现了布朗运动的马尔可夫性质 。
在高维欧几里得空间R n 中,热核等于
η
ε
=
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 }}},}
并且在作必要的修改后有着相同的物理解释。而且,当ε → 0 ,ηε → δ ,因此它也代表了一个初生δ 函数。
帕松核
帕松核
η
ε
(
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 }
是拉普拉斯方程 在上半平面中的基本解。它代表了一个半无限平板在边缘的电势 固定为δ 函数时的电势 。帕松核也和柯西分布 有紧密的联系。半群根据以下方程演变:
∂
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)}
其中的算子严谨地定义为傅里叶乘数
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 ).}
在波 动力学等物理范畴中,须要解不少双曲型偏微分方程 ,而这些方程有更多的奇异解。因此,从相关的柯西问题 的基本解所产生的初生δ 函数,一般都是震荡积分。例如,从穿音速 气体动力学 的欧拉-特里科米方程 的解,取得归一化艾里函数
ε
−
1
3
Ai
(
x
ε
−
1
3
)
.
{\displaystyle \varepsilon ^{-{\frac {1}{3}}}\operatorname {Ai} \left(x\varepsilon ^{-{\frac {1}{3}}}\right).}
虽然用的是傅里叶变换,但是不难看到,这从某种意义上产生了一个半群,只不过由于它不是绝对可积的,所以无法从上文更强的意义来定义半群。许多用振荡积分来建构的初生δ 函数,只能从分布意义上收敛(见下文的例子狄利克雷核 ),而不能从测度意义上收敛。
另一个例子是R 1+1 中的波动方程 的柯西问题:[ 48]
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}}}
解u 代表了一条具有无限弹性 的弦,一开始在原点处受到扰动后,距离平衡的位移程度。
其他同类的对单位元的近似还包括,广泛应用于电子和电信中的sinc函数
η
ε
(
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}
以及贝塞尔函数
η
ε
(
x
)
=
1
ε
J
1
ε
(
x
+
1
ε
)
.
{\displaystyle \eta _{\varepsilon }(x)={\frac {1}{\varepsilon }}J_{\frac {1}{\varepsilon }}\left({\frac {x+1}{\varepsilon }}\right).}
要对线性偏微分方程
L
[
u
]
=
f
{\displaystyle L[u]=f}
求解,其中L 是R n 上的一个微分算子 ,可以先取得基本解,即对以下方程求解:
L
[
u
]
=
δ
.
{\displaystyle L[u]=\delta .}
当L 相对简单的时候,通常直接利用傅里叶变换就可以求解(如上文提到的帕松核和热核)。如果算子比较复杂,可以先考虑更简单的方程
L
[
u
]
=
h
{\displaystyle L[u]=h}
其中h 是一个平面波 函数,意思是对于某矢量ξ,有
h
=
h
(
x
⋅
ξ
)
{\displaystyle h=h(x\cdot \xi )}
这样的方程可以用柯西-柯瓦列夫斯卡娅定理 求解(如果L 的系数是解析函数 ),或者用求积法求解(如果L 的系数是常数)。所以,如果δ 函数能够分解成平面波的话,理论上就能取得线性偏微分方程的解。
对δ 函数进行平面波分解的方法,最早是约翰·拉东 (Johann Radon)所发展的一套通用技巧之一,之后由弗瑞兹·约翰 进一步发展至这种形式(1955 )。[ d] 选择k ,使得n + k 是一个偶整数;对于实数s ,定义
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}}}
要取得δ ,对g (x · ξ ) 相对球体测度dω积分,再对积分施用拉普拉斯算子 的幂,其中ξ 属于单位球面 S n −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 }.}
此处的拉普拉斯算子理解为弱导数,所以上式的意思是,对于任何测试函数φ ,
φ
(
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 }.}
这一结果来自牛顿位势 公式(牛顿位势是帕松方程的基本解)。这实际上是拉东变换 的逆转公式,因为它能够从φ (x )在超曲面上的积分取回φ (x )的值。例如,若n 是奇数,而k = 1 ,则右边的积分等于
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}
其中Rφ (ξ , p ) 是φ 的拉东变换:
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.}
平面波分解的另一个等价表达式是(Gelfand & Shilov 1966–1968 ,I, §3.10 harvnb模板错误: 无指向目标: CITEREFGelfandShilov1966–1968 (帮助 ) )
δ
(
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 }}
对于偶数n ,且
δ
(
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 }}
对于奇数n 。
在对傅里叶级数 的研究当中,一个重要的问题是须要判断,和某周期函数 相关的傅里叶级数是否收敛到该函数,以及从何种意义上收敛。周期为2π的函数f 的傅里叶级数的第n 部分和,定义是与狄利克雷核 在区间[−π,π]上的卷积:
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)}}.}
因此
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}}
其中
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.}
傅里叶级数的一个基础结果说明,当N → ∞ ,狄利克雷核趋向于δ 函数的倍数。收敛指的是从分布意义上的收敛,即
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)}
对于一切紧支撑光滑函数f 。从而,在区间[−π,π]上有
δ
(
x
)
=
1
2
π
∑
n
=
−
∞
∞
e
i
n
x
.
{\displaystyle \delta (x)={\frac {1}{2\pi }}\sum _{n=-\infty }^{\infty }e^{inx}.}
尽管如此,但这并不对所有紧支撑连续函数成立,换言之,DN 从测度意义上并不弱收敛。鉴于傅里叶级数无法收敛,因此数学家建立了各种可和性方法 来达到收敛。从切萨罗求和法 发展出费耶核
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}.}
费耶核从一种更强的意义向δ 函数收敛:[ e]
∫
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)}
对于一切紧支撑连续函数f 。结果是,所有连续函数的傅里叶级数在每一点上都是切萨罗可和的,且和的值等于该函数的值。
狄拉克δ 分布是在包含所有平方可积函数 的希尔伯特空间 L 2 上所稠密定义的一个无界 线性泛函 。紧支撑光滑函数在L 2 中是一个稠密集 ,且δ 分布对于紧支撑光滑函数有良好定义。在许多应用中,可以对L 2 的某个子空间赋予更强的拓扑 ,使得δ 函数能够定义一个有界线性算子 。
索伯列夫空间
索伯列夫嵌入定理 应用在实数线R 上的索伯列夫空间 上时,意味着任何平方可积函数f ,只要满足
‖
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 }
就必定是连续的,而且满足
δ
[
f
]
=
|
f
(
0
)
|
<
C
‖
f
‖
H
1
.
{\displaystyle \delta [f]=|f(0)|<C\|f\|_{H^{1}}.}
从而δ 是一个在索伯列夫空间H 1 上的有界线性泛函。另一个等价的说法是,δ 是H 1 的连续对偶空间 H −1 的元素。更一般地说,在n 维中,有δ ∈ H −s (R n ) ,条件是s > n / 2 。
在复分析 中,δ 函数出现在柯西积分公式 中,公式说明,若D 是复平面 上一个具有光滑边缘的域,则
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}
对于一切在D 的闭包 内连续的全纯函数 。从而,对于此类全纯函数,δ 函数δ z 可以以柯西积分表示:
δ
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}}.}
更一般地来说,设H 2 (∂D )是一个哈代空间 ,它是所有在D 中直到D 的边缘都是连续的全纯函数,在L 2 (∂D ) 中的凸包。H 2 (∂D )中的函数可以唯一地延续成D 上的全纯函数,而且柯西积分公式仍然成立。特别是对于z ∈ D ,δ 函数δ z 是一个H 2 (∂D )上的连续线性泛函。这是多复数变量函数中的特殊情况:对于光滑域D ,塞格核 代替了柯西积分的角色。
在可分 希尔伯特空间 中,给定一个由函数{φ n }组成的标准正交基 (如一个紧自伴算子的归一化特征矢量 ),那么任何矢量f 都可以表达成:
f
=
∑
n
=
1
∞
α
n
φ
n
.
{\displaystyle f=\sum _{n=1}^{\infty }\alpha _{n}\varphi _{n}.}
系数{α n }可以如下得出:
α
n
=
⟨
φ
n
,
f
⟩
,
{\displaystyle \alpha _{n}=\langle \varphi _{n},f\rangle ,}
也可以写为
α
n
=
φ
n
†
f
,
{\displaystyle \alpha _{n}=\varphi _{n}^{\dagger }f,}
这是狄拉克符号 的一种。[ f] 在这种写法下,f 以并矢 方式展开:
f
=
∑
n
=
1
∞
φ
n
(
φ
n
†
f
)
.
{\displaystyle f=\sum _{n=1}^{\infty }\varphi _{n}\left(\varphi _{n}^{\dagger }f\right).}
设I 是该希尔伯特空间上的恒等算子 ,则表达式
I
=
∑
n
=
1
∞
φ
n
φ
n
†
{\displaystyle I=\sum _{n=1}^{\infty }\varphi _{n}\varphi _{n}^{\dagger }}
称为单位分解。当希尔伯特空间是L 2 (D ),包含所有在域D 上的平方可积函数,那么
φ
n
φ
n
†
{\displaystyle \varphi _{n}\varphi _{n}^{\dagger }}
就是一个积分算子,而f 可以重新表达为:
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 .}
右边的级数是在L 2 当中向f 收敛。就算f 是连续函数,点收敛极限也不一定存在。尽管如此,往往可以滥用符号,写
f
(
x
)
=
∫
δ
(
x
−
ξ
)
f
(
ξ
)
d
ξ
,
{\displaystyle f(x)=\int \,\delta (x-\xi )f(\xi )\,d\xi ,}
如此来表示δ 函数:
δ
(
x
−
ξ
)
=
∑
n
=
1
∞
φ
n
(
x
)
φ
n
∗
(
ξ
)
.
{\displaystyle \delta (x-\xi )=\sum _{n=1}^{\infty }\varphi _{n}(x)\varphi _{n}^{*}(\xi ).}
在适当的装备希尔伯特空间 (Φ, L 2 (D ), Φ*) 中,其中Φ ⊂ L 2 (D ) 包含所有紧支撑光滑函数,视乎基φ n 的性质,上方的级数有可能在Φ*中收敛。在大多数实际情况下,标准正交基来自于某个积分或微分算子,这时候级数会从分布的意义上收敛。
1827年柯西在若干论文中写下无限高、无限窄的函数时,用到了一个无限小数α ,使得函数δα 满足[ g]
∫
F
(
x
)
δ
α
(
x
)
=
F
(
0
)
.
{\displaystyle \int F(x)\delta _{\alpha }(x)=F(0).}
他(以及拉扎尔·卡诺 )把无限小数定义为一个趋向于零的序列。
非标准分析 能够严谨地处理无限小数。利用超实数 的语言,狄拉克δ 函数可以用含有无限小数的延伸实数来表达,详见Yamashita (2007) 论文所列出的相关书目。这样定义的δ 函数是真正意义上的函数,使得对于每个实函数F ,都有
∫
F
(
x
)
δ
α
(
x
)
=
F
(
0
)
,
{\displaystyle \int F(x)\delta _{\alpha }(x)=F(0),}
结果与傅里叶和柯西用别的语言所表达的一样。
狄拉克梳子是无限个相距T 的δ 函数
由一系列狄拉克测度组成的均匀脉冲串叫做狄拉克梳子 ,亦以西里尔字母Ш 相形称为Ш分布,是一个采样 函数,常用在数码信号处理 和离散时间信号分析中。狄拉克梳子是许多单个δ 函数的无限和,和的极限是分布意义上的极限:
Δ
(
x
)
=
∑
n
=
−
∞
∞
δ
(
x
−
n
)
,
{\displaystyle \Delta (x)=\sum _{n=-\infty }^{\infty }\delta (x-n),}
这可以理解为,在每个整数处都有一个单位质点。
狄拉克梳子是其自身的傅里叶变换(或乘以某个归一常数)。其重要性在于,若f 是速降函数 ,则f 在周期化后的结果以以下卷积表示:
(
f
∗
Δ
)
(
x
)
=
∑
n
=
−
∞
∞
f
(
x
−
n
)
.
{\displaystyle (f*\Delta )(x)=\sum _{n=-\infty }^{\infty }f(x-n).}
特别有,
(
f
∗
Δ
)
∧
=
f
^
Δ
^
=
f
^
Δ
,
{\displaystyle (f*\Delta )^{\wedge }={\hat {f}}{\widehat {\Delta }}={\hat {f}}\Delta ,}
这正是帕松求和公式。[ 54]
索霍茨基-魏尔斯特拉斯定理 是量子力学中重要的定理,它把δ 函数和分布p.v.1/x 联系起来,后者是函数1/x 的柯西主值 ,定义是
⟨
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.}
索霍茨基公式说明,
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),}
此处的极限是分布意义上的极限,就是说对于一切紧支撑光滑函数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.}
克罗内克δ 函数 δij 的定义是,
δ
i
j
=
{
1
i
=
j
0
i
≠
j
{\displaystyle \delta _{ij}={\begin{cases}1&i=j\\0&i\not =j\end{cases}}}
对于所有整数i 、j 。它满足以下的筛选性质:若
(
a
i
)
i
∈
Z
{\displaystyle (a_{i})_{i\in \mathbf {Z} }}
是一个两头无限的序列,则
∑
i
=
−
∞
∞
a
i
δ
i
k
=
a
k
.
{\displaystyle \sum _{i=-\infty }^{\infty }a_{i}\delta _{ik}=a_{k}.}
这和狄拉克δ 函数的筛选性质十分相似:对于任何R 上的实函数或复函数f ,有
∫
−
∞
∞
f
(
x
)
δ
(
x
−
x
0
)
d
x
=
f
(
x
0
)
.
{\displaystyle \int _{-\infty }^{\infty }f(x)\delta (x-x_{0})\,dx=f(x_{0}).}
也就是说,克罗内克δ 函数可以看作是与狄拉克δ 函数对应的离散函数。
在概率论 和统计学 中,狄拉克函数往往以概率密度函数 的身份,来代表一个离散分布 或部分离散、部分连续 的分布(概率密度函数一般只用作描述完全连续分布)。例如,设一组点x = {x 1 , ..., xn },对应概率为p 1 , ..., pn ;由这些点所组成的离散分布的概率密度函数可以写作
f
(
x
)
=
∑
i
=
1
n
p
i
δ
(
x
−
x
i
)
.
{\displaystyle f(x)=\sum _{i=1}^{n}p_{i}\delta (x-x_{i}).}
又举一例,设一个分布,其中十分之六的情况下输出标准正态分布 ,而十分之四的情况下输出单个数值3.5,这是一个部分连续、部分离散的混合分布。其密度函数是
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).}
也可以以完全不同的方法,用δ 函数表示扩散过程 (如布朗运动 )中的局部时 。一个随机过程 的局部时的表达式为
ℓ
(
x
,
t
)
=
∫
0
t
δ
(
x
−
B
(
s
)
)
d
s
{\displaystyle \ell (x,t)=\int _{0}^{t}\delta (x-B(s))\,ds}
这代表了该过程在某特定区间内,在点x 所花的时间。更准确地说,当只有一个维度时,上面的积分可以写成
ℓ
(
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}
其中1 [x −ε , x +ε ] 是区间[x −ε , x +ε ] 的指示函数 。
以下举一个例子,展示δ 函数如何在量子力学 中派上用场。一个粒子 的波函数 所给出的,是粒子出现在特定空间范围内的概率幅 。波函数假定属于希尔伯特空间L 2 (平方可积函数 空间),且粒子在某空间范围内出现的总概率,等于波函数的绝对值平方在该范围内的积分。一组波函数{φ n }叫做标准正交,如果
⟨
φ
n
|
φ
m
⟩
=
δ
n
m
{\displaystyle \langle \varphi _{n}|\varphi _{m}\rangle =\delta _{nm}}
其中δ 指的是克罗内克δ 函数,而不是狄拉克δ 函数。一组标准正交波函数叫做在平方可积函数空间中完备,如果任何波函数ψ 都可以表达为一些φ n 的线性组合:
ψ
=
∑
c
n
φ
n
,
{\displaystyle \psi =\sum c_{n}\varphi _{n},}
其中
c
n
=
⟨
φ
n
|
ψ
⟩
{\displaystyle c_{n}=\langle \varphi _{n}|\psi \rangle }
。量子力学中的哈密顿算符 量度的是(束缚态的)能级,而算符的所有本征函数 正正就组成了波函数完备标准正交系统,每个本征函数所对应的特征值等于能量值。这组能量值叫做这个哈密尔顿算符的光谱 。利用狄拉克符号 (如上 ),上式所表达的就是单位分解:
I
=
∑
|
φ
n
⟩
⟨
φ
n
|
.
{\displaystyle I=\sum |\varphi _{n}\rangle \langle \varphi _{n}|.}
此处,特征值都是离散的,但一个可观察量 的特征值也可以是连续的,就如位置算符 Qψ (x ) = xψ (x ) 。位置(在一维当中)的光谱是整条实数线,所以称作连续光谱 。不过,和哈密顿算符不同的是,位置算符并没有正式的本征函数。为了解决这一困局,通常会扩大所允许使用的函数,从普通的函数到所有分布。换言之,量子力学的希尔伯特空间要由合适的装备希尔伯特空间 取代。这样一来,位置算符就有了一套完备的本征分布,对应于实数线上的每个点y :
φ
y
(
x
)
=
δ
(
x
−
y
)
.
{\displaystyle \varphi _{y}(x)=\delta (x-y).}
位置的本征函数(分布)用狄拉克符号记作
φ
y
=
|
y
⟩
{\displaystyle \varphi _{y}=|y\rangle }
,亦称为位置本征态。
这种处理方法同样可以应用于动量算符 ,以及一切希尔伯特空间上的自伴无界算子 P ,前提是P 具有连续光谱,且不存在退化特征值。更确切地说,有一个实数的子集Ω(即算子的光谱)和一组分布φ y ,对应于Ω的每个元素y ,使得
P
φ
y
=
y
φ
y
.
{\displaystyle P\varphi _{y}=y\varphi _{y}.}
也就是说,φ y 是P 的特征矢量(本征分布)。如果这些特征矢量(作为分布)都满足归一化条件:
⟨
φ
y
,
φ
y
′
⟩
=
δ
(
y
−
y
′
)
,
{\displaystyle \langle \varphi _{y},\varphi _{y'}\rangle =\delta (y-y'),}
那么对于一切测试函数ψ,就有
ψ
(
x
)
=
∫
Ω
c
(
y
)
φ
y
(
x
)
d
y
{\displaystyle \psi (x)=\int _{\Omega }c(y)\varphi _{y}(x)\,dy}
其中
c
(
y
)
=
⟨
ψ
,
φ
y
⟩
.
{\displaystyle c(y)=\langle \psi ,\varphi _{y}\rangle .}
此处所得出的单位分解和离散的情况比较,有相似之处:
I
=
∫
Ω
|
φ
y
⟩
⟨
φ
y
|
d
y
{\displaystyle I=\int _{\Omega }|\varphi _{y}\rangle \,\langle \varphi _{y}|\,dy}
其中以算子为值的积分同样理解为弱积分。若P 的光谱同时含有连续和离散部分,则它的单位分解须包含跑遍所有离散态的和,再加上跑遍所有连续态的积分。
δ 函数在量子力学中还有众多特殊应用,例如δ 位势阱 。
在结构力学 中,δ 函数可以用来描述结构上的瞬时荷载或点荷载。一个谐振子 在t =0时突然受到冲量 为I 的力的冲击,其演变可以如下描述:
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),}
其中m 是质量,ξ是挠度,而k 是弹簧常数 。
根据欧拉﹣伯努力理论 ,一条细长的梁 的静力负荷挠度是
E
I
d
4
w
d
x
4
=
q
(
x
)
,
{\displaystyle EI{\frac {\mathrm {d} ^{4}w}{\mathrm {d} x^{4}}}=q(x),}
其中EI 是梁的弯曲刚度 ,w 是挠度 ,x 是空间坐标,而q (x )则是负荷分布。如果栋梁在x = x 0 处受到点力F 的负荷,那么负荷分布可以写作
q
(
x
)
=
F
δ
(
x
−
x
0
)
.
{\displaystyle q(x)=F\delta (x-x_{0}).}
由于δ 函数的积分是亥维赛阶跃函数 ,因此细长栋梁在多个点受到点力负荷时的静力负荷挠度,可以用一组分段多项式 来表示。
δ 函数还可以描述作用在一条梁上的点弯矩 。设两个相距d 的相反方向的点力F ,它们在栋梁上所产生的弯矩为M = Fd 。在保持M 不变的情况下,使d 趋向于零。假设所产生的弯矩位于x = 0,方向是顺时针,那么对栋梁的负荷分布就是
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}}}
因此点弯矩可以用δ 函数的导数来描述。对栋梁方程积分,得出的挠度一样是分段多项式。
要以狄拉克符号来表达这一节的内容,可见Levin 2002 ,Coordinate-space wave functions and completeness, pp.=109ff
Aratyn, Henrik; Rasinariu, Constantin, A short course in mathematical methods with Maple , World Scientific, 2006, ISBN 981-256-461-6
Arfken, G. B.; Weber, H. J., Mathematical Methods for Physicists 5th, Boston, Massachusetts: Academic Press, 2000, ISBN 978-0-12-059825-0
Bartle, Robert G. The elements of integration and Lebesgue measure . Wiley Classics Library. New York: John Wiley & Sons Inc. 1995. xii+179. ISBN 0-471-04222-6 . MR 1312157 .
Bracewell, R., The Fourier Transform and Its Applications 2nd, McGraw-Hill, 1986
Córdoba, A., La formule sommatoire de Poisson, C.R. Acad. Sci. Paris, Series I, 1988, 306 : 373–376
Courant, Richard ; Hilbert, David , Methods of Mathematical Physics, Volume II, Wiley-Interscience, 1962
Davis, Howard Ted; Thomson, Kendall T, Linear algebra and linear operators in engineering with applications in Mathematica , Academic Press, 2000, ISBN 0-12-206349-X
Debnath, Lokenath; Bhatta, Dambaru. Integral Transforms And Their Applications 2nd. CRC Press. 2007: 2 [2016-12-08 ] . ISBN 1-58488-575-0 . (原始内容存档 于2021-04-19).
Dennis, G. Zill; Warren, S. Wright. 4.5. Advanced Engineering Mathematics . 2009: 230 [2016-12-08 ] . ISBN 9780763782412 . (原始内容存档 于2021-04-19).
Dieudonné, Jean , Treatise on analysis. Vol. II, New York: Academic Press [Harcourt Brace Jovanovich Publishers], 1976, ISBN 978-0-12-215502-4 , MR 0530406
Dieudonné, Jean , Treatise on analysis. Vol. III, Boston, Massachusetts: Academic Press, 1972, MR 0350769
Dirac, Paul , The Principles of Quantum Mechanics 4th, Oxford at the Clarendon Press, 1958, ISBN 978-0-19-852011-5
Driggers, Ronald G., Encyclopedia of Optical Engineering, CRC Press, 2003, ISBN 978-0-8247-0940-2
Duistermaat, Hans; Kolk, Distributions: Theory and applications, Springer, 2010
Federer, Herbert, Geometric measure theory, Die Grundlehren der mathematischen Wissenschaften 153 , New York: Springer-Verlag: xiv+676, 1969, ISBN 978-3-540-60656-7 , MR 0257325
Fourier, J.B. The Analytical Theory of Heat English translation by Alexander Freeman, 1878. The University Press. 1822: 408 [2016-12-08 ] . (原始内容存档 于2021-04-19). , cf p 449 and pp 546–551. 法文原始文本请参见这里 (页面存档备份 ,存于互联网档案馆 )
Gelfand, I. M. ; Shilov, G. E., Generalized functions 1–5 , Academic Press, 1966
Grattan-Guinness, Ivor. Convolutions in French Mathematics, 1800–1840: From the Calculus and Mechanics to Mathematical Analysis and Mathematical Physics, Volume 2 . Birkhäuser. 2009: 653 [2016-12-08 ] . ISBN 3-7643-2238-1 . (原始内容存档 于2021-04-19).
Hartmann, William M., Signals, sound, and sensation , Springer, 1997 [2016-12-06 ] , ISBN 978-1-56396-283-7 , (原始内容存档 于2021-04-19)
Hewitt, E; Stromberg, K, Real and abstract analysis, Springer-Verlag, 1965
Hörmander, L. , The analysis of linear partial differential operators I, Grundl. Math. Wissenschaft. 256 , Springer, 1983, ISBN 3-540-12104-8 , MR 0717035
Isham, C. J., Lectures on quantum theory: mathematical and structural foundations, Imperial College Press, 1995, ISBN 978-81-7764-190-5
John, Fritz, Plane waves and spherical means applied to partial differential equations, Interscience Publishers, New York-London, 1955, MR 0075429
Komatsu, Hikosaburo. Fourier's hyperfunctions and Heaviside's pseudodifferential operators. Takahiro Kawai; Keiko Fujita (编). Microlocal Analysis and Complex Fourier Analysis . World Scientific. 2002: 200 [2016-12-08 ] . ISBN 981-238-161-9 . (原始内容存档 于2021-04-19).
Kracht, Manfred; Kreyszig, Erwin. On singular integral operators and generalizations. Themistocles M. Rassias (编). Topics in Mathematical Analysis: A Volume Dedicated to the Memory of A.L. Cauchy . World Scientific. 1989: 553 [2016-12-08 ] . ISBN 9971-5-0666-1 . (原始内容存档 于2021-04-19).
Lang, Serge , Undergraduate analysis, Undergraduate Texts in Mathematics 2nd, Berlin, New York: Springer-Verlag , 1997, ISBN 978-0-387-94841-6 , MR 1476913
Lange, Rutger-Jan, Potential theory, path integrals and the Laplacian of the indicator , Journal of High Energy Physics (Springer), 2012, 2012 (11): 29–30 [2014-10-30 ] , Bibcode:2012JHEP...11..032L , arXiv:1302.0864 , doi:10.1007/JHEP11(2012)032 , (原始内容存档 于2021-04-19)
Laugwitz, D., Definite values of infinite sums: aspects of the foundations of infinitesimal analysis around 1820, Arch. Hist. Exact Sci., 1989, 39 (3): 195–245, doi:10.1007/BF00329867
Levin, Frank S., Coordinate-space wave functions and completeness, An introduction to quantum theory , Cambridge University Press: 109ff , 2002, ISBN 0-521-59841-9
Li, Y. T.; Wong, R., Integral and series representations of the Dirac delta function, Commun. Pure Appl. Anal., 2008, 7 (2): 229–247, MR 2373214 , doi:10.3934/cpaa.2008.7.229
de la Madrid, R.; Bohm, A.; Gadella, M., Rigged Hilbert Space Treatment of Continuous Spectrum, Fortschr. Phys., 2002, 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., An Introduction to State Space, Quantum Mechanics Demystified, A Self-Teaching Guide , Demystified Series, New York: McGraw-Hill: 108, 2005-11-22 [2008-03-17 ] , ISBN 0-07-145546-9 , doi:10.1036/0071455469 , (原始内容存档 于2016-03-26)
Mitrović, Dragiša; Žubrinić, Darko. Fundamentals of Applied Functional Analysis: Distributions, Sobolev Spaces . CRC Press. 1998: 62 [2016-12-08 ] . ISBN 0-582-24694-6 . (原始内容存档 于2021-04-19).
Myint-U., Tyn; Debnath, Lokenath. Linear Partial Differential Equations for Scientists And Engineers 4th. Springer. 2007: 4 [2016-12-08 ] . ISBN 0-8176-4393-1 . (原始内容存档 于2021-04-02).
van der Pol, Balth.; Bremmer, H., Operational calculus 3rd, New York: Chelsea Publishing Co., 1987, ISBN 978-0-8284-0327-6 , MR 0904873
Rudin, W., Functional Analysis 2nd, McGraw-Hill, 1991, ISBN 0-07-054236-8
Vallée, Olivier; Soares, Manuel, Airy functions and applications to physics, London: Imperial College Press, 2004
Saichev, A I; Woyczyński, Wojbor Andrzej, Chapter1: Basic definitions and operations, Distributions in the Physical and Engineering Sciences: Distributional and fractal calculus, integral transforms, and wavelets , Birkhäuser, 1997, ISBN 0-8176-3924-1
Schwartz, L. , Théorie des distributions 1 , Hermann, 1950
Schwartz, L. , Théorie des distributions 2 , Hermann, 1951
Stein, Elias; Weiss, Guido, Introduction to Fourier Analysis on Euclidean Spaces, Princeton University Press, 1971, ISBN 0-691-08078-X
Strichartz, R., A Guide to Distribution Theory and Fourier Transforms, CRC Press, 1994, ISBN 0-8493-8273-4
Vladimirov, V. S., Equations of mathematical physics, Marcel Dekker, 1971, ISBN 0-8247-1713-9
Weisstein, Eric W. (编). Delta Function . at MathWorld --A Wolfram Web Resource. Wolfram Research, Inc. [2016-12-08 ] (英语) .
Weisstein, Eric W. (编). Doublet Function . at MathWorld --A Wolfram Web Resource. Wolfram Research, Inc. [2016-12-08 ] . (原始内容存档 于2021-04-19) (英语) .
Weisstein, Eric W. (编). Sifting Property . at MathWorld --A Wolfram Web Resource. Wolfram Research, Inc. [2016-12-08 ] . (原始内容存档 于2021-04-19) (英语) .
Yamashita, H., Pointwise analysis of scalar fields: A nonstandard approach, Journal of Mathematical Physics, 2006, 47 (9): 092301, Bibcode:2006JMP....47i2301Y , doi:10.1063/1.2339017
Yamashita, H., Comment on "Pointwise analysis of scalar fields: A nonstandard approach" [J. Math. Phys. 47, 092301 (2006)], Journal of Mathematical Physics, 2007, 48 (8): 084101, Bibcode:2007JMP....48h4101Y , doi:10.1063/1.2771422