Skip to content

留数定理と定積分への応用:孤立特異点の分類から実積分の計算まで

Prerequisite:正則関数の強力な性質:無限回微分可能性から一致の定理まで

Raw

This content is not available in your language yet.

  • 正則関数を孤立特異点 aa のまわりでローラン展開したとき、(za)1(z-a)^{-1} の係数 c1c_{-1} だけが aa を囲む閉曲線に沿う積分に生き残ります。この c1c_{-1} を留数と呼びます。
  • 孤立特異点は主要部の形で三つに分かれます。除去可能特異点は「近くで有界」、mm 位の極は「(za)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} を含む積分はジョルダンの補題つきの半円、実軸上に単純極があるときは極を小円でくぼませて主値を取ります。

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

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

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

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

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

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

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

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

記号:領域とは連結開集合のことです。D(a,r)={zC:za<r}D(a,r) = \{z \in \mathbb{C} : |z-a| < r\}D(a,r)={z:0<za<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 で定義します。円周 za=r|z-a| = r には断らない限り正の向き(反時計回り)を与え、閉曲線に沿う積分は \oint と書きます。

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

γf(z)dzγf(z)dzL(γ)maxzγ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)|

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

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

f(z)=n=cn(za)n,cn=12πiζa=ρf(ζ)(ζa)n+1dζ(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_nCC の内部にあり互いに他の外部にある正の向きの単純閉曲線とします。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

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

Definition 3.1孤立特異点

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

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

Definition 3.2主要部と特異点の三分類

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

Pa(z)=n1cn(za)nP_a(z) = \sum_{n \le -1} c_n (z-a)^n

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

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

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

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リーマンの除去可能特異点定理

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

  • (i) aaff の除去可能特異点である。
  • (ii) f(a)f(a) の値を適当に定めると、ffD(a,R)D(a,R) 全体で正則になる。
  • (iii) ある 0<rR0 < r \le RM>0M > 0 があって、0<za<r0 < |z-a| < rf(z)M|f(z)| \le M が成り立つ。
  • (iv) limza(za)f(z)=0\lim_{z \to a}(z-a)f(z) = 0
Proof(Theorem 3.3)

(i) \Rightarrow (ii):主要部が 00 なので D(a,R)D'(a,R)f(z)=n0cn(za)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)(za)f(z)Mza0|(z-a)f(z)| \le M|z-a| \to 0zaz \to a)です。

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

cn12π2πρmaxza=ρf(z)ρn+1=ρnmaxza=ρf(z)=ρn1η(ρ)|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)

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

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

Proposition 3.4極の特徴づけ

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

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

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

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

Proof(Proposition 3.4)

(1) (i) \Rightarrow (iii)mm 位の極なら D(a,R)D'(a,R)f(z)=nmcn(za)nf(z) = \sum_{n \ge -m}c_n(z-a)^n かつ cm0c_{-m} \ne 0 です。両辺に (za)m(z-a)^m を掛けると

(za)mf(z)=nmcn(za)n+m=k0ckm(za)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)=cm0g(a) = c_{-m} \ne 0 であり、0<za<R0 < |z-a| < Rf(z)=g(z)/(za)mf(z) = g(z)/(z-a)^m です。

(iii) \Rightarrow (ii)(za)mf(z)=g(z)g(a)0(z-a)^m f(z) = g(z) \to g(a) \ne 0zaz \to agg の連続性)です。

(ii) \Rightarrow (i)h(z):=(za)mf(z)h(z) := (z-a)^m f(z) とおくと、hhD(a,R)D'(a,R) で正則で zaz \to a のとき有限の極限 L0L \ne 0 を持つので、aa の近くで有界です。Theorem 3.3 の (iii) \Rightarrow (i),(ii) より aahh の除去可能特異点で、h(a)=Lh(a) = L と定めれば hhD(a,R)D(a,R) で正則になります。hh をテイラー展開して h(z)=k0bk(za)kh(z) = \sum_{k \ge 0}b_k(z-a)^kb0=L0b_0 = L \ne 0)とすると、0<za<R0 < |z-a| < R

f(z)=h(z)(za)m=k0bk(za)kmf(z) = \frac{h(z)}{(z-a)^m} = \sum_{k \ge 0} b_k (z-a)^{k-m}

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

(2) 十分性aamm 位の極なら (iii) より f=g/(za)mf = g/(z-a)^mg(a)0g(a) \ne 0 なので、f(z)=g(z)zamg(a)=|f(z)| = |g(z)|\,|z-a|^{-m} \to |g(a)| \cdot \infty = \infty です。

必要性f(z)|f(z)| \to \infty とします。ある r>0r > 0 があって 0<za<r0 < |z-a| < rf(z)1|f(z)| \ge 1、特に f(z)0f(z) \ne 0 です。そこで u:=1/fu := 1/f とおくと uuD(a,r)D'(a,r) で正則、u1|u| \le 1 と有界なので、Theorem 3.3 より aauu の除去可能特異点で、uuD(a,r)D(a,r) で正則に延長されます。f|f| \to \infty より u(a)=0u(a) = 0 です。uuD(a,r)D'(a,r)00 にならないので u≢0u \not\equiv 0 であり、Lemma 7.2[正則関数の強力な性質] で見たとおり、恒等的に 00 でない正則関数の零点は有限の位数を持ちます。その位数を mm とすれば u(z)=(za)mv(z)u(z) = (z-a)^m v(z)vvD(a,r)D(a,r') で正則、v(a)0v(a) \ne 0 と書けます。よって 0<za<r0 < |z-a| < r'f(z)=1/v(z)(za)mf(z) = \dfrac{1/v(z)}{(z-a)^m} となり、1/v1/vaa の近くで正則(v(a)0v(a)\ne 0 と連続性より近傍で v0v \ne 0)かつ aa00 でないので、(1)(iii) から aamm 位の極です。

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

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

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

Proof(Theorem 3.5)

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

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

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

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

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

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

Example 3.6特異点の分類の実例

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

(b) f(z)=1z2(z1)f(z) = \dfrac{1}{z^2(z-1)}z=0z = 0 では z2f(z)=1z110z^2 f(z) = \frac{1}{z-1} \to -1 \ne 0 なので Proposition 3.4 より 22 位の極、z=1z = 1 では (z1)f(z)=1z210(z-1)f(z) = \frac{1}{z^2} \to 1 \ne 0 なので単純極です。

(c) f(z)=1sinzf(z) = \dfrac{1}{\sin z}z=kπz = k\pikZk \in \mathbb{Z})。sin(kπ)=0\sin(k\pi) = 0 かつ ddzsinzz=kπ=cos(kπ)=(1)k0\frac{d}{dz}\sin z\big|_{z=k\pi} = \cos(k\pi) = (-1)^k \ne 0 なので kπk\pisinz\sin z11 位の零点、したがって

(zkπ)f(z)=zkπsinzsin(kπ)1cos(kπ)=(1)k0(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 = 0ew=n0wn/n!e^{w} = \sum_{n\ge0}w^n/n!w=1/zw = 1/z を代入して

e1/z=n01n!zn=1+1z+12!z2+(z0)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<z0 < |z| で収束します。負べきの項が無限個あるので 00 は真性特異点です。実際、値の散らばり方を直接確かめられます。w0w \ne 0 を任意に取り w=weiθw = |w|e^{i\theta} と書くと、zk=1logw+i(θ+2πk)z_k = \dfrac{1}{\log|w| + i(\theta + 2\pi k)}kNk \in \mathbb{N})は kk \to \infty00 に収束し、e1/zk=we^{1/z_k} = w です。つまり 00 のどんな近傍でも、ff00 以外のすべての複素数を無限回取ります。Theorem 3.5 の「稠密」よりずっと強い状況です。

Remark 3.7

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

Remark 3.8

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

Definition 4.1留数

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

Res(f,a):=c1\operatorname{Res}(f, a) := c_{-1}

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

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

でもある。

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

Proposition 4.2極における留数の計算法

aaff の孤立特異点とする。

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

Res(f,a)=1(m1)!limzadm1dzm1[(za)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)=limza(za)f(z)\operatorname{Res}(f,a) = \lim_{z \to a}(z-a)f(z)

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

Res(f,a)=g(a)h(a).\operatorname{Res}(f,a) = \frac{g(a)}{h'(a)}.
Proof(Proposition 4.2)

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

φ(z)=(za)mf(z)=k0ckm(za)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) で正則関数を定めます(Proposition 3.4 の証明と同じ議論です)。この級数の (za)m1(z-a)^{m-1} の係数は k=m1k = m-1 の項、すなわち c1c_{-1} です。一方、正則関数のテイラー係数の公式(Theorem 5.2[正則関数の強力な性質])から (za)m1(z-a)^{m-1} の係数は φ(m1)(a)/(m1)!\varphi^{(m-1)}(a)/(m-1)! です。両者を等しいとおけば結論を得ます。φ\varphiaa で正則なので φ(m1)\varphi^{(m-1)}aa で連続で、φ(m1)(a)=limzaφ(m1)(z)\varphi^{(m-1)}(a) = \lim_{z\to a}\varphi^{(m-1)}(z) と極限で書けます。

(2):(1) で m=1m = 1 とすると 0!=10! = 100 階微分は恒等写像なので、Res(f,a)=limza(za)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 より aahh11 位の零点で、aa の近傍で h(z)=(za)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)zaf(z) = \dfrac{g(z)/h_1(z)}{z-a} であり、g/h1g/h_1aa で正則、aa での値は g(a)/h(a)0g(a)/h'(a) \ne 0 ですから、Proposition 3.4 (1)(iii) より aa は単純極です。(2) を使って

Res(f,a)=limza(za)f(z)=limzag(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)}

を得ます。

Example 4.3留数の計算例

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

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

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

Res(f,2)=limz2ddz[zz1]=limz2(z1)z(z1)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 以上大きい有理関数では、後の Lemma 6.1 (i) と Theorem 5.1 を大きな円周に適用すると留数の総和が 00 と分かります。検算に使えます。

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

Res(f,0)=12!limz0d2dz2[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=z3+z2+12z1+16+e^z/z^3 = z^{-3} + z^{-2} + \frac{1}{2}z^{-1} + \frac{1}{6} + \cdotsz1z^{-1} の係数と一致します。

(c) f(z)=cotz=coszsinzf(z) = \cot z = \dfrac{\cos z}{\sin z}z=kπz = k\pig=cosg = \cosh=sinh = \sin とすると g(kπ)=(1)k0g(k\pi) = (-1)^k \ne 0h(kπ)=0h(k\pi) = 0h(kπ)=cos(kπ)=(1)k0h'(k\pi) = \cos(k\pi) = (-1)^k \ne 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.

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

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

z2e1/z=z2n01n!zn=n01n!z2nz^2 e^{1/z} = z^2\sum_{n\ge0}\frac{1}{n!}z^{-n} = \sum_{n \ge 0}\frac{1}{n!}z^{2-n}

で、z1z^{-1} が現れるのは 2n=12-n = -1、すなわち n=3n = 3 のときですから Res(f,0)=13!=16\operatorname{Res}(f,0) = \frac{1}{3!} = \frac{1}{6} です。留数は極の専売特許ではありません。

Theorem 5.1留数定理

DCD \subset \mathbb{C} を単連結領域とする。CCDD に含まれる区分的に C1C^1 級の単純閉曲線とし、正の向き(内部を左に見る向き)を与える。CC の内部(ジョルダンの曲線定理でいう有界成分)を intC\operatorname{int}C と書く。相異なる有限個の点 a1,,anintCa_1,\dots,a_n \in \operatorname{int}C を取り、ffD{a1,,an}D \setminus \{a_1,\dots,a_n\} で正則とする。このとき

Cf(z)dz=2πij=1nRes(f,aj).\oint_C f(z)\,dz = 2\pi i \sum_{j=1}^{n}\operatorname{Res}(f,a_j).
Proof(Theorem 5.1)

DD は単連結なので intCD\operatorname{int}C \subset D です(単連結領域内の閉曲線の内部は領域に含まれます)。したがって各 aja_jff の孤立特異点で、留数が定義されています。

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

ε<12minjkajak,ε<minjdist(aj,(intC))\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 は有限個で、intC\operatorname{int}C は開集合なので、このような ε\varepsilon が存在します)。γj\gamma_j を正の向きの円周 zaj=ε|z - a_j| = \varepsilon とすると、閉円板 D(aj,ε)\overline{D(a_j,\varepsilon)} は互いに交わらず、すべて intC\operatorname{int}C に含まれます。

積分路を分解するff は「CC の内側かつすべての γj\gamma_j の外側」の閉領域を含む開集合 D{a1,,an}D \setminus \{a_1,\dots,a_n\} で正則です。§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 を固定します。ffD(aj,Rj)D'(a_j, R_j) でローラン展開 f(z)=k=ck(j)(zaj)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 上で一様収束します。一様収束する連続関数列は長さ有限の曲線上で項別積分できる(Lemma 2.2[正則関数の強力な性質])ので

γjf(z)dz=k=ck(j)γj(zaj)kdz.\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(zaj)kdz=02πεkeikθiεeiθdθ=iεk+102π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 .

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

γjf(z)dz=2πic1(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)

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

Remark 5.2

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

ここから先は応用です。実積分を計算するときの手順はいつも同じで、次の三つに分かれます。

実積分の型使う輪郭余分な部分を消す道具
P(x)/Q(x)dx\int_{-\infty}^{\infty}P(x)/Q(x)\,dxdegQdegP+2\deg Q \ge \deg P + 2上半平面の半円Lemma 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上半平面の半円Lemma 6.1 (ii)(ジョルダンの補題)
実軸上に単純極があるとき極を小円でくぼませた輪郭Lemma 8.3 と主値
0xa1g(x)dx\int_0^{\infty}x^{a-1}g(x)\,dxaa は非整数)鍵穴(Appendix)分枝の指定と大小円の評価
a₁a₂−RROΓR実軸上の線分
有理関数の広義積分に使う輪郭。実軸上の線分と上半平面の半円弧をつなぎ、上半平面にある極(×)だけを内部に取り込みます。半径 R を大きくすると弧の寄与が消えます。

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

Lemma 6.1大円弧の寄与が消える条件

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

(i) limRRM(R)=0\lim_{R\to\infty} R\,M(R) = 0 ならば limRΓ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λzdzπλM(R)\left|\int_{\Gamma_R}f(z)e^{i\lambda z}\,dz\right| \le \frac{\pi}{\lambda}M(R)

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

Proof(Lemma 6.1)

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

ΓRf(z)dzπRM(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=Rdθ|dz| = R\,d\theta に注意して

ΓRf(z)eiλzdz0πM(R)eλRsinθRdθ=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) = 0u(θ)=cosθθsinθcosθ=θsinθ<0u'(\theta) = \cos\theta - \theta\sin\theta - \cos\theta = -\theta\sin\theta < 00<θπ/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θe2λRθ/πe^{-\lambda R\sin\theta} \le e^{-2\lambda R\theta/\pi} なので

0π/2eλRsinθdθ0π/2e2λRθ/πdθ=π2λR(1eλ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) を得ます。

Proposition 6.2有理関数の広義積分

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

  • degQdegP+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πiIma>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 にわたる。

Proof(Proposition 6.2)

p=degPp = \deg Pq=degQq = \deg Q、最高次係数を α,β\alpha, \beta とします。z|z|\to\infty のとき zqpP(z)/Q(z)α/βz^{\,q-p}P(z)/Q(z) \to \alpha/\beta なので、ある R11R_1 \ge 1K>0K > 0 があって

zR1    f(z)KzpqKz2|z| \ge R_1 \implies |f(z)| \le K|z|^{\,p-q} \le K|z|^{-2}

です(最後の不等号に qp2q - p \ge 2zR11|z| \ge R_1 \ge 1 を使いました)。したがって xR1f(x)dx2R1Kx2dx<\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 にならない)なので、広義積分は絶対収束します。また RR1R \ge R_1 に対し RM(R)K/R0R\,M(R) \le K/R \to 0 なので Lemma 6.1 (i) の仮定が満たされます。

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

RRf(x)dx+ΓRf(z)dz=2πiIma>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 によりません。RR \to \infty とすると、左辺第 1 項は絶対収束性より f(x)dx\int_{-\infty}^{\infty}f(x)dx に収束し、第 2 項は Lemma 6.1 (i) より 00 に収束します。

Example 6.3分母が 4 次の場合

dxx4+1\displaystyle\int_{-\infty}^{\infty}\frac{dx}{x^4+1} を計算します。P=1P = 1Q=z4+1Q = z^4+1degQdegP=42\deg Q - \deg P = 4 \ge 2Q(x)=x4+11>0Q(x) = x^4 + 1 \ge 1 > 0 より実軸上に零点はないので Proposition 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/43π/43\pi/4 の二つです。いずれも QQ の単純零点(Q=4z3Q' = 4z^300 にならない)なので、Proposition 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^3aa に直しました。この書き換えで計算がぐっと楽になります)。したがって

Res=14(eiπ/4+e3iπ/4)=14(22(1+i)+22(1+i))=142i\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=π22.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 .

被積分関数は正なので答えが正の実数になったのは整合的です。また x1|x| \ge 1 では 1x4+11x2+1\frac{1}{x^4+1} \le \frac{1}{x^2+1}x1|x| \le 1 では 1x4+11\frac{1}{x^4+1} \le 1 なので、積分値は 2+x1dxx2+1=2+π23.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+z12,sinθ=zz12i,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)izdz=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

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

Proposition 7.1a + b cosθ 型の積分

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

02πdθa+bcosθ=2πa2b2.\int_0^{2\pi}\frac{d\theta}{a + b\cos\theta} = \frac{2\pi}{\sqrt{a^2-b^2}} .
Proof(Proposition 7.1)

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

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

I=z=11a+b2(z+1z)dziz=z=12dzi(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±a2b2bz_{\pm} = \frac{-a \pm \sqrt{a^2-b^2}}{b}

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

z+<1|z_+| < 1 を確かめます。z+=a+a2b2bz_+ = \dfrac{-a+\sqrt{a^2-b^2}}{b}a2b2<a\sqrt{a^2-b^2} < a より負で、z+=aa2b2b|z_+| = \dfrac{a - \sqrt{a^2-b^2}}{b} です。z+<1    ab<a2b2|z_+| < 1 \iff a - b < \sqrt{a^2-b^2} であり、両辺とも正(a>ba>b)なので 2 乗して  (ab)2<(ab)(a+b)    ab<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)} とすると、Proposition 4.2 (2) より

Res(g,z+)=limzz+2(zz+)ib(zz+)(zz)=2ib(z+z)=2ib2a2b2b=1ia2b2.\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}} .

Theorem 5.1 より I=2πi1ia2b2=2πa2b2I = 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π2516=2π3\displaystyle\int_0^{2\pi}\frac{d\theta}{5+4\cos\theta} = \frac{2\pi}{\sqrt{25-16}} = \frac{2\pi}{3} です。次は分子に cos2θ\cos 2\theta が入る場合で、00 に多重極が現れます。

Example 7.2分子に cos 2θ がある積分

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

I=z=1(z2+z2)/2(2z2+5z2)/zdziz=z=1z4+12z2i(2z2+5z2)dz=i2z=1z4+1z2(2z1)(z2)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+5z2=(2z1)(z2)-2z^2+5z-2 = -(2z-1)(z-2)12i=i2\frac{1}{-2i} = \frac{i}{2} を使いました)。被積分関数を gg とすると、単位円の内部の特異点は z=0z=022 位の極)と z=1/2z=1/2(単純極)で、z=2z=2 は外部です。

z=1/2z = 1/2Proposition 4.2 (2) より

Res(g,12)=z4+1z22(z2)z=1/2=116+1142(122)=17/163/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 = 0m=2m=2 として φ(z)=z2g(z)=z4+12z25z+2\varphi(z) = z^2 g(z) = \dfrac{z^4+1}{2z^2-5z+2} とおくと

φ(z)=4z3(2z25z+2)(z4+1)(4z5)(2z25z+2)2,φ(0)=01(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 です。念のためローラン展開でも確かめます。125z+2z2=12115z2z22=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 となり、z1z^{-1} の係数は 5/45/4 で一致します。

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

I=i22π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λxdx\int_{-\infty}^{\infty}f(x)\cos\lambda x\,dx のような積分では、cosλz\cos\lambda z をそのまま複素平面に延長してはいけません。上半平面では cos(iy)=coshy\cos(iy) = \cosh yyy \to \infty で指数的に発散するからです。代わりに eiλze^{i\lambda z} を使います。これは上半平面で eiλz=eλImz1|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 なので、最後に実部と虚部を取れば目的の積分が得られます。

Example 8.1ローレンツ型のフーリエ変換

cosxx2+1dx\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}Example 6.3 と同じ半円輪郭 CRC_RR>1R>1)で積分します。FF の特異点は z=±iz = \pm i の単純極で、上半平面にあるのは z=iz = i だけです。Proposition 4.2 (3) を g(z)=eizg(z) = e^{iz}h(z)=z2+1h(z)=z^2+1h(z)=2zh'(z) = 2z に適用して

Res(F,i)=eii2i=e12i.\operatorname{Res}(F, i) = \frac{e^{i\cdot i}}{2i} = \frac{e^{-1}}{2i}.

Theorem 5.1 より RReixx2+1dx+ΓRF(z)dz=2πie12i=π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} です。

円弧の評価には Lemma 6.1 (ii) を使います。z=R>1|z| = R > 1 のとき f(z)=1z2+11R21|f(z)| = \dfrac{1}{|z^2+1|} \le \dfrac{1}{R^2-1}(三角不等式 z2+1z21|z^2+1| \ge |z|^2 - 1)なので M(R)1R210M(R) \le \dfrac{1}{R^2-1}\to0 であり、ΓRFπ11R210\left|\int_{\Gamma_R}F\right| \le \dfrac{\pi}{1}\cdot\dfrac{1}{R^2-1}\to0 です。よって

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

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

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

Definition 8.2コーシーの主値

cRc \in \mathbb{R} とし、ffR{c}\mathbb{R}\setminus\{c\} で連続とする。極限

p.v.f(x)dx:=limRε+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 留数」になる、と憶えると自然です。

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

ccff の単純極とし、ρ=Res(f,c)\rho = \operatorname{Res}(f,c) とおく。実数 α,β\alpha, \beta に対し、γε\gamma_\varepsilonz=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 なら時計回りである。)

Proof(Lemma 8.3)

cc は単純極なので、ローラン展開は f(z)=ρzc+n0cn(zc)nf(z) = \dfrac{\rho}{z-c} + \sum_{n \ge 0}c_n(z-c)^n の形です。第 2 項を g(z)g(z) とおくと、ggD(c,R)D(c,R) で正則なべき級数の和ですから、ある r>0r>0M>0M>0 に対し zcr|z-c| \le rg(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 より

γερzcdz=ραβ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βαεM0(ε+0).\left|\int_{\gamma_\varepsilon}g(z)\,dz\right| \le |\beta-\alpha|\,\varepsilon\, M \longrightarrow 0 \qquad (\varepsilon \to +0).

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

Example 8.4ディリクレ積分

0sinxxdx=π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:y0}\{iy : y \le 0\} を除いた単連結領域 DD に含まれ、ffDD で正則です(00DD に属しません)。よって Theorem 5.1 を特異点 00 個の場合に使って(これはコーシーの積分定理そのものです)

εReixxdx+ΓRfdz+Rεeixxdx+γεfdz=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εeitt(dt)=εReittdt\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 項と合わせて

εReixeixxdx=2iεRsinxxdx\int_\varepsilon^R\frac{e^{ix}-e^{-ix}}{x}\,dx = 2i\int_\varepsilon^R\frac{\sin x}{x}\,dx

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

以上を合わせて ε+0\varepsilon \to +0RR \to \infty とすると

2i0sinxxdxiπ=0,すなわち0sinxxdx=π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.cosxxdx=0\mathrm{p.v.}\displaystyle\int_{-\infty}^{\infty}\frac{\cos x}{x}dx = 0(奇関数なので当然です)、虚部から上の結果が出ます。

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

Exercise 9.1

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

Solution

特異点は分母の零点 z=0z=022 位)と z=1z=111 位)です。eze^z はこれらの点で 00 にならないので、位数は分母のとおりです。実際 z2f(z)=ezz110z^2f(z) = \frac{e^z}{z-1} \to -1 \ne 0(z1)f(z)=ezz2e0(z-1)f(z) = \frac{e^z}{z^2}\to e \ne 0 なので Proposition 3.4 より 0022 位の極、11 は単純極です。

z=1z=1Proposition 4.2 (2) より Res(f,1)=limz1ezz2=e\operatorname{Res}(f,1) = \lim_{z\to1}\frac{e^z}{z^2} = e

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

φ(z)=ez(z1)ez(z1)2=ez(z2)(z1)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 .

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

Exercise 9.2標準

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

Solution

被積分関数は偶関数なので 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 にならないので Proposition 6.2 が使えます。上半平面の極は単純極 z=iz=iz=2iz=2i です。

Res(f,i)=z2(z+i)(z2+4)z=i=12i3=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 を使いました。)和は i6i3=i6\frac{i}{6}-\frac{i}{3} = -\frac{i}{6} なので

f(x)dx=2πi(i6)=π3,0f(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=14A+B=04A+B=0、すなわち A=13A = -\frac13B=43B=\frac43 です。0dxx2+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}

で一致します。

Exercise 9.3標準

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

Solution

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

I=z=14z2(z2+4z+1)2dziz=4iz=1zdz(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+30.268z_0 = -2+\sqrt3 \approx -0.268z1=233.73z_1 = -2-\sqrt3 \approx -3.73 で、単位円の内部にあるのは z0z_0 だけです(z0z1=1z_0z_1 = 1 なので一方が内側、他方が外側)。分母は (zz0)2(zz1)2(z-z_0)^2(z-z_1)^2 なので z0z_022 位の極です。m=2m=2 として Proposition 4.2 (1) を使います。φ(z)=z(zz1)2\varphi(z) = \dfrac{z}{(z-z_1)^2} とすると

φ(z)=(zz1)2z2(zz1)(zz1)4=(zz1)2z(zz1)3=zz1(zz1)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}.

z0z1=23z_0 - z_1 = 2\sqrt3z0z1=(2+3)(23)=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}.

Theorem 5.1 より

I=4i2πi163=8π63=4π33=43π92.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 .

検算:Proposition 7.1 の式 02πdθa+bcosθ=2π(a2b2)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(a2b2)3/2-2\pi a(a^2-b^2)^{-3/2} です。a=2a=2b=1b=1 とすると 4π33/2=4π33\frac{4\pi}{3^{3/2}} = \frac{4\pi}{3\sqrt3} で一致します。

Exercise 9.4

a>0a > 0 とするとき 0xsinxx2+a2dx=π2ea\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} には Lemma 6.1 (i) が使えないことに注意してください。

Solution

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 だけで、Proposition 4.2 (3) を g(z)=zeizg(z) = ze^{iz}h(z)=z2+a2h(z) = z^2+a^2h(z)=2zh'(z)=2z に適用して

Res(F,ia)=iaeiia2ia=ea2.\operatorname{Res}(F, ia) = \frac{ia\,e^{i\cdot ia}}{2ia} = \frac{e^{-a}}{2}.

Theorem 5.1 より RRxeixx2+a2dx+ΓRFdz=2πiea2=iπea\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)RR2a2|f(z)| \le \dfrac{R}{R^2-a^2} なので M(R)0M(R) \to 0 ですが、RM(R)10R\,M(R) \to 1 \ne 0 なので Lemma 6.1 (i) は使えません。指数因子を持つ (ii)(ジョルダンの補題)なら λ=1\lambda = 1 として

ΓRF(z)dzπM(R)πRR2a20\left|\int_{\Gamma_R}F(z)\,dz\right| \le \pi M(R) \le \frac{\pi R}{R^2-a^2}\longrightarrow 0

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

xsinxx2+a2dx=πea\int_{-\infty}^{\infty}\frac{x\sin x}{x^2+a^2}\,dx = \pi e^{-a}

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

  • 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 に対して

0xa11+xdx=πsinπa\int_0^{\infty}\frac{x^{a-1}}{1+x}\,dx = \frac{\pi}{\sin \pi a}

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

分枝を決めるLogz=logz+iargz\operatorname{Log}z = \log|z| + i\arg zargz(0,2π)\arg z \in (0, 2\pi) と定め、za1:=e(a1)Logzz^{a-1} := e^{(a-1)\operatorname{Log}z} とします。これは C\mathbb{C} から非負実軸を除いた領域で正則です。f(z)=za11+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 だけです。

円周上の評価。大円上では za1=Ra1|z^{a-1}| = R^{a-1}1+zR1|1+z| \ge R-1 なので、ML 不等式より寄与は 2πRRa1R12πRa102\pi R\cdot\dfrac{R^{a-1}}{R-1} \sim 2\pi R^{a-1} \to 0a<1a<1)です。小円上では 2πεεa11ε2πεa02\pi\varepsilon\cdot\dfrac{\varepsilon^{a-1}}{1-\varepsilon} \sim 2\pi\varepsilon^{a}\to 0a>0a>0)です。0<a<10 < a < 1 という仮定は、この二つの評価を同時に成立させるためにあります。

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

留数と結論Res(f,1)=(1)a1=e(a1)iπ=eiπa\operatorname{Res}(f,-1) = (-1)^{a-1} = e^{(a-1)i\pi} = -e^{i\pi a} なので、Theorem 5.1 より (1e2π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πieiπa1e2πia=2πieiπaeiπa=2πi2isinπ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 と置換すれば

0dxx(1+x)=02tdtt(1+t2)=2[arctant]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)Γ(1a)=π/sinπa\Gamma(a)\Gamma(1-a) = \pi/\sin\pi a の複素解析による証明にもなっています。

Report an error in this article ・Operated by: Mugen Giken LLCPricingTermsLegal notice

© 2026 夢現技研合同会社 ・Feeding the text to an LLM is welcome. Code samples are MIT licensed.