指数積分

数学において、指数積分（しすうせきぶん、英: exponential integral）Ei は指数関数を含む積分によって定義される特殊関数の一つである。 実数 x≠0 に対し指数積分 Ei(x) は次のように定義される。 Ei ⁡ ( x ) = − p . v . ⁡ ∫ − x ∞ e − t t d t = p . v . ⁡ ∫ − ∞ x e t t d t {\displaystyle \operatorname {Ei} (x)=-\operatorname {p.\!v.} \int _{-x}^{\infty }{\frac {e^{-t}}{t}}\,\operatorname {d} \!t=\operatorname {p.\!v.} \int _{-\infty }^{x}{\frac {e^{t}}{t}}\,\operatorname {d} \!t} ただし p.v. はコーシーの主値を表す。この関数は初等関数でないことがリッシュのアルゴリズムによって示されている。 以下、本稿ではこれを Eireal(x) で表す。 E i r e a l ⁡ ( x ) = lim ϵ → + 0 ( − ∫ − x − ϵ e − t t d t − ∫ ϵ ∞ e − t t d t ) ( x > 0 ) E i r e a l ⁡ ( x ) = − ∫ − x ∞ e − t t d t ( x < 0 ) {\displaystyle {\begin{aligned}\operatorname {{Ei}^{real}} (x)&=\lim _{\epsilon \to +0}\left(-\int _{-x}^{-\epsilon }{\frac {e^{-t}}{t}}\,\operatorname {d} \!t-\int _{\epsilon }^{\infty }{\frac {e^{-t}}{t}}\,\operatorname {d} \!t\right)\quad &(x>0)\\\operatorname {{Ei}^{real}} (x)&=-\int _{-x}^{\infty }{\frac {e^{-t}}{t}}\,\operatorname {d} \!t\quad &(x<0)\end{aligned}}} 複素数 z に対し指数積分 Ei(z) は次のように定義される。 Ei ⁡ ( z ) = − π i + ∫ − ∞ − 0 i 1 − 0 i e t t d t + ∫ 1 z e t t d t {\displaystyle \operatorname {Ei} (z)=-\pi i+\int _{-\infty -0i}^{1-0i}{\frac {e^{t}}{t}}\,\operatorname {d} \!t+\int _{1}^{z}{\frac {e^{t}}{t}}\,\operatorname {d} \!t} これは多価関数であるが、本稿では負の実軸で分枝切断を行い正の実軸上で実数値をとるようにする。(文献によっては定義が異なる) Ei ⁡ ( x ± 0 i ) = E i r e a l ⁡ ( x ) ± π i ( x < 0 ) , Ei ⁡ ( x ) = E i r e a l ⁡ ( x ) ( x > 0 ) {\displaystyle \operatorname {Ei} (x\pm 0i)=\operatorname {{Ei}^{real}} (x)\pm \pi i\quad (x<0),\quad \operatorname {Ei} (x)=\operatorname {{Ei}^{real}} (x)\quad (x>0)} 複素関数 Ein(z) を次のように定める。 Ein ⁡ ( z ) = ∫ 0 z 1 − e − t t d t {\displaystyle \operatorname {Ein} (z)=\int _{0}^{z}{\frac {1-e^{-t}}{t}}\,\operatorname {d} \!t} これは複素平面全体で正則となり、 Ein ⁡ ( z ) − E 1 ( z ) − log ⁡ z = ∫ 0 z 1 − e − t t d t − ∫ z ∞ e − t t d t − ∫ 1 z 1 t d ⁡ t = ∫ 0 1 1 − e − t t d t − ∫ 1 ∞ e − t t d t = − ∫ 0 ∞ ( log ⁡ t ) e − t d t = γ {\displaystyle {\begin{aligned}\operatorname {Ein} (z)-E_{1}(z)-\log z&=\int _{0}^{z}{\frac {1-e^{-t}}{t}}\,\operatorname {d} \!t-\int _{z}^{\infty }{\frac {e^{-t}}{t}}\,\operatorname {d} \!t-\int _{1}^{z}{\frac {1}{t}}\,\operatorname {d} t\\&=\int _{0}^{1}{\frac {1-e^{-t}}{t}}\,\operatorname {d} \!t-\int _{1}^{\infty }{\frac {e^{-t}}{t}}\,\operatorname {d} \!t\\&=-\int _{0}^{\infty }(\log t)e^{-t}\,\operatorname {d} \!t\\&=\gamma \end{aligned}}} が成り立つ。ただしγはオイラーの定数である。これにより E1, Ei は E 1 ( z ) = − γ − log ⁡ z + Ein ⁡ ( z ) Ei ⁡ ( z ) = γ + log ⁡ z − Ein ⁡ ( − z ) {\displaystyle {\begin{aligned}E_{1}(z)&=-\gamma -\log z+\operatorname {Ein} (z)\\\operatorname {Ei} (z)&=\gamma +\log z-\operatorname {Ein} (-z)\end{aligned}}} と表され、多価性にまつわる問題を複素対数関数 log z に封じ込めることができる。 Ein(z) のテイラー展開は次のように与えられる。 Ein ⁡ ( z ) = ∫ 0 z ∑ k = 1 ∞ ( − t ) k − 1 k ! d t = − ∑ k = 1 ∞ ( − z ) k k k ! {\displaystyle {\begin{aligned}\operatorname {Ein} (z)&=\int _{0}^{z}\sum _{k=1}^{\infty }{\frac {(-t)^{k-1}}{k!}}\,\operatorname {d} \!t\\&=-\sum _{k=1}^{\infty }{\frac {(-z)^{k}}{k\;k!}}\end{aligned}}} これは複素平面全体で収束する。また次のような展開も可能である。 Ein ⁡ ( z ) = ∫ 0 z e − t ∑ k = 1 ∞ t k − 1 k ! d t = e − z ∑ n = 1 ∞ ( ∑ k = 1 n 1 k ) z n n ! Ein ⁡ ( z ) = ∫ 0 z e − t / 2 ∑ k = 0 ∞ ( t / 2 ) 2 k ( 2 k + 1 ) ! d t = e − z / 2 ∑ n = 1 ∞ ( ∑ k = 0 ⌊ ( n − 1 ) / 2 ⌋ 2 2 k + 1 ) ( z / 2 ) n n ! {\displaystyle {\begin{aligned}\operatorname {Ein} (z)&=\int _{0}^{z}e^{-t}\sum _{k=1}^{\infty }{\frac {t^{k-1}}{k!}}\,\operatorname {d} \!t\\&=e^{-z}\sum _{n=1}^{\infty }\left(\sum _{k=1}^{n}{\frac {1}{k}}\right){\frac {z^{n}}{n!}}\\\operatorname {Ein} (z)&=\int _{0}^{z}e^{-t/2}\sum _{k=0}^{\infty }{\frac {(t/2)^{2k}}{(2k+1)!}}\,\operatorname {d} \!t\\&=e^{-z/2}\sum _{n=1}^{\infty }\left(\sum _{k=0}^{\lfloor (n-1)/2\rfloor }{\frac {2}{2k+1}}\right){\frac {(z/2)^{n}}{n!}}\end{aligned}}} z の絶対値が十分大きいとき E1 は次のように近似できる。 E 1 ( z ) = − e − z { ∑ k = 1 n ( k − 1 ) ! ( − 1 z ) k + O ( 1 z n + 1 ) } {\displaystyle E_{1}(z)=-e^{-z}\left\{\sum _{k=1}^{n}(k-1)!\left(-{\frac {1}{z}}\right)^{k}+O\left({\frac {1}{z^{n+1}}}\right)\right\}} 右辺は n→∞ で発散するので適当な項数で打ち切って使用する。 指数積分は以下のように一般化できる。 E n ( z ) = z n − 1 ∫ z ∞ e − t t n d t = ∫ 1 ∞ e − z t t n dt {\displaystyle E_{n}(z)=z^{n-1}\int _{z}^{\infty }{\frac {e^{-t}}{t^{n}}}\,\operatorname {d} \!t=\int _{1}^{\infty }{\frac {e^{-zt}}{t^{n}}}\,\operatorname {dt} \!} これを n 次の指数積分と呼び、以下のように不完全ガンマ関数を用いて以下のように表せる。 　 E n ( x ) = x n − 1 Γ ( 1 − n , x ) {\displaystyle E_{n}(x)=x^{n-1}\Gamma (1-n,x)} また、以下の式を Ei(z) と記すこともある。 E 1 ( z ) = ∫ z ∞ e − t t d t {\displaystyle E_{1}(z)=\int _{z}^{\infty }{\frac {e^{-t}}{t}}\,\operatorname {d} \!t} このときは次のように分枝をとる。 Im ⁡ ( E n ( x ± 0 i ) ) = ∓ π ( − x ) n − 1 i ( x < 0 ) , Im ⁡ ( E n ( x ) ) = 0 ( x > 0 ) E 1 ( x ± 0 i ) = − E i r e a l ⁡ ( − x ) ∓ π i ( x < 0 ) , E 1 ( x ) = − E i r e a l ⁡ ( − x ) ( x > 0 ) {\displaystyle {\begin{aligned}&\operatorname {Im} (E_{n}(x\pm 0i))=\mp \pi (-x)^{n-1}i\quad (x<0),\quad \operatorname {Im} (E_{n}(x))=0\quad (x>0)\\&E_{1}(x\pm 0i)=-\operatorname {{Ei}^{real}} (-x)\mp \pi i\quad (x<0),\quad E_{1}(x)=-\operatorname {{Ei}^{real}} (-x)\quad (x>0)\end{aligned}}} 両者は次のような関係で結ばれる。 Ei ⁡ ( z ) = − E 1 ( − z ) + π i ( Im ⁡ ( z ) < 0 ) , Ei ⁡ ( z ) = − E 1 ( − z ) − π i ( Im ⁡ ( z ) > 0 ) {\displaystyle \operatorname {Ei} (z)=-E_{1}(-z)+\pi i\quad (\operatorname {Im} (z)<0),\quad \operatorname {Ei} (z)=-E_{1}(-z)-\pi i\quad (\operatorname {Im} (z)>0)} 指数積分は以下のような近似を持つ。 　 E 1 ( x ) = ( A − 7.7 + B ) − 0.13 {\displaystyle E_{1}(x)=(A^{-7.7}+B)^{-0.13}} ただし A = log ⁡ ( ( 0.56146 x + 0.65 ) ( 1 + x ) ) B = x 4 e 7.7 x ( 2 + x ) 3.7 {\displaystyle {\begin{aligned}A&=\log \left(\left({\frac {0.56146}{x}}+0.65\right)(1+x)\right)\\B&=x^{4}e^{7.7x}(2+x)^{3.7}\end{aligned}}} E 1 ( x ) = e − x x + 1 1 + 1 x + 2 1 + 2 x + 3 ⋱ {\displaystyle E_{1}(x)={\cfrac {e^{-x}}{x+{\cfrac {1}{1+{\cfrac {1}{x+{\cfrac {2}{1+{\cfrac {2}{x+{\cfrac {3}{\ddots }}}}}}}}}}}}} Ei ⁡ ( z ) = γ + log ⁡ z + z ⋅ 2 F 2 [ 1 , 1 2 , 2 ; z ] {\displaystyle {\begin{aligned}\operatorname {Ei} (z)&=\gamma +\log {z}+z\cdot {_{2}F_{2}}\left[{\begin{matrix}1,1\\2,2\end{matrix}};z\right]\\\end{aligned}}} 正弦積分 (sine integral) は正弦関数を含む積分によって定義される関数である。被積分関数は非正規化sinc関数という。 Si ⁡ ( z ) = ∫ 0 z sin ⁡ t t d t si ⁡ ( z ) = − ∫ z ∞ sin ⁡ t t d t = Si ⁡ ( z ) − π 2 {\displaystyle {\begin{aligned}&\operatorname {Si} (z)=\int _{0}^{z}{\frac {\sin {t}}{t}}\,\operatorname {d} \!t\\&\operatorname {si} (z)=-\int _{z}^{\infty }{\frac {\sin {t}}{t}}\,\operatorname {d} \!t=\operatorname {Si} (z)-{\frac {\pi }{2}}\end{aligned}}} 余弦積分 (cosine integral) は余弦関数を含む積分によって定義される関数である。 Ci ⁡ ( z ) = − ∫ z z + ∞ cos ⁡ t t d t {\displaystyle \operatorname {Ci} (z)=-\int _{z}^{z+\infty }{\frac {\cos {t}}{t}}\,\operatorname {d} \!t} 複素関数としての余弦積分は多価であるが、次のように複素対数関数と正則関数の和で表すことができる。 Ci ⁡ ( z ) = γ + log ⁡ z − Cin ⁡ ( z ) Cin ⁡ ( z ) = ∫ 0 z 1 − cos ⁡ t t d t {\displaystyle {\begin{aligned}\operatorname {Ci} (z)&=\gamma +\log {z}-\operatorname {Cin} (z)\\\operatorname {Cin} (z)&=\int _{0}^{z}{\frac {1-\cos {t}}{t}}\,\operatorname {d} \!t\end{aligned}}} 任意の複素数 z に対して次の関係が成り立つ。 Ein ⁡ ( ± i z ) = Cin ⁡ ( z ) ± i Si ⁡ ( z ) {\displaystyle \operatorname {Ein} (\pm iz)=\operatorname {Cin} (z)\pm i\operatorname {Si} (z)} 対数積分 (logarithmic integal) は対数関数の逆数の積分によって定義される関数である。詳しくは対数積分を参照。 Li ⁡ ( z ) = Ei ⁡ ( log ⁡ z ) − Ei ⁡ ( log ⁡ 2 ) = ∫ 2 z 1 log ⁡ t d t li ⁡ ( z ) = Ei ⁡ ( log ⁡ z ) = p . v . ⁡ ∫ 0 2 1 log ⁡ t d t + ∫ 2 z 1 log ⁡ t d t {\displaystyle {\begin{aligned}\operatorname {Li} (z)&=\operatorname {Ei} (\log {z})-\operatorname {Ei} (\log {2})=\int _{2}^{z}{\frac {1}{\log {t}}}\,\operatorname {d} \!t\\\operatorname {li} (z)&=\operatorname {Ei} (\log {z})=\operatorname {p.\!v.} \int _{0}^{2}{\frac {1}{\log {t}}}\,\operatorname {d} \!t+\int _{2}^{z}{\frac {1}{\log {t}}}\,\operatorname {d} \!t\end{aligned}}} ただし p.v. はコーシーの主値を表す。対数積分は素数の分布を表す公式(素数定理)に現れる。 ^ Wolfram Mathworld: Exponential Integral ^ SpringerLink: Integral exponential function Weisstein, Eric W. "Exponential Integral". mathworld.wolfram.com (英語). Weisstein, Eric W. "En-Function". mathworld.wolfram.com (英語).