MathLabs

解析学

特殊関数

初等関数を拡張し、数学と物理学のいたるところに現れるガンマ関数、ベータ関数、ゼータ関数、ベッセル関数。

直観階乗を超えて: 整数の間を埋める

階乗 n!=1⋅2⋯nn! = 1\cdot 2 \cdots n は整数に対してしか意味を持たないが、1,2,6,24,120,…1, 2, 6, 24, 120, \dots というように非常になめらかに増加するため、それらの点すべてを通り、実数さらには複素数の入力に対しても定義される曲線を求めたくなる。特殊関数とはまさにこれである: 自然な積分や無限和から生まれた小さな関数の族であり、初等的な算術のパターン(階乗、べき、和)を、確率論、数論、物理学が共有するなめらかで普遍的な仕組みへと補間し一般化するものである。

負の実軸に沿って並ぶ一連の極を示す複素平面のカラープロット。
負の実軸を含む領域のドメインカラーリング図。ガンマ関数はすべての非正整数においてはしご状の単純極を持つ。

大学ガンマ関数、ベータ関数、ゼータ関数

定義: ガンマ関数

Re⁡(z)>0\operatorname{Re}(z) > 0 に対して Γ(z)=∫0∞tz−1e−t dt\Gamma(z)=\displaystyle\int_0^\infty t^{z-1}e^{-t}\,dt と定義する。部分積分により Γ(z+1)=z Γ(z)\Gamma(z+1) = z\,\Gamma(z) が示され、Γ(1)=∫0∞e−t dt=1\Gamma(1) = \int_0^\infty e^{-t}\,dt = 1 であるから、帰納法により任意の非負整数 nn に対して Γ(n+1)=n!\Gamma(n+1)=n! が成り立つ: ガンマ関数は階乗を補間する。この漸化式はまた、Γ\Gamma を z=0,−1,−2,…z = 0, -1, -2, \dots を除く zz 全体へ有理型に拡張し、そこでは単純極を持つ。

Γ(z)=∫0∞tz−1e−t dt\Gamma(z)=\int_0^\infty t^{z-1}e^{-t}\,dt

ベータ関数 B(x,y)=∫01tx−1(1−t)y−1 dtB(x,y) = \displaystyle\int_0^1 t^{x-1}(1-t)^{y-1}\,dt はガンマ関数の比 B(x,y)=Γ(x)Γ(y)Γ(x+y)B(x,y) = \dfrac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)} としてまとめられ、ベータ確率分布の背後にある自然な正規化定数である。より深い恒等式が両者を対称的に結びつける: オイラーの反射公式 Γ(z)Γ(1−z)=πsin⁡(πz)\Gamma(z)\Gamma(1-z) = \dfrac{\pi}{\sin(\pi z)} は、整数でないすべての zz に対して成り立ち、Γ(z)\Gamma(z) の極を sin⁡(πz)\sin(\pi z) の零点と結びつける。

Γ(z)Γ(1−z)=πsin⁡(πz)\Gamma(z)\Gamma(1-z) = \dfrac{\pi}{\sin(\pi z)}

定義: リーマンゼータ関数とその解析接続

Re⁡(s)>1\operatorname{Re}(s) > 1 に対して ζ(s)=∑n=1∞1ns\zeta(s) = \displaystyle\sum_{n=1}^{\infty} \dfrac{1}{n^s} と定義する。この級数はそれ以外の場所では発散するが、ζ(s)\zeta(s) は関数等式 ζ(s)=2sπs−1sin⁡ ⁣(πs2)Γ(1−s) ζ(1−s)\zeta(s) = 2^s \pi^{s-1} \sin\!\left(\dfrac{\pi s}{2}\right) \Gamma(1-s)\, \zeta(1-s) を通じて C\mathbb{C} 全体上の、s=1s=1 にただ一つの単純極を持つ有理型関数へと拡張される。この関数等式は ζ(s)\zeta(s) を ζ(1−s)\zeta(1-s) とガンマ因子で表す。これにより ζ(−1)=−112\zeta(-1) = -\dfrac{1}{12} のような値を厳密に意味づけることができ、これは文字通り 1+2+3+⋯1+2+3+\cdots を足し合わせることとは何の関係もない。

ζ(s)=2sπs−1sin⁡ ⁣(πs2)Γ(1−s) ζ(1−s)\zeta(s) = 2^s \pi^{s-1} \sin\!\left(\dfrac{\pi s}{2}\right) \Gamma(1-s)\, \zeta(1-s)
4つの特殊関数の概観
関数定義式主要な性質
ガンマ, Γ(z)\Gamma(z)Γ(z)=∫0∞tz−1e−t dt\Gamma(z)=\displaystyle\int_0^\infty t^{z-1}e^{-t}\,dtn!n! を補間; 0,−1,−2,…0,-1,-2,\dots に単純極
ベータ, B(x,y)B(x,y)B(x,y)=∫01tx−1(1−t)y−1 dtB(x,y) = \displaystyle\int_0^1 t^{x-1}(1-t)^{y-1}\,dtB(x,y)=Γ(x)Γ(y)Γ(x+y)B(x,y) = \dfrac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)}
ゼータ, ζ(s)\zeta(s)ζ(s)=∑n=1∞1ns\zeta(s) = \displaystyle\sum_{n=1}^{\infty} \dfrac{1}{n^s}有理型接続;s=1s=1 に極が一つ
ベッセル, Jν(x)J_\nu(x)Jν(x)=∑k=0∞(−1)kk! Γ(k+ν+1)(x2)2k+νJ_\nu(x) = \displaystyle\sum_{k=0}^{\infty} \dfrac{(-1)^k}{k!\,\Gamma(k+\nu+1)} \left(\dfrac{x}{2}\right)^{2k+\nu}x2y′′+xy′+(x2−ν2)y=0x^2y''+xy'+(x^2-\nu^2)y=0 の解

発展主要な定理

任意の非負整数 nn に対して Γ(n+1)=n!\Gamma(n+1)=n! が成り立つ。ここで Γ\Gamma は Re⁡(z)>0\operatorname{Re}(z) > 0 に対して Γ(z)=∫0∞tz−1e−t dt\Gamma(z)=\displaystyle\int_0^\infty t^{z-1}e^{-t}\,dt で定義される。

なぜ正しいのか?

これは、ガンマ積分が恣意的な一般化ではなく、階乗の乗法構造を正しく再現する一意な滑らかな延長であることを示しており、だからこそ非整数を代入するのに正しい対象なのである。

証明

手順1(基底段階)。Γ(1)=∫0∞t0e−t dt=∫0∞e−t dt=[−e−t]0∞=0−(−1)=1=0!\Gamma(1) = \int_0^\infty t^0 e^{-t}\,dt = \int_0^\infty e^{-t}\,dt = \big[-e^{-t}\big]_0^\infty = 0-(-1) = 1 = 0! であるから、n=0n=0 で公式が成り立つ。

手順2(部分積分による漸化式)。Re⁡(z)>0\operatorname{Re}(z) > 0 に対して、u=tzu=t^z、dv=e−t dtdv=e^{-t}\,dt として Γ(z)=∫0∞tz−1e−t dt\Gamma(z)=\displaystyle\int_0^\infty t^{z-1}e^{-t}\,dt を部分積分すると Γ(z+1)=∫0∞tze−t dt=[−tze−t]0∞+z∫0∞tz−1e−t dt=0+z Γ(z)\Gamma(z+1) = \int_0^\infty t^z e^{-t}\,dt = \big[-t^z e^{-t}\big]_0^\infty + z\int_0^\infty t^{z-1}e^{-t}\,dt = 0 + z\,\Gamma(z) となる。ここで境界項が両端でゼロになることを用いた。これにより Γ(z+1)=z Γ(z)\Gamma(z+1) = z\,\Gamma(z) が示される。

手順3(帰納法)。ある非負整数 kk について Γ(k+1)=k!\Gamma(k+1) = k! が成り立つと仮定する。手順2で z=k+1z=k+1 とすると、Γ(k+2)=(k+1) Γ(k+1)=(k+1)⋅k!=(k+1)!\Gamma(k+2) = (k+1)\,\Gamma(k+1) = (k+1)\cdot k! = (k+1)! である。

手順4(結論)。基底段階 n=0n=0 が成り立ち(手順1)、帰納段階が n=kn=k から n=k+1n=k+1 へ等式を引き継ぐので(手順3)、数学的帰納法により任意の非負整数 nn に対して Γ(n+1)=n!\Gamma(n+1)=n! が成り立つ。

0<Re⁡(z)<10 < \operatorname{Re}(z) < 1 を満たす任意の zz に対して Γ(z)Γ(1−z)=πsin⁡(πz)\Gamma(z)\Gamma(1-z) = \dfrac{\pi}{\sin(\pi z)} が成り立つ。

なぜ正しいのか?

この公式は、一見純粋に代数的な対象(2つのガンマ値の比)を三角関数と結びつけ、非正整数における Γ\Gamma の極を、すべての整数における sin⁡(πz)\sin(\pi z) の零点と正確に一致させる恒等式である——周回積分による留数計算こそが、右辺の正弦を生み出すものである。

証明

手順1(ベータ関数に帰着する)。y=1−zy=1-z として B(x,y)=Γ(x)Γ(y)Γ(x+y)B(x,y) = \dfrac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)} を用いると、Γ(z)Γ(1−z)=B(z,1−z)=∫01tz−1(1−t)−z dt\Gamma(z)\Gamma(1-z) = B(z,1-z) = \int_0^1 t^{z-1}(1-t)^{-z}\,dt となる。u=t/(1−t)u = t/(1-t) と置換すると、0<Re⁡(z)<10<\operatorname{Re}(z)<1 に対してこれは ∫0∞uz−11+u du\displaystyle\int_0^\infty \dfrac{u^{z-1}}{1+u}\,du になる。

手順2(鍵穴型の積分路を設定する)。I=∫0∞uz−11+u duI=\int_0^\infty \dfrac{u^{z-1}}{1+u}\,du を評価するために、正の実軸に沿って切断された wz−1w^{z-1} の分枝を用い、wz−11+w\dfrac{w^{z-1}}{1+w} を鍵穴型の積分路——切断のすぐ上を外向きに、半径 RR の大円、切断のすぐ下を戻り、00 の周りの半径 ε\varepsilon の小円——に沿って積分する。

手順3(唯一の極における留数)。積分路の内部にある wz−11+w\dfrac{w^{z-1}}{1+w} の唯一の特異点は w=−1w=-1 における単純極であり、そこでは(偏角 π\pi の分枝を用いて)Res=(−1)z−1=eiπ(z−1)\mathrm{Res} = (-1)^{z-1} = e^{i\pi(z-1)} であるから、留数定理により鍵穴積分は 2πi eiπ(z−1)2\pi i\, e^{i\pi(z-1)} となる。

手順4(切断の両側を比較し R→∞R\to\infty、ε→0\varepsilon\to 0 とする)。0<Re⁡(z)<10<\operatorname{Re}(z)<1 のとき大円と小円の寄与は消え、2本の直線部分は wz−1w^{z-1} が切断を横切る際に得る位相だけ異なるので、I(1−e2πi(z−1))=2πi eiπ(z−1)I\big(1 - e^{2\pi i(z-1)}\big) = 2\pi i\, e^{i\pi(z-1)} が得られる。これを II について解き、sin⁡θ=eiθ−e−iθ2i\sin\theta = \frac{e^{i\theta}-e^{-i\theta}}{2i} を用いて指数関数を整理すると I=πsin⁡(πz)I = \dfrac{\pi}{\sin(\pi z)} となり、手順1と合わせて Γ(z)Γ(1−z)=πsin⁡(πz)\Gamma(z)\Gamma(1-z) = \dfrac{\pi}{\sin(\pi z)} が得られる。

発展実世界での応用と具体例

ガンマ関数は待ち時間やベイズ事前分布(ガンマ、ベータ、カイ二乗、スチューデントのt)をモデル化する確率分布を正規化する。ベッセル関数は円柱幾何における波動方程式や熱方程式の自然な解であるため、太鼓の膜の音響学、光学的回折パターン、アンテナの放射パターンに現れる。そしてゼータ関数は解析的整数論(素数の分布)の要であると同時に、カシミール効果や弦理論において発散する和に有限の値を割り当てるゼータ関数正規化を通じて物理学にも現れる。

例: カイ二乗分布の正規化

統計学者が k=5k=5 個の独立な標準正規変数の二乗和を、カイ二乗密度 f(x)=C xk/2−1e−x/2f(x) = C\, x^{k/2-1}e^{-x/2}(x>0x>0)でモデル化する。∫0∞f(x) dx=1\int_0^\infty f(x)\,dx = 1 となるような正規化定数 CC を求めよ。

解答

手順1(ガンマ積分に合わせて置換する)。t=x/2t = x/2 とおくと x=2tx=2t、dx=2 dtdx = 2\,dt であるから、∫0∞xk/2−1e−x/2 dx=∫0∞(2t)k/2−1e−t 2 dt=2k/2∫0∞tk/2−1e−t dt\int_0^\infty x^{k/2-1}e^{-x/2}\,dx = \int_0^\infty (2t)^{k/2-1}e^{-t}\, 2\,dt = 2^{k/2}\int_0^\infty t^{k/2-1}e^{-t}\,dt。

手順2(ガンマ関数を認識する)。残った積分は定義によりちょうど Γ(k/2)\Gamma(k/2) であるから、∫0∞xk/2−1e−x/2 dx=2k/2 Γ(k/2)\int_0^\infty x^{k/2-1}e^{-x/2}\,dx = 2^{k/2}\,\Gamma(k/2) である。

手順3(k=5k=5 を代入する)。ここで k/2=5/2k/2 = 5/2 であり、漸化式 Γ(5/2)=(3/2)(1/2)Γ(1/2)=34π\Gamma(5/2) = (3/2)(1/2)\Gamma(1/2) = \dfrac{3}{4}\sqrt{\pi} を用いると、25/2=422^{5/2}=4\sqrt{2} を簡約して積分は 42⋅34π=32π4\sqrt{2}\cdot\frac{3}{4}\sqrt{\pi}=3\sqrt{2\pi} となる。

手順4(CC を求める)。C⋅32π=1C \cdot 3\sqrt{2\pi} = 1 であるから、正規化定数は C=132πC = \dfrac{1}{3\sqrt{2\pi}} であり、一般のカイ二乗公式 C=12k/2Γ(k/2)C = \dfrac{1}{2^{k/2}\Gamma(k/2)} と一致する。

例: 円形の太鼓の膜の基本振動数

半径 aa の円形の太鼓が、動径方向に対称な振動モード u(r,t)=J0(kr)cos⁡(ωt)u(r,t) = J_0(kr)\cos(\omega t) を持ち、剛体の縁で u(a,t)=0u(a,t)=0(すべての tt について)を満たすとする。J0J_0 の最初の正の零点が j0,1≈2.405j_{0,1}\approx 2.405 であるとき、基本角振動数 ω1\omega_1 を aa と波の速さ cc(ただし ω=ck\omega = ck)を用いて表せ。

解答

手順1(境界条件を適用する)。u(a,t)=J0(ka)cos⁡(ωt)u(a,t)=J_0(ka)\cos(\omega t) がすべての tt についてゼロにならなければならないので、J0(ka)=0J_0(ka) = 0 が必要であり、kaka はベッセル関数 J0J_0 の零点でなければならない。

手順2(基本モードを選ぶ)。基本モードは最小の正の kk に対応し、したがって J0J_0 の最小の正の零点 j0,1≈2.405j_{0,1}\approx 2.405 に対応する。よって k1a=j0,1k_1 a = j_{0,1} であるから k1=j0,1/ak_1 = j_{0,1}/a である。

手順3(角振動数に変換する)。膜の波動方程式の分散関係 ω=ck\omega = ck を用いると、ω1=ck1=c j0,1a\omega_1 = c k_1 = \dfrac{c\, j_{0,1}}{a} となる。

手順4(解釈する)。数値を代入すると ω1≈2.405 ca\omega_1 \approx \dfrac{2.405\,c}{a} となる: 半径 aa を2倍にすると基本振動数は半分になり、これは大きな太鼓がより低い音を出すという物理的な直感とまさに一致する。

Γ(6)\Gamma(6) はいくらか。

Γ ⁣(12)\Gamma\!\left(\dfrac{1}{2}\right) はいくらか。

ζ(s)\zeta(s) の有理型接続はただ一つの極をどこに持ち、それはどのような種類か。

ある音響エンジニアが円形の太鼓の膜を設計し、中心で有限にとどまり剛体の縁でゼロになるモード形状を必要としている。これらの動径方向のモード形状を与える特殊関数の族はどれか。

参考文献

  1. Frank W. J. Olver, Ronald F. Boisvert, Daniel W. Lozier, Charles W. Clark (eds.) (2023). NIST Digital Library of Mathematical Functions, Chapter 5: Gamma Function
  2. Frank W. J. Olver, Ronald F. Boisvert, Daniel W. Lozier, Charles W. Clark (eds.) (2023). NIST Digital Library of Mathematical Functions, Chapter 25: Zeta and Related Functions
  3. David J. Platt, Timothy S. Trudgian (2021). The Riemann Hypothesis Is True Up to 3×10^12 · arXiv:2004.09765