未分類 kanrinin

ガウス積分を高校数学で計算する

ガウス積分

　\displaystyle \int_{-\infty}^{\infty}e^{-x^2}dx =\sqrt{\pi}  …①

の証明は通常、大学の微分積分の範囲で重積分の「極座標変換」の知識を使って計算されます。

 

計算したい積分をI とおくと、

　\displaystyle I^2=\left(\int_{-\infty}^{\infty}e^{-x^2}dx\right)\cdot \left(\int_{-\infty}^{\infty}e^{-y^2}dy \right)

　　\displaystyle =\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}e^{-x^2-y^2}dxdy

　　\displaystyle =\int_{0}^{2\pi} \int_{0}^{\infty}e^{-r^2}rdr d\theta

　　\displaystyle =\int_{0}^{2\pi} \left[-\frac{1}{2}e^{-r^2}\right]_{r=0}^{r=\infty}d\theta

　　\displaystyle =2\pi \cdot \frac12=\pi

 という計算から、 I=\sqrt{\pi}が示されます。

 

今回は上の計算を可能な限り高校数学の道具だけで計算してみようと思います。そもそも積分区間が\displaystyle -\infty\to \infty の広義積分である時点で高校数学の範囲を超えているわけですが、広義積分といえども有限区間の積分で最後に極限をとる

　\displaystyle \int_{-\infty}^{\infty}e^{-x^2}dx =\lim_{R\to \infty}\int_{-R}^{R}e^{-x^2}dx

ということですから高校数学の範ちゅうにしてしまえ、というノリで計算してしまいます。

 

\displaystyle e^{-x^2} は偶関数ですから、ガウス積分①を示すには

　\displaystyle \int_{0}^{\infty}e^{-x^2}dx=\frac{\sqrt{\pi}}{2} …②

を示せば十分です。さて、やや天下り的ですが、次の積分を計算してみましょう(道具があまり使えない代わりに巧みな工夫がいるのです！)

　\displaystyle \int_{0}^{\infty}\left\{\int_{0}^{\infty}xe^{-x^2(y^2+1)}dx \right\}dy …③

上③の積分を計算してみましょう。x で積分する際はy はただの定数扱いですから、

　\displaystyle \int_{0}^{\infty}\left[-\frac{1}{2(y^2+1)}e^{-x^2(y^2+1)}\right]_{x=0}^{x=\infty}dy

　\displaystyle =\int_{0}^{\infty}\frac{1}{2(y^2+1)}dy

ここで\displaystyle y=\tan\theta とおくと\displaystyle \frac{dy}{d\theta}=\frac{1}{\cos^2\theta} 、変数の対応は\displaystyle x が 0から\infty までうごくとき\theta は0 から\frac{\pi}{2} まで動くから

　\displaystyle \int_{0}^{\infty}\frac{1}{2(y^2+1)}dy =\frac12\int_{0}^{\frac{\pi}{2}}\frac{1}{\tan^2\theta +1 }\cdot \frac{1}{\cos^2\theta}d\theta

　\displaystyle \hspace{30mm}=\frac12\int_{0}^{\frac{\pi}{2}}d\theta =\frac{\pi}{4} …④

となります。さて、今度は③の計算をy から先に行ってみましょう。\displaystyle e^{-x^2(y^2+1)}=e^{-(xy)^2}\cdot e^{-x^2} と\displaystyle y で積分する際にはx は定数扱いであることに注意すると、

　\displaystyle \int_{0}^{\infty}\left\{\int_{0}^{\infty}xe^{-x^2(y^2+1)}dx \right\}dy

　\displaystyle =\int_{0}^{\infty}\left\{\int_{0}^{\infty}xe^{-x^2(y^2+1)}dy \right\}dx

　\displaystyle =\int_{0}^{\infty}xe^{-x^2}\left\{\int_{0}^{\infty}e^{-(xy)^2}dy\right\}dx

　(ここでy での積分においてxy=Y と置換すれば\displaystyle \frac{dy}{dY}=\frac{1}{x} により)

　\displaystyle =\int_{0}^{\infty}xe^{-x^2}\left\{\int_{0}^{\infty}e^{-Y^2}\cdot \frac{1}{x}dY\right\}dx

　\displaystyle =\int_{0}^{\infty}e^{-x^2}\left\{\int_{0}^{\infty}e^{-Y^2}dY\right\}dx

　\displaystyle =\left\{\int_{0}^{\infty}e^{-Y^2}dY\right\}\cdot \int_{0}^{\infty}e^{-x^2}dx

　\displaystyle =\left\{\int_{0}^{\infty}e^{-x^2}dx \right\}^2 …⑤

天下り的に定義した積分③の値は④より\displaystyle \frac{\pi}{4} であるが一方別のやり方で計算した⑤によるとこれはガウス積分(の積分区間が半分)の2乗\displaystyle \left(\int_{0}^{\infty}e^{-x^2}dx\right)^2 であることがわかりました。

よって②の式

　\displaystyle \int_{0}^{\infty}e^{-x^2}dx=\frac{\sqrt{\pi}}{2}

が成立することがわかります。

 

では。

 

 

 

 

 

RECOMMENDこちらの記事も人気です。