跳到主要内容

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

Author

思齐塾, 祭音Myyura

Description

(1) f(x,y)=xy(x2+y21)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+z21x^2 + y^2 + z^2 \leq 1 とし, g(x,y,z)=x2y2z2g(x,y,z) = x^2y^2z^2 とおく。

(a) ggBB における最大値を求めよ。

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

题目描述

(1) 设

f(x,y)=xy(x2+y21).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+z21,x^2+y^2+z^2\leq 1,

并令

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

(a) 求 ggBB 上的最大值。

(b) 计算

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

Kai

(1) (a)

f(x,y)=xy(x2+y21)=x3y+xy3xyf(x,y) = xy(x^2+y^2-1) = x^3y + xy^3 - xy
fx=3x2y+y3y\frac{\partial f}{\partial x} = 3x^2y + y^3 - y
fy=x3+3xy2x\frac{\partial f}{\partial y} = x^3 + 3xy^2 - x

ヤコビ行列は

Jf=(fxfy)=(3x2y+y3yx3+3xy2x)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}

ヘッセ行列は

2fx2=6xy\frac{\partial^2 f}{\partial x^2} = 6xy
2fxy=3x2+3y21\frac{\partial^2 f}{\partial x \partial y} = 3x^2 + 3y^2 - 1
2fy2=6xy\frac{\partial^2 f}{\partial y^2} = 6xy
Hf=(2fx22fxy2fyx2fy2)=(6xy3x2+3y213x2+3y216xy)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) fx=3x2y+y3y=y(3x2+y21)=0\frac{\partial f}{\partial x} = 3x^2y + y^3 - y = y(3x^2 + y^2 - 1) = 0 fy=x3+3xy2x=x(x2+3y21)=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(y21)=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(x21)=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: x0x \neq 0 and y0y \neq 0 . Then 3x2+y2=13x^2 + y^2 = 1 and x2+3y2=1x^2 + 3y^2 = 1 . Subtracting gives 2x22y2=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+141)=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+141)=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+141)=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+141)=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:

pointsfdetHfclassification(0,0), (0,±1), (±1,0)01 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+y21)z = f(x,y) = xy(x^2+y^2-1) 接平面の方程式は zf(1,1)=fx(1,1)(x1)+fy(1,1)(y1)z - f(1,1) = f_x(1,1)(x-1) + f_y(1,1)(y-1) f(1,1)=1(1+11)=1f(1,1) = 1(1+1-1) = 1 fx(1,1)=3(1)2(1)+(1)31=3+11=3f_x(1,1) = 3(1)^2(1) + (1)^3 - 1 = 3+1-1 = 3 fy(1,1)=(1)3+3(1)(1)21=1+31=3f_y(1,1) = (1)^3 + 3(1)(1)^2 - 1 = 1+3-1 = 3 Therefore, z1=3(x1)+3(y1)z - 1 = 3(x-1) + 3(y-1) z1=3x3+3y3z - 1 = 3x - 3 + 3y - 3 z=3x+3y5z = 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 0r10 \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=r6sin4θcos2θcos2ϕsin2ϕ=r6sin4θcos2θ(14sin22ϕ)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 0r10 \le r \le 1 , 0θπ0 \le \theta \le \pi , 0ϕ2π0 \le \phi \le 2\pi . Since r1r \le 1 , r6r^6 is maximized when r=1r = 1 . sin4θcos2θ\sin^4 \theta \cos^2 \theta is maximized when θ\theta is between 0 and π\pi . ddθ(sin4θcos2θ)=4sin3θcos3θ2sin5θcosθ=2sin3θcosθ(2cos2θsin2θ)=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 2cos2θ=sin2θ2 \cos^2 \theta = \sin^2 \theta . Since 2cos2θ=sin2θ2 \cos^2 \theta = \sin^2 \theta , 2cos2θ=1cos2θ2 \cos^2 \theta = 1 - \cos^2 \theta , so 3cos2θ=13 \cos^2 \theta = 1 , and cosθ=±13\cos \theta = \pm \frac{1}{\sqrt{3}} . Then sinθ=±23\sin \theta = \pm \sqrt{\frac{2}{3}} . Thus, the maximum occurs when θ=arccos(13)\theta = \arccos (\frac{1}{\sqrt{3}}) and sinθ=23\sin \theta = \sqrt{\frac{2}{3}} . Thus, the maximum value of sin4θcos2θ\sin^4 \theta \cos^2 \theta is (23)213=427(\frac{2}{3})^2 \frac{1}{3} = \frac{4}{27} . 14sin22ϕ\frac{1}{4} \sin^2 2\phi is maximized when sin2ϕ=±1\sin 2\phi = \pm 1 , which means 2ϕ=π2,3π22\phi = \frac{\pi}{2}, \frac{3\pi}{2} . So ϕ=π4,3π4\phi = \frac{\pi}{4}, \frac{3\pi}{4} . The maximum value of 14sin22ϕ\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 1642714=1271^6 \cdot \frac{4}{27} \cdot \frac{1}{4} = \frac{1}{27} .

(b) Bg(x,y,z)dxdydz=Bx2y2z2dxdydz=010π02π(rsinθcosϕ)2(rsinθsinϕ)2(rcosθ)2r2sinθdϕdθdr=01r8dr0πsin5θcos2θdθ02πcos2ϕsin2ϕdϕ=190πsin5θcos2θdθ02π14sin22ϕdϕ=190πsin5θcos2θdθ1402π1cos4ϕ2dϕ=190πsin5θcos2θdθ14122π=π360πsin5θcos2θ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πsin5θcos2θdθ=20π/2sin5θcos2θ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,

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