跳到主要内容

埼玉大学 理工学研究科 数理電子情報系専攻 情報システム工学コース 2018年2月実施 线性代数

Author

思齐塾, 祭音Myyura

Description

  1. 次の2次形式について考える。[ Consider the following quadratic form. ]
f(x,y,z)=3x2+3y2+3z2+2xy+2yz+2zxf(x, y, z) = 3x^2 + 3y^2 + 3z^2 + 2xy + 2yz + 2zx

以下の問に答えよ. [Solve the following problems. ]

(a) ベクトル x=(x,y,z)Tx = (x, y, z)^T と対称行列 AA を用いて、2次形式 f(x,y,z)=xTAxf(x,y,z) = x^T Ax と表すとき、 AA を求めよ. ここで上付添字 T^T は転置を表す. [ Find a symmetric matrix AA when the quadratic form is described as f(x,y,z)=xTAxf(x, y, z) = x^T Ax using a vector x=(x,y,z)Tx = (x, y, z)^T and the symmetric matrix AA . The superscript T^T denotes transposition. ]

(b) 前問 (a) で求めた対称行列 AA の固有値と固有ベクトルを求めよ. [ Find the eigenvalues and eigenvectors of the symmetric matrix AA obtained in the problem (a). ]

(c) 直交行列を用いて,前問 (a)で求めた対称行列 AA を対角化せよ. [Diagonalize the symmetric matrix AA obtained in the problem (a) using an orthogonal matrix. ]

(d) 変数 x,y,zx, y, zx2+y2+z2=1x^2 + y^2 + z^2 = 1 を満たすとき、2次形式 f(x,y,z)f(x,y,z) の最大値と最小値を求めよ. [ Find the maximum and minimum of the quadratic form f(x,y,z)f(x, y, z) , if the variables x,yx, y , and zz satisfy x2+y2+z2=1x^2 + y^2 + z^2 = 1 . ]

(e) 前問 (d) で求めた f(x,y,z)f(x,y,z) の最大値および最小値を与える x,y,zx,y,z の値を一組ずつ求めよ. [ Find the values of x,yx, y , and zz for the maximum and minimum of f(x,y,z)f (x, y, z) obtained in the problem (d), respectively. ]

题目描述

  1. 考虑二次型
f(x,y,z)=3x2+3y2+3z2+2xy+2yz+2zx.f(x,y,z)=3x^2+3y^2+3z^2+2xy+2yz+2zx.

回答下列问题。

(a) 使用向量

x=(x,y,z)T\mathbf{x}=(x,y,z)^T

与对称矩阵 AA,将二次型表示为

f(x,y,z)=xTAx.f(x,y,z)=\mathbf{x}^TA\mathbf{x}.

AA。这里上标 TT 表示转置。

(b) 求 (a) 中对称矩阵 AA 的特征值与特征向量。

(c) 使用正交矩阵将 (a) 中的对称矩阵 AA 对角化。

(d) 当变量 x,y,zx,y,z 满足

x2+y2+z2=1x^2+y^2+z^2=1

时,求二次型 f(x,y,z)f(x,y,z) 的最大值与最小值。

(e) 对 (d) 中的最大值与最小值,分别求一组使 f(x,y,z)f(x,y,z) 取得该值的 x,y,zx,y,z

Kai

(a) Solution:

Let x=(xyz)x = \begin{pmatrix} x \\ y \\ z \end{pmatrix} . Then f(x,y,z)=3x2+3y2+3z2+2xy+2yz+2zx=xTAxf(x, y, z) = 3x^2 + 3y^2 + 3z^2 + 2xy + 2yz + 2zx = x^T A x , where AA is a symmetric matrix.

xTAx=(xyz)(abcbdecef)(xyz)x^T A x = \begin{pmatrix} x & y & z \end{pmatrix} \begin{pmatrix} a & b & c \\ b & d & e \\ c & e & f \end{pmatrix} \begin{pmatrix} x \\ y \\ z \end{pmatrix}
=ax2+dy2+fz2+2bxy+2cxz+2eyz= ax^2 + dy^2 + fz^2 + 2bxy + 2cxz + 2eyz

Comparing the coefficients, we have: a=3,d=3,f=3,b=1,c=1,e=1a = 3, d = 3, f = 3, b = 1, c = 1, e = 1 Thus, the matrix AA is:

A=(311131113)A = \begin{pmatrix} 3 & 1 & 1 \\ 1 & 3 & 1 \\ 1 & 1 & 3 \end{pmatrix}

(b) Solution:

To find the eigenvalues of AA , we need to solve the characteristic equation AλI=0|A - \lambda I| = 0 .

AλI=3λ1113λ1113λ=(3λ)[(3λ)21]1[(3λ)1]+1[1(3λ)]=0|A - \lambda I| = \begin{vmatrix} 3-\lambda & 1 & 1 \\ 1 & 3-\lambda & 1 \\ 1 & 1 & 3-\lambda \end{vmatrix} = (3-\lambda)[(3-\lambda)^2 - 1] - 1[(3-\lambda)-1] + 1[1-(3-\lambda)] = 0
(3λ)((3λ)21)(2λ)(2λ)=0(3-\lambda)((3-\lambda)^2 - 1) - (2-\lambda) - (2-\lambda) = 0
(3λ)(96λ+λ21)4+2λ=0(3-\lambda)(9 - 6\lambda + \lambda^2 - 1) - 4 + 2\lambda = 0
(3λ)(λ26λ+8)4+2λ=0(3-\lambda)(\lambda^2 - 6\lambda + 8) - 4 + 2\lambda = 0
3λ218λ+24λ3+6λ28λ4+2λ=03\lambda^2 - 18\lambda + 24 - \lambda^3 + 6\lambda^2 - 8\lambda - 4 + 2\lambda = 0
λ3+9λ224λ+20=0-\lambda^3 + 9\lambda^2 - 24\lambda + 20 = 0
λ39λ2+24λ20=0\lambda^3 - 9\lambda^2 + 24\lambda - 20 = 0
(λ2)(λ27λ+10)=0(\lambda - 2)(\lambda^2 - 7\lambda + 10) = 0
(λ2)(λ2)(λ5)=0(\lambda - 2)(\lambda - 2)(\lambda - 5) = 0

So, the eigenvalues are λ1=2,λ2=2,λ3=5\lambda_1 = 2, \lambda_2 = 2, \lambda_3 = 5 .

For λ=2\lambda = 2 :

(A2I)v=0(111111111)(xyz)=(000)(A - 2I)v = 0 \Rightarrow \begin{pmatrix} 1 & 1 & 1 \\ 1 & 1 & 1 \\ 1 & 1 & 1 \end{pmatrix} \begin{pmatrix} x \\ y \\ z \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \\ 0 \end{pmatrix}
x+y+z=0z=xyx + y + z = 0 \Rightarrow z = -x - y

The eigenvectors are of the form v=(xyxy)=x(101)+y(011)v = \begin{pmatrix} x \\ y \\ -x-y \end{pmatrix} = x\begin{pmatrix} 1 \\ 0 \\ -1 \end{pmatrix} + y\begin{pmatrix} 0 \\ 1 \\ -1 \end{pmatrix} . Let v1=(101)v_1 = \begin{pmatrix} 1 \\ 0 \\ -1 \end{pmatrix} and v2=(011)v_2 = \begin{pmatrix} 0 \\ 1 \\ -1 \end{pmatrix} . These are linearly independent, but not orthogonal. We use Gram-Schmidt to orthogonalize them: u1=v1=(101)u_1 = v_1 = \begin{pmatrix} 1 \\ 0 \\ -1 \end{pmatrix} u2=v2v2u1u1u1u1=(011)12(101)=(1/211/2)u_2 = v_2 - \frac{v_2 \cdot u_1}{u_1 \cdot u_1} u_1 = \begin{pmatrix} 0 \\ 1 \\ -1 \end{pmatrix} - \frac{1}{2} \begin{pmatrix} 1 \\ 0 \\ -1 \end{pmatrix} = \begin{pmatrix} -1/2 \\ 1 \\ -1/2 \end{pmatrix} . We can take u2=(121)u_2 = \begin{pmatrix} -1 \\ 2 \\ -1 \end{pmatrix} . So the eigenvectors corresponding to λ=2\lambda = 2 are u1=(101)u_1 = \begin{pmatrix} 1 \\ 0 \\ -1 \end{pmatrix} and u2=(121)u_2 = \begin{pmatrix} -1 \\ 2 \\ -1 \end{pmatrix} .

For λ=5\lambda = 5 :

(A5I)v=0(211121112)(xyz)=(000)(A - 5I)v = 0 \Rightarrow \begin{pmatrix} -2 & 1 & 1 \\ 1 & -2 & 1 \\ 1 & 1 & -2 \end{pmatrix} \begin{pmatrix} x \\ y \\ z \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \\ 0 \end{pmatrix}

From the first equation, 2x+y+z=0-2x + y + z = 0 . From the second equation, x2y+z=0x - 2y + z = 0 . Subtracting the two equations gives 3x+3y=0x=y-3x + 3y = 0 \Rightarrow x = y . Substituting x=yx=y into the first equation: 2x+x+z=0z=x-2x + x + z = 0 \Rightarrow z = x . So the eigenvector is of the form v=(xxx)v = \begin{pmatrix} x \\ x \\ x \end{pmatrix} . We can take v3=(111)v_3 = \begin{pmatrix} 1 \\ 1 \\ 1 \end{pmatrix} .

(c) Solution:

We normalize the eigenvectors obtained in (b): u1=(101)12(101)u_1 = \begin{pmatrix} 1 \\ 0 \\ -1 \end{pmatrix} \Rightarrow \frac{1}{\sqrt{2}} \begin{pmatrix} 1 \\ 0 \\ -1 \end{pmatrix} u2=(121)16(121)u_2 = \begin{pmatrix} -1 \\ 2 \\ -1 \end{pmatrix} \Rightarrow \frac{1}{\sqrt{6}} \begin{pmatrix} -1 \\ 2 \\ -1 \end{pmatrix} v3=(111)13(111)v_3 = \begin{pmatrix} 1 \\ 1 \\ 1 \end{pmatrix} \Rightarrow \frac{1}{\sqrt{3}} \begin{pmatrix} 1 \\ 1 \\ 1 \end{pmatrix} The orthogonal matrix PP is formed by these normalized eigenvectors:

P=(12161302613121613)P = \begin{pmatrix} \frac{1}{\sqrt{2}} & \frac{-1}{\sqrt{6}} & \frac{1}{\sqrt{3}} \\ 0 & \frac{2}{\sqrt{6}} & \frac{1}{\sqrt{3}} \\ \frac{-1}{\sqrt{2}} & \frac{-1}{\sqrt{6}} & \frac{1}{\sqrt{3}} \end{pmatrix}

And PTAP=DP^T A P = D , where D=(200020005)D = \begin{pmatrix} 2 & 0 & 0 \\ 0 & 2 & 0 \\ 0 & 0 & 5 \end{pmatrix} .

(d) Solution:

Since f(x,y,z)=xTAxf(x, y, z) = x^T A x and we have diagonalized AA as PTAP=DP^T A P = D , let x=Pyx = Py , where y=(xyz)y = \begin{pmatrix} x' \\ y' \\ z' \end{pmatrix} . Then f(x,y,z)=yTDy=2x2+2y2+5z2f(x, y, z) = y^T D y = 2x'^2 + 2y'^2 + 5z'^2 . Also, x2+y2+z2=1(Py)T(Py)=1yTPTPy=1yTy=1x2+y2+z2=1x^2 + y^2 + z^2 = 1 \Rightarrow (Py)^T(Py) = 1 \Rightarrow y^T P^T P y = 1 \Rightarrow y^T y = 1 \Rightarrow x'^2 + y'^2 + z'^2 = 1 . We want to find the maximum and minimum of 2x2+2y2+5z22x'^2 + 2y'^2 + 5z'^2 subject to x2+y2+z2=1x'^2 + y'^2 + z'^2 = 1 . Since x2+y2=1z2x'^2 + y'^2 = 1 - z'^2 , f(x,y,z)=2(1z2)+5z2=2+3z2f(x, y, z) = 2(1 - z'^2) + 5z'^2 = 2 + 3z'^2 . Since 0z210 \leq z'^2 \leq 1 , the minimum value of f(x,y,z)f(x, y, z) is 22 (when z=0z' = 0 ) and the maximum value is 2+3=52 + 3 = 5 (when z=±1z' = \pm 1 ).

(e) Solution:

From (d), the maximum value of f(x,y,z)f(x, y, z) is 55 when z=±1z' = \pm 1 and x=y=0x' = y' = 0 . Then x=P(001)=(131313)x = P \begin{pmatrix} 0 \\ 0 \\ 1 \end{pmatrix} = \begin{pmatrix} \frac{1}{\sqrt{3}} \\ \frac{1}{\sqrt{3}} \\ \frac{1}{\sqrt{3}} \end{pmatrix} . Or x=(131313)x = \begin{pmatrix} -\frac{1}{\sqrt{3}} \\ -\frac{1}{\sqrt{3}} \\ -\frac{1}{\sqrt{3}} \end{pmatrix} . The minimum value of f(x,y,z)f(x, y, z) is 22 when z=0z' = 0 and x2+y2=1x'^2 + y'^2 = 1 . In this case, x=P(xy0)=(x2y62y6x2y6)x = P \begin{pmatrix} x' \\ y' \\ 0 \end{pmatrix} = \begin{pmatrix} \frac{x'}{\sqrt{2}} - \frac{y'}{\sqrt{6}} \\ \frac{2y'}{\sqrt{6}} \\ -\frac{x'}{\sqrt{2}} - \frac{y'}{\sqrt{6}} \end{pmatrix} . For example, when x=1,y=0x' = 1, y' = 0 , x=(12012)x = \begin{pmatrix} \frac{1}{\sqrt{2}} \\ 0 \\ -\frac{1}{\sqrt{2}} \end{pmatrix} . When x=0,y=1x'=0, y'=1 , x=(162616)x = \begin{pmatrix} -\frac{1}{\sqrt{6}} \\ \frac{2}{\sqrt{6}} \\ -\frac{1}{\sqrt{6}} \end{pmatrix} .