三次元対象物の複素積分表現(事例紹介)

本記事では, 複素解析の教科書ではあまり見られない,三次元対象物の複素積分による表現をいくつかの事例で紹介します. 従来と少し異なる視点を提供することにより, 複素解析を学ばれる方々の刺激になることを期待しています.

1. 基本的な式

ここでは, コーシーの積分公式を含む複素解析の基本的な式を取り上げる. 詳しい定義や導出等は複素解析の教科書をご参照願いたい.

さて, D は複素平面上の単連結領域(穴が開いていない領域)とし, C はそれを囲うある長さを持つ単純閉曲線(自身と交わらない閉じた曲線)とする. D の任意の一点 a において, 以下のコーシー・ポンペイウの公式(Cauchy-Pompeiu Formula)が成り立つ.

f(a) = \cfrac{1}{2\pi i} \displaystyle\oint_{C=\partial D} \cfrac{f(z)}{z-a} dz + \cfrac{1}{2\pi i} \displaystyle\int \displaystyle\int_{D} \cfrac{ \frac{\partial f(z)}{\partial \bar{z}}}{z-a} dz d \bar{z} \tag{1.1}

ここで, \bar{z} は, 複素数 z の複素共役(complex conjugate)である. また,

dz d \bar{z} = -2 i d D \tag{1.2}

であることから, 式(1.1)は二項目を書き変えて,

f(a) = \cfrac{1}{2\pi i} \displaystyle\oint_{C=\partial D} \cfrac{f(z)}{z-a} dz - \cfrac{1}{\pi} \displaystyle\int_{D} \cfrac{ \frac{\partial f(z)}{\partial \bar{z}}}{z-a} dD \tag{1.3}

とも表せる.

さて, f(z)D 上の正則関数(holomorphic function)であるとき, \cfrac{\partial f(z)}{\partial \bar{z}}=0 であるので, 式(1.1)あるいは式(1.3)は,

f(a) = \cfrac{1}{2\pi i} \displaystyle\oint_{C} \cfrac{f(z)}{z-a} dz \tag{1.4}

となる. これがコーシーの積分公式(Cauchy Integral Formula)と呼ばれるものである. また, 式(1.4)の特別な場合 f(z) \equiv 1 として, いわゆるコーシーの積分定理(Cauchy Integral Theorem)が成り立つ.

1 = \cfrac{1}{2\pi i} \displaystyle\oint_{C} \cfrac{1}{z-a} dz \tag{1.5}

そして, 式(1.4)と式(1.5)から次が成り立つ.

f(a) = \cfrac{\displaystyle\oint_{C} \cfrac{f(z)}{z-a} dz }{\displaystyle\oint_{C} \cfrac{1}{z-a}dz}  \tag{1.6}

なお, 式(1.1)において, f(z)=\pi (z-a)  \bar{z} (これは正則関数ではない)とおけば,

\cfrac{1}{2i} \displaystyle\oint_{C=\partial D} \bar{z} dz = \cfrac{1}{2i} \displaystyle\int \displaystyle\int_{D} d \bar{z} dz = \displaystyle\int_{D} dD = D \tag{1.7}

という D に関する基本的な関係式が得られる.

2. コーシーの積分公式の幾何学的解釈

三次元対象物の複素積分による表現に入る前に, 複素積分自体の幾何学的意味を見るために, ある変数変換により式(1.6)を書き換え, コーシーの積分公式の幾何学的な解釈を行ってみよう.

2.1 変数変換

以下の変数変換を考える.

z^{\prime} = \log(z-a) \tag{2.1}

ここで, \log(z) は自然対数である. 複素関数の対数は一般に多価性があるが, 本稿では1価に制限されているものとする.

ここで, z-a = re^{i \theta} , 0 \leq \theta \leq 2\pi とすると, この変数変換に伴い z^{\prime} = \log(r)+i \theta , 0 \leq \theta \leq 2\pi になり, 単純閉曲線 C は, 開いた曲線 C^{\prime} になる.

2.2 幾何学的解釈

式(1.6)は, \cfrac{d}{dz} \log(z) = \cfrac{1}{z} 及び変数変換(2.1)を用いると, 以下のように書き換えられる.

f(a) &= \cfrac{ \displaystyle\oint_{C} f(z) d( \log(z-a) ) }{ \displaystyle\oint_{C} d( \log(z-a) ) } \tag{2.2}\\  &= \cfrac{ \displaystyle\int_{C^{\prime}} f( e^{z^{\prime}} + a ) d z^{\prime} }{ \displaystyle\int_{C^{\prime}} d z^{\prime} } \tag{2.3}

式(2.3)によれば, f(a) は, (開いた)曲線 C^{\prime} に沿って z^{\prime} が動いた時の関数 g(z^{\prime}) \equiv f(e^{z^{\prime}}+ a) の平均値(あるいは重心)を与えていると解釈できる. それゆえ, 式(2.3)は, 平均値の定理(mean-value theorem)と呼ばれる.

2.3 解釈の整合性

実は, 上記の議論で,

\cfrac{1}{2\pi i} \displaystyle\oint_{C=\partial D} \cfrac{f(z)}{z-a} dz \tag{2.4}

という積分は, 変数変換(2.1)を行わなくてもそのまま, C 上を \cfrac{1}{z-a} という関数について z で積分するとき, f(z) という重みを与えて平均化している, とも解釈でき, しかもこの解釈自体は f(z) が正則か否かには関係ない. そのため, たとえば, 式(1.1)の右辺第一項にもこの解釈を適用可能である.

さて, 平均値(2.4)は, 平均値(2.4)自体を関数 \cfrac{1}{z-a}C にそって z で積分する合計値と一致するはずである. すなわち,

\displaystyle \oint_{C= \partial D} \cfrac{1}{z-a} \left( \cfrac{1}{2\pi i} \displaystyle \oint_{C=\partial D} \cfrac{f(w)}{w-a} dw \right) dz = \displaystyle \oint_{C=\partial D} \cfrac{f(z)}{z-a} dz \tag{2.5}

実際, ここで, 左辺の括弧内に式(1.1)を用いれば,

\displaystyle \oint_{C= \partial D} \cfrac{1}{z-a} \left( f(a) -  \cfrac{1}{2\pi i} \displaystyle\int \displaystyle\int_{D} \cfrac{ \frac{\partial f(w)}{\partial \bar{w}}}{w-a} dw d \bar{w} \right) dz = \displaystyle \oint_{C=\partial D} \cfrac{f(z)}{z-a} dz \tag{2.6}

であり, 左辺は,

2\pi i  f(a) -  \displaystyle\int \displaystyle\int_{D} \cfrac{ \frac{\partial f(w)}{\partial \bar{w}}}{w-a} dw d \bar{w} \tag{2.7}

であることから, 両辺を 2\pi i で割れば, コーシー・ポンペイウの公式が再現され, この公式と整合していることが確認される.

3. 三次元的な対象物を複素積分で表現する

筆者は, 中学の終わりごろから, 独学で微分積分学を学び, ついでベクトル解析を学び, n 次元球などの一般次元の空間の対象物を取り扱えるようになったあとで, 複素解析を学び始めた途端, 空間が突如二次元の世界に限定されてしまったような印象を持った. たとえば, せっかく習得したストークスの定理(Stokes' Theorem)などはどこへ行ってしまったのか, と思ったりした.

しかし, もちろん, 複素解析には本来そのような限定はない. 三次元以上の空間の対象と結び付けることが可能である. ここでは, 簡単な事例を挙げてそのことを示したい.

3.1 立体の体積

式(1.2)(または, 式(1.7))から,

D = \frac{1}{2i} \displaystyle\int \displaystyle\int_{D} d \bar{z} dz \tag{3.1}

である. ここで, D が時間的に変化する(つまり C=\partial D が時間的に変化する)としよう. すなわち, 各時点 t での複素平面というものを考えることにする.

D(t) = \frac{1}{2i} \displaystyle\int \displaystyle\int_{D(t)} d \bar{z} dz \tag{3.2}

立体の体積を複素積分で表現するために, 立体を一方向に平面でスライスしていく. このとき各平面が各時点の複素平面であるようにする. すると, 時刻 t=0 から 時刻 t=T までかけて D(t) は点から立体の断面になり, 立体の体積 V は, 以下のように表せる.

V &= \displaystyle\int_{V} dV \\&= \displaystyle\int_{0}^{T} D(t) dt \\&= \frac{1}{2i} \displaystyle\int_{0}^{T} \displaystyle\int \displaystyle\int_{D(t)} d \bar{z} dz dt\tag{3.3}

3.2 球の体積

ここで, 具体的な例として, 3次元の球を対象に考えてみよう. 球をある直径に沿って刻々とスライスしていく断面 D(t) を考える.時刻 t=0 から 時刻 t=T までかけて D(t) は点から半径 T の円盤になり, 時刻 t=T から 時刻 t=2T までかけて D(t) は再び点になるとする. 時刻 t のときの D(t) は,

D(t) = \left\{ z \,\, | \,\, |z| \leq \sqrt{2Tt-t^{2}} \right\} \tag{3.4}

となり, 時刻 t=0 から 時刻 t=2T まで厚み dt の円盤 D(t) を積分する形で球の体積が求まり,

\frac{4 \pi}{3} T^{3} &=  \displaystyle\int_{0}^{2T} D(t) dt  \\&= \displaystyle\int_{0}^{2T} \pi (\sqrt{2Tt-t^{2}})^2 dt  \\&= \frac{1}{2i} \displaystyle\int_{0}^{2T} \displaystyle\int \displaystyle\int_{D(t)} d \bar{z} dz dt  \tag{3.5}
という関係が得られる.
timeslicedcp.png

ところで, 式(3.5)では, 時刻 t の円盤(つまり2次元球) D(t) を足し上げて三次元球の体積を求めたわけだが, 同様にして三次元球を足し上げることで, 四次元球の体積を求めることができる. 時刻 t のときの三次元球の体積 B_{3}(t) は,

B_{3}(t) = \frac{4 \pi}{3} \bigl( \sqrt{2Tt-t^{2}} \bigr)^{3} \tag{3.6}

であり, 四次元球の体積は,

B_{4}= \displaystyle\int_{0}^{2T} B_{3}(t) dt = \frac{{\pi}^{2} T^{4}}{2} \tag{3.7}

となる. このことを踏まえ, 時刻をもう一つ増やして, 式(3.5)に類似した形で B_{4} について複素積分で表すと,

B_{4} = \frac{1}{2i} \displaystyle\int_{0}^{2T} \displaystyle\int_{0}^{2 \sqrt{2Tt_{2} - {t_{2}}^2}} \displaystyle\int \displaystyle\int_{D(t_{1})} d \bar{z} dz  dt_{1} dt_{2} \tag{3.8}

となる. このようにして, 複素積分を一般次元の球の体積と結び付けられる. なお, ここで, D(t_{1})=B_{2}(t_{1})=\pi (\sqrt{2T t_{1}-{t_{1}}^{2}})^2 である.

3.3 ストークスの定理

3.1項と同様に, 各時点の複素平面を考えることで三次元的な空間を作る. 座標としては, z=x+iyt を使って, 位置ベクトル \bm{r}=(x,iy,t) を考える. すると, 線素は d{\bm r}=(dx,idy,dt) , 面積要素は d \bm{S} = (idydt,dtdx,idxdy) = \frac{1}{2i}((dz-d \bar{z})dt, idt(dz+d \bar{z}), d \bar{z} dz  ) になる. ただし, ここで, dz=dx+idy , d\bar{z}=dx-idy である.

このような複素数を含んだベクトル表示における二つのベクトル \bm{A}= (a_{1}, a_{2}, a_{3}) , \bm{B}= (b_{1}, b_{2}, b_{3}) の内積及び外積を次のように定義することとする.

\bm{A} \cdot \bm{B} \equiv  \overline{a_{1}}b_{1} + \overline{a_{2}}b_{2} + \overline{a_{3}}b_{3} \bm{A} \times \bm{B} \equiv ( \overline{a_{2}}b_{3} -  \overline{a_{3}}b_{2}, \overline{a_{3}}b_{1} -  \overline{a_{1}}b_{3}, \overline{a_{1}}b_{2} -  \overline{a_{2}}b_{1} )

これらはそれぞれ成分が実数の場合の定義を包含している. なお,このとき,ベクトル \bm{A} の大きさ(ノルム)は, 成分が実数の場合と同様に \sqrt{ \bm{A} \cdot \bm{A}} で与えられる.

さて, ベクトル場 \bm{A}= (A_{1}, iA_{2}, A_{3}) に対し, 同三次元空間の単純閉曲線 L とそれを縁とする曲面 S について,

\displaystyle \oint_{L= \partial S} \bm{A} \cdot d \bm{r} &= \displaystyle \oint_{L= \partial S} A_{1} dx + (-i A_{2}) (i dy) +A_{3} dt \\  &= \displaystyle \oint_{L= \partial S} A_{1} dx +  A_{2} dy +A_{3} dt \\ &=  \displaystyle\int \displaystyle\int_{S} \Bigl( \frac{\partial A_{3}}{\partial y}- \frac{\partial A_{2}}{\partial t} \Bigr) dydt +  \Bigl( \frac{\partial A_{1}}{\partial t}- \frac{\partial A_{3}}{\partial x} \Bigr) dtdx + \Bigl( \frac{\partial A_{2}}{\partial x}- \frac{\partial A_{1}}{\partial y} \Bigr) dxdy \\ &=  \displaystyle\int \displaystyle\int_{S} \Bigl( \frac{1}{i}\frac{\partial A_{3}}{\partial y}- \frac{\partial (-i A_{2})}{\partial t} \Bigr) (idy)dt +  \Bigl( \frac{\partial A_{1}}{\partial t}- \frac{\partial A_{3}}{\partial x} \Bigr) dtdx + \Bigl( \frac{\partial (-i A_{2})}{\partial x}- \frac{1}{i}\frac{\partial A_{1}}{\partial y} \Bigr) dx(idy) \\ &=  \displaystyle\int \displaystyle\int_{S} \overline{ \Bigl( i \frac{\partial A_{3}}{\partial y}- \frac{\partial (i A_{2})}{\partial t} \Bigr)} (idy)dt + \overline{ \Bigl( \frac{\partial A_{1}}{\partial t}- \frac{\partial A_{3}}{\partial x} \Bigr)} dtdx + \overline{ \Bigl( \frac{\partial (i A_{2})}{\partial x}- \Bigl( i \frac{\partial A_{1}}{\partial y} \Bigr) \Bigr)} dx(idy) \\ &= \displaystyle\int_{S} \nabla \times \bm{A} \cdot d \bm{S} \tag{3.9}

であり, 実数解析のストークスの定理を利用することで, そのままストークスの定理(Stokes' Theorem)が成り立つ.

\displaystyle \oint_{L= \partial S} \bm{A} \cdot d \bm{r} = \displaystyle\int_{S} \nabla \times \bm{A} \cdot d \bm{S} \tag{3.10}

ただし, ここで, \nabla \equiv \Bigl( \frac{\partial}{\partial x}, \frac{1}{i} \frac{\partial}{\partial y},  \frac{\partial}{\partial t} \Bigr) である.

ガウスの定理(Gauss' Theorem)については,三次元空間のベクトル場 \bm{B}= (B_{1}, iB_{2}, B_{3}) を考えれば, 同三次元空間の単純閉曲面 S とそれを縁とする体積 V について,

\displaystyle \oint_{S= \partial V} \bm{B} \cdot d \bm{S} &= \, \, \, \displaystyle \int \! \! \! \!\displaystyle \int_{S= \partial V }\! \! \! \! \! \! \! \! \! \! \! \! \! \! \! \! \! \! \bigcirc \,\,\,\,\,\,\, B_{1} (idydt) +(- i B_{2}) dtdx +B_{3} (idxdy) \\ &= i \displaystyle \int \! \! \! \! \displaystyle \int_{S= \partial V }\! \! \! \! \! \! \! \! \! \! \! \! \! \! \! \! \! \! \bigcirc \,\,\,\,\,\,\, B_{1} dydt + (-B_{2}) dtdx + B_{3} dxdy \\ &= i \displaystyle \int \displaystyle \int \! \displaystyle \int_{V} \Bigl( \frac{\partial B_{1}}{\partial x} +  \frac{\partial (-B_{2})}{\partial y} + \frac{\partial B_{3}}{\partial z}  \Bigr) dxdydt \\ &=  \displaystyle \int \displaystyle \int \! \displaystyle \int_{V} \Bigl( \frac{\partial B_{1}}{\partial x} - \frac{1}{i} \frac{\partial (iB_{2})}{\partial y} + \frac{\partial B_{3}}{\partial z}  \Bigr) (idxdydt) \\ &=  \displaystyle \int_{V} \nabla \cdot \bm{B} dV  \tag{3.11}

であり, 実数解析のガウスの定理を利用することで, そのままガウスの定理が成り立つ.

\displaystyle \oint_{S= \partial V} \bm{B} \cdot d \bm{S} = \displaystyle\int_{V} \nabla \cdot \bm{B} d V \tag{3.12}

同様にして, ベクトル解析の諸公式を複素積分で表現することができる. ここでは詳しく展開できないが, 当然のことながら, 三次元の流体力学等を複素積分で表現することも可能である.

3.4 パップスの定理

3.3項で導入した 位置ベクトル \bm{r} , 線素 d \bm{r} 及び面積要素 d \bm{S} の表式を用いれば, 幾何学のパップス・ギュルダンの定理(Pappus-Guldinus theorem)(以下, パップスの定理)を複素積分で表現できる.

パップスの定理では, 断面上のすべての点が断面に垂直になるように(すなわち d{\bm r}\cdot d{\bm S}=dr dS となるように)断面 S を動かし, それが掃する体積 VS の重心の動いた道のり L と面積 S の積になる. 3.2項では, 直線方向に時点の異なる複素平面が並んだが, この並び方は回転してもいい. このようなことを利用して, たとえば, 半円盤を直径の周りに回転させて球を作り, その体積から半円盤の重心の位置を求めたり, これを高次化して, 半球を直径断面の周りに回転させて四次元球を作り, その体積から半球の重心の位置を求めたりすることができる. 重心の軌道のパラメータを t^{\prime} とすると, パップスの定理は一般式としては,

V = S \times \displaystyle\int_{L}  \Bigm\| \cfrac{d}{dt^{\prime}} \Bigl( \cfrac{\int_{S} \bm{r} dS }{\int_{S} dS} \Bigr) \Bigm\| d t^{\prime} \tag{3.13}

と表すことができる. ただし, 上で, dr= \sqrt{ d \bm{r} \cdot d \bm{r} } , dS= \sqrt{ d \bm{S} \cdot d \bm{S} } である. (パップスの定理について, 詳しくは本記事末の関連メモをご覧いただきたい.)

3.5 補足

多変数複素解析では, (z_{1}, z_{2}, \cdots ,z_{n}) \in \mathbb{C}^n を用いて, 2n 次元の空間 \mathbb{R}^{2n} 内の体積を扱うことができる.

dV = \Bigl( \frac{1}{2i} \Bigr) ^{n} d \overline{z_{1}} dz_{1} d \bar{z_{2}} dz_{2} \cdots d \overline{z_{n}} dz_{n} \tag{3.14}

4. おわりに

本記事では, 三次元対象物を複素積分で表現する事例をいくつか示しました. いわば直接見える対象物を直接は見えない世界(複素数の世界)に埋め込んでいる恰好になっています. 逆に, 直接は見えない複素数の世界を直接見えるこちら側に持ってこられるならば(理解とは結局そういうことなのかもしれませんが), もっと面白いことが分かってくるかもしれません.

The English version of this article is here .

On Generalizing The Theorem of Pappus is here2.