ディラックのデルタ関数





Dirac distribution PDF.png

数学におけるディラックのデルタ関数(デルタかんすう、英: delta function)、制御工学におけるインパルス関数 (インパルスかんすう、英: impulse function) とは、任意の実連続関数 f : RR に対し、


f(x)δ(x)dx=f(0){displaystyle int _{-infty }^{infty }f(x)delta (x),dx=f(0)}int _{{-infty }}^{{infty }}f(x)delta (x),dx=f(0)

を満たす実数値シュワルツ超関数 δ のことである。これはクロネッカーのデルタ


i=−+∞fiδi,j=fj{displaystyle sum _{i=-infty }^{+infty }f_{i}delta _{i,j}=f_{j}}sum _{{i=-infty }}^{{+infty }}f_{i}delta _{{i,j}}=f_{j}

の自然な拡張になっている。


ディラックのデルタ関数は、デルタ超関数 (delta distribution) あるいは単にディラックデルタ (Dirac's delta) とも呼ばれる。これを最初に定義して量子力学の定式化に用いた物理学者ポール・ディラックに因み、この名称が付いている。デルタ関数は古典的な意味での関数ではないシュワルツ超関数 (distribution) の最初の例になっている。


ディラックのデルタの「関数」としての性質は、形式的に次のように述べることができる。まず、f(x) として実直線上常に一定の値 1 をとる関数をとり、デルタ関数をデルタ関数自身と f(x) = 1 との積であると見ることにより


δ(x)dx=1{displaystyle int _{-infty }^{infty }delta (x)dx=1}int _{{-infty }}^{{infty }}delta (x)dx=1

である。一方、積分値が fx = 0 での値にしかよらないことから


δ(x)=0(x≠0){displaystyle delta (x)=0quad (xneq 0)}delta (x)=0quad (xneq 0)

でなければならないが、その上で積分値が 0 でない有限の値をとるためには


δ(0)=∞{displaystyle delta (0)=infty }delta (0)=infty

が満たされなければならない。




目次






  • 1 概要


  • 2 初等関数による近似


    • 2.1 正規分布の密度関数による近似


    • 2.2 Sinc関数による近似




  • 3 佐藤超関数としての定義


  • 4 ディラック測度


  • 5 演算子体の単位元


  • 6 関連項目


  • 7 参考文献





概要


上記のように特徴付けられるデルタ関数 δ(x) は、その名前にも現れているように、あたかも通常の関数であるかのように扱われることも珍しくないが、実際には通常の意味の関数と見なすことはできない。
例えば、デルタ関数を連続関数で表すことができないことは以下のようにして分かる。δ(x) が連続関数だったとして x = 0 でゼロでない値をとるならば x = 0 を含む小区間で非ゼロでなければならず、x ≠ 0δ(x) = 0 という条件を満たせない。したがって x ≠ 0δ(x) = 0 ならばそれは常に 0 の値をとる関数であり、他の関数と掛けて積分しても 0 以外の値をとることはない。点 x = 0 においてのみ不連続であることを認めても、デルタ関数の特徴付けに用いられている積分が、通常の関数の(広義)リーマン積分やルベーグ積分として理解されるならば、このような関数の積分は恒等的に 0 に等しい関数を積分するのと同じであり積分値は 0 になる。したがって、このような条件を満たすような通常の関数は存在しない。


しかし、通常の意味ではまったく関数ではないデルタ関数は、適当な枠組みの下では意味を持ち、例えばデルタ分布はヘヴィサイドの階段関数の弱微分(超関数の意味での微分)を与えている。



初等関数による近似


デルタ関数は実軸上滑らかで有界な関数の空間 Cb∞(R){displaystyle C_{b}^{infty }(mathbb {R} )}C_{b}^{infty }(mathbb{R} ) 上の汎関数になっているが、Cb∞(R){displaystyle C_{b}^{infty }(mathbb {R} )}C_{b}^{infty }(mathbb{R} ) の双対空間の中でデルタ関数に弱収束するような関数の族 φt、つまり



f(x)ϕt(x)dx→f(0)(t→0){displaystyle int _{-infty }^{infty }f(x)phi _{t}(x)dxrightarrow f(0)quad (trightarrow 0)}int _{{-infty }}^{infty }f(x)phi _{t}(x)dxrightarrow f(0)quad (trightarrow 0)



が任意の f∈Cb∞(R){displaystyle fin C_{b}^{infty }(mathbb {R} )}fin C_{b}^{infty }(mathbb{R} ) について成り立つような族 φt がいくつか知られている。同様にして、滑らかかつ有界とは別な条件を満たす関数の空間の上の汎関数としての弱収束の表示も与えられている。以下に代表的例を 2 つ挙げる。



正規分布の密度関数による近似


中心 μ, 分散 σ2 の正規分布の密度関数



(x)=12πσexp{−(x−μ)22σ2}{displaystyle f_{mu ,sigma }(x)={frac {1}{{sqrt {2pi }},sigma }}exp !left{-{frac {(x-mu )^{2}}{2sigma ^{2}}}right}}f_{{mu ,sigma }}(x)={frac  {1}{{sqrt  {2pi }},sigma }}exp !left{-{frac  {(x-mu )^{2}}{2sigma ^{2}}}right}



は、デルタ関数の満たすべき性質



(x)dx=1{displaystyle int _{-infty }^{infty }f_{mu ,sigma }(x),dx=1}int _{{-infty }}^{{infty }}f_{{mu ,sigma }}(x),dx=1



を満たす。さらに、μ = 0σ → 0 とすれば x = 0 の近傍の外で一様に fσ(x) → 0 (x ≠ 0) かつ fσ(0) → +∞ である。これは、σ → 0 とすることで、関数族 fσ が汎関数としてデルタ関数に近づくことを意味する。したがって、デルタ関数はある意味で正規分布の密度関数の極限と見なすことができ、



limσ012πσexp(−x22σ2)=δ(x){displaystyle lim _{sigma to 0}{frac {1}{{sqrt {2pi }},sigma }}exp !left(-{frac {x^{2}}{2sigma ^{2}}}right)=delta (x)}lim _{{sigma to 0}}{frac  {1}{{sqrt  {2pi }},sigma }}exp !left(-{frac  {x^{2}}{2sigma ^{2}}}right)=delta (x)



と表現される。デルタ関数の表現に正規分布を用いたが、このことから、デルタ関数は正規分布の一種であると考えることが可能である。デルタ関数は、特殊な確率分布の表現に有用である。



Sinc関数による近似


Sinc関数から変数変換とスケーリングによって得られる関数族



fk(x)=sin⁡kxπx(k∈R){displaystyle f_{k}(x)={frac {sin kx}{pi x}}quad (kin mathbb {R} )}f_{k}(x)={frac  {sin kx}{pi x}}quad (kin {mathbb  {R}})



は、デルタ関数の満たすべき条件



fk(x)dx=1{displaystyle int _{-infty }^{infty }f_{k}(x),dx=1}int _{{-infty }}^{{infty }}f_{k}(x),dx=1



を満たす。ただし、これは左辺を広義積分 lima→aa{displaystyle lim _{arightarrow infty }int _{-a}^{a}}lim _{{arightarrow infty }}int _{{-a}}^{a} として解釈した際に成立する等式である。上記の例と違ってこの関数族は k → ∞ としても各点収束しないが、任意のコンパクト台の滑らかな関数 g に対して



limk→fk(x)g(x)dx=g(0){displaystyle lim _{kto infty }int _{-infty }^{infty }f_{k}(x)g(x)dx=g(0)}lim _{{kto infty }}int _{{-infty }}^{{infty }}f_{k}(x)g(x)dx=g(0)



が成り立っている。これも弱収束の意味でデルタ関数を近似していると考えられ、



limk→sin⁡kxπx=δ(x){displaystyle lim _{kto infty }{frac {sin kx}{pi x}}=delta (x)}lim _{{kto infty }}{frac  {sin kx}{pi x}}=delta (x)



と表現される。



佐藤超関数としての定義


佐藤超関数の流儀では、ディラックのデルタ関数は複素領域から実軸への抽象的境界値



δ(x):=−12πi(1x+i0−1x−i0){displaystyle delta (x):={frac {-1}{2pi i}}left({frac {1}{x+i0}}-{frac {1}{x-i0}}right)}delta (x):={frac  {-1}{2pi i}}left({frac  {1}{x+i0}}-{frac  {1}{x-i0}}right)



と定義される。ここで抽象的境界値とは正則関数のある種の同値類を表すが、直感的には x ≠ 0 ならば



12πi(1x+i0−1x−i0)=−12πi(1x−1x)=0{displaystyle {frac {-1}{2pi i}}left({frac {1}{x+i0}}-{frac {1}{x-i0}}right)={frac {-1}{2pi i}}left({frac {1}{x}}-{frac {1}{x}}right)=0}{frac  {-1}{2pi i}}left({frac  {1}{x+i0}}-{frac  {1}{x-i0}}right)={frac  {-1}{2pi i}}left({frac  {1}{x}}-{frac  {1}{x}}right)=0



である。また、デルタ関数の最も重要な性質である



δ(x)f(x)dx=f(0){displaystyle int delta (x)f(x),dx=f(0)}int delta (x)f(x),dx=f(0)



は、複素解析学のコーシーの積分公式から導かれる。厳密な定義には層係数のコホモロジー論を必要とするが、1 変数の場合は比較的容易に理論展開できる。



ディラック測度


ディラック関数は以下のようにして定まるディラック測度 (Dirac measure) δ0 の非形式的な密度関数だと解釈することができる。実直線のボレル部分集合 A に対して、A0 を含む場合 δ0(A) = 1、そうでない場合 δ0(A) = 0 とすると、δ0σ-加法性を持っている。この測度に関する有界ボレル関数の積分は


f(x)dδ0(x)=f(0){displaystyle int f(x)ddelta _{0}(x)=f(0)}int f(x)ddelta _{0}(x)=f(0)

であり、形式的に 0(x) = δ(x)dx が成り立っている。



演算子体の単位元


ミクシンスキーの演算子法に従い、R≥0 = [0, ∞) 上の複素数値連続関数の全体 C = C([0, ∞); C) が畳み込み


(f∗g)(x):=∫0xf(x−ξ)g(ξ)dξ{displaystyle (f*g)(x):=int _{0}^{x}f(x-xi )g(xi ),dxi }(f*g)(x):=int _{0}^{x}f(x-xi )g(xi ),dxi

に関して零因子を持たないというティッチマーシュの定理(英語版)を用いて、(単位元を持たない可換な)整域としての C の商体 M を構築する Mティッチマーシュ・ミクシンスキー代数や、ミクシンスキー演算子(ヘヴィサイド演算子、—超関数)の体などと呼ばれる。M には C にはなかった乗法の単位元 δ = {δ(x)} が付加されているが、この δ(x) はしばしばデルタ関数と看做される。


実際 δ は、特に定数関数 1 に対応する積分作用素 l = {1} ∈ CM に対して = δl = l すなわち、形式上は任意の x に対して


0xδ)dξ=1{displaystyle int _{0}^{x}delta (xi ),dxi =1}int _{0}^{x}delta (xi ),dxi =1

を満たさなければならない(もし δR≥0 上の連続関数ならば、x = 0 とすれば左辺は 0 となるから、これを C の中だけで考えることはできない)。再び形式的な議論だが、この被積分関数を δ(x)[0, x] の指示関数との値ごとの積と見なすことで、無限区間でのデルタ関数の性質が満たされると考えることができる。一方で、十分小さな ε > 0 に対し


ε)dξ=∫0xδ)dξδ)dξ=0{displaystyle int _{varepsilon }^{x}delta (xi ),dxi =int _{0}^{x}delta (xi ),dxi -int _{0}^{varepsilon }delta (xi ),dxi =0}int _{varepsilon }^{x}delta (xi ),dxi =int _{0}^{x}delta (xi ),dxi -int _{0}^{varepsilon }delta (xi ),dxi =0

だから、x ≠ 0δ(x) = 0 が満たされていると考えることができる。



関連項目



  • ヘヴィサイドの階段関数

  • 超関数

  • インパルス応答

  • 再生核

  • クロネッカーのデルタ



参考文献



  • Weisstein, Eric W. "Delta Function". MathWorld(英語)..mw-parser-output cite.citation{font-style:inherit}.mw-parser-output .citation q{quotes:"""""""'""'"}.mw-parser-output .citation .cs1-lock-free a{background:url("//upload.wikimedia.org/wikipedia/commons/thumb/6/65/Lock-green.svg/9px-Lock-green.svg.png")no-repeat;background-position:right .1em center}.mw-parser-output .citation .cs1-lock-limited a,.mw-parser-output .citation .cs1-lock-registration a{background:url("//upload.wikimedia.org/wikipedia/commons/thumb/d/d6/Lock-gray-alt-2.svg/9px-Lock-gray-alt-2.svg.png")no-repeat;background-position:right .1em center}.mw-parser-output .citation .cs1-lock-subscription a{background:url("//upload.wikimedia.org/wikipedia/commons/thumb/a/aa/Lock-red-alt-2.svg/9px-Lock-red-alt-2.svg.png")no-repeat;background-position:right .1em center}.mw-parser-output .cs1-subscription,.mw-parser-output .cs1-registration{color:#555}.mw-parser-output .cs1-subscription span,.mw-parser-output .cs1-registration span{border-bottom:1px dotted;cursor:help}.mw-parser-output .cs1-ws-icon a{background:url("//upload.wikimedia.org/wikipedia/commons/thumb/4/4c/Wikisource-logo.svg/12px-Wikisource-logo.svg.png")no-repeat;background-position:right .1em center}.mw-parser-output code.cs1-code{color:inherit;background:inherit;border:inherit;padding:inherit}.mw-parser-output .cs1-hidden-error{display:none;font-size:100%}.mw-parser-output .cs1-visible-error{font-size:100%}.mw-parser-output .cs1-maint{display:none;color:#33aa33;margin-left:0.3em}.mw-parser-output .cs1-subscription,.mw-parser-output .cs1-registration,.mw-parser-output .cs1-format{font-size:95%}.mw-parser-output .cs1-kern-left,.mw-parser-output .cs1-kern-wl-left{padding-left:0.2em}.mw-parser-output .cs1-kern-right,.mw-parser-output .cs1-kern-wl-right{padding-right:0.2em}


  • Francois Treves (2006). Topological Vector Spaces, Distributions and Kernels. Dover Publications. 




Popular posts from this blog

How to reconfigure Docker Trusted Registry 2.x.x to use CEPH FS mount instead of NFS and other traditional...

is 'sed' thread safe

How to make a Squid Proxy server?