9

Γ, B関数のお話

525
1

はじめに

入試問題の背景にあることも時々ある, Γ(ガンマ)関数, B(ベータ)関数について, 簡単な説明をしていきます.

目次

  • Γ関数の定義と性質
  • B関数の定義と性質
  • Γ数とB関数の関係のお話
  • その応用
  • おわりに

Γ関数の定義と性質

正の実数xに対し, Γ関数Γ(x)を以下で定義します.

Γ(x)=0tx1etdt

ただし, 0f(x)dx=limR0Rf(x)dx という意味です.

このとき, 部分積分より,
Γ(x+1)=0txetdt=[txet]00xtx1etdt=x0tx1etdt=xΓ(x)

従って,

Γ(x+1)=xΓ(x)

が成り立ちます.

さらに, Γ(1)=0etdt=limR[et]0R=1なので, 自然数nに対し,
Γ(n+1)=nΓ(n)=n(n1)Γ(n1)==n(n1)21Γ(1)=n!

従って,

Γ(n+1)=n!

が成り立ちます.

以上より, Γ関数は, 階乗の概念の一般化と言うことができます.

階乗の一般化はいくつもあり得るのでは, と思う方もいるかもしれませんが, 実は, 自然数のときに階乗と一致してかつ良い性質(複素微分可能性と対数凸性などです)を持つような一般化は, このΓ関数しかないことが知られています.

B関数の定義と性質

正の実数x,yに対し, B関数B(x,y)を以下で定義します.

B(x,y)=01tx1(1t)y1dt

まず, 実数α,β(α<β)に対してs=α+t(βα)=(1t)α+tβと置換すると, t=sαβα, 1t=βsβα なので,

01tx1(1t)y1dt=αβ(sαβα)x1(βsβα)y1dsβα=1(βα)x+y1αβ(sα)x1(βs)y1ds

従って,

αβ(sα)x1(βs)y1ds=(βα)x+y1B(x,y)

となります.

これは俗に言う1/6公式などの一般化になっていることがわかると思います. (B関数の値の求め方は次節で紹介します.)

次に, 定義式において t=sin2θ と置換すると,
01tx1(1t)y1dt=0π2(sin2θ)x1(cos2θ)y12sinθcosθdθ=20π2sin2x1θcos2y1θdθ

従って,

0π2sin2x1θcos2y1θdθ=12B(x,y)

となります.

Γ関数とB関数の関係のお話

この節では, 綺麗な関係式

B(x,y)=Γ(x)Γ(y)Γ(x+y)

の証明をしていきます.

これには重積分が必要となるので, 多少, 厳密さを欠いた証明となってしまいますが許してください.

定義式より,
Γ(x+y)B(x,y)=0tx+y1etdt20π2sin2x1θcos2y1θdθ=20r2(x+y)1er2dr20π2sin2x1θcos2y1θdθ   (t=r2)=4r=0r=θ=0θ=π2(rcosθ)2x1(rsinθ)2y1er2rdθdr

ここで, u=rcosθ, v=rsinθ とおくと, r,θの動きに合わせて, (u,v)は, uv平面の第1象限, 即ち u>0,v>0 を動きます.

いま, この重積分の意味を考えてみると, rθ平面上の微小な長方形drdθ (2変数なので微小な区間ではなく微小な長方形となります. ) に対して, 関数の値と微小面積の積を足し合わせたものと捉えられます.

この長方形が u=rcosθ, v=rsinθ なる変数変換によってuv平面でどのような図形になるかを考えると, これは各辺がdr,rdθの長方形となります. (具体的には, r,θdr,dθ増えたときの点(u,v)の移動を考えると, ベクトル(cosθsinθ)dr, (rsinθrcosθ)dθ によって囲まれた部分を動くとわかります. )

従って, 微小面積はuv平面ではrdrdθとなることがわかりました. つまり rdrdθ=dudv となります.

これを用いて変数変換を行うと, u,vの積分区間はどちらも0であることにも注意して,

Γ(x+y)B(x,y)=4r=0r=θ=0θ=π2(rcosθ)2x1(rsinθ)2y1er2rdθdr=400u2x1v2y1e(u2+v2)dudv

となります. 後はこれを変形していけば,
Γ(x+y)B(x,y)=400u2x1v2y1e(u2+v2)dudv=20u2x1eu2du20v2y1ev2dv=0sx1esds0ty1etdt   (u2=s,v2=t)=Γ(x)Γ(y)

以上より, 示すことができました. お疲れ様でした.

その応用

では, Γ関数とB関数の関係式

B(x,y)=Γ(x)Γ(y)Γ(x+y)

を用いて, いろいろな計算をしていきましょう.

まずは, 1/6公式の一般化として紹介した式

αβ(sα)x1(βs)y1ds=(βα)x+y1B(x,y)

の右辺を具体的に計算してみます. 簡単のため自然数m,nを用いてx=m+1, y=n+1 とおくと,
αβ(sα)m(βs)nds=(βα)m+n+1B(m+1,n+1)=(βα)m+n+1Γ(m+1)Γ(n+1)Γ(m+n+2)=m!n!(m+n+1)!(βα)m+n+1

と計算できます. 見やすいように積分変数をxに変えて,

αβ(xα)m(βx)ndx=m!n!(m+n+1)!(βα)m+n+1

を得ます. 実際にm=n=1とすると1/6公式になっていることが分かると思います.

次に, 例の関係式で x=y=12 としてみます. すると,
B(12,12)=Γ(12)Γ(12)Γ(1)=Γ(12)2

一方, B関数の, 三角関数を用いた積分表示から

B(12,12)=20π2dθ=π

であるので, Γ(x)>0 と併せて以下を得ます.

Γ(12)=π

これは予想外の値になっていて面白いと思います.

最後に, B関数の, 三角関数の積分による表示

0π2sin2x1θcos2y1θdθ=12B(x,y)

において x=n+12, y=12 としてみます. すると,

0π2sin2nθdθ=12B(n+12,12)=12Γ(n+12)Γ(12)Γ(n+1)=n12nΓ(n12)Γ(n)π2==n12nn32n1121Γ(12)Γ(1)π2=2n12n2n32n212π2

のようになります. 奇数乗の場合も同様に求められます.

おわりに

ここまで読んでくださった方, 本当にありがとうございます.

簡単な紹介と言ったのに, いろいろ書いて長くなってしまってすみません.

この先の発展としては, Γ関数の無限積表示, 対数微分, 相反公式, 等があると思います. 興味がある方はお勉強してみると, とても面白いですよ.

では改めて, ありがとうございました.

投稿日:2020119
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

東大数理M1

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. はじめに
  2. 目次
  3. Γ関数の定義と性質
  4. B関数の定義と性質
  5. Γ関数とB関数の関係のお話
  6. その応用
  7. おわりに