この章で扱う留数定理は、この種の積分を次の手続きに置き換えます。被積分関数を複素平面に延長し、積分路を閉曲線になるように閉じ、閉曲線の内部にある特異点で「留数」という一つの数を計算して足す。それだけです。∫ − ∞ ∞ d x x 4 + 1 \int_{-\infty}^{\infty}\frac{dx}{x^4+1} ∫ − ∞ ∞ x 4 + 1 d x は、後で見るように二つの留数の足し算だけで π / 2 \pi/\sqrt{2} π / 2 と分かります。
記号 :領域とは連結開集合のことです。D ( a , r ) = { z ∈ C : ∣ z − a ∣ < r } D(a,r) = \{z \in \mathbb{C} : |z-a| < r\} D ( a , r ) = { z ∈ C : ∣ z − a ∣ < r } 、D ′ ( a , r ) = { z : 0 < ∣ z − a ∣ < r } D'(a,r) = \{z : 0 < |z-a| < r\} D ′ ( a , r ) = { z : 0 < ∣ z − a ∣ < r } (中心を抜いた円板、除外近傍)と書きます。曲線はすべて区分的に C 1 C^1 C 1 級とし、∫ γ f ( z ) d z = ∫ α β f ( γ ( t ) ) γ ′ ( t ) d t \int_\gamma f(z)\,dz = \int_\alpha^\beta f(\gamma(t))\gamma'(t)\,dt ∫ γ f ( z ) d z = ∫ α β f ( γ ( t )) γ ′ ( t ) d t で定義します。円周 ∣ z − a ∣ = r |z-a| = r ∣ z − a ∣ = r には断らない限り正の向き(反時計回り)を与え、閉曲線に沿う積分は ∮ \oint ∮ と書きます。
長さと最大値による評価(ML 不等式) :曲線 γ \gamma γ の長さを L ( γ ) L(\gamma) L ( γ ) とすると
∣ ∫ γ f ( z ) d z ∣ ≤ ∫ γ ∣ f ( z ) ∣ ∣ d z ∣ ≤ L ( γ ) ⋅ max z ∈ γ ∣ f ( z ) ∣ \left|\int_\gamma f(z)\,dz\right| \le \int_\gamma |f(z)|\,|dz| \le L(\gamma)\cdot\max_{z \in \gamma}|f(z)| ∫ γ f ( z ) d z ≤ ∫ γ ∣ f ( z ) ∣ ∣ d z ∣ ≤ L ( γ ) ⋅ z ∈ γ max ∣ f ( z ) ∣ が成り立ちます(Lemma 2.4[コーシーの積分定理と積分公式] )。この章の評価はほとんどすべてこの不等式に帰着します。
ローラン展開 :正則関数の強力な性質 の Theorem 6.1[正則関数の強力な性質] で示したとおり、f f f が円環 r < ∣ z − a ∣ < R r < |z-a| < R r < ∣ z − a ∣ < R で正則ならば、そこで
f ( z ) = ∑ n = − ∞ ∞ c n ( z − a ) n , c n = 1 2 π i ∮ ∣ ζ − a ∣ = ρ f ( ζ ) ( ζ − a ) n + 1 d ζ ( r < ρ < R ) f(z) = \sum_{n=-\infty}^{\infty}c_n (z-a)^n,
\qquad
c_n = \frac{1}{2\pi i}\oint_{|\zeta-a|=\rho}\frac{f(\zeta)}{(\zeta-a)^{n+1}}\,d\zeta
\quad (r < \rho < R) f ( z ) = n = − ∞ ∑ ∞ c n ( z − a ) n , c n = 2 π i 1 ∮ ∣ ζ − a ∣ = ρ ( ζ − a ) n + 1 f ( ζ ) d ζ ( r < ρ < R ) と一意的に展開でき、この級数は円環の任意のコンパクト部分集合上で一様収束します。係数 c n c_n c n は ρ \rho ρ の取り方によりません。
多重連結領域に対するコーシーの積分定理 :C C C を正の向きの単純閉曲線、γ 1 , … , γ n \gamma_1,\dots,\gamma_n γ 1 , … , γ n を C C C の内部にあり互いに他の外部にある正の向きの単純閉曲線とします。f f f が「C C C の内側かつすべての γ j \gamma_j γ j の外側」という閉領域を含むある開集合で正則ならば
∮ C f ( z ) d z = ∑ j = 1 n ∮ γ j f ( z ) d z \oint_C f(z)\,dz = \sum_{j=1}^{n}\oint_{\gamma_j} f(z)\,dz ∮ C f ( z ) d z = j = 1 ∑ n ∮ γ j f ( z ) d z が成り立ちます。証明は前章の Theorem 5.1[コーシーの積分定理と積分公式] と同じ要領で、C C C と各 γ j \gamma_j γ j を細い「橋」(切断線)でつないで単連結な部分に分け、各橋を往復する積分が符号違いで相殺することを使います。留数定理はこの定理の言い換えにすぎません。
Definition 3.1 (孤立特異点 )
a ∈ C a \in \mathbb{C} a ∈ C とする。ある R > 0 R > 0 R > 0 が存在して f f f が除外近傍 D ′ ( a , R ) = { z : 0 < ∣ z − a ∣ < R } D'(a,R) = \{z : 0 < |z-a| < R\} D ′ ( a , R ) = { z : 0 < ∣ z − a ∣ < R } で正則であるとき、a a a を f f f の孤立特異点 と呼ぶ。f f f は a a a 自身では定義されていなくてもよいし、定義されていても正則でなくてよい。
孤立特異点のまわりでは f f f は円環 0 < ∣ z − a ∣ < R 0 < |z-a| < R 0 < ∣ z − a ∣ < R で正則ですから、§2 のローラン展開が使えます。展開の負べきの部分が特異点の「悪さ」を測る量です。
Definition 3.2 (主要部と特異点の三分類 )
a a a を f f f の孤立特異点とし、D ′ ( a , R ) D'(a,R) D ′ ( a , R ) でのローラン展開を f ( z ) = ∑ n = − ∞ ∞ c n ( z − a ) n f(z) = \sum_{n=-\infty}^{\infty}c_n(z-a)^n f ( z ) = ∑ n = − ∞ ∞ c n ( z − a ) n とする。負べきの部分
P a ( z ) = ∑ n ≤ − 1 c n ( z − a ) n P_a(z) = \sum_{n \le -1} c_n (z-a)^n P a ( z ) = n ≤ − 1 ∑ c n ( z − a ) n を a a a における f f f の主要部 と呼ぶ。主要部の形により、a a a を次のように分類する。
c n = 0 c_n = 0 c n = 0 (すべての n ≤ − 1 n \le -1 n ≤ − 1 )、すなわち主要部が 0 0 0 のとき、a a a を除去可能特異点 という。
ある正の整数 m m m について c − m ≠ 0 c_{-m} \ne 0 c − m = 0 かつ c n = 0 c_n = 0 c n = 0 (すべての n < − m n < -m n < − m )のとき、a a a を m m m 位の極 という。m = 1 m = 1 m = 1 のときは単純極 ともいう。
c n ≠ 0 c_n \ne 0 c n = 0 となる n ≤ − 1 n \le -1 n ≤ − 1 が無限個あるとき、a a a を真性特異点 という。
ローラン展開は一意なので、この分類は展開の取り方によらず定まる。
flowchart TD
A["孤立特異点 a のまわりでローラン展開する"] --> B["主要部(負べきの項の全体)を見る"]
B --> C["負べきの項が無い"]
B --> D["負べきが有限個で、最低次が -m"]
B --> E["負べきが無限個ある"]
C --> C2["除去可能特異点:a まで正則に延長できる。近くでの有界性で判定できる"]
D --> D2["m 位の極:z が a に近づくと絶対値が無限大に発散する"]
E --> E2["真性特異点:どんな近傍の像も複素平面で稠密になる"] 孤立特異点の三分類。主要部の項数を数えるだけで型が決まり、それぞれに極限による判定条件が対応します。 分類の定義はローラン係数を全部調べよと言っているように見えますが、実際にはもっと安上がりな判定法があります。まず除去可能特異点です。
Theorem 3.3 (リーマンの除去可能特異点定理 )
a a a を f f f の孤立特異点とし、f f f は D ′ ( a , R ) D'(a,R) D ′ ( a , R ) で正則とする。次の 4 条件は同値である。
(i) a a a は f f f の除去可能特異点である。
(ii) f ( a ) f(a) f ( a ) の値を適当に定めると、f f f は D ( a , R ) D(a,R) D ( a , R ) 全体で正則になる。
(iii) ある 0 < r ≤ R 0 < r \le R 0 < r ≤ R と M > 0 M > 0 M > 0 があって、0 < ∣ z − a ∣ < r 0 < |z-a| < r 0 < ∣ z − a ∣ < r で ∣ f ( z ) ∣ ≤ M |f(z)| \le M ∣ f ( z ) ∣ ≤ M が成り立つ。
(iv) lim z → a ( z − a ) f ( z ) = 0 \lim_{z \to a}(z-a)f(z) = 0 lim z → a ( z − a ) f ( z ) = 0 。
Proof(Theorem 3.3) (i) ⇒ \Rightarrow ⇒ (ii) :主要部が 0 0 0 なので D ′ ( a , R ) D'(a,R) D ′ ( a , R ) で f ( z ) = ∑ n ≥ 0 c n ( z − a ) n f(z) = \sum_{n \ge 0}c_n(z-a)^n f ( z ) = ∑ n ≥ 0 c n ( z − a ) n です。右辺のべき級数は D ′ ( a , R ) D'(a,R) D ′ ( a , R ) で収束するので収束半径は R R R 以上であり、D ( a , R ) D(a,R) D ( a , R ) 全体で正則関数を定めます。f ( a ) : = c 0 f(a) := c_0 f ( a ) := c 0 と定めれば、f f f はこのべき級数と D ( a , R ) D(a,R) D ( a , R ) 全体で一致し、正則になります。
(ii) ⇒ \Rightarrow ⇒ (iii) :正則関数は連続なので、r = R / 2 r = R/2 r = R /2 とすればコンパクト集合 D ( a , r ) ‾ \overline{D(a,r)} D ( a , r ) 上で ∣ f ∣ |f| ∣ f ∣ は最大値 M M M を取ります。
(iii) ⇒ \Rightarrow ⇒ (iv) :∣ ( z − a ) f ( z ) ∣ ≤ M ∣ z − a ∣ → 0 |(z-a)f(z)| \le M|z-a| \to 0 ∣ ( z − a ) f ( z ) ∣ ≤ M ∣ z − a ∣ → 0 (z → a z \to a z → a )です。
(iv) ⇒ \Rightarrow ⇒ (i) :0 < ρ < R 0 < \rho < R 0 < ρ < R に対し η ( ρ ) : = max ∣ z − a ∣ = ρ ∣ ( z − a ) f ( z ) ∣ \eta(\rho) := \max_{|z-a| = \rho}|(z-a)f(z)| η ( ρ ) := max ∣ z − a ∣ = ρ ∣ ( z − a ) f ( z ) ∣ とおくと、仮定より η ( ρ ) → 0 \eta(\rho) \to 0 η ( ρ ) → 0 (ρ → 0 \rho \to 0 ρ → 0 )です。§2 のローラン係数の式と ML 不等式から、n ≤ − 1 n \le -1 n ≤ − 1 に対して
∣ c n ∣ ≤ 1 2 π ⋅ 2 π ρ ⋅ max ∣ z − a ∣ = ρ ∣ f ( z ) ∣ ρ n + 1 = ρ − n max ∣ z − a ∣ = ρ ∣ f ( z ) ∣ = ρ − n − 1 η ( ρ ) |c_n| \le \frac{1}{2\pi}\cdot 2\pi\rho \cdot \frac{\max_{|z-a|=\rho}|f(z)|}{\rho^{\,n+1}}
= \rho^{-n}\max_{|z-a|=\rho}|f(z)|
= \rho^{-n-1}\,\eta(\rho) ∣ c n ∣ ≤ 2 π 1 ⋅ 2 π ρ ⋅ ρ n + 1 max ∣ z − a ∣ = ρ ∣ f ( z ) ∣ = ρ − n ∣ z − a ∣ = ρ max ∣ f ( z ) ∣ = ρ − n − 1 η ( ρ ) が成り立ちます。ここで n ≤ − 1 n \le -1 n ≤ − 1 より − n − 1 ≥ 0 -n-1 \ge 0 − n − 1 ≥ 0 なので、ρ < 1 \rho < 1 ρ < 1 のとき ρ − n − 1 ≤ 1 \rho^{-n-1} \le 1 ρ − n − 1 ≤ 1 です。よって ∣ c n ∣ ≤ η ( ρ ) → 0 |c_n| \le \eta(\rho) \to 0 ∣ c n ∣ ≤ η ( ρ ) → 0 となり、c n c_n c n は ρ \rho ρ によらない定数ですから c n = 0 c_n = 0 c n = 0 です。これがすべての n ≤ − 1 n \le -1 n ≤ − 1 で成り立つので主要部は 0 0 0 です。
∎ 条件 (iii) が重要です。「近くで有界」という、ローラン係数を一つも計算しなくてよい条件だけで、特異点が消せると分かります。次は極です。
Proposition 3.4 (極の特徴づけ )
a a a を f f f の孤立特異点、m m m を正の整数とする。
(1) 次の 3 条件は同値である。
(i) a a a は f f f の m m m 位の極である。
(ii) 極限 lim z → a ( z − a ) m f ( z ) \lim_{z \to a}(z-a)^m f(z) lim z → a ( z − a ) m f ( z ) が存在し、その値は 0 0 0 でない。
(iii) ある r > 0 r > 0 r > 0 と D ( a , r ) D(a,r) D ( a , r ) で正則な関数 g g g で g ( a ) ≠ 0 g(a) \ne 0 g ( a ) = 0 を満たすものがあり、0 < ∣ z − a ∣ < r 0 < |z-a| < r 0 < ∣ z − a ∣ < r で f ( z ) = g ( z ) ( z − a ) m f(z) = \dfrac{g(z)}{(z-a)^m} f ( z ) = ( z − a ) m g ( z ) と書ける。
(2) a a a が(何らかの位数の)極であるための必要十分条件は lim z → a ∣ f ( z ) ∣ = ∞ \lim_{z \to a}|f(z)| = \infty lim z → a ∣ f ( z ) ∣ = ∞ である。
Proof(Proposition 3.4) (1) (i) ⇒ \Rightarrow ⇒ (iii) :m m m 位の極なら D ′ ( a , R ) D'(a,R) D ′ ( a , R ) で f ( z ) = ∑ n ≥ − m c n ( z − a ) n f(z) = \sum_{n \ge -m}c_n(z-a)^n f ( z ) = ∑ n ≥ − m c n ( z − a ) n かつ c − m ≠ 0 c_{-m} \ne 0 c − m = 0 です。両辺に ( z − a ) m (z-a)^m ( z − a ) m を掛けると
( z − a ) m f ( z ) = ∑ n ≥ − m c n ( z − a ) n + m = ∑ k ≥ 0 c k − m ( z − a ) k (z-a)^m f(z) = \sum_{n \ge -m} c_n (z-a)^{n+m} = \sum_{k \ge 0} c_{k-m}(z-a)^k ( z − a ) m f ( z ) = n ≥ − m ∑ c n ( z − a ) n + m = k ≥ 0 ∑ c k − m ( z − a ) k となり、右辺は D ′ ( a , R ) D'(a,R) D ′ ( a , R ) で収束するべき級数ですから収束半径は R R R 以上で、D ( a , R ) D(a,R) D ( a , R ) で正則関数 g g g を定めます。g ( a ) = c − m ≠ 0 g(a) = c_{-m} \ne 0 g ( a ) = c − m = 0 であり、0 < ∣ z − a ∣ < R 0 < |z-a| < R 0 < ∣ z − a ∣ < R で f ( z ) = g ( z ) / ( z − a ) m f(z) = g(z)/(z-a)^m f ( z ) = g ( z ) / ( z − a ) m です。
(iii) ⇒ \Rightarrow ⇒ (ii) :( z − a ) m f ( z ) = g ( z ) → g ( a ) ≠ 0 (z-a)^m f(z) = g(z) \to g(a) \ne 0 ( z − a ) m f ( z ) = g ( z ) → g ( a ) = 0 (z → a z \to a z → a 、g g g の連続性)です。
(ii) ⇒ \Rightarrow ⇒ (i) :h ( z ) : = ( z − a ) m f ( z ) h(z) := (z-a)^m f(z) h ( z ) := ( z − a ) m f ( z ) とおくと、h h h は D ′ ( a , R ) D'(a,R) D ′ ( a , R ) で正則で z → a z \to a z → a のとき有限の極限 L ≠ 0 L \ne 0 L = 0 を持つので、a a a の近くで有界です。Theorem 3.3 の (iii) ⇒ \Rightarrow ⇒ (i),(ii) より a a a は h h h の除去可能特異点で、h ( a ) = L h(a) = L h ( a ) = L と定めれば h h h は D ( a , R ) D(a,R) D ( a , R ) で正則になります。h h h をテイラー展開して h ( z ) = ∑ k ≥ 0 b k ( z − a ) k h(z) = \sum_{k \ge 0}b_k(z-a)^k h ( z ) = ∑ k ≥ 0 b k ( z − a ) k (b 0 = L ≠ 0 b_0 = L \ne 0 b 0 = L = 0 )とすると、0 < ∣ z − a ∣ < R 0 < |z-a| < R 0 < ∣ z − a ∣ < R で
f ( z ) = h ( z ) ( z − a ) m = ∑ k ≥ 0 b k ( z − a ) k − m f(z) = \frac{h(z)}{(z-a)^m} = \sum_{k \ge 0} b_k (z-a)^{k-m} f ( z ) = ( z − a ) m h ( z ) = k ≥ 0 ∑ b k ( z − a ) k − m となります。これは f f f のローラン展開そのもの(一意性より)で、最低次の項は b 0 ( z − a ) − m b_0(z-a)^{-m} b 0 ( z − a ) − m 、b 0 ≠ 0 b_0 \ne 0 b 0 = 0 ですから a a a は m m m 位の極です。
(2) 十分性 :a a a が m m m 位の極なら (iii) より f = g / ( z − a ) m f = g/(z-a)^m f = g / ( z − a ) m 、g ( a ) ≠ 0 g(a) \ne 0 g ( a ) = 0 なので、∣ f ( z ) ∣ = ∣ g ( z ) ∣ ∣ z − a ∣ − m → ∣ g ( a ) ∣ ⋅ ∞ = ∞ |f(z)| = |g(z)|\,|z-a|^{-m} \to |g(a)| \cdot \infty = \infty ∣ f ( z ) ∣ = ∣ g ( z ) ∣ ∣ z − a ∣ − m → ∣ g ( a ) ∣ ⋅ ∞ = ∞ です。
必要性 :∣ f ( z ) ∣ → ∞ |f(z)| \to \infty ∣ f ( z ) ∣ → ∞ とします。ある r > 0 r > 0 r > 0 があって 0 < ∣ z − a ∣ < r 0 < |z-a| < r 0 < ∣ z − a ∣ < r で ∣ f ( z ) ∣ ≥ 1 |f(z)| \ge 1 ∣ f ( z ) ∣ ≥ 1 、特に f ( z ) ≠ 0 f(z) \ne 0 f ( z ) = 0 です。そこで u : = 1 / f u := 1/f u := 1/ f とおくと u u u は D ′ ( a , r ) D'(a,r) D ′ ( a , r ) で正則、∣ u ∣ ≤ 1 |u| \le 1 ∣ u ∣ ≤ 1 と有界なので、Theorem 3.3 より a a a は u u u の除去可能特異点で、u u u は D ( a , r ) D(a,r) D ( a , r ) で正則に延長されます。∣ f ∣ → ∞ |f| \to \infty ∣ f ∣ → ∞ より u ( a ) = 0 u(a) = 0 u ( a ) = 0 です。u u u は D ′ ( a , r ) D'(a,r) D ′ ( a , r ) で 0 0 0 にならないので u ≢ 0 u \not\equiv 0 u ≡ 0 であり、Lemma 7.2[正則関数の強力な性質] で見たとおり、恒等的に 0 0 0 でない正則関数の零点は有限の位数を持ちます。その位数を m m m とすれば u ( z ) = ( z − a ) m v ( z ) u(z) = (z-a)^m v(z) u ( z ) = ( z − a ) m v ( z ) 、v v v は D ( a , r ′ ) D(a,r') D ( a , r ′ ) で正則、v ( a ) ≠ 0 v(a) \ne 0 v ( a ) = 0 と書けます。よって 0 < ∣ z − a ∣ < r ′ 0 < |z-a| < r' 0 < ∣ z − a ∣ < r ′ で f ( z ) = 1 / v ( z ) ( z − a ) m f(z) = \dfrac{1/v(z)}{(z-a)^m} f ( z ) = ( z − a ) m 1/ v ( z ) となり、1 / v 1/v 1/ v は a a a の近くで正則(v ( a ) ≠ 0 v(a)\ne 0 v ( a ) = 0 と連続性より近傍で v ≠ 0 v \ne 0 v = 0 )かつ a a a で 0 0 0 でないので、(1)(iii) から a a a は m m m 位の極です。
∎ 残るのが真性特異点です。Theorem 3.3 と Proposition 3.4 により、a a a が真性特異点であることは「z → a z \to a z → a で f ( z ) f(z) f ( z ) が有限の極限も持たず、∣ f ( z ) ∣ → ∞ |f(z)| \to \infty ∣ f ( z ) ∣ → ∞ でもない」ことと同値です。極限が存在しないだけでなく、値の散らばり方は極端です。
Theorem 3.5 (カゾラーティ・ワイエルシュトラスの定理 )
a a a を f f f の真性特異点とし、f f f は D ′ ( a , R ) D'(a,R) D ′ ( a , R ) で正則とする。このとき任意の 0 < r ≤ R 0 < r \le R 0 < r ≤ R に対して像 f ( D ′ ( a , r ) ) f(D'(a,r)) f ( D ′ ( a , r )) は C \mathbb{C} C で稠密である。すなわち、任意の w ∈ C w \in \mathbb{C} w ∈ C と任意の ε > 0 \varepsilon > 0 ε > 0 に対し、0 < ∣ z − a ∣ < r 0 < |z-a| < r 0 < ∣ z − a ∣ < r かつ ∣ f ( z ) − w ∣ < ε |f(z) - w| < \varepsilon ∣ f ( z ) − w ∣ < ε を満たす z z z が存在する。
Proof(Theorem 3.5) 稠密でないと仮定します。すると、ある w ∈ C w \in \mathbb{C} w ∈ C と δ > 0 \delta > 0 δ > 0 があって、0 < ∣ z − a ∣ < r 0 < |z-a| < r 0 < ∣ z − a ∣ < r のすべての z z z で ∣ f ( z ) − w ∣ ≥ δ |f(z) - w| \ge \delta ∣ f ( z ) − w ∣ ≥ δ となります。そこで
g ( z ) : = 1 f ( z ) − w ( 0 < ∣ z − a ∣ < r ) g(z) := \frac{1}{f(z) - w} \qquad (0 < |z-a| < r) g ( z ) := f ( z ) − w 1 ( 0 < ∣ z − a ∣ < r ) とおくと、分母は 0 0 0 にならないので g g g は D ′ ( a , r ) D'(a,r) D ′ ( a , r ) で正則で、∣ g ∣ ≤ 1 / δ |g| \le 1/\delta ∣ g ∣ ≤ 1/ δ と有界です。Theorem 3.3 より a a a は g g g の除去可能特異点で、g g g は D ( a , r ) D(a,r) D ( a , r ) で正則に延長されます。
g ( a ) ≠ 0 g(a) \ne 0 g ( a ) = 0 の場合:f ( z ) = w + 1 / g ( z ) f(z) = w + 1/g(z) f ( z ) = w + 1/ g ( z ) は a a a の近くで有界なので、再び Theorem 3.3 により a a a は f f f の除去可能特異点です。
g ( a ) = 0 g(a) = 0 g ( a ) = 0 の場合:g g g は D ′ ( a , r ) D'(a,r) D ′ ( a , r ) で 0 0 0 にならないので g ≢ 0 g \not\equiv 0 g ≡ 0 であり、a a a は g g g の有限位数 m m m の零点です。Proposition 3.4 の (1)(iii) の形に書けば f ( z ) − w f(z) - w f ( z ) − w は a a a で m m m 位の極を持ち、定数 w w w を足しても主要部は変わらないので f f f も a a a で m m m 位の極を持ちます。
どちらの場合も a a a が真性特異点であることに反します。
∎ Example 3.6 (特異点の分類の実例 )
(a) f ( z ) = sin z z f(z) = \dfrac{\sin z}{z} f ( z ) = z sin z の z = 0 z = 0 z = 0 。sin z = z − z 3 3 ! + z 5 5 ! − ⋯ \sin z = z - \frac{z^3}{3!} + \frac{z^5}{5!} - \cdots sin z = z − 3 ! z 3 + 5 ! z 5 − ⋯ を z z z で割ると f ( z ) = 1 − z 2 6 + z 4 120 − ⋯ f(z) = 1 - \frac{z^2}{6} + \frac{z^4}{120} - \cdots f ( z ) = 1 − 6 z 2 + 120 z 4 − ⋯ となり、負べきの項がありません。よって 0 0 0 は除去可能特異点で、f ( 0 ) : = 1 f(0) := 1 f ( 0 ) := 1 と定めれば整関数になります。
(b) f ( z ) = 1 z 2 ( z − 1 ) f(z) = \dfrac{1}{z^2(z-1)} f ( z ) = z 2 ( z − 1 ) 1 。z = 0 z = 0 z = 0 では z 2 f ( z ) = 1 z − 1 → − 1 ≠ 0 z^2 f(z) = \frac{1}{z-1} \to -1 \ne 0 z 2 f ( z ) = z − 1 1 → − 1 = 0 なので Proposition 3.4 より 2 2 2 位の極、z = 1 z = 1 z = 1 では ( z − 1 ) f ( z ) = 1 z 2 → 1 ≠ 0 (z-1)f(z) = \frac{1}{z^2} \to 1 \ne 0 ( z − 1 ) f ( z ) = z 2 1 → 1 = 0 なので単純極です。
(c) f ( z ) = 1 sin z f(z) = \dfrac{1}{\sin z} f ( z ) = sin z 1 の z = k π z = k\pi z = k π (k ∈ Z k \in \mathbb{Z} k ∈ Z )。sin ( k π ) = 0 \sin(k\pi) = 0 sin ( k π ) = 0 かつ d d z sin z ∣ z = k π = cos ( k π ) = ( − 1 ) k ≠ 0 \frac{d}{dz}\sin z\big|_{z=k\pi} = \cos(k\pi) = (-1)^k \ne 0 d z d sin z z = k π = cos ( k π ) = ( − 1 ) k = 0 なので k π k\pi k π は sin z \sin z sin z の 1 1 1 位の零点、したがって
( z − k π ) f ( z ) = z − k π sin z − sin ( k π ) ⟶ 1 cos ( k π ) = ( − 1 ) k ≠ 0 (z-k\pi)f(z) = \frac{z-k\pi}{\sin z - \sin(k\pi)} \longrightarrow \frac{1}{\cos(k\pi)} = (-1)^k \ne 0 ( z − k π ) f ( z ) = sin z − sin ( k π ) z − k π ⟶ cos ( k π ) 1 = ( − 1 ) k = 0 (微分係数の定義そのもの)となり、k π k\pi k π は単純極です。
(d) f ( z ) = e 1 / z f(z) = e^{1/z} f ( z ) = e 1/ z の z = 0 z = 0 z = 0 。e w = ∑ n ≥ 0 w n / n ! e^{w} = \sum_{n\ge0}w^n/n! e w = ∑ n ≥ 0 w n / n ! に w = 1 / z w = 1/z w = 1/ z を代入して
e 1 / z = ∑ n ≥ 0 1 n ! z − n = 1 + 1 z + 1 2 ! z 2 + ⋯ ( z ≠ 0 ) e^{1/z} = \sum_{n \ge 0}\frac{1}{n!}z^{-n} = 1 + \frac{1}{z} + \frac{1}{2!z^2} + \cdots
\qquad (z \ne 0) e 1/ z = n ≥ 0 ∑ n ! 1 z − n = 1 + z 1 + 2 ! z 2 1 + ⋯ ( z = 0 ) が 0 < ∣ z ∣ 0 < |z| 0 < ∣ z ∣ で収束します。負べきの項が無限個あるので 0 0 0 は真性特異点です。実際、値の散らばり方を直接確かめられます。w ≠ 0 w \ne 0 w = 0 を任意に取り w = ∣ w ∣ e i θ w = |w|e^{i\theta} w = ∣ w ∣ e i θ と書くと、z k = 1 log ∣ w ∣ + i ( θ + 2 π k ) z_k = \dfrac{1}{\log|w| + i(\theta + 2\pi k)} z k = log ∣ w ∣ + i ( θ + 2 π k ) 1 (k ∈ N k \in \mathbb{N} k ∈ N )は k → ∞ k \to \infty k → ∞ で 0 0 0 に収束し、e 1 / z k = w e^{1/z_k} = w e 1/ z k = w です。つまり 0 0 0 のどんな近傍でも、f f f は 0 0 0 以外のすべての複素数を無限回取ります。Theorem 3.5 の「稠密」よりずっと強い状況です。
Definition 4.1 (留数 )
a a a を f f f の孤立特異点とし、D ′ ( a , R ) D'(a,R) D ′ ( a , R ) でのローラン展開を f ( z ) = ∑ n = − ∞ ∞ c n ( z − a ) n f(z) = \sum_{n=-\infty}^{\infty}c_n(z-a)^n f ( z ) = ∑ n = − ∞ ∞ c n ( z − a ) n とする。係数 c − 1 c_{-1} c − 1 を f f f の a a a における留数 と呼び、
Res ( f , a ) : = c − 1 \operatorname{Res}(f, a) := c_{-1} Res ( f , a ) := c − 1 と書く。ローラン係数の積分表示(§2)で n = − 1 n = -1 n = − 1 とすれば、0 < ρ < R 0 < \rho < R 0 < ρ < R に対して
Res ( f , a ) = 1 2 π i ∮ ∣ z − a ∣ = ρ f ( z ) d z \operatorname{Res}(f,a) = \frac{1}{2\pi i}\oint_{|z-a| = \rho} f(z)\,dz Res ( f , a ) = 2 π i 1 ∮ ∣ z − a ∣ = ρ f ( z ) d z でもある。
定義そのものは「ローラン展開せよ」と言っていますが、極の場合には展開せずに留数を取り出す公式があります。実用上はこの公式だけで九割方の計算が片付きます。
Proof(Proposition 4.2) (1) :仮定より D ′ ( a , R ) D'(a,R) D ′ ( a , R ) で f ( z ) = ∑ n ≥ − m c n ( z − a ) n f(z) = \sum_{n \ge -m}c_n(z-a)^n f ( z ) = ∑ n ≥ − m c n ( z − a ) n なので、
φ ( z ) = ( z − a ) m f ( z ) = ∑ k ≥ 0 c k − m ( z − a ) k \varphi(z) = (z-a)^m f(z) = \sum_{k \ge 0}c_{k-m}(z-a)^k φ ( z ) = ( z − a ) m f ( z ) = k ≥ 0 ∑ c k − m ( z − a ) k であり、右辺は D ′ ( a , R ) D'(a,R) D ′ ( a , R ) で収束するべき級数、したがって収束半径は R R R 以上で D ( a , R ) D(a,R) D ( a , R ) で正則関数を定めます(Proposition 3.4 の証明と同じ議論です)。この級数の ( z − a ) m − 1 (z-a)^{m-1} ( z − a ) m − 1 の係数は k = m − 1 k = m-1 k = m − 1 の項、すなわち c − 1 c_{-1} c − 1 です。一方、正則関数のテイラー係数の公式(Theorem 5.2[正則関数の強力な性質] )から ( z − a ) m − 1 (z-a)^{m-1} ( z − a ) m − 1 の係数は φ ( m − 1 ) ( a ) / ( m − 1 ) ! \varphi^{(m-1)}(a)/(m-1)! φ ( m − 1 ) ( a ) / ( m − 1 )! です。両者を等しいとおけば結論を得ます。φ \varphi φ は a a a で正則なので φ ( m − 1 ) \varphi^{(m-1)} φ ( m − 1 ) は a a a で連続で、φ ( m − 1 ) ( a ) = lim z → a φ ( m − 1 ) ( z ) \varphi^{(m-1)}(a) = \lim_{z\to a}\varphi^{(m-1)}(z) φ ( m − 1 ) ( a ) = lim z → a φ ( m − 1 ) ( z ) と極限で書けます。
(2) :(1) で m = 1 m = 1 m = 1 とすると 0 ! = 1 0! = 1 0 ! = 1 、0 0 0 階微分は恒等写像なので、Res ( f , a ) = lim z → a ( z − a ) f ( z ) \operatorname{Res}(f,a) = \lim_{z\to a}(z-a)f(z) Res ( f , a ) = lim z → a ( z − a ) f ( z ) です。
(3) :h ( a ) = 0 h(a) = 0 h ( a ) = 0 かつ h ′ ( a ) ≠ 0 h'(a) \ne 0 h ′ ( a ) = 0 より a a a は h h h の 1 1 1 位の零点で、a a a の近傍で h ( z ) = ( z − a ) h 1 ( z ) h(z) = (z-a)h_1(z) h ( z ) = ( z − a ) h 1 ( z ) 、h 1 h_1 h 1 は正則で h 1 ( a ) = h ′ ( a ) ≠ 0 h_1(a) = h'(a) \ne 0 h 1 ( a ) = h ′ ( a ) = 0 と書けます。よって f ( z ) = g ( z ) / h 1 ( z ) z − a f(z) = \dfrac{g(z)/h_1(z)}{z-a} f ( z ) = z − a g ( z ) / h 1 ( z ) であり、g / h 1 g/h_1 g / h 1 は a a a で正則、a a a での値は g ( a ) / h ′ ( a ) ≠ 0 g(a)/h'(a) \ne 0 g ( a ) / h ′ ( a ) = 0 ですから、Proposition 3.4 (1)(iii) より a a a は単純極です。(2) を使って
Res ( f , a ) = lim z → a ( z − a ) f ( z ) = lim z → a g ( z ) h 1 ( z ) = g ( a ) h ′ ( a ) \operatorname{Res}(f,a) = \lim_{z\to a}(z-a)f(z) = \lim_{z\to a}\frac{g(z)}{h_1(z)} = \frac{g(a)}{h'(a)} Res ( f , a ) = z → a lim ( z − a ) f ( z ) = z → a lim h 1 ( z ) g ( z ) = h ′ ( a ) g ( a ) を得ます。
∎ Example 4.3 (留数の計算例 )
(a) f ( z ) = z ( z − 1 ) ( z + 2 ) 2 f(z) = \dfrac{z}{(z-1)(z+2)^2} f ( z ) = ( z − 1 ) ( z + 2 ) 2 z 。z = 1 z=1 z = 1 は単純極なので Proposition 4.2 (2) より
Res ( f , 1 ) = lim z → 1 z ( z + 2 ) 2 = 1 9 . \operatorname{Res}(f,1) = \lim_{z\to1}\frac{z}{(z+2)^2} = \frac{1}{9}. Res ( f , 1 ) = z → 1 lim ( z + 2 ) 2 z = 9 1 . z = − 2 z=-2 z = − 2 は 2 2 2 位の極なので (1) で m = 2 m=2 m = 2 とし、
Res ( f , − 2 ) = lim z → − 2 d d z [ z z − 1 ] = lim z → − 2 ( z − 1 ) − z ( z − 1 ) 2 = − 1 9 . \operatorname{Res}(f,-2) = \lim_{z\to-2}\frac{d}{dz}\left[\frac{z}{z-1}\right]
= \lim_{z\to-2}\frac{(z-1) - z}{(z-1)^2} = \frac{-1}{9}. Res ( f , − 2 ) = z → − 2 lim d z d [ z − 1 z ] = z → − 2 lim ( z − 1 ) 2 ( z − 1 ) − z = 9 − 1 . 二つの留数の和が 0 0 0 になっているのは偶然ではありません。分母の次数が分子より 2 2 2 以上大きい有理関数では、後の Lemma 6.1 (i) と Theorem 5.1 を大きな円周に適用すると留数の総和が 0 0 0 と分かります。検算に使えます。
(b) f ( z ) = e z z 3 f(z) = \dfrac{e^z}{z^3} f ( z ) = z 3 e z 。m = 3 m=3 m = 3 として
Res ( f , 0 ) = 1 2 ! lim z → 0 d 2 d z 2 [ e z ] = 1 2 . \operatorname{Res}(f,0) = \frac{1}{2!}\lim_{z\to0}\frac{d^2}{dz^2}\bigl[e^z\bigr] = \frac{1}{2}. Res ( f , 0 ) = 2 ! 1 z → 0 lim d z 2 d 2 [ e z ] = 2 1 . ローラン展開 e z / z 3 = z − 3 + z − 2 + 1 2 z − 1 + 1 6 + ⋯ e^z/z^3 = z^{-3} + z^{-2} + \frac{1}{2}z^{-1} + \frac{1}{6} + \cdots e z / z 3 = z − 3 + z − 2 + 2 1 z − 1 + 6 1 + ⋯ の z − 1 z^{-1} z − 1 の係数と一致します。
(c) f ( z ) = cot z = cos z sin z f(z) = \cot z = \dfrac{\cos z}{\sin z} f ( z ) = cot z = sin z cos z の z = k π z = k\pi z = k π 。g = cos g = \cos g = cos 、h = sin h = \sin h = sin とすると g ( k π ) = ( − 1 ) k ≠ 0 g(k\pi) = (-1)^k \ne 0 g ( k π ) = ( − 1 ) k = 0 、h ( k π ) = 0 h(k\pi) = 0 h ( k π ) = 0 、h ′ ( k π ) = cos ( k π ) = ( − 1 ) k ≠ 0 h'(k\pi) = \cos(k\pi) = (-1)^k \ne 0 h ′ ( k π ) = cos ( k π ) = ( − 1 ) k = 0 なので Proposition 4.2 (3) より
Res ( cot , k π ) = cos ( k π ) cos ( k π ) = 1. \operatorname{Res}(\cot, k\pi) = \frac{\cos(k\pi)}{\cos(k\pi)} = 1. Res ( cot , k π ) = cos ( k π ) cos ( k π ) = 1. すべての極で留数が 1 1 1 です。
(d) f ( z ) = z 2 e 1 / z f(z) = z^2 e^{1/z} f ( z ) = z 2 e 1/ z の z = 0 z = 0 z = 0 。真性特異点なので (1)〜(3) は使えません。定義に戻ってローラン展開します。
z 2 e 1 / z = z 2 ∑ n ≥ 0 1 n ! z − n = ∑ n ≥ 0 1 n ! z 2 − n z^2 e^{1/z} = z^2\sum_{n\ge0}\frac{1}{n!}z^{-n} = \sum_{n \ge 0}\frac{1}{n!}z^{2-n} z 2 e 1/ z = z 2 n ≥ 0 ∑ n ! 1 z − n = n ≥ 0 ∑ n ! 1 z 2 − n で、z − 1 z^{-1} z − 1 が現れるのは 2 − n = − 1 2-n = -1 2 − n = − 1 、すなわち n = 3 n = 3 n = 3 のときですから Res ( f , 0 ) = 1 3 ! = 1 6 \operatorname{Res}(f,0) = \frac{1}{3!} = \frac{1}{6} Res ( f , 0 ) = 3 ! 1 = 6 1 です。留数は極の専売特許ではありません。
Theorem 5.1 (留数定理 )
D ⊂ C D \subset \mathbb{C} D ⊂ C を単連結領域とする。C C C を D D D に含まれる区分的に C 1 C^1 C 1 級の単純閉曲線とし、正の向き(内部を左に見る向き)を与える。C C C の内部(ジョルダンの曲線定理でいう有界成分)を int C \operatorname{int}C int C と書く。相異なる有限個の点 a 1 , … , a n ∈ int C a_1,\dots,a_n \in \operatorname{int}C a 1 , … , a n ∈ int C を取り、f f f は D ∖ { a 1 , … , a n } D \setminus \{a_1,\dots,a_n\} D ∖ { a 1 , … , a n } で正則とする。このとき
∮ C f ( z ) d z = 2 π i ∑ j = 1 n Res ( f , a j ) . \oint_C f(z)\,dz = 2\pi i \sum_{j=1}^{n}\operatorname{Res}(f,a_j). ∮ C f ( z ) d z = 2 π i j = 1 ∑ n Res ( f , a j ) . Proof(Theorem 5.1) D D D は単連結なので int C ⊂ D \operatorname{int}C \subset D int C ⊂ D です(単連結領域内の閉曲線の内部は領域に含まれます)。したがって各 a j a_j a j は f f f の孤立特異点で、留数が定義されています。
小円を取る 。ε > 0 \varepsilon > 0 ε > 0 を
ε < 1 2 min j ≠ k ∣ a j − a k ∣ , ε < min j dist ( a j , ∂ ( int C ) ) \varepsilon < \tfrac{1}{2}\min_{j \ne k}|a_j - a_k|,
\qquad
\varepsilon < \min_{j}\operatorname{dist}(a_j, \partial(\operatorname{int}C)) ε < 2 1 j = k min ∣ a j − a k ∣ , ε < j min dist ( a j , ∂ ( int C )) を満たすように取ります(a j a_j a j は有限個で、int C \operatorname{int}C int C は開集合なので、このような ε \varepsilon ε が存在します)。γ j \gamma_j γ j を正の向きの円周 ∣ z − a j ∣ = ε |z - a_j| = \varepsilon ∣ z − a j ∣ = ε とすると、閉円板 D ( a j , ε ) ‾ \overline{D(a_j,\varepsilon)} D ( a j , ε ) は互いに交わらず、すべて int C \operatorname{int}C int C に含まれます。
積分路を分解する 。f f f は「C C C の内側かつすべての γ j \gamma_j γ j の外側」の閉領域を含む開集合 D ∖ { a 1 , … , a n } D \setminus \{a_1,\dots,a_n\} D ∖ { a 1 , … , a n } で正則です。§2 の多重連結領域に対するコーシーの積分定理より
∮ C f ( z ) d z = ∑ j = 1 n ∮ γ j f ( z ) d z . \oint_C f(z)\,dz = \sum_{j=1}^{n}\oint_{\gamma_j}f(z)\,dz . ∮ C f ( z ) d z = j = 1 ∑ n ∮ γ j f ( z ) d z . 小円上で項別積分する 。j j j を固定します。f f f は D ′ ( a j , R j ) D'(a_j, R_j) D ′ ( a j , R j ) でローラン展開 f ( z ) = ∑ k = − ∞ ∞ c k ( j ) ( z − a j ) k f(z) = \sum_{k=-\infty}^{\infty}c^{(j)}_k(z-a_j)^k f ( z ) = ∑ k = − ∞ ∞ c k ( j ) ( z − a j ) k を持ちます。必要なら ε \varepsilon ε をさらに小さくして ε < R j \varepsilon < R_j ε < R j としておきます。円周 γ j \gamma_j γ j は円環 D ′ ( a j , R j ) D'(a_j,R_j) D ′ ( a j , R j ) のコンパクト部分集合なので、この級数は γ j \gamma_j γ j 上で一様収束します。一様収束する連続関数列は長さ有限の曲線上で項別積分できる(Lemma 2.2[正則関数の強力な性質] )ので
∮ γ j f ( z ) d z = ∑ k = − ∞ ∞ c k ( j ) ∮ γ j ( z − a j ) k d z . \oint_{\gamma_j} f(z)\,dz = \sum_{k=-\infty}^{\infty}c^{(j)}_k \oint_{\gamma_j}(z-a_j)^k\,dz . ∮ γ j f ( z ) d z = k = − ∞ ∑ ∞ c k ( j ) ∮ γ j ( z − a j ) k d z . 右辺の積分を計算します。z = a j + ε e i θ z = a_j + \varepsilon e^{i\theta} z = a j + ε e i θ (0 ≤ θ ≤ 2 π 0 \le \theta \le 2\pi 0 ≤ θ ≤ 2 π )と置くと d z = i ε e i θ d θ dz = i\varepsilon e^{i\theta}d\theta d z = i ε e i θ d θ なので
∮ γ j ( z − a j ) k d z = ∫ 0 2 π ε k e i k θ ⋅ i ε e i θ d θ = i ε k + 1 ∫ 0 2 π e i ( k + 1 ) θ d θ . \oint_{\gamma_j}(z-a_j)^k dz
= \int_0^{2\pi}\varepsilon^k e^{ik\theta}\cdot i\varepsilon e^{i\theta}\,d\theta
= i\,\varepsilon^{\,k+1}\int_0^{2\pi}e^{i(k+1)\theta}\,d\theta . ∮ γ j ( z − a j ) k d z = ∫ 0 2 π ε k e ik θ ⋅ i ε e i θ d θ = i ε k + 1 ∫ 0 2 π e i ( k + 1 ) θ d θ . k ≠ − 1 k \ne -1 k = − 1 のとき ∫ 0 2 π e i ( k + 1 ) θ d θ = [ e i ( k + 1 ) θ i ( k + 1 ) ] 0 2 π = 0 \int_0^{2\pi}e^{i(k+1)\theta}d\theta = \left[\dfrac{e^{i(k+1)\theta}}{i(k+1)}\right]_0^{2\pi} = 0 ∫ 0 2 π e i ( k + 1 ) θ d θ = [ i ( k + 1 ) e i ( k + 1 ) θ ] 0 2 π = 0 (e 2 π i ( k + 1 ) = 1 e^{2\pi i(k+1)} = 1 e 2 π i ( k + 1 ) = 1 )です。k = − 1 k = -1 k = − 1 のときは被積分関数が 1 1 1 なので 2 π 2\pi 2 π となり、積分は i ⋅ ε 0 ⋅ 2 π = 2 π i i \cdot \varepsilon^0 \cdot 2\pi = 2\pi i i ⋅ ε 0 ⋅ 2 π = 2 π i です。よって
∮ γ j f ( z ) d z = 2 π i c − 1 ( j ) = 2 π i Res ( f , a j ) \oint_{\gamma_j}f(z)\,dz = 2\pi i\, c^{(j)}_{-1} = 2\pi i \operatorname{Res}(f,a_j) ∮ γ j f ( z ) d z = 2 π i c − 1 ( j ) = 2 π i Res ( f , a j ) となり(最後の等号は Definition 4.1 )、これを j j j について足せば結論を得ます。
∎ ここから先は応用です。実積分を計算するときの手順はいつも同じで、次の三つに分かれます。
実積分の型 使う輪郭 余分な部分を消す道具 ∫ − ∞ ∞ P ( x ) / Q ( x ) d x \int_{-\infty}^{\infty}P(x)/Q(x)\,dx ∫ − ∞ ∞ P ( x ) / Q ( x ) d x (deg Q ≥ deg P + 2 \deg Q \ge \deg P + 2 deg Q ≥ deg P + 2 )上半平面の半円 Lemma 6.1 (i)∫ 0 2 π R ( cos θ , sin θ ) d θ \int_0^{2\pi}R(\cos\theta,\sin\theta)\,d\theta ∫ 0 2 π R ( cos θ , sin θ ) d θ 単位円 ∣ z ∣ = 1 \lvert z\rvert = 1 ∣ z ∣ = 1 不要(輪郭が最初から閉じている) ∫ − ∞ ∞ f ( x ) e i λ x d x \int_{-\infty}^{\infty}f(x)e^{i\lambda x}dx ∫ − ∞ ∞ f ( x ) e iλ x d x (λ > 0 \lambda > 0 λ > 0 )上半平面の半円 Lemma 6.1 (ii)(ジョルダンの補題)実軸上に単純極があるとき 極を小円でくぼませた輪郭 Lemma 8.3 と主値∫ 0 ∞ x a − 1 g ( x ) d x \int_0^{\infty}x^{a-1}g(x)\,dx ∫ 0 ∞ x a − 1 g ( x ) d x (a a a は非整数)鍵穴(Appendix) 分枝の指定と大小円の評価
有理関数の広義積分に使う輪郭。実軸上の線分と上半平面の半円弧をつなぎ、上半平面にある極(×)だけを内部に取り込みます。半径 R を大きくすると弧の寄与が消えます。 輪郭を閉じるために付け足した半円弧の寄与は、半径を大きくすると消えてほしいわけです。それを保証するのが次の補題です。(i) と (ii) の仮定の違いに注意してください。(ii) のほうが f f f に対する要求がゆるく、そのぶん e i λ z e^{i\lambda z} e iλ z という因子が必要です。
Lemma 6.1 (大円弧の寄与が消える条件 )
R 0 > 0 R_0 > 0 R 0 > 0 とし、Γ R \Gamma_R Γ R を上半平面の半円弧 z = R e i θ z = Re^{i\theta} z = R e i θ (0 ≤ θ ≤ π 0 \le \theta \le \pi 0 ≤ θ ≤ π )とする。f f f は { z : ∣ z ∣ ≥ R 0 , Im z ≥ 0 } \{z : |z| \ge R_0,\ \operatorname{Im}z \ge 0\} { z : ∣ z ∣ ≥ R 0 , Im z ≥ 0 } で連続とし、M ( R ) : = max 0 ≤ θ ≤ π ∣ f ( R e i θ ) ∣ M(R) := \max_{0 \le \theta \le \pi}|f(Re^{i\theta})| M ( R ) := max 0 ≤ θ ≤ π ∣ f ( R e i θ ) ∣ とおく。
(i) lim R → ∞ R M ( R ) = 0 \lim_{R\to\infty} R\,M(R) = 0 lim R → ∞ R M ( R ) = 0 ならば lim R → ∞ ∫ Γ R f ( z ) d z = 0 \displaystyle\lim_{R\to\infty}\int_{\Gamma_R}f(z)\,dz = 0 R → ∞ lim ∫ Γ R f ( z ) d z = 0 。
(ii)(ジョルダンの補題) λ > 0 \lambda > 0 λ > 0 とする。R > R 0 R > R_0 R > R 0 のとき
∣ ∫ Γ R f ( z ) e i λ z d z ∣ ≤ π λ M ( R ) \left|\int_{\Gamma_R}f(z)e^{i\lambda z}\,dz\right| \le \frac{\pi}{\lambda}M(R) ∫ Γ R f ( z ) e iλ z d z ≤ λ π M ( R ) が成り立つ。特に lim R → ∞ M ( R ) = 0 \lim_{R\to\infty}M(R) = 0 lim R → ∞ M ( R ) = 0 ならば lim R → ∞ ∫ Γ R f ( z ) e i λ z d z = 0 \displaystyle\lim_{R\to\infty}\int_{\Gamma_R}f(z)e^{i\lambda z}dz = 0 R → ∞ lim ∫ Γ R f ( z ) e iλ z d z = 0 。
Proof(Lemma 6.1) (i) :Γ R \Gamma_R Γ R の長さは π R \pi R π R なので、§2 の ML 不等式より
∣ ∫ Γ R f ( z ) d z ∣ ≤ π R ⋅ M ( R ) ⟶ 0. \left|\int_{\Gamma_R}f(z)\,dz\right| \le \pi R \cdot M(R) \longrightarrow 0 . ∫ Γ R f ( z ) d z ≤ π R ⋅ M ( R ) ⟶ 0. (ii) :z = R e i θ = R cos θ + i R sin θ z = Re^{i\theta} = R\cos\theta + iR\sin\theta z = R e i θ = R cos θ + i R sin θ に対し
∣ e i λ z ∣ = ∣ e i λ R cos θ ∣ ⋅ ∣ e − λ R sin θ ∣ = e − λ R sin θ |e^{i\lambda z}| = |e^{i\lambda R\cos\theta}|\cdot|e^{-\lambda R\sin\theta}| = e^{-\lambda R\sin\theta} ∣ e iλ z ∣ = ∣ e iλ R c o s θ ∣ ⋅ ∣ e − λ R s i n θ ∣ = e − λ R s i n θ です(第一の因子は絶対値 1 1 1 )。∣ d z ∣ = R d θ |dz| = R\,d\theta ∣ d z ∣ = R d θ に注意して
∣ ∫ Γ R f ( z ) e i λ z d z ∣ ≤ ∫ 0 π M ( R ) e − λ R sin θ R d θ = 2 R M ( R ) ∫ 0 π / 2 e − λ R sin θ d θ \left|\int_{\Gamma_R}f(z)e^{i\lambda z}dz\right|
\le \int_0^{\pi}M(R)\,e^{-\lambda R \sin\theta}\,R\,d\theta
= 2RM(R)\int_0^{\pi/2}e^{-\lambda R\sin\theta}\,d\theta ∫ Γ R f ( z ) e iλ z d z ≤ ∫ 0 π M ( R ) e − λ R s i n θ R d θ = 2 R M ( R ) ∫ 0 π /2 e − λ R s i n θ d θ となります(sin ( π − θ ) = sin θ \sin(\pi-\theta) = \sin\theta sin ( π − θ ) = sin θ による対称性で [ 0 , π / 2 ] [0,\pi/2] [ 0 , π /2 ] の積分の 2 倍にしました)。
ここでジョルダンの不等式 sin θ ≥ 2 θ π \sin\theta \ge \dfrac{2\theta}{\pi} sin θ ≥ π 2 θ (0 ≤ θ ≤ π / 2 0 \le \theta \le \pi/2 0 ≤ θ ≤ π /2 )を使います。証明しておきます。ψ ( θ ) = sin θ / θ \psi(\theta) = \sin\theta/\theta ψ ( θ ) = sin θ / θ (θ > 0 \theta > 0 θ > 0 )とおくと ψ ′ ( θ ) = θ cos θ − sin θ θ 2 \psi'(\theta) = \dfrac{\theta\cos\theta - \sin\theta}{\theta^2} ψ ′ ( θ ) = θ 2 θ cos θ − sin θ です。分子を u ( θ ) = θ cos θ − sin θ u(\theta) = \theta\cos\theta - \sin\theta u ( θ ) = θ cos θ − sin θ とすると u ( 0 ) = 0 u(0) = 0 u ( 0 ) = 0 、u ′ ( θ ) = cos θ − θ sin θ − cos θ = − θ sin θ < 0 u'(\theta) = \cos\theta - \theta\sin\theta - \cos\theta = -\theta\sin\theta < 0 u ′ ( θ ) = cos θ − θ sin θ − cos θ = − θ sin θ < 0 (0 < θ ≤ π / 2 0 < \theta \le \pi/2 0 < θ ≤ π /2 )なので u < 0 u < 0 u < 0 、したがって ψ \psi ψ は ( 0 , π / 2 ] (0,\pi/2] ( 0 , π /2 ] で狭義単調減少です。よって ψ ( θ ) ≥ ψ ( π / 2 ) = 1 π / 2 = 2 π \psi(\theta) \ge \psi(\pi/2) = \dfrac{1}{\pi/2} = \dfrac{2}{\pi} ψ ( θ ) ≥ ψ ( π /2 ) = π /2 1 = π 2 、すなわち sin θ ≥ 2 θ / π \sin\theta \ge 2\theta/\pi sin θ ≥ 2 θ / π です。
これを使うと e − λ R sin θ ≤ e − 2 λ R θ / π e^{-\lambda R\sin\theta} \le e^{-2\lambda R\theta/\pi} e − λ R s i n θ ≤ e − 2 λ R θ / π なので
∫ 0 π / 2 e − λ R sin θ d θ ≤ ∫ 0 π / 2 e − 2 λ R θ / π d θ = π 2 λ R ( 1 − e − λ R ) ≤ π 2 λ R \int_0^{\pi/2}e^{-\lambda R\sin\theta}d\theta
\le \int_0^{\pi/2}e^{-2\lambda R\theta/\pi}d\theta
= \frac{\pi}{2\lambda R}\left(1 - e^{-\lambda R}\right)
\le \frac{\pi}{2\lambda R} ∫ 0 π /2 e − λ R s i n θ d θ ≤ ∫ 0 π /2 e − 2 λ R θ / π d θ = 2 λ R π ( 1 − e − λ R ) ≤ 2 λ R π となり、これを代入して 2 R M ( R ) ⋅ π 2 λ R = π λ M ( R ) 2RM(R)\cdot\dfrac{\pi}{2\lambda R} = \dfrac{\pi}{\lambda}M(R) 2 R M ( R ) ⋅ 2 λ R π = λ π M ( R ) を得ます。
∎ Proposition 6.2 (有理関数の広義積分 )
P , Q P, Q P , Q を複素係数の多項式とし、次を仮定する。
deg Q ≥ deg P + 2 \deg Q \ge \deg P + 2 deg Q ≥ deg P + 2 。
Q Q Q は実軸上に零点を持たない。
このとき f = P / Q f = P/Q f = P / Q に対し ∫ − ∞ ∞ f ( x ) d x \int_{-\infty}^{\infty}f(x)\,dx ∫ − ∞ ∞ f ( x ) d x は絶対収束し、
∫ − ∞ ∞ P ( x ) Q ( x ) d x = 2 π i ∑ Im a > 0 Res ( f , a ) \int_{-\infty}^{\infty}\frac{P(x)}{Q(x)}\,dx = 2\pi i \sum_{\operatorname{Im}a > 0}\operatorname{Res}(f, a) ∫ − ∞ ∞ Q ( x ) P ( x ) d x = 2 π i Im a > 0 ∑ Res ( f , a ) が成り立つ。和は上半平面にある Q Q Q のすべての零点 a a a にわたる。
Proof(Proposition 6.2) p = deg P p = \deg P p = deg P 、q = deg Q q = \deg Q q = deg Q 、最高次係数を α , β \alpha, \beta α , β とします。∣ z ∣ → ∞ |z|\to\infty ∣ z ∣ → ∞ のとき z q − p P ( z ) / Q ( z ) → α / β z^{\,q-p}P(z)/Q(z) \to \alpha/\beta z q − p P ( z ) / Q ( z ) → α / β なので、ある R 1 ≥ 1 R_1 \ge 1 R 1 ≥ 1 と K > 0 K > 0 K > 0 があって
∣ z ∣ ≥ R 1 ⟹ ∣ f ( z ) ∣ ≤ K ∣ z ∣ p − q ≤ K ∣ z ∣ − 2 |z| \ge R_1 \implies |f(z)| \le K|z|^{\,p-q} \le K|z|^{-2} ∣ z ∣ ≥ R 1 ⟹ ∣ f ( z ) ∣ ≤ K ∣ z ∣ p − q ≤ K ∣ z ∣ − 2 です(最後の不等号に q − p ≥ 2 q - p \ge 2 q − p ≥ 2 と ∣ z ∣ ≥ R 1 ≥ 1 |z| \ge R_1 \ge 1 ∣ z ∣ ≥ R 1 ≥ 1 を使いました)。したがって ∫ ∣ x ∣ ≥ R 1 ∣ f ( x ) ∣ d x ≤ 2 ∫ R 1 ∞ K x − 2 d x < ∞ \int_{|x| \ge R_1}|f(x)|dx \le 2\int_{R_1}^{\infty}Kx^{-2}dx < \infty ∫ ∣ x ∣ ≥ R 1 ∣ f ( x ) ∣ d x ≤ 2 ∫ R 1 ∞ K x − 2 d x < ∞ であり、[ − R 1 , R 1 ] [-R_1,R_1] [ − R 1 , R 1 ] 上では f f f が連続(Q Q Q は実軸上で 0 0 0 にならない)なので、広義積分は絶対収束します。また R ≥ R 1 R \ge R_1 R ≥ R 1 に対し R M ( R ) ≤ K / R → 0 R\,M(R) \le K/R \to 0 R M ( R ) ≤ K / R → 0 なので Lemma 6.1 (i) の仮定が満たされます。
Q Q Q の零点は有限個です。R > max ( R 1 , max a ∣ a ∣ ) R > \max(R_1, \max_a|a|) R > max ( R 1 , max a ∣ a ∣ ) を取り、C R C_R C R を線分 [ − R , R ] [-R,R] [ − R , R ] と半円弧 Γ R \Gamma_R Γ R をつないだ正の向きの単純閉曲線とします。C R C_R C R とその内部は単連結領域 C \mathbb{C} C に含まれ、f f f は Q Q Q の零点を除いて正則で、C R C_R C R の内部にある零点はちょうど上半平面の零点全部です(実軸上には零点がなく、R R R の取り方から下半平面の零点は内部に入りません)。Theorem 5.1 より
∫ − R R f ( x ) d x + ∫ Γ R f ( z ) d z = 2 π i ∑ Im a > 0 Res ( f , a ) \int_{-R}^{R}f(x)\,dx + \int_{\Gamma_R}f(z)\,dz = 2\pi i \sum_{\operatorname{Im}a > 0}\operatorname{Res}(f,a) ∫ − R R f ( x ) d x + ∫ Γ R f ( z ) d z = 2 π i Im a > 0 ∑ Res ( f , a ) です。右辺は R R R によりません。R → ∞ R \to \infty R → ∞ とすると、左辺第 1 項は絶対収束性より ∫ − ∞ ∞ f ( x ) d x \int_{-\infty}^{\infty}f(x)dx ∫ − ∞ ∞ f ( x ) d x に収束し、第 2 項は Lemma 6.1 (i) より 0 0 0 に収束します。
∎ Example 6.3 (分母が 4 次の場合 )
∫ − ∞ ∞ d x x 4 + 1 \displaystyle\int_{-\infty}^{\infty}\frac{dx}{x^4+1} ∫ − ∞ ∞ x 4 + 1 d x を計算します。P = 1 P = 1 P = 1 、Q = z 4 + 1 Q = z^4+1 Q = z 4 + 1 で deg Q − deg P = 4 ≥ 2 \deg Q - \deg P = 4 \ge 2 deg Q − deg P = 4 ≥ 2 、Q ( x ) = x 4 + 1 ≥ 1 > 0 Q(x) = x^4 + 1 \ge 1 > 0 Q ( x ) = x 4 + 1 ≥ 1 > 0 より実軸上に零点はないので Proposition 6.2 が使えます。
z 4 = − 1 = e i π z^4 = -1 = e^{i\pi} z 4 = − 1 = e iπ の解は z k = e i π ( 2 k + 1 ) / 4 z_k = e^{i\pi(2k+1)/4} z k = e iπ ( 2 k + 1 ) /4 (k = 0 , 1 , 2 , 3 k = 0,1,2,3 k = 0 , 1 , 2 , 3 )、すなわち e i π / 4 , e 3 i π / 4 , e 5 i π / 4 , e 7 i π / 4 e^{i\pi/4},\ e^{3i\pi/4},\ e^{5i\pi/4},\ e^{7i\pi/4} e iπ /4 , e 3 iπ /4 , e 5 iπ /4 , e 7 iπ /4 です。上半平面にあるのは偏角が π / 4 \pi/4 π /4 と 3 π / 4 3\pi/4 3 π /4 の二つです。いずれも Q Q Q の単純零点(Q ′ = 4 z 3 Q' = 4z^3 Q ′ = 4 z 3 は 0 0 0 にならない)なので、Proposition 4.2 (3) より、a a a をそのどちらかとして
Res ( f , a ) = 1 4 a 3 = a 4 a 4 = a 4 ⋅ ( − 1 ) = − a 4 \operatorname{Res}(f,a) = \frac{1}{4a^3} = \frac{a}{4a^4} = \frac{a}{4\cdot(-1)} = -\frac{a}{4} Res ( f , a ) = 4 a 3 1 = 4 a 4 a = 4 ⋅ ( − 1 ) a = − 4 a です(a 4 = − 1 a^4 = -1 a 4 = − 1 を使って a 3 a^3 a 3 を a a a に直しました。この書き換えで計算がぐっと楽になります)。したがって
∑ Res = − 1 4 ( e i π / 4 + e 3 i π / 4 ) = − 1 4 ( 2 2 ( 1 + i ) + 2 2 ( − 1 + i ) ) = − 1 4 ⋅ 2 i \sum \operatorname{Res} = -\frac{1}{4}\left(e^{i\pi/4} + e^{3i\pi/4}\right)
= -\frac{1}{4}\left(\frac{\sqrt2}{2}(1+i) + \frac{\sqrt2}{2}(-1+i)\right)
= -\frac{1}{4}\cdot\sqrt{2}\,i ∑ Res = − 4 1 ( e iπ /4 + e 3 iπ /4 ) = − 4 1 ( 2 2 ( 1 + i ) + 2 2 ( − 1 + i ) ) = − 4 1 ⋅ 2 i となり、
∫ − ∞ ∞ d x x 4 + 1 = 2 π i ⋅ ( − 2 i 4 ) = 2 π 2 4 = π 2 ≈ 2.2214. \int_{-\infty}^{\infty}\frac{dx}{x^4+1} = 2\pi i \cdot\left(-\frac{\sqrt2 i}{4}\right) = \frac{2\pi\sqrt2}{4} = \frac{\pi}{\sqrt2} \approx 2.2214 . ∫ − ∞ ∞ x 4 + 1 d x = 2 π i ⋅ ( − 4 2 i ) = 4 2 π 2 = 2 π ≈ 2.2214. 被積分関数は正なので答えが正の実数になったのは整合的です。また ∣ x ∣ ≥ 1 |x| \ge 1 ∣ x ∣ ≥ 1 では 1 x 4 + 1 ≤ 1 x 2 + 1 \frac{1}{x^4+1} \le \frac{1}{x^2+1} x 4 + 1 1 ≤ x 2 + 1 1 、∣ x ∣ ≤ 1 |x| \le 1 ∣ x ∣ ≤ 1 では 1 x 4 + 1 ≤ 1 \frac{1}{x^4+1} \le 1 x 4 + 1 1 ≤ 1 なので、積分値は 2 + ∫ ∣ x ∣ ≥ 1 d x x 2 + 1 = 2 + π 2 ≈ 3.57 2 + \int_{|x|\ge1}\frac{dx}{x^2+1} = 2 + \frac{\pi}{2} \approx 3.57 2 + ∫ ∣ x ∣ ≥ 1 x 2 + 1 d x = 2 + 2 π ≈ 3.57 以下と分かり、2.22 2.22 2.22 という値と矛盾しません。
θ \theta θ について周期 2 π 2\pi 2 π の三角有理式の積分 ∫ 0 2 π R ( cos θ , sin θ ) d θ \int_0^{2\pi}R(\cos\theta,\sin\theta)\,d\theta ∫ 0 2 π R ( cos θ , sin θ ) d θ は、実積分が最初から閉じた区間を一周しているので、輪郭を付け足す必要すらありません。z = e i θ z = e^{i\theta} z = e i θ と置けば、これはそのまま単位円周に沿う複素積分になります。オイラーの公式から
cos θ = z + z − 1 2 , sin θ = z − z − 1 2 i , d θ = d z i z \cos\theta = \frac{z + z^{-1}}{2},\qquad
\sin\theta = \frac{z - z^{-1}}{2i},\qquad
d\theta = \frac{dz}{iz} cos θ = 2 z + z − 1 , sin θ = 2 i z − z − 1 , d θ = i z d z です。最後の関係は、単位円周を γ ( θ ) = e i θ \gamma(\theta) = e^{i\theta} γ ( θ ) = e i θ (0 ≤ θ ≤ 2 π 0 \le \theta \le 2\pi 0 ≤ θ ≤ 2 π )で径数付けると γ ′ ( θ ) = i e i θ = i z \gamma'(\theta) = ie^{i\theta} = iz γ ′ ( θ ) = i e i θ = i z なので、任意の連続関数 F F F に対し
∮ ∣ z ∣ = 1 F ( z ) i z d z = ∫ 0 2 π F ( e i θ ) i e i θ i e i θ d θ = ∫ 0 2 π F ( e i θ ) d θ \oint_{|z|=1}\frac{F(z)}{iz}\,dz = \int_0^{2\pi}\frac{F(e^{i\theta})}{ie^{i\theta}}\,ie^{i\theta}\,d\theta = \int_0^{2\pi}F(e^{i\theta})\,d\theta ∮ ∣ z ∣ = 1 i z F ( z ) d z = ∫ 0 2 π i e i θ F ( e i θ ) i e i θ d θ = ∫ 0 2 π F ( e i θ ) d θ が成り立つ、という意味です。
Proposition 7.1 (a + b cosθ 型の積分 )
a , b a, b a , b を実数とし a > ∣ b ∣ ≥ 0 a > |b| \ge 0 a > ∣ b ∣ ≥ 0 とする。このとき a + b cos θ ≥ a − ∣ b ∣ > 0 a + b\cos\theta \ge a - |b| > 0 a + b cos θ ≥ a − ∣ b ∣ > 0 で被積分関数は連続であり、
∫ 0 2 π d θ a + b cos θ = 2 π a 2 − b 2 . \int_0^{2\pi}\frac{d\theta}{a + b\cos\theta} = \frac{2\pi}{\sqrt{a^2-b^2}} . ∫ 0 2 π a + b cos θ d θ = a 2 − b 2 2 π . Proof(Proposition 7.1) b = 0 b = 0 b = 0 のときは左辺 = 2 π / a = 2\pi/a = 2 π / a 、右辺 = 2 π / a 2 = 2 π / a = 2\pi/\sqrt{a^2} = 2\pi/a = 2 π / a 2 = 2 π / a (a > 0 a > 0 a > 0 )で成立します。b < 0 b < 0 b < 0 のときは θ ↦ θ + π \theta \mapsto \theta + \pi θ ↦ θ + π と置換すると cos ( θ + π ) = − cos θ \cos(\theta+\pi) = -\cos\theta cos ( θ + π ) = − cos θ で、被積分関数は周期 2 π 2\pi 2 π なので積分区間は [ 0 , 2 π ] [0,2\pi] [ 0 , 2 π ] のままでよく、b b b を − b -b − b に取り替えた式に帰着します。右辺は b b b の符号によらないので、以下 b > 0 b > 0 b > 0 とします。
z = e i θ z = e^{i\theta} z = e i θ と置く上の対応により
I = ∮ ∣ z ∣ = 1 1 a + b 2 ( z + 1 z ) ⋅ d z i z = ∮ ∣ z ∣ = 1 2 d z i ( b z 2 + 2 a z + b ) . I = \oint_{|z|=1}\frac{1}{a + \frac{b}{2}\left(z + \frac1z\right)}\cdot\frac{dz}{iz}
= \oint_{|z|=1}\frac{2\,dz}{i\left(bz^2 + 2az + b\right)} . I = ∮ ∣ z ∣ = 1 a + 2 b ( z + z 1 ) 1 ⋅ i z d z = ∮ ∣ z ∣ = 1 i ( b z 2 + 2 a z + b ) 2 d z . 分母の 2 次式 b z 2 + 2 a z + b bz^2+2az+b b z 2 + 2 a z + b の零点は
z ± = − a ± a 2 − b 2 b z_{\pm} = \frac{-a \pm \sqrt{a^2-b^2}}{b} z ± = b − a ± a 2 − b 2 で、a > b > 0 a > b > 0 a > b > 0 より a 2 − b 2 > 0 a^2 - b^2 > 0 a 2 − b 2 > 0 なのでどちらも実数、かつ相異なります。解と係数の関係から z + z − = b / b = 1 z_+z_- = b/b = 1 z + z − = b / b = 1 です。
∣ z + ∣ < 1 |z_+| < 1 ∣ z + ∣ < 1 を確かめます。z + = − a + a 2 − b 2 b z_+ = \dfrac{-a+\sqrt{a^2-b^2}}{b} z + = b − a + a 2 − b 2 は a 2 − b 2 < a \sqrt{a^2-b^2} < a a 2 − b 2 < a より負で、∣ z + ∣ = a − a 2 − b 2 b |z_+| = \dfrac{a - \sqrt{a^2-b^2}}{b} ∣ z + ∣ = b a − a 2 − b 2 です。∣ z + ∣ < 1 ⟺ a − b < a 2 − b 2 |z_+| < 1 \iff a - b < \sqrt{a^2-b^2} ∣ z + ∣ < 1 ⟺ a − b < a 2 − b 2 であり、両辺とも正(a > b a>b a > b )なので 2 乗して ( a − b ) 2 < ( a − b ) ( a + b ) ⟺ a − b < a + b ⟺ b > 0 \ (a-b)^2 < (a-b)(a+b) \iff a-b < a+b \iff b > 0 ( a − b ) 2 < ( a − b ) ( a + b ) ⟺ a − b < a + b ⟺ b > 0 。これは仮定です。よって ∣ z + ∣ < 1 |z_+|<1 ∣ z + ∣ < 1 、そして ∣ z − ∣ = 1 / ∣ z + ∣ > 1 |z_-| = 1/|z_+| > 1 ∣ z − ∣ = 1/∣ z + ∣ > 1 です。特に単位円周上に零点はありません。
したがって単位円の内部にある特異点は単純極 z + z_+ z + だけです。g ( z ) = 2 i ( b z 2 + 2 a z + b ) g(z) = \dfrac{2}{i(bz^2+2az+b)} g ( z ) = i ( b z 2 + 2 a z + b ) 2 とすると、Proposition 4.2 (2) より
Res ( g , z + ) = lim z → z + 2 ( z − z + ) i b ( z − z + ) ( z − z − ) = 2 i b ( z + − z − ) = 2 i b ⋅ 2 a 2 − b 2 b = 1 i a 2 − b 2 . \operatorname{Res}(g, z_+) = \lim_{z\to z_+}\frac{2(z-z_+)}{i\,b(z-z_+)(z-z_-)}
= \frac{2}{i\,b(z_+-z_-)}
= \frac{2}{i\,b\cdot\frac{2\sqrt{a^2-b^2}}{b}}
= \frac{1}{i\sqrt{a^2-b^2}} . Res ( g , z + ) = z → z + lim i b ( z − z + ) ( z − z − ) 2 ( z − z + ) = i b ( z + − z − ) 2 = i b ⋅ b 2 a 2 − b 2 2 = i a 2 − b 2 1 . Theorem 5.1 より I = 2 π i ⋅ 1 i a 2 − b 2 = 2 π a 2 − b 2 I = 2\pi i \cdot \dfrac{1}{i\sqrt{a^2-b^2}} = \dfrac{2\pi}{\sqrt{a^2-b^2}} I = 2 π i ⋅ i a 2 − b 2 1 = a 2 − b 2 2 π です。
∎ たとえば ∫ 0 2 π d θ 5 + 4 cos θ = 2 π 25 − 16 = 2 π 3 \displaystyle\int_0^{2\pi}\frac{d\theta}{5+4\cos\theta} = \frac{2\pi}{\sqrt{25-16}} = \frac{2\pi}{3} ∫ 0 2 π 5 + 4 cos θ d θ = 25 − 16 2 π = 3 2 π です。次は分子に cos 2 θ \cos 2\theta cos 2 θ が入る場合で、0 0 0 に多重極が現れます。
Example 7.2 (分子に cos 2θ がある積分 )
I = ∫ 0 2 π cos 2 θ 5 − 4 cos θ d θ \displaystyle I = \int_0^{2\pi}\frac{\cos2\theta}{5-4\cos\theta}\,d\theta I = ∫ 0 2 π 5 − 4 cos θ cos 2 θ d θ を計算します。5 − 4 cos θ ≥ 1 > 0 5 - 4\cos\theta \ge 1 > 0 5 − 4 cos θ ≥ 1 > 0 なので被積分関数は連続です。z = e i θ z = e^{i\theta} z = e i θ とすると cos 2 θ = z 2 + z − 2 2 \cos2\theta = \frac{z^2+z^{-2}}{2} cos 2 θ = 2 z 2 + z − 2 、5 − 4 cos θ = 5 − 2 ( z + z − 1 ) = − 2 z 2 + 5 z − 2 z 5-4\cos\theta = 5 - 2(z+z^{-1}) = \dfrac{-2z^2+5z-2}{z} 5 − 4 cos θ = 5 − 2 ( z + z − 1 ) = z − 2 z 2 + 5 z − 2 なので
I = ∮ ∣ z ∣ = 1 ( z 2 + z − 2 ) / 2 ( − 2 z 2 + 5 z − 2 ) / z ⋅ d z i z = ∮ ∣ z ∣ = 1 z 4 + 1 2 z 2 i ( − 2 z 2 + 5 z − 2 ) d z = i 2 ∮ ∣ z ∣ = 1 z 4 + 1 z 2 ( 2 z − 1 ) ( z − 2 ) d z I = \oint_{|z|=1}\frac{(z^2+z^{-2})/2}{(-2z^2+5z-2)/z}\cdot\frac{dz}{iz}
= \oint_{|z|=1}\frac{z^4+1}{2z^2\,i\,(-2z^2+5z-2)}\,dz
= \frac{i}{2}\oint_{|z|=1}\frac{z^4+1}{z^2(2z-1)(z-2)}\,dz I = ∮ ∣ z ∣ = 1 ( − 2 z 2 + 5 z − 2 ) / z ( z 2 + z − 2 ) /2 ⋅ i z d z = ∮ ∣ z ∣ = 1 2 z 2 i ( − 2 z 2 + 5 z − 2 ) z 4 + 1 d z = 2 i ∮ ∣ z ∣ = 1 z 2 ( 2 z − 1 ) ( z − 2 ) z 4 + 1 d z です(− 2 z 2 + 5 z − 2 = − ( 2 z − 1 ) ( z − 2 ) -2z^2+5z-2 = -(2z-1)(z-2) − 2 z 2 + 5 z − 2 = − ( 2 z − 1 ) ( z − 2 ) と 1 − 2 i = i 2 \frac{1}{-2i} = \frac{i}{2} − 2 i 1 = 2 i を使いました)。被積分関数を g g g とすると、単位円の内部の特異点は z = 0 z=0 z = 0 (2 2 2 位の極)と z = 1 / 2 z=1/2 z = 1/2 (単純極)で、z = 2 z=2 z = 2 は外部です。
z = 1 / 2 z = 1/2 z = 1/2 :Proposition 4.2 (2) より
Res ( g , 1 2 ) = z 4 + 1 z 2 ⋅ 2 ( z − 2 ) ∣ z = 1 / 2 = 1 16 + 1 1 4 ⋅ 2 ( 1 2 − 2 ) = 17 / 16 − 3 / 4 = − 17 12 . \operatorname{Res}(g,\tfrac12) = \frac{z^4+1}{z^2\cdot 2(z-2)}\bigg|_{z=1/2}
= \frac{\frac{1}{16}+1}{\frac14 \cdot 2\left(\frac12-2\right)}
= \frac{17/16}{-3/4} = -\frac{17}{12}. Res ( g , 2 1 ) = z 2 ⋅ 2 ( z − 2 ) z 4 + 1 z = 1/2 = 4 1 ⋅ 2 ( 2 1 − 2 ) 16 1 + 1 = − 3/4 17/16 = − 12 17 . z = 0 z = 0 z = 0 :m = 2 m=2 m = 2 として φ ( z ) = z 2 g ( z ) = z 4 + 1 2 z 2 − 5 z + 2 \varphi(z) = z^2 g(z) = \dfrac{z^4+1}{2z^2-5z+2} φ ( z ) = z 2 g ( z ) = 2 z 2 − 5 z + 2 z 4 + 1 とおくと
φ ′ ( z ) = 4 z 3 ( 2 z 2 − 5 z + 2 ) − ( z 4 + 1 ) ( 4 z − 5 ) ( 2 z 2 − 5 z + 2 ) 2 , φ ′ ( 0 ) = 0 − 1 ⋅ ( − 5 ) 4 = 5 4 , \varphi'(z) = \frac{4z^3(2z^2-5z+2) - (z^4+1)(4z-5)}{(2z^2-5z+2)^2},
\qquad
\varphi'(0) = \frac{0 - 1\cdot(-5)}{4} = \frac{5}{4}, φ ′ ( z ) = ( 2 z 2 − 5 z + 2 ) 2 4 z 3 ( 2 z 2 − 5 z + 2 ) − ( z 4 + 1 ) ( 4 z − 5 ) , φ ′ ( 0 ) = 4 0 − 1 ⋅ ( − 5 ) = 4 5 , なので Res ( g , 0 ) = φ ′ ( 0 ) / 1 ! = 5 / 4 \operatorname{Res}(g,0) = \varphi'(0)/1! = 5/4 Res ( g , 0 ) = φ ′ ( 0 ) /1 ! = 5/4 です。念のためローラン展開でも確かめます。1 2 − 5 z + 2 z 2 = 1 2 ⋅ 1 1 − 5 z − 2 z 2 2 = 1 2 ( 1 + 5 z 2 + O ( z 2 ) ) = 1 2 + 5 4 z + O ( z 2 ) \dfrac{1}{2-5z+2z^2} = \dfrac{1}{2}\cdot\dfrac{1}{1 - \frac{5z-2z^2}{2}} = \dfrac12\left(1 + \frac{5z}{2} + O(z^2)\right) = \frac12 + \frac{5}{4}z + O(z^2) 2 − 5 z + 2 z 2 1 = 2 1 ⋅ 1 − 2 5 z − 2 z 2 1 = 2 1 ( 1 + 2 5 z + O ( z 2 ) ) = 2 1 + 4 5 z + O ( z 2 ) なので、( z 4 + 1 ) (z^4+1) ( z 4 + 1 ) を掛けて z 2 z^2 z 2 で割ると g ( z ) = 1 2 z 2 + 5 4 z + ⋯ g(z) = \dfrac{1}{2z^2} + \dfrac{5}{4z} + \cdots g ( z ) = 2 z 2 1 + 4 z 5 + ⋯ となり、z − 1 z^{-1} z − 1 の係数は 5 / 4 5/4 5/4 で一致します。
したがって留数の和は − 17 12 + 5 4 = − 17 12 + 15 12 = − 1 6 -\frac{17}{12} + \frac{5}{4} = -\frac{17}{12}+\frac{15}{12} = -\frac16 − 12 17 + 4 5 = − 12 17 + 12 15 = − 6 1 で、Theorem 5.1 より
I = i 2 ⋅ 2 π i ⋅ ( − 1 6 ) = − π ⋅ ( − 1 6 ) = π 6 . I = \frac{i}{2}\cdot 2\pi i \cdot\left(-\frac16\right) = -\pi\cdot\left(-\frac16\right) = \frac{\pi}{6}. I = 2 i ⋅ 2 π i ⋅ ( − 6 1 ) = − π ⋅ ( − 6 1 ) = 6 π . ∫ − ∞ ∞ f ( x ) cos λ x d x \int_{-\infty}^{\infty}f(x)\cos\lambda x\,dx ∫ − ∞ ∞ f ( x ) cos λ x d x のような積分では、cos λ z \cos\lambda z cos λ z をそのまま複素平面に延長してはいけません。上半平面では cos ( i y ) = cosh y \cos(iy) = \cosh y cos ( i y ) = cosh y が y → ∞ y \to \infty y → ∞ で指数的に発散するからです。代わりに e i λ z e^{i\lambda z} e iλ z を使います。これは上半平面で ∣ e i λ z ∣ = e − λ Im z ≤ 1 |e^{i\lambda z}| = e^{-\lambda\operatorname{Im}z} \le 1 ∣ e iλ z ∣ = e − λ Im z ≤ 1 (λ > 0 \lambda>0 λ > 0 )と抑えられ、しかも実軸上では cos λ x + i sin λ x \cos\lambda x + i\sin\lambda x cos λ x + i sin λ x なので、最後に実部と虚部を取れば目的の積分が得られます。
Example 8.1 (ローレンツ型のフーリエ変換 )
∫ − ∞ ∞ cos x x 2 + 1 d x \displaystyle\int_{-\infty}^{\infty}\frac{\cos x}{x^2+1}\,dx ∫ − ∞ ∞ x 2 + 1 cos x d x を求めます。f ( z ) = 1 z 2 + 1 f(z) = \dfrac{1}{z^2+1} f ( z ) = z 2 + 1 1 、λ = 1 \lambda = 1 λ = 1 として F ( z ) = f ( z ) e i z F(z) = f(z)e^{iz} F ( z ) = f ( z ) e i z を Example 6.3 と同じ半円輪郭 C R C_R C R (R > 1 R>1 R > 1 )で積分します。F F F の特異点は z = ± i z = \pm i z = ± i の単純極で、上半平面にあるのは z = i z = i z = i だけです。Proposition 4.2 (3) を g ( z ) = e i z g(z) = e^{iz} g ( z ) = e i z 、h ( z ) = z 2 + 1 h(z)=z^2+1 h ( z ) = z 2 + 1 、h ′ ( z ) = 2 z h'(z) = 2z h ′ ( z ) = 2 z に適用して
Res ( F , i ) = e i ⋅ i 2 i = e − 1 2 i . \operatorname{Res}(F, i) = \frac{e^{i\cdot i}}{2i} = \frac{e^{-1}}{2i}. Res ( F , i ) = 2 i e i ⋅ i = 2 i e − 1 . Theorem 5.1 より ∫ − R R e i x x 2 + 1 d x + ∫ Γ R F ( z ) d z = 2 π i ⋅ e − 1 2 i = π e \displaystyle\int_{-R}^{R}\frac{e^{ix}}{x^2+1}dx + \int_{\Gamma_R}F(z)\,dz = 2\pi i\cdot\frac{e^{-1}}{2i} = \frac{\pi}{e} ∫ − R R x 2 + 1 e i x d x + ∫ Γ R F ( z ) d z = 2 π i ⋅ 2 i e − 1 = e π です。
円弧の評価には Lemma 6.1 (ii) を使います。∣ z ∣ = R > 1 |z| = R > 1 ∣ z ∣ = R > 1 のとき ∣ f ( z ) ∣ = 1 ∣ z 2 + 1 ∣ ≤ 1 R 2 − 1 |f(z)| = \dfrac{1}{|z^2+1|} \le \dfrac{1}{R^2-1} ∣ f ( z ) ∣ = ∣ z 2 + 1∣ 1 ≤ R 2 − 1 1 (三角不等式 ∣ z 2 + 1 ∣ ≥ ∣ z ∣ 2 − 1 |z^2+1| \ge |z|^2 - 1 ∣ z 2 + 1∣ ≥ ∣ z ∣ 2 − 1 )なので M ( R ) ≤ 1 R 2 − 1 → 0 M(R) \le \dfrac{1}{R^2-1}\to0 M ( R ) ≤ R 2 − 1 1 → 0 であり、∣ ∫ Γ R F ∣ ≤ π 1 ⋅ 1 R 2 − 1 → 0 \left|\int_{\Gamma_R}F\right| \le \dfrac{\pi}{1}\cdot\dfrac{1}{R^2-1}\to0 ∫ Γ R F ≤ 1 π ⋅ R 2 − 1 1 → 0 です。よって
∫ − ∞ ∞ e i x x 2 + 1 d x = π e . \int_{-\infty}^{\infty}\frac{e^{ix}}{x^2+1}\,dx = \frac{\pi}{e} . ∫ − ∞ ∞ x 2 + 1 e i x d x = e π . 実部を取ると ∫ − ∞ ∞ cos x x 2 + 1 d x = π e \displaystyle\int_{-\infty}^{\infty}\frac{\cos x}{x^2+1}dx = \frac{\pi}{e} ∫ − ∞ ∞ x 2 + 1 cos x d x = e π 、虚部を取ると ∫ − ∞ ∞ sin x x 2 + 1 d x = 0 \displaystyle\int_{-\infty}^{\infty}\frac{\sin x}{x^2+1}dx = 0 ∫ − ∞ ∞ x 2 + 1 sin x d x = 0 です。後者は被積分関数が奇関数であることからも分かり、検算になります。なお cos x x 2 + 1 \dfrac{\cos x}{x^2+1} x 2 + 1 cos x の絶対値は 1 x 2 + 1 \dfrac{1}{x^2+1} x 2 + 1 1 以下なので、この広義積分は絶対収束しています。
被積分関数が実軸上に極を持つと、広義積分は普通の意味では発散します。それでも「対称に切り取れば」極限が存在することがあります。
Definition 8.2 (コーシーの主値 )
c ∈ R c \in \mathbb{R} c ∈ R とし、f f f は R ∖ { c } \mathbb{R}\setminus\{c\} R ∖ { c } で連続とする。極限
p . v . ∫ − ∞ ∞ f ( x ) d x : = lim R → ∞ ε → + 0 ( ∫ − R c − ε f ( x ) d x + ∫ c + ε R f ( x ) d x ) \mathrm{p.v.}\int_{-\infty}^{\infty}f(x)\,dx
:= \lim_{\substack{R \to \infty \\ \varepsilon \to +0}}
\left(\int_{-R}^{c-\varepsilon}f(x)\,dx + \int_{c+\varepsilon}^{R}f(x)\,dx\right) p.v. ∫ − ∞ ∞ f ( x ) d x := R → ∞ ε → + 0 lim ( ∫ − R c − ε f ( x ) d x + ∫ c + ε R f ( x ) d x ) が存在するとき、その値を f f f のコーシーの主値 と呼ぶ。左右を同じ ε \varepsilon ε で切り落とすのが要点である。通常の広義積分が存在するときは、主値も存在して同じ値になる。
輪郭のほうは、実軸上の極を半径 ε \varepsilon ε の小円弧で迂回して避けます。その小円弧の寄与を測るのが次の補題です。「一周すれば 2 π i × 2\pi i \times 2 π i × 留数」だったものが、「角度 β − α \beta-\alpha β − α だけ回れば i ( β − α ) × i(\beta-\alpha)\times i ( β − α ) × 留数」になる、と憶えると自然です。
Lemma 8.3 (単純極まわりの小円弧(くぼみ補題) )
c c c を f f f の単純極とし、ρ = Res ( f , c ) \rho = \operatorname{Res}(f,c) ρ = Res ( f , c ) とおく。実数 α , β \alpha, \beta α , β に対し、γ ε \gamma_\varepsilon γ ε を z = c + ε e i θ z = c + \varepsilon e^{i\theta} z = c + ε e i θ (θ \theta θ が α \alpha α から β \beta β まで動く)で定まる円弧とする。このとき
lim ε → + 0 ∫ γ ε f ( z ) d z = i ( β − α ) ρ . \lim_{\varepsilon \to +0}\int_{\gamma_\varepsilon}f(z)\,dz = i(\beta - \alpha)\rho . ε → + 0 lim ∫ γ ε f ( z ) d z = i ( β − α ) ρ . (β > α \beta > \alpha β > α なら反時計回り、β < α \beta < \alpha β < α なら時計回りである。)
Proof(Lemma 8.3) c c c は単純極なので、ローラン展開は f ( z ) = ρ z − c + ∑ n ≥ 0 c n ( z − c ) n f(z) = \dfrac{\rho}{z-c} + \sum_{n \ge 0}c_n(z-c)^n f ( z ) = z − c ρ + ∑ n ≥ 0 c n ( z − c ) n の形です。第 2 項を g ( z ) g(z) g ( z ) とおくと、g g g は D ( c , R ) D(c,R) D ( c , R ) で正則なべき級数の和ですから、ある r > 0 r>0 r > 0 と M > 0 M>0 M > 0 に対し ∣ z − c ∣ ≤ r |z-c| \le r ∣ z − c ∣ ≤ r で ∣ g ( z ) ∣ ≤ M |g(z)| \le M ∣ g ( z ) ∣ ≤ M です(連続関数のコンパクト集合上の有界性)。
主要部の寄与は ε \varepsilon ε によらず計算できます。z = c + ε e i θ z = c+\varepsilon e^{i\theta} z = c + ε e i θ 、d z = i ε e i θ d θ dz = i\varepsilon e^{i\theta}d\theta d z = i ε e i θ d θ より
∫ γ ε ρ z − c d z = ρ ∫ α β i ε e i θ ε e i θ d θ = i ρ ( β − α ) . \int_{\gamma_\varepsilon}\frac{\rho}{z-c}\,dz
= \rho\int_\alpha^\beta \frac{i\varepsilon e^{i\theta}}{\varepsilon e^{i\theta}}\,d\theta
= i\rho(\beta-\alpha). ∫ γ ε z − c ρ d z = ρ ∫ α β ε e i θ i ε e i θ d θ = i ρ ( β − α ) . 残りは ML 不等式で消えます。γ ε \gamma_\varepsilon γ ε の長さは ∣ β − α ∣ ε |\beta-\alpha|\varepsilon ∣ β − α ∣ ε なので、ε ≤ r \varepsilon \le r ε ≤ r のとき
∣ ∫ γ ε g ( z ) d z ∣ ≤ ∣ β − α ∣ ε M ⟶ 0 ( ε → + 0 ) . \left|\int_{\gamma_\varepsilon}g(z)\,dz\right| \le |\beta-\alpha|\,\varepsilon\, M \longrightarrow 0
\qquad (\varepsilon \to +0). ∫ γ ε g ( z ) d z ≤ ∣ β − α ∣ ε M ⟶ 0 ( ε → + 0 ) . 両者を足して結論を得ます。
∎ Example 8.4 (ディリクレ積分 )
∫ 0 ∞ sin x x d x = π 2 \displaystyle\int_0^{\infty}\frac{\sin x}{x}\,dx = \frac{\pi}{2} ∫ 0 ∞ x sin x d x = 2 π を示します。f ( z ) = e i z z f(z) = \dfrac{e^{iz}}{z} f ( z ) = z e i z とし、0 < ε < R 0 < \varepsilon < R 0 < ε < R に対して次の四つをつないだ輪郭 C C C を考えます。線分 [ ε , R ] [\varepsilon, R] [ ε , R ] 、大円弧 Γ R \Gamma_R Γ R (反時計回りに 0 0 0 から π \pi π へ)、線分 [ − R , − ε ] [-R,-\varepsilon] [ − R , − ε ] 、小円弧 γ ε \gamma_\varepsilon γ ε (θ \theta θ が π \pi π から 0 0 0 へ、すなわち時計回り)。
C C C とその内部は、複素平面から半直線 { i y : y ≤ 0 } \{iy : y \le 0\} { i y : y ≤ 0 } を除いた単連結領域 D D D に含まれ、f f f は D D D で正則です(0 0 0 は D D D に属しません)。よって Theorem 5.1 を特異点 0 0 0 個の場合に使って(これはコーシーの積分定理そのものです)
∫ ε R e i x x d x + ∫ Γ R f d z + ∫ − R − ε e i x x d x + ∫ γ ε f d z = 0. \int_\varepsilon^R \frac{e^{ix}}{x}dx + \int_{\Gamma_R}f\,dz + \int_{-R}^{-\varepsilon}\frac{e^{ix}}{x}dx + \int_{\gamma_\varepsilon}f\,dz = 0 . ∫ ε R x e i x d x + ∫ Γ R f d z + ∫ − R − ε x e i x d x + ∫ γ ε f d z = 0. 第 3 項は x = − t x = -t x = − t と置換すると ∫ − R − ε e i x x d x = ∫ R ε e − i t − t ( − d t ) = − ∫ ε R e − i t t d t \displaystyle\int_{-R}^{-\varepsilon}\frac{e^{ix}}{x}dx = \int_{R}^{\varepsilon}\frac{e^{-it}}{-t}(-dt) = -\int_\varepsilon^R\frac{e^{-it}}{t}dt ∫ − R − ε x e i x d x = ∫ R ε − t e − i t ( − d t ) = − ∫ ε R t e − i t d t なので、第 1 項と合わせて
∫ ε R e i x − e − i x x d x = 2 i ∫ ε R sin x x d x \int_\varepsilon^R\frac{e^{ix}-e^{-ix}}{x}\,dx = 2i\int_\varepsilon^R\frac{\sin x}{x}\,dx ∫ ε R x e i x − e − i x d x = 2 i ∫ ε R x sin x d x になります。第 2 項は Lemma 6.1 (ii) で f ( z ) = 1 / z f(z) = 1/z f ( z ) = 1/ z 、λ = 1 \lambda = 1 λ = 1 とすると M ( R ) = 1 / R → 0 M(R) = 1/R \to 0 M ( R ) = 1/ R → 0 なので 0 0 0 に収束します。第 4 項は、0 0 0 が f f f の単純極で Res ( f , 0 ) = lim z → 0 z ⋅ e i z z = 1 \operatorname{Res}(f,0) = \lim_{z\to0}z\cdot\frac{e^{iz}}{z} = 1 Res ( f , 0 ) = lim z → 0 z ⋅ z e i z = 1 ですから、Lemma 8.3 を α = π \alpha = \pi α = π 、β = 0 \beta = 0 β = 0 として適用し i ( 0 − π ) ⋅ 1 = − i π i(0-\pi)\cdot1 = -i\pi i ( 0 − π ) ⋅ 1 = − iπ に収束します。
以上を合わせて ε → + 0 \varepsilon \to +0 ε → + 0 、R → ∞ R \to \infty R → ∞ とすると
2 i ∫ 0 ∞ sin x x d x − i π = 0 , すなわち ∫ 0 ∞ sin x x d x = π 2 . 2i\int_0^{\infty}\frac{\sin x}{x}\,dx - i\pi = 0,
\qquad\text{すなわち}\qquad
\int_0^{\infty}\frac{\sin x}{x}\,dx = \frac{\pi}{2}. 2 i ∫ 0 ∞ x sin x d x − iπ = 0 , すなわち ∫ 0 ∞ x sin x d x = 2 π . 同じ計算を実部・虚部に分けて読むと p . v . ∫ − ∞ ∞ e i x x d x = i π \mathrm{p.v.}\displaystyle\int_{-\infty}^{\infty}\frac{e^{ix}}{x}dx = i\pi p.v. ∫ − ∞ ∞ x e i x d x = iπ であり、実部から p . v . ∫ − ∞ ∞ cos x x d x = 0 \mathrm{p.v.}\displaystyle\int_{-\infty}^{\infty}\frac{\cos x}{x}dx = 0 p.v. ∫ − ∞ ∞ x cos x d x = 0 (奇関数なので当然です)、虚部から上の結果が出ます。
注意すべき点が二つあります。第一に、sin x x \dfrac{\sin x}{x} x sin x は x = 0 x=0 x = 0 で除去可能特異点を持つ(Example 3.6 (a))ので ε → 0 \varepsilon \to 0 ε → 0 の極限は問題なく、ε \varepsilon ε を入れたのは e i z / z e^{iz}/z e i z / z を扱うための便宜です。第二に、この広義積分は絶対収束しません(Proposition 2.2[ボールウェイン積分] )。実際 ∫ 0 ∞ ∣ sin x x ∣ d x = ∞ \int_0^{\infty}\left|\frac{\sin x}{x}\right|dx = \infty ∫ 0 ∞ x s i n x d x = ∞ であり、値は R → ∞ R \to \infty R → ∞ の極限としてのみ意味を持ちます。実解析の側からの別証明は Theorem 3.2[ボールウェイン積分] にあり、この繊細さが表に出る有名な例として ボールウェイン積分 があります。
Exercise 9.1 易
f ( z ) = e z z 2 ( z − 1 ) f(z) = \dfrac{e^z}{z^2(z-1)} f ( z ) = z 2 ( z − 1 ) e z のすべての孤立特異点を求め、それぞれの位数と留数を求めてください。
Solution 特異点は分母の零点 z = 0 z=0 z = 0 (2 2 2 位)と z = 1 z=1 z = 1 (1 1 1 位)です。e z e^z e z はこれらの点で 0 0 0 にならないので、位数は分母のとおりです。実際 z 2 f ( z ) = e z z − 1 → − 1 ≠ 0 z^2f(z) = \frac{e^z}{z-1} \to -1 \ne 0 z 2 f ( z ) = z − 1 e z → − 1 = 0 、( z − 1 ) f ( z ) = e z z 2 → e ≠ 0 (z-1)f(z) = \frac{e^z}{z^2}\to e \ne 0 ( z − 1 ) f ( z ) = z 2 e z → e = 0 なので Proposition 3.4 より 0 0 0 は 2 2 2 位の極、1 1 1 は単純極です。
z = 1 z=1 z = 1 :Proposition 4.2 (2) より Res ( f , 1 ) = lim z → 1 e z z 2 = e \operatorname{Res}(f,1) = \lim_{z\to1}\frac{e^z}{z^2} = e Res ( f , 1 ) = lim z → 1 z 2 e z = e 。
z = 0 z=0 z = 0 :同 (1) で m = 2 m=2 m = 2 とし、φ ( z ) = e z z − 1 \varphi(z) = \dfrac{e^z}{z-1} φ ( z ) = z − 1 e z に対し
φ ′ ( z ) = e z ( z − 1 ) − e z ( z − 1 ) 2 = e z ( z − 2 ) ( z − 1 ) 2 , Res ( f , 0 ) = φ ′ ( 0 ) = 1 ⋅ ( − 2 ) 1 = − 2. \varphi'(z) = \frac{e^z(z-1) - e^z}{(z-1)^2} = \frac{e^z(z-2)}{(z-1)^2},
\qquad
\operatorname{Res}(f,0) = \varphi'(0) = \frac{1\cdot(-2)}{1} = -2 . φ ′ ( z ) = ( z − 1 ) 2 e z ( z − 1 ) − e z = ( z − 1 ) 2 e z ( z − 2 ) , Res ( f , 0 ) = φ ′ ( 0 ) = 1 1 ⋅ ( − 2 ) = − 2. したがって留数の和は e − 2 e - 2 e − 2 です。
Exercise 9.2 標準
∫ 0 ∞ x 2 ( x 2 + 1 ) ( x 2 + 4 ) d x \displaystyle\int_0^{\infty}\frac{x^2}{(x^2+1)(x^2+4)}\,dx ∫ 0 ∞ ( x 2 + 1 ) ( x 2 + 4 ) x 2 d x を留数定理で計算し、部分分数分解による初等的な計算と一致することを確かめてください。
Solution 被積分関数は偶関数なので ∫ 0 ∞ = 1 2 ∫ − ∞ ∞ \displaystyle\int_0^{\infty} = \frac12\int_{-\infty}^{\infty} ∫ 0 ∞ = 2 1 ∫ − ∞ ∞ です。f ( z ) = z 2 ( z 2 + 1 ) ( z 2 + 4 ) f(z) = \dfrac{z^2}{(z^2+1)(z^2+4)} f ( z ) = ( z 2 + 1 ) ( z 2 + 4 ) z 2 は deg \deg deg の差が 2 2 2 、分母は実軸上で 0 0 0 にならないので Proposition 6.2 が使えます。上半平面の極は単純極 z = i z=i z = i と z = 2 i z=2i z = 2 i です。
Res ( f , i ) = z 2 ( z + i ) ( z 2 + 4 ) ∣ z = i = − 1 2 i ⋅ 3 = − 1 6 i = i 6 , \operatorname{Res}(f,i) = \frac{z^2}{(z+i)(z^2+4)}\bigg|_{z=i} = \frac{-1}{2i\cdot3} = \frac{-1}{6i} = \frac{i}{6}, Res ( f , i ) = ( z + i ) ( z 2 + 4 ) z 2 z = i = 2 i ⋅ 3 − 1 = 6 i − 1 = 6 i , Res ( f , 2 i ) = z 2 ( z 2 + 1 ) ( z + 2 i ) ∣ z = 2 i = − 4 ( − 3 ) ⋅ 4 i = 1 3 i = − i 3 . \operatorname{Res}(f,2i) = \frac{z^2}{(z^2+1)(z+2i)}\bigg|_{z=2i} = \frac{-4}{(-3)\cdot4i} = \frac{1}{3i} = -\frac{i}{3}. Res ( f , 2 i ) = ( z 2 + 1 ) ( z + 2 i ) z 2 z = 2 i = ( − 3 ) ⋅ 4 i − 4 = 3 i 1 = − 3 i . (1 i = − i \frac1i = -i i 1 = − i を使いました。)和は i 6 − i 3 = − i 6 \frac{i}{6}-\frac{i}{3} = -\frac{i}{6} 6 i − 3 i = − 6 i なので
∫ − ∞ ∞ f ( x ) d x = 2 π i ( − i 6 ) = π 3 , ∫ 0 ∞ f ( x ) d x = π 6 . \int_{-\infty}^{\infty}f(x)\,dx = 2\pi i\left(-\frac{i}{6}\right) = \frac{\pi}{3},
\qquad
\int_0^{\infty}f(x)\,dx = \frac{\pi}{6}. ∫ − ∞ ∞ f ( x ) d x = 2 π i ( − 6 i ) = 3 π , ∫ 0 ∞ f ( x ) d x = 6 π . 初等的な確認:x 2 ( x 2 + 1 ) ( x 2 + 4 ) = A x 2 + 1 + B x 2 + 4 \dfrac{x^2}{(x^2+1)(x^2+4)} = \dfrac{A}{x^2+1}+\dfrac{B}{x^2+4} ( x 2 + 1 ) ( x 2 + 4 ) x 2 = x 2 + 1 A + x 2 + 4 B とおくと x 2 = A ( x 2 + 4 ) + B ( x 2 + 1 ) x^2 = A(x^2+4)+B(x^2+1) x 2 = A ( x 2 + 4 ) + B ( x 2 + 1 ) より A + B = 1 A+B=1 A + B = 1 、4 A + B = 0 4A+B=0 4 A + B = 0 、すなわち A = − 1 3 A = -\frac13 A = − 3 1 、B = 4 3 B=\frac43 B = 3 4 です。∫ 0 ∞ d x x 2 + a 2 = π 2 a \int_0^{\infty}\frac{dx}{x^2+a^2} = \frac{\pi}{2a} ∫ 0 ∞ x 2 + a 2 d x = 2 a π (a > 0 a>0 a > 0 )を使うと
− 1 3 ⋅ π 2 + 4 3 ⋅ π 4 = − π 6 + π 3 = π 6 -\frac13\cdot\frac{\pi}{2} + \frac43\cdot\frac{\pi}{4} = -\frac{\pi}{6}+\frac{\pi}{3} = \frac{\pi}{6} − 3 1 ⋅ 2 π + 3 4 ⋅ 4 π = − 6 π + 3 π = 6 π で一致します。
Exercise 9.3 標準
∫ 0 2 π d θ ( 2 + cos θ ) 2 \displaystyle\int_0^{2\pi}\frac{d\theta}{(2+\cos\theta)^2} ∫ 0 2 π ( 2 + cos θ ) 2 d θ を計算してください。
Solution z = e i θ z = e^{i\theta} z = e i θ と置くと 2 + cos θ = 2 + z + z − 1 2 = z 2 + 4 z + 1 2 z 2+\cos\theta = 2 + \dfrac{z+z^{-1}}{2} = \dfrac{z^2+4z+1}{2z} 2 + cos θ = 2 + 2 z + z − 1 = 2 z z 2 + 4 z + 1 なので
I = ∮ ∣ z ∣ = 1 4 z 2 ( z 2 + 4 z + 1 ) 2 ⋅ d z i z = 4 i ∮ ∣ z ∣ = 1 z d z ( z 2 + 4 z + 1 ) 2 . I = \oint_{|z|=1}\frac{4z^2}{(z^2+4z+1)^2}\cdot\frac{dz}{iz}
= \frac{4}{i}\oint_{|z|=1}\frac{z\,dz}{(z^2+4z+1)^2}. I = ∮ ∣ z ∣ = 1 ( z 2 + 4 z + 1 ) 2 4 z 2 ⋅ i z d z = i 4 ∮ ∣ z ∣ = 1 ( z 2 + 4 z + 1 ) 2 z d z . z 2 + 4 z + 1 = 0 z^2+4z+1 = 0 z 2 + 4 z + 1 = 0 の解は z 0 = − 2 + 3 ≈ − 0.268 z_0 = -2+\sqrt3 \approx -0.268 z 0 = − 2 + 3 ≈ − 0.268 と z 1 = − 2 − 3 ≈ − 3.73 z_1 = -2-\sqrt3 \approx -3.73 z 1 = − 2 − 3 ≈ − 3.73 で、単位円の内部にあるのは z 0 z_0 z 0 だけです(z 0 z 1 = 1 z_0z_1 = 1 z 0 z 1 = 1 なので一方が内側、他方が外側)。分母は ( z − z 0 ) 2 ( z − z 1 ) 2 (z-z_0)^2(z-z_1)^2 ( z − z 0 ) 2 ( z − z 1 ) 2 なので z 0 z_0 z 0 は 2 2 2 位の極です。m = 2 m=2 m = 2 として Proposition 4.2 (1) を使います。φ ( z ) = z ( z − z 1 ) 2 \varphi(z) = \dfrac{z}{(z-z_1)^2} φ ( z ) = ( z − z 1 ) 2 z とすると
φ ′ ( z ) = ( z − z 1 ) 2 − z ⋅ 2 ( z − z 1 ) ( z − z 1 ) 4 = ( z − z 1 ) − 2 z ( z − z 1 ) 3 = − z − z 1 ( z − z 1 ) 3 . \varphi'(z) = \frac{(z-z_1)^2 - z\cdot2(z-z_1)}{(z-z_1)^4} = \frac{(z-z_1)-2z}{(z-z_1)^3} = \frac{-z-z_1}{(z-z_1)^3}. φ ′ ( z ) = ( z − z 1 ) 4 ( z − z 1 ) 2 − z ⋅ 2 ( z − z 1 ) = ( z − z 1 ) 3 ( z − z 1 ) − 2 z = ( z − z 1 ) 3 − z − z 1 . z 0 − z 1 = 2 3 z_0 - z_1 = 2\sqrt3 z 0 − z 1 = 2 3 、− z 0 − z 1 = − ( − 2 + 3 ) − ( − 2 − 3 ) = 4 -z_0-z_1 = -(-2+\sqrt3)-(-2-\sqrt3) = 4 − z 0 − z 1 = − ( − 2 + 3 ) − ( − 2 − 3 ) = 4 なので
Res = φ ′ ( z 0 ) = 4 ( 2 3 ) 3 = 4 24 3 = 1 6 3 . \operatorname{Res} = \varphi'(z_0) = \frac{4}{(2\sqrt3)^3} = \frac{4}{24\sqrt3} = \frac{1}{6\sqrt3}. Res = φ ′ ( z 0 ) = ( 2 3 ) 3 4 = 24 3 4 = 6 3 1 . Theorem 5.1 より
I = 4 i ⋅ 2 π i ⋅ 1 6 3 = 8 π 6 3 = 4 π 3 3 = 4 3 π 9 ≈ 2.418. I = \frac{4}{i}\cdot2\pi i\cdot\frac{1}{6\sqrt3} = \frac{8\pi}{6\sqrt3} = \frac{4\pi}{3\sqrt3} = \frac{4\sqrt3\,\pi}{9} \approx 2.418 . I = i 4 ⋅ 2 π i ⋅ 6 3 1 = 6 3 8 π = 3 3 4 π = 9 4 3 π ≈ 2.418. 検算:Proposition 7.1 の式 ∫ 0 2 π d θ a + b cos θ = 2 π ( a 2 − b 2 ) − 1 / 2 \int_0^{2\pi}\frac{d\theta}{a+b\cos\theta} = 2\pi(a^2-b^2)^{-1/2} ∫ 0 2 π a + b c o s θ d θ = 2 π ( a 2 − b 2 ) − 1/2 を a a a で微分すると、左辺は − ∫ 0 2 π d θ ( a + b cos θ ) 2 -\int_0^{2\pi}\frac{d\theta}{(a+b\cos\theta)^2} − ∫ 0 2 π ( a + b c o s θ ) 2 d θ 、右辺は − 2 π a ( a 2 − b 2 ) − 3 / 2 -2\pi a(a^2-b^2)^{-3/2} − 2 π a ( a 2 − b 2 ) − 3/2 です。a = 2 a=2 a = 2 、b = 1 b=1 b = 1 とすると 4 π 3 3 / 2 = 4 π 3 3 \frac{4\pi}{3^{3/2}} = \frac{4\pi}{3\sqrt3} 3 3/2 4 π = 3 3 4 π で一致します。
Exercise 9.4 難
a > 0 a > 0 a > 0 とするとき ∫ 0 ∞ x sin x x 2 + a 2 d x = π 2 e − a \displaystyle\int_0^{\infty}\frac{x\sin x}{x^2+a^2}\,dx = \frac{\pi}{2}e^{-a} ∫ 0 ∞ x 2 + a 2 x sin x d x = 2 π e − a を示してください。z z 2 + a 2 \dfrac{z}{z^2+a^2} z 2 + a 2 z には Lemma 6.1 (i) が使えないことに注意してください。
Solution f ( z ) = z z 2 + a 2 f(z) = \dfrac{z}{z^2+a^2} f ( z ) = z 2 + a 2 z 、F ( z ) = f ( z ) e i z F(z) = f(z)e^{iz} F ( z ) = f ( z ) e i z とし、半径 R > a R > a R > a の半円輪郭 C R C_R C R を取ります。F F F の上半平面の特異点は単純極 z = i a z = ia z = ia だけで、Proposition 4.2 (3) を g ( z ) = z e i z g(z) = ze^{iz} g ( z ) = z e i z 、h ( z ) = z 2 + a 2 h(z) = z^2+a^2 h ( z ) = z 2 + a 2 、h ′ ( z ) = 2 z h'(z)=2z h ′ ( z ) = 2 z に適用して
Res ( F , i a ) = i a e i ⋅ i a 2 i a = e − a 2 . \operatorname{Res}(F, ia) = \frac{ia\,e^{i\cdot ia}}{2ia} = \frac{e^{-a}}{2}. Res ( F , ia ) = 2 ia ia e i ⋅ ia = 2 e − a . Theorem 5.1 より ∫ − R R x e i x x 2 + a 2 d x + ∫ Γ R F d z = 2 π i ⋅ e − a 2 = i π e − a \displaystyle\int_{-R}^{R}\frac{xe^{ix}}{x^2+a^2}dx + \int_{\Gamma_R}F\,dz = 2\pi i\cdot\frac{e^{-a}}{2} = i\pi e^{-a} ∫ − R R x 2 + a 2 x e i x d x + ∫ Γ R F d z = 2 π i ⋅ 2 e − a = iπ e − a です。
円弧の評価:∣ z ∣ = R > a |z|=R>a ∣ z ∣ = R > a のとき ∣ f ( z ) ∣ ≤ R R 2 − a 2 |f(z)| \le \dfrac{R}{R^2-a^2} ∣ f ( z ) ∣ ≤ R 2 − a 2 R なので M ( R ) → 0 M(R) \to 0 M ( R ) → 0 ですが、R M ( R ) → 1 ≠ 0 R\,M(R) \to 1 \ne 0 R M ( R ) → 1 = 0 なので Lemma 6.1 (i) は使えません。指数因子を持つ (ii)(ジョルダンの補題)なら λ = 1 \lambda = 1 λ = 1 として
∣ ∫ Γ R F ( z ) d z ∣ ≤ π M ( R ) ≤ π R R 2 − a 2 ⟶ 0 \left|\int_{\Gamma_R}F(z)\,dz\right| \le \pi M(R) \le \frac{\pi R}{R^2-a^2}\longrightarrow 0 ∫ Γ R F ( z ) d z ≤ π M ( R ) ≤ R 2 − a 2 π R ⟶ 0 と結論できます。よって lim R → ∞ ∫ − R R x e i x x 2 + a 2 d x = i π e − a \displaystyle\lim_{R\to\infty}\int_{-R}^{R}\frac{xe^{ix}}{x^2+a^2}dx = i\pi e^{-a} R → ∞ lim ∫ − R R x 2 + a 2 x e i x d x = iπ e − a 。虚部を取ると
∫ − ∞ ∞ x sin x x 2 + a 2 d x = π e − a \int_{-\infty}^{\infty}\frac{x\sin x}{x^2+a^2}\,dx = \pi e^{-a} ∫ − ∞ ∞ x 2 + a 2 x sin x d x = π e − a であり、被積分関数は偶関数(奇関数どうしの積)なので ∫ 0 ∞ = π 2 e − a \displaystyle\int_0^{\infty} = \frac{\pi}{2}e^{-a} ∫ 0 ∞ = 2 π e − a です。なおこの積分は絶対収束せず、R → ∞ R \to \infty R → ∞ の対称極限として存在します。
L. V. Ahlfors, Complex Analysis , 3rd ed., McGraw-Hill, 1979 — 第 4 章に留数計算と定積分への応用がまとまっています。真性特異点の値の分布(ピカールの定理)は第 8 章です。
E. M. Stein and R. Shakarchi, Complex Analysis (Princeton Lectures in Analysis II), Princeton University Press, 2003 — 第 3 章「Meromorphic Functions and the Logarithm」が留数公式と本記事で扱った型の実積分を扱っています。
W. Rudin, Real and Complex Analysis , 3rd ed., McGraw-Hill, 1987 — 第 10 章に孤立特異点の分類と、巻き数を使った一般形の留数定理があります。
神保道夫『複素関数入門』岩波書店、2003 — 日本語で読める入門書。留数と実積分への応用を丁寧に扱っています。
A.-L. Cauchy, Exercices de mathématiques , 1826 — 「résidu(留数)」という語と計算法が導入された原典です。
扱う積分 。0 < a < 1 0 < a < 1 0 < a < 1 に対して
∫ 0 ∞ x a − 1 1 + x d x = π sin π a \int_0^{\infty}\frac{x^{a-1}}{1+x}\,dx = \frac{\pi}{\sin \pi a} ∫ 0 ∞ 1 + x x a − 1 d x = sin π a π を示します。x a − 1 x^{a-1} x a − 1 は a a a が整数でないと一価な複素関数に延長できないので、これまでの輪郭は使えません。切断(分枝の切れ目)を正の実軸に置き、その両岸を往復する「鍵穴輪郭」を使います。
分枝を決める 。Log z = log ∣ z ∣ + i arg z \operatorname{Log}z = \log|z| + i\arg z Log z = log ∣ z ∣ + i arg z 、arg z ∈ ( 0 , 2 π ) \arg z \in (0, 2\pi) arg z ∈ ( 0 , 2 π ) と定め、z a − 1 : = e ( a − 1 ) Log z z^{a-1} := e^{(a-1)\operatorname{Log}z} z a − 1 := e ( a − 1 ) Log z とします。これは C \mathbb{C} C から非負実軸を除いた領域で正則です。f ( z ) = z a − 1 1 + z f(z) = \dfrac{z^{a-1}}{1+z} f ( z ) = 1 + z z a − 1 の特異点は z = − 1 = e i π z = -1 = e^{i\pi} z = − 1 = e iπ の単純極だけです。
輪郭 。0 < ε < 1 < R 0 < \varepsilon < 1 < R 0 < ε < 1 < R とし、切断の上岸を ε \varepsilon ε から R R R まで進み、大円 ∣ z ∣ = R |z|=R ∣ z ∣ = R を反時計回りに一周し、下岸を R R R から ε \varepsilon ε まで戻り、小円 ∣ z ∣ = ε |z|=\varepsilon ∣ z ∣ = ε を時計回りに一周して閉じます。厳密には上岸・下岸を実軸から δ \delta δ だけずらした線分に取り、最後に δ → 0 \delta \to 0 δ → 0 とします。この閉曲線の内部にある特異点は z = − 1 z=-1 z = − 1 だけです。
円周上の評価 。大円上では ∣ z a − 1 ∣ = R a − 1 |z^{a-1}| = R^{a-1} ∣ z a − 1 ∣ = R a − 1 、∣ 1 + z ∣ ≥ R − 1 |1+z| \ge R-1 ∣1 + z ∣ ≥ R − 1 なので、ML 不等式より寄与は 2 π R ⋅ R a − 1 R − 1 ∼ 2 π R a − 1 → 0 2\pi R\cdot\dfrac{R^{a-1}}{R-1} \sim 2\pi R^{a-1} \to 0 2 π R ⋅ R − 1 R a − 1 ∼ 2 π R a − 1 → 0 (a < 1 a<1 a < 1 )です。小円上では 2 π ε ⋅ ε a − 1 1 − ε ∼ 2 π ε a → 0 2\pi\varepsilon\cdot\dfrac{\varepsilon^{a-1}}{1-\varepsilon} \sim 2\pi\varepsilon^{a}\to 0 2 π ε ⋅ 1 − ε ε a − 1 ∼ 2 π ε a → 0 (a > 0 a>0 a > 0 )です。0 < a < 1 0 < a < 1 0 < a < 1 という仮定は、この二つの評価を同時に成立させるためにあります。
両岸の差 。上岸では arg z → 0 \arg z \to 0 arg z → 0 なので z a − 1 → x a − 1 z^{a-1} \to x^{a-1} z a − 1 → x a − 1 、下岸では arg z → 2 π \arg z \to 2\pi arg z → 2 π なので z a − 1 → x a − 1 e 2 π i ( a − 1 ) = x a − 1 e 2 π i a z^{a-1}\to x^{a-1}e^{2\pi i(a-1)} = x^{a-1}e^{2\pi ia} z a − 1 → x a − 1 e 2 π i ( a − 1 ) = x a − 1 e 2 π ia です。下岸は逆向きに進むので、二つの線分の寄与の和は I : = ∫ 0 ∞ x a − 1 1 + x d x I := \int_0^{\infty}\frac{x^{a-1}}{1+x}dx I := ∫ 0 ∞ 1 + x x a − 1 d x を使って ( 1 − e 2 π i a ) I \left(1 - e^{2\pi i a}\right)I ( 1 − e 2 π ia ) I と書けます。切断があるおかげで両岸が相殺せず、求めたい積分が残るのがこの輪郭の要点です。
留数と結論 。Res ( f , − 1 ) = ( − 1 ) a − 1 = e ( a − 1 ) i π = − e i π a \operatorname{Res}(f,-1) = (-1)^{a-1} = e^{(a-1)i\pi} = -e^{i\pi a} Res ( f , − 1 ) = ( − 1 ) a − 1 = e ( a − 1 ) iπ = − e iπ a なので、Theorem 5.1 より ( 1 − e 2 π i a ) I = 2 π i ( − e i π a ) \left(1-e^{2\pi ia}\right)I = 2\pi i\left(-e^{i\pi a}\right) ( 1 − e 2 π ia ) I = 2 π i ( − e iπ a ) 、すなわち
I = − 2 π i e i π a 1 − e 2 π i a = − 2 π i e − i π a − e i π a = − 2 π i − 2 i sin π a = π sin π a . I = \frac{-2\pi i\,e^{i\pi a}}{1 - e^{2\pi i a}}
= \frac{-2\pi i}{e^{-i\pi a} - e^{i\pi a}}
= \frac{-2\pi i}{-2i\sin\pi a}
= \frac{\pi}{\sin\pi a}. I = 1 − e 2 π ia − 2 π i e iπ a = e − iπ a − e iπ a − 2 π i = − 2 i sin π a − 2 π i = sin π a π . 検算 。a = 1 / 2 a = 1/2 a = 1/2 とすると I = π / sin ( π / 2 ) = π I = \pi/\sin(\pi/2) = \pi I = π / sin ( π /2 ) = π です。一方 x = t 2 x = t^2 x = t 2 と置換すれば
∫ 0 ∞ d x x ( 1 + x ) = ∫ 0 ∞ 2 t d t t ( 1 + t 2 ) = 2 [ arctan t ] 0 ∞ = π \int_0^{\infty}\frac{dx}{\sqrt{x}(1+x)} = \int_0^{\infty}\frac{2t\,dt}{t(1+t^2)} = 2\bigl[\arctan t\bigr]_0^{\infty} = \pi ∫ 0 ∞ x ( 1 + x ) d x = ∫ 0 ∞ t ( 1 + t 2 ) 2 t d t = 2 [ arctan t ] 0 ∞ = π で一致します。この積分はベータ関数・ガンマ関数の相反公式 Γ ( a ) Γ ( 1 − a ) = π / sin π a \Gamma(a)\Gamma(1-a) = \pi/\sin\pi a Γ ( a ) Γ ( 1 − a ) = π / sin π a の複素解析による証明にもなっています。