跳到主要内容

京都大学 情報学研究科 システム科学専攻 2019年8月実施 数学 II

Author

思齐塾, 祭音Myyura

Description

日本語版

問1

xyzxyz -空間上の有界な閉集合

D={(x,y,z):x2a2+y2b2z,x0,y0,0z1}D = \left\{ (x, y, z) : \frac{x^2}{a^2} + \frac{y^2}{b^2} \leq z, x \geq 0, y \geq 0, 0 \leq z \leq 1 \right\}

について考える. ただし, a,ba, b は正の定数であるとする. 以下の設問に答えよ.

(i) 次のように変数 (r,θ,s)(r, \theta, s) を変数 (x,y,z)(x, y, z) に移す写像のヤコビ行列式を求めよ.

x=arcosθ,y=brsinθ,z=sx = ar \cos\theta, \quad y = br \sin\theta, \quad z = s

ただし, r0,0θπ/2r \geq 0, 0 \leq \theta \leq \pi/2 であるとする.

(ii) 設問 (i) の写像によって xyzxyz -空間上の集合 DD に移される rθsr\theta s -空間上の集合は

E={(r,θ,s):0r,0θπ/2,0s1}E = \{ (r, \theta, s) : 0 \leq r \leq \boxed{\quad}, 0 \leq \theta \leq \pi/2, 0 \leq s \leq 1 \}

と表される. このとき, 空欄 \boxed{\quad} に入る式を書け.

(iii) Ddxdydz\displaystyle\iiint_D dxdydz を求めよ.

(iv) 以下の積分が有限となるような正の整数 l,m,nl, m, n の組のうち, l+mnl+m \leq n を満たすものを全て求めよ.

Dxlymzndxdydz\iiint_D \frac{x^l y^m}{z^n} dxdydz

問2

数列 {xn:n=1,2,...}\{x_n : n = 1, 2, ...\} を次のように定義する。

x1=1,xn=sinxn1,n=2,3,...x_1 = 1, \quad x_n = \sin x_{n-1}, \quad n = 2, 3, ...

また、2つの実数値関数 ffgg が, limx0f(x)/g(x)=0\lim_{x \to 0} f(x)/g(x) = 0 を満たすとき、ランダウの記法によって、 f(x)=o(g(x))f(x) = o(g(x)) と表記する.ただし, f(x)=o(1)f(x) = o(1) は, limx0f(x)=0\lim_{x \to 0} f(x) = 0 を意味するものとする。以下の設問に答えよ。

(i) 次の不等式が成り立つことを(帰納法などにより)示せ。

0<xn+1<xn,n=1,2,...0 < x_{n+1} < x_n, \quad n = 1, 2, ...

(ii) limnxn\lim_{n \to \infty} x_n の値を求めよ。

(iii) sinx\sin x のマクローリン展開を書け。ただし、 xx の5次以上の項は、ランダウの記法を利用してまとめて表記せよ。

(iv) 次式を満たす定数 aabb を求めよ。

1(sinx)2=ax2+b+o(1)\frac{1}{(\sin x)^2} = \frac{a}{x^2} + b + o(1)

(v) 設問 (ii) と (iv) の結果を利用して, limn(1xn+121xn2)\lim_{n \to \infty} \left(\frac{1}{x_{n+1}^2} - \frac{1}{x_n^2}\right) の値を導け。

(vi) 設問 (v) の結果を利用して、次式が成り立つことを示せ。

limnnxn2=3\lim_{n \to \infty} nx_n^2 = 3

题目描述

問1

a,ba,b 为正常数,考虑 xyzxyz 空间中的有界闭集

D={(x,y,z) | x2a2+y2b2z,x0, y0, 0z1}.D= \left\{ (x,y,z) \ \middle|\ \frac{x^2}{a^2}+\frac{y^2}{b^2}\leq z,\quad x\geq0,\ y\geq0,\ 0\leq z\leq1 \right\}.

回答:

  1. 对从 (r,θ,s)(r,\theta,s)(x,y,z)(x,y,z) 的变量变换
x=arcosθ,y=brsinθ,z=s,x=ar\cos\theta, \qquad y=br\sin\theta, \qquad z=s,

其中 r0r\geq00θπ/20\leq\theta\leq\pi/2,求 Jacobian 行列式

(x,y,z)(r,θ,s).\frac{\partial(x,y,z)} {\partial(r,\theta,s)}.
  1. 在上述映射下,被映到 DDrθsr\theta s 空间区域可写为
E={(r,θ,s) | 0rs, 0θπ2, 0s1}.E= \left\{ (r,\theta,s) \ \middle|\ 0\leq r\leq\boxed{\phantom{\sqrt{s}}}, \ 0\leq\theta\leq\frac{\pi}{2}, \ 0\leq s\leq1 \right\}.

写出方框中应填的表达式。 3. 计算 DD 的体积

Ddxdydz.\iiint_D dx\,dy\,dz.
  1. 在使广义积分
Dxlymzndxdydz\iiint_D \frac{x^l y^m}{z^n}\,dx\,dy\,dz

有限的正整数三元组 (l,m,n)(l,m,n) 中,列出所有还满足 l+mnl+m\leq n 的三元组。

問2

定义数列 {xn:n=1,2,}\{x_n:n=1,2,\ldots\}

x1=1,xn=sinxn1(n=2,3,).x_1=1, \qquad x_n=\sin x_{n-1} \quad(n=2,3,\ldots).

对两个实值函数 f,gf,g,若

limx0f(x)g(x)=0,\lim_{x\to0}\frac{f(x)}{g(x)}=0,

则用 Landau 记号写作 f(x)=o(g(x))f(x)=o(g(x));特别地, f(x)=o(1)f(x)=o(1) 表示 limx0f(x)=0\lim_{x\to0}f(x)=0。回答:

  1. 证明(可使用归纳法)
0<xn+1<xn,n=1,2,.0<x_{n+1}<x_n, \qquad n=1,2,\ldots.
limnxn.\lim_{n\to\infty}x_n.
  1. 写出 sinx\sin x 的 Maclaurin 展开,并用上述 Landau 记号统一表示五次及更高次的项。
  2. 求常数 a,ba,b,使
1(sinx)2=ax2+b+o(1)(x0).\frac{1}{(\sin x)^2} = \frac{a}{x^2}+b+o(1) \qquad(x\to0).
  1. 利用第 2、4 问的结果求
limn(1xn+121xn2).\lim_{n\to\infty} \left( \frac{1}{x_{n+1}^2} -\frac{1}{x_n^2} \right).
  1. 利用第 5 问证明
limnnxn2=3.\lim_{n\to\infty}nx_n^2=3.

Kai

問1

(i) ヤコビ行列式の計算

与えられた変数変換は

x=arcosθ,y=brsinθ,z=sx = ar \cos\theta, \quad y = br \sin\theta, \quad z = s

ヤコビ行列 JJ は、各変数の偏導関数から構成される行列です。

J=(x,y,z)(r,θ,s)=(xrxθxsyryθyszrzθzs)J = \frac{\partial(x, y, z)}{\partial(r, \theta, s)} = \begin{pmatrix} \frac{\partial x}{\partial r} & \frac{\partial x}{\partial \theta} & \frac{\partial x}{\partial s} \\ \frac{\partial y}{\partial r} & \frac{\partial y}{\partial \theta} & \frac{\partial y}{\partial s} \\ \frac{\partial z}{\partial r} & \frac{\partial z}{\partial \theta} & \frac{\partial z}{\partial s} \end{pmatrix}

各偏導関数を計算します。

xr=acosθ,xθ=arsinθ,xs=0\frac{\partial x}{\partial r} = a \cos\theta, \quad \frac{\partial x}{\partial \theta} = -ar \sin\theta, \quad \frac{\partial x}{\partial s} = 0
yr=bsinθ,yθ=brcosθ,ys=0\frac{\partial y}{\partial r} = b \sin\theta, \quad \frac{\partial y}{\partial \theta} = br \cos\theta, \quad \frac{\partial y}{\partial s} = 0
zr=0,zθ=0,zs=1\frac{\partial z}{\partial r} = 0, \quad \frac{\partial z}{\partial \theta} = 0, \quad \frac{\partial z}{\partial s} = 1

よって、ヤコビ行列は

J=(acosθarsinθ0bsinθbrcosθ0001)J = \begin{pmatrix} a \cos\theta & -ar \sin\theta & 0 \\ b \sin\theta & br \cos\theta & 0 \\ 0 & 0 & 1 \end{pmatrix}

ヤコビ行列式(ヤコビアン)は、この行列の行列式です。

det(J)=1det(acosθarsinθbsinθbrcosθ)=(acosθ)(brcosθ)(arsinθ)(bsinθ)\det(J) = 1 \cdot \det \begin{pmatrix} a \cos\theta & -ar \sin\theta \\ b \sin\theta & br \cos\theta \end{pmatrix} = (a \cos\theta)(br \cos\theta) - (-ar \sin\theta)(b \sin\theta)
=abrcos2θ+abrsin2θ=abr(cos2θ+sin2θ)=abr= abr \cos^2\theta + abr \sin^2\theta = abr(\cos^2\theta + \sin^2\theta) = abr

よって、ヤコビ行列式は abrabr です。

(ii) 新しい積分領域 E の決定

集合 DD の定義域を (r,θ,s)(r, \theta, s) 座標系に変換します。 DD の条件は以下の通りです。

  1. x2a2+y2b2z\frac{x^2}{a^2} + \frac{y^2}{b^2} \leq z
  2. x0,y0x \geq 0, y \geq 0
  3. 0z10 \leq z \leq 1

変数変換 x=arcosθ,y=brsinθ,z=sx = ar \cos\theta, y = br \sin\theta, z = s を代入します。

  1. (arcosθ)2a2+(brsinθ)2b2s    a2r2cos2θa2+b2r2sin2θb2s    r2(cos2θ+sin2θ)s    r2s\frac{(ar \cos\theta)^2}{a^2} + \frac{(br \sin\theta)^2}{b^2} \leq s \implies \frac{a^2 r^2 \cos^2\theta}{a^2} + \frac{b^2 r^2 \sin^2\theta}{b^2} \leq s \implies r^2(\cos^2\theta + \sin^2\theta) \leq s \implies r^2 \leq s . r0r \geq 0 なので, 0rs0 \leq r \leq \sqrt{s} となります。
  2. x=arcosθ0x = ar \cos\theta \geq 0y=brsinθ0y = br \sin\theta \geq 0 . a,b,ra, b, r は非負なので, cosθ0\cos\theta \geq 0 かつ sinθ0\sin\theta \geq 0 が必要です。これは与えられた条件 0θπ/20 \leq \theta \leq \pi/2 を満たします。
  3. 0z1    0s10 \leq z \leq 1 \implies 0 \leq s \leq 1 .

これらの条件をまとめると、 rθsr\theta s -空間上の集合 EE

E={(r,θ,s):0rs,0θπ/2,0s1}E = \{ (r, \theta, s) : 0 \leq r \leq \sqrt{s}, 0 \leq \theta \leq \pi/2, 0 \leq s \leq 1 \}

したがって、空欄に入る式は s\sqrt{s} です。

(iii) 体積の計算

求める体積は、領域 DD 上での 11 の三重積分です。変数変換を用いて計算します。

Ddxdydz=Edet(J)drdθds=Eabrdrdθds\iiint_D dxdydz = \iiint_E |\det(J)| dr d\theta ds = \iiint_E abr dr d\theta ds

a,b>0,r0a,b>0, r \ge 0 なので abr=abr|abr| = abr です。積分範囲を用いて計算します。

V=010π/20sabrdrdθdsV = \int_0^1 \int_0^{\pi/2} \int_0^{\sqrt{s}} abr \, dr \, d\theta \, ds

まず rr で積分します。

0sabrdr=ab[r22]0s=ab(s)22=ab2s\int_0^{\sqrt{s}} abr \, dr = ab \left[ \frac{r^2}{2} \right]_0^{\sqrt{s}} = ab \frac{(\sqrt{s})^2}{2} = \frac{ab}{2}s

次に θ\theta で積分します。

0π/2ab2sdθ=ab2s[θ]0π/2=ab2sπ2=πab4s\int_0^{\pi/2} \frac{ab}{2}s \, d\theta = \frac{ab}{2}s [\theta]_0^{\pi/2} = \frac{ab}{2}s \cdot \frac{\pi}{2} = \frac{\pi ab}{4}s

最後に ss で積分します。

01πab4sds=πab4[s22]01=πab412=πab8\int_0^1 \frac{\pi ab}{4}s \, ds = \frac{\pi ab}{4} \left[ \frac{s^2}{2} \right]_0^1 = \frac{\pi ab}{4} \cdot \frac{1}{2} = \frac{\pi ab}{8}

よって、体積は πab8\frac{\pi ab}{8} です。

(iv) 積分の収束条件

与えられた積分を II とします。

I=DxlymzndxdydzI = \iiint_D \frac{x^l y^m}{z^n} dxdydz

変数変換を行うと、

I=E(arcosθ)l(brsinθ)msn(abr)drdθdsI = \iiint_E \frac{(ar \cos\theta)^l (br \sin\theta)^m}{s^n} (abr) \, dr d\theta ds
I=al+1bm+1Erl+m+1coslθsinmθsndrdθdsI = a^{l+1}b^{m+1} \iiint_E \frac{r^{l+m+1} \cos^l\theta \sin^m\theta}{s^n} \, dr d\theta ds

積分を分離します。

I=al+1bm+1(0π/2coslθsinmθdθ)(011sn(0srl+m+1dr)ds)I = a^{l+1}b^{m+1} \left( \int_0^{\pi/2} \cos^l\theta \sin^m\theta \, d\theta \right) \left( \int_0^1 \frac{1}{s^n} \left( \int_0^{\sqrt{s}} r^{l+m+1} \, dr \right) ds \right)

l,ml, m は正の整数なので、 θ\theta の積分は有限値になります。 rr の積分を計算します。

0srl+m+1dr=[rl+m+2l+m+2]0s=(s)l+m+2l+m+2=sl+m+22l+m+2\int_0^{\sqrt{s}} r^{l+m+1} \, dr = \left[ \frac{r^{l+m+2}}{l+m+2} \right]_0^{\sqrt{s}} = \frac{(\sqrt{s})^{l+m+2}}{l+m+2} = \frac{s^{\frac{l+m+2}{2}}}{l+m+2}

これを ss の積分に代入します。

011snsl+m+22l+m+2ds=1l+m+201sl+m+22nds\int_0^1 \frac{1}{s^n} \frac{s^{\frac{l+m+2}{2}}}{l+m+2} \, ds = \frac{1}{l+m+2} \int_0^1 s^{\frac{l+m+2}{2} - n} \, ds

この広義積分が有限(収束)となるためには、指数が 1-1 より大きい必要があります。

l+m+22n>1\frac{l+m+2}{2} - n > -1
l+m+2>2(n1)=2n2l+m+2 > 2(n-1) = 2n - 2
l+m+4>2n    n<l+m+42=2+l+m2l+m+4 > 2n \implies n < \frac{l+m+4}{2} = 2 + \frac{l+m}{2}

問題の条件として、 l,m,nl, m, n は正の整数であり、 l+mnl+m \leq n を満たす必要があります。 したがって、以下の連立不等式を満たす正の整数の組 (l,m,n)(l,m,n) を求めます。

l+mn<2+l+m2l+m \leq n < 2 + \frac{l+m}{2}

この不等式から、

l+m<2+l+m2l+m < 2 + \frac{l+m}{2}
l+m2<2    l+m<4\frac{l+m}{2} < 2 \implies l+m < 4

l,ml, m は正の整数なので、 l1,m1l \geq 1, m \geq 1 であり、 l+m2l+m \geq 2 です。したがって、 l+ml+m の取りうる値は 2233 です。

場合 1: l+m=2l+m=2 l,ml, m は正の整数なので、 (l,m)=(1,1)(l,m) = (1,1) のみです。 このとき、 nn の条件は、

2n<2+22    2n<32 \leq n < 2 + \frac{2}{2} \implies 2 \leq n < 3

nn は整数なので、 n=2n=2 となります。 よって、組 (1,1,2)(1,1,2) が一つの解です。

場合 2: l+m=3l+m=3 l,ml, m は正の整数なので、 (l,m)(l,m)(1,2)(1,2) または (2,1)(2,1) です。 このとき、 nn の条件は、

3n<2+32    3n<3.53 \leq n < 2 + \frac{3}{2} \implies 3 \leq n < 3.5

nn は整数なので、 n=3n=3 となります。 よって、組 (1,2,3)(1,2,3)(2,1,3)(2,1,3) が解となります。

以上をまとめると、求める正の整数の組 (l,m,n)(l, m, n) は以下の3つです。

(1,1,2),(1,2,3),(2,1,3)(1, 1, 2), \quad (1, 2, 3), \quad (2, 1, 3)

問2

(i) 帰納法による証明: n=1n = 1 のとき, x1=1x_1 = 1 であり, x2=sinx1=sin1x_2 = \sin x_1 = \sin 1 . 0<sin1<10 < \sin 1 < 1 なので, 0<x2<x10 < x_2 < x_1 が成り立つ。

n=kn = k のとき, 0<xk+1<xk0 < x_{k+1} < x_k が成り立つと仮定する。 xk+2=sinxk+1x_{k+2} = \sin x_{k+1} である。 0<xk+1<xk0 < x_{k+1} < x_k より, 0<sinxk+1<xk+10 < \sin x_{k+1} < x_{k+1} が成り立つ。 したがって, 0<xk+2<xk+10 < x_{k+2} < x_{k+1} が成り立つ。 よって, すべての nn に対して, 0<xn+1<xn0 < x_{n+1} < x_n が成り立つ。

(ii) 0<xn+1<xn0 < x_{n+1} < x_n より, {xn}\{x_n\} は単調減少数列であり, 下に有界であるから, 極限が存在する。 limnxn=L\lim_{n \to \infty} x_n = L とすると, limnxn+1=limnsinxn\lim_{n \to \infty} x_{n+1} = \lim_{n \to \infty} \sin x_n より, L=sinLL = \sin L となる。 この方程式の解は L=0L = 0 である。 したがって, limnxn=0\lim_{n \to \infty} x_n = 0 である。

(iii) sinx\sin x のマクローリン展開は次の通りである。

sinx=xx33!+o(x4)=xx36+o(x4).\sin x =x-\frac{x^3}{3!}+o(x^4) =x-\frac{x^3}{6}+o(x^4).

(iv)

sinx=x(1x26+o(x3))\sin x =x\left(1-\frac{x^2}{6}+o(x^3)\right)
1(sinx)2=1x2(1x26+o(x3))2=1x2(1+x23+o(x2))=1x2+13+o(1).\frac{1}{(\sin x)^2} =\frac{1}{x^2} \left(1-\frac{x^2}{6}+o(x^3)\right)^{-2} =\frac{1}{x^2} \left(1+\frac{x^2}{3}+o(x^2)\right) =\frac{1}{x^2}+\frac{1}{3}+o(1).

したがって, a=1a = 1 , b=13b = \frac{1}{3} である。

(v) xn+1=sinxn=xnxn36+o(xn4)x_{n+1} = \sin x_n = x_n - \frac{x_n^3}{6} + o(x_n^4) である。

1xn+12=1sin2xn=1xn2+13+o(1)\frac{1}{x_{n+1}^2} = \frac{1}{\sin^2 x_n} = \frac{1}{x_n^2} + \frac{1}{3} + o(1)
limn(1xn+121xn2)=limn(1sin2xn1xn2)=13\lim_{n \to \infty} \left(\frac{1}{x_{n+1}^2} - \frac{1}{x_n^2}\right) = \lim_{n \to \infty} \left(\frac{1}{\sin^2 x_n} - \frac{1}{x_n^2}\right) = \frac{1}{3}

(vi)

1xn+121xn2=13+o(1)\frac{1}{x_{n+1}^2} - \frac{1}{x_n^2} = \frac{1}{3} + o(1)
k=1n1(1xk+121xk2)=k=1n1(13+o(1))\sum_{k=1}^{n-1} \left(\frac{1}{x_{k+1}^2} - \frac{1}{x_k^2}\right) = \sum_{k=1}^{n-1} \left(\frac{1}{3} + o(1)\right)
1xn21x12=n13+o(n)\frac{1}{x_n^2} - \frac{1}{x_1^2} = \frac{n-1}{3} + o(n)
1xn2=n3+o(n)\frac{1}{x_n^2} = \frac{n}{3} + o(n)
nxn2=nn3+o(n)=113+o(1)nx_n^2 = \frac{n}{\frac{n}{3} + o(n)} = \frac{1}{\frac{1}{3} + o(1)}
limnnxn2=3\lim_{n \to \infty} nx_n^2 = 3