跳到主要内容

北海道大学 情報科学院 情報科学専攻 システム情報科学コース 2022年8月実施 専門科目1 問1 (応用数学I)

Author

Miyake

Description

題意の要約。出典:大学公開原卷の保存版(PDF 3ページ)

1-1)

x=(x1,x2,x3)T\boldsymbol x=(x_1,x_2,x_3)^T および二次形式

Q=2x12+5x22+2x322x1x22x2x3+4x3x1Q=2x_1^2+5x_2^2+2x_3^2-2x_1x_2-2x_2x_3+4x_3x_1

について、次を求める。

(a) Q=xTAxQ=\boldsymbol x^TA\boldsymbol x を満たす対称行列 AA

(b) AA の全固有値と、それぞれに属する単位固有ベクトル。

(c) xTx=1\boldsymbol x^T\boldsymbol x=1 の下での QQ の最大値と、その理由。

1-2)

定数 a,b,ca,b,c に対し、デカルト座標の単位ベクトルを i,j,k\boldsymbol i,\boldsymbol j,\boldsymbol k として B=axi+byj+czk\boldsymbol B=ax\boldsymbol i+by\boldsymbol j+cz\boldsymbol k を考える。

(a) divB\operatorname{div}\boldsymbol BcurlB\operatorname{curl}\boldsymbol B を求める。

(b) 原点中心・半径 rr の上半球面 SSz0z\ge0、底面を含まない)について、原点側を内側とした外向き単位法線を点 (x,y,z)(x,y,z) で表す。

(c) (b) の面 SS に対する通量 SBndA\iint_S\boldsymbol B\cdot\boldsymbol n\,dA を求める。dAdA は面積要素である。

题目描述

1-1. 二次型与特征值:令 x=(x1,x2,x3)\boldsymbol x=(x_1,x_2,x_3)^\top,考虑

Q=2x12+5x22+2x322x1x22x2x3+4x3x1.Q=2x_1^2+5x_2^2+2x_3^2-2x_1x_2-2x_2x_3+4x_3x_1.
  1. 求满足 Q=xAxQ=\boldsymbol x^\top A\boldsymbol x 的对称矩阵 AA
  2. AA 的全部特征值和相应单位特征向量。
  3. xx=1\boldsymbol x^\top\boldsymbol x=1 的条件下求 QQ 的最大值,并说明理由。

1-2. 向量分析:给定常数 a,b,ca,b,c,令 B=(ax,by,cz)\boldsymbol B=(ax,by,cz)

  1. divB\operatorname{div}\boldsymbol BcurlB\operatorname{curl}\boldsymbol B
  2. 对以原点为球心、半径为 rr 的上半球面 SSz0z\ge0,不含底面),写出点 (x,y,z)(x,y,z) 处的外向单位法向量。
  3. 求穿过该半球面的通量 SBndA\iint_S\boldsymbol B\cdot\boldsymbol n\,dA

Kai

1-1)

(a)

A=(212151212) \begin{align} A = \begin{pmatrix} 2 & -1 & 2 \\ -1 & 5 & -1 \\ 2 & -1 & 2 \end{pmatrix} \end{align}

(b)

AA の固有値を aa とすると、

0=det(2a1215a1212a)=a3+9a218a=a(a3)(a6)  a=0,3,6\begin{align} 0 &= \det \begin{pmatrix} 2-a & -1 & 2 \\ -1 & 5-a & -1 \\ 2 & -1 & 2-a \end{pmatrix} \\ &= -a^3 + 9a^2 - 18a \\ &= -a(a-3)(a-6) \\ \therefore \ \ a &= 0, 3, 6 \end{align}

である。 固有値 a=0,3,6a=0,3,6 に属する大きさ 11 の固有ベクトルは、それぞれ、

12(101),  13(111),  16(121)\begin{align} \frac{1}{\sqrt{2}} \begin{pmatrix} 1 \\ 0 \\ -1 \end{pmatrix} , \ \ \frac{1}{\sqrt{3}} \begin{pmatrix} 1 \\ 1 \\ 1 \end{pmatrix} , \ \ \frac{1}{\sqrt{6}} \begin{pmatrix} 1 \\ -2 \\ 1 \end{pmatrix} \end{align}

である。

(c)

AA の最大の固有値は 66 であるから、 xTx=1\boldsymbol{x}^T\boldsymbol{x}=1 のときの QQ の最大値は 66 である。実際、上記の正規直交固有ベクトルを u0,u3,u6\boldsymbol u_0,\boldsymbol u_3,\boldsymbol u_6 として x=λdλuλ\boldsymbol x=\sum_\lambda d_\lambda\boldsymbol u_\lambda と書けば、λdλ2=1\sum_\lambda d_\lambda^2=1 であり、Q=3d32+6d626Q=3d_3^2+6d_6^2\le6 となる。x=±u6\boldsymbol x=\pm\boldsymbol u_6 で等号が成立する。

1-2)

(a)

divB=a+b+c,curlB=0.\operatorname{div}\boldsymbol B=a+b+c,\qquad \operatorname{curl}\boldsymbol B=\boldsymbol0.

(b)

球面上では x2+y2+z2=r2x^2+y^2+z^2=r^2 なので、

n=1r(x,y,z).\boldsymbol n=\frac1r(x,y,z).

(c)

上半球内部を VV とし、赤道の円板 DD を加えて閉曲面を作る。DD 上では z=0z=0 であり、下向き法線に対する通量は cz=0-cz=0 である。発散定理より、

SBndA=V(a+b+c)dV=2πr33(a+b+c).\iint_S\boldsymbol B\cdot\boldsymbol n\,dA =\iiint_V(a+b+c)\,dV =\boxed{\frac{2\pi r^3}{3}(a+b+c)}.