ディガンマ関数

数学において、ディガンマ関数（でぃがんまかんすう、英: digamma function）とはガンマ関数の対数微分で定義される特殊関数。ポリガンマ関数の一種である。 ガンマ関数 Γ(z) に対し、その対数微分 ψ ( z ) = d d z ln ⁡ Γ ( z ) = Γ ′ ( z ) Γ ( z ) {\displaystyle \psi (z)={\frac {d}{dz}}\ln {\Gamma (z)}={\frac {\Gamma '(z)}{\Gamma (z)}}} をディガンマ関数と呼ぶ。 ディガンマ関数はz = 0, −1, −2, ... で1位の極をもち，それらの点を除く全複素平面では解析的になる。 ガンマ関数のワイエルシュトラスの無限乗積表示 1 Γ ( z ) = lim n → ∞ z ( z + 1 ) ⋯ ( z + n ) n z n ! {\displaystyle {\frac {1}{\Gamma (z)}}=\lim _{n\to \infty }{\frac {z(z+1)\cdots (z+n)}{n^{z}n!}}} を対数微分することで、ディガンマ関数における ψ ( z ) = lim n → ∞ { ln ⁡ n − 1 z − ∑ k = 1 n 1 z + k } {\displaystyle \psi (z)=\lim _{n\to \infty }\left\{\ln {n}-{\frac {1}{z}}-\sum _{k=1}^{n}{\frac {1}{z+k}}\right\}} という表示を得る。特にz　=1とすれば、特殊値 ψ ( 1 ) = lim n → ∞ { ln ⁡ n − ∑ k = 1 n 1 k } = − γ {\displaystyle \psi (1)=\lim _{n\to \infty }\left\{\ln {n}-\sum _{k=1}^{n}{\frac {1}{k}}\right\}=-\gamma } を得る。但し、γ =0.5772...はオイラーの定数である。 また、ディガンマ関数は次の漸化式を満たす。 ψ ( z + 1 ) = ψ ( z ) + 1 z {\displaystyle \psi (z+1)=\psi (z)+{\frac {1}{z}}} この関係式から、一般に ψ ( z + n ) = ψ ( z ) + ∑ k = 1 n 1 z + k − 1 {\displaystyle \psi (z+n)=\psi (z)+\sum _{k=1}^{n}{\frac {1}{z+k-1}}} であり、特に z =1とすれば、特殊値 ψ ( n + 1 ) = − γ + ∑ k = 1 n 1 k {\displaystyle \psi (n+1)=-\gamma +\sum _{k=1}^{n}{\frac {1}{k}}} が得られる。 ディガンマ関数とその導関数はz ≠0, -1, -2, -3...で次の級数表示を持つ。 ψ ( z ) = − γ − ∑ n = 0 ∞ ( 1 z + n − 1 n + 1 ) = − γ + ∑ n = 0 ∞ z − 1 ( n + 1 ) ( z + n ) {\displaystyle \psi (z)=-\gamma -\sum _{n=0}^{\infty }{\biggl (}{\frac {1}{z+n}}-{\frac {1}{n+1}}{\biggr )}=-\gamma +\sum _{n=0}^{\infty }{\frac {z-1}{(n+1)(z+n)}}} ψ ( k ) ( z ) = ( − 1 ) k + 1 k ! ∑ n = 0 ∞ 1 ( z + n ) k + 1 {\displaystyle \psi ^{(k)}(z)=(-1)^{k+1}k!\sum _{n=0}^{\infty }{\frac {1}{(z+n)^{k+1}}}} これらの級数は、ガンマ関数のワイエルシュトラスの無限乗積表示 1 Γ ( z ) = z e γ z ∏ n = 1 ∞ ( 1 + z n ) e − z / n {\displaystyle {\frac {1}{\Gamma (z)}}=ze^{\gamma z}\prod _{n=1}^{\infty }{\biggl (}1+{\frac {z}{n}}{\biggr )}e^{-z/n}} の対数微分から導かれるものである、 また、z =0でのテイラー展開により、|z |<1の領域で次のように級数表示される。 ψ ( z + 1 ) = − γ + ∑ n = 2 ∞ ( − 1 ) n ζ ( n ) z n − 1 {\displaystyle \psi (z+1)=-\gamma +\sum _{n=2}^{\infty }(-1)^{n}\zeta (n)z^{n-1}} ただし、ζ(n )はリーマンゼータ関数を表す。 Rez > 0のとき、ディガンマ関数は次の積分表示を持つ。 ψ ( z ) = ∫ 0 ∞ ( e − s − 1 ( 1 + s ) z ) d s s {\displaystyle \psi (z)=\int _{0}^{\infty }{\biggl (}e^{-s}-{\frac {1}{(1+s)^{z}}}{\biggr )}{\frac {ds}{s}}} ψ ( z ) = ∫ 0 ∞ ( e − s s − e − z s 1 − e − s ) d s {\displaystyle \psi (z)=\int _{0}^{\infty }{\biggl (}{\frac {e^{-s}}{s}}-{\frac {e^{-zs}}{1-e^{-s}}}{\biggr )}ds} ψ ( z ) = − γ + ∫ 1 ∞ s z − 1 − 1 s z ( s − 1 ) d s {\displaystyle \psi (z)=-\gamma +\int _{1}^{\infty }{\frac {s^{z-1}-1}{s^{z}(s-1)}}ds} ψ ( z + 1 ) = ln ⁡ z − 1 2 z − ∫ 0 ∞ ( 1 2 coth ⁡ s 2 − 1 s ) e − z s d s {\displaystyle \psi (z+1)=\ln {z}-{\frac {1}{2z}}-\int _{0}^{\infty }{\biggl (}{\frac {1}{2}}\operatorname {coth} {\frac {s}{2}}-{\frac {1}{s}}{\biggr )}e^{-zs}ds} 但し、coths/2 は双曲線余接関数を表す。 また、ディガンマ関数同士の差について、以下が成り立つ。 ψ ( y ) − ψ ( x ) = ∫ 0 1 u x − 1 − u y − 1 1 − u d u {\displaystyle \psi (y)-\psi (x)=\int _{0}^{1}{\frac {u^{x-1}-u^{y-1}}{1-u}}du} ガンマ関数の相反公式に対し、対数微分をとることで次の関係式が導かれる。 ψ ( 1 − z ) − ψ ( z ) = π cot ⁡ π z {\displaystyle \psi (1-z)-\psi (z)=\pi \operatorname {cot} \pi z} 但し、cot πz は余接関数を表す。 z →∞　(|argz | < π)のとき、ディガンマ関数は次の漸近展開をもつ。 ψ ( z ) ∼ ln ⁡ z − 1 2 z − ∑ n = 1 ∞ B 2 n 2 n z 2 n = ln ⁡ z − 1 2 z − 1 12 z 2 + 1 120 z 4 − 1 252 z 6 + ⋯ {\displaystyle {\begin{aligned}\psi (z)&\sim \ln {z}-{\frac {1}{2z}}-\sum _{n=1}^{\infty }{\frac {B_{2n}}{2nz^{2n}}}\\&=\ln {z}-{\frac {1}{2z}}-{\frac {1}{12z^{2}}}+{\frac {1}{120z^{4}}}-{\frac {1}{252z^{6}}}+\cdots \end{aligned}}} 但し、B2nはベルヌーイ数である。 ディガンマ関数は、正の整数において、次の値をとる。 ψ ( 1 ) = − γ {\displaystyle \psi (1)=-\gamma } ψ ( n ) = − γ + ∑ k = 1 n − 1 1 k = − γ + H n − 1 ( n = 2 , 3 , 4 , ⋯ ) {\displaystyle \psi (n)=-\gamma +\sum _{k=1}^{n-1}{\frac {1}{k}}=-\gamma +H_{n-1}\qquad (n=2,3,4,\cdots )} 但し、Hn-1は調和数を表す。 また、正の半整数において、次の値をとる。 ψ ( 1 / 2 ) = − γ − 2 ln ⁡ 2 {\displaystyle \psi (1/2)=-\gamma -2\ln {2}} ψ ( n + 1 / 2 ) = − γ − 2 ln ⁡ 2 + 2 ∑ k = 0 n − 1 1 2 k + 1 ( n = 1 , 2 , 3 , ⋯ ) {\displaystyle \psi (n+1/2)=-\gamma -2\ln {2}+2\sum _{k=0}^{n-1}{\frac {1}{2k+1}}\qquad (n=1,2,3,\cdots )} Milton Abramowitz and Irene A. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables. Dover (1965) ISBN 978-0486612720 E. T. Whittaker and G. N. Watson, A Course of Modern Analysis. Cambridge University Press (1927; reprinted 1996) ISBN 978-0521588072 George B. Arfken and Hans J. Weber, Mathematical Methods for Physicists, Academic Press; ジョージ．ブラウン・アルフケン、ハンス．J・ウェーバー (著)、 権平健一郎、神原武志、小山直人 (翻訳) 『基礎物理数学第4版Vol．3　特殊関数』 講談社 (2001) ISBN 978-4061539792 ガンマ関数 ポリガンマ関数