概述
UA OPTI512R 傅立叶光学导论6 Dirac函数及其性质
- Dirac函数
- Comb函数
- 涉及Dirac函数的微分与积分
Dirac函数
在工程科学中,Dirac函数通常用来表示脉冲,所以也被称为脉冲函数,这一讲讨论Dirac函数的正式定义与常用性质。
定义 引入
s
i
n
c
sinc
sinc函数,
s
i
n
c
(
x
)
=
sin
π
x
π
x
sinc(x)=frac{sin pi x}{ pi x}
sinc(x)=πxsinπx
定义Dirac函数为
δ
(
x
)
=
lim
b
→
0
1
∣
b
∣
s
i
n
c
(
x
b
)
delta(x)=lim_{b to 0} frac{1}{|b|}sincleft( frac{x}{b} right)
δ(x)=b→0lim∣b∣1sinc(bx)
这个定义非常有用,当需要推导Dirac函数的相关结论时,我们总是可以从这个定义出发,比如因为
s
i
n
c
sinc
sinc函数是偶函数,所以Dirac函数也是偶函数,
δ
(
−
x
)
=
δ
(
x
)
delta(-x)=delta(x)
δ(−x)=δ(x)
另外,虽然这里的定义用的是sinc函数,但实际上也可以用rect、gaus等随着scale趋于0,函数值全部集中在原点处的函数去定义Dirac函数。
性质1(sifting property)
对任意well-behaved function
f
(
x
)
f(x)
f(x),
∫
x
1
x
2
f
(
x
)
δ
(
x
−
x
0
)
d
x
=
f
(
x
0
)
χ
(
x
1
,
x
2
)
(
x
0
)
int_{x_1}^{x_2} f(x)delta(x-x_0)dx =f(x_0)chi_{(x_1,x_2)}(x_0)
∫x1x2f(x)δ(x−x0)dx=f(x0)χ(x1,x2)(x0)
x
0
∈
(
x
1
,
x
2
)
x_0 in (x_1,x_2)
x0∈(x1,x2)时
χ
(
x
1
,
x
2
)
(
x
0
)
=
1
chi_{(x_1,x_2)}(x_0)=1
χ(x1,x2)(x0)=1,否则
χ
(
x
1
,
x
2
)
(
x
0
)
=
0
chi_{(x_1,x_2)}(x_0)=0
χ(x1,x2)(x0)=0。如果
f
(
x
)
f(x)
f(x)恒等于1,则
∫
x
1
x
2
f
(
x
)
δ
(
x
−
x
0
)
d
x
=
∫
x
1
x
2
δ
(
x
−
x
0
)
d
x
=
1
int_{x_1}^{x_2} f(x)delta(x-x_0)dx = int_{x_1}^{x_2} delta(x-x_0)dx=1
∫x1x2f(x)δ(x−x0)dx=∫x1x2δ(x−x0)dx=1
这说明Dirac函数围成的面积为1。并且由于
∫
x
1
x
2
f
(
x
)
δ
(
x
−
x
0
)
d
x
=
f
(
x
0
)
χ
(
x
1
,
x
2
)
(
x
0
)
=
∫
x
1
x
2
f
(
x
0
)
δ
(
x
−
x
0
)
d
x
int_{x_1}^{x_2} f(x)delta(x-x_0)dx =f(x_0)chi_{(x_1,x_2)}(x_0)=int_{x_1}^{x_2} f(x_0)delta(x-x_0)dx
∫x1x2f(x)δ(x−x0)dx=f(x0)χ(x1,x2)(x0)=∫x1x2f(x0)δ(x−x0)dx
与Dirac函数的乘积可以被简化为
f
(
x
)
δ
(
x
−
x
0
)
=
f
(
x
0
)
δ
(
x
−
x
0
)
f(x)delta(x-x_0)=f(x_0)delta(x-x_0)
f(x)δ(x−x0)=f(x0)δ(x−x0)
比如
x
δ
(
x
−
x
0
)
=
x
0
δ
(
x
−
x
0
)
δ
(
x
)
δ
(
x
−
x
0
)
=
δ
(
x
0
)
δ
(
x
−
x
0
)
xdelta(x-x_0)=x_0delta(x-x_0) \ delta(x)delta(x-x_0)=delta(x_0)delta(x-x_0)
xδ(x−x0)=x0δ(x−x0)δ(x)δ(x−x0)=δ(x0)δ(x−x0)
但是还是要谨慎一点,因为 δ ( x 0 ) δ ( x − x 0 ) delta(x_0)delta(x-x_0) δ(x0)δ(x−x0)不管怎么想都有点违和,当 x 0 ≠ 0 x_0ne 0 x0=0时, δ ( x 0 ) δ ( x − x 0 ) = 0 delta(x_0)delta(x-x_0)=0 δ(x0)δ(x−x0)=0,而当 x 0 = 0 x_0=0 x0=0时,Dirac函数是不定型 ∞ ⋅ ∞ infty cdot infty ∞⋅∞,所以这个等式并不成立。
性质2(scaling property)
对任意well-behaved function
f
(
x
)
f(x)
f(x),
∫
x
1
x
2
f
(
x
)
δ
(
x
−
x
0
b
)
d
x
=
u
=
x
/
b
∫
x
1
x
2
f
(
b
u
)
δ
(
u
−
x
0
b
)
d
(
∣
b
∣
u
)
=
∣
b
∣
f
(
b
⋅
x
0
b
)
=
∣
b
∣
f
(
x
0
)
int_{x_1}^{x_2} f(x)delta left( frac{x-x_0}{b}right)dx=_{u=x/b}int_{x_1}^{x_2} f(bu)delta left( u-frac{x_0}{b}right)d(|b|u) \ = |b|fleft( b cdot frac{x_0}{b} right)=|b|f(x_0)
∫x1x2f(x)δ(bx−x0)dx=u=x/b∫x1x2f(bu)δ(u−bx0)d(∣b∣u)=∣b∣f(b⋅bx0)=∣b∣f(x0)
根据性质1,
∫
x
1
x
2
f
(
x
)
δ
(
x
−
x
0
b
)
d
x
=
∫
x
1
x
2
f
(
x
)
[
∣
b
∣
δ
(
x
−
x
0
)
]
d
x
int_{x_1}^{x_2} f(x)delta left( frac{x-x_0}{b}right)dx=int_{x_1}^{x_2} f(x)[|b|delta(x-x_0)]dx
∫x1x2f(x)δ(bx−x0)dx=∫x1x2f(x)[∣b∣δ(x−x0)]dx
所以
δ
(
x
−
x
0
b
)
=
∣
b
∣
δ
(
x
−
x
0
)
delta left( frac{x-x_0}{b}right)=|b|delta(x-x_0)
δ(bx−x0)=∣b∣δ(x−x0)
scaling换一种形式来写也是一样的,比如
a
>
0
a>0
a>0,
δ
(
a
x
−
x
0
)
=
δ
(
x
−
x
0
a
a
)
=
1
∣
a
∣
δ
(
x
−
x
0
a
)
delta(ax-x_0)=deltaleft( frac{x-frac{x_0}{a}}{a} right)=frac{1}{|a|}delta left( x-frac{x_0}{a}right)
δ(ax−x0)=δ(ax−ax0)=∣a∣1δ(x−ax0)
Comb函数
定义Comb函数为
c
o
m
b
(
x
)
=
∑
n
=
−
∞
+
∞
δ
(
x
−
n
)
comb(x) = sum_{n=-infty}^{+infty}delta(x-n)
comb(x)=n=−∞∑+∞δ(x−n)
根据Dirac函数的性质2,
c
o
m
b
(
x
b
)
=
∣
b
∣
∑
n
=
−
∞
+
∞
δ
(
x
−
n
b
)
combleft( frac{x}{b} right)=|b| sum_{n=-infty}^{+infty}delta(x-nb)
comb(bx)=∣b∣n=−∞∑+∞δ(x−nb)
Comb函数一个很重要的作用是采样,假设一个连续波形为
f
(
x
)
f(x)
f(x),计算
f
(
x
)
(
1
∣
b
∣
c
o
m
b
(
x
b
)
)
=
∑
n
=
−
∞
+
∞
f
(
x
)
δ
(
x
−
n
b
)
=
∑
n
=
−
∞
+
∞
f
(
n
b
)
δ
(
x
−
n
b
)
f(x) left( frac{1}{|b|} combleft( frac{x}{b} right)right)=sum_{n=-infty}^{+infty} f(x)delta(x-nb) \ = sum_{n=-infty}^{+infty} f(nb)delta(x-nb)
f(x)(∣b∣1comb(bx))=n=−∞∑+∞f(x)δ(x−nb)=n=−∞∑+∞f(nb)δ(x−nb)
这样就从这个波形中采样得到了一个离散信号 { f ( n b ) } {f(nb)} {f(nb)}。
涉及Dirac函数的微分与积分
先考虑Dirac函数的原函数,
u
(
x
)
=
∫
−
∞
x
δ
(
α
)
d
α
=
χ
(
0
,
+
∞
)
(
x
)
u(x)= int_{-infty}^x delta(alpha)dalpha=chi_{(0,+infty)}(x)
u(x)=∫−∞xδ(α)dα=χ(0,+∞)(x)
这其实就是上一讲定义的step function,阶跃函数,但是阶跃函数在
x
=
0
x=0
x=0处的取值实际上可以取
[
0
,
1
]
[0,1]
[0,1]上的任意值,所以阶跃函数的含义比
u
(
x
)
u(x)
u(x)更宽泛,但是不管阶跃函数在原点处取值是什么,它的导数总是等于Dirac函数的
δ
(
x
)
=
d
d
x
s
t
e
p
(
x
)
=
d
d
x
u
(
x
)
delta(x)=frac{d}{dx}step(x)=frac{d}{dx}u(x)
δ(x)=dxdstep(x)=dxdu(x)
接下来讨论Dirac函数的导数,先把记号叙述一下,
d
d
x
δ
(
x
)
=
δ
′
(
x
)
,
d
n
d
x
n
δ
(
x
)
=
δ
(
n
)
(
x
)
frac{d}{dx}delta(x)=delta'(x),frac{d^n}{dx^n} delta(x)=delta^{(n)}(x)
dxdδ(x)=δ′(x),dxndnδ(x)=δ(n)(x)
它们分别代表Dirac函数的一阶导数和
n
n
n阶导数,虽然可以根据
δ
(
x
)
delta(x)
δ(x)的定义写出
δ
′
(
x
)
delta'(x)
δ′(x)的定义
δ
′
(
x
)
=
lim
b
→
0
d
d
x
(
1
∣
b
∣
s
i
n
c
(
x
b
)
)
delta'(x)=lim_{b to 0} frac{d}{dx} left( frac{1}{|b|} sinc left( frac{x}{b} right) right)
δ′(x)=b→0limdxd(∣b∣1sinc(bx))
但它们有什么意义呢?我们从积分的角度理解一下,假设
f
(
x
)
f(x)
f(x)是一个良定义的有界可微函数,
x
0
<
∞
x_0<infty
x0<∞
∫
−
∞
+
∞
f
(
x
)
δ
′
(
x
−
x
0
)
d
x
=
∫
−
∞
+
∞
f
(
x
)
d
δ
(
x
−
x
0
)
=
f
(
x
)
δ
(
x
−
x
0
)
∣
−
∞
+
∞
−
∫
−
∞
+
∞
f
′
(
x
)
δ
(
x
−
x
0
)
d
x
=
−
∫
−
∞
+
∞
f
′
(
x
)
δ
(
x
−
x
0
)
d
x
=
−
f
′
(
x
0
)
int_{-infty}^{+infty}f(x)delta'(x-x_0)dx=int_{-infty}^{+infty}f(x)ddelta(x-x_0) \ = f(x)delta(x-x_0) |_{-infty}^{+infty}-int_{-infty}^{+infty}f'(x)delta(x-x_0)dx \ = -int_{-infty}^{+infty}f'(x)delta(x-x_0)dx=-f'(x_0)
∫−∞+∞f(x)δ′(x−x0)dx=∫−∞+∞f(x)dδ(x−x0)=f(x)δ(x−x0)∣−∞+∞−∫−∞+∞f′(x)δ(x−x0)dx=−∫−∞+∞f′(x)δ(x−x0)dx=−f′(x0)
可以自行验证
∫
−
∞
+
∞
f
(
x
)
δ
(
n
)
(
x
−
x
0
)
d
x
=
(
−
1
)
n
f
(
n
)
(
x
)
int_{-infty}^{+infty} f(x)delta^{(n)}(x-x_0)dx = (-1)^n f^{(n)}(x)
∫−∞+∞f(x)δ(n)(x−x0)dx=(−1)nf(n)(x)
最后
以上就是隐形宝贝为你收集整理的UA OPTI512R 傅立叶光学导论6 Dirac函数及其性质的全部内容,希望文章能够帮你解决UA OPTI512R 傅立叶光学导论6 Dirac函数及其性质所遇到的程序开发问题。
如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。
发表评论 取消回复