跳到主要内容

京都大学 情報学研究科 システム科学専攻 2020年8月実施 微積分

Author

思齐塾, 祭音Myyura

Description

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

题目描述

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 的三元组。

Kai

(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)