Excel VBA 数学教室ではアフィリエイトプログラムを利用して商品を紹介しています。

超関数とテスト関数・ディラックのデルタ関数

このページでは冒頭で超関数について簡単に触れ、その中のひとつであるディラックのデルタ関数について解説します。 

超関数とテスト関数

超関数の性質を調べるときには、テスト関数 ϕ(x) というものを多用します。少し難しい表現になってしまいますが、テスト関数とは「台がコンパクトであり、何回でも微分可能な連続関数」のことです。「台」とは関数値が 0 でない定義域のことです。「コンパクト」は厳密に考えると難しいのですが、とりあえず「有界閉区間の集合である」と理解しておけば十分です。つまり「台がコンパクトである」というのは、ϕ(±)=0 であることを意味しています。大雑把に次のような関数をイメージしてください。
 
コンパクトな関数の数学グラフ
このようなテスト関数 ϕ(x) に対して
 (A)ϕ(x)f(x)dx
が有限の値をとるような関数 f(x)超関数と定義します。超関数は普通の関数を拡張定義したものなので、二次関数や指数関数のような、ごく普通の関数も(上の積分は有限値をとるので)やはり超関数に含まれます。ただし、後述するデルタ関数のように、普通の関数とは明らかに異なる性質をもつ関数もその範疇に含まれることになります。そして
 (B)ϕ(x)f(x)dx=ϕ(x)g(x)dx
が成り立つことを「超関数 f(x)g(x) が等しい」と定義することにします。つまり f(x)g(x) の各点での値は必ずしも一致しません(つまりグラフが重なりません)が、全区間で積分した値が同じであれば、それを f(x)=g(x) と書くことにします。

ディラックのデルタ関数

デルタ関数はイギリスの物理学者 P.A.M.Dirac によって考案された関数です。図のように x 軸上に原点を中心として 2m の幅の中に物質が密度 1/2m で分布しているような状況を考えます。
 
Excel分布関数のグラフ

この分布関数 p(x,m)
 p(x,m)={12m(mxm)0(x<m,m<x)
のように表すことができます。p(x,m)x の全区間で積分すると
 p(x,m)dx=mm12mdx=1
となります(つまり総質量が 1 です)。この分布関数において m をどんどん小さくしていくと下の図のように原点付近に圧縮された密度の高い状態となります。
 
Excel分布関数のグラフ②

もちろん先ほどの積分計算の結果は同じですから総質量は 1 に保たれています。そして m0 という極限を考えると
 limm0p(x,m)={(x=0)0(x0)
という関数が現れます。このような特異な関数をデルタ関数とよび、δ(x) と書きます。デルタ関数を正確に図に表すことはできませんが、限りなく細長くしたような棒線をイメージしてください。ここまでをまとめると、デルタ関数は
 (C)δ(x)={(x=0)0(x0)(D)δ(x)dx=1
といった性質をもつ関数です。ここで分布関数にテスト関数 ϕ(x) をかけて全区間で積分してみます。
 ϕ(x)p(x,m)dx=mmϕ(x)2mdx
ここで積分の平均値の定理(面積を適当な高さの長方形に均すことができるという定理)によって、
 mmϕ(x)2mdx=ϕ(c)2m2m=ϕ(c)
となるような ϕ(c) が存在しています。m0 の極限をとれば、明らかに ϕ(c)ϕ(0) となります。すなわち
 (E)ϕ(x)δ(x)dx=ϕ(0)
となります。ここで改めて、デルタ関数は上式 (E) を満たすような関数であると定義します。

デルタ関数の基本公式についてまとめておきます。
 (F)f(x)δ(x)=f(0)δ(x)(G)xδ(x)=0(H)δ(ax)=δ(x)a(I)δ(x)=δ(x)
【基本公式[F]と[G]の証明】先に述べたように、超関数 f(x)g(x) が等しいというのは、
 ϕ(x)f(x)dx=ϕ(x)g(x)dx
が成り立つということです。そこで左辺 f(x)δ(x)ϕ(x) をかけて積分すると
 {f(x)δ(x)}ϕ(x)dx=δ(x){f(x)ϕ(x)}dx
ここでデルタ関数の基本的な性質である
 δ(x)ϕ(x)dx=ϕ(0)
を用いると
 {f(x)δ(x)}ϕ(x)dx=f(0)ϕ(0)
となります。また右辺 f(0)δ(x) にテスト関数 ϕ(x) をかけて積分すると
 f(0)δ(x)ϕ(x)dx=f(0)ϕ(0)
となります。よって
 f(x)δ(x)=f(0)δ(x)
が証明されました。この公式で f(x)=x とおけば
 xδ(x)=0
は明らかです。f(x)=x に限らず、f(x)=x2 など、原点で 0 になる関数を δ(x) に掛ければ 0 になります。そうでなければ公式 [F] より y 切片の値に δ(x) を乗じた関数が得られます。グラフを思い描いてみると、公式というより、ほとんど自明の事実です。
【基本公式[H]と[I]の証明】δ(ax) にテスト関数 ϕ(x) をかけて積分します。a>0 のとき、ax=t とおくと、
 ϕ(x)δ(ax)dx=1aϕ(ta)δ(t)dt=ϕ(0)a=ϕ(x)δ(x)adx
a<0 のとき、a=b(b>0),ax=bx=t とおくと、
 ϕ(x)δ(ax)dx=1bϕ(tb)δ(t)dt=ϕ(0)b=ϕ(0)a=ϕ(x)δ(x)adx
よって被積分関数同士を比較すると
 δ(ax)=δ(x)a
が得られます。この公式で a=1 とすると
 δ(x)=δ(x)
すなわちデルタ関数が偶関数であることがわかります。

エクセルや数学に関するコメントをお寄せください