コンテンツにスキップ

留数定理と定積分への応用:孤立特異点の分類から実積分・鍵穴輪郭・級数の和まで

対象読者: 本テーマ第 1〜4 章(コーシーの積分公式、テイラー展開、一致の定理)を学んだ学部生。広義積分と三角関数の積分の計算に慣れていることを仮定します。

前提となる記事:正則関数の強力な性質:無限回微分可能性から一致の定理まで

生 Markdown

複素数の世界で関数を調べていると、「ある点でだけ関数の値が決まらない」ことがよくあります。そうした点のまわりで関数がどんな暴れ方をするかを調べ、そこから一つの数を取り出すと、実数の範囲では計算が大変な積分(たとえば ∫−∞∞dxx4+1\int_{-\infty}^{\infty}\frac{dx}{x^4+1})が足し算だけで求まります。この記事では、その数の意味と、それを使った計算の流れを、例を通して順に見ていきます。

  • 正則関数を孤立特異点 aa のまわりでローラン展開したとき、(z−a)−1(z-a)^{-1} の係数 c−1c_{-1} だけが aa を囲む閉曲線に沿う積分に生き残ります。この c−1c_{-1} を留数と呼びます。
  • 孤立特異点は主要部の形で三つに分かれます。除去可能特異点は「近くで有界」、mm 位の極は「(z−a)mf(z)(z-a)^m f(z) が 00 でない有限の極限を持つ」、真性特異点は「どんな小さい近傍の像も複素平面で稠密(複素平面のどの点のそばにも像の点がある、ということ)」と、いずれも極限の言葉に翻訳できます。
  • 留数定理は「単純閉曲線に沿う積分 =2πi×= 2\pi i \times(内部の留数の和)」です。証明は、多重連結領域に対するコーシーの積分定理で積分路を各特異点まわりの小円に分解し、小円上でローラン級数を項別積分するだけです。
  • 実積分への応用は三段構えです。(1) 被積分関数を複素平面に延長して積分路を閉じる、(2) 付け足した部分の寄与が 00 に消えることを不等式で示す、(3) 内部の留数を足す。
  • 型ごとに使う輪郭が決まっています。有理関数の (−∞,∞)(-\infty,\infty) 積分は上半平面の半円、周期 2π2\pi の三角関数の積分は単位円、eiλxe^{i\lambda x} を含む積分はジョルダンの補題(半円弧上の積分が半径を大きくすると 00 に近づくことを保証する不等式)つきの半円、実軸上に単純極があるときは極を小円でくぼませて主値(極の両側を同じ幅で切り落として作る積分値)を取ります。
  • 補遺として、xa−1x^{a-1} を含む積分は正の実軸に切れ目を入れた鍵穴輪郭で 2πi1−e2πia∑Res⁡\frac{2\pi i}{1-e^{2\pi ia}}\sum\operatorname{Res} に帰着し(§9)、整数にわたる級数の和 ∑nf(n)\sum_nf(n) は πcot⁡πz\pi\cot\pi z をかけて正方形の周で積分すると極の留数の和に帰着します(§10。ζ(2)=π2/6\zeta(2)=\pi^2/6 や ∑1n2+a2=πacoth⁡πa\sum\frac1{n^2+a^2}=\frac\pi a\coth\pi a が出ます)。

1. 動機:実積分を複素平面に持ち上げる

Section titled “1. 動機:実積分を複素平面に持ち上げる”

微分積分学で ∫−∞∞dxx2+1\int_{-\infty}^{\infty}\frac{dx}{x^2+1} を計算するのは簡単です。原始関数 arctan⁡x\arctan x が分かっているので、lim⁡R→∞[arctan⁡x]−RR=π\lim_{R\to\infty}[\arctan x]_{-R}^{R} = \pi で終わります。ところが分母を x4+1x^4+1 に変えただけで景色が一変します。実数の範囲では x4+1=(x2+2x+1)(x2−2x+1)x^4+1 = (x^2+\sqrt{2}x+1)(x^2-\sqrt{2}x+1) と因数分解し、部分分数に分け、四つの項をそれぞれ arctan⁡\arctan と log⁡\log で積分し、最後に極限を取る、という長い計算になります。∫0∞sin⁡xx dx\int_0^{\infty}\frac{\sin x}{x}\,dx に至っては、初等関数で書ける原始関数がそもそも存在しないので、パラメータ微分などの技巧なしには実数の世界で手が出ません。

この章で扱う留数定理は、この種の積分を次の手続きに置き換えます。被積分関数を複素平面に延長し、積分路を閉曲線になるように閉じ、閉曲線の内部にある特異点で「留数」という一つの数を計算して足す。それだけです。∫−∞∞dxx4+1\int_{-\infty}^{\infty}\frac{dx}{x^4+1} は、後で見るように二つの留数の足し算だけで π/2\pi/\sqrt{2} と分かります。

なぜそんなことができるのでしょうか。出発点は コーシーの積分定理(定理 4.3[コーシーの積分定理と積分公式])です。単連結領域で正則な関数は閉曲線に沿う積分が 00 になります。つまり、閉曲線に沿う積分が 00 でないとしたら、その値は内部にある特異点だけが原因です。では特異点はどれだけ寄与するのでしょうか。中心 aa の小さい円周 ∣z−a∣=r|z-a| = r 上で (z−a)n(z-a)^n を積分してみます。n≠−1n \ne -1 のとき (z−a)n(z-a)^n は C∖{a}\mathbb{C}\setminus\{a\} 全体で一価(各点で値が一つに決まる)な原始関数 (z−a)n+1n+1\frac{(z-a)^{n+1}}{n+1} を持つので、閉曲線に沿う積分は 00 です。ところが n=−1n = -1 のときだけ、原始関数の候補は log⁡(z−a)\log(z-a) という多価関数しかなく、円周を一周すると偏角が 2π2\pi 増えて、積分は 2πi2\pi i を返します。

したがって、ローラン展開 f(z)=∑n=−∞∞cn(z−a)nf(z) = \sum_{n=-\infty}^{\infty}c_n(z-a)^n を円周上で項別に積分すると、無限個の項のうち c−1(z−a)−1c_{-1}(z-a)^{-1} ただ一つだけが生き残ります。この「取り残される数」がフランス語で résidu(残りもの)、日本語で留数です。コーシーは 1826 年の『Exercices de mathématiques』でこの語を導入し、実積分の計算に系統的に用いました。この章は、その計算術を最初から最後まで具体的に組み立てます。

2. 準備:記号とすでに証明した道具

Section titled “2. 準備:記号とすでに証明した道具”

記号:領域とは連結開集合のことです。D(a,r)={z∈C:∣z−a∣<r}D(a,r) = \{z \in \mathbb{C} : |z-a| < r\}、D′(a,r)={z:0<∣z−a∣<r}D'(a,r) = \{z : 0 < |z-a| < r\}(中心を抜いた円板、除外近傍)と書きます。曲線はすべて区分的に C1C^1 級とし、∫γf(z) dz=∫αβf(γ(t))γ′(t) dt\int_\gamma f(z)\,dz = \int_\alpha^\beta f(\gamma(t))\gamma'(t)\,dt で定義します。円周 ∣z−a∣=r|z-a| = r には断らない限り正の向き(反時計回り)を与え、閉曲線に沿う積分は ∮\oint と書きます。

長さと最大値による評価(ML 不等式):曲線 γ\gamma の長さを L(γ)L(\gamma) とすると

∣∫γf(z) dz∣≤∫γ∣f(z)∣ ∣dz∣≤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)|

が成り立ちます(補題 2.4[コーシーの積分定理と積分公式])。この章の評価はほとんどすべてこの不等式に帰着します。

ローラン展開:正則関数の強力な性質 の 定理 6.1[正則関数の強力な性質] で示したとおり、ff が円環 r<∣z−a∣<Rr < |z-a| < R で正則ならば、そこで

f(z)=∑n=−∞∞cn(z−a)n,cn=12π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)

と一意的に展開でき、この級数は円環の任意のコンパクト部分集合上で一様収束します。係数 cnc_n は ρ\rho の取り方によりません。

多重連結領域に対するコーシーの積分定理:CC を正の向きの単純閉曲線、γ1,…,γn\gamma_1,\dots,\gamma_n を CC の内部にあり互いに他の外部にある正の向きの単純閉曲線とします。ff が「CC の内側かつすべての γj\gamma_j の外側」という閉領域を含むある開集合で正則ならば

∮Cf(z) dz=∑j=1n∮γjf(z) dz\oint_C f(z)\,dz = \sum_{j=1}^{n}\oint_{\gamma_j} f(z)\,dz

が成り立ちます。証明はコーシーの積分定理と積分公式の 定理 5.1[コーシーの積分定理と積分公式] と同じ要領で、CC と各 γj\gamma_j を細い「橋」(切断線)でつないで単連結な部分に分け、各橋を往復する積分が符号違いで相殺することを使います。留数定理はこの定理の言い換えにすぎません。

定義 3.1(孤立特異点)↗

a∈Ca \in \mathbb{C} とする。ある R>0R > 0 が存在して ff が除外近傍 D′(a,R)={z:0<∣z−a∣<R}D'(a,R) = \{z : 0 < |z-a| < R\} で正則であるとき、aa を ff の孤立特異点と呼ぶ。ff は aa 自身では定義されていなくてもよいし、定義されていても正則でなくてよい。

孤立特異点のまわりでは ff は円環 0<∣z−a∣<R0 < |z-a| < R で正則ですから、§2 のローラン展開が使えます。展開の負べきの部分が特異点の「悪さ」を測る量です。

定義 3.2(主要部と特異点の三分類)↗

aa を ff の孤立特異点とし、D′(a,R)D'(a,R) でのローラン展開を f(z)=∑n=−∞∞cn(z−a)nf(z) = \sum_{n=-\infty}^{\infty}c_n(z-a)^n とする。負べきの部分

Pa(z)=∑n≤−1cn(z−a)nP_a(z) = \sum_{n \le -1} c_n (z-a)^n

を aa における ff の主要部と呼ぶ。主要部の形により、aa を次のように分類する。

  • cn=0c_n = 0(すべての n≤−1n \le -1)、すなわち主要部が 00 のとき、aa を除去可能特異点という。
  • ある正の整数 mm について c−m≠0c_{-m} \ne 0 かつ cn=0c_n = 0(すべての n<−mn < -m)のとき、aa を mm 位の極という。m=1m = 1 のときは単純極ともいう。
  • cn≠0c_n \ne 0 となる n≤−1n \le -1 が無限個あるとき、aa を真性特異点という。

ローラン展開は一意なので、この分類は展開の取り方によらず定まる。

孤立特異点の三分類。主要部の項数を数えるだけで型が決まり、それぞれに極限による判定条件が対応します。

分類の定義はローラン係数を全部調べよと言っているように見えますが、実際にはもっと安上がりな判定法があります。まず除去可能特異点です。

定理 3.3(リーマンの除去可能特異点定理)✓

aa を ff の孤立特異点とし、ff は D′(a,R)D'(a,R) で正則とする。次の 4 条件は同値である。

  • (i) aa は ff の除去可能特異点である。
  • (ii) f(a)f(a) の値を適当に定めると、ff は D(a,R)D(a,R) 全体で正則になる。
  • (iii) ある 0<r≤R0 < r \le R と M>0M > 0 があって、0<∣z−a∣<r0 < |z-a| < r で ∣f(z)∣≤M|f(z)| \le M が成り立つ。
  • (iv) lim⁡z→a(z−a)f(z)=0\lim_{z \to a}(z-a)f(z) = 0。
証明(定理 3.3)

(i) ⇒\Rightarrow (ii):主要部が 00 なので D′(a,R)D'(a,R) で f(z)=∑n≥0cn(z−a)nf(z) = \sum_{n \ge 0}c_n(z-a)^n です。右辺のべき級数は D′(a,R)D'(a,R) で収束するので収束半径は RR 以上であり、D(a,R)D(a,R) 全体で正則関数を定めます。f(a):=c0f(a) := c_0 と定めれば、ff はこのべき級数と D(a,R)D(a,R) 全体で一致し、正則になります。

(ii) ⇒\Rightarrow (iii):正則関数は連続なので、r=R/2r = R/2 とすればコンパクト集合 D(a,r)‾\overline{D(a,r)} 上で ∣f∣|f| は最大値 MM を取ります。

(iii) ⇒\Rightarrow (iv):∣(z−a)f(z)∣≤M∣z−a∣→0|(z-a)f(z)| \le M|z-a| \to 0(z→az \to a)です。

(iv) ⇒\Rightarrow (i):0<ρ<R0 < \rho < R に対し η(ρ):=max⁡∣z−a∣=ρ∣(z−a)f(z)∣\eta(\rho) := \max_{|z-a| = \rho}|(z-a)f(z)| とおくと、仮定より η(ρ)→0\eta(\rho) \to 0(ρ→0\rho \to 0)です。§2 のローラン係数の式と ML 不等式から、n≤−1n \le -1 に対して

∣cn∣≤12π⋅2πρ⋅max⁡∣z−a∣=ρ∣f(z)∣ρ n+1=ρ−nmax⁡∣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)

が成り立ちます。ここで n≤−1n \le -1 より −n−1≥0-n-1 \ge 0 なので、ρ<1\rho < 1 のとき ρ−n−1≤1\rho^{-n-1} \le 1 です。よって ∣cn∣≤η(ρ)→0|c_n| \le \eta(\rho) \to 0 となり、cnc_n は ρ\rho によらない定数ですから cn=0c_n = 0 です。これがすべての n≤−1n \le -1 で成り立つので主要部は 00 です。

条件 (iii) が重要です。「近くで有界」という、ローラン係数を一つも計算しなくてよい条件だけで、特異点が消せると分かります。次は極です。

命題 3.4(極の特徴づけ)✓

aa を ff の孤立特異点、mm を正の整数とする。

(1) 次の 3 条件は同値である。

  • (i) aa は ff の mm 位の極である。
  • (ii) 極限 lim⁡z→a(z−a)mf(z)\lim_{z \to a}(z-a)^m f(z) が存在し、その値は 00 でない。
  • (iii) ある r>0r > 0 と D(a,r)D(a,r) で正則な関数 gg で g(a)≠0g(a) \ne 0 を満たすものがあり、0<∣z−a∣<r0 < |z-a| < r で f(z)=g(z)(z−a)mf(z) = \dfrac{g(z)}{(z-a)^m} と書ける。

(2) aa が(何らかの位数の)極であるための必要十分条件は lim⁡z→a∣f(z)∣=∞\lim_{z \to a}|f(z)| = \infty である。

証明(命題 3.4)

(1) (i) ⇒\Rightarrow (iii):mm 位の極なら D′(a,R)D'(a,R) で f(z)=∑n≥−mcn(z−a)nf(z) = \sum_{n \ge -m}c_n(z-a)^n かつ c−m≠0c_{-m} \ne 0 です。両辺に (z−a)m(z-a)^m を掛けると

(z−a)mf(z)=∑n≥−mcn(z−a)n+m=∑k≥0ck−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

となり、右辺は D′(a,R)D'(a,R) で収束するべき級数ですから収束半径は RR 以上で、D(a,R)D(a,R) で正則関数 gg を定めます。g(a)=c−m≠0g(a) = c_{-m} \ne 0 であり、0<∣z−a∣<R0 < |z-a| < R で f(z)=g(z)/(z−a)mf(z) = g(z)/(z-a)^m です。

(iii) ⇒\Rightarrow (ii):(z−a)mf(z)=g(z)→g(a)≠0(z-a)^m f(z) = g(z) \to g(a) \ne 0(z→az \to a、gg の連続性)です。

(ii) ⇒\Rightarrow (i):h(z):=(z−a)mf(z)h(z) := (z-a)^m f(z) とおくと、hh は D′(a,R)D'(a,R) で正則で z→az \to a のとき有限の極限 L≠0L \ne 0 を持つので、aa の近くで有界です。定理 3.3 の (iii) ⇒\Rightarrow (i),(ii) より aa は hh の除去可能特異点で、h(a)=Lh(a) = L と定めれば hh は D(a,R)D(a,R) で正則になります。hh をテイラー展開して h(z)=∑k≥0bk(z−a)kh(z) = \sum_{k \ge 0}b_k(z-a)^k(b0=L≠0b_0 = L \ne 0)とすると、0<∣z−a∣<R0 < |z-a| < R で

f(z)=h(z)(z−a)m=∑k≥0bk(z−a)k−mf(z) = \frac{h(z)}{(z-a)^m} = \sum_{k \ge 0} b_k (z-a)^{k-m}

となります。これは ff のローラン展開そのもの(一意性より)で、最低次の項は b0(z−a)−mb_0(z-a)^{-m}、b0≠0b_0 \ne 0 ですから aa は mm 位の極です。

(2) 必要性:aa が mm 位の極なら (iii) より f=g/(z−a)mf = g/(z-a)^m、g(a)≠0g(a) \ne 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)∣→∞|f(z)| \to \infty とします。ある r>0r > 0 があって 0<∣z−a∣<r0 < |z-a| < r で ∣f(z)∣≥1|f(z)| \ge 1、特に f(z)≠0f(z) \ne 0 です。そこで u:=1/fu := 1/f とおくと uu は D′(a,r)D'(a,r) で正則、∣u∣≤1|u| \le 1 と有界なので、定理 3.3 より aa は uu の除去可能特異点で、uu は D(a,r)D(a,r) で正則に延長されます。∣f∣→∞|f| \to \infty より u(a)=0u(a) = 0 です。uu は D′(a,r)D'(a,r) で 00 にならないので u≢0u \not\equiv 0 であり、補題 7.2[正則関数の強力な性質] で見たとおり、恒等的に 00 でない正則関数の零点は有限の位数を持ちます。その位数を mm とすれば u(z)=(z−a)mv(z)u(z) = (z-a)^m v(z)、vv は D(a,r′)D(a,r') で正則、v(a)≠0v(a) \ne 0 と書けます。よって 0<∣z−a∣<r′0 < |z-a| < r' で f(z)=1/v(z)(z−a)mf(z) = \dfrac{1/v(z)}{(z-a)^m} となり、1/v1/v は aa の近くで正則(v(a)≠0v(a)\ne 0 と連続性より近傍で v≠0v \ne 0)かつ aa で 00 でないので、(1)(iii) から aa は mm 位の極です。

残るのが真性特異点です。定理 3.3 と 命題 3.4 により、aa が真性特異点であることは「z→az \to a で f(z)f(z) が有限の極限も持たず、∣f(z)∣→∞|f(z)| \to \infty でもない」ことと同値です。極限が存在しないだけでなく、値の散らばり方は極端です。

次の定理は、真性特異点ではまわりの値が極端に散らばることを述べています。稠密とは、複素平面のどの点のそばにも像の点が見つかる、ということです。証明は背理法で、もし像があるひとつの点 ww のそばを避けているなら、1/(f−w)1/(f-w) が特異点の近くで有界になり、いま示した除去可能特異点の定理で正則に延長できます。すると ff 自身が除去可能特異点か極になってしまい、真性特異点であることに矛盾する、という流れです。

定理 3.5(カゾラーティ・ワイエルシュトラスの定理)✓

aa を ff の真性特異点とし、ff は D′(a,R)D'(a,R) で正則とする。このとき任意の 0<r≤R0 < r \le R に対して像 f(D′(a,r))f(D'(a,r)) は C\mathbb{C} で稠密である。すなわち、任意の w∈Cw \in \mathbb{C} と任意の ε>0\varepsilon > 0 に対し、0<∣z−a∣<r0 < |z-a| < r かつ ∣f(z)−w∣<ε|f(z) - w| < \varepsilon を満たす zz が存在する。

証明(定理 3.5)

稠密でないと仮定します。すると、ある w∈Cw \in \mathbb{C} と δ>0\delta > 0 があって、0<∣z−a∣<r0 < |z-a| < r のすべての zz で ∣f(z)−w∣≥δ|f(z) - w| \ge \delta となります。そこで

g(z):=1f(z)−w(0<∣z−a∣<r)g(z) := \frac{1}{f(z) - w} \qquad (0 < |z-a| < r)

とおくと、分母は 00 にならないので gg は D′(a,r)D'(a,r) で正則で、∣g∣≤1/δ|g| \le 1/\delta と有界です。定理 3.3 より aa は gg の除去可能特異点で、gg は D(a,r)D(a,r) で正則に延長されます。

g(a)≠0g(a) \ne 0 の場合:f(z)=w+1/g(z)f(z) = w + 1/g(z) は aa の近くで有界なので、再び 定理 3.3 により aa は ff の除去可能特異点です。

g(a)=0g(a) = 0 の場合:gg は D′(a,r)D'(a,r) で 00 にならないので g≢0g \not\equiv 0 であり、aa は gg の有限位数 mm の零点です。命題 3.4 の (1)(iii) の形に書けば f(z)−wf(z) - w は aa で mm 位の極を持ち、定数 ww を足しても主要部は変わらないので ff も aa で mm 位の極を持ちます。

どちらの場合も aa が真性特異点であることに反します。

例 3.6(特異点の分類の実例)

(a) f(z)=sin⁡zzf(z) = \dfrac{\sin z}{z} の z=0z = 0。sin⁡z=z−z33!+z55!−⋯\sin z = z - \frac{z^3}{3!} + \frac{z^5}{5!} - \cdots を zz で割ると f(z)=1−z26+z4120−⋯f(z) = 1 - \frac{z^2}{6} + \frac{z^4}{120} - \cdots となり、負べきの項がありません。よって 00 は除去可能特異点で、f(0):=1f(0) := 1 と定めれば整関数になります。

(b) f(z)=1z2(z−1)f(z) = \dfrac{1}{z^2(z-1)}。z=0z = 0 では z2f(z)=1z−1→−1≠0z^2 f(z) = \frac{1}{z-1} \to -1 \ne 0 なので 命題 3.4 より 22 位の極、z=1z = 1 では (z−1)f(z)=1z2→1≠0(z-1)f(z) = \frac{1}{z^2} \to 1 \ne 0 なので単純極です。

(c) f(z)=1sin⁡zf(z) = \dfrac{1}{\sin z} の z=kπz = k\pi(k∈Zk \in \mathbb{Z})。sin⁡(kπ)=0\sin(k\pi) = 0 かつ ddzsin⁡z∣z=kπ=cos⁡(kπ)=(−1)k≠0\frac{d}{dz}\sin z\big|_{z=k\pi} = \cos(k\pi) = (-1)^k \ne 0 なので kπk\pi は sin⁡z\sin z の 11 位の零点、したがって

(z−kπ)f(z)=z−kπsin⁡z−sin⁡(kπ)⟶1cos⁡(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

(微分係数の定義そのもの)となり、kπk\pi は単純極です。

(d) f(z)=e1/zf(z) = e^{1/z} の z=0z = 0。ew=∑n≥0wn/n!e^{w} = \sum_{n\ge0}w^n/n! に w=1/zw = 1/z を代入して

e1/z=∑n≥01n!z−n=1+1z+12!z2+⋯(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)

が 0<∣z∣0 < |z| で収束します。負べきの項が無限個あるので 00 は真性特異点です。実際、値の散らばり方を直接確かめられます。w≠0w \ne 0 を任意に取り w=∣w∣eiθw = |w|e^{i\theta} と書くと、zk=1log⁡∣w∣+i(θ+2πk)z_k = \dfrac{1}{\log|w| + i(\theta + 2\pi k)}(k∈Nk \in \mathbb{N})は k→∞k \to \infty で 00 に収束し、e1/zk=we^{1/z_k} = w です。つまり 00 のどんな近傍でも、ff は 00 以外のすべての複素数を無限回取ります。定理 3.5 の「稠密」よりずっと強い状況です。

注意 3.7

例 3.6 (d) の現象は一般に成り立ちます。ピカールの大定理によれば、真性特異点のどんな除外近傍でも、ff は高々 1 個の例外値を除くすべての複素数値を無限回取ります。証明はこの章の範囲を超えるので、Ahlfors『Complex Analysis』の第 8 章(Picard の定理を扱う章)を参照してください。

注意 3.8

すべての特異点が孤立しているわけではありません。f(z)=1sin⁡(1/z)f(z) = \dfrac{1}{\sin(1/z)} は z=1/(kπ)z = 1/(k\pi)(kk は 00 でない整数)で極を持ちますが、これらは k→∞k \to \infty で 00 に集積します。したがって 00 のどんな除外近傍にも ff が正則でない点が含まれ、00 は 定義 3.1 の意味の孤立特異点ではありません。ローラン展開も留数も定義できないので、この章の道具は 00 には使えません。

定義 4.1(留数)↗

aa を ff の孤立特異点とし、D′(a,R)D'(a,R) でのローラン展開を f(z)=∑n=−∞∞cn(z−a)nf(z) = \sum_{n=-\infty}^{\infty}c_n(z-a)^n とする。係数 c−1c_{-1} を ff の aa における留数と呼び、

Res⁡(f,a):=c−1\operatorname{Res}(f, a) := c_{-1}

と書く。ローラン係数の積分表示(§2)で n=−1n = -1 とすれば、0<ρ<R0 < \rho < R に対して

Res⁡(f,a)=12πi∮∣z−a∣=ρf(z) dz\operatorname{Res}(f,a) = \frac{1}{2\pi i}\oint_{|z-a| = \rho} f(z)\,dz

でもある。

定義そのものは「ローラン展開せよ」と言っていますが、極の場合には展開せずに留数を取り出す公式があります。実用上はこの公式だけで九割方の計算が片付きます。

命題 4.2(極における留数の計算法)✓

aa を ff の孤立特異点とする。

(1) aa が ff の高々 mm 位の極である(すなわちローラン展開で n<−mn < -m の係数がすべて 00 である)とき、φ(z):=(z−a)mf(z)\varphi(z) := (z-a)^m f(z) は aa で正則に延長され、

Res⁡(f,a)=1(m−1)!lim⁡z→ad m−1dz m−1[(z−a)mf(z)].\operatorname{Res}(f,a) = \frac{1}{(m-1)!}\lim_{z \to a}\frac{d^{\,m-1}}{dz^{\,m-1}}\Bigl[(z-a)^m f(z)\Bigr].

(2) 特に aa が単純極(m=1m=1)のとき Res⁡(f,a)=lim⁡z→a(z−a)f(z)\operatorname{Res}(f,a) = \lim_{z \to a}(z-a)f(z)。

(3) f=g/hf = g/h とし、g,hg,h は aa の近傍で正則、g(a)≠0g(a) \ne 0、h(a)=0h(a) = 0、h′(a)≠0h'(a) \ne 0 とする。このとき aa は ff の単純極で

Res⁡(f,a)=g(a)h′(a).\operatorname{Res}(f,a) = \frac{g(a)}{h'(a)}.
証明(命題 4.2)

(1):仮定より D′(a,R)D'(a,R) で f(z)=∑n≥−mcn(z−a)nf(z) = \sum_{n \ge -m}c_n(z-a)^n なので、

φ(z)=(z−a)mf(z)=∑k≥0ck−m(z−a)k\varphi(z) = (z-a)^m f(z) = \sum_{k \ge 0}c_{k-m}(z-a)^k

であり、右辺は D′(a,R)D'(a,R) で収束するべき級数、したがって収束半径は RR 以上で D(a,R)D(a,R) で正則関数を定めます(命題 3.4 の証明と同じ議論です)。この級数の (z−a)m−1(z-a)^{m-1} の係数は k=m−1k = m-1 の項、すなわち c−1c_{-1} です。一方、正則関数のテイラー係数の公式(定理 5.2[正則関数の強力な性質])から (z−a)m−1(z-a)^{m-1} の係数は φ(m−1)(a)/(m−1)!\varphi^{(m-1)}(a)/(m-1)! です。両者を等しいとおけば結論を得ます。φ\varphi は aa で正則なので φ(m−1)\varphi^{(m-1)} は aa で連続で、φ(m−1)(a)=lim⁡z→aφ(m−1)(z)\varphi^{(m-1)}(a) = \lim_{z\to a}\varphi^{(m-1)}(z) と極限で書けます。

(2):(1) で m=1m = 1 とすると 0!=10! = 1、00 階微分は恒等写像なので、Res⁡(f,a)=lim⁡z→a(z−a)f(z)\operatorname{Res}(f,a) = \lim_{z\to a}(z-a)f(z) です。

(3):h(a)=0h(a) = 0 かつ h′(a)≠0h'(a) \ne 0 より aa は hh の 11 位の零点で、aa の近傍で h(z)=(z−a)h1(z)h(z) = (z-a)h_1(z)、h1h_1 は正則で h1(a)=h′(a)≠0h_1(a) = h'(a) \ne 0 と書けます。よって f(z)=g(z)/h1(z)z−af(z) = \dfrac{g(z)/h_1(z)}{z-a} であり、g/h1g/h_1 は aa で正則、aa での値は g(a)/h′(a)≠0g(a)/h'(a) \ne 0 ですから、命題 3.4 (1)(iii) より aa は単純極です。(2) を使って

Res⁡(f,a)=lim⁡z→a(z−a)f(z)=lim⁡z→ag(z)h1(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)}

を得ます。

例 4.3(留数の計算例)

(a) f(z)=z(z−1)(z+2)2f(z) = \dfrac{z}{(z-1)(z+2)^2}。z=1z=1 は単純極なので 命題 4.2 (2) より

Res⁡(f,1)=lim⁡z→1z(z+2)2=19.\operatorname{Res}(f,1) = \lim_{z\to1}\frac{z}{(z+2)^2} = \frac{1}{9}.

z=−2z=-2 は 22 位の極なので (1) で m=2m=2 とし、

Res⁡(f,−2)=lim⁡z→−2ddz[zz−1]=lim⁡z→−2(z−1)−z(z−1)2=−19.\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}.

二つの留数の和が 00 になっているのは偶然ではありません。分母の次数が分子より 22 以上大きい有理関数では、後の 補題 6.1 (i) と 定理 5.1 を大きな円周に適用すると留数の総和が 00 と分かります。検算に使えます。

(b) f(z)=ezz3f(z) = \dfrac{e^z}{z^3}。m=3m=3 として

Res⁡(f,0)=12!lim⁡z→0d2dz2[ez]=12.\operatorname{Res}(f,0) = \frac{1}{2!}\lim_{z\to0}\frac{d^2}{dz^2}\bigl[e^z\bigr] = \frac{1}{2}.

ローラン展開 ez/z3=z−3+z−2+12z−1+16+⋯e^z/z^3 = z^{-3} + z^{-2} + \frac{1}{2}z^{-1} + \frac{1}{6} + \cdots の z−1z^{-1} の係数と一致します。

(c) f(z)=cot⁡z=cos⁡zsin⁡zf(z) = \cot z = \dfrac{\cos z}{\sin z} の z=kπz = k\pi。g=cos⁡g = \cos、h=sin⁡h = \sin とすると g(kπ)=(−1)k≠0g(k\pi) = (-1)^k \ne 0、h(kπ)=0h(k\pi) = 0、h′(kπ)=cos⁡(kπ)=(−1)k≠0h'(k\pi) = \cos(k\pi) = (-1)^k \ne 0 なので 命題 4.2 (3) より

Res⁡(cot⁡,kπ)=cos⁡(kπ)cos⁡(kπ)=1.\operatorname{Res}(\cot, k\pi) = \frac{\cos(k\pi)}{\cos(k\pi)} = 1.

すべての極で留数が 11 です。

(d) f(z)=z2e1/zf(z) = z^2 e^{1/z} の z=0z = 0。真性特異点なので (1)〜(3) は使えません。定義に戻ってローラン展開します。

z2e1/z=z2∑n≥01n!z−n=∑n≥01n!z2−nz^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−1z^{-1} が現れるのは 2−n=−12-n = -1、すなわち n=3n = 3 のときですから Res⁡(f,0)=13!=16\operatorname{Res}(f,0) = \frac{1}{3!} = \frac{1}{6} です。留数は極の専売特許ではありません。

定理 5.1(留数定理)✓

D⊂CD \subset \mathbb{C} を開集合とし、CC を区分的に C1C^1 級の単純閉曲線で、正の向き(内部を左に見る向き)を与えたものとする。CC の内部(ジョルダンの曲線定理でいう有界成分)を int⁡C\operatorname{int}C と書き、CC と int⁡C\operatorname{int}C はともに DD に含まれるとする(DD が単連結領域なら、C⊂DC\subset D だけからこれが従う)。相異なる有限個の点 a1,…,an∈int⁡Ca_1,\dots,a_n \in \operatorname{int}C を取り、ff は D∖{a1,…,an}D \setminus \{a_1,\dots,a_n\} で正則とする。このとき

∮Cf(z) dz=2πi∑j=1nRes⁡(f,aj).\oint_C f(z)\,dz = 2\pi i \sum_{j=1}^{n}\operatorname{Res}(f,a_j).

証明の方針は、§1 の考え方をそのまま厳密にすることです。まず大きな曲線 CC に沿う積分を、各特異点のまわりの小さな円周に沿う積分の和に置き換えます。小円の上ではローラン展開を項ごとに積分でき、(z−a)−1(z-a)^{-1} の項だけが残ります。途中の領域の連結性や境界の向きの確認は、この置き換えが正しいことを確かめるための作業です。

証明(定理 5.1)

仮定から int⁡C⊂D\operatorname{int}C \subset D です(DD が単連結領域の場合に C⊂DC\subset D からこれが従うことは、定理 4.3[コーシーの積分定理と積分公式] の証明の後半で示しました)。したがって各 aja_j は ff の孤立特異点で、留数が定義されています。

小円を取る。ε>0\varepsilon > 0 を

ε<12min⁡j≠k∣aj−ak∣,ε<min⁡jdist⁡(aj,∂(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))

を満たすように取ります(aja_j は有限個で、int⁡C\operatorname{int}C は開集合なので、このような ε\varepsilon が存在します)。γj\gamma_j を正の向きの円周 ∣z−aj∣=ε|z - a_j| = \varepsilon とすると、閉円板 D(aj,ε)‾\overline{D(a_j,\varepsilon)} は互いに交わらず、すべて int⁡C\operatorname{int}C に含まれます。

積分路を分解する。Ω:=int⁡C∖⋃jD(aj,ε)‾\Omega := \operatorname{int}C \setminus \bigcup_j \overline{D(a_j,\varepsilon)} とおきます。Ω\Omega は有界領域です(ε\varepsilon より少し大きい半径の円環も int⁡C\operatorname{int}C に含まれるので、int⁡C\operatorname{int}C の中の道が閉円板に入る部分はその円環の中の円弧で迂回でき、連結性が保たれます)。その境界は CC と γ1,…,γn\gamma_1,\dots,\gamma_n の n+1n+1 本の単純閉曲線からなり、Ω\Omega を左に見る向きは CC については正の向き、各 γj\gamma_j については逆向きです。ff は Ω‾\overline{\Omega} を含む開集合 D∖{a1,…,an}D \setminus \{a_1,\dots,a_n\} で正則なので、系 3.3[正則関数の強力な性質] により f=u+ivf = u+iv の u,vu,v は C1C^1 級で、コーシー・リーマンの関係式(定理 4.1[正則関数とコーシー・リーマンの関係式])を満たします。補題 4.2[コーシーの積分定理と積分公式] を Ω\Omega に適用すると被積分関数が恒等的に 00 なので、境界に沿う積分は 00 で、命題 2.3[コーシーの積分定理と積分公式] (2)(3) で向きと和を整理すると(§2 の多重連結領域に対するコーシーの積分定理)

∮Cf(z) dz=∑j=1n∮γjf(z) dz.\oint_C f(z)\,dz = \sum_{j=1}^{n}\oint_{\gamma_j}f(z)\,dz .

小円上で項別積分する。jj を固定します。ff は D′(aj,Rj)D'(a_j, R_j) でローラン展開 f(z)=∑k=−∞∞ck(j)(z−aj)kf(z) = \sum_{k=-\infty}^{\infty}c^{(j)}_k(z-a_j)^k を持ちます。必要なら ε\varepsilon をさらに小さくして ε<Rj\varepsilon < R_j としておきます。円周 γj\gamma_j は円環 D′(aj,Rj)D'(a_j,R_j) のコンパクト部分集合なので、この級数は γj\gamma_j 上で一様収束します。一様収束する連続関数列は長さ有限の曲線上で項別積分できる(補題 2.2[正則関数の強力な性質])ので

∮γjf(z) dz=∑k=−∞∞ck(j)∮γj(z−aj)k dz.\oint_{\gamma_j} f(z)\,dz = \sum_{k=-\infty}^{\infty}c^{(j)}_k \oint_{\gamma_j}(z-a_j)^k\,dz .

右辺の積分を計算します。z=aj+εeiθz = a_j + \varepsilon e^{i\theta}(0≤θ≤2π0 \le \theta \le 2\pi)と置くと dz=iεeiθdθdz = i\varepsilon e^{i\theta}d\theta なので

∮γj(z−aj)kdz=∫02πεkeikθ⋅iεeiθ dθ=i ε k+1∫02πei(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 .

k≠−1k \ne -1 のとき ∫02πei(k+1)θdθ=[ei(k+1)θi(k+1)]02π=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(e2πi(k+1)=1e^{2\pi i(k+1)} = 1)です。k=−1k = -1 のときは被積分関数が 11 なので 2π2\pi となり、積分は i⋅ε0⋅2π=2πii \cdot \varepsilon^0 \cdot 2\pi = 2\pi i です。よって

∮γjf(z) dz=2πi c−1(j)=2πiRes⁡(f,aj)\oint_{\gamma_j}f(z)\,dz = 2\pi i\, c^{(j)}_{-1} = 2\pi i \operatorname{Res}(f,a_j)

となり(最後の等号は 定義 4.1)、これを jj について足せば結論を得ます。

注意 5.2

向きは結果の符号を決めるので常に確認してください。CC を逆向きに回れば右辺は −2πi∑Res⁡-2\pi i \sum \operatorname{Res} になります。また、特異点が CC の上にある場合はこの定理は適用できません。実軸上に極がある実積分を扱う §8 では、輪郭を変形して極を避ける必要が出てきます。

ここから先は応用です。実積分を計算するときの手順はいつも同じで、積分の型ごとに使う輪郭と、余分な部分を消す道具が決まっています。

実積分の型使う輪郭余分な部分を消す道具
∫−∞∞P(x)/Q(x) dx\int_{-\infty}^{\infty}P(x)/Q(x)\,dx(deg⁡Q≥deg⁡P+2\deg Q \ge \deg P + 2)上半平面の半円補題 6.1 (i)
∫02πR(cos⁡θ,sin⁡θ) dθ\int_0^{2\pi}R(\cos\theta,\sin\theta)\,d\theta単位円 ∣z∣=1\lvert z\rvert = 1不要(輪郭が最初から閉じている)
∫−∞∞f(x)eiλxdx\int_{-\infty}^{\infty}f(x)e^{i\lambda x}dx(λ>0\lambda > 0)上半平面の半円補題 6.1 (ii)(ジョルダンの補題)
実軸上に単純極があるとき極を小円でくぼませた輪郭補題 8.3 と主値
∫0∞xa−1g(x) dx\int_0^{\infty}x^{a-1}g(x)\,dx(aa は非整数)鍵穴輪郭(§9。付録に最初の例)分枝の指定と大小円の評価(定理 9.1)
∑n∈Zf(n)\sum_{n\in\mathbb{Z}}f(n)(ff は有理関数)πcot⁡πz\pi\cot\pi z をかけた正方形の周(§10)補題 10.1 の一様評価(定理 10.2)
a₁a₂−RROΓR実軸上の線分
有理関数の広義積分に使う輪郭。実軸上の線分と上半平面の半円弧をつなぎ、上半平面にある極(×)だけを内部に取り込みます。半径 R を大きくすると弧の寄与が消えます。

輪郭を閉じるために付け足した半円弧の寄与は、半径を大きくすると消えてほしいわけです。それを保証するのが次の補題です。(i) と (ii) の仮定の違いに注意してください。(ii) のほうが ff に対する要求がゆるく、そのぶん eiλze^{i\lambda z} という因子が必要です。

補題 6.1(大円弧の寄与が消える条件)✓

R0>0R_0 > 0 とし、ΓR\Gamma_R を上半平面の半円弧 z=Reiθz = Re^{i\theta}(0≤θ≤π0 \le \theta \le \pi)とする。ff は {z:∣z∣≥R0, Im⁡z≥0}\{z : |z| \ge R_0,\ \operatorname{Im}z \ge 0\} で連続とし、M(R):=max⁡0≤θ≤π∣f(Reiθ)∣M(R) := \max_{0 \le \theta \le \pi}|f(Re^{i\theta})| とおく。

(i) lim⁡R→∞R M(R)=0\lim_{R\to\infty} R\,M(R) = 0 ならば lim⁡R→∞∫ΓRf(z) dz=0\displaystyle\lim_{R\to\infty}\int_{\Gamma_R}f(z)\,dz = 0。

(ii)(ジョルダンの補題) λ>0\lambda > 0 とする。R>R0R > R_0 のとき

∣∫ΓRf(z)eiλz dz∣≤πλM(R)\left|\int_{\Gamma_R}f(z)e^{i\lambda z}\,dz\right| \le \frac{\pi}{\lambda}M(R)

が成り立つ。特に lim⁡R→∞M(R)=0\lim_{R\to\infty}M(R) = 0 ならば lim⁡R→∞∫ΓRf(z)eiλzdz=0\displaystyle\lim_{R\to\infty}\int_{\Gamma_R}f(z)e^{i\lambda z}dz = 0。

証明(補題 6.1)

(i):ΓR\Gamma_R の長さは πR\pi R なので、§2 の ML 不等式より

∣∫ΓRf(z) dz∣≤πR⋅M(R)⟶0.\left|\int_{\Gamma_R}f(z)\,dz\right| \le \pi R \cdot M(R) \longrightarrow 0 .

(ii):z=Reiθ=Rcos⁡θ+iRsin⁡θz = Re^{i\theta} = R\cos\theta + iR\sin\theta に対し

∣eiλz∣=∣eiλRcos⁡θ∣⋅∣e−λRsin⁡θ∣=e−λRsin⁡θ|e^{i\lambda z}| = |e^{i\lambda R\cos\theta}|\cdot|e^{-\lambda R\sin\theta}| = e^{-\lambda R\sin\theta}

です(第一の因子は絶対値 11)。∣dz∣=R dθ|dz| = R\,d\theta に注意して

∣∫ΓRf(z)eiλzdz∣≤∫0πM(R) e−λRsin⁡θ R dθ=2RM(R)∫0π/2e−λRsin⁡θ 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

となります(sin⁡(π−θ)=sin⁡θ\sin(\pi-\theta) = \sin\theta による対称性で [0,π/2][0,\pi/2] の積分の 2 倍にしました)。

ここでジョルダンの不等式 sin⁡θ≥2θπ\sin\theta \ge \dfrac{2\theta}{\pi}(0≤θ≤π/20 \le \theta \le \pi/2)を使います。証明しておきます。ψ(θ)=sin⁡θ/θ\psi(\theta) = \sin\theta/\theta(θ>0\theta > 0)とおくと ψ′(θ)=θcos⁡θ−sin⁡θθ2\psi'(\theta) = \dfrac{\theta\cos\theta - \sin\theta}{\theta^2} です。分子を u(θ)=θcos⁡θ−sin⁡θu(\theta) = \theta\cos\theta - \sin\theta とすると u(0)=0u(0) = 0、u′(θ)=cos⁡θ−θsin⁡θ−cos⁡θ=−θsin⁡θ<0u'(\theta) = \cos\theta - \theta\sin\theta - \cos\theta = -\theta\sin\theta < 0(0<θ≤π/20 < \theta \le \pi/2)なので u<0u < 0、したがって ψ\psi は (0,π/2](0,\pi/2] で狭義単調減少です。よって ψ(θ)≥ψ(π/2)=1π/2=2π\psi(\theta) \ge \psi(\pi/2) = \dfrac{1}{\pi/2} = \dfrac{2}{\pi}、すなわち sin⁡θ≥2θ/π\sin\theta \ge 2\theta/\pi です。

これを使うと e−λRsin⁡θ≤e−2λRθ/πe^{-\lambda R\sin\theta} \le e^{-2\lambda R\theta/\pi} なので

∫0π/2e−λRsin⁡θdθ≤∫0π/2e−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}

となり、これを代入して 2RM(R)⋅π2λR=πλM(R)2RM(R)\cdot\dfrac{\pi}{2\lambda R} = \dfrac{\pi}{\lambda}M(R) を得ます。

命題 6.2(有理関数の広義積分)✓

P,QP, Q を複素係数の多項式とし、次を仮定する。

  • deg⁡Q≥deg⁡P+2\deg Q \ge \deg P + 2。
  • QQ は実軸上に零点を持たない。

このとき f=P/Qf = P/Q に対し ∫−∞∞f(x) dx\int_{-\infty}^{\infty}f(x)\,dx は絶対収束し、

∫−∞∞P(x)Q(x) dx=2πi∑Im⁡a>0Res⁡(f,a)\int_{-\infty}^{\infty}\frac{P(x)}{Q(x)}\,dx = 2\pi i \sum_{\operatorname{Im}a > 0}\operatorname{Res}(f, a)

が成り立つ。和は上半平面にある QQ のすべての零点 aa にわたる。

証明(命題 6.2)

p=deg⁡Pp = \deg P、q=deg⁡Qq = \deg Q、最高次係数を α,β\alpha, \beta とします。∣z∣→∞|z|\to\infty のとき z q−pP(z)/Q(z)→α/βz^{\,q-p}P(z)/Q(z) \to \alpha/\beta なので、ある R1≥1R_1 \ge 1 と K>0K > 0 があって

∣z∣≥R1  ⟹  ∣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}

です(最後の不等号に q−p≥2q - p \ge 2 と ∣z∣≥R1≥1|z| \ge R_1 \ge 1 を使いました)。したがって ∫∣x∣≥R1∣f(x)∣dx≤2∫R1∞Kx−2dx<∞\int_{|x| \ge R_1}|f(x)|dx \le 2\int_{R_1}^{\infty}Kx^{-2}dx < \infty であり、[−R1,R1][-R_1,R_1] 上では ff が連続(QQ は実軸上で 00 にならない)なので、広義積分は絶対収束します。また R≥R1R \ge R_1 に対し R M(R)≤K/R→0R\,M(R) \le K/R \to 0 なので 補題 6.1 (i) の仮定が満たされます。

QQ の零点は有限個です。R>max⁡(R1,max⁡a∣a∣)R > \max(R_1, \max_a|a|) を取り、CRC_R を線分 [−R,R][-R,R] と半円弧 ΓR\Gamma_R をつないだ正の向きの単純閉曲線とします。CRC_R とその内部は単連結領域 C\mathbb{C} に含まれ、ff は QQ の零点を除いて正則で、CRC_R の内部にある零点はちょうど上半平面の零点全部です(実軸上には零点がなく、RR の取り方から下半平面の零点は内部に入りません)。定理 5.1 より

∫−RRf(x) dx+∫ΓRf(z) dz=2πi∑Im⁡a>0Res⁡(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)

です。右辺は RR によりません。R→∞R \to \infty とすると、左辺第 1 項は絶対収束性より ∫−∞∞f(x)dx\int_{-\infty}^{\infty}f(x)dx に収束し、第 2 項は 補題 6.1 (i) より 00 に収束します。

例 6.3(分母が 4 次の場合)

∫−∞∞dxx4+1\displaystyle\int_{-\infty}^{\infty}\frac{dx}{x^4+1} を計算します。P=1P = 1、Q=z4+1Q = z^4+1 で deg⁡Q−deg⁡P=4≥2\deg Q - \deg P = 4 \ge 2、Q(x)=x4+1≥1>0Q(x) = x^4 + 1 \ge 1 > 0 より実軸上に零点はないので 命題 6.2 が使えます。

z4=−1=eiπz^4 = -1 = e^{i\pi} の解は zk=eiπ(2k+1)/4z_k = e^{i\pi(2k+1)/4}(k=0,1,2,3k = 0,1,2,3)、すなわち eiπ/4, e3iπ/4, e5iπ/4, e7iπ/4e^{i\pi/4},\ e^{3i\pi/4},\ e^{5i\pi/4},\ e^{7i\pi/4} です。上半平面にあるのは偏角が π/4\pi/4 と 3π/43\pi/4 の二つです。いずれも QQ の単純零点(Q′=4z3Q' = 4z^3 は 00 にならない)なので、命題 4.2 (3) より、aa をそのどちらかとして

Res⁡(f,a)=14a3=a4a4=a4⋅(−1)=−a4\operatorname{Res}(f,a) = \frac{1}{4a^3} = \frac{a}{4a^4} = \frac{a}{4\cdot(-1)} = -\frac{a}{4}

です(a4=−1a^4 = -1 を使って a3a^3 を aa に直しました。この書き換えで計算がぐっと楽になります)。したがって

∑Res⁡=−14(eiπ/4+e3iπ/4)=−14(22(1+i)+22(−1+i))=−14⋅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

となり、

∫−∞∞dxx4+1=2πi⋅(−2i4)=2π24=π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∣≥1|x| \ge 1 では 1x4+1≤1x2+1\frac{1}{x^4+1} \le \frac{1}{x^2+1}、∣x∣≤1|x| \le 1 では 1x4+1≤1\frac{1}{x^4+1} \le 1 なので、積分値は 2+∫∣x∣≥1dxx2+1=2+π2≈3.572 + \int_{|x|\ge1}\frac{dx}{x^2+1} = 2 + \frac{\pi}{2} \approx 3.57 以下と分かり、2.222.22 という値と矛盾しません。

7. 三角関数の周期積分:単位円を回る

Section titled “7. 三角関数の周期積分:単位円を回る”

θ\theta について周期 2π2\pi の三角有理式の積分 ∫02πR(cos⁡θ,sin⁡θ) dθ\int_0^{2\pi}R(\cos\theta,\sin\theta)\,d\theta は、実積分が最初から閉じた区間を一周しているので、輪郭を付け足す必要すらありません。z=eiθz = e^{i\theta} と置けば、これはそのまま単位円周に沿う複素積分になります。オイラーの公式から

cos⁡θ=z+z−12,sin⁡θ=z−z−12i,dθ=dziz\cos\theta = \frac{z + z^{-1}}{2},\qquad \sin\theta = \frac{z - z^{-1}}{2i},\qquad d\theta = \frac{dz}{iz}

です。最後の関係は、単位円周を γ(θ)=eiθ\gamma(\theta) = e^{i\theta}(0≤θ≤2π0 \le \theta \le 2\pi)で径数付けると γ′(θ)=ieiθ=iz\gamma'(\theta) = ie^{i\theta} = iz なので、任意の連続関数 FF に対し

∮∣z∣=1F(z)iz dz=∫02πF(eiθ)ieiθ ieiθ dθ=∫02πF(eiθ) 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

が成り立つ、という意味です。

命題 7.1(a + b cosθ 型の積分)✓

a,ba, b を実数とし a>∣b∣≥0a > |b| \ge 0 とする。このとき a+bcos⁡θ≥a−∣b∣>0a + b\cos\theta \ge a - |b| > 0 で被積分関数は連続であり、

∫02πdθa+bcos⁡θ=2πa2−b2.\int_0^{2\pi}\frac{d\theta}{a + b\cos\theta} = \frac{2\pi}{\sqrt{a^2-b^2}} .
証明(命題 7.1)

b=0b = 0 のときは左辺 =2π/a= 2\pi/a、右辺 =2π/a2=2π/a= 2\pi/\sqrt{a^2} = 2\pi/a(a>0a > 0)で成立します。b<0b < 0 のときは θ↦θ+π\theta \mapsto \theta + \pi と置換すると cos⁡(θ+π)=−cos⁡θ\cos(\theta+\pi) = -\cos\theta で、被積分関数は周期 2π2\pi なので積分区間は [0,2π][0,2\pi] のままでよく、bb を −b-b に取り替えた式に帰着します。右辺は bb の符号によらないので、以下 b>0b > 0 とします。

z=eiθz = e^{i\theta} と置く上の対応により

I=∮∣z∣=11a+b2(z+1z)⋅dziz=∮∣z∣=12 dzi(bz2+2az+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)} .

分母の 2 次式 bz2+2az+bbz^2+2az+b の零点は

z±=−a±a2−b2bz_{\pm} = \frac{-a \pm \sqrt{a^2-b^2}}{b}

で、a>b>0a > b > 0 より a2−b2>0a^2 - b^2 > 0 なのでどちらも実数、かつ相異なります。解と係数の関係から z+z−=b/b=1z_+z_- = b/b = 1 です。

∣z+∣<1|z_+| < 1 を確かめます。z+=−a+a2−b2bz_+ = \dfrac{-a+\sqrt{a^2-b^2}}{b} は a2−b2<a\sqrt{a^2-b^2} < a より負で、∣z+∣=a−a2−b2b|z_+| = \dfrac{a - \sqrt{a^2-b^2}}{b} です。∣z+∣<1  ⟺  a−b<a2−b2|z_+| < 1 \iff a - b < \sqrt{a^2-b^2} であり、両辺とも正(a>ba>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。これは仮定です。よって ∣z+∣<1|z_+|<1、そして ∣z−∣=1/∣z+∣>1|z_-| = 1/|z_+| > 1 です。特に単位円周上に零点はありません。

したがって単位円の内部にある特異点は単純極 z+z_+ だけです。g(z)=2i(bz2+2az+b)g(z) = \dfrac{2}{i(bz^2+2az+b)} とすると、命題 4.2 (2) より

Res⁡(g,z+)=lim⁡z→z+2(z−z+)i b(z−z+)(z−z−)=2i b(z+−z−)=2i b⋅2a2−b2b=1ia2−b2.\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}} .

定理 5.1 より I=2πi⋅1ia2−b2=2πa2−b2I = 2\pi i \cdot \dfrac{1}{i\sqrt{a^2-b^2}} = \dfrac{2\pi}{\sqrt{a^2-b^2}} です。

たとえば ∫02πdθ5+4cos⁡θ=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} です。次は分子に cos⁡2θ\cos 2\theta が入る場合で、00 に多重極が現れます。

例 7.2(分子に cos 2θ がある積分)

I=∫02πcos⁡2θ5−4cos⁡θ dθ\displaystyle I = \int_0^{2\pi}\frac{\cos2\theta}{5-4\cos\theta}\,d\theta を計算します。5−4cos⁡θ≥1>05 - 4\cos\theta \ge 1 > 0 なので被積分関数は連続です。z=eiθz = e^{i\theta} とすると cos⁡2θ=z2+z−22\cos2\theta = \frac{z^2+z^{-2}}{2}、5−4cos⁡θ=5−2(z+z−1)=−2z2+5z−2z5-4\cos\theta = 5 - 2(z+z^{-1}) = \dfrac{-2z^2+5z-2}{z} なので

I=∮∣z∣=1(z2+z−2)/2(−2z2+5z−2)/z⋅dziz=∮∣z∣=1z4+12z2 i (−2z2+5z−2) dz=i2∮∣z∣=1z4+1z2(2z−1)(z−2) dzI = \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

です(−2z2+5z−2=−(2z−1)(z−2)-2z^2+5z-2 = -(2z-1)(z-2) と 1−2i=i2\frac{1}{-2i} = \frac{i}{2} を使いました)。被積分関数を gg とすると、単位円の内部の特異点は z=0z=0(22 位の極)と z=1/2z=1/2(単純極)で、z=2z=2 は外部です。

z=1/2z = 1/2:命題 4.2 (2) より

Res⁡(g,12)=z4+1z2⋅2(z−2)∣z=1/2=116+114⋅2(12−2)=17/16−3/4=−1712.\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}.

z=0z = 0:m=2m=2 として φ(z)=z2g(z)=z4+12z2−5z+2\varphi(z) = z^2 g(z) = \dfrac{z^4+1}{2z^2-5z+2} とおくと

φ′(z)=4z3(2z2−5z+2)−(z4+1)(4z−5)(2z2−5z+2)2,φ′(0)=0−1⋅(−5)4=54,\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},

なので Res⁡(g,0)=φ′(0)/1!=5/4\operatorname{Res}(g,0) = \varphi'(0)/1! = 5/4 です。念のためローラン展開でも確かめます。12−5z+2z2=12⋅11−5z−2z22=12(1+5z2+O(z2))=12+54z+O(z2)\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) なので、(z4+1)(z^4+1) を掛けて z2z^2 で割ると g(z)=12z2+54z+⋯g(z) = \dfrac{1}{2z^2} + \dfrac{5}{4z} + \cdots となり、z−1z^{-1} の係数は 5/45/4 で一致します。

したがって留数の和は −1712+54=−1712+1512=−16-\frac{17}{12} + \frac{5}{4} = -\frac{17}{12}+\frac{15}{12} = -\frac16 で、定理 5.1 より

I=i2⋅2πi⋅(−16)=−π⋅(−16)=π6.I = \frac{i}{2}\cdot 2\pi i \cdot\left(-\frac16\right) = -\pi\cdot\left(-\frac16\right) = \frac{\pi}{6}.

∫−∞∞f(x)cos⁡λx dx\int_{-\infty}^{\infty}f(x)\cos\lambda x\,dx のような積分では、cos⁡λz\cos\lambda z をそのまま複素平面に延長してはいけません。上半平面では cos⁡(iy)=cosh⁡y\cos(iy) = \cosh y が y→∞y \to \infty で指数的に発散するからです。代わりに eiλze^{i\lambda z} を使います。これは上半平面で ∣eiλz∣=e−λIm⁡z≤1|e^{i\lambda z}| = e^{-\lambda\operatorname{Im}z} \le 1(λ>0\lambda>0)と抑えられ、しかも実軸上では cos⁡λx+isin⁡λx\cos\lambda x + i\sin\lambda x なので、最後に実部と虚部を取れば目的の積分が得られます。

例 8.1(ローレンツ型のフーリエ変換)

∫−∞∞cos⁡xx2+1 dx\displaystyle\int_{-\infty}^{\infty}\frac{\cos x}{x^2+1}\,dx を求めます。f(z)=1z2+1f(z) = \dfrac{1}{z^2+1}、λ=1\lambda = 1 として F(z)=f(z)eizF(z) = f(z)e^{iz} を 例 6.3 と同じ半円輪郭 CRC_R(R>1R>1)で積分します。FF の特異点は z=±iz = \pm i の単純極で、上半平面にあるのは z=iz = i だけです。命題 4.2 (3) を g(z)=eizg(z) = e^{iz}、h(z)=z2+1h(z)=z^2+1、h′(z)=2zh'(z) = 2z に適用して

Res⁡(F,i)=ei⋅i2i=e−12i.\operatorname{Res}(F, i) = \frac{e^{i\cdot i}}{2i} = \frac{e^{-1}}{2i}.

定理 5.1 より ∫−RReixx2+1dx+∫ΓRF(z) dz=2πi⋅e−12i=π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} です。

円弧の評価には 補題 6.1 (ii) を使います。∣z∣=R>1|z| = R > 1 のとき ∣f(z)∣=1∣z2+1∣≤1R2−1|f(z)| = \dfrac{1}{|z^2+1|} \le \dfrac{1}{R^2-1}(三角不等式 ∣z2+1∣≥∣z∣2−1|z^2+1| \ge |z|^2 - 1)なので M(R)≤1R2−1→0M(R) \le \dfrac{1}{R^2-1}\to0 であり、∣∫ΓRF∣≤π1⋅1R2−1→0\left|\int_{\Gamma_R}F\right| \le \dfrac{\pi}{1}\cdot\dfrac{1}{R^2-1}\to0 です。よって

∫−∞∞eixx2+1 dx=πe.\int_{-\infty}^{\infty}\frac{e^{ix}}{x^2+1}\,dx = \frac{\pi}{e} .

実部を取ると ∫−∞∞cos⁡xx2+1dx=πe\displaystyle\int_{-\infty}^{\infty}\frac{\cos x}{x^2+1}dx = \frac{\pi}{e}、虚部を取ると ∫−∞∞sin⁡xx2+1dx=0\displaystyle\int_{-\infty}^{\infty}\frac{\sin x}{x^2+1}dx = 0 です。後者は被積分関数が奇関数であることからも分かり、検算になります。なお cos⁡xx2+1\dfrac{\cos x}{x^2+1} の絶対値は 1x2+1\dfrac{1}{x^2+1} 以下なので、この広義積分は絶対収束しています。

被積分関数が実軸上に極を持つと、広義積分は普通の意味では発散します。それでも「対称に切り取れば」極限が存在することがあります。

定義 8.2(コーシーの主値)↗

c∈Rc \in \mathbb{R} とし、ff は R∖{c}\mathbb{R}\setminus\{c\} で連続とする。極限

p.v.∫−∞∞f(x) dx:=lim⁡R→∞ε→+0(∫−Rc−εf(x) dx+∫c+εRf(x) dx)\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)

が存在するとき、その値を ff のコーシーの主値と呼ぶ。左右を同じ ε\varepsilon で切り落とすのが要点である。通常の広義積分が存在するときは、主値も存在して同じ値になる。

輪郭のほうは、実軸上の極を半径 ε\varepsilon の小円弧で迂回して避けます。その小円弧の寄与を測るのが次の補題です。「一周すれば 2πi×2\pi i \times 留数」だったものが、「角度 β−α\beta-\alpha だけ回れば i(β−α)×i(\beta-\alpha)\times 留数」になる、と憶えると自然です。

補題 8.3(単純極まわりの小円弧(くぼみ補題))✓

cc を ff の単純極とし、ρ=Res⁡(f,c)\rho = \operatorname{Res}(f,c) とおく。実数 α,β\alpha, \beta に対し、γε\gamma_\varepsilon を z=c+εeiθz = c + \varepsilon e^{i\theta}(θ\theta が α\alpha から β\beta まで動く)で定まる円弧とする。このとき

lim⁡ε→+0∫γεf(z) dz=i(β−α)ρ.\lim_{\varepsilon \to +0}\int_{\gamma_\varepsilon}f(z)\,dz = i(\beta - \alpha)\rho .

(β>α\beta > \alpha なら反時計回り、β<α\beta < \alpha なら時計回りである。)

証明(補題 8.3)

cc は単純極なので、ローラン展開は f(z)=ρz−c+∑n≥0cn(z−c)nf(z) = \dfrac{\rho}{z-c} + \sum_{n \ge 0}c_n(z-c)^n の形です。第 2 項を g(z)g(z) とおくと、gg は D(c,R)D(c,R) で正則なべき級数の和ですから、ある r>0r>0 と M>0M>0 に対し ∣z−c∣≤r|z-c| \le r で ∣g(z)∣≤M|g(z)| \le M です(連続関数のコンパクト集合上の有界性)。

主要部の寄与は ε\varepsilon によらず計算できます。z=c+εeiθz = c+\varepsilon e^{i\theta}、dz=iεeiθdθdz = i\varepsilon e^{i\theta}d\theta より

∫γερz−c dz=ρ∫αβiεeiθεeiθ 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).

残りは ML 不等式で消えます。γε\gamma_\varepsilon の長さは ∣β−α∣ε|\beta-\alpha|\varepsilon なので、ε≤r\varepsilon \le r のとき

∣∫γεg(z) dz∣≤∣β−α∣ ε M⟶0(ε→+0).\left|\int_{\gamma_\varepsilon}g(z)\,dz\right| \le |\beta-\alpha|\,\varepsilon\, M \longrightarrow 0 \qquad (\varepsilon \to +0).

両者を足して結論を得ます。

次の例では、原点を小さく避ける輪郭を使います。向きは、実軸上の二本の線分は左から右、大きな半円弧は反時計回り、原点を避ける小円弧は時計回りです。小円弧が時計回りなので、くぼみ補題の寄与は負の符号(−iπ×-i\pi\times 留数)で出ます。

例 8.4(ディリクレ積分)

∫0∞sin⁡xx dx=π2\displaystyle\int_0^{\infty}\frac{\sin x}{x}\,dx = \frac{\pi}{2} を示します。f(z)=eizzf(z) = \dfrac{e^{iz}}{z} とし、0<ε<R0 < \varepsilon < R に対して次の四つをつないだ輪郭 CC を考えます。線分 [ε,R][\varepsilon, R]、大円弧 ΓR\Gamma_R(反時計回りに 00 から π\pi へ)、線分 [−R,−ε][-R,-\varepsilon]、小円弧 γε\gamma_\varepsilon(θ\theta が π\pi から 00 へ、すなわち時計回り)。

CC とその内部は、複素平面から半直線 {iy:y≤0}\{iy : y \le 0\} を除いた単連結領域 DD に含まれ、ff は DD で正則です(00 は DD に属しません)。よって 定理 5.1 を特異点 00 個の場合に使って(これはコーシーの積分定理そのものです)

∫εReixxdx+∫ΓRf dz+∫−R−εeixxdx+∫γεf dz=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 .

第 3 項は x=−tx = -t と置換すると ∫−R−εeixxdx=∫Rεe−it−t(−dt)=−∫εRe−ittdt\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 なので、第 1 項と合わせて

∫εReix−e−ixx dx=2i∫εRsin⁡xx dx\int_\varepsilon^R\frac{e^{ix}-e^{-ix}}{x}\,dx = 2i\int_\varepsilon^R\frac{\sin x}{x}\,dx

になります。第 2 項は 補題 6.1 (ii) で f(z)=1/zf(z) = 1/z、λ=1\lambda = 1 とすると M(R)=1/R→0M(R) = 1/R \to 0 なので 00 に収束します。第 4 項は、00 が ff の単純極で Res⁡(f,0)=lim⁡z→0z⋅eizz=1\operatorname{Res}(f,0) = \lim_{z\to0}z\cdot\frac{e^{iz}}{z} = 1 ですから、補題 8.3 を α=π\alpha = \pi、β=0\beta = 0 として適用し i(0−π)⋅1=−iπi(0-\pi)\cdot1 = -i\pi に収束します。

以上を合わせて ε→+0\varepsilon \to +0、R→∞R \to \infty とすると

2i∫0∞sin⁡xx dx−iπ=0,すなわち∫0∞sin⁡xx dx=π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}.

同じ計算を実部・虚部に分けて読むと p.v.∫−∞∞eixxdx=iπ\mathrm{p.v.}\displaystyle\int_{-\infty}^{\infty}\frac{e^{ix}}{x}dx = i\pi であり、実部から p.v.∫−∞∞cos⁡xxdx=0\mathrm{p.v.}\displaystyle\int_{-\infty}^{\infty}\frac{\cos x}{x}dx = 0(奇関数なので当然です)、虚部から上の結果が出ます。

注意すべき点が二つあります。第一に、sin⁡xx\dfrac{\sin x}{x} は x=0x=0 で除去可能特異点を持つ(例 3.6 (a))ので ε→0\varepsilon \to 0 の極限は問題なく、ε\varepsilon を入れたのは eiz/ze^{iz}/z を扱うための便宜です。第二に、この広義積分は絶対収束しません(命題 2.2[ボールウェイン積分])。実際 ∫0∞∣sin⁡xx∣dx=∞\int_0^{\infty}\left|\frac{\sin x}{x}\right|dx = \infty であり、値は R→∞R \to \infty の極限としてのみ意味を持ちます。実解析の側からの別証明は 定理 3.2[ボールウェイン積分] にあり、この繊細さが表に出る有名な例として ボールウェイン積分 があります。

9. 鍵穴輪郭と分数べきを含む積分

Section titled “9. 鍵穴輪郭と分数べきを含む積分”

xa−1x^{a-1}(aa が整数でない)は、原点のまわりを一周すると値が変わる多価関数です。そこで正の実軸に切れ目を入れ、その両岸を往復する鍵穴輪郭を使います。付録の計算(∫0∞xa−11+xdx\int_0^\infty\frac{x^{a-1}}{1+x}dx)を、有理関数 gg 一般に拡張します。

C∖[0,∞)\mathbb{C}\setminus[0,\infty) の上で Log⁡z=ln⁡∣z∣+iarg⁡z\operatorname{Log}z=\ln|z|+i\arg z(arg⁡z∈(0,2π)\arg z\in(0,2\pi))と定め、za−1=e(a−1)Log⁡zz^{a-1}=e^{(a-1)\operatorname{Log}z} とします。Log⁡\operatorname{Log} は C∖[0,∞)\mathbb{C}\setminus[0,\infty) で正則で、正の実軸の上岸では arg⁡→0\arg\to0、下岸では arg⁡→2π\arg\to2\pi です。

z₁z₂|z| = R|z| = ε上岸:arg z → 0下岸:arg z → 2π
鍵穴輪郭。切れ目(正の実軸)の上岸を右へ、大円を反時計回りに、下岸を左へ、小円を時計回りに回ります。見やすいように、切れ目との隙間(δ)は大きく描いています。×は被積分関数の極です。

定理 9.1(鍵穴輪郭による分数べきの積分)✓

gg を有理関数で、[0,∞)[0,\infty) 上に極をもたず、ある整数 m≥1m\ge1 と定数 CC、R0>0R_0>0 について ∣g(z)∣≤C∣z∣−m|g(z)|\le C|z|^{-m}(∣z∣≥R0|z|\ge R_0)を満たすものとする。aa は 0<a<m0<a<m を満たす整数でない実数とし、gg の極を z1,…,zpz_1,\dots,z_p とする。このとき広義積分は絶対収束し

∫0∞xa−1g(x) dx=2πi1−e2πia∑k=1pRes⁡(za−1g(z), zk)\int_0^\infty x^{a-1}g(x)\,dx=\frac{2\pi i}{1-e^{2\pi ia}}\sum_{k=1}^p\operatorname{Res}\bigl(z^{a-1}g(z),\,z_k\bigr)

が成り立つ。

証明の方針は次のとおりです。切れ目の上岸と下岸では za−1z^{a-1} の偏角が 00 と 2π2\pi で、値が e2πiae^{2\pi ia} 倍だけ違うので、両岸の積分は打ち消し合わず (1−e2πia)J(1-e^{2\pi ia})J として残ります。輪郭は切れ目から少し(δ\delta)ずらして閉じ、先に δ→0\delta\to0、次に ε→0\varepsilon\to0 と R→∞R\to\infty の極限を取ります。

証明(定理 9.1)

f(z)=za−1g(z)f(z)=z^{a-1}g(z) は C∖[0,∞)\mathbb{C}\setminus[0,\infty) で zkz_k を除いて正則です。絶対収束:gg は 00 の近くで有界(∣g∣≤M|g|\le M)なので ∣xa−1g∣≤Mxa−1|x^{a-1}g|\le Mx^{a-1} は a>0a>0 で可積分、xx が大きいところでは ∣xa−1g∣≤Cxa−1−m|x^{a-1}g|\le Cx^{a-1-m} が a<ma<m で可積分です。

輪郭. 0<ε<min⁡k∣zk∣0<\varepsilon<\min_k|z_k|、R>max⁡(R0,max⁡k∣zk∣)R>\max(R_0,\max_k|z_k|) とし、δ>0\delta>0 を、すべての zkz_k の偏角が (δ,2π−δ)(\delta,2\pi-\delta) に入るほど小さく取ります。次の 4 つをつないだ閉曲線 Γ\Gamma を考えます:(1) 半直線 arg⁡z=δ\arg z=\delta 上を εeiδ\varepsilon e^{i\delta} から ReiδRe^{i\delta} へ、(2) 円 ∣z∣=R|z|=R を偏角 δ\delta から 2π−δ2\pi-\delta まで反時計回りに(CRC_R)、(3) 半直線 arg⁡z=2π−δ\arg z=2\pi-\delta 上を Rei(2π−δ)Re^{i(2\pi-\delta)} から εei(2π−δ)\varepsilon e^{i(2\pi-\delta)} へ、(4) 円 ∣z∣=ε|z|=\varepsilon を偏角 2π−δ2\pi-\delta から δ\delta まで時計回りに(CεC_\varepsilon)。これは区分的に C1C^1 な単純閉曲線で、その内部(円環の切れた扇形)と Γ\Gamma は D=C∖[0,∞)D=\mathbb{C}\setminus[0,\infty) に含まれ、内部にすべての zkz_k を含みます。定理 5.1 により

∮Γf dz=2πi∑kRes⁡(f,zk)\oint_\Gamma f\,dz=2\pi i\sum_k\operatorname{Res}(f,z_k)

です。∮Γf dz=I1(δ)+IR(δ)+I3(δ)+Iε(δ)\oint_\Gamma f\,dz=I_1(\delta)+I_R(\delta)+I_3(\delta)+I_\varepsilon(\delta)((1)〜(4) の寄与)とします。

δ→0\delta\to0. f(reiδ)eiδ=ra−1eiaδg(reiδ)f(re^{i\delta})e^{i\delta}=r^{a-1}e^{ia\delta}g(re^{i\delta}) は、gg が [ε,R][\varepsilon,R] を含む開集合で連続(極がない)ことから、r∈[ε,R]r\in[\varepsilon,R] について ra−1g(r)r^{a-1}g(r) に一様収束します。よって I1(δ)→J:=∫εRra−1g(r) drI_1(\delta)\to J:=\int_\varepsilon^Rr^{a-1}g(r)\,dr です。同様に f(rei(2π−δ))=ra−1e(a−1)i(2π−δ)g(re−iδ)→e2πi(a−1)ra−1g(r)=e2πiara−1g(r)f(re^{i(2\pi-\delta)})=r^{a-1}e^{(a-1)i(2\pi-\delta)}g(re^{-i\delta})\to e^{2\pi i(a-1)}r^{a-1}g(r)=e^{2\pi ia}r^{a-1}g(r) (r∈[ε,R]r\in[\varepsilon,R] について一様)で、dz=ei(2π−δ)drdz=e^{i(2\pi-\delta)}dr の因子は 11 に収束し、向きが逆なので I3(δ)→−e2πiaJI_3(\delta)\to-e^{2\pi ia}J です。円弧は δ\delta によらない評価ができます。∣za−1∣=∣z∣a−1|z^{a-1}|=|z|^{a-1} なので ML 不等式から ∣IR(δ)∣≤2πR⋅Ra−1⋅CR−m=2πCRa−m|I_R(\delta)|\le2\pi R\cdot R^{a-1}\cdot CR^{-m}=2\pi CR^{a-m}、∣Iε(δ)∣≤2πε⋅εa−1M=2πMεa|I_\varepsilon(\delta)|\le2\pi\varepsilon\cdot\varepsilon^{a-1}M=2\pi M\varepsilon^a です。したがって

∣2πi∑kRes⁡(f,zk)−(1−e2πia)J∣≤2πCRa−m+2πMεa\Bigl|2\pi i\sum_k\operatorname{Res}(f,z_k)-\bigl(1-e^{2\pi ia}\bigr)J\Bigr|\le2\pi CR^{a-m}+2\pi M\varepsilon^a

です。ここで ε→0\varepsilon\to0、R→∞R\to\infty とすると、a>0a>0 と a<ma<m から右辺は 00 に収束し、J→∫0∞xa−1g dxJ\to\int_0^\infty x^{a-1}g\,dx です。a∉Za\notin\mathbb{Z} なので 1−e2πia≠01-e^{2\pi ia}\ne0 で、主張を得ます。

例 9.2(π/sin πa と π/(2 sin(πa/2)))

(a) g(z)=11+zg(z)=\frac1{1+z}(m=1m=1)、0<a<10<a<1。極 −1=eiπ-1=e^{i\pi} での za−1gz^{a-1}g の留数は eiπ(a−1)=−eiπae^{i\pi(a-1)}=-e^{i\pi a} なので、定理 9.1 から ∫0∞xa−11+xdx=2πi1−e2πia(−eiπa)=πsin⁡πa\int_0^\infty\frac{x^{a-1}}{1+x}dx=\frac{2\pi i}{1-e^{2\pi ia}}(-e^{i\pi a})=\frac{\pi}{\sin\pi a} です(付録の計算と同じ)。

(b) g(z)=11+z2g(z)=\frac1{1+z^2}(m=2m=2)、0<a<20<a<2、a≠1a\ne1。極は i=eiπ/2i=e^{i\pi/2} と −i=e3iπ/2-i=e^{3i\pi/2} で、留数は eiπ(a−1)/22i\frac{e^{i\pi(a-1)/2}}{2i} と e3iπ(a−1)/2−2i\frac{e^{3i\pi(a-1)/2}}{-2i} です。θ=πa2\theta=\frac{\pi a}2 とおくと eiπ(a−1)/2=−ieiθe^{i\pi(a-1)/2}=-ie^{i\theta}、e3iπ(a−1)/2=ie3iθe^{3i\pi(a-1)/2}=ie^{3i\theta} なので、留数の和は 12i(−ieiθ−ie3iθ)=−e2iθcos⁡θ\frac1{2i}(-ie^{i\theta}-ie^{3i\theta})=-e^{2i\theta}\cos\theta、1−e2πia=1−e4iθ=−2ie2iθsin⁡2θ1-e^{2\pi ia}=1-e^{4i\theta}=-2ie^{2i\theta}\sin2\theta です。よって

∫0∞xa−11+x2dx=2πi (−e2iθcos⁡θ)−2ie2iθsin⁡2θ=πcos⁡θsin⁡2θ=π2sin⁡(πa/2)\int_0^\infty\frac{x^{a-1}}{1+x^2}dx=\frac{2\pi i\,(-e^{2i\theta}\cos\theta)}{-2ie^{2i\theta}\sin2\theta}=\frac{\pi\cos\theta}{\sin2\theta}=\frac{\pi}{2\sin(\pi a/2)}

です。a=1a=1 では直接 ∫0∞dx1+x2=π2\int_0^\infty\frac{dx}{1+x^2}=\frac\pi2 で、この式の a→1a\to1 の値と一致します。

aa が整数のときは 1−e2πia=01-e^{2\pi ia}=0 で上の公式が使えません。このとき、対数 Log⁡z\operatorname{Log}z を被積分関数にかけます。

命題 9.3(対数を使う鍵穴積分)✓

gg を有理関数で、[0,∞)[0,\infty) 上に極をもたず、∣g(z)∣≤C∣z∣−2|g(z)|\le C|z|^{-2}(∣z∣≥R0|z|\ge R_0)を満たすものとし、極を z1,…,zpz_1,\dots,z_p とする。このとき

∫0∞g(x) dx=−∑k=1pRes⁡(g(z)Log⁡z, zk)\int_0^\infty g(x)\,dx=-\sum_{k=1}^p\operatorname{Res}\bigl(g(z)\operatorname{Log}z,\ z_k\bigr)

が成り立つ(Log⁡z\operatorname{Log}z は arg⁡z∈(0,2π)\arg z\in(0,2\pi) の分枝)。

証明(命題 9.3)

定理 9.1 と同じ輪郭を F(z)=g(z)Log⁡zF(z)=g(z)\operatorname{Log}z に使います。上の半直線では Log⁡→ln⁡r\operatorname{Log}\to\ln r、下の半直線では ln⁡r+2πi\ln r+2\pi i なので、δ→0\delta\to0 で 2 本の線分の寄与の和は ∫εRg(r)ln⁡r dr−∫εRg(r)(ln⁡r+2πi) dr=−2πi∫εRg(r) dr\int_\varepsilon^Rg(r)\ln r\,dr-\int_\varepsilon^Rg(r)(\ln r+2\pi i)\,dr=-2\pi i\int_\varepsilon^Rg(r)\,dr に収束します。円弧は ∣Log⁡z∣≤∣ln⁡∣z∣∣+2π|\operatorname{Log}z|\le|\ln|z||+2\pi から、∣IR∣≤2πR⋅CR−2(ln⁡R+2π)→0|I_R|\le2\pi R\cdot CR^{-2}(\ln R+2\pi)\to0、∣Iε∣≤2πεM(∣ln⁡ε∣+2π)→0|I_\varepsilon|\le2\pi\varepsilon M(|\ln\varepsilon|+2\pi)\to0 です。定理 5.1 から −2πi∫0∞g dx=2πi∑kRes⁡(F,zk)-2\pi i\int_0^\infty g\,dx=2\pi i\sum_k\operatorname{Res}(F,z_k) です。

例 9.4(∫dx/(x³+1))

g(z)=1z3+1g(z)=\frac1{z^3+1} の極は zk=eiθkz_k=e^{i\theta_k}、θk=π3,π,5π3\theta_k=\frac\pi3,\pi,\frac{5\pi}3 で(すべて偏角が (0,2π)(0,2\pi) に入る)、Res⁡(gLog⁡z,zk)=Log⁡zk3zk2=iθk3e−2iθk\operatorname{Res}(g\operatorname{Log}z,z_k)=\frac{\operatorname{Log}z_k}{3z_k^2}=\frac{i\theta_k}3e^{-2i\theta_k} です。和は i3[π3e−2πi/3+π+5π3e2πi/3]\frac i3\bigl[\frac\pi3e^{-2\pi i/3}+\pi+\frac{5\pi}3e^{2\pi i/3}\bigr] で、e±2πi/3=−12±32ie^{\pm2\pi i/3}=-\frac12\pm\frac{\sqrt3}2i から括弧の中は 23π3i\frac{2\sqrt3\pi}{3}i となり、和は −23π9-\frac{2\sqrt3\pi}9 です。したがって 命題 9.3 より ∫0∞dxx3+1=23π9=2π33\int_0^\infty\frac{dx}{x^3+1}=\frac{2\sqrt3\pi}9=\frac{2\pi}{3\sqrt3} です。

留数定理は、整数にわたる和を、複素平面の極の寄与に置き換えるのにも使えます。πcot⁡πz\pi\cot\pi z は、すべての整数に留数 11 の単純極をもつ関数です。

整数すべてに極をもつ πcot⁡πz\pi\cot\pi z に有理関数をかけ、大きな正方形の周で積分します。そのためには、周上で ∣cot⁡πz∣|\cot\pi z| が NN によらない定数で抑えられている必要があります。頂点を (N+12)(±1±i)(N+\frac12)(\pm1\pm i) に取るのは、縦の辺が整数点(極)を通らないようにするためです。

補題 10.1(cot πz の一様な評価)✓

N≥1N\ge1 とし、CNC_N を 4 点 (N+12)(±1±i)(N+\frac12)(\pm1\pm i) を頂点とする正方形の周とする。CNC_N 上で ∣cot⁡πz∣≤coth⁡π2|\cot\pi z|\le\coth\frac\pi2 である(右辺は NN によらない定数)。

証明(補題 10.1)

cot⁡πz=i eiπz+e−iπzeiπz−e−iπz\cot\pi z=i\,\frac{e^{i\pi z}+e^{-i\pi z}}{e^{i\pi z}-e^{-i\pi z}} です。縦の辺 z=±(N+12)+iyz=\pm(N+\frac12)+iy:cot⁡\cot は周期 1 なので z=N+12+iyz=N+\frac12+iy で調べればよく、u=eiπ(N+1/2)=±iu=e^{i\pi(N+1/2)}=\pm i(u−1=−uu^{-1}=-u)として eiπz=ue−πye^{i\pi z}=ue^{-\pi y}、e−iπz=−ueπye^{-i\pi z}=-ue^{\pi y} なので cot⁡πz=ie−πy−eπye−πy+eπy=−itanh⁡πy\cot\pi z=i\frac{e^{-\pi y}-e^{\pi y}}{e^{-\pi y}+e^{\pi y}}=-i\tanh\pi y で、∣cot⁡πz∣≤1|\cot\pi z|\le1 です。横の辺 z=x+iyz=x+iy、∣y∣=N+12≥12|y|=N+\frac12\ge\frac12:∣eiπz∣=e−πy|e^{i\pi z}|=e^{-\pi y}、∣e−iπz∣=eπy|e^{-i\pi z}|=e^{\pi y} なので、y>0y>0 なら ∣cot⁡πz∣≤eπy+e−πyeπy−e−πy=coth⁡πy≤coth⁡π2|\cot\pi z|\le\frac{e^{\pi y}+e^{-\pi y}}{e^{\pi y}-e^{-\pi y}}=\coth\pi y\le\coth\frac\pi2(coth⁡\coth は (0,∞)(0,\infty) で減少)、y<0y<0 も同様です。coth⁡π2>1\coth\frac\pi2>1 なので、全体の上界は coth⁡π2\coth\frac\pi2 です。

定理 10.2(留数による級数の和)✓

f=P/Qf=P/Q を有理関数で deg⁡Q≥deg⁡P+2\deg Q\ge\deg P+2 を満たすものとし、ff の極を z1,…,zpz_1,\dots,z_p とする。このとき

∑n∈Zn∉{z1,…,zp}f(n)=−∑k=1pRes⁡(πcot⁡(πz) f(z), zk)\sum_{\substack{n\in\mathbb{Z}\\ n\notin\{z_1,\dots,z_p\}}}f(n)=-\sum_{k=1}^p\operatorname{Res}\bigl(\pi\cot(\pi z)\,f(z),\,z_k\bigr)

が成り立つ(左辺は絶対収束する)。

証明(定理 10.2)

h(z)=πcot⁡(πz)f(z)h(z)=\pi\cot(\pi z)f(z) とおきます。∣f(z)∣≤K∣z∣−2|f(z)|\le K|z|^{-2}(∣z∣|z| が大きいとき)なので左辺は絶対収束します。πcot⁡πz=πcos⁡πzsin⁡πz\pi\cot\pi z=\pi\frac{\cos\pi z}{\sin\pi z} は整数 nn で単純極をもち、sin⁡πz=(−1)nπ(z−n)+O((z−n)3)\sin\pi z=(-1)^n\pi(z-n)+O((z-n)^3)、cos⁡πn=(−1)n\cos\pi n=(-1)^n から πcot⁡πz=1z−n+(正則)\pi\cot\pi z=\frac1{z-n}+(\text{正則})、留数 11 です。したがって ff の極でない整数 nn では、hh は単純極で Res⁡(h,n)=f(n)\operatorname{Res}(h,n)=f(n) です。NN を、すべての zkz_k が CNC_N の内部にあり、CNC_N 上に ff の極がないほど大きく取ります(補題 10.1 の正方形)。開正方形 D={∣Re⁡z∣<N+1, ∣Im⁡z∣<N+1}D=\{|\operatorname{Re}z|<N+1,\ |\operatorname{Im}z|<N+1\} は CNC_N とその内部を含み、hh は DD で、CNC_N の内部にある有限個の点(∣n∣≤N|n|\le N の整数と zkz_k)を除いて正則です。定理 5.1 により

∮CNh dz=2πi[∑∣n∣≤Nn∉{zk}f(n)+∑k=1pRes⁡(h,zk)]\oint_{C_N}h\,dz=2\pi i\Bigl[\sum_{\substack{|n|\le N\\ n\notin\{z_k\}}}f(n)+\sum_{k=1}^p\operatorname{Res}(h,z_k)\Bigr]

です。補題 10.1 から CNC_N 上で ∣h∣≤πcoth⁡π2⋅K(N+12)−2|h|\le\pi\coth\frac\pi2\cdot K(N+\tfrac12)^{-2}、CNC_N の長さは 8(N+12)8(N+\frac12) なので、ML 不等式から ∣∮CNh∣≤8πKcoth⁡(π/2)N+1/2→0\bigl|\oint_{C_N}h\bigr|\le\frac{8\pi K\coth(\pi/2)}{N+1/2}\to0 です。N→∞N\to\infty として結論を得ます。

例 10.3(ζ(2) = π²/6)

f(z)=1z2f(z)=\frac1{z^2} は極 00(整数)のみをもち、定理 10.2 から ∑n≠01n2=−Res⁡(πcot⁡πzz2,0)\sum_{n\ne0}\frac1{n^2}=-\operatorname{Res}\bigl(\frac{\pi\cot\pi z}{z^2},0\bigr) です。sin⁡πz=πz(1−π2z26+⋯ )\sin\pi z=\pi z\bigl(1-\frac{\pi^2z^2}6+\cdots\bigr)、cos⁡πz=1−π2z22+⋯\cos\pi z=1-\frac{\pi^2z^2}2+\cdots から πcot⁡πz=1z(1−π2z22)(1+π2z26)+⋯=1z−π23z+O(z3)\pi\cot\pi z=\frac1z\bigl(1-\frac{\pi^2z^2}2\bigr)\bigl(1+\frac{\pi^2z^2}6\bigr)+\cdots=\frac1z-\frac{\pi^2}3z+O(z^3) です。したがって πcot⁡πzz2=1z3−π23z+⋯\frac{\pi\cot\pi z}{z^2}=\frac1{z^3}-\frac{\pi^2}{3z}+\cdots の留数は −π23-\frac{\pi^2}3 で、∑n≠01n2=π23\sum_{n\ne0}\frac1{n^2}=\frac{\pi^2}3、すなわち ζ(2)=∑n≥11n2=π26\zeta(2)=\sum_{n\ge1}\frac1{n^2}=\frac{\pi^2}6 です。Fourier 級数の Parseval の等式(例 5.2[フーリエ級数])による証明と同じ値です。

例 10.4(Σ 1/(n² + a²) = (π/a) coth πa)

a>0a>0 とし f(z)=1z2+a2f(z)=\frac1{z^2+a^2} とします。極 ±ia\pm ia は整数でなく、cot⁡(iπa)=−icoth⁡πa\cot(i\pi a)=-i\coth\pi a(cot⁡(ix)=cos⁡ixsin⁡ix=cosh⁡xisinh⁡x=−icoth⁡x\cot(ix)=\frac{\cos ix}{\sin ix}=\frac{\cosh x}{i\sinh x}=-i\coth x)から Res⁡(πcot⁡(πz)f, ia)=πcot⁡(iπa)2ia=−πcoth⁡πa2a\operatorname{Res}(\pi\cot(\pi z)f,\,ia)=\frac{\pi\cot(i\pi a)}{2ia}=-\frac{\pi\coth\pi a}{2a}、同様に −ia-ia で πcot⁡(−iπa)−2ia=−πcoth⁡πa2a\frac{\pi\cot(-i\pi a)}{-2ia}=-\frac{\pi\coth\pi a}{2a} です。定理 10.2 より

∑n∈Z1n2+a2=πacoth⁡πa,したがって∑n≥11n2+a2=πacoth⁡πa−12a2\sum_{n\in\mathbb{Z}}\frac1{n^2+a^2}=\frac\pi a\coth\pi a,\qquad\text{したがって}\qquad\sum_{n\ge1}\frac1{n^2+a^2}=\frac{\pi a\coth\pi a-1}{2a^2}

です。a→0a\to0 とすると、πacoth⁡πa=1+(πa)23+⋯\pi a\coth\pi a=1+\frac{(\pi a)^2}3+\cdots から右辺は π26\frac{\pi^2}6 に収束し、ζ(2)=π26\zeta(2)=\frac{\pi^2}6 と整合します。

演習 11.1易

f(z)=ezz2(z−1)f(z) = \dfrac{e^z}{z^2(z-1)} のすべての孤立特異点を求め、それぞれの位数と留数を求めてください。

解答

特異点は分母の零点 z=0z=0(22 位)と z=1z=1(11 位)です。eze^z はこれらの点で 00 にならないので、位数は分母のとおりです。実際 z2f(z)=ezz−1→−1≠0z^2f(z) = \frac{e^z}{z-1} \to -1 \ne 0、(z−1)f(z)=ezz2→e≠0(z-1)f(z) = \frac{e^z}{z^2}\to e \ne 0 なので 命題 3.4 より 00 は 22 位の極、11 は単純極です。

z=1z=1:命題 4.2 (2) より Res⁡(f,1)=lim⁡z→1ezz2=e\operatorname{Res}(f,1) = \lim_{z\to1}\frac{e^z}{z^2} = e。

z=0z=0:同 (1) で m=2m=2 とし、φ(z)=ezz−1\varphi(z) = \dfrac{e^z}{z-1} に対し

φ′(z)=ez(z−1)−ez(z−1)2=ez(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 .

したがって留数の和は e−2e - 2 です。

演習 11.2標準

∫0∞x2(x2+1)(x2+4) dx\displaystyle\int_0^{\infty}\frac{x^2}{(x^2+1)(x^2+4)}\,dx を留数定理で計算し、部分分数分解による初等的な計算と一致することを確かめてください。

解答

被積分関数は偶関数なので ∫0∞=12∫−∞∞\displaystyle\int_0^{\infty} = \frac12\int_{-\infty}^{\infty} です。f(z)=z2(z2+1)(z2+4)f(z) = \dfrac{z^2}{(z^2+1)(z^2+4)} は deg⁡\deg の差が 22、分母は実軸上で 00 にならないので 命題 6.2 が使えます。上半平面の極は単純極 z=iz=i と z=2iz=2i です。

Res⁡(f,i)=z2(z+i)(z2+4)∣z=i=−12i⋅3=−16i=i6,\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,2i)=z2(z2+1)(z+2i)∣z=2i=−4(−3)⋅4i=13i=−i3.\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}.

(1i=−i\frac1i = -i を使いました。)和は i6−i3=−i6\frac{i}{6}-\frac{i}{3} = -\frac{i}{6} なので

∫−∞∞f(x) dx=2πi(−i6)=π3,∫0∞f(x) dx=π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}.

初等的な確認:x2(x2+1)(x2+4)=Ax2+1+Bx2+4\dfrac{x^2}{(x^2+1)(x^2+4)} = \dfrac{A}{x^2+1}+\dfrac{B}{x^2+4} とおくと x2=A(x2+4)+B(x2+1)x^2 = A(x^2+4)+B(x^2+1) より A+B=1A+B=1、4A+B=04A+B=0、すなわち A=−13A = -\frac13、B=43B=\frac43 です。∫0∞dxx2+a2=π2a\int_0^{\infty}\frac{dx}{x^2+a^2} = \frac{\pi}{2a}(a>0a>0)を使うと

−13⋅π2+43⋅π4=−π6+π3=π6-\frac13\cdot\frac{\pi}{2} + \frac43\cdot\frac{\pi}{4} = -\frac{\pi}{6}+\frac{\pi}{3} = \frac{\pi}{6}

で一致します。

演習 11.3標準

∫02πdθ(2+cos⁡θ)2\displaystyle\int_0^{2\pi}\frac{d\theta}{(2+\cos\theta)^2} を計算してください。

解答

z=eiθz = e^{i\theta} と置くと 2+cos⁡θ=2+z+z−12=z2+4z+12z2+\cos\theta = 2 + \dfrac{z+z^{-1}}{2} = \dfrac{z^2+4z+1}{2z} なので

I=∮∣z∣=14z2(z2+4z+1)2⋅dziz=4i∮∣z∣=1z dz(z2+4z+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}.

z2+4z+1=0z^2+4z+1 = 0 の解は z0=−2+3≈−0.268z_0 = -2+\sqrt3 \approx -0.268 と z1=−2−3≈−3.73z_1 = -2-\sqrt3 \approx -3.73 で、単位円の内部にあるのは z0z_0 だけです(z0z1=1z_0z_1 = 1 なので一方が内側、他方が外側)。分母は (z−z0)2(z−z1)2(z-z_0)^2(z-z_1)^2 なので z0z_0 は 22 位の極です。m=2m=2 として 命題 4.2 (1) を使います。φ(z)=z(z−z1)2\varphi(z) = \dfrac{z}{(z-z_1)^2} とすると

φ′(z)=(z−z1)2−z⋅2(z−z1)(z−z1)4=(z−z1)−2z(z−z1)3=−z−z1(z−z1)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}.

z0−z1=23z_0 - z_1 = 2\sqrt3、−z0−z1=−(−2+3)−(−2−3)=4-z_0-z_1 = -(-2+\sqrt3)-(-2-\sqrt3) = 4 なので

Res⁡=φ′(z0)=4(23)3=4243=163.\operatorname{Res} = \varphi'(z_0) = \frac{4}{(2\sqrt3)^3} = \frac{4}{24\sqrt3} = \frac{1}{6\sqrt3}.

定理 5.1 より

I=4i⋅2πi⋅163=8π63=4π33=43 π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 .

検算:命題 7.1 の式 ∫02πdθa+bcos⁡θ=2π(a2−b2)−1/2\int_0^{2\pi}\frac{d\theta}{a+b\cos\theta} = 2\pi(a^2-b^2)^{-1/2} を aa で微分すると、左辺は −∫02πdθ(a+bcos⁡θ)2-\int_0^{2\pi}\frac{d\theta}{(a+b\cos\theta)^2}、右辺は −2πa(a2−b2)−3/2-2\pi a(a^2-b^2)^{-3/2} です。a=2a=2、b=1b=1 とすると 4π33/2=4π33\frac{4\pi}{3^{3/2}} = \frac{4\pi}{3\sqrt3} で一致します。

演習 11.4難

a>0a > 0 とするとき ∫0∞xsin⁡xx2+a2 dx=π2e−a\displaystyle\int_0^{\infty}\frac{x\sin x}{x^2+a^2}\,dx = \frac{\pi}{2}e^{-a} を示してください。zz2+a2\dfrac{z}{z^2+a^2} には 補題 6.1 (i) が使えないことに注意してください。

解答

f(z)=zz2+a2f(z) = \dfrac{z}{z^2+a^2}、F(z)=f(z)eizF(z) = f(z)e^{iz} とし、半径 R>aR > a の半円輪郭 CRC_R を取ります。FF の上半平面の特異点は単純極 z=iaz = ia だけで、命題 4.2 (3) を g(z)=zeizg(z) = ze^{iz}、h(z)=z2+a2h(z) = z^2+a^2、h′(z)=2zh'(z)=2z に適用して

Res⁡(F,ia)=ia ei⋅ia2ia=e−a2.\operatorname{Res}(F, ia) = \frac{ia\,e^{i\cdot ia}}{2ia} = \frac{e^{-a}}{2}.

定理 5.1 より ∫−RRxeixx2+a2dx+∫ΓRF dz=2πi⋅e−a2=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} です。

円弧の評価:∣z∣=R>a|z|=R>a のとき ∣f(z)∣≤RR2−a2|f(z)| \le \dfrac{R}{R^2-a^2} なので M(R)→0M(R) \to 0 ですが、R M(R)→1≠0R\,M(R) \to 1 \ne 0 なので 補題 6.1 (i) は使えません。指数因子を持つ (ii)(ジョルダンの補題)なら λ=1\lambda = 1 として

∣∫ΓRF(z) dz∣≤πM(R)≤πRR2−a2⟶0\left|\int_{\Gamma_R}F(z)\,dz\right| \le \pi M(R) \le \frac{\pi R}{R^2-a^2}\longrightarrow 0

と結論できます。よって lim⁡R→∞∫−RRxeixx2+a2dx=iπe−a\displaystyle\lim_{R\to\infty}\int_{-R}^{R}\frac{xe^{ix}}{x^2+a^2}dx = i\pi e^{-a}。虚部を取ると

∫−∞∞xsin⁡xx2+a2 dx=πe−a\int_{-\infty}^{\infty}\frac{x\sin x}{x^2+a^2}\,dx = \pi e^{-a}

であり、被積分関数は偶関数(奇関数どうしの積)なので ∫0∞=π2e−a\displaystyle\int_0^{\infty} = \frac{\pi}{2}e^{-a} です。なおこの積分は絶対収束せず、R→∞R \to \infty の対称極限として存在します。

演習 11.5標準

定理 10.2 で f(z)=1z4f(z)=\frac1{z^4} を使い、ζ(4)=π490\zeta(4)=\frac{\pi^4}{90} を示せ。

解答

cot⁡w\cot w の級数を求めます。sin⁡w=w(1−w26+w4120)+⋯\sin w=w\bigl(1-\frac{w^2}6+\frac{w^4}{120}\bigr)+\cdots、cos⁡w=1−w22+w424−⋯\cos w=1-\frac{w^2}2+\frac{w^4}{24}-\cdots で、u=w26−w4120u=\frac{w^2}6-\frac{w^4}{120} とすると 11−u=1+u+u2+⋯=1+w26+(−1120+136)w4+⋯=1+w26+7w4360+⋯\frac1{1-u}=1+u+u^2+\cdots=1+\frac{w^2}6+\bigl(-\frac1{120}+\frac1{36}\bigr)w^4+\cdots=1+\frac{w^2}6+\frac{7w^4}{360}+\cdots です。これに 1−w22+w4241-\frac{w^2}2+\frac{w^4}{24} をかけると w4w^4 の係数は 7360−112+124=−145\frac7{360}-\frac1{12}+\frac1{24}=-\frac1{45} で、cot⁡w=1w−w3−w345−⋯\cot w=\frac1w-\frac w3-\frac{w^3}{45}-\cdots です。w=πzw=\pi z として πcot⁡πz=1z−π2z3−π4z345−⋯\pi\cot\pi z=\frac1z-\frac{\pi^2z}3-\frac{\pi^4z^3}{45}-\cdots なので、πcot⁡πzz4\frac{\pi\cot\pi z}{z^4} の z−1z^{-1} の係数(留数)は −π445-\frac{\pi^4}{45} です。定理 10.2 より ∑n≠01n4=π445\sum_{n\ne0}\frac1{n^4}=\frac{\pi^4}{45}、ζ(4)=π490\zeta(4)=\frac{\pi^4}{90} です。

演習 11.6標準

定理 9.1 を使って ∫0∞x(1+x)2 dx=π2\int_0^\infty\frac{\sqrt x}{(1+x)^2}\,dx=\frac\pi2 を示せ。

解答

a=32a=\frac32、g(z)=1(1+z)2g(z)=\frac1{(1+z)^2}(m=2m=2、0<a<20<a<2、整数でない)です。極 −1=eiπ-1=e^{i\pi} は 2 位で、Res⁡(z1/2(1+z)2,−1)=ddzz1/2∣z=−1=12(−1)−1/2=12e−iπ/2=−i2\operatorname{Res}\bigl(\frac{z^{1/2}}{(1+z)^2},-1\bigr)=\frac{d}{dz}z^{1/2}\big|_{z=-1}=\frac12(-1)^{-1/2}=\frac12e^{-i\pi/2}=-\frac i2 です(z1/2=e12Log⁡zz^{1/2}=e^{\frac12\operatorname{Log}z}、arg⁡(−1)=π\arg(-1)=\pi)。1−e2πia=1−e3πi=21-e^{2\pi ia}=1-e^{3\pi i}=2 なので、定理 9.1 から ∫0∞x1/2(1+x)2dx=2πi2⋅(−i2)=π2\int_0^\infty\frac{x^{1/2}}{(1+x)^2}dx=\frac{2\pi i}{2}\cdot\bigl(-\frac i2\bigr)=\frac\pi2 です。ここで xa−1=x1/2x^{a-1}=x^{1/2} です。

  • 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(留数)」という語と計算法が導入された原典です。

Appendix: 分数べきを含む積分と鍵穴輪郭

Section titled “Appendix: 分数べきを含む積分と鍵穴輪郭”

扱う積分。0<a<10 < a < 1 に対して

∫0∞xa−11+x dx=πsin⁡πa\int_0^{\infty}\frac{x^{a-1}}{1+x}\,dx = \frac{\pi}{\sin \pi a}

を示します。xa−1x^{a-1} は aa が整数でないと一価な複素関数に延長できないので、これまでの輪郭は使えません。切断(分枝の切れ目)を正の実軸に置き、その両岸を往復する「鍵穴輪郭」を使います。

分枝を決める。Log⁡z=log⁡∣z∣+iarg⁡z\operatorname{Log}z = \log|z| + i\arg z、arg⁡z∈(0,2π)\arg z \in (0, 2\pi) と定め、za−1:=e(a−1)Log⁡zz^{a-1} := e^{(a-1)\operatorname{Log}z} とします。これは C\mathbb{C} から非負実軸を除いた領域で正則です。f(z)=za−11+zf(z) = \dfrac{z^{a-1}}{1+z} の特異点は z=−1=eiπz = -1 = e^{i\pi} の単純極だけです。

輪郭。0<ε<1<R0 < \varepsilon < 1 < R とし、切断の上岸を ε\varepsilon から RR まで進み、大円 ∣z∣=R|z|=R を反時計回りに一周し、下岸を RR から ε\varepsilon まで戻り、小円 ∣z∣=ε|z|=\varepsilon を時計回りに一周して閉じます。厳密には上岸・下岸を実軸から δ\delta だけずらした線分に取り、最後に δ→0\delta \to 0 とします。この閉曲線の内部にある特異点は z=−1z=-1 だけです。

円周上の評価。大円上では ∣za−1∣=Ra−1|z^{a-1}| = R^{a-1}、∣1+z∣≥R−1|1+z| \ge R-1 なので、ML 不等式より寄与は 2πR⋅Ra−1R−1∼2πRa−1→02\pi R\cdot\dfrac{R^{a-1}}{R-1} \sim 2\pi R^{a-1} \to 0(a<1a<1)です。小円上では 2πε⋅εa−11−ε∼2πεa→02\pi\varepsilon\cdot\dfrac{\varepsilon^{a-1}}{1-\varepsilon} \sim 2\pi\varepsilon^{a}\to 0(a>0a>0)です。0<a<10 < a < 1 という仮定は、この二つの評価を同時に成立させるためにあります。

両岸の差。上岸では arg⁡z→0\arg z \to 0 なので za−1→xa−1z^{a-1} \to x^{a-1}、下岸では arg⁡z→2π\arg z \to 2\pi なので za−1→xa−1e2πi(a−1)=xa−1e2πiaz^{a-1}\to x^{a-1}e^{2\pi i(a-1)} = x^{a-1}e^{2\pi ia} です。下岸は逆向きに進むので、二つの線分の寄与の和は I:=∫0∞xa−11+xdxI := \int_0^{\infty}\frac{x^{a-1}}{1+x}dx を使って (1−e2πia)I\left(1 - e^{2\pi i a}\right)I と書けます。切断があるおかげで両岸が相殺せず、求めたい積分が残るのがこの輪郭の要点です。

留数と結論。Res⁡(f,−1)=(−1)a−1=e(a−1)iπ=−eiπa\operatorname{Res}(f,-1) = (-1)^{a-1} = e^{(a-1)i\pi} = -e^{i\pi a} なので、定理 5.1 より (1−e2πia)I=2πi(−eiπa)\left(1-e^{2\pi ia}\right)I = 2\pi i\left(-e^{i\pi a}\right)、すなわち

I=−2πi eiπa1−e2πia=−2πie−iπa−eiπa=−2πi−2isin⁡π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}.

検算。a=1/2a = 1/2 とすると I=π/sin⁡(π/2)=πI = \pi/\sin(\pi/2) = \pi です。一方 x=t2x = t^2 と置換すれば

∫0∞dxx(1+x)=∫0∞2t dtt(1+t2)=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

で一致します。この積分はベータ関数・ガンマ関数の相反公式 Γ(a)Γ(1−a)=π/sin⁡πa\Gamma(a)\Gamma(1-a) = \pi/\sin\pi a の複素解析による証明にもなっています。

このテーマの 100 本ノック(100 問)

執筆: Claude Code(Opus 5.5)、Claude Code(Sonnet 5.5)・監査 5 回(Claude Code)・最終 2026-10-12

この記事の誤りを報告する ・運営: 夢現技研合同会社 ・料金 ・利用条件 ・特定商取引法に基づく表記

© 2026 夢現技研合同会社 ・コード例は MIT ライセンスです。