コンテンツにスキップ

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

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

生 Markdown

この年度は一般教育科目として数学と物理が 1 冊の問題冊子にまとめられ、数学 3 問・物理 3 問の合計 6 問すべてが必答でした(平成10年8月25日、9時00分〜12時00分)。本記事は数学の 3 問を扱います。3 問とも計算量は小さく、代わりに「示せ」「証明せよ」が全体の半分以上を占めるため、論理を閉じる書き方ができるかどうかで点差がつく構成です。

問題分野主題
第1問複素解析・フーリエ解析ディリクレ積分と絶対値関数のフーリエ級数
第2問線形代数特性多項式の係数と複素固有値をもつ条件
第3問特殊関数重み e−x2e^{-x^2} に関する直交多項式

第1問では α\alpha の符号による場合分けと、留数定理をそのまま使えない(極が積分路上にある)ことへの対処が要点です。第2問は 3 次方程式の判別式を自力で作る問題に帰着します。第3問は「nn 次多項式は HℓH_\ell で展開できる」という問題文中の注意が 3 つの設問すべての鍵になっており、これを補題として最初に立てておくと証明が短くなります。

第1問 ディリクレ積分と絶対値関数のフーリエ級数

Section titled “第1問 ディリクレ積分と絶対値関数のフーリエ級数”

設問(1)は実数 α\alpha に対する

I(α)=∫0∞dx sin⁡αxxI(\alpha)=\int_0^{\infty} dx\,\frac{\sin\alpha x}{x}

を複素積分で求める問題です。設問(2)は区間 [−π,π][-\pi,\pi] 上の f(x)=∣x∣f(x)=|x| を

f(x)=a02+∑n=1∞(ancos⁡nx+bnsin⁡nx)f(x)=\frac{a_0}{2}+\sum_{n=1}^{\infty}\left(a_n\cos nx+b_n\sin nx\right)

と展開したときの係数を求め、その結果から奇数の平方の逆数和が π2/8\pi^2/8 になることを示す問題です。

まず積分が存在することを確認します。被積分関数は x→0x\to 0 で sin⁡αx/x→α\sin\alpha x/x\to\alpha と有界なので原点は特異点ではありません。無限遠側は、α≠0\alpha\neq 0 として部分積分すれば

∫1Rsin⁡αxx dx=[−cos⁡αxαx]1R−1α∫1Rcos⁡αxx2 dx\int_1^{R}\frac{\sin\alpha x}{x}\,dx=\left[-\frac{\cos\alpha x}{\alpha x}\right]_1^{R}-\frac{1}{\alpha}\int_1^{R}\frac{\cos\alpha x}{x^2}\,dx

となり、右辺第 2 項は ∣cos⁡αx/x2∣≤x−2|\cos\alpha x/x^2|\le x^{-2} より R→∞R\to\infty で絶対収束します。よって I(α)I(\alpha) は広義リーマン積分(上端を RR で打ち切った積分の R→∞R\to\infty の極限として定める積分)として収束します(∣sin⁡αx∣/x|\sin\alpha x|/x 自体は可積分でないので、絶対収束ではありません)。

α↦−α\alpha\mapsto-\alpha で被積分関数は符号を変えるので I(−α)=−I(α)I(-\alpha)=-I(\alpha) であり、I(0)=0I(0)=0 です。したがって α>0\alpha\gt 0 の場合を計算すれば十分です。

以下 α>0\alpha\gt 0 とし、C∖{0}\mathbb{C}\setminus\{0\} で正則な

g(z)=eiαzzg(z)=\frac{e^{i\alpha z}}{z}

を、次の閉曲線に沿って積分します。正の実軸上の線分 [ε,R][\varepsilon,R]、上半平面の大半円 ΓR:z=Reiθ (θ:0→π)\Gamma_R:z=Re^{i\theta}\ (\theta:0\to\pi)、負の実軸上の線分 [−R,−ε][-R,-\varepsilon]、そして原点の上側を回る小半円 γε:z=εeiθ (θ:π→0)\gamma_\varepsilon:z=\varepsilon e^{i\theta}\ (\theta:\pi\to 0) をこの順につないだ経路です(0<ε<R0\lt\varepsilon\lt R)。原点はこの閉曲線の内部から小半円によって除かれているので、囲まれた領域で gg は正則です。コーシーの積分定理より

∫εReiαxxdx+∫ΓRg(z)dz+∫−R−εeiαxxdx+∫γεg(z)dz=0\int_{\varepsilon}^{R}\frac{e^{i\alpha x}}{x}dx+\int_{\Gamma_R}g(z)dz+\int_{-R}^{-\varepsilon}\frac{e^{i\alpha x}}{x}dx+\int_{\gamma_\varepsilon}g(z)dz=0

が成り立ちます。

大半円の寄与を評価します。z=Reiθz=Re^{i\theta} では dz/z=i dθdz/z=i\,d\theta かつ ∣eiαz∣=e−αRsin⁡θ|e^{i\alpha z}|=e^{-\alpha R\sin\theta} なので、[0,π/2][0,\pi/2] で成り立つジョルダンの不等式 sin⁡θ≥2θ/π\sin\theta\ge 2\theta/\pi(sin⁡\sin の凹性から従います)を用いて

∣∫ΓRg(z)dz∣=∣∫0πeiαReiθ i dθ∣≤∫0πe−αRsin⁡θdθ=2∫0π/2e−αRsin⁡θdθ≤2∫0π/2e−2αRθ/πdθ=παR(1−e−αR)≤παR\begin{aligned} \left\lvert\int_{\Gamma_R}g(z)dz\right\rvert &=\left\lvert\int_0^{\pi}e^{i\alpha Re^{i\theta}}\,i\,d\theta\right\rvert \le\int_0^{\pi}e^{-\alpha R\sin\theta}d\theta =2\int_0^{\pi/2}e^{-\alpha R\sin\theta}d\theta\\ &\le 2\int_0^{\pi/2}e^{-2\alpha R\theta/\pi}d\theta =\frac{\pi}{\alpha R}\left(1-e^{-\alpha R}\right)\le\frac{\pi}{\alpha R} \end{aligned}

と押さえられます。α>0\alpha\gt 0 なので R→∞R\to\infty でこれは 00 に収束します。ここで α\alpha の符号が本質的に効いていて、α<0\alpha\lt 0 なら上半平面では被積分関数が発散し、この評価は使えません。

小半円の寄与は、z=εeiθz=\varepsilon e^{i\theta} で dz/z=i dθdz/z=i\,d\theta、向きは θ:π→0\theta:\pi\to 0 なので

∫γεg(z)dz=i∫π0eiαεeiθdθ=−iπ+i∫π0(eiαεeiθ−1)dθ\int_{\gamma_\varepsilon}g(z)dz=i\int_{\pi}^{0}e^{i\alpha\varepsilon e^{i\theta}}d\theta =-i\pi+i\int_{\pi}^{0}\left(e^{i\alpha\varepsilon e^{i\theta}}-1\right)d\theta

であり、∣w∣≤αε|w|\le\alpha\varepsilon に対する ∣eiw−1∣≤∣w∣e∣w∣|e^{iw}-1|\le|w|e^{|w|} から残差は παεeαε\pi\alpha\varepsilon e^{\alpha\varepsilon} 以下です。よって ε→0\varepsilon\to 0 で ∫γεg dz→−iπ\int_{\gamma_\varepsilon}g\,dz\to-i\pi となります。1/z1/z の原点まわりの半周が全周の半分の 2πi/22\pi i/2 に符号を付けたものになっている、という標準的な結果です。

実軸上の 2 つの線分は 1 つにまとめられます。負側で x=−ux=-u と置換すると

∫−R−εeiαxxdx=−∫εRe−iαuudu\int_{-R}^{-\varepsilon}\frac{e^{i\alpha x}}{x}dx=-\int_{\varepsilon}^{R}\frac{e^{-i\alpha u}}{u}du

なので、2 つの和は

∫εReiαx−e−iαxxdx=2i∫εRsin⁡αxxdx\int_{\varepsilon}^{R}\frac{e^{i\alpha x}-e^{-i\alpha x}}{x}dx=2i\int_{\varepsilon}^{R}\frac{\sin\alpha x}{x}dx

です。以上を積分定理の式に戻し、R→∞R\to\infty、ε→0\varepsilon\to 0 とすると

2i∫0∞sin⁡αxxdx=iπ2i\int_0^{\infty}\frac{\sin\alpha x}{x}dx=i\pi

を得ます。したがって α>0\alpha\gt 0 で I(α)=π/2I(\alpha)=\pi/2 です。奇関数性と I(0)=0I(0)=0 を合わせて、答えは

∫0∞dx sin⁡αxx=π2 sgn α={  π/2(α>0)  0(α=0)−π/2(α<0)\int_0^{\infty}dx\,\frac{\sin\alpha x}{x}=\frac{\pi}{2}\,\mathrm{sgn}\,\alpha= \begin{cases} \ \ \pi/2 & (\alpha\gt 0)\\ \ \ 0 & (\alpha=0)\\ -\pi/2 & (\alpha\lt 0) \end{cases}

です。値が ∣α∣|\alpha| に依存しないことは、α>0\alpha\gt 0 で y=αxy=\alpha x と置換すると α\alpha が完全に消えることから独立に確認できます。なお実部を取れば P ⁣∫−∞∞cos⁡(αx)/x dx=0\mathrm{P}\!\int_{-\infty}^{\infty}\cos(\alpha x)/x\,dx=0 で、これは被積分関数が奇関数であることと整合します。

f(x)=∣x∣f(x)=|x| は偶関数です。f(x)sin⁡nxf(x)\sin nx は奇関数で、対称区間での積分は 00 になるので

bn=1π∫−ππ∣x∣sin⁡nx dx=0(n=1,2,⋯ )b_n=\frac{1}{\pi}\int_{-\pi}^{\pi}|x|\sin nx\,dx=0\qquad(n=1,2,\cdots)

です。余弦係数は偶関数性から

an=1π∫−ππ∣x∣cos⁡nx dx=2π∫0πxcos⁡nx dxa_n=\frac{1}{\pi}\int_{-\pi}^{\pi}|x|\cos nx\,dx=\frac{2}{\pi}\int_0^{\pi}x\cos nx\,dx

となります。n=0n=0 では

a0=2π∫0πx dx=2π⋅π22=πa_0=\frac{2}{\pi}\int_0^{\pi}x\,dx=\frac{2}{\pi}\cdot\frac{\pi^2}{2}=\pi

です。n≥1n\ge 1 では部分積分して

∫0πxcos⁡nx dx=[xsin⁡nxn]0π−1n∫0πsin⁡nx dx=0+[cos⁡nxn2]0π=(−1)n−1n2\int_0^{\pi}x\cos nx\,dx=\left[\frac{x\sin nx}{n}\right]_0^{\pi}-\frac{1}{n}\int_0^{\pi}\sin nx\,dx =0+\left[\frac{\cos nx}{n^2}\right]_0^{\pi}=\frac{(-1)^n-1}{n^2}

となるので

an=2π⋅(−1)n−1n2={0(n even)−4πn2(n odd)a_n=\frac{2}{\pi}\cdot\frac{(-1)^n-1}{n^2}= \begin{cases} 0 & (n\ \text{even})\\[2pt] -\dfrac{4}{\pi n^2} & (n\ \text{odd}) \end{cases}

です。答えは a0=πa_0=\pi、偶数 n≥2n\ge 2 で an=0a_n=0、奇数 nn で an=−4/(πn2)a_n=-4/(\pi n^2)、そしてすべての nn で bn=0b_n=0 です。フーリエ級数は

∣x∣=π2−4π∑k=0∞cos⁡(2k+1)x(2k+1)2(x∈[−π,π])|x|=\frac{\pi}{2}-\frac{4}{\pi}\sum_{k=0}^{\infty}\frac{\cos(2k+1)x}{(2k+1)^2} \qquad(x\in[-\pi,\pi])

の形になります。偶数次の余弦が落ちる理由は、∣x∣−π/2|x|-\pi/2 が x→π−xx\to\pi-x の下で符号を変える(半区間ごとの反対称性をもつ)ことに対応しています。

まず級数が x=0x=0 で f(0)=0f(0)=0 に収束することを確かめます。係数は ∑k4/{π(2k+1)2}<∞\sum_{k}4/\{\pi(2k+1)^2\}\lt\infty と絶対総和可能なので、ワイエルシュトラスの M 判定法により右辺の級数は R\mathbb{R} 上一様収束し、その和 g(x)g(x) は連続関数です。他方 ff を 2π2\pi 周期に拡張したものは f(−π)=f(π)=πf(-\pi)=f(\pi)=\pi より連続で、フェイェールの定理(連続な周期関数のフーリエ級数は、部分和の平均であるチェザロ平均をとると各点でもとの関数に収束する、という定理)からそのフーリエ級数のチェザロ平均は各点で ff に収束します。収束する級数のチェザロ平均は同じ和に収束するので g=fg=f です。区分的に滑らかな連続関数に対するディリクレの定理を使って直接 g=fg=f と結論しても同じです(f′(x)=sgn xf'(x)=\mathrm{sgn}\,x は x=0,±πx=0,\pm\pi を除いて連続)。

そこで上で得た展開に x=0x=0 を代入します。f(0)=0f(0)=0、cos⁡0=1\cos 0=1 なので

0=π2−4π∑k=0∞1(2k+1)2,すなわち∑k=0∞1(2k+1)2=π4⋅π2=π280=\frac{\pi}{2}-\frac{4}{\pi}\sum_{k=0}^{\infty}\frac{1}{(2k+1)^2}, \qquad\text{すなわち}\qquad \sum_{k=0}^{\infty}\frac{1}{(2k+1)^2}=\frac{\pi}{4}\cdot\frac{\pi}{2}=\frac{\pi^2}{8}

が従います。これが示すべき式です。

検算を 3 つ挙げます。第 1 に、既知の ζ(2)=π2/6\zeta(2)=\pi^2/6(ζ(s)=∑n≥1n−s\zeta(s)=\sum_{n\ge1}n^{-s} はゼータ関数)を使うと、奇数項の和は ∑n oddn−2=ζ(2)−∑m≥1(2m)−2=ζ(2)(1−1/4)=π2/8\sum_{n\ \mathrm{odd}}n^{-2}=\zeta(2)-\sum_{m\ge1}(2m)^{-2}=\zeta(2)(1-1/4)=\pi^2/8 で一致します。第 2 に、同じ展開に x=πx=\pi を入れると右辺は π/2+(4/π)(π2/8)=π\pi/2+(4/\pi)(\pi^2/8)=\pi となり f(π)=πf(\pi)=\pi に一致します。第 3 に、パーセバルの等式

1π∫−ππf(x)2dx=a022+∑n=1∞(an2+bn2)\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)^2dx=\frac{a_0^2}{2}+\sum_{n=1}^{\infty}\left(a_n^2+b_n^2\right)

の左辺は 2π2/32\pi^2/3、右辺は π2/2+(16/π2)∑k(2k+1)−4\pi^2/2+(16/\pi^2)\sum_k(2k+1)^{-4} なので ∑k(2k+1)−4=π4/96\sum_k(2k+1)^{-4}=\pi^4/96 となり、ζ(4)=π4/90\zeta(4)=\pi^4/90 から得られる ζ(4)(1−1/16)=π4/96\zeta(4)(1-1/16)=\pi^4/96 と一致します。

第2問 特性多項式の係数と複素固有値をもつ条件

Section titled “第2問 特性多項式の係数と複素固有値をもつ条件”

33 次の実行列 A={aij}A=\{a_{ij}\} の特性多項式を

F(λ)≡det⁡[A−λI]=−λ3+Pλ2+Qλ+RF(\lambda)\equiv\det[A-\lambda I]=-\lambda^3+P\lambda^2+Q\lambda+R

と置きます(II は単位行列)。設問(1)は P=Tr⁡AP=\operatorname{Tr}A、Q=12(Tr⁡A2−(Tr⁡A)2)Q=\frac{1}{2}\left(\operatorname{Tr}A^2-(\operatorname{Tr}A)^2\right)、R=det⁡AR=\det A を示す問題、設問(2)は P=0P=0 のときに AA が複素共役な(実でない)固有値の対をもつための条件を Q,RQ,R で表す問題です。

RR は最も簡単で、λ=0\lambda=0 を代入すれば

R=F(0)=det⁡[A−0⋅I]=det⁡AR=F(0)=\det[A-0\cdot I]=\det A

です。

PP と QQ は行列式の置換展開

det⁡[A−λI]=∑σ∈S3sgn(σ)∏i=13(A−λI)iσ(i)\det[A-\lambda I]=\sum_{\sigma\in S_3}\mathrm{sgn}(\sigma)\prod_{i=1}^{3}(A-\lambda I)_{i\sigma(i)}

から読み取ります。λ\lambda は対角成分 (A−λI)ii=aii−λ(A-\lambda I)_{ii}=a_{ii}-\lambda にしか現れません。恒等置換でない σ∈S3\sigma\in S_3 は少なくとも 2 つの添字を動かすので、固定点は高々 1 つ、すなわちその項は λ\lambda について高々 1 次です。したがって λ3\lambda^3 と λ2\lambda^2 の係数は恒等置換の項

(a11−λ)(a22−λ)(a33−λ)=−λ3+(a11+a22+a33)λ2−(a11a22+a22a33+a33a11)λ+a11a22a33(a_{11}-\lambda)(a_{22}-\lambda)(a_{33}-\lambda) =-\lambda^3+(a_{11}+a_{22}+a_{33})\lambda^2-(a_{11}a_{22}+a_{22}a_{33}+a_{33}a_{11})\lambda+a_{11}a_{22}a_{33}

だけから来ます。λ2\lambda^2 の係数を比較して

P=a11+a22+a33=Tr⁡AP=a_{11}+a_{22}+a_{33}=\operatorname{Tr}A

です。

λ1\lambda^1 の係数には、恒等置換の項に加えて、固定点をちょうど 1 つもつ置換、すなわち 3 つの互換の項が寄与します(3-サイクルは固定点をもたず λ\lambda を含みません)。たとえば互換 σ=(2 3)\sigma=(2\,3) の項は符号 −1-1 で −(a11−λ)a23a32-(a_{11}-\lambda)a_{23}a_{32} となり、λ\lambda の係数として +a23a32+a_{23}a_{32} を与えます。他の 2 つの互換も同様なので

Q=(a12a21+a23a32+a31a13)−(a11a22+a22a33+a33a11)Q=\left(a_{12}a_{21}+a_{23}a_{32}+a_{31}a_{13}\right)-\left(a_{11}a_{22}+a_{22}a_{33}+a_{33}a_{11}\right)

です。一方、(A2)ii=∑jaijaji(A^2)_{ii}=\sum_j a_{ij}a_{ji} より

Tr⁡A2=∑i,jaijaji=∑iaii2+2(a12a21+a23a32+a31a13)\operatorname{Tr}A^2=\sum_{i,j}a_{ij}a_{ji}=\sum_i a_{ii}^2+2\left(a_{12}a_{21}+a_{23}a_{32}+a_{31}a_{13}\right)

であり、また

(Tr⁡A)2=∑iaii2+2(a11a22+a22a33+a33a11)(\operatorname{Tr}A)^2=\sum_i a_{ii}^2+2\left(a_{11}a_{22}+a_{22}a_{33}+a_{33}a_{11}\right)

なので、差の半分をとると

12(Tr⁡A2−(Tr⁡A)2)=(a12a21+a23a32+a31a13)−(a11a22+a22a33+a33a11)=Q\frac{1}{2}\left(\operatorname{Tr}A^2-(\operatorname{Tr}A)^2\right) =\left(a_{12}a_{21}+a_{23}a_{32}+a_{31}a_{13}\right)-\left(a_{11}a_{22}+a_{22}a_{33}+a_{33}a_{11}\right)=Q

となり、3 式すべてが示されました。

検算として固有値表示と突き合わせます。固有値を λ1,λ2,λ3\lambda_1,\lambda_2,\lambda_3(重複を込めて、複素数の範囲で)とすると F(λ)=−(λ−λ1)(λ−λ2)(λ−λ3)F(\lambda)=-(\lambda-\lambda_1)(\lambda-\lambda_2)(\lambda-\lambda_3) なので P=∑iλiP=\sum_i\lambda_i、Q=−∑i<jλiλjQ=-\sum_{i\lt j}\lambda_i\lambda_j、R=λ1λ2λ3R=\lambda_1\lambda_2\lambda_3 です。Tr⁡A=∑iλi\operatorname{Tr}A=\sum_i\lambda_i、Tr⁡A2=∑iλi2=(∑iλi)2−2∑i<jλiλj\operatorname{Tr}A^2=\sum_i\lambda_i^2=(\sum_i\lambda_i)^2-2\sum_{i\lt j}\lambda_i\lambda_j を使うと上の 3 式が再現され、これはニュートンの恒等式(べき和と基本対称式を結ぶ関係式)の p2=e12−2e2p_2=e_1^2-2e_2 にほかなりません。

P=0P=0 のとき、固有値は

g(λ)≡−F(λ)=λ3−Qλ−R=0g(\lambda)\equiv-F(\lambda)=\lambda^3-Q\lambda-R=0

の根です。gg は実係数なので、実でない根は必ず複素共役の対で現れます。3 次ですから、根の配置は「3 個とも実」か「実根 1 個+複素共役対」の二択であり、複素共役の固有値をもつことは「gg の根がすべて実、ではない」ことと同値です。これを gg の増減から判定します。g′(λ)=3λ2−Qg'(\lambda)=3\lambda^2-Q です。

Q>0Q\gt 0 の場合。gg は λ=±m\lambda=\pm m(m≡Q/3m\equiv\sqrt{Q/3})に極値をもち、λ=−m\lambda=-m が極大、λ=+m\lambda=+m が極小です。3 次関数のグラフを考えると、根が 3 個とも実であることは、極大値が 00 以上かつ極小値が 00 以下であること、すなわち g(−m)≥0≥g(m)g(-m)\ge 0\ge g(m) と同値です(等号のときは極値点が重根になります)。極値は m3=Q3mm^3=\frac{Q}{3}m を使って

g(±m)=±m3∓Qm−R=∓2Q3m−Rg(\pm m)=\pm m^3\mp Qm-R=\mp\frac{2Q}{3}m-R

と書け、つねに g(−m)−g(m)=4Q3m>0g(-m)-g(m)=\frac{4Q}{3}m\gt 0 なので、「3 実根でない」ことは 2 つの極値が同符号であること、つまり

g(m) g(−m)=R2−4Q29m2=R2−4Q327>0g(m)\,g(-m)=R^2-\frac{4Q^2}{9}m^2=R^2-\frac{4Q^3}{27}\gt 0

と同値です。

Q<0Q\lt 0 の場合。g′(λ)=3λ2−Q>0g'(\lambda)=3\lambda^2-Q\gt 0 なので gg は狭義単調増加で、実根はちょうど 1 個、しかも単根です(重根なら g′=0g'=0 も必要ですが g′>0g'\gt 0)。よって残り 2 根は必ず複素共役対です。このとき 27R2≥0>4Q327R^2\ge 0\gt 4Q^3 なので、上の不等式 R2>4Q3/27R^2\gt 4Q^3/27 も自動的に成り立っています。

Q=0Q=0 の場合。g(λ)=λ3−Rg(\lambda)=\lambda^3-R です。R=0R=0 なら λ=0\lambda=0 の三重根ですべて実、R≠0R\neq 0 なら実根は λ=R1/3\lambda=R^{1/3} の 1 個だけで(狭義単調増加)、複素共役対をもちます。これも 27R2>0=4Q327R^2\gt 0=4Q^3 という判定と一致します。

以上を合わせると、どの場合も同じ 1 本の不等式にまとまります。答えは

27R2>4Q3(同値に 4Q3−27R2<0)27R^2\gt 4Q^3\qquad\left(\text{同値に}\ 4Q^3-27R^2\lt 0\right)

のとき、かつそのときに限り AA は複素共役の固有値をもちます。左辺 4Q3−27R24Q^3-27R^2 は λ3−Qλ−R\lambda^3-Q\lambda-R の判別式 ∏i<j(λi−λj)2\prod_{i\lt j}(\lambda_i-\lambda_j)^2 にほかならず、根が a, u±iva,\ u\pm iv(v≠0v\neq 0)なら (λ1−λ2)(λ1−λ3)=(a−u)2+v2>0(\lambda_1-\lambda_2)(\lambda_1-\lambda_3)=(a-u)^2+v^2\gt 0、(λ2−λ3)2=−4v2<0(\lambda_2-\lambda_3)^2=-4v^2\lt 0 で判別式は負、全根実なら非負、という一般論とも整合します。

具体例で確かめます。固有値 ±i, 0\pm i,\ 0 をもつ実行列(回転生成子に第 3 行・列として零を足したもの)では Tr⁡A=0\operatorname{Tr}A=0、Tr⁡A2=−2\operatorname{Tr}A^2=-2 より Q=−1Q=-1、R=0R=0 で、27R2=0>−4=4Q327R^2=0\gt-4=4Q^3 となり複素固有値ありと正しく判定されます。A=diag(1,−1,0)A=\mathrm{diag}(1,-1,0) では Q=1Q=1、R=0R=0 で 0>40\gt 4 は偽、実際すべて実固有値です。境界の例として固有値 1,1,−21,1,-2 では Q=3Q=3、R=−2R=-2 で 27R2=108=4Q327R^2=108=4Q^3 と等号になり、重根(すべて実)の場合が等号に対応することが確認できます。

第3問 重み e−x2e^{-x^2} に関する直交多項式

Section titled “第3問 重み e−x2e^{-x^2}e−x2 に関する直交多項式”

Hn(x)H_n(x)(n=0,1,2,…n=0,1,2,\ldots)は最高次係数 1 の nn 次多項式で、

Hn(x)=xn+(n より低次の多項式),∫−∞∞dx e−x2Hn(x)Hm(x)=0(n≠m)H_n(x)=x^n+(n\ \text{より低次の多項式}),\qquad \int_{-\infty}^{\infty}dx\,e^{-x^2}H_n(x)H_m(x)=0\quad(n\neq m)

を満たすものとします。以下、内積を ⟨f,g⟩=∫−∞∞e−x2f(x)g(x) dx\langle f,g\rangle=\int_{-\infty}^{\infty}e^{-x^2}f(x)g(x)\,dx、ノルムの 2 乗を hn=⟨Hn,Hn⟩h_n=\langle H_n,H_n\rangle と書きます。重み e−x2e^{-x^2} は正で Hn≢0H_n\not\equiv 0 なので hn>0h_n\gt 0 です。与えられた積分公式はモーメント μk=∫−∞∞xke−x2dx\mu_k=\int_{-\infty}^{\infty}x^k e^{-x^2}dx の言葉で

μ2m=π (2m)!4mm!,μ2m+1=0\mu_{2m}=\sqrt{\pi}\,\frac{(2m)!}{4^m m!},\qquad \mu_{2m+1}=0

(奇数次は被積分関数が奇関数なので 00)となり、具体的には μ0=π\mu_0=\sqrt{\pi}、μ2=12π\mu_2=\frac{1}{2}\sqrt{\pi}、μ4=34π\mu_4=\frac{3}{4}\sqrt{\pi} です。

問題文の注意「nn 次多項式 P(x)P(x) は P=∑ℓ=0ncℓHℓP=\sum_{\ell=0}^{n}c_\ell H_\ell と展開できる」を、次の形の補題として全設問で使います。H0,…,HnH_0,\ldots,H_n は次数が互いに異なるので一次独立で、nn 次以下の多項式の空間(次元 n+1n+1)の基底をなし、展開係数は直交性から cℓ=⟨P,Hℓ⟩/hℓc_\ell=\langle P,H_\ell\rangle/h_\ell と一意に定まります。ここから 2 つの系が出ます。第 1 に、HnH_n は次数 n−1n-1 以下の任意の多項式と直交します(相手を H0,…,Hn−1H_0,\ldots,H_{n-1} で展開すればよい)。第 2 に、最高次係数 1 の nn 次多項式で H0,…,Hn−1H_0,\ldots,H_{n-1} のすべてと直交するものは HnH_n に限ります。実際そのような KK に対し K−HnK-H_n は次数 n−1n-1 以下で全 HℓH_\ell(ℓ<n\ell\lt n)と直交するので、展開係数がすべて 00、つまり K=HnK=H_n です。この一意性を以下「補題(B)」と呼びます。

H0=1H_0=1 です(最高次係数 1 の 0 次式)。Hn=xn+(低次)H_n=x^n+(\text{低次}) と置き、1,x,…,xn−11,x,\ldots,x^{n-1} との直交(上の系より H0,…,Hn−1H_0,\ldots,H_{n-1} との直交と同値)を課します。

n=1n=1: H1=x+cH_1=x+c に対し ⟨H1,1⟩=μ1+cμ0=cπ=0\langle H_1,1\rangle=\mu_1+c\mu_0=c\sqrt{\pi}=0 より c=0c=0 で、

H1(x)=xH_1(x)=x

です。

n=2n=2: H2=x2+bx+cH_2=x^2+bx+c に対し ⟨H2,1⟩=μ2+cμ0=π(12+c)=0\langle H_2,1\rangle=\mu_2+c\mu_0=\sqrt{\pi}\left(\tfrac{1}{2}+c\right)=0 より c=−12c=-\tfrac{1}{2}、⟨H2,x⟩=μ3+bμ2+cμ1=bμ2=0\langle H_2,x\rangle=\mu_3+b\mu_2+c\mu_1=b\mu_2=0 より b=0b=0 で、

H2(x)=x2−12H_2(x)=x^2-\frac{1}{2}

です。

n=3n=3: H3=x3+ax2+bx+cH_3=x^3+ax^2+bx+c に対し、⟨H3,x⟩=μ4+bμ2=π(34+b2)=0\langle H_3,x\rangle=\mu_4+b\mu_2=\sqrt{\pi}\left(\tfrac{3}{4}+\tfrac{b}{2}\right)=0 より b=−32b=-\tfrac{3}{2} です。残る 2 条件は

⟨H3,1⟩=aμ2+cμ0=0,⟨H3,x2⟩=aμ4+cμ2=0\langle H_3,1\rangle=a\mu_2+c\mu_0=0,\qquad \langle H_3,x^2\rangle=a\mu_4+c\mu_2=0

で(奇数次モーメントはすべて落としました)、この (a,c)(a,c) の連立 1 次方程式の係数行列式は μ22−μ0μ4=π4−3π4=−π2≠0\mu_2^2-\mu_0\mu_4=\tfrac{\pi}{4}-\tfrac{3\pi}{4}=-\tfrac{\pi}{2}\neq 0 なので a=c=0a=c=0 です。よって

H3(x)=x3−32xH_3(x)=x^3-\frac{3}{2}x

です。検算として、物理で標準のエルミート多項式 H2phys=4x2−2H_2^{\mathrm{phys}}=4x^2-2、H3phys=8x3−12xH_3^{\mathrm{phys}}=8x^3-12x を最高次係数で割ると x2−12x^2-\tfrac{1}{2}、x3−32xx^3-\tfrac{3}{2}x となり一致します。また ⟨H1,H3⟩=μ4−32μ2=(34−34)π=0\langle H_1,H_3\rangle=\mu_4-\tfrac{3}{2}\mu_2=\left(\tfrac{3}{4}-\tfrac{3}{4}\right)\sqrt{\pi}=0 が直接確認できます。

Kn(x)≡(−1)nHn(−x)K_n(x)\equiv(-1)^n H_n(-x) と置き、すべての nn で Kn=HnK_n=H_n となることを nn についての帰納法で示します。これが示されれば Hn(−x)−(−1)nHn(x)=(−1)n{Kn(x)−Hn(x)}=0H_n(-x)-(-1)^nH_n(x)=(-1)^n\{K_n(x)-H_n(x)\}=0 です。

まず KnK_n の最高次項は (−1)n(−x)n=xn(-1)^n(-x)^n=x^n なので、KnK_n はつねに最高次係数 1 の nn 次多項式です。n=0n=0 では K0=H0=1K_0=H_0=1 で成立します(設問(1)より K1=−H1(−x)=x=H1K_1=-H_1(-x)=x=H_1 も直接確かめられます)。

m<nm\lt n なるすべての mm で Km=HmK_m=H_m、すなわち Hm(−x)=(−1)mHm(x)H_m(-x)=(-1)^mH_m(x) が成り立つと仮定します。m<nm\lt n に対し、置換 x=−tx=-t(dx=−dtdx=-dt、積分範囲は ∞→−∞\infty\to-\infty が反転して元に戻り、e−x2e^{-x^2} は偶関数)を使うと

⟨Kn,Hm⟩=(−1)n∫−∞∞e−x2Hn(−x)Hm(x) dx=(−1)n∫−∞∞e−t2Hn(t)Hm(−t) dt=(−1)n+m∫−∞∞e−t2Hn(t)Hm(t) dt=(−1)n+m⟨Hn,Hm⟩=0\begin{aligned} \langle K_n,H_m\rangle &=(-1)^n\int_{-\infty}^{\infty}e^{-x^2}H_n(-x)H_m(x)\,dx =(-1)^n\int_{-\infty}^{\infty}e^{-t^2}H_n(t)H_m(-t)\,dt\\ &=(-1)^{n+m}\int_{-\infty}^{\infty}e^{-t^2}H_n(t)H_m(t)\,dt =(-1)^{n+m}\langle H_n,H_m\rangle=0 \end{aligned}

となります(2 行目で帰納法の仮定を使いました)。つまり KnK_n は最高次係数 1 の nn 次多項式で H0,…,Hn−1H_0,\ldots,H_{n-1} のすべてと直交するので、補題(B)により Kn=HnK_n=H_n です。帰納法が閉じて、すべての nn で

Hn(−x)=(−1)nHn(x)H_n(-x)=(-1)^nH_n(x)

が示されました。設問(1)の H2=x2−12H_2=x^2-\tfrac{1}{2}(偶)、H3=x3−32xH_3=x^3-\tfrac{3}{2}x(奇)と整合します。

以下 n≥1n\ge 1 とします(Hn−1H_{n-1} が現れるため)。また境界項について、多項式 pp に対し e−x2p(x)→0e^{-x^2}p(x)\to 0(x→±∞x\to\pm\infty)となることを部分積分で繰り返し使います。

第 1 式。D(x)≡xHn(x)−Hn+1(x)D(x)\equiv xH_n(x)-H_{n+1}(x) は、最高次 xn+1x^{n+1} が打ち消し合うので次数 nn 以下の多項式です。補題により D=∑ℓ=0ncℓHℓD=\sum_{\ell=0}^{n}c_\ell H_\ell、cℓhℓ=⟨D,Hℓ⟩c_\ell h_\ell=\langle D,H_\ell\rangle と展開できます。ℓ≤n\ell\le n では ⟨Hn+1,Hℓ⟩=0\langle H_{n+1},H_\ell\rangle=0 なので cℓhℓ=⟨xHn,Hℓ⟩=⟨Hn,xHℓ⟩c_\ell h_\ell=\langle xH_n,H_\ell\rangle=\langle H_n,xH_\ell\rangle です(重み付き内積の中で xx を移しただけ)。ℓ≤n−2\ell\le n-2 なら xHℓxH_\ell は次数 ℓ+1≤n−1\ell+1\le n-1 の多項式なので HnH_n と直交し、cℓ=0c_\ell=0 です。ℓ=n\ell=n では

⟨xHn,Hn⟩=∫−∞∞e−x2x Hn(x)2 dx=0\langle xH_n,H_n\rangle=\int_{-\infty}^{\infty}e^{-x^2}x\,H_n(x)^2\,dx=0

です。設問(2)より Hn(x)2H_n(x)^2 は偶関数で、xx 倍した被積分関数が奇関数になるからです。残るのは ℓ=n−1\ell=n-1 だけで、

xHn(x)=Hn+1(x)+anHn−1(x),an=⟨xHn,Hn−1⟩hn−1xH_n(x)=H_{n+1}(x)+a_nH_{n-1}(x),\qquad a_n=\frac{\langle xH_n,H_{n-1}\rangle}{h_{n-1}}

が示されました。係数 ana_n は具体的に求まります。xe−x2=−12(e−x2)′xe^{-x^2}=-\tfrac{1}{2}\left(e^{-x^2}\right)' を使って部分積分すると

⟨xHn,Hn−1⟩=−12∫−∞∞(e−x2)′HnHn−1 dx=12∫−∞∞e−x2(Hn′Hn−1+HnHn−1′)dx\langle xH_n,H_{n-1}\rangle =-\frac{1}{2}\int_{-\infty}^{\infty}\left(e^{-x^2}\right)'H_nH_{n-1}\,dx =\frac{1}{2}\int_{-\infty}^{\infty}e^{-x^2}\left(H_n'H_{n-1}+H_nH_{n-1}'\right)dx

です。Hn−1′H_{n-1}' は次数 n−2n-2 なので第 2 項は 00、第 1 項は Hn′=nxn−1+(低次)H_n'=nx^{n-1}+(\text{低次}) より Hn′−nHn−1H_n'-nH_{n-1} が次数 n−2n-2 以下で Hn−1H_{n-1} と直交するため ⟨Hn′,Hn−1⟩=n hn−1\langle H_n',H_{n-1}\rangle=n\,h_{n-1} となります。よって ⟨xHn,Hn−1⟩=n2hn−1\langle xH_n,H_{n-1}\rangle=\tfrac{n}{2}h_{n-1}、すなわち

an=n2a_n=\frac{n}{2}

です。設問(1)で確かめると、xH1=x2=H2+12H0xH_1=x^2=H_2+\tfrac{1}{2}H_0 で a1=12a_1=\tfrac{1}{2}、xH2=x3−x2=H3+H1xH_2=x^3-\tfrac{x}{2}=H_3+H_1 で a2=1a_2=1 となり、an=n/2a_n=n/2 と一致します。なお途中の式は hn=⟨Hn,xHn−1⟩=n2hn−1h_n=\langle H_n,xH_{n-1}\rangle=\tfrac{n}{2}h_{n-1}(第 1 式の形から ⟨xHn−1,Hn⟩=⟨Hn,Hn⟩\langle xH_{n-1},H_n\rangle=\langle H_n,H_n\rangle)も与え、hn=n!2nπh_n=\tfrac{n!}{2^n}\sqrt{\pi} が従います。h1=μ2=12πh_1=\mu_2=\tfrac{1}{2}\sqrt{\pi}、h2=μ4−μ2+14μ0=12πh_2=\mu_4-\mu_2+\tfrac{1}{4}\mu_0=\tfrac{1}{2}\sqrt{\pi} と直接計算しても合います。

第 2 式。Hn′(x)H_n'(x) は次数 n−1n-1、最高次項 nxn−1nx^{n-1} の多項式なので、E(x)≡Hn′(x)−nHn−1(x)E(x)\equiv H_n'(x)-nH_{n-1}(x) は次数 n−2n-2 以下です。ℓ≤n−2\ell\le n-2 に対する展開係数は ⟨E,Hℓ⟩=⟨Hn′,Hℓ⟩−n⟨Hn−1,Hℓ⟩=⟨Hn′,Hℓ⟩\langle E,H_\ell\rangle=\langle H_n',H_\ell\rangle-n\langle H_{n-1},H_\ell\rangle=\langle H_n',H_\ell\rangle で決まります。部分積分すると

⟨Hn′,Hℓ⟩=[e−x2HnHℓ]−∞∞−∫−∞∞Hn(e−x2Hℓ)′dx=∫−∞∞e−x2Hn(2xHℓ−Hℓ′)dx\langle H_n',H_\ell\rangle =\left[e^{-x^2}H_nH_\ell\right]_{-\infty}^{\infty} -\int_{-\infty}^{\infty}H_n\left(e^{-x^2}H_\ell\right)'dx =\int_{-\infty}^{\infty}e^{-x^2}H_n\left(2xH_\ell-H_\ell'\right)dx

となり、2xHℓ−Hℓ′2xH_\ell-H_\ell' は次数 ℓ+1≤n−1\ell+1\le n-1 の多項式なので HnH_n と直交して 00 です。よって EE の展開係数はすべて 00、E≡0E\equiv 0 で

ddxHn(x)=nHn−1(x)\frac{d}{dx}H_n(x)=nH_{n-1}(x)

が示されました。設問(1)では H1′=1=1⋅H0H_1'=1=1\cdot H_0、H2′=2x=2H1H_2'=2x=2H_1、H3′=3x2−32=3(x2−12)=3H2H_3'=3x^2-\tfrac{3}{2}=3\left(x^2-\tfrac{1}{2}\right)=3H_2 となり、すべて一致します。

出典: 東京大学大学院理学系研究科 物理学専攻 平成11年度 修士課程 入学試験問題 数学(一般教育科目、数学・物理合冊のうち数学)。問題文は要約して引用しています。

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

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

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