指数積分

数学において、指数積分(しすうせきぶん、: exponential integralEi指数関数を含む積分によって定義される特殊関数の一つである。

定義

実関数としての指数積分

実数 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}

これは多価関数であるが、本稿では負の実軸で分枝切断を行い正の実軸上で実数値をとるようにする。[1][2](文献によっては定義が異なる)

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)}

近似

指数積分は以下のような近似を持つ。

SwameeとOhijaの近似

  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.コーシーの主値を表す。対数積分は素数の分布を表す公式(素数定理)に現れる。

出典

  1. ^ Wolfram Mathworld: Exponential Integral
  2. ^ SpringerLink: Integral exponential function

外部リンク

  • Weisstein, Eric W. "Exponential Integral". mathworld.wolfram.com (英語).
  • Weisstein, Eric W. "En-Function". mathworld.wolfram.com (英語).