京都大学 情報学研究科 通信情報システム専攻 2022年8月実施 専門基礎A [A-1]
Author
SUN , 祭音Myyura
Description
Answer all the following questions.
(1)
Find all the local maxima and minima, and corresponding x x x and y y y with respect to function f ( x , y ) f(x,y) f ( x , y ) . Let x x x and y y y be real numbers.
f ( x , y ) = x 3 − x 2 y + x y 2 − x f(x,y)=x^3-x^2y+xy^2-x f ( x , y ) = x 3 − x 2 y + x y 2 − x
(2)
Let D D D be a domain bounded by y = x y=\sqrt{x} y = x and y = x y=x y = x , where x ≥ 0 x\ge 0 x ≥ 0 . Compute the following integral I I I .
I = ∬ D e − y d x d y I=\iint_D e^{-y}\,dx\,dy I = ∬ D e − y d x d y
(3)
Find the length of the curve given as follows.
y = x 3 / 2 ( 0 ≤ x ≤ 4 3 ) y=x^{3/2}\quad \left(0\le x\le \frac{4}{3}\right) y = x 3/2 ( 0 ≤ x ≤ 3 4 )
(4)
Find the eigenvectors of matrix A A A , and show the conditions on which they become orthogonal to each other. Let x x x be a real number.
A = ( x a b x ) A=\begin{pmatrix}
x & a\\
b & x
\end{pmatrix} A = ( x b a x )
题目描述
回答全部问题。
对实数 x , y x,y x , y ,求函数
f ( x , y ) = x 3 − x 2 y + x y 2 − x f(x,y)=x^3-x^2y+xy^2-x f ( x , y ) = x 3 − x 2 y + x y 2 − x
的所有局部极大值、局部极小值及其对应的 ( x , y ) (x,y) ( x , y ) 。
设 D D D 为 x ≥ 0 x\ge0 x ≥ 0 时曲线 y = x y=\sqrt x y = x 与 y = x y=x y = x 围成的区域,计算
I = ∬ D e − y d x d y . I=\iint_D e^{-y}\,dx\,dy. I = ∬ D e − y d x d y .
求曲线
y = x 3 / 2 , 0 ≤ x ≤ 4 3 y=x^{3/2},\qquad 0\le x\le\frac43 y = x 3/2 , 0 ≤ x ≤ 3 4
的弧长。
对矩阵
A = ( x a b x ) , A=\begin{pmatrix}x&a\\b&x\end{pmatrix}, A = ( x b a x ) ,
求其特征向量,并给出这些特征向量彼此正交的条件,其中 x x x 为实数。
Kai
(1)
To find critical points, set partial derivatives to zero:
{ f x = 3 x 2 − 2 y x + y 2 − 1 = 0 f y = − x 2 + 2 x y = x ( 2 y − x ) = 0 \begin{cases}
f_x = 3x^2 - 2yx + y^2 - 1 = 0 \\
f_y = -x^2 + 2xy = x(2y - x) = 0
\end{cases} { f x = 3 x 2 − 2 y x + y 2 − 1 = 0 f y = − x 2 + 2 x y = x ( 2 y − x ) = 0
From f y = 0 f_y = 0 f y = 0 , we have x = 0 x = 0 x = 0 or x = 2 y x = 2y x = 2 y .
If x = 0 x = 0 x = 0 , then y 2 − 1 = 0 ⇒ y = ± 1 y^2 - 1 = 0 \Rightarrow y = \pm 1 y 2 − 1 = 0 ⇒ y = ± 1 . Points: ( 0 , 1 ) , ( 0 , − 1 ) (0, 1), (0, -1) ( 0 , 1 ) , ( 0 , − 1 ) .
If x = 2 y x = 2y x = 2 y , then 3 ( 2 y ) 2 − 2 y ( 2 y ) + y 2 − 1 = 0 ⇒ 9 y 2 = 1 ⇒ y = ± 1 3 3(2y)^2 - 2y(2y) + y^2 - 1 = 0 \Rightarrow 9y^2 = 1 \Rightarrow y = \pm \frac{1}{3} 3 ( 2 y ) 2 − 2 y ( 2 y ) + y 2 − 1 = 0 ⇒ 9 y 2 = 1 ⇒ y = ± 3 1 . Points: ( 2 3 , 1 3 ) , ( − 2 3 , − 1 3 ) (\frac{2}{3}, \frac{1}{3}), (-\frac{2}{3}, -\frac{1}{3}) ( 3 2 , 3 1 ) , ( − 3 2 , − 3 1 ) .
Second-order partial derivatives:
f x x = 6 x − 2 y , f x y = − 2 x + 2 y , f y y = 2 x f_{xx} = 6x - 2y, \quad f_{xy} = -2x + 2y, \quad f_{yy} = 2x f xx = 6 x − 2 y , f x y = − 2 x + 2 y , f yy = 2 x
Hessian determinant H = f x x f y y − ( f x y ) 2 H = f_{xx}f_{yy} - (f_{xy})^2 H = f xx f yy − ( f x y ) 2 :
For ( 0 , 1 ) : H = ( − 2 ) ( 0 ) − ( 2 ) 2 = − 4 < 0 (0, 1): H = ( -2)(0) - (2)^2 = -4 < 0 ( 0 , 1 ) : H = ( − 2 ) ( 0 ) − ( 2 ) 2 = − 4 < 0 (Saddle point)
For ( 0 , − 1 ) : H = ( 2 ) ( 0 ) − ( − 2 ) 2 = − 4 < 0 (0, -1): H = (2)(0) - (-2)^2 = -4 < 0 ( 0 , − 1 ) : H = ( 2 ) ( 0 ) − ( − 2 ) 2 = − 4 < 0 (Saddle point)
For ( 2 3 , 1 3 ) : H = ( 10 3 ) ( 4 3 ) − ( − 2 3 ) 2 = 4 > 0 (\frac{2}{3}, \frac{1}{3}): H = (\frac{10}{3})(\frac{4}{3}) - (-\frac{2}{3})^2 = 4 > 0 ( 3 2 , 3 1 ) : H = ( 3 10 ) ( 3 4 ) − ( − 3 2 ) 2 = 4 > 0 . Since f x x = 10 3 > 0 f_{xx} = \frac{10}{3} > 0 f xx = 3 10 > 0 , it is a local minimum:
f ( 2 3 , 1 3 ) = ( 2 3 ) 3 − ( 2 3 ) 2 ( 1 3 ) + ( 2 3 ) ( 1 3 ) 2 − 2 3 = − 4 9 f(\frac{2}{3}, \frac{1}{3}) = (\frac{2}{3})^3 - (\frac{2}{3})^2(\frac{1}{3}) + (\frac{2}{3})(\frac{1}{3})^2 - \frac{2}{3} = -\frac{4}{9} f ( 3 2 , 3 1 ) = ( 3 2 ) 3 − ( 3 2 ) 2 ( 3 1 ) + ( 3 2 ) ( 3 1 ) 2 − 3 2 = − 9 4
For ( − 2 3 , − 1 3 ) : H = ( − 10 3 ) ( − 4 3 ) − ( 2 3 ) 2 = 4 > 0 (-\frac{2}{3}, -\frac{1}{3}): H = (-\frac{10}{3})(-\frac{4}{3}) - (\frac{2}{3})^2 = 4 > 0 ( − 3 2 , − 3 1 ) : H = ( − 3 10 ) ( − 3 4 ) − ( 3 2 ) 2 = 4 > 0 . Since f x x = − 10 3 < 0 f_{xx} = -\frac{10}{3} < 0 f xx = − 3 10 < 0 , it is a local maximum:
f ( − 2 3 , − 1 3 ) = ( − 2 3 ) 3 − ( − 2 3 ) 2 ( − 1 3 ) + ( − 2 3 ) ( − 1 3 ) 2 − ( − 2 3 ) = 4 9 \begin{aligned}
f(-\frac{2}{3}, -\frac{1}{3}) &= (-\frac{2}{3})^3 - (-\frac{2}{3})^2(-\frac{1}{3}) \\
&+ (-\frac{2}{3})(-\frac{1}{3})^2 - (-\frac{2}{3}) \\
&= \frac{4}{9}
\end{aligned} f ( − 3 2 , − 3 1 ) = ( − 3 2 ) 3 − ( − 3 2 ) 2 ( − 3 1 ) + ( − 3 2 ) ( − 3 1 ) 2 − ( − 3 2 ) = 9 4
(2)
The domain D D D is bounded by y = x y = \sqrt{x} y = x and y = x y = x y = x , which intersect at ( 0 , 0 ) (0, 0) ( 0 , 0 ) and ( 1 , 1 ) (1, 1) ( 1 , 1 ) .
For 0 ≤ y ≤ 1 0 \le y \le 1 0 ≤ y ≤ 1 , the range of x x x is y 2 ≤ x ≤ y y^2 \le x \le y y 2 ≤ x ≤ y .
I = ∬ D e − y d A = ∫ 0 1 ∫ y 2 y e − y d x d y = ∫ 0 1 ( y − y 2 ) e − y d y I = \iint_D e^{-y} dA = \int_0^1 \int_{y^2}^y e^{-y} dx dy = \int_0^1 (y - y^2) e^{-y} dy I = ∬ D e − y d A = ∫ 0 1 ∫ y 2 y e − y d x d y = ∫ 0 1 ( y − y 2 ) e − y d y
Using integration by parts:
∫ ( y − y 2 ) e − y d y = − ( y − y 2 ) e − y + ∫ ( 1 − 2 y ) e − y d y = ( y 2 − y ) e − y − ( 1 − 2 y ) e − y + ∫ ( − 2 ) e − y d y = ( y 2 − y − 1 + 2 y + 2 ) e − y = ( y 2 + y + 1 ) e − y \begin{aligned}
\int (y - y^2) e^{-y} dy &= -(y - y^2)e^{-y} + \int (1 - 2y) e^{-y} dy\\
&= (y^2 - y)e^{-y} - (1 - 2y)e^{-y} + \int (-2) e^{-y} dy\\
&= (y^2 - y - 1 + 2y + 2)e^{-y} \\
&= (y^2 + y + 1)e^{-y}
\end{aligned} ∫ ( y − y 2 ) e − y d y = − ( y − y 2 ) e − y + ∫ ( 1 − 2 y ) e − y d y = ( y 2 − y ) e − y − ( 1 − 2 y ) e − y + ∫ ( − 2 ) e − y d y = ( y 2 − y − 1 + 2 y + 2 ) e − y = ( y 2 + y + 1 ) e − y
Evaluating from 0 to 1:
I = [ ( y 2 + y + 1 ) e − y ] 0 1 = 3 e − 1 − 1 I = \left[ (y^2 + y + 1)e^{-y} \right]_0^1 = 3e^{-1} - 1 I = [ ( y 2 + y + 1 ) e − y ] 0 1 = 3 e − 1 − 1
(3)
The arc length L L L for y = x 3 / 2 y = x^{3/2} y = x 3/2 from x = 0 x=0 x = 0 to x = 4 / 3 x=4/3 x = 4/3 :
y ′ = 3 2 x 1 / 2 ⇒ 1 + ( y ′ ) 2 = 1 + 9 4 x y' = \frac{3}{2}x^{1/2} \Rightarrow 1 + (y')^2 = 1 + \frac{9}{4}x y ′ = 2 3 x 1/2 ⇒ 1 + ( y ′ ) 2 = 1 + 4 9 x
L = ∫ 0 4 / 3 1 + 9 4 x d x L = \int_0^{4/3} \sqrt{1 + \frac{9}{4}x} dx L = ∫ 0 4/3 1 + 4 9 x d x
Let u = 1 + 9 4 x u = 1 + \frac{9}{4}x u = 1 + 4 9 x , then d u = 9 4 d x du = \frac{9}{4}dx d u = 4 9 d x .
When x = 0 , u = 1 x=0, u=1 x = 0 , u = 1 ; when x = 4 / 3 , u = 4 x=4/3, u=4 x = 4/3 , u = 4 .
L = ∫ 1 4 u ⋅ 4 9 d u = 4 9 [ 2 3 u 3 / 2 ] 1 4 = 8 27 ( 4 3 / 2 − 1 3 / 2 ) = 8 27 ( 8 − 1 ) = 56 27 \begin{aligned}
L &= \int_1^4 \sqrt{u} \cdot \frac{4}{9} du = \frac{4}{9} \left[ \frac{2}{3}u^{3/2} \right]_1^4 \\
&= \frac{8}{27}(4^{3/2} - 1^{3/2}) = \frac{8}{27}(8 - 1) \\
&= \frac{56}{27}
\end{aligned} L = ∫ 1 4 u ⋅ 9 4 d u = 9 4 [ 3 2 u 3/2 ] 1 4 = 27 8 ( 4 3/2 − 1 3/2 ) = 27 8 ( 8 − 1 ) = 27 56
(4)
The characteristic equation is
det ( A − λ I ) = ∣ x − λ a b x − λ ∣ = ( x − λ ) 2 − a b = 0. \det(A-\lambda I)
=\begin{vmatrix}x-\lambda&a\\b&x-\lambda\end{vmatrix}
=(x-\lambda)^2-ab=0. det ( A − λ I ) = x − λ b a x − λ = ( x − λ ) 2 − ab = 0.
If a b ≠ 0 ab\ne0 ab = 0 , put s = a b s=\sqrt{ab} s = ab (over C \mathbb C C ). Then
λ ± = x ± s . \lambda_\pm=x\pm s. λ ± = x ± s .
For λ + = x + s \lambda_+=x+s λ + = x + s ,
( − s a b − s ) ( v 1 v 2 ) = 0 , − s v 1 + a v 2 = 0 , \begin{pmatrix}-s&a\\b&-s\end{pmatrix}
\begin{pmatrix}v_1\\v_2\end{pmatrix}=0,
\qquad -sv_1+av_2=0, ( − s b a − s ) ( v 1 v 2 ) = 0 , − s v 1 + a v 2 = 0 ,
and for λ − = x − s \lambda_-=x-s λ − = x − s ,
( s a b s ) ( v 1 v 2 ) = 0 , s v 1 + a v 2 = 0. \begin{pmatrix}s&a\\b&s\end{pmatrix}
\begin{pmatrix}v_1\\v_2\end{pmatrix}=0,
\qquad sv_1+av_2=0. ( s b a s ) ( v 1 v 2 ) = 0 , s v 1 + a v 2 = 0.
Therefore
E λ + = span { [ a s ] } , E λ − = span { [ a − s ] } . E_{\lambda_+}=\operatorname{span}\left\{\begin{bmatrix}a\\s\end{bmatrix}\right\},
\qquad
E_{\lambda_-}=\operatorname{span}\left\{\begin{bmatrix}a\\-s\end{bmatrix}\right\}. E λ + = span { [ a s ] } , E λ − = span { [ a − s ] } .
If a = b = 0 a=b=0 a = b = 0 , the eigenspace for λ = x \lambda=x λ = x is R 2 \mathbb R^2 R 2 . If a b = 0 ab=0 ab = 0 but ( a , b ) ≠ ( 0 , 0 ) (a,b)\ne(0,0) ( a , b ) = ( 0 , 0 ) , the only eigenvalue is x x x and its eigenspace is span { ( 1 , 0 ) T } \operatorname{span}\{(1,0)^T\} span {( 1 , 0 ) T } for a ≠ 0 a\ne0 a = 0 , or span { ( 0 , 1 ) T } \operatorname{span}\{(0,1)^T\} span {( 0 , 1 ) T } for b ≠ 0 b\ne0 b = 0 .
For two real eigenvectors we need a b > 0 ab>0 ab > 0 . Their inner product is
[ a s ] T [ a − s ] = a 2 − a b = a ( a − b ) , \begin{bmatrix}a\\s\end{bmatrix}^{\!T}
\begin{bmatrix}a\\-s\end{bmatrix}
=a^2-ab=a(a-b), [ a s ] T [ a − s ] = a 2 − ab = a ( a − b ) ,
so they are orthogonal exactly when a = b ≠ 0 a=b\ne0 a = b = 0 . When a = b = 0 a=b=0 a = b = 0 , any orthogonal basis is an eigenbasis. Hence A A A has an orthogonal real eigenbasis if and only if a = b a=b a = b .