ガンマ関数

提供: testwiki
ナビゲーションに移動 検索に移動
テンプレート:Math のグラフ
テンプレート:Math の絶対値
(グラフ中「Re」は テンプレート:Mvar に相当、「Im」は テンプレート:Mvar に相当)

テンプレート:読み仮名 ruby不使用とは、数学において階乗の概念を複素数全体に拡張した特殊関数。複素階乗とも。一般に Γ(z) と表記される。

自然数 n に対しては、ガンマ関数と n の階乗との間では次の関係式が成り立つ:

n!=Γ(n+1), Γ(n)=(n1)!.

1729年に数学者レオンハルト・オイラーによって無限乗積の形で、最初に導入された[1]Γ という記号は、1814年にルジャンドルが導入した[1]。また、それ以前にガウスが得ており Π などと表記していた(ただし、Π(z)=Γ(z+1) であった)。

定義

実部が正となる複素数 z に対して、次の広域積分で定義される複素関数:

Γ(z)=0tz1etdt(z>0,)

ガンマ関数と呼ぶ[2]。この積分表示は第二種オイラー積分とも呼ばれる。

一般の複素数 z に対しては解析接続もしくは次の極限で定義される。

Γ(z)=limnnzn!k=0n(z+k).

他にも互いに同値となるいくつかの定義が存在する。

基本的性質

0 または負の整数でない、かつ実部が正の任意の複素数 z に対して、

Γ(z+1)=0ettzdt=[ettz]0+z0ettz1dt=zΓ(z)([ettz]0=0),

となることから、Γ(z+1)=zΓ(z) が成り立つ。またさらに、

Γ(1)=0etdt=[et]0=limt(et+1)=1,

である。これらの性質から、任意の正の整数 n に対して、

Γ(n+1)=nΓ(n)=n(n1)Γ(n1)==n!Γ(1)=n!,

より Γ(n+1)=n! が成り立つ。その意味でガンマ関数は階乗の定義域を複素平面に拡張したものとなっている。

歴史的には、ガンマ関数は「階乗の複素数への拡張となるもの」(複素階乗)の実例として、オイラーにより考案された。階乗の複素数への拡張となる関数は無数に存在するが「正の実軸上で対数凸である解析関数」という条件を付ければ、それは一意に定まりガンマ関数に他ならない(ボーア・モレルップの定理)。

右半平面においてオイラー積分で定義されたガンマ関数は全平面に有理型解析接続する。

ガンマ関数は零点を持たず、原点と負の整数に一位の極を持つ。その留数は、

Res(Γ,n)=(1)nn!.

であるテンプレート:Sfn

また、1/2 に対するガンマ関数の値は、ガウス積分の結果に一致する。

Γ(12)=π.

これより、自然数 n に対して、

Γ(12+n)=(2n1)!!2nπ,

が成立することがわかる。ここで !!二重階乗を表す。この性質を利用して高次元の球の体積と表面積を求めることができる。また、

Γ(12n)=(2)n(2n1)!!π.

定義の整合性

定義の積分表示と極限表示が一致することを示す。 テンプレート:Indent とすれば

limn(1tn)n=et

であるから直感的には

limnGn(z)=0tz1etdt

である。(厳密にははさみうちの原理によって証明される)テンプレート:Math の置換により テンプレート:Indent となる.テンプレート:Mvar を除く部分を テンプレート:Math として テンプレート:Indent これにより テンプレート:Indent を得る。故に テンプレート:Indent である。

ワイエルシュトラスの乗積表示

オイラーの乗積表示からオイラーの定数

γ=limn(k=1n1klogn)

を括り出すとワイエルシュトラスの乗積表示が得られる。ワイエルシュトラスはガンマ関数が負の整数にを持つことを嫌って逆数を用いたテンプレート:要出典。ガンマ関数の逆数は複素平面全体で正則であるテンプレート:Sfnテンプレート:Indent

ハンケルの積分表示

ガンマ関数は次の周回積分で表される[3]。積分経路は正の無限大から実軸の上側に沿って原点に至り、原点を正の向きに回り、実軸の下側に沿って無限大に戻るものとする。但し、その偏角はπarg(t)π,0arg(s)2πとする。 テンプレート:Indent これをハンケルの積分表示と呼ぶ。このハンケルの積分表示は、積分経路を適当に変形し、数値積分でガンマ関数の値を求めるために使われることがある[4]

ハンケルの積分表示の導出

極座標表示 (t)=reiθ を用いると、実軸の上側に沿う部分は θ=πr= から r=δ まで、原点を回る部分は r=δθ=π から θ=π まで、実軸の下側に沿う部分は θ=πr=δ から r= までとなる。 テンプレート:Indent z>0 とすると δ0δz0 であるから テンプレート:Indent である。しかし、左辺の被積分関数は z が有界であるかぎり正則であるから、左辺は複素平面全体に解析接続する。従って、 テンプレート:Indent である。 s=reiθ とすれば、同様にして テンプレート:Indent を得る。また、相反公式により、 テンプレート:Indent を得る。

スターリングの公式

テンプレート:Main zでの漸近展開として、ガンマ関数はスターリングの公式で近似される。この漸近近似は複素平面全体(負の実数を除く)で成立するが、|argz|=πに近づくにつれ近似の誤差が大きくなる(極限の収束が遅くなる)ため、応用上は相反公式などを用いて|argz|π/2程度に制限することもある。 テンプレート:Indent

相反公式

次の恒等式相反公式(reflection formula)という[5]テンプレート:Indent 相補公式とも呼ばれるテンプレート:Sfn。 この恒等式はオイラーの乗積表示から得られる。 テンプレート:Indent この分母は正弦関数の無限乗積展開であるから、 テンプレート:Indent である。相反公式にz=12を代入すれば テンプレート:Indent となり テンプレート:Indent を得る。

ルジャンドルの倍数公式

次の恒等式をルジャンドルの倍数公式と呼ぶ。これはガウスの乗法公式の特別な場合である。 テンプレート:Indent

証明

ベータ関数は以下のように表される。

B(z1,z2)=Γ(z1)Γ(z2)Γ(z1+z2)=01tz11(1t)z21dt

ここで z1=z2=z とおくと、

Γ2(z)Γ(2z)=01tz1(1t)z1dt

t=1+x2 とおくと

Γ2(z)Γ(2z)=122z111(1x2)z1dx

(1x2)z1 は偶関数なので

22z1Γ2(z)=2Γ(2z)01(1x2)z1dx

ここで

B(12,z)=01t121(1t)z1dt,t=s2

とすると

B(12,z)=201(1s2)z1ds=201(1x2)z1dx

よって

22z1Γ2(z)=Γ(2z)B(12,z)

よって

B(12,z)=Γ(12)Γ(z)Γ(z+12),Γ(12)=π

よって以下の式が成り立つ。

Γ(z)Γ(z+12)=212zπΓ(2z).

乗法公式

次の恒等式をガウスの乗法公式(multiplication formula)という。 テンプレート:Indent

証明

両辺の比をf(z)とすると テンプレート:Indent 故に、任意に大きな自然数mについてf(z+m)=f(z)が成立する。スターリングの公式により テンプレート:Indent 途中で テンプレート:Indent を適用した。 テンプレート:Indent であり、故に テンプレート:Indent が成立する。

微分方程式

(x, y, y1, , yn)を変数とする多項式F(x, y, y1, , yn)に対し、 テンプレート:Indent の形で表される微分方程式を代数的微分方程式という。ガンマ関数はいかなる代数的微分方程式も満たさないことが知られている[5]。このことはヘルダーが1887年に最初に証明を与えた後 [6]テンプレート:仮リンク[7]テンプレート:仮リンク[8][9]テンプレート:仮リンク[10]ハウスドルフ[11]により、別証明や一般化がなされた。

いくつかの特殊値

テンプレート:See also テンプレート:Indent

ポリガンマ関数

テンプレート:Main ガンマ関数の対数微分

ψ(z)=ddzlogΓ(z)

ディガンマ関数(Digamma function)と呼ぶ。同様の対数微分を繰り返した関数

ψ(n)(z)=dn+1dzn+1logΓ(z)

を、ポリガンマ関数(Polygamma function)と呼ぶ。

脚注

テンプレート:脚注ヘルプ テンプレート:Reflist

参考文献

関連項目

外部リンク

テンプレート:Integral テンプレート:Normdaten

  1. 1.0 1.1 E. T. Whittaker and G. N. Watson (1927), Chapter XII, §12.1
  2. Wolfram mathworld: Gamma Function
  3. Springer Online Reference Works: Gamma-function
  4. Schmelzer & Trefethen (2007), Computing the Gamma function using contour integrals and rational approximations
  5. 5.0 5.1 小松 (2004)、第2章
  6. Otto Ludwig Hölder, "Über die Eigenschaft der Gammafunction keiner algebraischen Differentialgleichung zu genügen," Math. Ann., 28, (1887) pp. 1–13. テンプレート:Doi
  7. Eliakim Hastings Moore, "Concerning transcendentally transcendental functions," Math. Ann., 48 (1897), pp. 49–74. テンプレート:Doi
  8. A. Ostrowski, "Neuer Beweis der Hölderschen Satzes, dass die Gammafunktion keiner algebraischen Differntialgleichung genügt." Math. Ann. 79 (1919), pp. 286–288. テンプレート:Doi
  9. A. Ostrowski, "Zum Hölderschen Satz über Γ(x). Math. Ann. 94 (1925), pp. 248–251. テンプレート:Doi
  10. E. W. Barnes, "The theory of the Gamma function," Messenger of Math. 29 (1900), pp. 64–128.
  11. F. Hausdorff, "Zum Hölderschen Satz über Γ(x)," Math. Ann. 94 (1925), pp. 244–247. テンプレート:Doi