我是靠谱客的博主 隐形宝贝,最近开发中收集的这篇文章主要介绍UA OPTI512R 傅立叶光学导论6 Dirac函数及其性质,觉得挺不错的,现在分享给大家,希望可以做个参考。

概述

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)=b0limb1sinc(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)δ(xx0)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)δ(xx0)dx=x1x2δ(xx0)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)δ(xx0)dx=f(x0)χ(x1,x2)(x0)=x1x2f(x0)δ(xx0)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)δ(xx0)=f(x0)δ(xx0)

比如
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δ(xx0)=x0δ(xx0)δ(x)δ(xx0)=δ(x0)δ(xx0)

但是还是要谨慎一点,因为 δ ( x 0 ) δ ( x − x 0 ) delta(x_0)delta(x-x_0) δ(x0)δ(xx0)不管怎么想都有点违和,当 x 0 ≠ 0 x_0ne 0 x0=0时, δ ( x 0 ) δ ( x − x 0 ) = 0 delta(x_0)delta(x-x_0)=0 δ(x0)δ(xx0)=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)δ(bxx0)dx=u=x/bx1x2f(bu)δ(ubx0)d(bu)=bf(bbx0)=bf(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)δ(bxx0)dx=x1x2f(x)[bδ(xx0)]dx

所以
δ ( x − x 0 b ) = ∣ b ∣ δ ( x − x 0 ) delta left( frac{x-x_0}{b}right)=|b|delta(x-x_0) δ(bxx0)=bδ(xx0)

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) δ(axx0)=δ(axax0)=a1δ(xax0)

Comb函数

定义Comb函数为
c o m b ( x ) = ∑ n = − ∞ + ∞ δ ( x − n ) comb(x) = sum_{n=-infty}^{+infty}delta(x-n) comb(x)=n=+δ(xn)

根据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)=bn=+δ(xnb)

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)(b1comb(bx))=n=+f(x)δ(xnb)=n=+f(nb)δ(xnb)

这样就从这个波形中采样得到了一个离散信号 { 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)=b0limdxd(b1sinc(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)δ(xx0)dx=+f(x)dδ(xx0)=f(x)δ(xx0)++f(x)δ(xx0)dx=+f(x)δ(xx0)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)(xx0)dx=(1)nf(n)(x)

最后

以上就是隐形宝贝为你收集整理的UA OPTI512R 傅立叶光学导论6 Dirac函数及其性质的全部内容,希望文章能够帮你解决UA OPTI512R 傅立叶光学导论6 Dirac函数及其性质所遇到的程序开发问题。

如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。

本图文内容来源于网友提供,作为学习参考使用,或来自网络收集整理,版权属于原作者所有。
点赞(56)

评论列表共有 0 条评论

立即
投稿
返回
顶部