平成19年度 東大院 物理学専攻 修士 数学 解答

出典と、この解答の位置づけ

本記事が扱う問題は、東京大学大学院理学系研究科物理学専攻の大学院入試問題からの引用です（引用元: https://www.phys.s.u-tokyo.ac.jp/）。

• 掲載している解答・解説は東京大学が公表したものではなく、夢現技研合同会社が独自に作成したものです。誤りが含まれる可能性があります。
• 本記事に問題文は転載していません。問題は上記の専攻公式サイトでご確認ください。
• 本記事の内容について、東京大学および同専攻へのお問い合わせはご遠慮ください。誤りのご指摘はお問い合わせまでお願いします。

大問2題で、2問すべてが必答です（試験時間は60分）。第1問は正規行列とユニタリ対角化という前提を与えたうえで、2次の回転行列を題材に直交性・行列指数関数・複素対角化をたどり、後半でファンデルモンド行列式に移ります。第2問はコーシー–リーマン関係式だけを出発点に、調和性、ストークスの定理による線積分の消滅、コーシーの積分定理、そして $\int e^{-x^2}\cos(2px)\,dx$ の評価まで一本道で誘導する構成です。計算量は軽く、微分の交換可能性の根拠、ストークスの定理を適用できる条件、長方形経路の側辺の寄与が消えることの評価といった論理の詰めが配点の中心になります。

問題	分野	主題
第1問	線形代数	回転行列の直交性・行列指数関数・ファンデルモンド行列式
第2問	複素解析・フーリエ解析	コーシー–リーマン関係式、コーシーの積分定理、ガウス型積分

第1問 回転行列の直交性・行列指数関数・ファンデルモンド行列式

設定

前提として次の事実が与えられます。正方行列 $A$ に対し $U^{-1}AU$ が対角行列となるユニタリ行列 $U$ が存在するための必要十分条件は、$A$ が正規行列であること、すなわち

$$A^{*}A=AA^{*}$$

が成り立つことです。ここで $A^{*}$ は $A$ を転置して複素共役をとった行列（随伴行列）です。実成分のユニタリ行列を直交行列と呼びます。

前半では 2 次の行列

$$B=\begin{pmatrix} \cos\theta & -\sin\theta \\ \sin\theta & \cos\theta \end{pmatrix} \qquad (\theta \in \mathbb{R})$$

を扱い、後半では変数の組 $\{x_1,\dots,x_n\}$ に対する $n$ 次の行列

$$C=\begin{pmatrix} 1 & 1 & \cdots & 1 \\ x_1 & x_2 & \cdots & x_n \\ \vdots & \vdots & & \vdots \\ x_1^{\,n-1} & x_2^{\,n-1} & \cdots & x_n^{\,n-1} \end{pmatrix}, \qquad C_{ij}=x_j^{\,i-1}$$

を扱います。$\lvert C \rvert$ は $C$ の行列式を表します。

設問1

$B$ の成分は $\theta$ が実数なのですべて実数です。したがって複素共役をとる操作は何もせず、$B^{*}=B^{\mathsf{T}}$ が成り立ちます。転置は

$$B^{\mathsf{T}}=\begin{pmatrix} \cos\theta & \sin\theta \\ -\sin\theta & \cos\theta \end{pmatrix}$$

です。積を計算すると

$$\begin{aligned} B^{\mathsf{T}}B&=\begin{pmatrix} \cos\theta & \sin\theta \\ -\sin\theta & \cos\theta \end{pmatrix} \begin{pmatrix} \cos\theta & -\sin\theta \\ \sin\theta & \cos\theta \end{pmatrix}\\[2pt] &=\begin{pmatrix} \cos^2\theta+\sin^2\theta & -\cos\theta\sin\theta+\sin\theta\cos\theta \\ -\sin\theta\cos\theta+\cos\theta\sin\theta & \sin^2\theta+\cos^2\theta \end{pmatrix} =\begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix}=E \end{aligned}$$

となり、同じ計算で

$$BB^{\mathsf{T}}=\begin{pmatrix} \cos^2\theta+\sin^2\theta & \cos\theta\sin\theta-\sin\theta\cos\theta \\ \sin\theta\cos\theta-\cos\theta\sin\theta & \sin^2\theta+\cos^2\theta \end{pmatrix}=E$$

も得られます。使ったのは $\cos^2\theta+\sin^2\theta=1$ だけです。

よって $B^{*}B=BB^{*}=E$、つまり $B$ はユニタリ行列であり、しかも成分が実数なので定義により直交行列です。答えは、$B^{\mathsf{T}}B=BB^{\mathsf{T}}=E$ が実数成分のまま成立するので $B$ は直交行列である、ということになります。

なお $B^{*}B=BB^{*}=E$ は $B$ が正規行列であることも意味しており、設定で述べられた定理から $B$ がユニタリ行列で対角化できることが保証されます。これを実行するのが設問3です。また $\lvert B \rvert=\cos^2\theta+\sin^2\theta=1$ なので $B$ は回転（$SO(2)$ の元）であり、鏡映は含みません。

設問2

$A$ を回転の生成子で書き直してから偶数次と奇数次に分けます。

$$A=\begin{pmatrix} 0 & -\theta \\ \theta & 0 \end{pmatrix}=\theta J, \qquad J\equiv\begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix}$$

とおくと

$$J^2=\begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix}\begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix} =\begin{pmatrix} -1 & 0 \\ 0 & -1 \end{pmatrix}=-E$$

です。これを繰り返せば、$k=0,1,2,\dots$ に対して

$$J^{2k}=(-1)^k E,\qquad J^{2k+1}=(-1)^k J$$

が成り立ちます（$k$ についての帰納法。$J^{2(k+1)}=J^{2k}J^2=(-1)^k E\cdot(-E)=(-1)^{k+1}E$）。

$A^m=\theta^m J^m$ なので、定義式は

$$\exp A=\sum_{m=0}^{\infty}\frac{\theta^m J^m}{m!}$$

となります。この級数は絶対収束します。実際、行列ノルムを劣乗法的なもの（たとえば作用素ノルム）にとれば $\lVert A^m/m!\rVert \le \lVert A\rVert^m/m!$ であり、$\sum_m \lVert A\rVert^m/m!=e^{\lVert A\rVert}<\infty$ です。$2\times 2$ 行列の空間は有限次元なので完備であり、絶対収束する級数は収束します。絶対収束していれば項の並べ替えが自由にできるので、$m$ を偶数 $m=2k$ と奇数 $m=2k+1$ に分けて

$$\begin{aligned} \exp A&=\sum_{k=0}^{\infty}\frac{\theta^{2k}}{(2k)!}J^{2k}+\sum_{k=0}^{\infty}\frac{\theta^{2k+1}}{(2k+1)!}J^{2k+1}\\[2pt] &=\left(\sum_{k=0}^{\infty}\frac{(-1)^k\theta^{2k}}{(2k)!}\right)E +\left(\sum_{k=0}^{\infty}\frac{(-1)^k\theta^{2k+1}}{(2k+1)!}\right)J\\[2pt] &=(\cos\theta)\,E+(\sin\theta)\,J \end{aligned}$$

と書けます。最後の等号は $\cos$ と $\sin$ のテイラー級数（実数 $\theta$ 全体で収束）です。成分に戻すと

$$\exp A=\cos\theta\begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix} +\sin\theta\begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix} =\begin{pmatrix} \cos\theta & -\sin\theta \\ \sin\theta & \cos\theta \end{pmatrix}=B$$

となり、示すべき式が得られました。

検算として $\theta=0$ を入れると $A=O$、$\exp O=E$ で $B=E$ と一致します。また $\theta$ で微分すると $\dfrac{d}{d\theta}\exp(\theta J)=-\sin\theta\,E+\cos\theta\,J=J\exp(\theta J)$ で、行列指数関数が満たすべき微分方程式と整合します。設問3の対角化を使った別の導出は下で述べます。

設問3

固有値は特性方程式から求めます。

$$\lvert B-\lambda E \rvert=(\cos\theta-\lambda)^2+\sin^2\theta=\lambda^2-2\lambda\cos\theta+1=0$$

です。判別式は $4\cos^2\theta-4=-4\sin^2\theta\le 0$ なので、$\theta$ が実数の範囲では根は

$$\lambda_{\pm}=\cos\theta\pm i\sin\theta=e^{\pm i\theta}$$

の2つです（$\sin\theta$ の符号がどちらでも、複号がこの2根を尽くします）。$\sin\theta\neq 0$ のとき固有値は実数ではないので、対角化行列は実行列ではとれません。

固有ベクトルを求めます。$\lambda_{+}=e^{i\theta}$ に対する条件の第1行は $(\cos\theta-e^{i\theta})w_1-\sin\theta\,w_2=-i\sin\theta\,w_1-\sin\theta\,w_2=0$、すなわち $w_2=-iw_1$ です。同様に $\lambda_{-}=e^{-i\theta}$ では $w_2=+iw_1$ となります。そこで規格化して

$$\boldsymbol{w}_{\pm}=\frac{1}{\sqrt{2}}\begin{pmatrix} 1 \\ \mp i \end{pmatrix}$$

とおきます。この $\boldsymbol{w}_\pm$ が固有ベクトルであることは、$\sin\theta=0$ の場合も含めて直接確認できます。

$$B\begin{pmatrix} 1 \\ -i \end{pmatrix} =\begin{pmatrix} \cos\theta+i\sin\theta \\ \sin\theta-i\cos\theta \end{pmatrix} =e^{i\theta}\begin{pmatrix} 1 \\ -i \end{pmatrix}, \qquad B\begin{pmatrix} 1 \\ i \end{pmatrix} =\begin{pmatrix} \cos\theta-i\sin\theta \\ \sin\theta+i\cos\theta \end{pmatrix} =e^{-i\theta}\begin{pmatrix} 1 \\ i \end{pmatrix}$$

（第2成分は $-i(\cos\theta+i\sin\theta)=\sin\theta-i\cos\theta$、$i(\cos\theta-i\sin\theta)=\sin\theta+i\cos\theta$ から従います。）これらを列に並べて

$$U=\frac{1}{\sqrt{2}}\begin{pmatrix} 1 & 1 \\ -i & i \end{pmatrix}$$

とおくと、$U$ はユニタリです。

$$U^{*}U=\frac{1}{2}\begin{pmatrix} 1 & i \\ 1 & -i \end{pmatrix}\begin{pmatrix} 1 & 1 \\ -i & i \end{pmatrix} =\frac{1}{2}\begin{pmatrix} 2 & 0 \\ 0 & 2 \end{pmatrix}=E$$

したがって $U^{-1}=U^{*}$ であり、固有値関係 $BU=U\,\mathrm{diag}(e^{i\theta},e^{-i\theta})$ から

$$U^{-1}BU=U^{*}BU=\begin{pmatrix} e^{i\theta} & 0 \\ 0 & e^{-i\theta} \end{pmatrix}$$

が答えです。$\theta$ が $\pi$ の整数倍のときは $B=\pm E$ ですでに対角ですが、そのときも $e^{i\theta}=e^{-i\theta}=\pm1$ なので上の式はそのまま成り立ちます。つまり $\theta$ の値による場合分けは不要です。

固有値の絶対値が $\lvert e^{\pm i\theta}\rvert=1$ であることは $B$ がユニタリ行列であることと整合し、積 $e^{i\theta}e^{-i\theta}=1$ は $\lvert B \rvert=1$ と、和 $e^{i\theta}+e^{-i\theta}=2\cos\theta$ はトレースと一致します。

同じ $U$ は設問2の $A=\theta J$ も対角化します。$J\boldsymbol{w}_\pm=\mp i\boldsymbol{w}_\pm$ なので $U^{-1}AU=\mathrm{diag}(i\theta,-i\theta)$ であり、対角行列に対しては指数関数が成分ごとの指数関数になるので

$$U^{-1}(\exp A)U=\exp\!\left(U^{-1}AU\right)=\begin{pmatrix} e^{i\theta} & 0 \\ 0 & e^{-i\theta} \end{pmatrix}=U^{-1}BU$$

となり、$\exp A=B$ が再現されます。設問2の別解であり、検算にもなっています（$U^{-1}(\exp A)U=\exp(U^{-1}AU)$ は級数の各項に $U^{-1}\cdot U$ を挟めば従います）。

設問4

$D_n(x_1,\dots,x_n)\equiv\lvert C \rvert$ とおき、$n$ についての帰納法で

$$D_n=\prod_{1\le i<j\le n}(x_j-x_i)$$

を示します。比例係数がちょうど 1 になることまで出ます。

$n=1$ のとき $C=(1)$ で $D_1=1$、右辺は空積なので 1 で一致します。

$n\ge 2$ とし、$D_{n-1}$ について主張が成り立つと仮定します。行番号 $k$ を $k=n,n-1,\dots,2$ の順に（大きい方から）とり、第 $k$ 行から第 $k-1$ 行の $x_1$ 倍を引きます。降順に処理するので、引く側の第 $k-1$ 行はまだ変更されていない元の行です。行に他の行の定数倍を加える操作は行列式を変えないので $D_n$ は不変です。第 $k$ 行 $(k\ge2)$ の第 $j$ 成分は

$$x_j^{\,k-1}-x_1 x_j^{\,k-2}=x_j^{\,k-2}\,(x_j-x_1)$$

になります。とくに $j=1$ では $x_1^{\,k-2}(x_1-x_1)=0$ です。第1行は $(1,1,\dots,1)$ のまま変わりません。したがって第1列は $(1,0,\dots,0)^{\mathsf{T}}$ となり、第1列に沿う余因子展開で第1項だけが残ります（符号は $(-1)^{1+1}=+1$）。

$$D_n=\det M,\qquad M_{k-1,\,j-1}=x_j^{\,k-2}(x_j-x_1)\quad (k,j=2,\dots,n)$$

$M$ は $(n-1)\times(n-1)$ 行列で、その第 $j-1$ 列は共通因子 $(x_j-x_1)$ を持ちます。行列式は各列について1次なので、列ごとに因子をくくり出せて

$$D_n=\left\{\prod_{j=2}^{n}(x_j-x_1)\right\}\det\left[x_j^{\,k-2}\right]_{k,j=2,\dots,n} =\left\{\prod_{j=2}^{n}(x_j-x_1)\right\}D_{n-1}(x_2,\dots,x_n)$$

となります。最後の行列式は、行の指数が $0,1,\dots,n-2$、列の変数が $x_2,\dots,x_n$ という、まさに $n-1$ 変数の同じ形の行列式です。帰納法の仮定を使えば

$$D_n=\prod_{j=2}^{n}(x_j-x_1)\cdot\prod_{2\le i<j\le n}(x_j-x_i)=\prod_{1\le i<j\le n}(x_j-x_i)$$

が得られ、主張が示されました。$i=1$ の因子と $i\ge2$ の因子を合わせるとちょうどすべての組 $(i,j)$、$i<j$ を一度ずつ拾うことに注意します。

別の見方でも同じ結論が出ます。$D_n$ は $x_1,\dots,x_n$ の多項式で、$x_i=x_j$（$i\neq j$）とおくと第 $i$ 列と第 $j$ 列が一致して $D_n=0$ になるので、因数定理から各 $(x_j-x_i)$ で割り切れます。これらは互いに素な既約多項式なので積 $\prod_{i<j}(x_j-x_i)$ で割り切れます。一方 $D_n$ の全次数は各項が行ごとに $0,1,\dots,n-1$ 次の因子を1つずつ拾うので $\sum_{k=0}^{n-1}k=n(n-1)/2$ で、積の次数も組の個数 $\binom{n}{2}=n(n-1)/2$ に等しい。よって商は定数です。定数は対角成分の積に対応する単項式 $x_2x_3^2\cdots x_n^{\,n-1}$ の係数を両辺で比べれば 1 と決まります。

小さい $n$ で確かめると、$n=2$ では $\lvert C \rvert=x_2-x_1$、$n=3$ では $\lvert C \rvert=(x_2-x_1)(x_3-x_1)(x_3-x_2)$ で、いずれも公式と一致します。乱数を入れた数値計算でも $n=2$ から $6$ まで両辺が一致します。

設問5

設問4の結果 $\lvert C \rvert=\prod_{1\le i<j\le n}(x_j-x_i)$ から、$\lvert C \rvert\neq 0$ であることは、積のどの因子も 0 でないこと、すなわちすべての組 $i<j$ について $x_i\neq x_j$ であることと同値です。したがって $C$ が正則である（逆行列をもつ）のは、$x_1,\dots,x_n$ が互いにすべて異なる場合であり、そのときに限ります。

逆向きも直接見えます。もし $x_i=x_j$ となる $i\neq j$ があれば $C$ の第 $i$ 列と第 $j$ 列が完全に一致するので、列が線形従属になり $\lvert C \rvert=0$ で $C$ は正則ではありません。なお $n=1$ のときは条件が空で、$\lvert C \rvert=1$ なので常に正則です。

この条件の意味づけを付けておきます。$C$ の転置 $C^{\mathsf{T}}$ は、係数 $(a_0,a_1,\dots,a_{n-1})$ の多項式 $P(x)=\sum_{k=0}^{n-1}a_k x^k$ を値の組 $(P(x_1),\dots,P(x_n))$ に写す行列です。$C$ が正則であることは、与えられた $n$ 個の値を $n$ 点 $x_1,\dots,x_n$ で取る $n-1$ 次以下の多項式がただ一つ存在すること（ラグランジュ補間の一意可解性）と同じ内容です。点が重複していれば、その点で異なる値を要求すると解がなく、同じ値を要求すると解が一意でなくなります。互いに異なる、という条件はこの事情をそのまま表しています。

第2問 コーシー–リーマン関係式、コーシーの積分定理、ガウス型積分

設定

2つの実変数 $x,y$ の実関数 $u(x,y)$、$v(x,y)$ は、平面上のいたるところで有限な値をとり、何回でも微分可能で、さらに

$$\frac{\partial u}{\partial x}=\frac{\partial v}{\partial y}, \qquad \frac{\partial u}{\partial y}=-\frac{\partial v}{\partial x}$$

を満たすとします。これがコーシー–リーマン関係式で、以下ではこの2式を式(1)と呼びます。無限回微分可能という仮定は平面全体で成り立つとします。

この仮定だけを出発点に、$u,v$ が調和関数であること、$xy$ 平面上の任意の閉経路まわりの2つの線積分が消えること（式(4)）、$f(z)=u+iv$ に対する $\oint_{\mathrm{C}}f(z)dz=0$（式(5)）を順に導き、最後に $f(z)=e^{-z^2}$ に適用して定積分を求めます。

設問1

式(1)の第1式を $x$ で偏微分し、第2式を $y$ で偏微分します。

$$\frac{\partial^2 u}{\partial x^2}=\frac{\partial}{\partial x}\frac{\partial v}{\partial y}=\frac{\partial^2 v}{\partial x\,\partial y}, \qquad \frac{\partial^2 u}{\partial y^2}=\frac{\partial}{\partial y}\left(-\frac{\partial v}{\partial x}\right)=-\frac{\partial^2 v}{\partial y\,\partial x}$$

$v$ は任意回微分可能で、とくに2階偏導関数が連続なので、シュワルツ（クレロー）の定理により偏微分の順序を交換できます。

$$\frac{\partial^2 v}{\partial x\,\partial y}=\frac{\partial^2 v}{\partial y\,\partial x}$$

したがって

$$\frac{\partial^2 u}{\partial x^2}+\frac{\partial^2 u}{\partial y^2} =\frac{\partial^2 v}{\partial x\,\partial y}-\frac{\partial^2 v}{\partial y\,\partial x}=0$$

です。$v$ についても同様で、式(1)の第2式を $x$ で、第1式を $y$ で偏微分すると

$$\frac{\partial^2 v}{\partial x^2}=\frac{\partial}{\partial x}\left(-\frac{\partial u}{\partial y}\right)=-\frac{\partial^2 u}{\partial x\,\partial y}, \qquad \frac{\partial^2 v}{\partial y^2}=\frac{\partial}{\partial y}\frac{\partial u}{\partial x}=\frac{\partial^2 u}{\partial y\,\partial x}$$

となり、今度は $u$ の2階偏導関数の連続性から交換できて

$$\frac{\partial^2 v}{\partial x^2}+\frac{\partial^2 v}{\partial y^2} =-\frac{\partial^2 u}{\partial x\,\partial y}+\frac{\partial^2 u}{\partial y\,\partial x}=0$$

を得ます。以上で式(2)が示されました。使った仮定は、式(1)と、$u,v$ が $C^2$ 級であること（無限回微分可能性から従う）の2つだけです。$u,v$ はともにラプラス方程式を満たす調和関数です。

設問2

$xy$ 平面内の経路なので、面 $\mathrm{S}$ も $xy$ 平面内にとれます。このとき法線ベクトルは $\boldsymbol{n}=\boldsymbol{e}_z$（経路 $\mathrm{C}$ を反時計回りに一周する向きに合わせた右ねじの向き）です。$z$ 成分をもたず $x,y$ にしか依存しないベクトル場 $\boldsymbol{A}=(A_x(x,y),A_y(x,y),0)$ については

$$\boldsymbol{A}\cdot d\boldsymbol{x}=A_x\,dx+A_y\,dy, \qquad \left[\nabla\times\boldsymbol{A}\right]\cdot\boldsymbol{n}=\frac{\partial A_y}{\partial x}-\frac{\partial A_x}{\partial y}$$

となります。あとは求める線積分の被積分形になるように $\boldsymbol{A}$ を選ぶだけです。$u,v$ は平面全体で無限回微分可能なので、以下のどちらの $\boldsymbol{A}$ も微分可能で、除外すべき特異点はなく、$\mathrm{C}$ が囲む領域 $\mathrm{S}$ の全体でストークスの定理を使えます。

(ア) 第1式。$\boldsymbol{A}=(u,-v,0)$ とおくと $\boldsymbol{A}\cdot d\boldsymbol{x}=u\,dx-v\,dy$ です。回転の $z$ 成分は

$$\frac{\partial(-v)}{\partial x}-\frac{\partial u}{\partial y} =-\left(\frac{\partial v}{\partial x}+\frac{\partial u}{\partial y}\right)=0$$

で、最後の等号は式(1)の第2式 $\partial u/\partial y=-\partial v/\partial x$ そのものです。よってストークスの定理から

$$\oint_{\mathrm{C}}\left[u(x,y)dx-v(x,y)dy\right]=\iint_{\mathrm{S}}0\,dS=0$$

となります。

(イ) 第2式。$\boldsymbol{A}=(v,u,0)$ とおくと $\boldsymbol{A}\cdot d\boldsymbol{x}=v\,dx+u\,dy$ です。回転の $z$ 成分は

$$\frac{\partial u}{\partial x}-\frac{\partial v}{\partial y}=0$$

で、これは式(1)の第1式です。よって

$$\oint_{\mathrm{C}}\left[u(x,y)dy+v(x,y)dx\right]=\iint_{\mathrm{S}}0\,dS=0$$

となり、式(4)の2式がともに示されました。

面積分の被積分関数が恒等的に 0 なので、$\mathrm{S}$ の取り方や $\boldsymbol{n}$ の向き（$\pm\boldsymbol{e}_z$ のどちらにとるか、それに応じて $\mathrm{C}$ の向きを合わせる）には結果が依存しません。また、経路が自己交差する一般の閉じた経路の場合は、有限個の単純閉曲線に分解して各々に上の議論を適用すればよく、和も 0 です。別の言い方をすると、$\nabla\times\boldsymbol{A}$ が単連結な平面全体で恒等的に 0 なので $\boldsymbol{A}$ は保存場（あるスカラー関数の勾配）であり、その周回積分はどんな閉経路でも 0 になります。式(4)が「任意の」閉経路について成り立つ、という主張はこの意味で正しいことになります。

設問3

$z=x+iy$ なので $dz=dx+i\,dy$、また $f(z)=u(x,y)+iv(x,y)$ です。被積分形を実部と虚部に分けます。

$$f(z)\,dz=(u+iv)(dx+i\,dy) =\left(u\,dx-v\,dy\right)+i\left(v\,dx+u\,dy\right)$$

実部では $u\,dx$ と $i^2 v\,dy=-v\,dy$ が、虚部では $v\,dx$ と $u\,dy$ が集まります。経路 $\mathrm{C}$ 上で積分すると

$$\oint_{\mathrm{C}}f(z)\,dz =\oint_{\mathrm{C}}\left[u\,dx-v\,dy\right]+i\oint_{\mathrm{C}}\left[v\,dx+u\,dy\right]$$

となり、右辺の2つの周回積分は設問2で示した式(4)によりどちらも 0 です。よって

$$\oint_{\mathrm{C}}f(z)\,dz=0+i\cdot 0=0$$

が任意の閉経路 $\mathrm{C}$ について成り立ちます。これがコーシーの積分定理で、成立の根拠は式(1)（コーシー–リーマン関係式）と、$u,v$ が平面全体で滑らかで特異点をもたないことです。$f$ が正則でない領域を経路が囲む場合には、この論法は使えません。

設問4

$z=x+iy$ に対し $z^2=x^2-y^2+2ixy$ なので

$$e^{-z^2}=e^{-(x^2-y^2)}e^{-2ixy}=e^{\,y^2-x^2}\left[\cos(2xy)-i\sin(2xy)\right]$$

です。実部と虚部を読み取ると

$$u(x,y)=e^{\,y^2-x^2}\cos(2xy), \qquad v(x,y)=-e^{\,y^2-x^2}\sin(2xy)$$

となります。どちらも各点で有限な値をとり、初等関数の合成なので無限回微分可能です。以下 $E(x,y)\equiv e^{\,y^2-x^2}$ と略記します。$\partial E/\partial x=-2xE$、$\partial E/\partial y=2yE$ です。

$$\begin{aligned} \frac{\partial u}{\partial x}&=-2xE\cos(2xy)+E\cdot\left(-2y\sin(2xy)\right) =-2E\left[x\cos(2xy)+y\sin(2xy)\right]\\[2pt] \frac{\partial v}{\partial y}&=-2yE\sin(2xy)-E\cdot 2x\cos(2xy) =-2E\left[y\sin(2xy)+x\cos(2xy)\right] \end{aligned}$$

2つは一致するので $\partial u/\partial x=\partial v/\partial y$ です。次に

$$\begin{aligned} \frac{\partial u}{\partial y}&=2yE\cos(2xy)+E\cdot\left(-2x\sin(2xy)\right) =2E\left[y\cos(2xy)-x\sin(2xy)\right]\\[2pt] \frac{\partial v}{\partial x}&=2xE\sin(2xy)-E\cdot 2y\cos(2xy) =-2E\left[y\cos(2xy)-x\sin(2xy)\right] \end{aligned}$$

より $\partial u/\partial y=-\partial v/\partial x$ です。以上で $f(z)=e^{-z^2}$ の実部と虚部が式(1)を満たすことが示されました。数値微分でも $(x,y)=(0.3,0.7)$、$(-1.2,0.4)$、$(2.0,-1.1)$ の3点で両式が成り立つことを確認しています。

したがって設問1から設問3の結論がこの $f$ にすべて適用できます。とくに $u,v$ は調和で、任意の閉経路について $\oint_{\mathrm{C}}e^{-z^2}dz=0$ が成り立ちます。

設問5

まず $p>0$ とします。$f(z)=e^{-z^2}$ に対し、複素平面上の長方形の閉経路

$$\mathrm{C}_R:\quad -R\to R\to R+ip\to -R+ip\to -R$$

（$R>0$、反時計回り）をとります。設問4より $f$ の実部・虚部は式(1)を満たすので、設問3から

$$\oint_{\mathrm{C}_R}e^{-z^2}dz=0$$

です。この周回積分を4辺に分けます。

下辺（$z=x$、$x$ は $-R$ から $R$）は

$$I_1=\int_{-R}^{R}e^{-x^2}dx$$

上辺（$z=x+ip$、$x$ は $R$ から $-R$）は、$-(x+ip)^2=-x^2-2ipx+p^2$ より

$$I_3=\int_{R}^{-R}e^{p^2}e^{-x^2}e^{-2ipx}dx=-e^{p^2}\int_{-R}^{R}e^{-x^2}e^{-2ipx}dx$$

右辺（$z=R+iy$、$y$ は $0$ から $p$、$dz=i\,dy$）と左辺（$z=-R+iy$、$y$ は $p$ から $0$）については、$-(\pm R+iy)^2=-(R^2-y^2)\mp 2iRy$ なので被積分関数の絶対値は

$$\left\lvert e^{-(\pm R+iy)^2}\right\rvert=e^{\,y^2-R^2}\le e^{\,p^2-R^2} \qquad (0\le y\le p)$$

と評価できます。積分区間の長さは $p$ なので

$$\lvert I_2 \rvert\le p\,e^{p^2}e^{-R^2}, \qquad \lvert I_4 \rvert\le p\,e^{p^2}e^{-R^2}$$

であり、$p$ を固定して $R\to\infty$ とすると $e^{-R^2}\to0$ なのでどちらも 0 に収束します。ここで $p$ は $R$ に依らない定数であることが効いています。

$R\to\infty$ の極限をとります。$I_1$ は式(7)により $\sqrt{\pi}$ に収束し、$I_3$ の積分も $\lvert e^{-x^2}e^{-2ipx}\rvert=e^{-x^2}$ が可積分なので絶対収束します。$I_1+I_2+I_3+I_4=0$ で $I_2,I_4\to0$ だから

$$\sqrt{\pi}-e^{p^2}\int_{-\infty}^{\infty}e^{-x^2}e^{-2ipx}dx=0, \qquad\text{すなわち}\qquad \int_{-\infty}^{\infty}e^{-x^2}e^{-2ipx}dx=\sqrt{\pi}\,e^{-p^2}$$

を得ます。$e^{-2ipx}=\cos(2px)-i\sin(2px)$ と分けると、虚部の積分は $e^{-x^2}\sin(2px)$ が $x$ の奇関数で（絶対収束もしているので）0 です。実部を比べて

$$\int_{-\infty}^{\infty}e^{-x^2}\cos(2px)\,dx=\sqrt{\pi}\,e^{-p^2}$$

となります。

残りの場合を埋めます。$p=0$ のときは式(7)そのもので、上式も $\sqrt{\pi}$ を与えるので一致します。$p<0$ のときは $\cos(2px)$ が $p$ について偶関数なので左辺は $\lvert p \rvert$ での値に等しく、右辺 $\sqrt{\pi}e^{-p^2}$ も $p$ の偶関数です。したがって

$$\int_{-\infty}^{\infty}e^{-x^2}\cos(2px)\,dx=\sqrt{\pi}\,e^{-p^2}$$

がすべての実数 $p$ について成り立ちます。これが答えです。$p<0$ を長方形経路で直接扱うなら、$-\lvert p \rvert$ 側、つまり下半平面に長方形をとれば同じ計算になります。

検算を2つ付けます。ひとつは微分方程式による別解です。$I(p)\equiv\int_{-\infty}^{\infty}e^{-x^2}\cos(2px)dx$ とおくと、$p$ による偏導関数の絶対値が $\lvert 2xe^{-x^2}\sin(2px)\rvert\le 2\lvert x \rvert e^{-x^2}$ で $p$ に依らない可積分関数に押さえられるので、積分と微分を交換できます。$2xe^{-x^2}=-\frac{d}{dx}e^{-x^2}$ を使って部分積分すると

$$I'(p)=-\int_{-\infty}^{\infty}2x e^{-x^2}\sin(2px)dx =\left[e^{-x^2}\sin(2px)\right]_{-\infty}^{\infty}-2p\int_{-\infty}^{\infty}e^{-x^2}\cos(2px)dx=-2p\,I(p)$$

です（境界項は $e^{-x^2}\to0$ で消えます）。$I(0)=\sqrt{\pi}$ の初期条件で解くと $I(p)=\sqrt{\pi}e^{-p^2}$ で、上の結果と一致します。もうひとつは数値積分です。$p=0.5,\,1.0,\,2.0,\,-1.5$ で数値積分すると $1.38039$、$0.65205$、$0.032464$、$0.186815$ となり、$\sqrt{\pi}e^{-p^2}$ の値と有効数字5桁以上で一致します。

得られた式は、ガウス関数のフーリエ変換がまたガウス関数になるという事実にほかなりません。$p$ が大きくなるほど振動が細かくなって打ち消し合いが効き、積分値は $e^{-p^2}$ で急速に減衰します。

出典: 東京大学大学院理学系研究科 物理学専攻 平成19年度 修士課程 入学試験問題 数学。問題文は要約して引用しています。