tecktom.com

ウォリス積分とウォリス積

2024年度 高1発展2026年10月6日12分で読めます塾講師木下積分

いよいよ高一数Ⅲも最後ですね。今回は定積分と漸化式に関係する「ウォリス積分」および「ウォリス積」というものについて紹介します。

定理(ウォリス積分・ウォリス積)

nnを非負の整数とし、

In:=∫0π/2sin⁡nxdxI_n:=\int_0^{\pi/2}\sin ^n x dx

とする。

  • 1. ∫0π/2cos⁡nxdx=In\int_0^{\pi/2} \cos ^n x dx =I_nである。

  • 2. In+2=n+1n+2InI_{n+2}=\dfrac{n+1}{n+2}I_nが成立する。

  • 3. 次が成立する。(ウォリス積分)

    In={(n−1)!!n!!(nが奇数)(n−1)!!n!!π2(nが偶数)I_n= \begin{cases} \dfrac{\left(n-1\right)!!}{n!!}\quad \text{(}n\text{が奇数)}\\ \dfrac{\left(n-1\right)!!}{n!!}\dfrac{\pi}{2}\quad \text{(}n\text{が偶数)} \end{cases}

  • 4. ∏k=1∞(2k)2(2k−1)(2k+1)=π2\prod_{k=1}^{\infty}\dfrac{\left(2k\right)^2}{\left(2k-1\right)\left(2k+1\right)}=\dfrac{\pi}{2} が成立する。(ウォリス積)

はじめに、みなさんに馴染みがないと思われる記号や概念の紹介をしましょう。

定義(二重階乗)

非負の整数mmに対して、

m!!:={m⋅(m−2)⋯3⋅1(mが奇数)m⋅(m−2)⋯4⋅2(mが偶数)m!!:= \begin{cases} m\cdot (m-2)\cdots 3\cdot 1\quad \text{(}m\text{が奇数)}\\ m\cdot (m-2)\cdots 4\cdot 2\quad \text{(}m\text{が偶数)} \end{cases}

ただし、0!!=10!!=1、(−1)!!=1(-1)!!=1と約束する(定理のn=0,1n=0,1の場合に使う)。

まずは“!!!!”の記号から。これは「二重階乗」と呼ばれ、上の定義からも分かる通り階乗の二個飛ばしバージョンです。例えば、

{5!!=5⋅3⋅1=158!!=8⋅6⋅4⋅2=384\begin{cases} 5!! =5\cdot 3\cdot 1=15\\ 8!! =8\cdot 6\cdot 4\cdot 2=384 \end{cases}

みたいな感じです。あんまむずくないしかっこいいっしょ。

定義(総乗記号)

数列{am}\left\{a_m\right\}に対して、

∏k=lnak:=al⋅al+1⋅al+2⋯an\prod_{k=l}^n a_k :=a_l\cdot a_{l+1}\cdot a_{l+2}\cdots a_n

続いて“∏\prod”の記号です。これは、総乗記号と呼ばれ、総和記号∑\sumの掛け算バージョンです。上の定義をみても分かるように、例えば、

{∏k=14k=1⋅2⋅3⋅4=24∏k=13k2=12⋅22⋅32=36\begin{cases} \prod_{k=1}^{4} k =1\cdot 2\cdot 3\cdot 4=24\\ \prod_{k=1}^{3} k^2 =1^2\cdot 2^2\cdot 3^2=36 \end{cases}

のようになります。これも見た目こそかっこいいですが難しいことをやっているわけではないですし、なんなら知っておくと記述とかで楽できる場面が結構ある気がします。(ちなみに∏\prodはギリシャ文字のπ\piの大文字です。英語で積を表す“Product”の頭文字“P”に対応しています。)

定理(はさみうちの原理)

十分大きいnnでいつでもan≦bn≦cna_n\leqq b_n\leqq c_nが成り立つ数列{am}\left\{a_m\right\}、{bm}\left\{b_m\right\}、{cm}\left\{c_m\right\}について、

{lim⁡k→∞ak=αlim⁡k→∞ck=α\begin{aligned}&\begin{cases} \lim_{k\to \infty}a_k =\alpha\\ \lim_{k\to \infty}c_k =\alpha \end{cases}\end{aligned}

が成り立つとき、

lim⁡k→∞bk=α\begin{aligned}\lim_{k\to \infty}b_k =\alpha\end{aligned}

が成立する。

上記の定理をはさみうちの原理と言います。大学数学で厳密に極限を定義するとそこから導かれる定理となりますが、高校数学のうちは言及さえすれば証明なしで用いることのできる定理です。直感的にはまあ当たり前よね。

さて、準備も済んだので、いよいよ冒頭の定理を示しに行ってみましょう。

InI_nの定義から、要は

∫0π/2cos⁡nxdx=∫0π/2sin⁡nxdx\int_0^{\pi/2} \cos ^n x dx =\int_0^{\pi/2} \sin ^n x dx

ということを言っている。これには、置換積分を使ってやるのが楽。

以前の補足で登場した“King Property”1右辺で t=a+b−xt=a+b-x と置換すると、x=ax=a のとき t=bt=b、x=bx=b のとき t=at=a、dx=−dtdx=-dt なので左辺に戻ります。この性質はKing Property の記事でも扱っています。下の脚注へ:

定理(King Property)

∫abf(x)dx=∫abf(a+b−x)dx\int_a^b f(x)dx=\int_a^b f(a+b-x)dx

を用いれば、

∫0π/2sin⁡nxdx=∫0π/2sin⁡n(π2−x)dx=∫0π/2cos⁡nxdx(∵sin⁡(π2−x)=cos⁡x)\begin{aligned}&\int_0^{\pi/2} \sin ^n x dx\\&=\int_0^{\pi/2} \sin ^n \left(\dfrac{\pi}{2}-x\right) dx\\&=\int_0^{\pi/2} \cos ^n x dx\qquad \left(\because \sin\left(\dfrac{\pi}{2}-x\right)=\cos x\right)\end{aligned}

となる。よって、

∫0π/2sin⁡nxdx=∫0π/2cos⁡nxdx\int_0^{\pi/2} \sin ^n x dx =\int_0^{\pi/2} \cos ^n x dx

□

※この解法は有名なので覚えておいてもいいかも。

In+2=∫0π/2sin⁡n+2xdx=∫0π/2sin⁡xsin⁡n+1xdx=[−cos⁡xsin⁡n+1x]0π/2+∫0π/2(n+1)sin⁡nxcos⁡2xdx(部分積分)=∫0π/2(n+1)sin⁡nx(1−sin⁡2x)dx(∵[−cos⁡xsin⁡n+1x]0π/2=0)=(n+1)∫0π/2sin⁡nxdx−(n+1)∫0π/2sin⁡n+2xdx=(n+1){In−In+2}\begin{aligned}&I_{n+2}=\int_0^{\pi/2}\sin ^{n+2} x dx\\ &=\int_0^{\pi/2}\sin x\sin ^{n+1}x dx\\ &=\left[ -\cos x \sin^{n+1}x\right]_0^{\pi/2}\\ &\qquad+\int_0^{\pi/2}\left(n+1\right)\sin^n x\cos^2 xdx\quad\text{(部分積分)}\\ &=\int_0^{\pi/2}\left(n+1\right)\sin^n x\left(1-\sin^2 x \right) dx\\ &\qquad\left(\because\left[ -\cos x \sin^{n+1}x\right]_0^{\pi/2}=0\right)\\ &=(n+1)\int_0^{\pi/2}\sin^n xdx\\ &\qquad-(n+1)\int_0^{\pi/2}\sin^{n+2} xdx\\ &=(n+1)\left\{I_n-I_{n+2}\right\}\end{aligned}

となる。つまり、

In+2=(n+1){In−In+2}⇔(n+2)In+2=(n+1)In⇔In+2=n+1n+2In\begin{aligned}&I_{n+2}=(n+1)\left\{I_n-I_{n+2}\right\}\\ &\Leftrightarrow (n+2)I_{n+2}=(n+1)I_n\\ &\Leftrightarrow I_{n+2}=\dfrac{n+1}{n+2}I_n\end{aligned}

□

定理-2.を繰り返し用いると、

In=n−1nIn−2=n−1nn−3n−2In−4=n−1nn−3n−2n−5n−4In−6=⋯={(n−1)!!n!!I1(nが奇数)(n−1)!!n!!I0(nが偶数)\begin{aligned}&I_n=\dfrac{n-1}{n}I_{n-2}\\ &=\dfrac{n-1}{n}\dfrac{n-3}{n-2}I_{n-4}\\ &=\dfrac{n-1}{n}\dfrac{n-3}{n-2}\dfrac{n-5}{n-4}I_{n-6}\\ &=\cdots\\ &=\begin{cases} \dfrac{\left(n-1\right)!!}{n!!}I_1\quad \text{(}n\text{が奇数)}\\ \dfrac{\left(n-1\right)!!}{n!!}I_0\quad \text{(}n\text{が偶数)} \end{cases}\end{aligned}

となる。(二重階乗の定義を確認せよ!)ここで、

なので、

In={(n−1)!!n!!(nが奇数)(n−1)!!n!!π2(nが偶数)I_n= \begin{cases} \dfrac{\left(n-1\right)!!}{n!!}\quad \text{(}n\text{が奇数)}\\ \dfrac{\left(n-1\right)!!}{n!!}\dfrac{\pi}{2}\quad \text{(}n\text{が偶数)} \end{cases}

□

まず、[0,π2]\left[0,\dfrac{\pi}{2}\right]で0≦sin⁡x≦10\leqq\sin x\leqq 1より、自然数nnに対して、

sin⁡2n+2x≦sin⁡2n+1x≦sin⁡2nx\begin{aligned}&\sin^{2n+2}x\leqq\sin^{2n+1}x\leqq\sin^{2n}x\end{aligned}

であるから、これを積分しても大小関係は同じく、

I2n+2≦I2n+1≦I2n⋯♢\begin{aligned}&I_{2n+2}\leqq I_{2n+1}\leqq I_{2n}\quad \cdots\diamondsuit\end{aligned}

である。ここで、定理-3.を用いると、

♢⇔(2n+1)!!(2n+2)!!π2≦(2n)!!(2n+1)!!≦(2n−1)!!(2n)!!π2\begin{aligned}&\diamondsuit\Leftrightarrow \dfrac{\left(2n+1\right)!!}{(2n+2)!!}\dfrac{\pi}{2} \leqq \dfrac{\left(2n\right)!!}{(2n+1)!!}\\ &\qquad\leqq \dfrac{\left(2n-1\right)!!}{(2n)!!}\dfrac{\pi}{2}\end{aligned}

この全辺に(2n)!!(2n−1)!!(>0)\dfrac{(2n)!!}{(2n-1)!!}(>0)をかけてやって、

(2n+1)!!(2n+2)!!(2n)!!(2n−1)!!π2≦{(2n)!!}2(2n+1)!!(2n−1)!!≦π2⋯♡\begin{aligned}&\dfrac{\left(2n+1\right)!!}{(2n+2)!!}\dfrac{(2n)!!}{(2n-1)!!}\dfrac{\pi}{2}\\ &\qquad\leqq \dfrac{\left\{\left(2n\right)!!\right\}^2}{(2n+1)!!(2n-1)!!} \leqq \dfrac{\pi}{2}\quad\cdots\heartsuit\end{aligned}

を得る。ここで、二重階乗の定義から、

(2n+1)!!(2n+2)!!(2n)!!(2n−1)!!=(2n+1)(2n+2)(2n−1)!!(2n)!!(2n)!!(2n−1)!!=(2n+1)(2n+2)\begin{aligned}&\dfrac{\left(2n+1\right)!!}{(2n+2)!!}\dfrac{(2n)!!}{(2n-1)!!}\\ &=\dfrac{(2n+1)}{(2n+2)}\dfrac{\left(2n-1\right)!!}{(2n)!!}\dfrac{(2n)!!}{(2n-1)!!}\\ &=\dfrac{(2n+1)}{(2n+2)}\end{aligned}

となるから、

♡⇔(2n+1)(2n+2)π2≦{(2n)!!}2(2n+1)!!(2n−1)!!≦π2\begin{aligned}&\heartsuit \Leftrightarrow \dfrac{(2n+1)}{(2n+2)}\dfrac{\pi}{2}\\ &\qquad\leqq \dfrac{\left\{\left(2n\right)!!\right\}^2}{(2n+1)!!(2n-1)!!} \leqq \dfrac{\pi}{2}\end{aligned}

となる。

lim⁡n→∞(2n+1)(2n+2)π2=lim⁡n→∞(1+12n)(1+22n)π2=π2\begin{aligned}&\lim_{n\to \infty}\dfrac{(2n+1)}{(2n+2)}\dfrac{\pi}{2}\\ &=\lim_{n\to \infty}\dfrac{(1+\frac{1}{2n})}{(1+\frac{2}{2n})}\dfrac{\pi}{2}=\dfrac{\pi}{2}\end{aligned}

なので、はさみうちの原理から、

lim⁡n→∞{(2n)!!}2(2n+1)!!(2n−1)!!=π2\begin{aligned}&\lim_{n\to \infty}\dfrac{\left\{\left(2n\right)!!\right\}^2}{(2n+1)!!(2n-1)!!}=\dfrac{\pi}{2}\end{aligned}

が分かる。二重階乗と総乗記号の定義から、

{(2n)!!}2(2n+1)!!(2n−1)!!=∏k=1n(2k)2(2k−1)(2k+1)\begin{aligned}&\dfrac{\left\{\left(2n\right)!!\right\}^2}{(2n+1)!!(2n-1)!!}\\ &=\prod_{k=1}^{n}\dfrac{\left(2k\right)^2}{\left(2k-1\right)\left(2k+1\right)}\end{aligned}

であるので、

∏k=1∞(2k)2(2k−1)(2k+1)=π2\begin{aligned}&\prod_{k=1}^{\infty}\dfrac{\left(2k\right)^2}{\left(2k-1\right)\left(2k+1\right)}=\dfrac{\pi}{2}\end{aligned}

□

なかなかハードな計算でしたね。ただ、このウォリス積は、バラしてみると

∏k=1∞(2k)2(2k−1)(2k+1)=2⋅21⋅3×4⋅43⋅5×6⋅65⋅7×⋯\begin{aligned}&\prod_{k=1}^{\infty}\dfrac{\left(2k\right)^2}{\left(2k-1\right)\left(2k+1\right)}\\ &=\dfrac{2\cdot 2}{1\cdot 3}\times\dfrac{4\cdot 4}{3\cdot 5}\times\dfrac{6\cdot 6}{5\cdot 7}\times\cdots\end{aligned}

であり、つまり各項が有理数のものを掛け合わせていくとπ\piが出てくるんですね。なんとも不思議で美しい公式です。

また、一般には三角関数の高次(4次以上)の累乗の積分は

{ひとつ分離して置換積分で多項式の積分に帰着(奇数乗)半角公式でばらしまくる(偶数乗)\begin{aligned}&\begin{cases} \text{ひとつ分離して置換積分で多項式の積分に帰着(奇数乗)}\\ \text{半角公式でばらしまくる(偶数乗)} \end{cases}\end{aligned}

が定石で、その計算は結構大変だったりすることが多いんですが、定理-3.で示されたように、積分区間が0→π/20\to \pi/2であるような定積分に限ってはこのウォリス積分を用いれば一瞬で計算できてしまいます2たとえば ∫0π/2sin⁡6x dx=5⋅3⋅16⋅4⋅2⋅π2=5π32\displaystyle\int_0^{\pi/2}\sin^6 x\,dx=\frac{5\cdot3\cdot1}{6\cdot4\cdot2}\cdot\frac{\pi}{2}=\frac{5\pi}{32}、∫0π/2cos⁡5x dx=4⋅25⋅3⋅1=815\displaystyle\int_0^{\pi/2}\cos^5 x\,dx=\frac{4\cdot2}{5\cdot3\cdot1}=\frac{8}{15} です(cos⁡\cos は定理-1. により sin⁡\sin と同じ値です)。下の脚注へ。(僕はよくこれを使ってラクしてました。)

1 右辺で t=a+b−xt=a+b-x と置換すると、x=ax=a のとき t=bt=b、x=bx=b のとき t=at=a、dx=−dtdx=-dt なので左辺に戻ります。この性質はKing Property の記事でも扱っています。

2 たとえば ∫0π/2sin⁡6x dx=5⋅3⋅16⋅4⋅2⋅π2=5π32\displaystyle\int_0^{\pi/2}\sin^6 x\,dx=\frac{5\cdot3\cdot1}{6\cdot4\cdot2}\cdot\frac{\pi}{2}=\frac{5\pi}{32}、∫0π/2cos⁡5x dx=4⋅25⋅3⋅1=815\displaystyle\int_0^{\pi/2}\cos^5 x\,dx=\frac{4\cdot2}{5\cdot3\cdot1}=\frac{8}{15} です(cos⁡\cos は定理-1. により sin⁡\sin と同じ値です)。

塾の授業で配った補足資料をもとに書き直しました。

← 塾講師木下の記事一覧へ