跳到主要内容

東京大学 総合文化研究科 広域科学専攻 相関基礎科学系 2024年度 総合科目 第1問(数学)

Author​

祭音Myyura (Based on Miyake's answer refined with GPT 6 Astra)

Description​

以下の I〜V に答えよ。

I​

実変数 xx の実関数 y=y(x)y=y(x) に対する次の常微分方程式から一つを選び、選んだ番号を明記して一般解を初等関数で書け。ただし、初等関数の範囲で解けないものも含まれる。

(a)y′′−xy=0,(b)2xyy′+y2+x2=0,(c)y′+x2y+x5=0,(d)y′+y2+2x=0.\begin{array}{ll} \mathrm{(a)}&y''-xy=0,\\ \mathrm{(b)}&2xyy'+y^2+x^2=0,\\ \mathrm{(c)}&y'+x^2y+x^5=0,\\ \mathrm{(d)}&y'+y^2+2x=0. \end{array}

II​

次の定積分を n=0,1,2n=0,1,2 について求めよ。

In=∫0∞(log⁡x)nx2+1 dx.I_n=\int_0^\infty\frac{(\log x)^n}{x^2+1}\,dx.

III​

次の命題が真なら証明し、偽なら反例を挙げてそれが反例であることを示せ。

「実数から実数への微分可能な関数 f(x)f(x) が x→∞x\to\infty で収束するならば、lim⁡x→∞f′(x)=0\lim_{x\to\infty}f'(x)=0 である。」

IV​

nn を正の整数とする。nn 次エルミート行列 AA の固有値がすべて正の実数であるとき、AA は正定値であるという。正定値エルミート行列には A=C2A=C^2 を満たす可逆なエルミート行列 CC が存在することを用いてよい。以下を示せ。

  1. AA が正定値であるための必要十分条件は、任意の非零ベクトル v∈Cnv\in\mathbb C^n に対して (v,Av)>0(v,Av)>0 である。ただし (v,w)=∑i=1nv‾iwi(v,w)=\sum_{i=1}^n\overline v_iw_i。
  2. AA が正定値エルミート行列ならば、任意の可逆行列 XX に対して X†AXX^\dagger AX も正定値エルミート行列である。
  3. 一般の nn 次正方行列 X,YX,Y について、XX が可逆ならば XYXY と YXYX の固有値は一致する。
  4. A,BA,B が正定値エルミート行列ならば、ABAB の固有値はすべて正の実数である。

V​

t≥0t\geq0 における x∈Rx\in\mathbb R の確率分布関数 P(x,t)P(x,t) が

∂P∂t=∂2P∂x2+2∂P∂x\frac{\partial P}{\partial t}=\frac{\partial^2P}{\partial x^2}+2\frac{\partial P}{\partial x}

を満たすとする。

  1. 初期分布 P(x,0)=Ce−3x2P(x,0)=Ce^{-3x^2} が確率分布となるように CC を求めよ。
  2. P(x,0)P(x,0) の xx に関するフーリエ変換 P^(k,0)\widehat P(k,0) を求めよ。
  3. t>0t>0 における P(x,t)P(x,t) を求めよ。

题目描述​

I. 从上述四个常微分方程中任选一个,用初等函数写出通解,并注明选择的编号;其中有的方程不能用初等函数求解。

II. 求 In=∫0∞(log⁡x)n/(1+x2) dxI_n=\int_0^\infty (\log x)^n/(1+x^2)\,dx,其中 n=0,1,2n=0,1,2。

III. 判断“定义在全体实数上的可微实函数若在正无穷处收敛,则其导数趋于零”是否成立,并给出证明或反例。

IV. 证明正定厄米矩阵的二次型判据、可逆合同变换保持正定、可逆 XX 时 XYXY 与 YXYX 同谱,以及两个正定厄米矩阵之积的全部特征值均为正实数。可使用正定厄米矩阵存在可逆厄米平方根这一事实。

V. 对扩散方程 Pt=Pxx+2PxP_t=P_{xx}+2P_x,将初值 Ce−3x2Ce^{-3x^2} 归一化,求其傅里叶变换以及任意 t>0t>0 的概率密度。

Kai​

I​

(c) を選ぶ。積分因子 ex3/3e^{x^3/3} を掛けると

ddx(yex3/3)=−x5ex3/3.\frac{d}{dx}\left(ye^{x^3/3}\right)=-x^5e^{x^3/3}.

u=x3/3u=x^3/3 と置いて積分すれば ∫x5ex3/3dx=(x3−3)ex3/3\int x^5e^{x^3/3}dx=(x^3-3)e^{x^3/3} なので、任意定数 C1∈RC_1\in\mathbb R を用いて

y=C1e−x3/3−x3+3.\boxed{y=C_1e^{-x^3/3}-x^3+3}.

II​

−1<a<1-1<a<1 に対して、t=x2t=x^2 と置くと

J(a)=∫0∞xa1+x2 dx=12∫0∞t(a−1)/21+t dt=π2cos⁡(πa/2).J(a)=\int_0^\infty\frac{x^a}{1+x^2}\,dx =\frac12\int_0^\infty\frac{t^{(a-1)/2}}{1+t}\,dt =\frac{\pi}{2\cos(\pi a/2)}.

最後の等式はベータ積分と Γ(z)Γ(1−z)=π/sin⁡(πz)\Gamma(z)\Gamma(1-z)=\pi/\sin(\pi z) による。a=0a=0 の近傍で積分内微分できるため、

I0=J(0)=π2,I1=J′(0)=0,I2=J′′(0)=π38.\boxed{I_0=J(0)=\frac\pi2,\qquad I_1=J'(0)=0,\qquad I_2=J''(0)=\frac{\pi^3}{8}}.

III​

命題は偽である。反例として

f(x)={sin⁡(x2)/x,x≠0,0,x=0f(x)=\begin{cases}\sin(x^2)/x,&x\ne0,\\0,&x=0\end{cases}

をとる。f′(0)=lim⁡h→0sin⁡(h2)/h2=1f'(0)=\lim_{h\to0}\sin(h^2)/h^2=1 であり、全実数上で微分可能である。また ∣f(x)∣≤1/∣x∣|f(x)|\leq1/|x| より f(x)→0f(x)\to0。一方、x≠0x\ne0 では

f′(x)=2cos⁡(x2)−sin⁡(x2)x2.f'(x)=2\cos(x^2)-\frac{\sin(x^2)}{x^2}.

xn=2πnx_n=\sqrt{2\pi n} では f′(xn)=2f'(x_n)=2、yn=(2n+1)πy_n=\sqrt{(2n+1)\pi} では f′(yn)=−2f'(y_n)=-2 なので、導関数は 00 に収束しない。

IV​

(1)​

スペクトル定理より A=Udiag⁡(λ1,…,λn)U†A=U\operatorname{diag}(\lambda_1,\ldots,\lambda_n)U^\dagger と書ける。w=U†vw=U^\dagger v とすれば

v†Av=∑i=1nλi∣wi∣2.v^\dagger Av=\sum_{i=1}^n\lambda_i|w_i|^2.

全固有値が正なら、v≠0v\ne0 に対して右辺は正である。逆に固有ベクトル vi≠0v_i\ne0 を代入すると λi∥vi∥2>0\lambda_i\|v_i\|^2>0 なので全固有値が正になる。

(2)​

(X†AX)†=X†AX(X^\dagger AX)^\dagger=X^\dagger AX。XX の可逆性により v≠0v\ne0 なら Xv≠0Xv\ne0 であり、

v†X†AXv=(Xv)†A(Xv)>0.v^\dagger X^\dagger AXv=(Xv)^\dagger A(Xv)>0.

(1) より正定値である。

(3)​

X−1(XY)X=YXX^{-1}(XY)X=YX より両者は相似であり、重複度を含めて固有値が一致する。

(4)​

A=C2A=C^2 を満たす可逆エルミート行列 CC をとると

C−1(AB)C=CBC=C†BC.C^{-1}(AB)C=CBC=C^\dagger BC.

右辺は (2) より正定値エルミート行列である。相似変換は固有値を変えないので ABAB の全固有値は正の実数である。

V​

(1)​

∫−∞∞e−3x2dx=π/3\int_{-\infty}^{\infty}e^{-3x^2}dx=\sqrt{\pi/3} より

C=3/π.\boxed{C=\sqrt{3/\pi}}.

(2)​

変換を P^(k,t)=∫RP(x,t)e−ikxdx\widehat P(k,t)=\int_{\mathbb R}P(x,t)e^{-ikx}dx と定義する。ガウス積分より

P^(k,0)=e−k2/12.\boxed{\widehat P(k,0)=e^{-k^2/12}}.

(3)​

方程式をフーリエ変換すると ∂tP^=(−k2+2ik)P^\partial_t\widehat P=(-k^2+2ik)\widehat P。したがって

P^(k,t)=exp⁡[−(t+112)k2+2ikt].\widehat P(k,t)=\exp\left[-\left(t+\frac1{12}\right)k^2+2ikt\right].

逆変換により

P(x,t)=3π(1+12t)exp⁡[−3(x+2t)21+12t].\boxed{P(x,t)=\sqrt{\frac{3}{\pi(1+12t)}}\exp\left[-\frac{3(x+2t)^2}{1+12t}\right]}.

これは平均 −2t-2t、分散 1/6+2t1/6+2t の正規分布であり、初期条件と規格化を満たす。

Reference​