1

むかしわからなかった積分の解法②

209
0
$$$$

 自分が加入・所属する二つの団体の今年度アドベント・カレンダー用に寄稿しました、今回は “Victor Adamchik” というとある工学の博士の方により示された美しいある広義積分の一般形についてご紹介します。
 なお、式中$W$につきましては「ランベルトのW関数(/対数積関数)」というものでして定義等の確認は次のMathlogの記事 『LambertのW関数』 がオススメです。

$$\int_{-\infty}^{+\infty}\frac{a^2}{(e^x-(|a|x-b))^2+(a\pi)^2}\,dx=\frac{1}{W(e^{b/|a|}/|a|)+1}\quad\text{for }a,b\in\mathbb{R}.$$

 最初に、$x\mapsto\ln x$として置換して積分区間を半分にする。そして、次に部分分数分解を行う。
$$ \begin{align} \int_{-\infty}^{+\infty}\frac{a^2}{(e^x-(|a|x-b))^2+(a\pi)^2}\,dx&=\int_0^\infty\frac{a^2}{(x-(|a|\ln x-b))^2+(a\pi)^2}\cdot\frac{dx}{x} \\ &=\int_0^\infty\frac{|a|^2}{((x-(|a|\ln x-b))+i|a|\pi)((x-(|a|\ln x-b))-i|a|\pi)}\cdot\frac{dx}{x} \\ &=\int_0^\infty\frac{1}{2i\pi}\left(-\frac{|a|}{(x-(|a|\ln x-b))+i|a|\pi}+\frac{|a|}{(x-(|a|\ln x-b))-i|a|\pi}\right)\frac{dx}{x} \\ &=\frac{1}{2\pi i}\left(-\int_0^\infty\frac{|a|}{x-|a|(\ln x-i\pi)+b}\cdot\frac{dx}{x}+\int_0^\infty\frac{|a|}{x-|a|(\ln x+i\pi)+b}\cdot\frac{dx}{x}\right) \\ &=\frac{1}{2\pi i}\left(-\int_0^\infty\frac{|a|}{x-|a|\ln(e^{-i\pi}x)+b}\cdot\frac{dx}{x}+\int_0^\infty\frac{|a|}{x-|a|\ln(e^{+i\pi}x)+b}\cdot\frac{dx}{x}\right) \\ &=\frac{1}{2\pi i}\left(+\int_0^\infty\frac{|a|}{|a|\ln(e^{-i\pi}x)-x-b}\cdot\frac{dx}{x}-\int_0^\infty\frac{|a|}{|a|\ln(e^{+i\pi}x)-x-b}\cdot\frac{dx}{x}\right). \end{align} $$
 ここで、自然対数の分枝を考慮しつつオイラーの等式$e^{+i\pi}=e^{-i\pi}=-1$を考慮して、次の複素積分を考察する。続くその次の図は積分経路$\varGamma$の図示。
$$ \begin{align} \oint_\varGamma\frac{|a|}{|a|\ln(-z)-z-b}\cdot\frac{dz}{z}&\coloneqq\left(\underline{\int_{\ell_+\mathstrut}}_{(\mathrm{I})}+\underline{\int_{C_R\mathstrut}}_{(\mathrm{I\hspace{-.1em}I})}+\underline{\int_{\ell_-\mathstrut}}_{(\mathrm{I\hspace{-.1em}I\hspace{-.1em}I})}+\underline{\int_{c_\varepsilon\mathstrut}}_{(\mathrm{I\hspace{-.1em}V})}\right)\frac{|a|}{|a|\ln(-z)-z-b}\cdot\frac{dz}{z} \\ &\coloneqq\left(\underline{\int_{\varepsilon\mathstrut}^R}_{(\mathrm{I})}+\underline{\int_{z=Re^{+i\theta}\colon0\leq\theta\leq2\pi\mathstrut}}_{(\mathrm{I\hspace{-.1em}I})}+\underline{\int_{R\mathstrut}^\varepsilon}_{(\mathrm{I\hspace{-.1em}I\hspace{-.1em}I})}+\underline{\int_{z=\varepsilon e^{-i\theta}\colon0\leq\theta\leq2\pi\mathstrut}}_{(\mathrm{I\hspace{-.1em}V})}\right)\frac{|a|}{|a|\ln(-z)-z-b}\cdot\frac{dz}{z}.\tag*{$\cdots(\mathrm{V})$} \end{align} $$
積分経路!FORMULA[7][2074314156][0]と特異点達!FORMULA[8][1751739959][0] 積分経路$\varGamma$と特異点達$\Large\times$
原点$O$以外の孤立特異点の計算過程を次に示しておく、非自明の場合を考えているので無論以降も$a\neq0$である。
$$ \begin{align} |a|\ln(-z)-z-b=0&\Longleftrightarrow|a|(\ln(-z/|a|)+(-z/|a|))+|a|(\ln|a|-b/|a|)=0 \\ &\Longleftrightarrow|a|\ln((-z/|a|)e^{-z/|a|})+|a|\ln{}(|a|/e^{b/|a|})=0 \\ &\Longleftrightarrow\ln((-z/|a|)e^{-z/|a|})+\ln{}(|a|/e^{b/|a|})=0 \\ &\Longleftrightarrow\ln((-z/|a|)e^{-z/|a|})=\ln{}(e^{b/|a|}/|a|) \\ &\Longleftrightarrow(-z/|a|)e^{-z/|a|}=e^{b/|a|}/|a| \\ &\Longleftrightarrow(-z/|a|)=W(e^{b/|a|}/|a|) \\ &\Longleftrightarrow z=-|a|\cdot W(e^{b/|a|}/|a|). \end{align} $$
 最左辺の周回積分はコーシーの留数定理が適用可能、式途中の極限計算ではロピタルの定理を使用。
$$ \begin{align} \oint_\varGamma\frac{|a|}{|a|\ln(-z)-z-b}\cdot\frac{dz}{z}=2\pi i&\cdot\mathop{\mathrm{Res}}_{z=-|a|\cdot W(e^{b/|a|}/|a|)}\left(\frac{|a|}{|a|\ln(-z)-z-b}\cdot\frac{1}{z}\right) \\ =2\pi i\cdot\frac{1}{(1-1)!}&\cdot\lim_{z\to-|a|\cdot W(e^{b/|a|}/|a|)}\left(\frac{d^{1-1}}{dz^{1-1}}\left((z-(-|a|\cdot W(e^{b/|a|}/|a|)))^1\cdot\left(\frac{|a|}{|a|\ln(-z)-z-b}\cdot\frac{1}{z}\right)\right)\right) \\ =2\pi i&\cdot\lim_{z\to-|a|\cdot W(e^{b/|a|}/|a|)}\left((z+|a|\cdot W(e^{b/|a|}/|a|))\cdot\frac{|a|}{|a|\ln(-z)-z-b}\cdot\frac{1}{z}\right) \\ =2\pi i&\cdot\lim_{z\to-|a|\cdot W(e^{b/|a|}/|a|)}\left(\frac{z+|a|\cdot W(e^{b/|a|}/|a|)}{|a|\ln(-z)-z-b}\right)\cdot\left(-\frac{1}{W(e^{b/|a|}/|a|)}\right) \\ =2\pi i&\cdot\lim_{z\to-|a|\cdot W(e^{b/|a|}/|a|)}\left(\frac{\displaystyle\frac{d}{dz}(z+|a|\cdot W(e^{b/|a|}/|a|))}{\displaystyle\frac{d}{dz}(|a|\ln(-z)-z-b)}\right)\cdot\left(-\frac{1}{W(e^{b/|a|}/|a|)}\right) \\ =2\pi i&\cdot\left(-\frac{W(e^{b/|a|}/|a|)}{W(e^{b/|a|}/|a|)+1}\right)\cdot\left(-\frac{1}{W(e^{b/|a|}/|a|)}\right) \\ =2\pi i&\cdot\frac{1}{W(e^{b/|a|}/|a|)+1}. \end{align} $$
 よって、$R\to\infty,\varepsilon\to0$による積分の極限値が存在することとなる。
$$\lim_{R\to\infty,\varepsilon\to0}\oint_\varGamma\frac{|a|}{|a|\ln(-z)-z-b}\cdot\frac{dz}{z}=2\pi i\cdot\frac{1}{W(e^{b/|a|}/|a|)+1}.$$
 また、最右辺については第$2$$\mathrm{I\hspace{-.1em}I}$と第$4$$\mathrm{I\hspace{-.1em}V}$の積分の絶対値は次のとおり上から評価される。
$$ \begin{align} 0\leq\left|\int_{|z|=r\neq1}\frac{|a|}{|a|\ln(-z)-z-b}\cdot\frac{dz}{z}\right|&\leq\int_{|z|=r\neq1}\left|\frac{|a|}{|a|\ln(-z)-z-b}\right|\cdot\left|\frac{dz}{z}\right| \\ &\leq\int_{|z|=r\neq1}\left|\frac{dz}{z}\right|\cdot\sup_{|z|=r\neq1}\left|\frac{|a|}{|a|\ln(-z)-z-b}\right| \\ &=\int_{|z|=r\neq1}\left|\frac{dz}{z}\right|\cdot\sup_{|z|=r\neq1}\left(\frac{|a|}{||a|\ln z+z-b|}\right) \\ &\leq\int_{|z|=r\neq1}\left|\frac{dz}{z}\right|\cdot\sup_{|z|=r\neq1}\left(\frac{|a|}{|||a|\ln z+z|-|b||}\right) \\ &\leq\int_{|z|=r\neq1}\left|\frac{dz}{z}\right|\cdot\sup_{|z|=r\neq1}\left(\frac{|a|}{|||a||\ln z|-|z||-|b||}\right) \\ &=2\pi\cdot{}\sup_{\!\!|\theta|\leq\pi}\left(\frac{|a|}{\left|\left||a|\sqrt{\ln\!{}^2\,r+\theta^2\mathstrut}-r\right|-|b|\right|}\right) \\ &= \left\{ \begin{aligned} &2\pi\cdot{}\sup_{\!\!|\theta|\leq\pi}\left(\frac{|a|/r}{\left|\left||a|\sqrt{((\ln r)/r)^2+(\theta/r)^2\mathstrut}-1\right|-|b|/r\right|}\right)&, \\ &2\pi\cdot{}\sup_{\!\!|\theta|\leq\pi}\left(\frac{|a|/|\ln r|}{\left|\left||a|\sqrt{1+(\theta/(\ln r))^2\mathstrut}-r/|\ln r|\right|-|b|/|\ln r|\right|}\right) \end{aligned} \right\} \xrightarrow[\smash{\enspace r\to0\hspace{0.5em}\enspace}]{\enspace r\to\infty\enspace}0. \end{align} $$
以上の評価より、以下の極限がえられる。
$$\lim_{R\to\infty,\varepsilon\to0}\int_{C_R}\frac{|a|}{|a|\ln(-z)-z-b}\cdot\frac{dz}{z}=\lim_{R\to\infty,\varepsilon\to0}\int_{c_\varepsilon}\frac{|a|}{|a|\ln(-z)-z-b}\cdot\frac{dz}{z}=0.$$
 これらより、式$\mathrm{V}$から次が帰結される。
$$ \begin{align} &\phantom{\lim_{R\to\infty,\varepsilon\to0}}\oint_\varGamma\frac{|a|}{|a|\ln(-z)-z-b}\cdot\frac{dz}{z}\hspace{-2em}&&=\phantom{\lim_{R\to\infty,\varepsilon\to0}}\int_{\ell_+\sqcup\ell_-}\frac{|a|}{|a|\ln(-z)-z-b}\cdot\frac{dz}{z}+\phantom{\lim_{R\to\infty,\varepsilon\to0}}\int_{C_R\sqcup c_\varepsilon}\frac{|a|}{|a|\ln(-z)-z-b}\cdot\frac{dz}{z} \\ \Longrightarrow&\lim_{R\to\infty,\varepsilon\to0}\oint_\varGamma\frac{|a|}{|a|\ln(-z)-z-b}\cdot\frac{dz}{z}\hspace{-2em}&&=\lim_{R\to\infty,\varepsilon\to0}\int_{\ell_+\sqcup\ell_-}\frac{|a|}{|a|\ln(-z)-z-b}\cdot\frac{dz}{z}+\lim_{R\to\infty,\varepsilon\to0}\int_{C_R\sqcup c_\varepsilon}\frac{|a|}{|a|\ln(-z)-z-b}\cdot\frac{dz}{z} \\ \Longleftrightarrow&\lim_{R\to\infty,\varepsilon\to0}\oint_\varGamma\frac{|a|}{|a|\ln(-z)-z-b}\cdot\frac{dz}{z}\hspace{-2em}&&=\lim_{R\to\infty,\varepsilon\to0}\int_{\ell_+\sqcup\ell_-}\frac{|a|}{|a|\ln(-z)-z-b}\cdot\frac{dz}{z}+\lim_{R\to\infty,\varepsilon\to0}\int_{C_R\sqcup c_\varepsilon}\frac{|a|}{|a|\ln(-z)-z-b}\cdot\frac{dz}{z} \\ \Longleftrightarrow&\mathop{}2\pi i\cdot\frac{1}{W(e^{b/|a|}/|a|)+1}\hspace{-2em}&&=\lim_{R\to\infty,\varepsilon\to0}\int_{\ell_+\sqcup\ell_-}\frac{|a|}{|a|\ln(-z)-z-b}\cdot\frac{dz}{z}+0+0 \\ \Longleftrightarrow&\mathop{}2\pi i\cdot\frac{1}{W(e^{b/|a|}/|a|)+1}\hspace{-2em}&&=\lim_{R\to\infty,\varepsilon\to0}\left(\int_{\ell_+}\frac{|a|}{|a|\ln(-z)-z-b}\cdot\frac{dz}{z}+\int_{\ell_-}\frac{|a|}{|a|\ln(-z)-z-b}\cdot\frac{dz}{z}\right) \\ \Longleftrightarrow&\mathop{}2\pi i\cdot\frac{1}{W(e^{b/|a|}/|a|)+1}\hspace{-2em}&&=\lim_{R\to\infty,\varepsilon\to0}\left(\int_\varepsilon^R\frac{|a|}{|a|\ln(-z)-z-b}\cdot\frac{dz}{z}+\int_R^\varepsilon\frac{|a|}{|a|\ln(-e^{2\pi i}z)-e^{2\pi i}z-b}\cdot\frac{dz}{z}\right) \\ \Longleftrightarrow&\mathop{}2\pi i\cdot\frac{1}{W(e^{b/|a|}/|a|)+1}\hspace{-2em}&&=\phantom{\lim_{R\to\infty,\varepsilon\to0}}\int_0^\infty\frac{|a|}{|a|\ln(-z)-z-b}\cdot\frac{dz}{z}+\int_\infty^0\frac{|a|}{|a|\ln(-e^{2\pi i}z)-e^{2\pi i}z-b}\cdot\frac{dz}{z} \\ \Longleftrightarrow&\mathop{}2\pi i\cdot\frac{1}{W(e^{b/|a|}/|a|)+1}\hspace{-2em}&&=\phantom{\lim_{R\to\infty,\varepsilon\to0}}\underline{\underline{+\int_0^\infty\frac{|a|}{|a|\ln(e^{-i\pi}x)-x-b}\cdot\frac{dx}{x}-\int_0^\infty\frac{|a|}{|a|\ln(e^{+i\pi}x)-x-b}\cdot\frac{dx}{x}}}_{(\mathrm{V\hspace{-.1em}I})}. \end{align} $$
 ゆえに、最後に式$\mathrm{V\hspace{-.1em}I}$の部分を代入すれば所望の積分をえる。
$$ \begin{align} \int_{-\infty}^{+\infty}\frac{a^2}{(e^x-(|a|x-b))^2+(a\pi)^2}\,dx&=\frac{1}{2\pi i}\left(+\int_0^\infty\frac{|a|}{|a|\ln(e^{-i\pi}x)-x-b}\cdot\frac{dx}{x}-\int_0^\infty\frac{|a|}{|a|\ln(e^{+i\pi}x)-x-b}\cdot\frac{dx}{x}\right) \\ &=\frac{1}{2\pi i}\times\left(2\pi i\cdot\frac{1}{W(e^{b/|a|}/|a|)+1}\right) \\ &=\frac{1}{W(e^{b/|a|}/|a|)+1}. \end{align} $$
$\square$

 早速ですけど、パラメータ$a=1,b=0$を引き渡してみましょう。すると、綺麗な結果が出現しますがこれが冒頭の積分でした。

$$\int_{-\infty}^{+\infty}\frac{dx}{(e^x-x)^2+\pi^2}=\frac{1}{W(1)+1}.\tag*{$\cdots(0)$}$$

 折角ですので、所有PCのPythonで近似値を計算してみます。次に示すコードブロックは、Windows PowerShellで実際に入力した画面の再現です。最終行の値が今求めている近似値です、SciPyパッケージを利用すれば初めからランベルトのW関数が実装されていて便利でした。

      PS> python -m pip install scipy

Installing collected packages: scipy
Successfully installed scipy-1.14.1

>>> import numpy as np
>>> from scipy.special import lambertw
>>> value = 1 / (lambertw(1) + 1)
>>> print(value)

(0.6381037433651108+0j)
    

 さらに、一般形からただちに多くの系がえられます。ただし、$n\in\mathbb{N}$と約束します。

$$ \begin{align} &\int_{-\infty}^{+\infty}\frac{e^x}{(e^x-x)^2+\pi^2}\,dx&&=&&\frac{1}{W(1)+1}&&,\tag*{$\cdots(1)$} \\ &\int_{-\infty}^{+\infty}\frac{dx}{(e^x-x+1)^2+\pi^2}&&=&&\frac{1}{2}&&,\tag*{$\cdots(2)$} \\ &\int_{-\infty}^{+\infty}\frac{e^x}{(e^x-x+1)^2+\pi^2}\,dx&&=&&\frac{1}{2}&&,\tag*{$\cdots(3)$} \\ &\int_{-\infty}^{+\infty}\frac{dx}{(e^x-nx+1)^2+(n\pi)^2}&&=&&\frac{1}{n}-\frac{1}{n+1}&&,\tag*{$\cdots(4)$} \\ &\int_{-\infty}^{+\infty}\frac{e^x}{(e^x-nx+1)^2+(n\pi)^2}\,dx&&=&&\frac{1}{n+1}&&,\tag*{$\cdots(5)$} \\ &\int_{-\infty}^{+\infty}\frac{dx}{(e^x-x+e+1)^2+\pi^2}&&=&&\frac{1}{e+1}&&,\tag*{$\cdots(6)$} \\ &\int_{-\infty}^{+\infty}\frac{e^x}{(e^x-x+e+1)^2+\pi^2}\,dx&&=&&\frac{1}{e+1}&&,\tag*{$\cdots(7)$} \\ &\int_{-\infty}^{+\infty}\frac{dx}{(e^x-ex+e^2+2e)^2+(e\pi)^2}&&=&&\frac{1}{e^2}-\frac{1}{e}+\frac{1}{e+1}&&,\tag*{$\cdots(8)$} \\ &\int_{-\infty}^{+\infty}\frac{e^x}{(e^x-ex+e^2+2e)^2+(e\pi)^2}\,dx&&=&&\frac{1}{e}-\frac{1}{e+1}&&,\tag*{$\cdots(9)$} \\ &\int_{-\infty}^{+\infty}\frac{2e^x-x}{((e^x-x+1)^2+\pi^2)^2}\,dx&&=&&\frac{1}{16}&&.\tag*{$\cdots(10)$} \end{align} $$

 じつは、他にも同様の解法で別の形の積分がえられます。

$$ \begin{align} &\int_{-\infty}^{+\infty}\frac{dx}{(e^x+x+1)^2+\pi^2}&&=&&\frac{2}{3}&&,\tag*{$\cdots(11)$} \\ &\int_{-\infty}^{+\infty}\frac{e^x}{(e^x+x+1)^2+\pi^2}\,dx&&=&&\frac{1}{3}&&.\tag*{$\cdots(12)$} \end{align} $$

投稿日:15日前
OptHub AI Competition

この記事を高評価した人

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

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

バッジはありません。

投稿者

現在は20代エンジニアです、、、数学のなかでもとくに函数解析に興味がありますが、大学では可換環論を勉強していました。みかけによらず中国語が苦手です。。。

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中