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

問題

次の積分を求めよ。

  1. ∫−∞∞dxx2+1\displaystyle \int_{-\infty}^{\infty} \dfrac{dx}{x^2+1}
  2. p.v.∫−∞∞dxx4−1\displaystyle p.v. \int_{-\infty}^{\infty} \dfrac{dx}{x^4-1}
  3. ∫0∞dxx3+1\displaystyle \int_0^{\infty} \dfrac{dx}{x^3+1}

なお p.v.∫p.v. \int は主値積分を意味する。(補足を参照)

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

有理式 R(z)R(z) は R(z)=P(z)Q(z)R(z) = \dfrac{P(z)}{Q(z)} と2つの多項式 P(z),Q(z)P(z) , Q(z) を用いて表されるとしましょう。

  • deg⁡P+2≦deg⁡Q\deg P + 2 \leqq \deg Q
  • R(z)R(z) の極の位数が高々 11 となる

とき,半径が十分大きな円を含む積分経路を考えると,(−∞,∞)(-\infty,\infty) での積分を計算できます。早速見ていきましょう。

1.分母の多項式に実数解がないパターン

留数定理 の記事の具体例と同じパターンです。

(arctan⁡x)′=1x2+1(\arctan x)' = \dfrac{1}{x^2+1} を用いると簡単に計算できますが,留数定理を用いても積分できます。

解

実数 R>0R > 0 を十分大きくとる。

C1C_1 を (−R)(-R) から (R)(R) を結ぶ直線,C2C_2 を ∣z∣=R|z|=R の Im(z)≧0\mathrm{Im} (z) \geqq 0 の部分とする。なお,向きは反時計回りにつける。C=C1+C2C = C_1 + C_2 とする。

pic1

1z2+1\dfrac{1}{z^2+1} の特異点は ±i\pm i で,各点での位数は 11 である。

RR を十分大きく取っているため,ii は CC で囲まれた領域に入る。

よって留数定理から

∮Cdzz2+1=2πi  Res(1z2+1,i) \oint_{C} \dfrac{dz}{z^2+1} = 2 \pi i \; \mathrm{Res} \left( \dfrac{1}{z^2+1} , i \right)

である。

Res(1z2+1,i)=lim⁡z→i(z−i)⋅1z2+1=lim⁡z→i1z+i=12i\begin{aligned} \mathrm{Res} \left( \dfrac{1}{z^2+1} , i \right) &= \lim_{z \to i} (z-i) \cdot \dfrac{1}{z^2+1}\\ &= \lim_{z\to i} \dfrac{1}{z+i}\\ &=\dfrac{1}{2i} \end{aligned}

より, ∮Cdzz2+1=2πi  Res(1z2+1,i)=2πi12i=π\begin{aligned} \oint_{C} \dfrac{dz}{z^2+1} &= 2 \pi i \; \mathrm{Res} \left( \dfrac{1}{z^2+1} , i \right)\\ &= 2 \pi i \dfrac{1}{2i}\\ &= \pi \end{aligned}

上の結果は RR に寄らないため,R→∞R \to \infty としても積分値は変わらない。すなわち lim⁡R→∞∮Cdzz2+1=π \lim_{R \to \infty} \oint_{C} \dfrac{dz}{z^2+1} = \pi となる。

次に C1,C2C_1,C_2 での積分を計算する。

lim⁡R→∞∫C1dzz2+1=∫−∞∞dxx2+1 \lim_{R \to \infty} \int_{C_1} \dfrac{dz}{z^2+1} = \int_{-\infty}^{\infty} \dfrac{dx}{x^2+1} は求めたい値になる。一方 C2C_2 については,z=Reiθz=Re^{i\theta} と置換して三角不等式 ∣a+b∣≧∣a∣−∣b∣|a+b|\geqq |a|-|b| を使うと, ∣lim⁡R→∞∫C2dzz2+1∣≦lim⁡R→∞∫0π∣Rieiθ∣dθ∣R2e2iθ+1∣≦lim⁡R→∞∫0πRR2−1  dθ=0\begin{aligned} \left| \lim_{R \to \infty} \int_{C_2} \dfrac{dz}{z^2+1} \right| &\leqq \lim_{R \to \infty} \int_{0}^{\pi} \dfrac{|Rie^{i\theta}|d\theta}{|R^2e^{2i\theta}+1|}\\ &\leqq \lim_{R \to \infty} \int_0^{\pi} \dfrac{R}{R^2-1} \; d\theta\\ &= 0 \end{aligned} より lim⁡R→∞∫C2dzz2+1=0\displaystyle \lim_{R \to \infty} \int_{C_2} \dfrac{dz}{z^2+1} = 0 である。

以上より ∫−∞∞dxx2+1=lim⁡R→∞∫C1dzz2+1=lim⁡R→∞∫C1dzz2+1+∫C2dzz2+1=lim⁡R→∞∮Cdzz2+1=π\begin{aligned} \int_{-\infty}^{\infty} \dfrac{dx}{x^2+1} &= \lim_{R \to \infty} \int_{C_1} \dfrac{dz}{z^2+1}\\ &= \lim_{R \to \infty} \int_{C_1} \dfrac{dz}{z^2+1} + \int_{C_2} \dfrac{dz}{z^2+1}\\ &= \lim_{R \to \infty} \oint_{C} \dfrac{dz}{z^2+1}\\ &= \pi \end{aligned} である。

2.分母の多項式に実数解があるパターン

有理式 R(z)R(z) は R(z)=P(z)Q(z)R(z) = \dfrac{P(z)}{Q(z)} と2つの多項式 P(z),Q(z)P(z) , Q(z) を用いて表されるとしましょう。

deg⁡P+2≦deg⁡Q\deg P + 2 \leqq \deg Q かつ R(z)R(z) の極の位数が高々 11 となるとき,次に示す方法で (−∞,∞)(-\infty,\infty) での積分ができます。

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

R>0R>0 を,Δ(0,R)\Delta(0,R) が Δ(±1,δ)\Delta (\pm 1 , \delta) を含むように大きくとる。

  • CR={z∣∣z∣=R,Im(z)≧0}C_R = \{ z \mid |z| = R , \mathrm{Im} (z) \geqq 0 \}
  • LR,δL_{R,\delta}:実軸上の区間 [−R,−1−δ][-R,-1-\delta],[−1+δ,1−δ][-1+\delta,1-\delta],[1+δ,R][1+\delta,R] の和
  • Lε,δL_{\varepsilon,\delta}:実軸上の区間 [−1−δ,−1−ε][-1-\delta,-1-\varepsilon],[−1+ε,−1+δ][-1+\varepsilon,-1+\delta],[1−δ,1−ε][1-\delta,1-\varepsilon],[1+ε,1+δ][1+\varepsilon,1+\delta] の和
  • cε(1)={z∣∣z−1∣=ε,Im(z)≦0}c_{\varepsilon} (1) = \{ z \mid |z-1| = \varepsilon , \mathrm{Im} (z) \leqq 0 \}
  • cε(−1)={z∣∣z+1∣=ε,Im(z)≦0}c_{\varepsilon} (-1) = \{ z \mid |z+1| = \varepsilon , \mathrm{Im} (z) \leqq 0 \}

これらの和を CC とおく。このとき向きは反時計回りにとる。

pic2

CC は ±1,i\pm 1,i を含むため,留数定理から ∮Cdzz4−1=2πi(Res(1z4−1,1)+Res(1z4−1,−1)+Res(1z4−1,i))=2πi(1(z+1)(z2+1)∣z=1+1(z−1)(z2+1)∣z=−1+1(z+i)(z2−1)∣z=i)=2πi(14−14+i4)=−π2\begin{aligned} &\oint_C \dfrac{dz}{z^4-1}\\ &= 2\pi i \left( \mathrm{Res} \left( \dfrac{1}{z^4-1} , 1\right) + \mathrm{Res} \left( \dfrac{1}{z^4-1} , -1\right) + \mathrm{Res} \left( \dfrac{1}{z^4-1} , i\right) \right)\\ &= 2\pi i \left( \left. \dfrac{1}{(z+1)(z^2+1)} \right|_{z=1} + \left. \dfrac{1}{(z-1)(z^2+1)} \right|_{z=-1} +\left. \dfrac{1}{(z+i)(z^2-1)} \right|_{z=i} \right)\\ &= 2 \pi i \left( \dfrac{1}{4} - \dfrac{1}{4} + \dfrac{i}{4} \right)\\ &= -\dfrac{\pi}{2} \end{aligned} である。

CRC_R での積分は例1と同様に計算できる: lim⁡R→∞∣∫CRdzz4−1∣≦lim⁡R→∞∫0πRR4−1  dθ=lim⁡R→∞πRR4−1=0\begin{aligned} \lim_{R \to \infty} \left| \int_{C_R} \dfrac{dz}{z^4-1} \right| &\leqq \lim_{R \to \infty} \int_{0}^{\pi} \dfrac{R}{R^4-1} \; d\theta\\ &= \lim_{R \to \infty} \dfrac{\pi R}{R^4-1} = 0 \end{aligned}

残りの経路での積分を考える。

LR,δL_{R,\delta} と Lε,δL_{\varepsilon,\delta} が求めたい積分である:

lim⁡R→∞lim⁡ε→0∫LR,δ+Lε,δdzz4−1=lim⁡R→∞∫LR,δdzz4−1+lim⁡ε→0∫Lε,δdzz4−1=p.v.∫−∞∞dzz4−1\begin{aligned} &\lim_{R \to \infty} \lim_{\varepsilon \to 0} \int_{L_{R,\delta}+L_{\varepsilon,\delta}} \dfrac{dz}{z^4-1}\\ &=\lim_{R \to \infty} \int_{L_{R,\delta}} \dfrac{dz}{z^4-1} +\lim_{\varepsilon \to 0} \int_{L_{\varepsilon,\delta}} \dfrac{dz}{z^4-1}\\ &= p.v. \int_{-\infty}^{\infty} \dfrac{dz}{z^4-1} \end{aligned}

残りは,cε(1)c_{\varepsilon} (1) と cε(−1)c_{\varepsilon}(-1) での積分である。cε(1)c_{\varepsilon} (1) について考える。

1z4−1\dfrac{1}{z^4-1} の 11 周りのローラン展開を c−1z−1+h(z)\dfrac{c_{-1}}{z-1} + h(z)(h(z)h(z) は次数が非負の部分)ととる。このとき h(z)h(z) は cε(1)c_{\varepsilon} (1) 上で正則である。よって ∣h(z)∣|h(z)| は cε(1)c_{\varepsilon} (1) 上で最大値を取る。

こうして lim⁡ε→0∣∫cε(1)h(z)  dz∣≦lim⁡ε→0∫cε(1)∣h(z)∣  dz=lim⁡ε→0πεmax⁡z∈cε(1)(∣h(z)∣)=0\begin{aligned} \lim_{\varepsilon \to 0} \left| \int_{c_{\varepsilon} (1)} h(z) \; dz \right| &\leqq \lim_{\varepsilon \to 0} \int_{c_{\varepsilon} (1)} |h(z)| \; dz\\ &= \lim_{\varepsilon \to 0} \pi \varepsilon \max_{z \in c_{\varepsilon} (1)} (|h(z)|)\\ &= 0 \end{aligned} と評価される。一方 ∫cε(1)c−1z−1  dz=∫π2πic−1  dz=πic−1=πi  Res(1z4−1,1)=πi4\begin{aligned} \int_{c_{\varepsilon} (1)} \dfrac{c_{-1}}{z-1} \; dz &= \int_{\pi}^{2\pi} i c_{-1} \; dz\\ &= \pi i c_{-1}\\ &= \pi i \; \mathrm{Res} \left( \dfrac{1}{z^4-1} , 1 \right)\\ &= \dfrac{\pi i}{4} \end{aligned} である。こうして lim⁡ε→0∫cε(1)dzz4−1=πi4 \lim_{\varepsilon \to 0} \int_{c_{\varepsilon} (1)} \dfrac{dz}{z^4-1} = \dfrac{\pi i}{4} である。同様に lim⁡ε→0∫cε(−1)dzz4−1=−πi4 \lim_{\varepsilon \to 0} \int_{c_{\varepsilon} (-1)} \dfrac{dz}{z^4-1} = -\dfrac{\pi i}{4} となる。ゆえに lim⁡ε→0∫cε(−1)+cε(1)dzz4−1=πi4−πi4=0 \lim_{\varepsilon \to 0} \int_{c_{\varepsilon} (-1) + c_{\varepsilon} (1)} \dfrac{dz}{z^4-1} = \dfrac{\pi i}{4}-\dfrac{\pi i}{4} = 0

以上を合わせると p.v.∫−∞∞dzz4−1=lim⁡R→∞∫LR,δdzz4−1+lim⁡ε→0∫cε(−1)+cε(1)+Lε,δdzz4−1=lim⁡R→∞lim⁡ε→0∮Cdzz4−1=−π2\begin{aligned} &p.v. \int_{-\infty}^{\infty} \dfrac{dz}{z^4-1}\\ &= \lim_{R \to \infty} \int_{L_{R,\delta}} \dfrac{dz}{z^4-1} + \lim_{\varepsilon \to 0} \int_{c_{\varepsilon} (-1) + c_{\varepsilon} (1) + L_{\varepsilon,\delta}} \dfrac{dz}{z^4-1}\\ &= \lim_{R \to \infty} \lim_{\varepsilon \to 0} \oint_{C} \dfrac{dz}{z^4-1}\\ &= - \dfrac{\pi}{2} \end{aligned}

が得られる。

3.特殊な経路を使うパターン

一般に m,nm,n が整数で 0<m<n0 < m < n となるとき ∫0∞xm−1xn+1dx=πnsin⁡mπn \int_0^{\infty} \dfrac{x^{m-1}}{x^n+1} dx = \dfrac{\pi}{n \sin \frac{m\pi}{n}} と計算されます。

積分経路は 扇形 になるようにとります。

解

CC を偏角が 00 から 2nπ\dfrac{2}{n}\pi で半径が RR の扇形の周とする。向きは反時計回りにつける。また下図のように CC のうち実軸上の部分を C1C_1,円弧の部分を C2C_2,00 と Re2πniRe^{\frac{2\pi}{n}i} を結ぶ直線を C3C_3 とする。

pic3

zm−1zn+1\dfrac{z^{m-1}}{z^n+1} は CC 内に1位の特異点 z=eπniz = e^{\frac{\pi}{n}i} を持つ。よって留数定理より ∮Czm−1zn+1dz=2πi  Res(zm−1zn+1,eπni) \oint_C \dfrac{z^{m-1}}{z^n+1} dz = 2\pi i \; \mathrm{Res} \left( \dfrac{z^{m-1}}{z^n+1} , e^{\frac{\pi}{n}i} \right) である。右辺を計算しよう。 Res(zm−1zn+1,eπni)=lim⁡z→eπni(z−eπni)zm−1zn+1=lim⁡z→eπnizm−1(z−e3πni)(z−e5πni)⋯(z−e(2n−1)πni)=(eπni)m−1(eπni−e3πni)(eπni−e5πni)⋯(eπni−e(2n−1)πni)=e(m−1)πnieπni⋅(n−1)(1−e2πni)(1−e4πni)⋯(1−e2(n−1)πni)\begin{aligned} &\mathrm{Res} \left( \dfrac{z^{m-1}}{z^n+1} , e^{\frac{\pi}{n}i} \right)\\ &= \lim_{z \to e^{\frac{\pi}{n}i}} (z-e^{\frac{\pi}{n}i}) \dfrac{z^{m-1}}{z^n+1}\\ &= \lim_{z \to e^{\frac{\pi}{n}i}} \dfrac{z^{m-1}}{(z-e^{\frac{3\pi}{n}i})(z-e^{\frac{5\pi}{n}i}) \cdots (z-e^{\frac{(2n-1)\pi}{n}i})}\\ &= \dfrac{(e^{\frac{\pi}{n}i})^{m-1}}{(e^{\frac{\pi}{n}i}-e^{\frac{3\pi}{n}i})(e^{\frac{\pi}{n}i}-e^{\frac{5\pi}{n}i}) \cdots (e^{\frac{\pi}{n}i}-e^{\frac{(2n-1)\pi}{n}i})}\\ &=\dfrac{e^{\frac{(m-1)\pi}{n}i}}{e^{\frac{\pi}{n}i \cdot (n-1)} (1-e^{\frac{2\pi}{n}i})(1-e^{\frac{4\pi}{n}i}) \cdots (1-e^{\frac{2(n-1)\pi}{n}i})} \end{aligned} と計算されるが,e2kπni  (k=1,⋯n)e^{\frac{2k\pi}{n}i}\; (k = 1, \cdots n) は zn=1z^n=1 の解であるため, (zn−1)=(z−e2πni)⋯(z−e2nπni)=(z−1)(z−e2πni)⋯(z−e2(n−1)πni)\begin{aligned} (z^n-1) &= (z-e^{\frac{2\pi}{n}i}) \cdots (z- e^{\frac{2n\pi}{n}i})\\ &= (z-1) (z-e^{\frac{2\pi}{n}i}) \cdots (z-e^{\frac{2(n-1)\pi}{n}i}) \end{aligned} となる。(zn−1)=(z−1)(zn−1+zn−2+⋯1)(z^n-1) = (z-1)(z^{n-1} + z^{n-2} + \cdots 1) と合わせて zn−1+zn−2+⋯1=(z−e2πni)(z−e4πni)⋯(z−e2(n−1)πni) z^{n-1} + z^{n-2} + \cdots 1 = (z-e^{\frac{2\pi}{n}i})(z-e^{\frac{4\pi}{n}i}) \cdots (z-e^{\frac{2(n-1)\pi}{n}i}) である。こうして Res(zm−1zn+1,eπni)=e(m−1)πnieπni⋅(n−1)(1−e2πni)(1−e4πni)⋯(1−e2(n−1)πni)=e(m−1)πnie(n−1)πnin=−emπnin\begin{aligned} &\mathrm{Res} \left( \dfrac{z^{m-1}}{z^n+1} , e^{\frac{\pi}{n}i} \right)\\ &= \dfrac{e^{\frac{(m-1)\pi}{n}i}}{e^{\frac{\pi}{n}i \cdot (n-1)} (1-e^{\frac{2\pi}{n}i})(1-e^{\frac{4\pi}{n}i}) \cdots (1-e^{\frac{2(n-1)\pi}{n}i})}\\ &= \dfrac{e^{\frac{(m-1)\pi}{n}i}}{e^{\frac{(n-1)\pi}{n}i}n }\\ &=- \dfrac{e^{\frac{m\pi}{n}i}}{n} \end{aligned} を得る。

C1C_1 上の積分は lim⁡R→∞∫C1zm−1zn+1dz=∫0∞xm−1xn+1dx \lim_{R \to \infty} \int_{C_1} \dfrac{z^{m-1}}{z^n+1} dz = \int_0^{\infty} \dfrac{x^{m-1}}{x^n+1} dx となり求めたい積分。

C2C_2 上の積分はこれまでの例と同様に lim⁡R→∞∣∫C2zm−1zn+1dz∣≦lim⁡R→∞∫02nπRmRn−1dθ=0\begin{aligned} \lim_{R \to \infty} \left| \int_{C_2} \dfrac{z^{m-1}}{z^n+1} dz \right| &\leqq \lim_{R \to \infty} \int_0^{\frac{2}{n}\pi}\dfrac{R^{m}}{R^n - 1} d\theta\\ &= 0 \end{aligned} と計算できる。

C3C_3 上の積分を考える。z=e2πnixz= e^{\frac{2\pi}{n}i}x と置換すれば次のように計算される。 lim⁡R→∞∫C3zm−1zn+1dz=lim⁡R→∞∫R0e2π(m−1)nixm−1e2πni⋅nxn+1e2πnidx=−e2πmnilim⁡R→∞∫0Rxm−1xn+1dx=−e2πmni∫0∞xm−1xn+1dx\begin{aligned} \lim_{R \to \infty} \int_{C_3} \dfrac{z^{m-1}}{z^n+1} dz &= \lim_{R \to \infty} \int_{R}^0 \dfrac{e^{\frac{2\pi(m-1)}{n}i} x^{m-1}}{e^{\frac{2\pi}{n} i \cdot n} x^n + 1} e^{\frac{2\pi}{n}i} dx\\ &= -e^{\frac{2\pi m}{n}i} \lim_{R \to \infty} \int_0^R \dfrac{x^{m-1}}{x^n+1} dx\\ &= -e^{\frac{2\pi m}{n}i} \int_0^{\infty} \dfrac{x^{m-1}}{x^n+1} dx \end{aligned}

これらをまとめると lim⁡R→∞∮Czm−1zn+1dz=(1−e2πmni)∫0∞xm−1xn+1dx=−2πi⋅emπnin\begin{aligned} \lim_{R \to \infty} \oint_C \dfrac{z^{m-1}}{z^n+1} dz &= (1-e^{\frac{2\pi m}{n}i}) \int_0^{\infty} \dfrac{x^{m-1}}{x^n+1} dx\\ &=- 2\pi i \cdot \dfrac{e^{\frac{m\pi}{n}i}}{n} \end{aligned} を得る。こうして ∫0∞xm−1xn+1dx=2πie2mπni−1⋅emπnin=2iemπni−e−mπni⋅πn=πnsin⁡mπn\begin{aligned} \int_0^{\infty} \dfrac{x^{m-1}}{x^n+1} dx &= \dfrac{2\pi i}{e^{\frac{2m\pi}{n}i} - 1} \cdot \dfrac{e^{\frac{m\pi}{n}i}}{n}\\ &= \dfrac{2i}{e^{\frac{m\pi}{n}i}-e^{-\frac{m\pi}{n}i}} \cdot \dfrac{\pi}{n} \\ &= \dfrac{\pi}{n \sin \frac{m\pi}{n}} \end{aligned} となる。

補足:主値積分

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

複素関数 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.∫r−δr+δR(z)  dz=lim⁡ε→0(∫r−δr−εR(z)  dz+∫r+εr+δR(z)  dz) p.v. \int_{r - \delta}^{r+\delta} R(z) \; dz = \lim_{\varepsilon \to 0} \left( \int_{r-\delta}^{r-\varepsilon} R(z) \; dz + \int_{r+\varepsilon}^{r+\delta} R(z) \; dz \right) と定めます。

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

3番目の扇形を用いた積分は美しいですね。