跳到主要内容

明治大学 先端数理科学研究科 現象数理学専攻 2019年8月実施 微积分

Author​

思齐塾, 祭音Myyura

Description​

(1) f(x,y)=xy(x2+y2−1)f(x,y) = xy(x^2 + y^2 - 1) とする。

(a) ff のヤコビ行列とヘッセ行列を求めよ。

(b) ff の極値をすべて求めよ。

(c) ff のグラフ z=f(x,y)z = f(x,y) の (1,1,1)(1,1,1) における接平面の方程式を求めよ。

(2) BB を単位球 x2+y2+z2≤1x^2 + y^2 + z^2 \leq 1 とし, g(x,y,z)=x2y2z2g(x,y,z) = x^2y^2z^2 とおく。

(a) gg の BB における最大値を求めよ。

(b) ∭Bg(x,y,z)dxdydz\iiint_B g(x,y,z) dx dy dz の値を求めよ。

题目描述​

(1) 设

f(x,y)=xy(x2+y2−1).f(x,y)=xy(x^2+y^2-1).

(a) 求 ff 的 Jacobian 矩阵与 Hessian 矩阵。

(b) 求 ff 的全部极值。

(c) 求曲面 z=f(x,y)z=f(x,y) 在点 (1,1,1)(1,1,1) 处的切平面方程。

(2) 设 BB 为单位球体

x2+y2+z2≤1,x^2+y^2+z^2\leq 1,

并令

g(x,y,z)=x2y2z2.g(x,y,z)=x^2y^2z^2.

(a) 求 gg 在 BB 上的最大值。

(b) 计算

∭Bg(x,y,z) dx dy dz.\iiint_B g(x,y,z)\,dx\,dy\,dz.

Kai​

(1) (a)

f(x,y)=xy(x2+y2−1)=x3y+xy3−xyf(x,y) = xy(x^2+y^2-1) = x^3y + xy^3 - xy
∂f∂x=3x2y+y3−y\frac{\partial f}{\partial x} = 3x^2y + y^3 - y
∂f∂y=x3+3xy2−x\frac{\partial f}{\partial y} = x^3 + 3xy^2 - x

ヤコビ行列は

Jf=(∂f∂x∂f∂y)=(3x2y+y3−yx3+3xy2−x)J_f = \begin{pmatrix} \frac{\partial f}{\partial x} & \frac{\partial f}{\partial y} \end{pmatrix} = \begin{pmatrix} 3x^2y + y^3 - y & x^3 + 3xy^2 - x \end{pmatrix}

ヘッセ行列は

∂2f∂x2=6xy\frac{\partial^2 f}{\partial x^2} = 6xy
∂2f∂x∂y=3x2+3y2−1\frac{\partial^2 f}{\partial x \partial y} = 3x^2 + 3y^2 - 1
∂2f∂y2=6xy\frac{\partial^2 f}{\partial y^2} = 6xy
Hf=(∂2f∂x2∂2f∂x∂y∂2f∂y∂x∂2f∂y2)=(6xy3x2+3y2−13x2+3y2−16xy)H_f = \begin{pmatrix} \frac{\partial^2 f}{\partial x^2} & \frac{\partial^2 f}{\partial x \partial y} \\ \frac{\partial^2 f}{\partial y \partial x} & \frac{\partial^2 f}{\partial y^2} \end{pmatrix} = \begin{pmatrix} 6xy & 3x^2 + 3y^2 - 1 \\ 3x^2 + 3y^2 - 1 & 6xy \end{pmatrix}

(b) ∂f∂x=3x2y+y3−y=y(3x2+y2−1)=0\frac{\partial f}{\partial x} = 3x^2y + y^3 - y = y(3x^2 + y^2 - 1) = 0 ∂f∂y=x3+3xy2−x=x(x2+3y2−1)=0\frac{\partial f}{\partial y} = x^3 + 3xy^2 - x = x(x^2 + 3y^2 - 1) = 0

Case 1: x=0x=0 . Then y(y2−1)=0y(y^2 - 1) = 0 , so y=0,1,−1y = 0, 1, -1 . Critical points: (0,0),(0,1),(0,−1)(0,0), (0,1), (0,-1) . Case 2: y=0y=0 . Then x(x2−1)=0x(x^2 - 1) = 0 , so x=0,1,−1x = 0, 1, -1 . Critical points: (0,0),(1,0),(−1,0)(0,0), (1,0), (-1,0) . Case 3: x≠0x \neq 0 and y≠0y \neq 0 . Then 3x2+y2=13x^2 + y^2 = 1 and x2+3y2=1x^2 + 3y^2 = 1 . Subtracting gives 2x2−2y2=02x^2 - 2y^2 = 0 , so x2=y2x^2 = y^2 , and x=±yx = \pm y . Substituting x=yx = y into 3x2+y2=13x^2 + y^2 = 1 gives 4x2=14x^2 = 1 , so x=±12x = \pm \frac{1}{2} , and y=±12y = \pm \frac{1}{2} . Critical points: (12,12),(12,−12),(−12,12),(−12,−12)(\frac{1}{2}, \frac{1}{2}), (\frac{1}{2}, -\frac{1}{2}), (-\frac{1}{2}, \frac{1}{2}), (-\frac{1}{2}, -\frac{1}{2}) .

Critical points: (0,0),(0,1),(0,−1),(1,0),(−1,0),(12,12),(12,−12),(−12,12),(−12,−12)(0,0), (0,1), (0,-1), (1,0), (-1,0), (\frac{1}{2}, \frac{1}{2}), (\frac{1}{2}, -\frac{1}{2}), (-\frac{1}{2}, \frac{1}{2}), (-\frac{1}{2}, -\frac{1}{2}) . Evaluate ff at these points: f(0,0)=0f(0,0) = 0 , f(0,1)=0f(0,1) = 0 , f(0,−1)=0f(0,-1) = 0 , f(1,0)=0f(1,0) = 0 , f(−1,0)=0f(-1,0) = 0 . f(12,12)=14(14+14−1)=14(−12)=−18f(\frac{1}{2}, \frac{1}{2}) = \frac{1}{4}(\frac{1}{4} + \frac{1}{4} - 1) = \frac{1}{4}(-\frac{1}{2}) = -\frac{1}{8} . f(12,−12)=−14(14+14−1)=18f(\frac{1}{2}, -\frac{1}{2}) = -\frac{1}{4}(\frac{1}{4} + \frac{1}{4} - 1) = \frac{1}{8} . f(−12,12)=−14(14+14−1)=18f(-\frac{1}{2}, \frac{1}{2}) = -\frac{1}{4}(\frac{1}{4} + \frac{1}{4} - 1) = \frac{1}{8} . f(−12,−12)=14(14+14−1)=−18f(-\frac{1}{2}, -\frac{1}{2}) = \frac{1}{4}(\frac{1}{4} + \frac{1}{4} - 1) = -\frac{1}{8} .

The Hessian classification is as follows:

pointsfdet⁡Hfclassification(0,0), (0,±1), (±1,0)0−1 or −4saddle(12,12), (−12,−12)−182local minimum(12,−12), (−12,12)182local maximum\begin{array}{c|c|c|c} \text{points} & f & \det H_f & \text{classification}\\ \hline (0,0),\ (0,\pm1),\ (\pm1,0) & 0 & -1\ \text{or}\ -4 & \text{saddle}\\ (\frac12,\frac12),\ (-\frac12,-\frac12) & -\frac18 & 2 & \text{local minimum}\\ (\frac12,-\frac12),\ (-\frac12,\frac12) & \frac18 & 2 & \text{local maximum} \end{array}

Indeed, at the two minimum points the Hessian eigenvalues are 1,21,2 , while at the two maximum points they are −1,−2-1,-2 . Thus these four points give all the extrema.

(c) z=f(x,y)=xy(x2+y2−1)z = f(x,y) = xy(x^2+y^2-1) 接平面の方程式は z−f(1,1)=fx(1,1)(x−1)+fy(1,1)(y−1)z - f(1,1) = f_x(1,1)(x-1) + f_y(1,1)(y-1) f(1,1)=1(1+1−1)=1f(1,1) = 1(1+1-1) = 1 fx(1,1)=3(1)2(1)+(1)3−1=3+1−1=3f_x(1,1) = 3(1)^2(1) + (1)^3 - 1 = 3+1-1 = 3 fy(1,1)=(1)3+3(1)(1)2−1=1+3−1=3f_y(1,1) = (1)^3 + 3(1)(1)^2 - 1 = 1+3-1 = 3 Therefore, z−1=3(x−1)+3(y−1)z - 1 = 3(x-1) + 3(y-1) z−1=3x−3+3y−3z - 1 = 3x - 3 + 3y - 3 z=3x+3y−5z = 3x + 3y - 5

(2) (a) Using spherical coordinates x=rsin⁡θcos⁡ϕx = r \sin \theta \cos \phi , y=rsin⁡θsin⁡ϕy = r \sin \theta \sin \phi , z=rcos⁡θz = r \cos \theta , where 0≤r≤10 \le r \le 1 , 0≤θ≤π0 \le \theta \le \pi , 0≤ϕ≤2π0 \le \phi \le 2\pi .

g(x,y,z)=x2y2z2=(rsin⁡θcos⁡ϕ)2(rsin⁡θsin⁡ϕ)2(rcos⁡θ)2=r6sin⁡4θcos⁡2θcos⁡2ϕsin⁡2ϕ=r6sin⁡4θcos⁡2θ(14sin⁡22ϕ)g(x,y,z) = x^2y^2z^2 = (r \sin \theta \cos \phi)^2(r \sin \theta \sin \phi)^2(r \cos \theta)^2 = r^6 \sin^4 \theta \cos^2 \theta \cos^2 \phi \sin^2 \phi = r^6 \sin^4 \theta \cos^2 \theta (\frac{1}{4} \sin^2 2\phi)

We want to maximize this expression subject to 0≤r≤10 \le r \le 1 , 0≤θ≤π0 \le \theta \le \pi , 0≤ϕ≤2π0 \le \phi \le 2\pi . Since r≤1r \le 1 , r6r^6 is maximized when r=1r = 1 . sin⁡4θcos⁡2θ\sin^4 \theta \cos^2 \theta is maximized when θ\theta is between 0 and π\pi . ddθ(sin⁡4θcos⁡2θ)=4sin⁡3θcos⁡3θ−2sin⁡5θcos⁡θ=2sin⁡3θcos⁡θ(2cos⁡2θ−sin⁡2θ)=0\frac{d}{d\theta} (\sin^4 \theta \cos^2 \theta) = 4\sin^3 \theta \cos^3 \theta - 2\sin^5 \theta \cos \theta = 2 \sin^3 \theta \cos \theta (2 \cos^2 \theta - \sin^2 \theta) = 0 . So either sin⁡θ=0\sin \theta = 0 , cos⁡θ=0\cos \theta = 0 , or 2cos⁡2θ=sin⁡2θ2 \cos^2 \theta = \sin^2 \theta . Since 2cos⁡2θ=sin⁡2θ2 \cos^2 \theta = \sin^2 \theta , 2cos⁡2θ=1−cos⁡2θ2 \cos^2 \theta = 1 - \cos^2 \theta , so 3cos⁡2θ=13 \cos^2 \theta = 1 , and cos⁡θ=±13\cos \theta = \pm \frac{1}{\sqrt{3}} . Then sin⁡θ=2/3\sin\theta=\sqrt{2/3} because 0≤θ≤π0\leq\theta\leq\pi. The two maximizing polar angles are θ=arccos⁡(1/3)\theta=\arccos(1/\sqrt3) and θ=π−arccos⁡(1/3)\theta=\pi-\arccos(1/\sqrt3). Thus, the maximum value of sin⁡4θcos⁡2θ\sin^4 \theta \cos^2 \theta is (23)213=427(\frac{2}{3})^2 \frac{1}{3} = \frac{4}{27} . 14sin⁡22ϕ\frac{1}{4} \sin^2 2\phi is maximized when sin⁡2ϕ=±1\sin 2\phi = \pm 1 , which gives ϕ=π/4, 3π/4, 5π/4, 7π/4\phi=\pi/4,\ 3\pi/4,\ 5\pi/4,\ 7\pi/4 on [0,2π][0,2\pi]. The maximum value of 14sin⁡22ϕ\frac{1}{4} \sin^2 2\phi is 14\frac{1}{4} . Therefore, the maximum value of g(x,y,z)g(x,y,z) is 16⋅427⋅14=1271^6 \cdot \frac{4}{27} \cdot \frac{1}{4} = \frac{1}{27} .

(b) ∭Bg(x,y,z)dxdydz=∭Bx2y2z2dxdydz=∫01∫0π∫02π(rsin⁡θcos⁡ϕ)2(rsin⁡θsin⁡ϕ)2(rcos⁡θ)2r2sin⁡θdϕdθdr=∫01r8dr∫0πsin⁡5θcos⁡2θdθ∫02πcos⁡2ϕsin⁡2ϕdϕ=19∫0πsin⁡5θcos⁡2θdθ∫02π14sin⁡22ϕdϕ=19∫0πsin⁡5θcos⁡2θdθ14∫02π1−cos⁡4ϕ2dϕ=19∫0πsin⁡5θcos⁡2θdθ14122π=π36∫0πsin⁡5θcos⁡2θdθ\iiint_B g(x,y,z) dx dy dz = \iiint_B x^2y^2z^2 dx dy dz = \int_0^1 \int_0^{\pi} \int_0^{2\pi} (r \sin \theta \cos \phi)^2(r \sin \theta \sin \phi)^2(r \cos \theta)^2 r^2 \sin \theta d\phi d\theta dr = \int_0^1 r^8 dr \int_0^{\pi} \sin^5 \theta \cos^2 \theta d\theta \int_0^{2\pi} \cos^2 \phi \sin^2 \phi d\phi = \frac{1}{9} \int_0^{\pi} \sin^5 \theta \cos^2 \theta d\theta \int_0^{2\pi} \frac{1}{4} \sin^2 2\phi d\phi = \frac{1}{9} \int_0^{\pi} \sin^5 \theta \cos^2 \theta d\theta \frac{1}{4} \int_0^{2\pi} \frac{1 - \cos 4\phi}{2} d\phi = \frac{1}{9} \int_0^{\pi} \sin^5 \theta \cos^2 \theta d\theta \frac{1}{4} \frac{1}{2} 2\pi = \frac{\pi}{36} \int_0^{\pi} \sin^5 \theta \cos^2 \theta d\theta Using the beta integral,

∫0πsin⁡5θcos⁡2θ dθ=2∫0π/2sin⁡5θcos⁡2θ dθ=B(3,32)=Γ(3)Γ(3/2)Γ(9/2)=16105.\begin{aligned} \int_0^{\pi}\sin^5\theta\cos^2\theta\,d\theta &=2\int_0^{\pi/2}\sin^5\theta\cos^2\theta\,d\theta\\ &=B\left(3,\frac32\right) =\frac{\Gamma(3)\Gamma(3/2)}{\Gamma(9/2)} =\frac{16}{105}. \end{aligned}

Therefore,

∭Bx2y2z2 dx dy dz=π3616105=4π945.\boxed{ \iiint_Bx^2y^2z^2\,dx\,dy\,dz =\frac{\pi}{36}\frac{16}{105} =\frac{4\pi}{945} }.