留数定理を用いた三角関数の積分

問題

次の積分を求めよ。

  1. ∫02πdt2+cos⁡t\displaystyle \int_0^{2\pi} \dfrac{dt}{2+\cos t}
  2. ∫02πcos⁡2θ1−2acos⁡θ+a2dθ(∣a∣<1)\displaystyle \int_0^{2\pi} \dfrac{\cos 2\theta}{1-2 a \cos \theta +a^2} d\theta \quad (|a| < 1)
  3. ∫0∞sin⁡xx  dx\displaystyle \int_0^{\infty} \dfrac{\sin x}{x} \; dx
  4. ∫−∞∞sin⁡x2  dx\displaystyle \int_{-\infty}^{\infty} \sin x^2 \; dx
  5. ∫0∞e−x2cos⁡(ax)  dx(a>0)\displaystyle \int_0^{\infty} e^{-x^2} \cos (ax) \; dx \quad (a>0)

この記事では 留数定理 の応用として,三角関数を含む実積分の計算方法を紹介します。

1.三角関数の有理式

z=eitz=e^{it} とすると cos⁡t=12(z+1z),sin⁡t=12i(z−1z)\cos t = \dfrac{1}{2} \left(z+\dfrac{1}{z}\right) , \sin t = \dfrac{1}{2i} \left(z-\dfrac{1}{z}\right) となります。これを用いて被積分関数を zz の有理式に変形し積分をします。

解

z=eitz=e^{it} とおく。このとき dz=ieit  dtdz = ie^{it} \;dt である。よってdt=dzizdt = \dfrac{dz}{iz} である。また cos⁡t=12(z+1z)\cos t = \dfrac{1}{2} \left(z+\dfrac{1}{z}\right) である。

この変数変換によって ∫02πdt2+cos⁡t=∮∣z∣=112+12(z+1z)⋅dziz=2i∮∣z∣=1dzz2+4z+1\begin{aligned} \int_0^{2\pi} \dfrac{dt}{2+\cos t} &= \oint_{|z|=1} \dfrac{1}{2+ \frac{1}{2}\left(z+\frac{1}{z}\right)} \cdot \dfrac{dz}{iz}\\ &= \dfrac{2}{i} \oint_{|z|=1} \dfrac{dz}{z^2+4z+1} \end{aligned} が得られる。z2+4z+1=0z^2+4z+1=0 の解は z=−2±3z = -2 \pm \sqrt{3} である。このうち ∣z∣<1|z| < 1 に含まれるものは z=−2+3z=-2+\sqrt{3} である。こうして留数定理から ∮∣z∣=1dzz2+4z+1=2πi  Res(1z2+4z+1,−2+3)=2πilim⁡z→−2+31z+(2+3)=πi3\begin{aligned} \oint_{|z|=1} \dfrac{dz}{z^2+4z+1} &= 2 \pi i \; \mathrm{Res} \left( \dfrac{1}{z^2+4z+1} , -2+\sqrt{3}\right)\\ &= 2\pi i \lim_{z \to -2 + \sqrt{3}} \dfrac{1}{z+ (2+\sqrt{3})}\\ &= \dfrac{\pi i}{\sqrt{3}} \end{aligned} を得る。

以上より ∫02πdt2+cos⁡t=2i∮∣z∣=1dzz2+4z+1=23π \int_0^{2\pi} \dfrac{dt}{2+\cos t} = \dfrac{2}{i} \oint_{|z|=1} \dfrac{dz}{z^2+4z+1} = \dfrac{2}{\sqrt{3}} \pi

となる。

2.三角関数の有理式パート2

まずは先ほどと同じように計算してみましょう。

解

z=eiθz=e^{i\theta} とおくと dθ=dzizd\theta = \dfrac{dz}{iz} である。また cos⁡2θ=12(z2+1z2)\cos 2\theta = \dfrac{1}{2} \left(z^2+\dfrac{1}{z^2}\right) である。

この変数変換によって ∫02πcos⁡2θ1−2acos⁡θ+a2dθ=∮∣z∣=1z4+12z211−a(z+1z)+a2dziz=∮∣z∣=1z4+1z2(z−a)(1−az)dz2i\begin{aligned} &\int_0^{2\pi} \dfrac{\cos 2\theta}{1-2 a \cos \theta +a^2} d\theta\\ &= \oint_{|z|=1} \dfrac{z^4+1}{2z^2} \dfrac{1}{1-a (z+\frac{1}{z})+a^2} \dfrac{dz}{iz}\\ &= \oint_{|z|=1} \dfrac{z^4+1}{z^2 (z-a)(1-az)} \dfrac{dz}{2i} \end{aligned} となる。∣a∣<1|a| < 1 より ∣z∣=1|z| = 1 に囲われた領域における被積分関数の特異点は z=0,az=0,a となる。

Res(z4+1z2(z−a)(1−az),0)=lim⁡z→0ddzz4+1(z−a)(1−az)=lim⁡z→0−(z4+1)(1+a2−2az)+4z3(z−a)(1−az)(z−a)2(1−az)2=−1+a2a2Res(z4+1z2(z−a)(1−az),a)=lim⁡z→az4+1z2(1−az)=a4+1a2(1−a2)\begin{aligned} &\mathrm{Res} \left( \dfrac{z^4+1}{z^2 (z-a)(1-az)} , 0 \right)\\ &= \lim_{z\to 0} \dfrac{d}{dz} \dfrac{z^4+1}{(z-a)(1-az)}\\ &= \lim_{z \to 0} \dfrac{-(z^4+1) (1+a^2-2az) + 4z^3 (z-a)(1-az)}{(z-a)^2(1-az)^2}\\ &= -\dfrac{1+a^2}{a^2}\\ &\mathrm{Res} \left( \dfrac{z^4+1}{z^2 (z-a)(1-az)} , a \right)\\ &= \lim_{z \to a} \dfrac{z^4+1}{z^2 (1-az)}\\ &= \dfrac{a^4+1}{a^2 (1-a^2)} \end{aligned} となるため ∫02πcos⁡2θ1−2acos⁡θ+a2dθ=∮∣z∣=1z4+1z2(z−a)(1−az)dz2i=2πi⋅12i(Res(z4+1z2(z−a)(1−az),0)+Res(z4+1z2(z−a)(1−az),a))=π(−1+a2a2+a4+1a2(1−a2))=2πa21−a2\begin{aligned} &\int_0^{2\pi} \dfrac{\cos 2\theta}{1-2 a \cos \theta +a^2} d\theta\\ &= \oint_{|z|=1} \dfrac{z^4+1}{z^2 (z-a)(1-az)} \dfrac{dz}{2i}\\ &= 2\pi i \cdot \dfrac{1}{2i} \left(\mathrm{Res} \left( \dfrac{z^4+1}{z^2 (z-a)(1-az)} , 0 \right) + \mathrm{Res} \left( \dfrac{z^4+1}{z^2 (z-a)(1-az)} , a \right) \right)\\ &= \pi \left(-\dfrac{1+a^2}{a^2} + \dfrac{a^4+1}{a^2 (1-a^2)}\right)\\ &= \dfrac{2\pi a^2}{1-a^2} \end{aligned} と計算される。

別解

先ほどの方法では計算がやや煩雑になります。同様の方法で分子の cos⁡2θ\cos 2\theta を cos⁡nθ\cos n\theta に変えた積分を計算しようとすると,z=0z=0 で nn 次の極を持つことになり計算がさらに難化しそうです。

実はこの問題は オイラーの公式 を用いると簡単に計算できます。この方法により,分子の cos⁡2θ\cos 2\theta を cos⁡nθ\cos n\theta に変えた積分も簡単に計算できます。

別解

In=∫02πcos⁡nθ1−2acos⁡θ+a2dθJn=∫02πsin⁡nθ1−2acos⁡θ+a2dθ I_n = \int_0^{2\pi} \dfrac{\cos n\theta}{1-2 a \cos \theta +a^2} d\theta\\ J_n = \int_0^{2\pi} \dfrac{\sin n\theta}{1-2 a \cos \theta +a^2} d\theta とおく。このとき In+iJn=∫02πcos⁡nθ+isin⁡nθ1−2acos⁡θ+a2dθ=∫02πeinθ1−2acos⁡θ+a2dθ\begin{aligned} I_n + i J_n &= \int_0^{2\pi} \dfrac{\cos n\theta + i\sin n\theta}{1-2 a \cos \theta +a^2} d\theta\\ &= \int_0^{2\pi} \dfrac{e^{in\theta}}{1-2 a \cos \theta +a^2} d\theta \end{aligned} となる。

z=eiθz=e^{i\theta} とおくことで (In+iJn)=∮∣z∣=1zn(z−a)(1−az)dzi=1i∮∣z∣=1zn(z−a)(1−az)dz=2πi⋅1i⋅Res(zn(z−a)(1−az),a)=2πlim⁡z→azn1−az=2πan1−a2\begin{aligned} (I_n + i J_n) &= \oint_{|z|=1} \dfrac{z^n}{(z-a)(1-az)} \dfrac{dz}{i}\\ &= \dfrac{1}{i} \oint_{|z|=1} \dfrac{z^{n}}{(z-a)(1-az)} dz\\ &= 2\pi i \cdot \dfrac{1}{i} \cdot \mathrm{Res} \left( \dfrac{z^{n}}{(z-a)(1-az)} , a \right)\\ &= 2\pi \lim_{z \to a} \dfrac{z^{n}}{1-az}\\ &= \dfrac{2\pi a^n}{1-a^2} \end{aligned} と計算できる。よって In=∫02πcos⁡nθ1−2acos⁡θ+a2dθ=2πan1−a2 I_n = \int_0^{2\pi} \dfrac{\cos n\theta}{1-2 a \cos \theta +a^2} d\theta = \dfrac{2\pi a^n}{1-a^2} が得られる。

必要な留数計算が減り,一気に簡単になったとわかるでしょう。三角関数を含む積分ではオイラーの公式を積極的に使うようにしましょう。実際に次の問題もオイラーの公式を用いて計算することになります。

3.三角関数と有理式の積

有理式を R(z)=P(z)/Q(z)R(z) = P(z) / Q(z) と表したとき

  • deg⁡P+1≦deg⁡Q\deg P + 1 \leqq \deg Q
  • Q(z)=0Q(z) = 0 の実数解は重解ではない

を満たすときは次のように積分をしましょう。

  1. 次のように閉曲線を取る
    • 実部が正の特異点をすべて含むように,十分大きな RR を半径とする半円を取る
    • 実軸上の特異点の周囲 ε\varepsilon を避けるように −R-R から RR へ直線を取る
    • 特異点周りには(上で取った直線とつながるように)小さな半円を取る
  2. 閉曲線周りの積分値を留数定理から求める
  3. R→∞R \to \infty,ε→0\varepsilon \to 0 の評価をする

計算では,オイラーの公式により eize^{iz} の積分をすることを考えます。

解

十分小さい δ>0\delta > 0 を取る。RR を δ\delta より大になるように取る。

  • CR={z∈C∣∣z∣=R,Imz≧0}C_R = \{ z \in \mathbb{C} \mid |z| = R , \mathrm{Im} z \geqq 0 \}(下図青色)
  • LR,δL_{R,\delta}:実軸上の区間 [−R,−δ][-R,-\delta] と [δ,R][\delta , R] の和(下図赤色)
  • Lε,δL_{\varepsilon,\delta}:実軸上の区間 [−δ,−ε][-\delta,-\varepsilon] と [ε,δ][\varepsilon , \delta] の和(下図橙色)
  • Cε={z∈C∣∣z∣=ε,Imz≦0}C_{\varepsilon} = \{ z \in \mathbb{C} \mid |z| = \varepsilon , \mathrm{Im} z \leqq 0 \}(下図水色)

と定め,CC を上記4曲線の和とする。このとき CC は単純閉曲線である。

pic3-1

オイラーの公式より eizz=cos⁡zz+isin⁡zz \dfrac{e^{iz}}{z} = \dfrac{\cos z}{z} + i \dfrac{\sin z}{z} である。eizz\dfrac{e^{iz}}{z} を CC で積分して,その虚部を調べればよい。

eizz\dfrac{e^{iz}}{z} の孤立特異点は z=0z=0 で位数は1であるため,留数定理から ∮Ceizzdz=2πi  Res(eizz,0)=2πilim⁡z→0eiz=2πi\begin{aligned} \oint_C \dfrac{e^{iz}}{z} dz &= 2\pi i \; \mathrm{Res} \left(\dfrac{e^{iz}}{z} , 0\right)\\ &= 2\pi i\lim_{z \to 0} e^{iz}\\ &= 2\pi i \end{aligned} である。

CRC_R での積分を評価する。

lim⁡R→∞∣∫CReizzdz∣≦lim⁡R→∞∫CR∣eizzdz∣=lim⁡R→∞∫0π∣eiReiθReiθiReiθ∣dθ=lim⁡R→∞∫0π∣eiRcos⁡θe−Rsin⁡θ∣dθ=lim⁡R→∞∫0πe−Rsin⁡θdθ\begin{aligned} \lim_{R \to \infty} \left| \int_{C_R}\dfrac{e^{iz}}{z} dz \right| &\leqq \lim_{R \to \infty} \int_{C_R} \left| \dfrac{e^{iz}}{z} dz \right|\\ &= \lim_{R \to \infty} \int_0^{\pi} \left| \dfrac{e^{iRe^{i\theta}}}{Re^{i\theta}} iR e^{i\theta} \right| d\theta\\ &= \lim_{R \to \infty} \int_0^{\pi} \left| e^{iR\cos \theta} e^{-R\sin \theta} \right| d\theta\\ &= \lim_{R \to \infty} \int_0^{\pi} e^{-R\sin \theta} d\theta \end{aligned} と計算される。

最後の積分を詳しく評価する。0≦θ≦π20 \leqq \theta \leqq \dfrac{\pi}{2} で sin⁡θ≧2πθ\sin \theta \geqq \dfrac{2}{\pi}\theta である。

これを用いて式を変形すると ∫0πe−Rsin⁡θdθ=2∫0π2e−Rsin⁡θdθ≦2∫0π2e−R2πθdθ=πR−πRe−R→0(R→∞)\begin{aligned} \int_0^{\pi} e^{-R\sin \theta} d\theta &= 2 \int_0^{\frac{\pi}{2}} e^{-R\sin \theta} d\theta\\ &\leqq 2 \int_0^{\frac{\pi}{2}} e^{-R \frac{2}{\pi}\theta} d\theta\\ &= \dfrac{\pi}{R} - \dfrac{\pi}{R} e^{-R} \to 0 \quad (R \to \infty) \end{aligned} であるため, lim⁡R→∞∫CReizzdz=0 \lim_{R \to \infty} \int_{C_R} \dfrac{e^{iz}}{z} dz = 0 である。

ε→0\varepsilon \to 0 のときの CεC_{\varepsilon} 上の積分を評価する。eizz\dfrac{e^{iz}}{z} は z=0z=0 で1位の極を持つ。ゆえに正則関数 h(z)h(z) を用いると,z=0z=0 近傍で eizz\dfrac{e^{iz}}{z} をローラン展開したものを c−1z+h(z)\dfrac{c_{-1}}{z} + h(z) と表すことができる。h(z)h(z) は正則であるため ∣h(z)∣|h(z)| は CεC_{\varepsilon} 上で最大値を取る。

こうして lim⁡ε→0∣∫Cεh(z)  dz∣≦lim⁡ε→0∫Cε∣h(z)∣  dz=lim⁡ε→0πεmax⁡z∈Cε(∣h(z)∣)=0\begin{aligned} \lim_{\varepsilon \to 0} \left| \int_{C_{\varepsilon}} h(z) \; dz \right| &\leqq \lim_{\varepsilon \to 0} \int_{C_{\varepsilon}} |h(z)| \; dz\\ &= \lim_{\varepsilon \to 0} \pi \varepsilon \max_{z \in C_{\varepsilon}} (|h(z)|)\\ &= 0 \end{aligned} である。一方 ∫Cεc−1zdz=∫π2πc−1εeiθεieiθdθ=πic−1=πi  Res(eizz,0)=πilim⁡z→0zeizz=πi\begin{aligned} \int_{C_{\varepsilon}} \dfrac{c_{-1}}{z} dz &= \int_{\pi}^{2\pi} \dfrac{c_{-1}}{\varepsilon e^{i\theta}} \varepsilon i e^{i\theta} d\theta\\ &= \pi i c_{-1}\\ &= \pi i \; \mathrm{Res} \left( \dfrac{e^{iz}}{z} , 0 \right)\\ &= \pi i \lim_{z \to 0} z \dfrac{e^{iz}}{z} = \pi i \end{aligned} である。よって lim⁡ε→0∫Cεeizzdz=πi \lim_{\varepsilon \to 0} \int_{C_{\varepsilon}} \dfrac{e^{iz}}{z} dz = \pi i

LR,δ,Lε,δL_{R,\delta},L_{\varepsilon,\delta} での積分は,LR,δL_{R,\delta} は ε\varepsilonに,Lε,δL_{\varepsilon,\delta} は RR に依らないことから各々の極限を考えれば良く, lim⁡R→∞lim⁡ε→0∫LR,δ+Lε,δeizzdz=lim⁡R→∞∫LR,δeizzdz+lim⁡ε→0∫Lε,δeizzdz=p.v.∫−∞∞eixxdx\begin{aligned} \lim_{R \to \infty} \lim_{\varepsilon \to 0} \int_{L_{R,\delta}+L_{\varepsilon,\delta}} \dfrac{e^{iz}}{z} dz &= \lim_{R \to \infty} \int_{L_{R,\delta}} \dfrac{e^{iz}}{z} dz + \lim_{\varepsilon \to 0} \int_{L_{\varepsilon,\delta}} \dfrac{e^{iz}}{z} dz\\ &= p.v. \int_{-\infty}^{\infty} \dfrac{e^{ix}}{x} dx \end{aligned} となる。

なお p.v.∫p.v. \int に関しては,後述する主値積分を参照してください。広義積分の一種と考えて問題ありません。

以上をまとめると lim⁡R→∞lim⁡ε→0∮Ceizzdz=lim⁡R→∞lim⁡ε→0∫LR,δ+Lε,δeizzdz+lim⁡R→∞∫CReizzdz+lim⁡ε→0∫Cεeizzdz=p.v.∫−∞∞eixxdx+πi\begin{aligned} &\lim_{R \to \infty} \lim_{\varepsilon \to 0} \oint_{C} \dfrac{e^{iz}}{z} dz\\ & = \lim_{R \to \infty} \lim_{\varepsilon \to 0} \int_{L_{R,\delta}+L_{\varepsilon,\delta}} \dfrac{e^{iz}}{z} dz + \lim_{R \to \infty} \int_{C_R} \dfrac{e^{iz}}{z} dz + \lim_{\varepsilon \to 0} \int_{C_{\varepsilon}} \dfrac{e^{iz}}{z} dz\\ &= p.v. \int_{-\infty}^{\infty} \dfrac{e^{ix}}{x} dx + \pi i \end{aligned} 及び lim⁡R→∞lim⁡ε→0∮Ceizzdz=2πi \lim_{R \to \infty} \lim_{\varepsilon \to 0} \oint_{C} \dfrac{e^{iz}}{z} dz = 2\pi i であるため, p.v.∫−∞∞eixxdx=πi p.v. \int_{-\infty}^{\infty} \dfrac{e^{ix}}{x} dx = \pi i である。 p.v.∫−∞∞eixxdx=p.v.∫−∞∞cos⁡xxdx+p.v.∫−∞∞isin⁡xxdx p.v. \int_{-\infty}^{\infty} \dfrac{e^{ix}}{x} dx = p.v. \int_{-\infty}^{\infty} \dfrac{\cos x}{x} dx + p.v. \int_{-\infty}^{\infty} i\dfrac{\sin x}{x} dx より,実部と虚部を比較して ∫0∞sin⁡xxdx=12p.v.∫−∞∞sin⁡xxdx=π2 \int_0^{\infty} \dfrac{\sin x}{x} dx = \dfrac{1}{2} p.v. \int_{-\infty}^{\infty} \dfrac{\sin x}{x} dx = \dfrac{\pi}{2} を得る。

有理式の仮定がどのように効いているかを確認します。

CRC_R 上の積分の評価のとき,有理式の次数と eize^{iz} を被積分関数に含むことが効いて,R→0R \to 0 で積分が 00 に収束することがわかります。

また Q(z)Q(z) が重解を持たない仮定により CεC_{\varepsilon} 上の積分の評価で見通しよく計算できます。

sinc 関数

sin⁡zz\dfrac{\sin z}{z} は sinc 関数(シンク関数) と呼ばれ,信号処理などの分野で非常に重要な関数となります。

sinc 関数は z=0z=0 を孤立特異点に持ちますが,この特異点は可除特異点になるので z=0z=0 でも正則であるとみなすことができます。よって本質的には,sinc 関数の実数全体の積分は主値積分と扱わなくても良いです。

一方 eizz\dfrac{e^{iz}}{z} は z=0z=0 で極を持つので主値積分として考える必要があります。

発展例題

応用例題

次の積分を求めよ。 ∫−∞∞(sin⁡xx)3dx \int_{-\infty}^{\infty} \left( \dfrac{\sin x}{x} \right)^3 dx

三倍角の公式を使えば sin⁡zz3\dfrac{\sin z}{z^3} の形に帰着できます。

しかし下手をすると CεC_{\varepsilon} での評価が難しくなるかもしれません。

証明

三倍角の公式により sin⁡3x=3sin⁡x−sin⁡3x4\sin^3 x = \dfrac{3\sin x - \sin 3x}{4} であるため ∫−∞∞3eiz−e3iz4z3dz \int_{-\infty}^{\infty} \dfrac{3e^{iz} - e^{3iz}}{4z^3} dz を計算すればよい。

積分路を前の問題3と同じように取ると,同様の議論から lim⁡R→∞∫CR3eiz−e3iz4z3dz=0 \lim_{R \to \infty} \int_{C_R} \dfrac{3e^{iz} - e^{3iz}}{4z^3} dz = 0 であるため ∫−∞∞(sin⁡xx)3dx=−lim⁡ε→0∫Cε3eiz−e3iz4z3dz \int_{-\infty}^{\infty}\left( \dfrac{\sin x}{x} \right)^3 dx = - \lim_{\varepsilon \to 0} \int_{C_{\varepsilon}} \dfrac{3e^{iz} - e^{3iz}}{4z^3} dz を得る。

一方 z=0z=0 近傍でローラン展開すると,被積分関数は 12z3+34z+i−1316z+⋯ \dfrac{1}{2z^3} + \dfrac{3}{4z} + i - \dfrac{13}{16} z + \cdots である。00 次以上の項は前の例題同様 ε→0\varepsilon \to 0 で消える。

−3-3 次の項は ∫Cεdzz3=∫−π0ε−2ie−2iθdθ=[e−2iθ−2ε2]−π0=0\begin{aligned} &\int_{C_{\varepsilon}} \dfrac{dz}{z^3}\\ &= \int_{-\pi}^0 \varepsilon^{-2} i e^{-2i\theta} d\theta\\ &= \Big[ \dfrac{e^{-2i\theta}}{-2\varepsilon^2} \Big]_{-\pi}^0\\ &= 0 \end{aligned} となり消える。

こうして lim⁡ε→0∫Cε3eiz−e3iz4z3dz=−34πi \lim_{\varepsilon \to 0} \int_{C_{\varepsilon}} \dfrac{3e^{iz} - e^{3iz}}{4z^3} dz = -\dfrac{3}{4}\pi i であるため ∫−∞∞(sin⁡xx)3dx=34π \int_{-\infty}^{\infty}\left( \dfrac{\sin x}{x} \right)^3 dx = \dfrac{3}{4} \pi となる。

∫−∞∞sin⁡3xx3dx=27∫−∞∞sin⁡xx3dx \int_{-\infty}^{\infty} \dfrac{\sin 3x}{x^3} dx = 27 \int_{-\infty}^{\infty} \dfrac{\sin x}{x^3} dx であることから sin⁡xx3\dfrac{\sin x}{x^3} の積分だけを計算すればいいと考えてしまうと,CεC_{\varepsilon} で積分が消えません。

実はわざわざ分けない方が上手くいくという計算例でした。

また,わざわざ三倍角の公式を使わなくとも,そのままの形からローラン展開を用いて計算することもできます。

4.フレネル積分(ガウス積分を用いた計算)

∫−∞∞sin⁡x2  dx=π2\displaystyle \int_{-\infty}^{\infty} \sin x^2 \; dx = \sqrt{\dfrac{\pi}{2}} となります。これをフレネル積分といいます。オイラーの公式とガウス積分を用いて計算できます。

シンプルですが思いつきにくい積分経路を取ります。

解

sin⁡x2\sin x^2 は偶関数であるため ∫−∞∞sin⁡x2  dx=2∫0∞sin⁡x2  dx \int_{-\infty}^{\infty} \sin x^2 \; dx = 2 \int_0^{\infty} \sin x^2 \; dx である。

eix2=cos⁡x2+isin⁡x2e^{ix^2} = \cos x^2 + i \sin x^2 であるため,∫0∞eix2  dx\displaystyle \int_{0}^{\infty} e^{ix^2} \; dx を計算し虚部を見ればよい。

次のように経路を定める。

  • C1C_1:00 と RR を結ぶ線分(下図赤色)
  • C2={z∣∣z∣=R,0<arg⁡z<π4}C_2 = \{ z \mid |z| = R , 0 < \arg z < \frac{\pi}{4} \}(下図青色)
  • C3C_3:R(12+i2)R\left(\dfrac{1}{\sqrt{2}} + \dfrac{i}{\sqrt{2}}\right) と 00 を結ぶ線分(下図緑色)

C=C1+C2+C3C = C_1 + C_2 + C_3 とおく。向きは反時計回りにつける。

pic4-1

CC が囲う領域で eix2e^{ix^2} は正則なので ∮Ceiz2dz=0 \oint_C e^{iz^2} dz = 0 である。

計算すると lim⁡R→∞∫C1eiz2dz=∫0∞eix2dxlim⁡R→∞∫C3eiz2  dz=lim⁡R→∞∫R0eix2(12+i2)2(12+i2)dx=(12+i2)∫∞0e−x2x=−(12+i2)π2\begin{aligned} \lim_{ R \to \infty} \int_{C_1} e^{iz^2} dz &= \int_0^{\infty} e^{ix^2} dx\\ \lim_{R \to \infty} \int_{C_3} e^{iz^2} \; dz &= \lim_{R \to \infty} \int_{R}^0 e^{ix^2 (\frac{1}{\sqrt{2}} + \frac{i}{\sqrt{2}} )^2} \left( \dfrac{1}{\sqrt{2}} + \dfrac{i}{\sqrt{2}} \right) dx \\ &= \left( \dfrac{1}{\sqrt{2}} + \dfrac{i}{\sqrt{2}} \right) \int_{\infty}^0 e^{-x^2} x\\ &= -\left( \dfrac{1}{\sqrt{2}} + \dfrac{i}{\sqrt{2}} \right) \dfrac{\sqrt{\pi}}{2} \end{aligned} となる。一方 lim⁡R→∞∣∫C2eiz2dz∣=lim⁡R→∞∣∫0π4eiR2e2iθRieiθdθ∣=lim⁡R→∞∣∫0π4eiR2(cos⁡2θ+isin⁡2θ)Rieiθdθ∣≦lim⁡R→∞∫0π4∣eiR2(cos⁡2θ+isin⁡2θ)Rieiθ∣dθ=lim⁡R→∞∫0π4Re−R2sin⁡2θdθ=0\begin{aligned} &\lim_{R \to \infty} \left| \int_{C_2} e^{iz^2} dz \right| \\&= \lim_{R \to \infty} \left| \int_{0}^{\frac{\pi}{4}} e^{i R^2 e^{2i\theta}} R i e^{i\theta} d\theta \right|\\ &= \lim_{R \to \infty} \left| \int_{0}^{\frac{\pi}{4}} e^{i R^2 (\cos 2\theta + i \sin 2\theta)} R i e^{i\theta} d\theta \right|\\ &\leqq \lim_{R \to \infty} \int_{0}^{\frac{\pi}{4}} \left| e^{i R^2 (\cos 2\theta + i \sin 2\theta)} R i e^{i\theta} \right|d\theta\\ &= \lim_{R \to \infty} \int_0^{\frac{\pi}{4}} Re^{- R^2 \sin 2\theta} d\theta\\ &= 0 \end{aligned} である。こうして lim⁡R→∞∫Ceiz2dz=∫0∞eiz2−(π8+iπ8)=0 \lim_{R \to \infty} \int_C e^{iz^2} dz = \int_0^{\infty} e^{iz^2} - \left(\sqrt{\dfrac{\pi}{8}} + i \sqrt{\dfrac{\pi}{8}} \right) = 0 であるため ∫−∞∞sin⁡x2  dx=π2 \int_{-\infty}^{\infty} \sin x^2 \; dx = \sqrt{\dfrac{\pi}{2}} を得る。同時に ∫−∞∞cos⁡x2  dx=π2 \int_{-\infty}^{\infty} \cos x^2 \; dx = \sqrt{\dfrac{\pi}{2}} となる。

5.ガウス積分を用いた計算2

オイラーの公式とガウス積分を組み合わせて計算をしましょう。積分経路が少々思いつきにくいかもしれません。

まず,オイラーの公式を用いて e−x2−iax=e−x2cos⁡ax−ie−x2sin⁡ax e^{-x^2-iax} = e^{-x^2} \cos ax - i e^{-x^2} \sin ax という式を得ます。こうして得られた左辺は e−x2−iax=e−a24e−x2−iax+a24=e−14a2e−(x+a2i)2 e^{-x^2-iax} = e^{-\frac{a^2}{4}} e^{-x^2-iax+\frac{a^2}{4}} = e^{-\frac{1}{4} a^2} e^{-\left( x+\frac{a}{2} i \right)^2} と変形されます。よって e−(x+a2i)2e^{-\left( x+\frac{a}{2} i \right)^2} を積分した値がわかれば良いですね。xx を R\mathbb{R} で動かすとき,積分値は e−z2e^{-z^2} を Im(z)=a2\mathrm{Im} (z) = \dfrac{a}{2} なる直線上で積分した値となります。これら観察から解答に示す経路を取ることになります。

解

CC を −R<Re(z)<R,0<Im(z)<a2-R < \mathrm{Re} (z) < R , 0 < \mathrm{Im} (z) < \dfrac{a}{2} で定まる領域を囲う閉曲線とする。

pic5-1

e−z2e^{-z^2} は C\mathbb{C} 全体で正則であるため ∮Ce−z2dz=0 \oint_C e^{-z^2} dz = 0 である。

Im(z)=0\mathrm{Im} (z) = 0 上の積分を考えると,ガウス積分より lim⁡R→∞∫{Im(z)=0}∩Ce−z2dz=∫−∞∞e−x2dx=π \lim_{R \to \infty} \int_{\{\mathrm{Im} (z) = 0 \} \cap C} e^{-z^2}dz = \int_{-\infty}^{\infty} e^{-x^2}dx = \sqrt{\pi} となる。

Im(z)=a2\mathrm{Im} (z) = \dfrac{a}{2} 上の積分を考える。 ∫{Im(z)=a2}∩Ce−z2dz=−∫−RRe−(x+a2i)2dx=−∫−RRe−x2e−iaxea24dx=−ea24∫−RRe−x2(cos⁡ax−isin⁡ax)  dx=−ea24∫−RRe−x2cos⁡ax  dx\begin{aligned} \int_{\left\{\mathrm{Im} (z) = \frac{a}{2} \right\} \cap C} e^{-z^2} dz &= -\int_{-R}^R e^{-(x+\frac{a}{2}i)^2} dx\\ &= -\int_{-R}^R e^{-x^2} e^{-iax} e^{\frac{a^2}{4}} dx\\ &= -e^{\frac{a^2}{4}} \int_{-R}^R e^{-x^2} (\cos ax - i \sin ax) \; dx\\ &= -e^{\frac{a^2}{4}} \int_{-R}^R e^{-x^2} \cos ax \; dx \end{aligned} である。途中 e−x2sin⁡axe^{-x^2} \sin ax が奇関数であることを用いた。

よって lim⁡R→∞∫{Im(z)=a2}∩Ce−z2dz=−ea24∫−∞∞e−x2cos⁡ax  dx \lim_{R \to \infty} \int_{\left\{\mathrm{Im} (z) = \frac{a}{2} \right\} \cap C} e^{-z^2} dz = -e^{\frac{a^2}{4}} \int_{-\infty}^{\infty} e^{-x^2} \cos ax \; dx である。

Re(z)=R\mathrm{Re} (z) = R 上の積分を考える。 ∫{Re(z)=R}∩Ce−z2dz=∫0ae−(R+iy)2dy=∫0ae−R2e−2iRy−y2dy=e−R2∫0ae−2iRy−y2dy\begin{aligned} \int_{\{\mathrm{Re} (z) = R\} \cap C} e^{-z^2} dz &= \int_0^a e^{-(R+iy)^2} dy\\ &=\int_0^a e^{-R^2} e^{-2iRy-y^2} dy\\ &= e^{-R^2} \int_0^a e^{-2iRy-y^2} dy \end{aligned} となる。e−2iRy−y2e^{-2iRy-y^2} は積分域で有界であるため ∫0ae−2iRy−y2dy=M<+∞\displaystyle \int_0^a e^{-2iRy-y^2} dy = M < + \infty である。よって lim⁡R→∞∫{Re(z)=R}∩Ce−z2dz=lim⁡R→∞e−R2M=0 \lim_{R \to \infty} \int_{\{\mathrm{Re} (z) = R\} \cap C} e^{-z^2} dz = \lim_{R \to \infty} e^{-R^2} M =0 である。同様に Re(z)=−R\mathrm{Re} (z) = -R 上の積分も R→∞R \to \infty で 00 に収束する。

向きに注意して以上をまとめると lim⁡R→∞∮Ce−z2dz=π−ea24∫−∞∞e−x2cos⁡ax  dx=0\begin{aligned} \lim_{R\to \infty} \oint_C e^{-z^2} dz &= \sqrt{\pi} - e^{\frac{a^2}{4}} \int_{-\infty}^{\infty} e^{-x^2} \cos ax \; dx\\ &=0 \end{aligned} であるため ∫0∞e−x2cos⁡ax  dx=π2e−a24 \int_{0}^{\infty} e^{-x^2} \cos ax \; dx = \dfrac{\sqrt{\pi}}{2} e^{-\frac{a^2}{4}} となる。

補足:主値積分

実数全体の積分を考える際,実軸上に極があるケースを考えるときがあります。その時,そのまま積分することができないため,次のような方法で表記をします。

複素関数 R(z)R(z) は多項式 P,QP,Q を用いて R(z)=P(z)Q(z)R(z) = \dfrac{P(z)}{Q(z)} と表されるとしましょう。

R(z)R(z) が r∈Rr \in \mathbb{R} で高々1位の極を持つとします。この場合,rr を含む実積分は定義されないため,代わりに十分小さい δ>0\delta > 0 により [r−δ,r+δ][r - \delta , r+\delta] 上の積分値を考えます。

Q(r)=0Q(r) = 0 なる rr まわりの主値積分を,固定した小さな δ>0\delta > 0 に大して p.v.∫b−δb+δR(z)  dz=lim⁡ε→0(∫r−δr−εR(z)  dz+∫r+δr+εR(z)  dz) p.v. \int_{b - \delta}^{b+\delta} R(z) \; dz = \lim_{\varepsilon \to 0} \left( \int_{r-\delta}^{r-\varepsilon} R(z) \; dz + \int_{r+\delta}^{r+\varepsilon} R(z) \; dz \right) と定めます。

実数全体での主値積分を,Q(z)=0Q(z)=0 の各実数解 rr のまわりでは主値積分を考えるものとして p.v.∫−∞∞R(z)  dz p.v. \int_{-\infty}^{\infty} R(z) \; dz と書きます。

今後も様々な実積分を取り扱っていきます。