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

> 平成8年度 修士課程 数学（全3問）の解答。非線形振動子の摂動展開と永年項、3行2列行列の特異値分解、ビュフォンの針（長針の場合）を扱います。
> https://rikai.mugen-giken.com/exams/utokyo-physics/master-h8-math

:::caution[出典と、この解答の位置づけ]
本記事が扱う問題は、**東京大学大学院理学系研究科物理学専攻の大学院入試問題**からの引用です（引用元: [https://www.phys.s.u-tokyo.ac.jp/](https://www.phys.s.u-tokyo.ac.jp/)）。

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

平成8年度（平成7年8月実施）の数学は全3問で、3問すべてに解答する形式です。第1問は非線形振動子の摂動展開で、永年項の出現と展開の破綻という摂動論の定番テーマを問います。第2問は3行2列行列の特異値分解を具体的な行列で実行させる問題、第3問はビュフォンの針で、棒の長さが線の間隔以上という「長針」の場合を端点の位置から積み上げて計算させます。いずれも方針は素直で、計算の正確さと議論の詰めが得点を分けます。

なお、この問題冊子には物理学の問題も同綴されていますが、本記事では数学の3問のみを扱います。

| 問題 | 分野 | 主題 |
|---|---|---|
| 第1問 | 微分方程式 | 非線形振動子の摂動展開と永年項 |
| 第2問 | 線形代数 | 3行2列行列の特異値分解 |
| 第3問 | 確率統計 | ビュフォンの針（長針の場合） |

## 第1問 非線形振動子の摂動展開と永年項

### 設定

時間 $t$ の関数 $y(t)$ に対する非線形常微分方程式

$$
\frac{d^2y}{dt^2}+y=\varepsilon y^3
$$

を、初期条件 $y(0)=1$, $\dfrac{dy}{dt}(0)=0$ のもとで考えます。$\varepsilon$ は微小な正の定数です。必要なら恒等式 $(\cos\theta+i\sin\theta)^3=\cos 3\theta+i\sin 3\theta$ を用いてよい、とあります。

### 設問(a)

解を $\varepsilon$ のべき級数

$$
y(t;\varepsilon)=\sum_{n=0}^{\infty}\varepsilon^n y_n(t)
$$

と展開して方程式に代入します。右辺は

$$
\varepsilon y^3=\varepsilon\,(y_0+\varepsilon y_1+\cdots)^3=\varepsilon y_0^3+O(\varepsilon^2)
$$

なので、$\varepsilon$ の各次数を比較すると

$$
\begin{aligned}
O(\varepsilon^0):&\quad \frac{d^2y_0}{dt^2}+y_0=0,\\
O(\varepsilon^1):&\quad \frac{d^2y_1}{dt^2}+y_1=y_0^3
\end{aligned}
$$

が得られます。初期条件 $y(0)=1$, $y'(0)=0$ はすべての $\varepsilon$ で成り立つので、次数ごとに

$$
y_0(0)=1,\quad \frac{dy_0}{dt}(0)=0,\qquad y_1(0)=0,\quad \frac{dy_1}{dt}(0)=0
$$

です。以上が答えです。

### 設問(b)

$O(\varepsilon^0)$ の方程式は単振動で、初期条件から

$$
y_0(t)=\cos t
$$

です。次に $y_0^3=\cos^3 t$ を線形化します。与えられた恒等式の実部をとると $\cos 3\theta=\cos^3\theta-3\cos\theta\sin^2\theta=4\cos^3\theta-3\cos\theta$ なので

$$
\cos^3 t=\frac{3}{4}\cos t+\frac{1}{4}\cos 3t
$$

です。よって $y_1$ の方程式は

$$
\frac{d^2y_1}{dt^2}+y_1=\frac{3}{4}\cos t+\frac{1}{4}\cos 3t
$$

となります。右辺第1項は固有振動数と共鳴する強制項なので、特解に $t$ に比例する項が現れます。特解を $y_p=\alpha\,t\sin t+\beta\cos 3t$ とおいて代入すると、$\dfrac{d^2}{dt^2}(t\sin t)+t\sin t=2\cos t$、$(-9+1)\beta\cos 3t=\dfrac{1}{4}\cos 3t$ より

$$
2\alpha=\frac{3}{4},\quad -8\beta=\frac{1}{4}
\qquad\Longrightarrow\qquad
\alpha=\frac{3}{8},\quad \beta=-\frac{1}{32}
$$

です。一般解 $y_1=C\cos t+D\sin t+\dfrac{3}{8}t\sin t-\dfrac{1}{32}\cos 3t$ に初期条件を課すと、$y_1(0)=C-\dfrac{1}{32}=0$ より $C=\dfrac{1}{32}$、また

$$
\frac{dy_1}{dt}=-C\sin t+D\cos t+\frac{3}{8}(\sin t+t\cos t)+\frac{3}{32}\sin 3t
$$

の $t=0$ での値は $D$ なので $D=0$ です。答えは

$$
y_0(t)=\cos t,\qquad
y_1(t)=\frac{1}{32}\left(\cos t-\cos 3t\right)+\frac{3}{8}\,t\sin t
$$

です。検算として $y_1(0)=\dfrac{1}{32}-\dfrac{1}{32}=0$、$y_1'(0)=0$ を満たし、代入すれば $y_1''+y_1=\dfrac{3}{4}\cos t+\dfrac{1}{4}\cos 3t$ が確かめられます。

### 設問(c)

ここまでの結果から、2次までの近似解は

$$
y(t;\varepsilon)\simeq\cos t+\varepsilon\left[\frac{1}{32}(\cos t-\cos 3t)+\frac{3}{8}\,t\sin t\right]
$$

です。補正項に $t$ に比例して振幅が増大する項（永年項）$\dfrac{3}{8}\varepsilon t\sin t$ が含まれる点が問題です。展開が「$\varepsilon y_1$ は $y_0$ に対する小さな補正」であるためには $\varepsilon t\ll 1$、すなわち $t\ll 1/\varepsilon$ が必要で、$t$ が $1/\varepsilon$ 程度以上に大きくなるとこの展開は近似として破綻します。

一方、もとの方程式の真の解は $t\to\infty$ で発散しません。実際、エネルギー

$$
E=\frac{1}{2}\left(\frac{dy}{dt}\right)^2+\frac{1}{2}y^2-\frac{\varepsilon}{4}y^4
$$

は保存し、初期条件での値は $E=\dfrac{1}{2}-\dfrac{\varepsilon}{4}$ です。ポテンシャル $V(y)=\dfrac{1}{2}y^2-\dfrac{\varepsilon}{4}y^4$ は $y=\pm 1/\sqrt{\varepsilon}$ に高さ $\dfrac{1}{4\varepsilon}$ の障壁を持ち、$\varepsilon$ が微小なら初期エネルギーはこれより十分低いので、解は井戸の中の有界な周期振動にとどまります。永年項は真の解の性質ではなく、展開のとり方の産物です。実際、真の解は振幅に依存して振動数が $\omega=1-\dfrac{3}{8}\varepsilon+O(\varepsilon^2)$ とずれた周期解であり、これを $\varepsilon$ で素朴に展開すると

$$
\cos\omega t=\cos t+\frac{3}{8}\varepsilon t\sin t+O(\varepsilon^2)
$$

となって、まさに設問(b)の永年項が現れます。つまり本来は振動数の補正として現れるべき効果を固定振動数 $\cos t$ のまわりで展開したために、見かけ上の発散項が生じたのです。したがってこの素朴な展開は $t\ll 1/\varepsilon$ の範囲でのみ有効で、長時間の挙動を記述するには不適当です（振動数もあわせて展開するリントシュテット・ポアンカレ法などを使えば一様に有効な近似が得られます）。

## 第2問 3行2列行列の特異値分解

### 設定

3行2列の実行列 $A$ を

$$
A=U\Lambda V^{\mathrm{T}}
$$

の形に分解します（式(1)と呼びます）。$U$, $\Lambda$, $V$ はそれぞれ 3行2列、2行2列、2行2列の実行列で、

$$
V^{\mathrm{T}}V=VV^{\mathrm{T}}=U^{\mathrm{T}}U=\begin{pmatrix}1&0\\0&1\end{pmatrix},\qquad
\Lambda=\begin{pmatrix}\lambda_1&0\\0&\lambda_2\end{pmatrix},\qquad
\lambda_1\ge\lambda_2\ge 0
$$

を仮定します。単位ベクトル $\vec{e}_1=\begin{pmatrix}1\\0\end{pmatrix}$, $\vec{e}_2=\begin{pmatrix}0\\1\end{pmatrix}$ を用いてよい、とあります。

### 設問(a)

式(1)の分解ができたとして、$A^{\mathrm{T}}A$ を計算します。$U^{\mathrm{T}}U=I$（2次単位行列）と $\Lambda^{\mathrm{T}}=\Lambda$ より

$$
A^{\mathrm{T}}A=(U\Lambda V^{\mathrm{T}})^{\mathrm{T}}(U\Lambda V^{\mathrm{T}})
=V\Lambda\,U^{\mathrm{T}}U\,\Lambda V^{\mathrm{T}}
=V\Lambda^2 V^{\mathrm{T}},\qquad
\Lambda^2=\begin{pmatrix}\lambda_1^2&0\\0&\lambda_2^2\end{pmatrix}
$$

です。$V$ の第 $j$ 列は $V\vec{e}_j$ と書けます。$V^{\mathrm{T}}V=I$ を使うと

$$
A^{\mathrm{T}}A\,(V\vec{e}_j)=V\Lambda^2 V^{\mathrm{T}}V\vec{e}_j=V\Lambda^2\vec{e}_j=\lambda_j^2\,(V\vec{e}_j)\qquad(j=1,2)
$$

となります。すなわち $A^{\mathrm{T}}A$ の固有値は $\lambda_1^2$ と $\lambda_2^2$ であり、対応する固有ベクトルはそれぞれ $V\vec{e}_1$（$V$ の第1列）と $V\vec{e}_2$（$V$ の第2列）です。$A^{\mathrm{T}}A$ は2次行列なので固有値はこの2つで尽くされています（$V\vec{e}_1$, $V\vec{e}_2$ は直交する単位ベクトルなので1次独立です）。なお $U$ は $A^{\mathrm{T}}A$ の固有値・固有ベクトルには現れません。

### 設問(b)

与えられた行列は

$$
A=\begin{pmatrix}
\dfrac{3+\sqrt{2}}{2} & \dfrac{3-\sqrt{2}}{2}\\[2mm]
\dfrac{3-\sqrt{2}}{2} & \dfrac{3+\sqrt{2}}{2}\\[2mm]
\dfrac{3\sqrt{2}}{2} & \dfrac{3\sqrt{2}}{2}
\end{pmatrix}
$$

です。$A^{\mathrm{T}}A$ の成分は $A$ の列ベクトル同士の内積です。$p=\dfrac{3+\sqrt{2}}{2}$, $q=\dfrac{3-\sqrt{2}}{2}$, $r=\dfrac{3\sqrt{2}}{2}$ とおくと

$$
p^2+q^2=\frac{(11+6\sqrt{2})+(11-6\sqrt{2})}{4}=\frac{11}{2},\qquad
r^2=\frac{9}{2},\qquad
pq=\frac{9-2}{4}=\frac{7}{4}
$$

なので、対角成分は $p^2+q^2+r^2=\dfrac{11}{2}+\dfrac{9}{2}=10$、非対角成分は $2pq+r^2=\dfrac{7}{2}+\dfrac{9}{2}=8$ となり

$$
A^{\mathrm{T}}A=\begin{pmatrix}10&8\\8&10\end{pmatrix}
$$

です。固有方程式は $(10-\mu)^2-64=0$、すなわち $\mu=10\pm 8$ で、固有値は $\mu_1=18$ と $\mu_2=2$ です。固有ベクトルは、$\mu_1=18$ に対して $(10-18)v_1+8v_2=0$ より $v_1=v_2$、$\mu_2=2$ に対して $8v_1+8v_2=0$ より $v_1=-v_2$ です。長さ1で第1要素が非負という条件から、答えは

$$
\mu_1=18:\ \vec{v}_1=\frac{1}{\sqrt{2}}\begin{pmatrix}1\\1\end{pmatrix},\qquad
\mu_2=2:\ \vec{v}_2=\frac{1}{\sqrt{2}}\begin{pmatrix}1\\-1\end{pmatrix}
$$

です。

### 設問(c)

設問(a)より、$\lambda_1^2$, $\lambda_2^2$ は $A^{\mathrm{T}}A$ の固有値で、$\lambda_1\ge\lambda_2\ge 0$ から

$$
\lambda_1=\sqrt{18}=3\sqrt{2},\qquad \lambda_2=\sqrt{2},\qquad
\Lambda=\begin{pmatrix}3\sqrt{2}&0\\0&\sqrt{2}\end{pmatrix}
$$

です。$V$ の列は対応する固有ベクトルで、$V$ の第1行の各要素が非負という条件は設問(b)の $\vec{v}_1$, $\vec{v}_2$（どちらも第1要素が $1/\sqrt{2}>0$）で満たされるので

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

ととれます（この $V$ は対称かつ直交行列で、$V^{\mathrm{T}}V=VV^{\mathrm{T}}=I$ を満たします）。$U$ は式(1)を $V$ で右から解いて求めます。$\lambda_1,\lambda_2>0$ なので $\Lambda$ は正則で、$A=U\Lambda V^{\mathrm{T}}$ より $U=AV\Lambda^{-1}$、すなわち $U$ の第 $j$ 列は $\vec{u}_j=\dfrac{1}{\lambda_j}A\vec{v}_j$ です。

$$
A\vec{v}_1=\frac{1}{\sqrt{2}}\begin{pmatrix}p+q\\q+p\\2r\end{pmatrix}
=\frac{1}{\sqrt{2}}\begin{pmatrix}3\\3\\3\sqrt{2}\end{pmatrix}
\ \Longrightarrow\
\vec{u}_1=\frac{1}{3\sqrt{2}}\cdot\frac{1}{\sqrt{2}}\begin{pmatrix}3\\3\\3\sqrt{2}\end{pmatrix}
=\begin{pmatrix}1/2\\1/2\\1/\sqrt{2}\end{pmatrix}
$$

$$
A\vec{v}_2=\frac{1}{\sqrt{2}}\begin{pmatrix}p-q\\q-p\\0\end{pmatrix}
=\frac{1}{\sqrt{2}}\begin{pmatrix}\sqrt{2}\\-\sqrt{2}\\0\end{pmatrix}
=\begin{pmatrix}1\\-1\\0\end{pmatrix}
\ \Longrightarrow\
\vec{u}_2=\frac{1}{\sqrt{2}}\begin{pmatrix}1\\-1\\0\end{pmatrix}
$$

よって答えは

$$
U=\begin{pmatrix}
\dfrac{1}{2} & \dfrac{1}{\sqrt{2}}\\[2mm]
\dfrac{1}{2} & -\dfrac{1}{\sqrt{2}}\\[2mm]
\dfrac{1}{\sqrt{2}} & 0
\end{pmatrix},\qquad
\Lambda=\begin{pmatrix}3\sqrt{2}&0\\0&\sqrt{2}\end{pmatrix},\qquad
V=\frac{1}{\sqrt{2}}\begin{pmatrix}1&1\\1&-1\end{pmatrix}
$$

です。検算します。$\vec{u}_1\cdot\vec{u}_1=\dfrac{1}{4}+\dfrac{1}{4}+\dfrac{1}{2}=1$、$\vec{u}_2\cdot\vec{u}_2=\dfrac{1}{2}+\dfrac{1}{2}=1$、$\vec{u}_1\cdot\vec{u}_2=\dfrac{1}{2\sqrt{2}}-\dfrac{1}{2\sqrt{2}}=0$ なので $U^{\mathrm{T}}U=I$ が成り立ちます。さらに $U\Lambda V^{\mathrm{T}}$ の第1行は $\lambda_1 u_{11}\vec{v}_1^{\mathrm{T}}+\lambda_2 u_{12}\vec{v}_2^{\mathrm{T}}=\dfrac{3}{2}(1,\,1)+\dfrac{1}{\sqrt{2}}(1,\,-1)=\left(\dfrac{3+\sqrt{2}}{2},\,\dfrac{3-\sqrt{2}}{2}\right)$ となり、第2行・第3行も同様に $A$ の各行を再現します。

## 第3問 ビュフォンの針（長針の場合）

### 設定

床一面に間隔 $a$ の平行線が引かれていて、そこへ長さ $l\ (\ge a)$ の棒を上から落とします。棒が平行線と交わる確率 $P(l)$ を、まず棒の一端の位置を固定した条件付き確率 $p(x)$ から求め、次に端点の位置について平均する、という2段構えで計算します。棒の太さは無視し、方位角は床面内で一様に分布します。いわゆるビュフォンの針で、通常よく扱われる $l\le a$（短針）ではなく $l\ge a$（長針）の場合です。

### 設問(a)

棒の一端（端点 O と呼びます）が最も近い平行線から距離 $x\ (0\le x\le a/2)$ にあるとします。平行線に垂直な方向の座標をとり、最も近い線を高さ $0$、端点 O を高さ $x$ に置くと、隣の線は高さ $a$ にあります（O が線の下側にある配置は、方位角の分布が上下対称なので同じ結果を与えます）。棒の方位角 $\varphi$ を平行線の方向から測ると、$\varphi$ は $[0,2\pi)$ で一様分布し、もう一方の端点の高さは $x+l\sin\varphi$ です。

棒は線と交わるのは、下の線（高さ $0$）をまたぐ場合

$$
x+l\sin\varphi\le 0
\quad\Longleftrightarrow\quad
\sin\varphi\le-\frac{x}{l}
$$

または上の線（高さ $a$）をまたぐ場合

$$
x+l\sin\varphi\ge a
\quad\Longleftrightarrow\quad
\sin\varphi\ge\frac{a-x}{l}
$$

のいずれかです。$0\le x\le a/2$ と $l\ge a$ から $0\le\dfrac{x}{l}<\dfrac{a-x}{l}\le 1$ なので、2つの事象は排反で、しかもどちらの閾値も $[0,1]$ に収まります（さらに下や上の線をまたげば2本以上と交わりますが、「交わる」事象としては上の2条件で尽くされています）。

一様分布する $\varphi$ に対して、$0\le c\le 1$ のとき $\sin\varphi\ge c$ となる $\varphi$ の範囲は $[\arcsin c,\ \pi-\arcsin c]$ で、その測度は $\pi-2\arcsin c$ です。よって

$$
\mathrm{Prob}(\sin\varphi\ge c)=\frac{\pi-2\arcsin c}{2\pi}=\frac{1}{\pi}\arccos c
$$

です。$\sin\varphi\le-c$ の確率も対称性から同じ値です。したがって

$$
p(x)=\frac{1}{\pi}\left[\arccos\frac{x}{l}+\arccos\frac{a-x}{l}\right]
$$

が答えです。検算として、$x=0$（端点が線上）では $\arccos 0=\dfrac{\pi}{2}$ より $p(0)=\dfrac{1}{2}+\dfrac{1}{\pi}\arccos\dfrac{a}{l}\ge\dfrac{1}{2}$ となり、線上に端点があれば少なくとも下向き半分の方位で交わるという直観と合います。また $l\to\infty$ では両方の $\arccos$ が $\dfrac{\pi}{2}$ に近づき $p(x)\to 1$ で、十分長い棒はほぼ確実に線と交わることに対応します。

### 設問(b)

端点 O の最も近い線からの距離 $x$ は $[0,\ a/2]$ で一様に分布するので、密度は $\dfrac{2}{a}$ です。求める確率は

$$
P(l)=\frac{2}{a}\int_0^{a/2}p(x)\,dx
=\frac{2}{\pi a}\int_0^{a/2}\left[\arccos\frac{x}{l}+\arccos\frac{a-x}{l}\right]dx
$$

です。第2項で $u=a-x$ と置換すると積分区間は $[a/2,\ a]$ に移り、2つの積分がつながって

$$
P(l)=\frac{2}{\pi a}\int_0^{a}\arccos\frac{x}{l}\,dx
$$

となります。原始関数は

$$
\int\arccos\frac{x}{l}\,dx=x\arccos\frac{x}{l}-\sqrt{l^2-x^2}+C
$$

です（微分すると $\arccos\dfrac{x}{l}-\dfrac{x}{\sqrt{l^2-x^2}}+\dfrac{x}{\sqrt{l^2-x^2}}=\arccos\dfrac{x}{l}$ で確かめられます）。$l\ge a$ なので積分区間全体で $\arccos$ の引数は $[0,1]$ に収まり、

$$
\int_0^{a}\arccos\frac{x}{l}\,dx
=\left[x\arccos\frac{x}{l}-\sqrt{l^2-x^2}\right]_0^{a}
=a\arccos\frac{a}{l}+l-\sqrt{l^2-a^2}
$$

です。答えは

$$
P(l)=\frac{2}{\pi a}\left(a\arccos\frac{a}{l}+l-\sqrt{l^2-a^2}\right)
=\frac{2}{\pi}\arccos\frac{a}{l}+\frac{2}{\pi a}\left(l-\sqrt{l^2-a^2}\right)
$$

です。検算します。$l=a$ では $\arccos 1=0$、$\sqrt{l^2-a^2}=0$ より $P(a)=\dfrac{2}{\pi}$ となり、短針の公式 $P=\dfrac{2l}{\pi a}$ を $l=a$ で評価した値と連続につながります。$l\to\infty$ では $\arccos\dfrac{a}{l}\to\dfrac{\pi}{2}$、$l-\sqrt{l^2-a^2}=\dfrac{a^2}{l+\sqrt{l^2-a^2}}\to 0$ なので $P(l)\to 1$ となり、十分長い棒は確率1で線と交わるという当然の極限を再現します。また $P(l)$ は $p(x)\le 1$ の平均なので $1$ を超えることはなく、$l$ について単調増加で $\dfrac{2}{\pi}$ から $1$ へ向かいます。

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