Skip to content

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

Raw

This content is not available in your language yet.

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

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

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

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

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

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

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

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

f(x)=a02+n=1(ancosnx+bnsinnx)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 になることを示す問題です。

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

1Rsinαxxdx=[cosαxαx]1R1α1Rcosαxx2dx\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/x2x2|\cos\alpha x/x^2|\le x^{-2} より RR\to\infty で絶対収束します。よって I(α)I(\alpha) は広義リーマン積分として収束します(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=idθ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/\pisin\sin の凹性から従います)を用いて

ΓRg(z)dz=0πeiαReiθidθ0πeαRsinθdθ=20π/2eαRsinθdθ20π/2e2αRθ/πdθ=παR(1eα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 なので RR\to\infty でこれは 00 に収束します。ここで α\alpha の符号が本質的に効いていて、α<0\alpha\lt 0 なら上半平面では被積分関数が発散し、この評価は使えません。

小半円の寄与は、z=εeiθz=\varepsilon e^{i\theta}dz/z=idθ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 に対する eiw1wew|e^{iw}-1|\le|w|e^{|w|} から残差は παεeαε\pi\alpha\varepsilon e^{\alpha\varepsilon} 以下です。よって ε0\varepsilon\to 0γεgdziπ\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=εReiαuudu\int_{-R}^{-\varepsilon}\frac{e^{i\alpha x}}{x}dx=-\int_{\varepsilon}^{R}\frac{e^{-i\alpha u}}{u}du

なので、2 つの和は

εReiαxeiα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

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

2i0sinαxxdx=iπ2i\int_0^{\infty}\frac{\sin\alpha x}{x}dx=i\pi

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

0dxsinαxx=π2sgnα={  π/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 0y=αxy=\alpha x と置換すると α\alpha が完全に消えることから独立に確認できます。なお実部を取れば P ⁣cos(αx)/xdx=0\mathrm{P}\!\int_{-\infty}^{\infty}\cos(\alpha x)/x\,dx=0 で、これは被積分関数が奇関数であることと整合します。

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

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

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

an=1πππxcosnxdx=2π0πxcosnxdxa_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πxdx=2ππ22=πa_0=\frac{2}{\pi}\int_0^{\pi}x\,dx=\frac{2}{\pi}\cdot\frac{\pi^2}{2}=\pi

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

0πxcosnxdx=[xsinnxn]0π1n0πsinnxdx=0+[cosnxn2]0π=(1)n1n2\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)n1n2={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、偶数 n2n\ge 2an=0a_n=0、奇数 nnan=4/(πn2)a_n=-4/(\pi n^2)、そしてすべての nnbn=0b_n=0 です。フーリエ級数は

x=π24πk=0cos(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/2xπxx\to\pi-x の下で符号を変える(半区間ごとの反対称性をもつ)ことに対応しています。

まず級数が x=0x=0f(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) は連続関数です。他方 ff2π2\pi 周期に拡張したものは f(π)=f(π)=πf(-\pi)=f(\pi)=\pi より連続で、フェイェールの定理からそのフーリエ級数のチェザロ平均は各点で ff に収束します。収束する級数のチェザロ平均は同じ和に収束するので g=fg=f です。区分的に滑らかな連続関数に対するディリクレの定理を使って直接 g=fg=f と結論しても同じです(f(x)=sgnxf'(x)=\mathrm{sgn}\,xx=0,±πx=0,\pm\pi を除いて連続)。

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

0=π24πk=01(2k+1)2,すなわちk=01(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 を使うと、奇数項の和は n oddn2=ζ(2)m1(2m)2=ζ(2)(11/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)(11/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=TrAP=\operatorname{Tr}AQ=12(TrA2(TrA)2)Q=\frac{1}{2}\left(\operatorname{Tr}A^2-(\operatorname{Tr}A)^2\right)R=detAR=\det A を示す問題、設問(2)は P=0P=0 のときに AA が複素共役な(実でない)固有値の対をもつための条件を Q,RQ,R で表す問題です。

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

R=F(0)=det[A0I]=detAR=F(0)=\det[A-0\cdot I]=\det A

です。

PPQQ は行列式の置換展開

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=TrAP=a_{11}+a_{22}+a_{33}=\operatorname{Tr}A

です。

λ1\lambda^1 の係数には、恒等置換の項に加えて、固定点をちょうど 1 つもつ置換、すなわち 3 つの互換の項が寄与します(3-サイクルは固定点をもたず λ\lambda を含みません)。たとえば互換 σ=(23)\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} より

TrA2=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)

であり、また

(TrA)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(TrA2(TrA)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_iQ=i<jλiλjQ=-\sum_{i\lt j}\lambda_i\lambda_jR=λ1λ2λ3R=\lambda_1\lambda_2\lambda_3 です。TrA=iλi\operatorname{Tr}A=\sum_i\lambda_iTrA2=iλi2=(iλi)22i<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=e122e2p_2=e_1^2-2e_2 にほかなりません。

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

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

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

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

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

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

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

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

のとき、かつそのときに限り AA は複素共役の固有値をもちます。左辺 4Q327R24Q^3-27R^2λ3Qλ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 ivv0v\neq 0)なら (λ1λ2)(λ1λ3)=(au)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 行・列として零を足したもの)では TrA=0\operatorname{Tr}A=0TrA2=2\operatorname{Tr}A^2=-2 より Q=1Q=-1R=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=1R=0R=00>40\gt 4 は偽、実際すべて実固有値です。境界の例として固有値 1,1,21,1,-2 では Q=3Q=3R=2R=-227R2=108=4Q327R^2=108=4Q^3 と等号になり、重根(すべて実)の場合が等号に対応することが確認できます。

第3問 重み ex2e^{-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 より低次の多項式),dxex2Hn(x)Hm(x)=0(nm)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=ex2f(x)g(x)dx\langle f,g\rangle=\int_{-\infty}^{\infty}e^{-x^2}f(x)g(x)\,dx、ノルムの 2 乗を hn=Hn,Hnh_n=\langle H_n,H_n\rangle と書きます。重み ex2e^{-x^2} は正で Hn≢0H_n\not\equiv 0 なので hn>0h_n\gt 0 です。与えられた積分公式はモーメント μk=xkex2dx\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==0ncHP=\sum_{\ell=0}^{n}c_\ell H_\ell と展開できる」を、次の形の補題として全設問で使います。H0,,HnH_0,\ldots,H_n は次数が互いに異なるので一次独立で、nn 次以下の多項式の空間(次元 n+1n+1)の基底をなし、展開係数は直交性から c=P,H/hc_\ell=\langle P,H_\ell\rangle/h_\ell と一意に定まります。ここから 2 つの系が出ます。第 1 に、HnH_n は次数 n1n-1 以下の任意の多項式と直交します(相手を H0,,Hn1H_0,\ldots,H_{n-1} で展開すればよい)。第 2 に、最高次係数 1 の nn 次多項式で H0,,Hn1H_0,\ldots,H_{n-1} のすべてと直交するものは HnH_n に限ります。実際そのような KK に対し KHnK-H_n は次数 n1n-1 以下で全 HH_\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,,xn11,x,\ldots,x^{n-1} との直交(上の系より H0,,Hn1H_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)=x212H_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=π43π4=π20\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)=x332xH_3(x)=x^3-\frac{3}{2}x

です。検算として、物理で標準のエルミート多項式 H2phys=4x22H_2^{\mathrm{phys}}=4x^2-2H3phys=8x312xH_3^{\mathrm{phys}}=8x^3-12x を最高次係数で割ると x212x^2-\tfrac{1}{2}x332xx^3-\tfrac{3}{2}x となり一致します。また H1,H3=μ432μ2=(3434)π=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) と置き、すべての nnKn=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 なるすべての mmKm=HmK_m=H_m、すなわち Hm(x)=(1)mHm(x)H_m(-x)=(-1)^mH_m(x) が成り立つと仮定します。m<nm\lt n に対し、置換 x=tx=-tdx=dtdx=-dt、積分範囲は \infty\to-\infty が反転して元に戻り、ex2e^{-x^2} は偶関数)を使うと

Kn,Hm=(1)nex2Hn(x)Hm(x)dx=(1)net2Hn(t)Hm(t)dt=(1)n+met2Hn(t)Hm(t)dt=(1)n+mHn,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,,Hn1H_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=x212H_2=x^2-\tfrac{1}{2}(偶)、H3=x332xH_3=x^3-\tfrac{3}{2}x(奇)と整合します。

以下 n1n\ge 1 とします(Hn1H_{n-1} が現れるため)。また境界項について、多項式 pp に対し ex2p(x)0e^{-x^2}p(x)\to 0x±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==0ncHD=\sum_{\ell=0}^{n}c_\ell H_\ellch=D,Hc_\ell h_\ell=\langle D,H_\ell\rangle と展開できます。n\ell\le n では Hn+1,H=0\langle H_{n+1},H_\ell\rangle=0 なので ch=xHn,H=Hn,xHc_\ell h_\ell=\langle xH_n,H_\ell\rangle=\langle H_n,xH_\ell\rangle です(重み付き内積の中で xx を移しただけ)。n2\ell\le n-2 なら xHxH_\ell は次数 +1n1\ell+1\le n-1 の多項式なので HnH_n と直交し、c=0c_\ell=0 です。=n\ell=n では

xHn,Hn=ex2xHn(x)2dx=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 倍した被積分関数が奇関数になるからです。残るのは =n1\ell=n-1 だけで、

xHn(x)=Hn+1(x)+anHn1(x),an=xHn,Hn1hn1xH_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 は具体的に求まります。xex2=12(ex2)xe^{-x^2}=-\tfrac{1}{2}\left(e^{-x^2}\right)' を使って部分積分すると

xHn,Hn1=12(ex2)HnHn1dx=12ex2(HnHn1+HnHn1)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

です。Hn1H_{n-1}' は次数 n2n-2 なので第 2 項は 00、第 1 項は Hn=nxn1+(低次)H_n'=nx^{n-1}+(\text{低次}) より HnnHn1H_n'-nH_{n-1} が次数 n2n-2 以下で Hn1H_{n-1} と直交するため Hn,Hn1=nhn1\langle H_n',H_{n-1}\rangle=n\,h_{n-1} となります。よって xHn,Hn1=n2hn1\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_0a1=12a_1=\tfrac{1}{2}xH2=x3x2=H3+H1xH_2=x^3-\tfrac{x}{2}=H_3+H_1a2=1a_2=1 となり、an=n/2a_n=n/2 と一致します。なお途中の式は hn=Hn,xHn1=n2hn1h_n=\langle H_n,xH_{n-1}\rangle=\tfrac{n}{2}h_{n-1}(第 1 式の形から xHn1,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) は次数 n1n-1、最高次項 nxn1nx^{n-1} の多項式なので、E(x)Hn(x)nHn1(x)E(x)\equiv H_n'(x)-nH_{n-1}(x) は次数 n2n-2 以下です。n2\ell\le n-2 に対する展開係数は E,H=Hn,HnHn1,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=[ex2HnH]Hn(ex2H)dx=ex2Hn(2xHH)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

となり、2xHH2xH_\ell-H_\ell' は次数 +1n1\ell+1\le n-1 の多項式なので HnH_n と直交して 00 です。よって EE の展開係数はすべて 00E0E\equiv 0

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

が示されました。設問(1)では H1=1=1H0H_1'=1=1\cdot H_0H2=2x=2H1H_2'=2x=2H_1H3=3x232=3(x212)=3H2H_3'=3x^2-\tfrac{3}{2}=3\left(x^2-\tfrac{1}{2}\right)=3H_2 となり、すべて一致します。

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

Report an error in this article ・Operated by: Mugen Giken LLCPricingTermsLegal notice

© 2026 夢現技研合同会社 ・Feeding the text to an LLM is welcome. Code samples are MIT licensed.