コンテンツにスキップ

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

対象読者: 学部の物理数学(線形代数、微分積分、常微分方程式、複素関数論、フーリエ解析)を一通り学び、東大物理学専攻の院試問題を自分で解いた、または解こうとしている人。問題文は載せていないので、出典の問題を手元に置いて読みます。

生 Markdown

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

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

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

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

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

A∗A=AA∗A^{*}A=AA^{*}

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

前半では 2 次の行列

B=(cos⁡θ−sin⁡θsin⁡θcos⁡θ)(θ∈R)B=\begin{pmatrix} \cos\theta & -\sin\theta \\ \sin\theta & \cos\theta \end{pmatrix} \qquad (\theta \in \mathbb{R})

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

C=(11⋯1x1x2⋯xn⋮⋮⋮x1 n−1x2 n−1⋯xn n−1),Cij=xj i−1C=\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}

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

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

BT=(cos⁡θsin⁡θ−sin⁡θcos⁡θ)B^{\mathsf{T}}=\begin{pmatrix} \cos\theta & \sin\theta \\ -\sin\theta & \cos\theta \end{pmatrix}

です。積を計算すると

BTB=(cos⁡θsin⁡θ−sin⁡θcos⁡θ)(cos⁡θ−sin⁡θsin⁡θcos⁡θ)=(cos⁡2θ+sin⁡2θ−cos⁡θsin⁡θ+sin⁡θcos⁡θ−sin⁡θcos⁡θ+cos⁡θsin⁡θsin⁡2θ+cos⁡2θ)=(1001)=E\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}

となり、同じ計算で

BBT=(cos⁡2θ+sin⁡2θcos⁡θsin⁡θ−sin⁡θcos⁡θsin⁡θcos⁡θ−cos⁡θsin⁡θsin⁡2θ+cos⁡2θ)=EBB^{\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θ+sin⁡2θ=1\cos^2\theta+\sin^2\theta=1 だけです。

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

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

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

A=(0−θθ0)=θJ,J≡(0−110)A=\begin{pmatrix} 0 & -\theta \\ \theta & 0 \end{pmatrix}=\theta J, \qquad J\equiv\begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix}

とおくと

J2=(0−110)(0−110)=(−100−1)=−EJ^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,…k=0,1,2,\dots に対して

J2k=(−1)kE,J2k+1=(−1)kJJ^{2k}=(-1)^k E,\qquad J^{2k+1}=(-1)^k J

が成り立ちます(kk についての帰納法。J2(k+1)=J2kJ2=(−1)kE⋅(−E)=(−1)k+1EJ^{2(k+1)}=J^{2k}J^2=(-1)^k E\cdot(-E)=(-1)^{k+1}E)。

Am=θmJmA^m=\theta^m J^m なので、定義式は

exp⁡A=∑m=0∞θmJmm!\exp A=\sum_{m=0}^{\infty}\frac{\theta^m J^m}{m!}

となります。この級数は絶対収束します。実際、行列ノルムを劣乗法的なもの(たとえば作用素ノルム。劣乗法的とは ∥XY∥≤∥X∥∥Y∥\lVert XY\rVert\le\lVert X\rVert\lVert Y\rVert が成り立つことで、作用素ノルムは行列が長さ 1 のベクトルを最大で何倍に伸ばすかを表す量のことです)にとれば ∥Am/m!∥≤∥A∥m/m!\lVert A^m/m!\rVert \le \lVert A\rVert^m/m! であり、∑m∥A∥m/m!=e∥A∥<∞\sum_m \lVert A\rVert^m/m!=e^{\lVert A\rVert}<\infty です。2×22\times 2 行列の空間は有限次元なので完備であり、絶対収束する級数は収束します。絶対収束していれば項の並べ替えが自由にできるので、mm を偶数 m=2km=2k と奇数 m=2k+1m=2k+1 に分けて

exp⁡A=∑k=0∞θ2k(2k)!J2k+∑k=0∞θ2k+1(2k+1)!J2k+1=(∑k=0∞(−1)kθ2k(2k)!)E+(∑k=0∞(−1)kθ2k+1(2k+1)!)J=(cos⁡θ) E+(sin⁡θ) J\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⁡\cos と sin⁡\sin のテイラー級数(実数 θ\theta 全体で収束)です。成分に戻すと

exp⁡A=cos⁡θ(1001)+sin⁡θ(0−110)=(cos⁡θ−sin⁡θsin⁡θcos⁡θ)=B\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

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

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

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

∣B−λE∣=(cos⁡θ−λ)2+sin⁡2θ=λ2−2λcos⁡θ+1=0\lvert B-\lambda E \rvert=(\cos\theta-\lambda)^2+\sin^2\theta=\lambda^2-2\lambda\cos\theta+1=0

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

λ±=cos⁡θ±isin⁡θ=e±iθ\lambda_{\pm}=\cos\theta\pm i\sin\theta=e^{\pm i\theta}

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

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

w±=12(1∓i)\boldsymbol{w}_{\pm}=\frac{1}{\sqrt{2}}\begin{pmatrix} 1 \\ \mp i \end{pmatrix}

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

B(1−i)=(cos⁡θ+isin⁡θsin⁡θ−icos⁡θ)=eiθ(1−i),B(1i)=(cos⁡θ−isin⁡θsin⁡θ+icos⁡θ)=e−iθ(1i)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⁡θ+isin⁡θ)=sin⁡θ−icos⁡θ-i(\cos\theta+i\sin\theta)=\sin\theta-i\cos\theta、i(cos⁡θ−isin⁡θ)=sin⁡θ+icos⁡θi(\cos\theta-i\sin\theta)=\sin\theta+i\cos\theta から従います。)これらを列に並べて

U=12(11−ii)U=\frac{1}{\sqrt{2}}\begin{pmatrix} 1 & 1 \\ -i & i \end{pmatrix}

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

U∗U=12(1i1−i)(11−ii)=12(2002)=EU^{*}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∗U^{-1}=U^{*} であり、固有値関係 BU=U diag(eiθ,e−iθ)BU=U\,\mathrm{diag}(e^{i\theta},e^{-i\theta}) から

U−1BU=U∗BU=(eiθ00e−iθ)U^{-1}BU=U^{*}BU=\begin{pmatrix} e^{i\theta} & 0 \\ 0 & e^{-i\theta} \end{pmatrix}

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

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

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

U−1(exp⁡A)U=exp⁡ ⁣(U−1AU)=(eiθ00e−iθ)=U−1BUU^{-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\exp A=B が再現されます。設問2の別解であり、検算にもなっています(U−1(exp⁡A)U=exp⁡(U−1AU)U^{-1}(\exp A)U=\exp(U^{-1}AU) は級数の各項に U−1⋅UU^{-1}\cdot U を挟めば従います)。

Dn(x1,…,xn)≡∣C∣D_n(x_1,\dots,x_n)\equiv\lvert C \rvert とおき、nn についての帰納法で

Dn=∏1≤i<j≤n(xj−xi)D_n=\prod_{1\le i<j\le n}(x_j-x_i)

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

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

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

xj k−1−x1xj k−2=xj k−2 (xj−x1)x_j^{\,k-1}-x_1 x_j^{\,k-2}=x_j^{\,k-2}\,(x_j-x_1)

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

Dn=det⁡M,Mk−1, j−1=xj k−2(xj−x1)(k,j=2,…,n)D_n=\det M,\qquad M_{k-1,\,j-1}=x_j^{\,k-2}(x_j-x_1)\quad (k,j=2,\dots,n)

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

Dn={∏j=2n(xj−x1)}det⁡[xj k−2]k,j=2,…,n={∏j=2n(xj−x1)}Dn−1(x2,…,xn)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,…,n−20,1,\dots,n-2、列の変数が x2,…,xnx_2,\dots,x_n という、まさに n−1n-1 変数の同じ形の行列式です。帰納法の仮定を使えば

Dn=∏j=2n(xj−x1)⋅∏2≤i<j≤n(xj−xi)=∏1≤i<j≤n(xj−xi)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=1i=1 の因子と i≥2i\ge2 の因子を合わせるとちょうどすべての組 (i,j)(i,j)、i<ji<j を一度ずつ拾うことに注意します。

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

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

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

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

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

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

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

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

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

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

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

式(1)の第1式を xx で偏微分し、第2式を yy で偏微分します。

∂2u∂x2=∂∂x∂v∂y=∂2v∂x ∂y,∂2u∂y2=∂∂y(−∂v∂x)=−∂2v∂y ∂x\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}

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

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

したがって

∂2u∂x2+∂2u∂y2=∂2v∂x ∂y−∂2v∂y ∂x=0\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

です。vv についても同様で、式(1)の第2式を xx で、第1式を yy で偏微分すると

∂2v∂x2=∂∂x(−∂u∂y)=−∂2u∂x ∂y,∂2v∂y2=∂∂y∂u∂x=∂2u∂y ∂x\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}

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

∂2v∂x2+∂2v∂y2=−∂2u∂x ∂y+∂2u∂y ∂x=0\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,vu,v が C2C^2 級であること(無限回微分可能性から従う)の2つだけです。u,vu,v はともにラプラス方程式を満たす調和関数です。

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

A⋅dx=Ax dx+Ay dy,[∇×A]⋅n=∂Ay∂x−∂Ax∂y\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}

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

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

∂(−v)∂x−∂u∂y=−(∂v∂x+∂u∂y)=0\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式 ∂u/∂y=−∂v/∂x\partial u/\partial y=-\partial v/\partial x そのものです。よってストークスの定理から

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

となります。

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

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

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

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

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

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

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

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

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

∮Cf(z) dz=∮C[u dx−v dy]+i∮C[v dx+u dy]\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 です。よって

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

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

z=x+iyz=x+iy に対し z2=x2−y2+2ixyz^2=x^2-y^2+2ixy なので

e−z2=e−(x2−y2)e−2ixy=e y2−x2[cos⁡(2xy)−isin⁡(2xy)]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 y2−x2cos⁡(2xy),v(x,y)=−e y2−x2sin⁡(2xy)u(x,y)=e^{\,y^2-x^2}\cos(2xy), \qquad v(x,y)=-e^{\,y^2-x^2}\sin(2xy)

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

∂u∂x=−2xEcos⁡(2xy)+E⋅(−2ysin⁡(2xy))=−2E[xcos⁡(2xy)+ysin⁡(2xy)]∂v∂y=−2yEsin⁡(2xy)−E⋅2xcos⁡(2xy)=−2E[ysin⁡(2xy)+xcos⁡(2xy)]\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つは一致するので ∂u/∂x=∂v/∂y\partial u/\partial x=\partial v/\partial y です。次に

∂u∂y=2yEcos⁡(2xy)+E⋅(−2xsin⁡(2xy))=2E[ycos⁡(2xy)−xsin⁡(2xy)]∂v∂x=2xEsin⁡(2xy)−E⋅2ycos⁡(2xy)=−2E[ycos⁡(2xy)−xsin⁡(2xy)]\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}

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

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

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

CR:−R→R→R+ip→−R+ip→−R\mathrm{C}_R:\quad -R\to R\to R+ip\to -R+ip\to -R

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

∮CRe−z2dz=0\oint_{\mathrm{C}_R}e^{-z^2}dz=0

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

下辺(z=xz=x、xx は −R-R から RR)は

I1=∫−RRe−x2dxI_1=\int_{-R}^{R}e^{-x^2}dx

上辺(z=x+ipz=x+ip、xx は RR から −R-R)は、−(x+ip)2=−x2−2ipx+p2-(x+ip)^2=-x^2-2ipx+p^2 より

I3=∫R−Rep2e−x2e−2ipxdx=−ep2∫−RRe−x2e−2ipxdxI_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+iyz=R+iy、yy は 00 から pp、dz=i dydz=i\,dy)と左辺(z=−R+iyz=-R+iy、yy は pp から 00)については、−(±R+iy)2=−(R2−y2)∓2iRy-(\pm R+iy)^2=-(R^2-y^2)\mp 2iRy なので被積分関数の絶対値は

∣e−(±R+iy)2∣=e y2−R2≤e p2−R2(0≤y≤p)\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)

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

∣I2∣≤p ep2e−R2,∣I4∣≤p ep2e−R2\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}

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

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

π−ep2∫−∞∞e−x2e−2ipxdx=0,すなわち∫−∞∞e−x2e−2ipxdx=π e−p2\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)−isin⁡(2px)e^{-2ipx}=\cos(2px)-i\sin(2px) と分けると、虚部の積分は e−x2sin⁡(2px)e^{-x^2}\sin(2px) が xx の奇関数で(絶対収束もしているので)0 です。実部を比べて

∫−∞∞e−x2cos⁡(2px) dx=π e−p2\int_{-\infty}^{\infty}e^{-x^2}\cos(2px)\,dx=\sqrt{\pi}\,e^{-p^2}

となります。

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

∫−∞∞e−x2cos⁡(2px) dx=π e−p2\int_{-\infty}^{\infty}e^{-x^2}\cos(2px)\,dx=\sqrt{\pi}\,e^{-p^2}

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

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

I′(p)=−∫−∞∞2xe−x2sin⁡(2px)dx=[e−x2sin⁡(2px)]−∞∞−2p∫−∞∞e−x2cos⁡(2px)dx=−2p I(p)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−x2→0e^{-x^2}\to0 で消えます)。I(0)=πI(0)=\sqrt{\pi} の初期条件で解くと I(p)=πe−p2I(p)=\sqrt{\pi}e^{-p^2} で、上の結果と一致します。もうひとつは数値積分です。p=0.5, 1.0, 2.0, −1.5p=0.5,\,1.0,\,2.0,\,-1.5 で数値積分すると 1.380391.38039、0.652050.65205、0.0324640.032464、0.1868150.186815 となり、πe−p2\sqrt{\pi}e^{-p^2} の値と有効数字5桁以上で一致します。

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

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

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

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

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