院試hub

東京大学 院試 過去問 解答例

東大 数理科学研究科 数理科学専攻 専門科目B 2013年度 院試 解答例・解説

東京大学 数理科学研究科 数理科学専攻 専門科目B 2013年度の院試 過去問について、設問ごとの解法方針と確認点を解説。全18問収録の解答・解説PDFと併用できます。問題本文は含みません。

最終更新:

1 — 有限体上の2次一般線形群

標数 pp の落とし穴

複素数上なら pp 乗して1になる固有値は pp 個あるが,Fp\mathbb F_p 上では Xp1X^p-1 が重根 (X1)p(X-1)^p になる。この一点で,対象は 「半単純な巡回行列」ではなく「単位行列に冪零行列を足したもの」に変わる。

数え上げ

2次では非零冪零行列のJordan型は一種類だけである。したがって,共役類の大きさを 中心化群で割って数えるのが最短である。最後に N=0N=0 を忘れると p21p^2-1 で止まってしまう。

2013年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

(1)

gp=eg^p=e とする。代数閉包上で gg の固有値を α\alpha とすれば αp=1\alpha^p=1 である。標数 pp では Xp1=(X1)p X^p-1=(X-1)^p なので,固有値は重複度も込めてすべて 11 である。

(2)

従って g=I+Ng=I+N と書け,NN は2次の冪零行列である。逆に,任意の 冪零行列 NN について N2=0N^2=0 であるから (I+N)p=I+Np=I (I+N)^p=I+N^p=I が成り立つ。よって個数は,2次正方行列の冪零行列の個数に等しい。

非零冪零行列はすべて J=(0100) J=\begin{pmatrix}0&1\\0&0\end{pmatrix} に共役である。GL2(Fp)GL_2(\mathbb F_p) における JJ の中心化群は {aI+bJaFp×, bFp} \{aI+bJ\mid a\in \mathbb F_p^\times,\ b\in \mathbb F_p\} で,その位数は p(p1)p(p-1) である。一方 GL2(Fp)=(p21)(p2p)=p(p1)2(p+1) |GL_2(\mathbb F_p)|=(p^2-1)(p^2-p)=p(p-1)^2(p+1) だから,非零冪零行列は GL2(Fp)p(p1)=p21 \frac{|GL_2(\mathbb F_p)|}{p(p-1)}=p^2-1 個である。零行列も加えると,条件を満たす ggp2 p^2 個である。

最終答

固有値はすべて 11。また gp=eg^p=e となる元の個数は p2 \boxed{p^2} である。

2 — 3本の直線の有限被覆

幾何で見る

BB は3本の座標軸を原点で貼り合わせたものの座標環である。 T=x+y+zT=x+y+z は各枝でその枝の座標になる。一般点 T0T\ne0 では3枝が分離し, それが BAKK3B\otimes_AK\simeq K^3 という分解である。

商の次元

BA3B\subset A^3 は「原点で3つの値が等しい」という1条件ではなく, 3つの値を1つにそろえる条件である。したがって自由な3次元の値から対角1次元を 残すので,ずれの自由度は2次元である。

2013年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

(1)

BB の元を,3本の座標軸上の関数の組として見る。すなわち B{(u,v,w)C[X]C[Y]C[Z]u(0)=v(0)=w(0)} B\simeq \{(u,v,w)\in \mathbb C[X]\oplus\mathbb C[Y]\oplus\mathbb C[Z] \mid u(0)=v(0)=w(0)\} である。この表示では TT は各成分で座標関数として作用する。 また x=(T,0,0),y=(0,T,0),z=(0,0,T) x=(T,0,0),\qquad y=(0,T,0),\qquad z=(0,0,T) と書ける。

任意の元は a(T)+b(T)x+c(T)y a(T)+b(T)x+c(T)y と一意に表せる。実際,第3成分から a(T)a(T) が決まり,第1,第2成分から b(T),c(T)b(T),c(T) が決まる。従って BB は階数3の自由 AA-加群である。

(2)

K=C(T)K=\mathbb C(T) 上に移すと BAKKexKeyKez,ex=xT,ey=yT,ez=zT B\otimes_A K\simeq K e_x\oplus K e_y\oplus K e_z,\qquad e_x=\frac{x}{T},\quad e_y=\frac{y}{T},\quad e_z=\frac{z}{T} となる。ここで ex,ey,eze_x,e_y,e_z は互いに直交する原始冪等元で, ex+ey+ez=1e_x+e_y+e_z=1 である。よって冪等元は iIei(I{x,y,z}) \sum_{i\in I} e_i\qquad (I\subset \{x,y,z\}) の8個である。

(3)

これらの冪等元を BB に付け加えて生成される部分環は C=AexAeyAezA3 C=Ae_x\oplus Ae_y\oplus Ae_z\simeq A^3 である。上の組表示では B={(u,v,w)A3uvw(modT)},C=A3. B=\{(u,v,w)\in A^3\mid u\equiv v\equiv w \pmod T\},\qquad C=A^3. したがって C/BC/BT=0T=0 での3成分から対角成分を割ったものだから dimCC/B=31=2. \dim_{\mathbb C} C/B=3-1=2.

(4)

最後に単元群を比べる。BB は連結な3本のアフィン直線の座標環なので B×=C× B^\times=\mathbb C^\times である。一方 C×=(A×)3=(C×)3. C^\times=(A^\times)^3=(\mathbb C^\times)^3. 包含写像は対角埋め込みに対応するため,同型ではない。

最終答

BAAxAy,rankAB=3 \boxed{B\simeq A\oplus Ax\oplus Ay,\quad \operatorname{rank}_A B=3} Idem(BAK)={iIeiI{x,y,z}} \boxed{\operatorname{Idem}(B\otimes_AK)=\left\{\sum_{i\in I}e_i\mid I\subset\{x,y,z\}\right\}} dimCC/B=2,B×C× は同型でない。 \boxed{\dim_{\mathbb C} C/B=2,\qquad B^\times\to C^\times\text{ は同型でない。}}

3 — Nakayamaの補題と単項曲線

前半の要点

rrNN で非零因子であることを使う場所は,rϕ(m)=0r\phi(m)=0 から ϕ(m)=0\phi(m)=0 を出す一箇所である。その後は有限生成性を使ってNakayamaの補題を 適用する。

二項式の見つけ方

単項曲線の核では,同じ tt-次数をもつ単項式の差が関係式になる。 ここでは 26=5+7,45=6+27,37=35+6 2\cdot6=5+7,\qquad 4\cdot5=6+2\cdot7,\qquad 3\cdot7=3\cdot5+6 から上の3本が自然に出る。生成性は「どの単項式も有限個の代表に落とせる」ことを 確認するのが実戦的である。

2013年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

(1)

K=kerϕK=\ker\phi とおく。自然写像 K/rKM/rMK/rK\to M/rM の核を調べる。 kKk\in KM/rMM/rM で0に写るとすると,ある mMm\in M により k=rmk=rm と書ける。このとき 0=ϕ(k)=rϕ(m) 0=\phi(k)=r\phi(m) であり,rrNN の非零元を殺さないから ϕ(m)=0\phi(m)=0 である。 従って mKm\in K,すなわち krKk\in rK である。よって単射である。

(2)

次に ϕˉ:M/rMN/rN\bar\phi:M/rM\to N/rN が同型とする。単射性から K/rK=0K/rK=0 なので,Nakayamaの補題により K=0K=0 である。さらに余核 C=cokerϕC=\operatorname{coker}\phi について C/rC=0C/rC=0 となるので,再び Nakayamaの補題から C=0C=0 である。従って ϕ\phi は同型である。

(3)

最後に f:C[[x,y,z]]C[[t]],xt5,yt6,zt7 f:\mathbb C[[x,y,z]]\longrightarrow \mathbb C[[t]],\qquad x\mapsto t^5,\quad y\mapsto t^6,\quad z\mapsto t^7 を考える。次の3つの二項式は明らかに核に入る。 y2xz,x4yz2,z3x3y. y^2-xz,\qquad x^4-yz^2,\qquad z^3-x^3y. これらが生成するイデアルを II とおく。

商環では y2=xzy^2=xz により yy の次数を高々1にでき,さらに x4=yz2, z3=x3yx^4=yz^2,\ z^3=x^3y により高い冪を落とせる。したがって任意の単項式は, C[[x]]\mathbb C[[x]]-加群として 1, y, z, yz, z2, yz2 1,\ y,\ z,\ yz,\ z^2,\ yz^2 の線形結合に整理できる。tt-進値を見ると,これらの値はそれぞれ 0, 6, 7, 13, 14, 20 0,\ 6,\ 7,\ 13,\ 14,\ 20 を代表し,半群 5,6,7\langle5,6,7\rangle の剰余類を尽くす。この正規形の像が0なら 各係数が0になるので,核は II に一致する。

最終答

kerf=(y2xz, x4yz2, z3x3y) \boxed{\ker f=(\,y^2-xz,\ x^4-yz^2,\ z^3-x^3y\,)} である。

4 — 3次Heisenberg型の不変体

軌道と最小多項式

固定体の問題では,まず作用群を明確にする。XX は単に X,Y,ZX,Y,Z に動くだけでなく,同時スカラー倍も受けるため,軌道の大きさは3ではなく 9になる。

中間体の数え方

中間体の共役類は部分群の共役類で数える。ここで現れる群は3次のHeisenberg群で, 中心を除いた位数3部分群は,商 F32\mathbb F_3^2 の直線に従って4種類に分かれる。

2013年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

(1)

ρ\rhoρ(X)=ωX,ρ(Y)=ωY,ρ(Z)=ωZ \rho(X)=\omega X,\qquad \rho(Y)=\omega Y,\qquad \rho(Z)=\omega Z で定まる自己同型とする。計算すると στ=ρ2τσ \sigma\tau=\rho^2\tau\sigma であり,ρ\rho は中心的である。従って σ,τ\sigma,\tau が生成する群 GG は 位数 2727 の非可換な3群であり,K=LGK=L^G である。

XXGG-軌道は {ωiX, ωiY, ωiZi=0,1,2} \{\omega^iX,\ \omega^iY,\ \omega^iZ\mid i=0,1,2\} である。よって XXKK 上の最小多項式は i=02(UωiX)(UωiY)(UωiZ)=(U3X3)(U3Y3)(U3Z3) \prod_{i=0}^{2}(U-\omega^iX)(U-\omega^iY)(U-\omega^iZ) =(U^3-X^3)(U^3-Y^3)(U^3-Z^3) であり,次数は9である。したがって [K(X):K]=9. [K(X):K]=9.

(2)

GGLL に忠実に作用する有限群なので [L:K]=G=27. [L:K]=|G|=27. (3)

可換な最大商は G/[G,G]G/[G,G] である。ここで交換子群は中心 ρ\langle\rho\rangle で,位数3である。従って,KK 上アーベルとなる最大中間体 の次数は G/[G,G]=9. |G/[G,G]|=9.

(4)

最後に [L:M]=3[L:M]=3 となる中間体は,GG の位数3の部分群に対応する。中心 ρ\langle\rho\rangle は1つの共役類を与える。非中心的な位数3部分群は G/ρF32G/\langle\rho\rangle\simeq \mathbb F_3^2 の1次元部分空間ごとに1つの共役類に まとまる。1次元部分空間は 32131=4 \frac{3^2-1}{3-1}=4 個であるから,共役類の総数は 1+4=5 1+4=5 である。

最終答

[K(X):K]=9,[L:K]=27 \boxed{[K(X):K]=9,\qquad [L:K]=27} 最大アーベル中間体の次数は 9,[L:M]=3 の共役類は 5 個 \boxed{\text{最大アーベル中間体の次数は }9,\qquad [L:M]=3\text{ の共役類は }5\text{ 個}} である。

5 — 外微分と楕円体上の積分

向きの決め方

閉曲面の積分では,まず内部領域を決め,その境界として外向きに向きを入れるのが 最も安全である。この向きなら Vω=Vdω\int_{\partial V}\omega=\int_V d\omega と符号を 迷わず使える。

体積に落とす

この設問の核心は,微分形式そのものを曲面上で直接積分しないことにある。 dωd\omega が標準体積形式になるよう ff を選ぶため,曲面積分は楕円体の体積に 変換される。

2013年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

(1)

ω=xdydz2zf(y)dxdy+yf(y)dzdx \omega=x\,dy\wedge dz-2zf(y)\,dx\wedge dy+yf(y)\,dz\wedge dx とおく。外微分を計算すると dω={1f(y)+yf(y)}dxdydz. d\omega=\{1-f(y)+yf'(y)\}\,dx\wedge dy\wedge dz. 従って dω=dxdydzd\omega=dx\wedge dy\wedge dz となる条件は yf(y)f(y)=0 yf'(y)-f(y)=0 である。C1C^1 性と f(1)=1f(1)=1 から f(y)=y f(y)=y が得られる。

(2)

曲面 SSx2+y2+z2+xy+yz+zx=12 x^2+y^2+z^2+xy+yz+zx=\frac12 で与えられる楕円体である。内部を V={(x,y,z)x2+y2+z2+xy+yz+zx12} V=\left\{(x,y,z)\mid x^2+y^2+z^2+xy+yz+zx\le \frac12\right\} とし,S=VS=\partial V には外向き法線で向きを入れる。この向きでStokesの定理を 用いると Sω=Vdω=vol(V). \int_S\omega=\int_V d\omega=\operatorname{vol}(V).

二次形式の行列は A=(11/21/21/211/21/21/21) A=\begin{pmatrix} 1&1/2&1/2\\ 1/2&1&1/2\\ 1/2&1/2&1 \end{pmatrix} で,detA=1/2\det A=1/2 である。したがって vol(V)=4π3(1/2)3/2detA=2π3. \operatorname{vol}(V) =\frac{4\pi}{3}\frac{(1/2)^{3/2}}{\sqrt{\det A}} =\frac{2\pi}{3}.

最終答

f(y)=y,Sω=2π3 \boxed{f(y)=y},\qquad \boxed{\displaystyle \int_S\omega=\frac{2\pi}{3}} である。ただし SS は楕円体内部の境界として外向きに向きを入れる。

6 — 回転面の測地線とClairaut関係

保存量の由来

回転面では角度 θ\theta が循環座標である。従って f(x)2θ˙f(x)^2\dot\theta が保存される。これを幾何的な量に直すと rsinwr\sin w になり,これはClairautの関係と呼ばれる。

指数型回転面での判定

f(x)=exf(x)=e^x では回転半径 rrxx とともに単調に小さくなる。 Clairaut定数が正なら rr はその定数より小さくなれないため,下方へ無限には 進めない。下に抜けられるのは定数が0,すなわち純粋な子午線方向だけである。

2013年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

(1)

回転面を X(x,θ)=(x,f(x)cosθ,f(x)sinθ) X(x,\theta)=(x,f(x)\cos\theta,f(x)\sin\theta) でパラメータ表示する。このとき第一基本形式は ds2=(1+f(x)2)dx2+f(x)2dθ2 ds^2=(1+f'(x)^2)\,dx^2+f(x)^2\,d\theta^2 である。

まず θ=α\theta=\alpha が定数の曲線を考える。この曲線の速度は XxX_x 方向であり, 加速度の接成分は,曲線を弧長で取り直したとき消える。これは子午線が回転面の 反射対称面との交わりであり,その対称性から測地線になる,と言ってもよい。

(2)

次に一般の測地線を弧長 ss で表す。Lagrangian L=12{(1+f(x)2)x˙2+f(x)2θ˙2} L=\frac12\{(1+f'(x)^2)\dot x^2+f(x)^2\dot\theta^2\} θ\theta を含まないので,共役運動量 f(x)2θ˙ f(x)^2\dot\theta は一定である。速度ベクトルと子午線方向のなす角を ww とすると,回転方向成分の 大きさは f(x)θ˙f(x)|\dot\theta| である。弧長パラメータでは sinw=f(x)θ˙ \sin w=f(x)|\dot\theta| だから rsinw=f(x)sinw=f(x)2θ˙ r\sin w=f(x)\sin w=f(x)^2|\dot\theta| は一定である。

(3)

特に f(x)=exf(x)=e^x とする。Clairaut定数を c=rsinwc=r\sin w と書く。もし c>0c>0 なら r=excr=e^x\ge c でなければならないから xlogc x\ge \log c であり,xx 座標は下に有界である。従って xx 座標が下に有界でない測地線では c=0c=0 である。このとき sinw=0\sin w=0 なので,測地線は子午線方向に進み, パラメータの取り替えを除いて θ=α\theta=\alpha 一定の形である。

最終答

子午線は測地線であり,任意の測地線について r(t)sinw(t)=一定 \boxed{r(t)\sin w(t)=\text{一定}} が成り立つ。f(x)=exf(x)=e^x の場合,xx 座標が下に有界でない測地線は, パラメータの取り替えを除き子午線またはその一部に限られる。

7 — 5平面配置のホモロジー

包含排除がそのまま効く

平面,交線,三重点はいずれも可縮で,4重交点がないため,オイラー標数は 交わりの個数だけで計算できる。

神経の見方

この配置は「5頂点のうち3頂点までの交わりはあるが,4頂点の交わりはない」 という構造を持つ。したがって神経は4単体の2骨格であり,2次元の穴が4個残る。

2013年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

(1)

5個の平面を H1,,H5H_1,\dots,H_5 とする。各 HiH_i は可縮で,2平面の交わりは 直線,3平面の交わりは1点,4平面の交わりは空である。包含排除で χ(X)=(51)(52)+(53)=510+10=5. \chi(X)=\binom51-\binom52+\binom53=5-10+10=5.

(2)

単連結性は,平面を順に貼り合わせて示す。既に貼った和集合に新しい平面を加える とき,共通部分はその平面内の有限本の直線の和であり,連結である。平面自身は単連結 であり,共通部分も連結なので,van Kampenの定理により基本群は増えない。最初の 平面から始めて帰納的に π1(X)=0 \pi_1(X)=0 が従う。

(3)

またこの配置の良い近傍被覆を取ると,神経は5頂点の単体の2骨格になる。したがって XX は2次元CW複体と同じホモトピー型を持つ。よって Hi(X)=0H_i(X)=0i3i\ge3 で成り立ち,単連結性から H1(X;Z)=0H_1(X;\mathbb Z)=0 である。 従って χ(X)=rankH0rankH1+rankH2=1+rankH2. \chi(X)=\operatorname{rank}H_0-\operatorname{rank}H_1+\operatorname{rank}H_2 =1+\operatorname{rank}H_2. ゆえに H2(X;Z)Z4. H_2(X;\mathbb Z)\simeq \mathbb Z^4.

最終答

χ(X)=5,π1(X)=0 \boxed{\chi(X)=5},\qquad \boxed{\pi_1(X)=0} また整数係数ホモロジー群は H0(X;Z)Z,H1(X;Z)=0,H2(X;Z)Z4,Hi(X;Z)=0 (i3). \boxed{ H_0(X;\mathbb Z)\simeq \mathbb Z,\quad H_1(X;\mathbb Z)=0,\quad H_2(X;\mathbb Z)\simeq \mathbb Z^4,\quad H_i(X;\mathbb Z)=0\ (i\ge3). }

8 — 随伴軌道と不変な複素構造

不変性の使いどころ

存在しないことを示すには,全点を調べる必要はない。1点 A0A_0 の安定化群に 注目し,その線形作用と可換な JA0J_{A_0} があり得るかだけを見れば十分である。

対角作用との可換性

安定化群は接空間の2方向を異なる重みで伸縮する。このような作用と可換な写像は 2方向を混ぜられないため,複素構造に必要な90度回転型の写像を作れない。

2013年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

(1)

跡0の行列を A=(xyzx) A=\begin{pmatrix}x&y\\ z&-x\end{pmatrix} と書くと detA=x2yz. \det A=-x^2-yz. 従って X={(x,y,z)R3x2+yz=1} X=\{(x,y,z)\in\mathbb R^3\mid x^2+yz=1\} である。関数 F(x,y,z)=x2+yzF(x,y,z)=x^2+yz の勾配はこの集合上で消えないので, XXM2(R)M_2(\mathbb R) の部分多様体である。

(2)

gSL2(R)g\in SL_2(\mathbb R) に対して AgAg1A\mapsto gAg^{-1} は跡と行列式を保つ。 逆写像は Ag1AgA\mapsto g^{-1}Ag なので,これは XX から XX への微分同相写像で ある。

(3)

最後に,条件を満たす JJ が存在すると仮定して矛盾を導く。 A0=(1001) A_0=\begin{pmatrix}1&0\\0&-1\end{pmatrix} とする。接空間は TA0X={(0uv0)|u,vR}=RERF T_{A_0}X=\left\{ \begin{pmatrix}0&u\\v&0\end{pmatrix}\middle| u,v\in\mathbb R \right\} =\mathbb R E\oplus\mathbb R F である。ただし E=(0100),F=(0010). E=\begin{pmatrix}0&1\\0&0\end{pmatrix},\qquad F=\begin{pmatrix}0&0\\1&0\end{pmatrix}. A0A_0 を固定する対角行列 gs=(s00s1)(s>0) g_s=\begin{pmatrix}s&0\\0&s^{-1}\end{pmatrix}\qquad (s>0) の作用は Es2E,Fs2F E\mapsto s^2E,\qquad F\mapsto s^{-2}F である。条件(III)から JA0J_{A_0} はこの対角作用と可換でなければならない。 s41s^4\ne1 となる ss を取ると,可換な線形写像は E,FE,F 方向を別々に保つ 対角型に限られる。しかし実2次元空間上の対角型線形写像の二乗が I-I になることは ない。これは条件(II)に反する。

最終答

XX は2次元部分多様体であり,SL2(R)SL_2(\mathbb R) の共役作用は XX の微分同相を定める。ただし条件(I)--(III)をすべて満たす J:TXTX \boxed{J:TX\to TX} は存在しない。

9 — 単射正則関数と長さの単調性

平方根を作る理由

長さには f|f'| が現れる。f=g2f'=g^2 と書けると f=g2|f'|=|g|^2 になり, 円周平均をParsevalの等式で係数表示できる。

等号成立条件

平均値公式と三角不等式の等号条件を見るのが最短である。等号が成り立つなら, 円周上の値が同じ向きにそろう。正則関数が開写像であることと両立するのは定数の場合 だけである。

2013年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

(1)

単射正則関数では ff' は円板内で零点を持たない。単位円板は単連結なので, 正則関数 logf\log f' が存在する。従って g=exp(12logf) g=\exp\left(\frac12\log f'\right) とおけば g2=fg^2=f' である。

(2)

g(z)=n=0anzn g(z)=\sum_{n=0}^\infty a_nz^n と書く。曲線 γs(t)=f(seit)\gamma_s(t)=f(se^{it}) の速度は γs(t)=iseitf(seit) \gamma_s'(t)=ise^{it}f'(se^{it}) なので L(s)=s02πf(seit)dt=s02πg(seit)2dt. L(s)=s\int_0^{2\pi}|f'(se^{it})|\,dt =s\int_0^{2\pi}|g(se^{it})|^2\,dt. (3)

Parsevalの等式から L(s)=2πn=0an2s2n+1. L(s)=2\pi\sum_{n=0}^{\infty}|a_n|^2s^{2n+1}. gg は零でないので少なくとも1つの ana_n は非零であり,右辺は 0<s<10<s<1 で狭義単調増加である。

(4)

最後に,正則関数 hh について h(0)=12π02πh(reiθ)dθ h(0)=\frac1{2\pi}\int_0^{2\pi}h(re^{i\theta})\,d\theta である。仮定はこの平均の三角不等式で等号が成立していることを意味する。 従って各 0<r<δ0<r<\delta について,ほとんどすべての θ\thetah(reiθ)h(re^{i\theta}) は同じ偏角を持つ。連続性から円周上で hh の値は同一直線上に 入り,開写像定理により hhz<δ|z|<\delta で定数である。恒等定理より 円板全体で定数である。

最終答

L(s)=2πn=0an2s2n+1 \boxed{L(s)=2\pi\sum_{n=0}^{\infty}|a_n|^2s^{2n+1}} であり,これは ss の狭義単調増加関数である。また仮定を満たす正則関数 hh定数 \boxed{\text{定数}} である。

10 — 変分法とカテナリー

境界条件の使い方

φX0\varphi\in X_0 なので φ(±1)=0\varphi(\pm1)=0 である。第一変分の φ\varphi' を含む項は部分積分で φ\varphi の係数に移せる。

候補の見つけ方

この汎関数は回転面の面積に対応し,停留曲線はカテナリーになる。 λ=cosh1\lambda=\cosh1 と与えられているため,coshxcosh1\cosh x-\cosh1 が端点条件を ぴったり満たす。

2013年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

(1)

F[u]=11(u+λ)1+(u)2dx F[u]=\int_{-1}^{1}(u+\lambda)\sqrt{1+(u')^2}\,dx とおく。u,φX0u,\varphi\in X_0 に対して微分すると ddhF[u+hφ]h=0=11{φ1+(u)2+(u+λ)uφ1+(u)2}dx. \left.\frac{d}{dh}F[u+h\varphi]\right|_{h=0} =\int_{-1}^{1} \left\{\varphi\sqrt{1+(u')^2} +(u+\lambda)\frac{u'\varphi'}{\sqrt{1+(u')^2}}\right\}\,dx. 従って求める極限は上の式で与えられる。

(2)

uuC2C^2 級で,任意の φX0\varphi\in X_0 に対して第一変分が0になるとする。 端点で φ=0\varphi=0 なので部分積分により 11[1+(u)2ddx{(u+λ)u1+(u)2}]φdx=0. \int_{-1}^{1} \left[ \sqrt{1+(u')^2} -\frac{d}{dx}\left\{(u+\lambda)\frac{u'}{\sqrt{1+(u')^2}}\right\} \right]\varphi\,dx=0. よってEuler方程式は ddx{(u+λ)u1+(u)2}=1+(u)2. \frac{d}{dx}\left\{(u+\lambda)\frac{u'}{\sqrt{1+(u')^2}}\right\} =\sqrt{1+(u')^2}.

(3)

λ=cosh1=(e+e1)/2\lambda=\cosh1=(e+e^{-1})/2 であることに注意する。 u(x)=coshxcosh1 u(x)=\cosh x-\cosh1 とおけば u(±1)=0u(\pm1)=0 であり,u+λ=coshxu+\lambda=\cosh xu=sinhxu'=\sinh x である。 したがって 1+(u)2=coshx,(u+λ)u1+(u)2=sinhx \sqrt{1+(u')^2}=\cosh x,\qquad (u+\lambda)\frac{u'}{\sqrt{1+(u')^2}}=\sinh x となり,Euler方程式を満たす。

最終答

第一変分は 11{φ1+(u)2+(u+λ)uφ1+(u)2}dx \boxed{\int_{-1}^{1} \left\{\varphi\sqrt{1+(u')^2} +(u+\lambda)\frac{u'\varphi'}{\sqrt{1+(u')^2}}\right\}\,dx} である。停留関数は ddx{(u+λ)u1+(u)2}=1+(u)2 \boxed{ \frac{d}{dx}\left\{(u+\lambda)\frac{u'}{\sqrt{1+(u')^2}}\right\} =\sqrt{1+(u')^2}} を満たす。例として u(x)=coshxcosh1 \boxed{u(x)=\cosh x-\cosh1} が取れる。

11 — 長方形のDirichlet固有値

無理性の役割

b2b^2 が無理数であるため,異なる (m,n)(m,n) が同じ固有値を与えることがない。 固有関数の形そのものは通常の長方形のDirichlet問題と同じである。

漸近の見方

固有値の数え上げは楕円内の格子点数に帰着する。主項は境界ではなく面積で決まるので, 第1象限の楕円面積だけを計算すればよい。

2013年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

(1)

変数分離 u(x,y)=v(x)w(y)u(x,y)=v(x)w(y) を代入すると vvww=λ -\frac{v''}{v}-\frac{w''}{w}=\lambda である。境界条件から v(0)=v(1)=0,w(0)=w(b)=0 v(0)=v(1)=0,\qquad w(0)=w(b)=0 なので vm(x)=sin(mπx),wn(y)=sin(nπyb) v_m(x)=\sin(m\pi x),\qquad w_n(y)=\sin\left(\frac{n\pi y}{b}\right) が得られる。対応する固有値は λm,n=π2(m2+n2b2)(m,n1) \lambda_{m,n}=\pi^2\left(m^2+\frac{n^2}{b^2}\right)\qquad (m,n\ge1) である。

(2)

一般の固有関数 uu は,xx 方向の正弦級数で u(x,y)=m1um(y)sin(mπx) u(x,y)=\sum_{m\ge1}u_m(y)\sin(m\pi x) と展開できる。方程式に代入すると各 umu_mum+π2m2um=λum,um(0)=um(b)=0 -u_m''+\pi^2m^2u_m=\lambda u_m,\qquad u_m(0)=u_m(b)=0 を満たす。従って非零成分は上の wnw_n に限られる。b2b^2 が無理数なので λm,n\lambda_{m,n} の重複も起きず,固有関数はこれらで尽くされる。

(3)

数え上げ関数は N(μ)=#{(m,n)N2π2(m2+n2b2)μ} N(\mu)=\#\left\{(m,n)\in\mathbb N^2\mid \pi^2\left(m^2+\frac{n^2}{b^2}\right)\le\mu\right\} である。これは第1象限内の楕円 m2+n2b2μπ2 m^2+\frac{n^2}{b^2}\le \frac{\mu}{\pi^2} の格子点数であり,面積近似から N(μ)14πμπbμπ=b4πμ. N(\mu)\sim \frac14\pi\cdot\frac{\sqrt\mu}{\pi}\cdot\frac{b\sqrt\mu}{\pi} =\frac{b}{4\pi}\mu.

最終答

um,n(x,y)=sin(mπx)sin(nπyb) \boxed{u_{m,n}(x,y)=\sin(m\pi x)\sin\left(\frac{n\pi y}{b}\right)} λm,n=π2(m2+n2b2)(m,n1) \boxed{\lambda_{m,n}=\pi^2\left(m^2+\frac{n^2}{b^2}\right)\quad (m,n\ge1)} で尽くされる。また α=1,limμN(μ)μ=b4π \boxed{\alpha=1,\qquad \lim_{\mu\to\infty}\frac{N(\mu)}{\mu}=\frac{b}{4\pi}} である。

12 — 重み付き2乗可積分条件の閉性

有界閉球ではない点

L2L^2 ノルムだけでなく,xux uL2L^2 ノルムも制御している。 しかし L2L^2 収束から部分列のほとんど至る所収束を取り,重み付きの非負関数に Fatouを使えば,下半連続性がそのまま得られる。

2013年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

MMM={uL2(R)Ru2dμ+Rx2u2dμ1} M=\left\{u\in L^2(\mathbb R)\mid \int_{\mathbb R}|u|^2\,d\mu+\int_{\mathbb R}x^2|u|^2\,d\mu\le1\right\} とする。unMu_n\in Munuu_n\to u in L2(R)L^2(\mathbb R) と仮定する。 部分列を取れば unj(x)u(x)u_{n_j}(x)\to u(x) がほとんど至る所で成り立つ。

非負関数 (1+x2)unj(x)2 (1+x^2)|u_{n_j}(x)|^2 にFatouの補題を適用すると R(1+x2)u(x)2dμlim infjR(1+x2)unj(x)2dμ1. \int_{\mathbb R}(1+x^2)|u(x)|^2\,d\mu \le \liminf_{j\to\infty} \int_{\mathbb R}(1+x^2)|u_{n_j}(x)|^2\,d\mu \le1. 従って uMu\in M である。

最終答

M は L2(R) の閉部分集合である。 \boxed{M\text{ は }L^2(\mathbb R)\text{ の閉部分集合である。}}

13 — 感染症型常微分方程式

R0R_0 の意味

R0R_0 は,感染成分 (y,z)(y,z) の線形化で増殖が減衰に勝つかを表す量である。 R0>1R_0>1 なら正の平衡点が現れ,R0<1R_0<1 なら病気のない平衡点だけが残る。

消滅証明の重み

yyzz をそのまま足すと zz の項が消えにくい。係数 β2λ/(με)\beta_2\lambda/(\mu\varepsilon) を選ぶと zz の正項を相殺でき, Wγ(R01)yW'\le \gamma(R_0-1)y という単純な評価になる。

2013年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

(1)

Ω\Omega の各境界でベクトル場の向きを調べる。x=0x=0 では x=λ>0, x'=\lambda>0, y=0y=0 では y=0y'=0z=0z=0 では z=δy0z'=\delta y\ge0 である。また (x+y)=λμxγyλμ(x+y) (x+y)'=\lambda-\mu x-\gamma y\le \lambda-\mu(x+y) であり,γ>μ\gamma>\mu から x+y=λ/μx+y=\lambda/\mu の面では内向きである。 さらに z=δλ/(με)z=\delta\lambda/(\mu\varepsilon) では z=δyεzδλμδλμ=0. z'=\delta y-\varepsilon z\le \frac{\delta\lambda}{\mu} -\frac{\delta\lambda}{\mu}=0. 従って Ω\Omega は正に不変である。

(2)

平衡点を求める。y=z=0y=z=0 なら E0=(λμ,0,0) E_0=\left(\frac{\lambda}{\mu},0,0\right) を得る。内部平衡点では z=0z' =0 から z=(δ/ε)yz=(\delta/\varepsilon)yy=0y'=0 から x(β1+δβ2ε)=γ. x\left(\beta_1+\frac{\delta\beta_2}{\varepsilon}\right)=\gamma. A=β1+δβ2ε A=\beta_1+\frac{\delta\beta_2}{\varepsilon} とおけば x=γA,y=λμxγ,z=δεy. x_*=\frac{\gamma}{A},\qquad y_*=\frac{\lambda-\mu x_*}{\gamma},\qquad z_*=\frac{\delta}{\varepsilon}y_*. y>0y_*>0 となる条件は λAμγ>1, \frac{\lambda A}{\mu\gamma}>1, すなわち R0>1R_0>1 である。従って R0>1R_0>1 なら内部平衡点はただ1つである。

(3)

R0<1R_0<1 では内部平衡点はなく,境界平衡点は E0E_0 のみである。 (y,z)(y,z)-成分に関する E0E_0 での線形化は (β1λ/μγβ2λ/μδε). \begin{pmatrix} \beta_1\lambda/\mu-\gamma & \beta_2\lambda/\mu\\ \delta&-\varepsilon \end{pmatrix}. この行列の行列式は γε(1R0)>0 \gamma\varepsilon(1-R_0)>0 で,跡は負である。従って E0E_0 は局所漸近安定である。

(4)

最後に大域的な消滅を示す。x(t)λ/μx(t)\le\lambda/\mu なので W(t)=y(t)+β2λμεz(t) W(t)=y(t)+\frac{\beta_2\lambda}{\mu\varepsilon}z(t) とおくと W(t)γ(R01)y(t)0. W'(t)\le \gamma(R_0-1)y(t)\le0. よって WW は減少し,0y(t)dt<\int_0^\infty y(t)\,dt<\infty である。 方程式から yy' は有界なので y(t)0y(t)\to0。さらに z=δyεz z'=\delta y-\varepsilon z より z(t)0z(t)\to0 である。

最終答

Ω\Omega は正に不変である。R0>1R_0>1 では内部平衡点 x=γβ1+δβ2/ε,y=λμxγ,z=δεy \boxed{ x_*=\frac{\gamma}{\beta_1+\delta\beta_2/\varepsilon},\quad y_*=\frac{\lambda-\mu x_*}{\gamma},\quad z_*=\frac{\delta}{\varepsilon}y_*} がただ1つ存在する。R0<1R_0<1 では境界平衡点 (λ/μ,0,0) \boxed{\left(\lambda/\mu,0,0\right)} のみが存在し,これは局所漸近安定で,さらに limty(t)=0,limtz(t)=0 \boxed{\lim_{t\to\infty}y(t)=0,\qquad \lim_{t\to\infty}z(t)=0} である。

14 — 有限遷移系の到達性

AFの読み替え

有限遷移系では「どの無限経路もいつか XX に入る」は「XX を避け続ける 有向閉路がない」と同じである。この読み替えをすると,問題は有向閉路を壊す頂点集合 の問題になる。

下界と構成

m1m_1 の下界は4本の行閉路で出る。上界は対角的に4点を取って,残りに明示的な トポロジカル順序を与えればよい。第2問は最大値そのものではなく 88 以上を示せば 十分なので,8点の具体例を作るのが速い。

2013年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

(1)

まず AF(X)=QAF(X)=Q となることは,QXQ\setminus X の中だけを通る無限遷移列が 存在しないことと同値である。有限グラフでは,これは QXQ\setminus X が有向閉路を 含まないことと同値である。

固定した mm に対し,4点 (m,0),(m,1),(m,2),(m,3) (m,0),(m,1),(m,2),(m,3) は有向閉路を作る。4つの行は互いに交わらないので,すべての有向閉路を壊すには 少なくとも4点を XX に入れる必要がある。

一方 X0={(0,3),(1,2),(2,1),(3,0)} X_0=\{(0,3),(1,2),(2,1),(3,0)\} とおくと,QX0Q\setminus X_0 は次の順序で並べられる。 (0,0),(1,1),(2,2),(3,3),(1,0),(0,1),(3,2),(2,3),(3,1),(0,2),(2,0),(1,3). \begin{gathered} (0,0),(1,1),(2,2),(3,3),\\ (1,0),(0,1),(3,2),(2,3),\\ (3,1),(0,2),(2,0),(1,3). \end{gathered} この順序で各遷移は右へ進むか,または X0X_0 に入る。従って QX0Q\setminus X_0 には有向閉路がない。よって m1=4. m_1=4.

(2)

次に X1={(0,0),(0,2),(0,3),(1,1),(1,2),(1,3),(2,1),(2,2)} X_1=\{(0,0),(0,2),(0,3),(1,1),(1,2),(1,3),(2,1),(2,2)\} を取る。このとき AF(X1)=X1AF(X_1)=X_1 であり,X1X_1 内の遷移は (0,2)(0,3)(0,0),(1,3)(1,2)(1,1),(2,1)(2,2) (0,2)\to(0,3)\to(0,0),\quad (1,3)\to(1,2)\to(1,1),\quad (2,1)\to(2,2) およびこれらから外へ出るものだけで,有向閉路を持たない。従って EG(AF(X1))=. EG(AF(X_1))=\varnothing. この X1X_1 は8点からなるので M28. M_2\ge8.

最終答

m1=4,M28 \boxed{m_1=4},\qquad \boxed{M_2\ge8} である。

15 — 有限離散測度の直交多項式

行列式表示の意味

この QnQ_n はGram--Schmidtで作った直交多項式を,正規化せず行列式で書いたもの である。最後の列を xmx^m と内積に入れると,同じ列が2本できるため直交性が出る。

有限台の影響

台が N+1N+1 点しかないので,次数 N+1N+1 以上では通常の無限直交多項式列とは 違う挙動をする。QN+1Q_{N+1} は台の全点で消える多項式で,それより先は行列式が 恒等的に消える。

2013年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

(1)

p,q=k=0Np(k)q(k)wk\langle p,q\rangle=\sum_{k=0}^{N}p(k)q(k)w_k とする。 行列式で定義された QQ_\ell の最後の列を xmx^m と内積に入れると, m<m<\ell では最後の列が既存のモーメント列と一致する。従って行列式が0になり, Q,xm=0(m<) \langle Q_\ell,x^m\rangle=0\qquad (m<\ell) である。

(2)

Δn\Delta_n1,x,,xn1,x,\dots,x^n のGram行列式である。Gram--Schmidtの標準公式 から,QQ_\ell の最高次係数は Δ1\Delta_{\ell-1} であり, Q,Q=ΔΔ1(1N) \langle Q_\ell,Q_\ell\rangle=\Delta_\ell\Delta_{\ell-1}\qquad (1\le \ell\le N) となる。また =0\ell=0 では Q0,Q0=Δ0=1\langle Q_0,Q_0\rangle=\Delta_0=1 である。

(3)

xQnxQ_n は次数 n+1n+1 なので,直交性により xQn=αnQn+1+βnQn+γnQn1 xQ_n=\alpha_nQ_{n+1}+\beta_nQ_n+\gamma_nQ_{n-1} と書ける。最高次係数を比べると αn=Δn1Δn. \alpha_n=\frac{\Delta_{n-1}}{\Delta_n}. さらに Qn1Q_{n-1} と内積を取ると γnQn1,Qn1=xQn,Qn1=Qn,xQn1=αn1Qn,Qn. \gamma_n\langle Q_{n-1},Q_{n-1}\rangle =\langle xQ_n,Q_{n-1}\rangle =\langle Q_n,xQ_{n-1}\rangle =\alpha_{n-1}\langle Q_n,Q_n\rangle. 従って γn=ΔnΔn1. \gamma_n=\frac{\Delta_n}{\Delta_{n-1}}.

(4)

最後に,測度は 0,1,,N0,1,\dots,NN+1N+1 点だけに台を持つ。したがって k=0N(xk) \prod_{k=0}^{N}(x-k) は全ての多項式と内積0である。QN+1Q_{N+1} はこの多項式の非零定数倍であるため, その零点は 0,1,,N 0,1,\dots,N である。さらに 2\ell\ge2 ではGram行列の階数が足りず, QN+Q_{N+\ell} は零多項式になる。

最終答

Q,xm=0(m<) \boxed{\langle Q_\ell,x^m\rangle=0\quad(m<\ell)} Q,Q=ΔΔ1(1N) \boxed{\langle Q_\ell,Q_\ell\rangle=\Delta_\ell\Delta_{\ell-1}\quad(1\le \ell\le N)} γn=ΔnΔn1 \boxed{\gamma_n=\frac{\Delta_n}{\Delta_{n-1}}} また QN+1Q_{N+1} の零点は 0,1,,N\boxed{0,1,\dots,N} であり, 2\ell\ge2 では QN+Q_{N+\ell} は零多項式である。

16 — ランダムな添字付き和の極限定理

非零になる添字の組

XjX_j は平均0で独立なので,積の期待値が残るには各添字が少なくとも2回現れる 必要がある。jkj\ne k では,対応する2通りの入れ替わりだけが寄与する。

2つの極限の違い

qn0q_n\to0 ではほとんど全ての参照先が0に集中し,共通因子 X0X_0 が残る。 一方 qn=1n3q_n=1-n^{-3} では参照先が遠方に逃げ,しかも衝突しないため,通常の 独立和の中心極限定理に戻る。

2013年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

(1)

jkj\ne k とする。Njn,NknN_j^n,N_k^n の値で条件付けると E[XjXNjnXkXNkn] E[X_jX_{N_j^n}X_kX_{N_k^n}] が非零になり得るのは,添字が (Njn,Nkn)=(k,j)または(j,k) (N_j^n,N_k^n)=(k,j)\quad\text{または}\quad (j,k) となる場合だけである。従って E ⁣[XjXNjnXkXNkn1{Njn0, Nkn0}]=2v2(1qn)2qnj+k. E\!\left[X_jX_{N_j^n}X_kX_{N_k^n} 1_{\{N_j^n\ne0,\ N_k^n\ne0\}}\right] =2v^2(1-q_n)^2q_n^{j+k}.

(2)

qn0q_n\to0 とする。このとき Njn=0N_j^n=0 である確率が1に近づくので YnX0nj=1nXj Y_n-\frac{X_0}{\sqrt n}\sum_{j=1}^{n}X_j は確率収束で0に近づく。中心極限定理より 1nj=1nXjZ,ZN(0,v), \frac1{\sqrt n}\sum_{j=1}^{n}X_j\Rightarrow Z,\qquad Z\sim N(0,v), かつ ZZX0X_0 と独立に取れる。よって YnX0Z. Y_n\Rightarrow X_0Z.

(3)

次に qn=1n3q_n=1-n^{-3} とする。 P(Njnn)=1qnn+1n+1n3 P(N_j^n\le n)=1-q_n^{n+1}\sim \frac{n+1}{n^3} だから,和集合評価により P(min1jnNjnn)n{1qnn+1}0. P\left(\min_{1\le j\le n}N_j^n\le n\right) \le n\{1-q_n^{n+1}\}\longrightarrow0.

(4)

この場合,高確率で全ての NjnN_j^nnn より大きい。さらに P(Njn=Nkn)=m0(1qn)2qn2m=1qn1+qn=O(n3) P(N_j^n=N_k^n)=\sum_{m\ge0}(1-q_n)^2q_n^{2m} =\frac{1-q_n}{1+q_n}=O(n^{-3}) なので,1j<kn1\le j<k\le n の衝突確率の総和は O(n1)O(n^{-1}) で0に収束する。 従って高確率で XjXNjn(j=1,,n) X_jX_{N_j^n}\qquad (j=1,\dots,n) は互いに独立で,同分布,平均0,分散 v2v^2 の変数として扱える。中心極限定理から YnN(0,v2). Y_n\Rightarrow N(0,v^2).

最終答

E[XjXNjnXkXNkn1{Njn0,Nkn0}]=2v2(1qn)2qnj+k(jk) \boxed{ E[X_jX_{N_j^n}X_kX_{N_k^n}1_{\{N_j^n\ne0,N_k^n\ne0\}}] =2v^2(1-q_n)^2q_n^{j+k}\quad(j\ne k)} qn0 なら YnX0Z,ZN(0,v) \boxed{q_n\to0\text{ なら }Y_n\Rightarrow X_0Z,\quad Z\sim N(0,v)} qn=1n3 なら min1jnNjn>n が高確率で起こり,YnN(0,v2) \boxed{q_n=1-n^{-3}\text{ なら }\min_{1\le j\le n}N_j^n>n\text{ が高確率で起こり,} Y_n\Rightarrow N(0,v^2)} である。

17 — 対称化された有理式の双対性

対称化の係数

分割 I,JI,J を固定すると,II の中の並べ替えが m!m! 通り,JJ の中の 並べ替えが n!n! 通りある。これが係数 m!n!/(m+n)!m!n!/(m+n)! の由来である。

基本恒等式から全体へ

この種の問題は,まず最小の非自明な交換 h1,2=h2,1h_{1,2}=h_{2,1} を直接確認し, それを大きな対称化和の中に埋め込んで使う。全体の双対性は,多数の変数を一度に 動かすのではなく,1個ずつブロックの反対側へ移すことで得られる。

2013年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

(1)

定義から hm,n=1(m+n)!σSm+n1im<m+1jm+nλ(xσ(i),xσ(j)). h_{m,n}=\frac1{(m+n)!}\sum_{\sigma\in S_{m+n}} \prod_{1\le i\le m<m+1\le j\le m+n} \lambda(x_{\sigma(i)},x_{\sigma(j)}). 順列 σ\sigma が,前半 mm 個の添字集合 II と後半 nn 個の添字集合 JJ を定めると,同じ分割 (I,J)(I,J)m!n!m!n! 回現れる。従って hm,n=m!n!(m+n)!IJ={1,,m+n}I=m, J=niI, jJλ(xi,xj). h_{m,n} =\frac{m!n!}{(m+n)!} \sum_{\substack{I\sqcup J=\{1,\dots,m+n\}\\ |I|=m,\ |J|=n}} \prod_{i\in I,\ j\in J}\lambda(x_i,x_j).

(2)

次に λ(x,y)=λ(ay,x) \lambda(x,y)=\lambda(ay,x) を使って h1,2h_{1,2}h2,1h_{2,1} を共通分母で比較する。差の分子は (ab)(1a)(1b)(x1x2)(x2x3)(x3x1) (a-b)(1-a)(1-b)(x_1-x_2)(x_2-x_3)(x_3-x_1) を因子に持つ交代式である。一方,差全体は対称式であるから,この交代因子を持つ 対称式は零でなければならない。従って h1,2=h2,1 h_{1,2}=h_{2,1} である。

(3)

この基本恒等式を3変数ずつに適用すると,1個の変数を右側から左側へ移す操作が 可能になる。具体的には,対称化の和を,選んだ3変数と残りの変数に分け, 選んだ3変数に h1,2=h2,1h_{1,2}=h_{2,1} を適用する。これにより h1,n=hn,1(n3) h_{1,n}=h_{n,1}\qquad (n\ge3) が従う。

(4)

一般の m,nm,n については,隣り合うブロックサイズの交換を繰り返す。上の 11nn の交換を対称化の中で適用すれば hm,n=hm1,n+1 h_{m,n}=h_{m-1,n+1} 型の移動ができるので,これを有限回繰り返して hm,n=hn,m h_{m,n}=h_{n,m} を得る。

最終答

分割表示は hm,n=m!n!(m+n)!IJ={1,,m+n}I=m, J=niI, jJλ(xi,xj) \boxed{ h_{m,n} =\frac{m!n!}{(m+n)!} \sum_{\substack{I\sqcup J=\{1,\dots,m+n\}\\ |I|=m,\ |J|=n}} \prod_{i\in I,\ j\in J}\lambda(x_i,x_j)} である。また h1,2=h2,1,h1,n=hn,1 (n3),hm,n=hn,m \boxed{h_{1,2}=h_{2,1}},\qquad \boxed{h_{1,n}=h_{n,1}\ (n\ge3)},\qquad \boxed{h_{m,n}=h_{n,m}} が成り立つ。

18 — 半線形熱方程式の陽的差分法

安定性の入口

0<λ1/20<\lambda\le1/2 により,線形熱方程式部分は隣接点の凸結合になる。 そのため最大ノルムで増幅しない。この安定性に非線形項と打ち切り誤差を足して 評価する。

非線形項の扱い

U2u2=(Uu)(U+u)U^2-u^2=(U-u)(U+u) と分解し,uK|u|\le K を使うと Ek(Ek+2K)E_k(E_k+2K) が出る。あとは離散不等式を常微分方程式で上から比較すれば, 収束は y(t)=O(h2)y(t)=O(h^2) から従う。

2013年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

(1)

厳密解を格子点でTaylor展開する。時間方向と空間方向の打ち切り誤差は,定義された RR により u(xi,tk+1)u(xi,tk)τu(xi1,tk)2u(xi,tk)+u(xi+1,tk)h2u(xi,tk)2h2R \left| \frac{u(x_i,t_{k+1})-u(x_i,t_k)}{\tau} -\frac{u(x_{i-1},t_k)-2u(x_i,t_k)+u(x_{i+1},t_k)}{h^2} -u(x_i,t_k)^2 \right| \le h^2R と評価できる。

差分解との差を取ると ei,k+1=λei1,k+(12λ)ei,k+λei+1,k +τ{Ui,k2u(xi,tk)2}+τh2ρi,k e_{i,k+1} =\lambda e_{i-1,k}+(1-2\lambda)e_{i,k}+\lambda e_{i+1,k} \ +\tau\{U_{i,k}^2-u(x_i,t_k)^2\}+\tau h^2\rho_{i,k} となる。ただし ρi,kR|\rho_{i,k}|\le R である。0<λ1/20<\lambda\le1/2 なので線形部分は 凸結合であり, Ui,k2ui2ei,k{ei,k+2ui}Ek(Ek+2K) |U_{i,k}^2-u_i^2| \le |e_{i,k}|\{|e_{i,k}|+2|u_i|\} \le E_k(E_k+2K) である。従って Ek+1Ek+τEk(Ek+2K)+τh2R,E0=0. E_{k+1}\le E_k+\tau E_k(E_k+2K)+\tau h^2R,\qquad E_0=0.

(2)

常微分方程式 y=y(y+2K)+h2R,y(0)=0 y'=y(y+2K)+h^2R,\qquad y(0)=0 を考える。0yK0\le y\le K の範囲では y3Ky+h2R. y'\le 3K y+h^2R. 比較から y(t)h2R3K(e3Kt1) y(t)\le \frac{h^2R}{3K}(e^{3Kt}-1) であり,仮定 e2KTh2R<K2e^{2KT}h^2R<K^2 の下で少なくとも 0<t<T0<t<T では爆発せず, 一意解を持つ。

(3)

関数 F(y)=y(y+2K)+h2RF(y)=y(y+2K)+h^2R は増加関数である。E0=y(0)E_0=y(0) であり,前段の 不等式はEuler型の上方比較になっているので,帰納的に Eky(tk)(0kM) E_k\le y(t_k)\qquad (0\le k\le M) が得られる。

(4)

最後に h0h\to0 とすると,右辺の常微分方程式の解は一様に0へ収束する。 実際,0tT0\le t\le T0y(t)CTh2R 0\le y(t)\le C_T h^2R となる定数 CTC_T が取れる。従って max1iN, 1kMei,kmax0kMy(tk)0. \max_{1\le i\le N,\ 1\le k\le M}|e_{i,k}| \le \max_{0\le k\le M}y(t_k)\longrightarrow0.

最終答

誤差は Ek+1Ek+τEk(Ek+2K)+τh2R,E0=0 \boxed{E_{k+1}\le E_k+\tau E_k(E_k+2K)+\tau h^2R,\qquad E_0=0} を満たす。比較方程式 y=y(y+2K)+h2R,y(0)=0 \boxed{y'=y(y+2K)+h^2R,\qquad y(0)=0} は仮定の下で 0<t<T0<t<T に一意解を持ち, Eky(tk)(0kM) \boxed{E_k\le y(t_k)\quad(0\le k\le M)} である。従って limh0max1iN1kMei,k=0 \boxed{\displaystyle \lim_{h\to0}\max_{\substack{1\le i\le N\\1\le k\le M}}|e_{i,k}|=0} が成り立つ。

東京大学 専門科目B — 他の年度