跳到主要内容

筑波大学 理工情報生命学術院 システム情報工学研究群 社会工学学位プログラム 2020年8月実施 线性代数

Author

思齐塾, 祭音Myyura

Description

数学

問題IとIIの両方に答えよ。問題ごとに別々の解答用紙を使用せよ。 以下では,実数全体の集合を R\mathbb{R} とする。

I. nn 次元実ベクトル aRna \in \mathbb{R}^n , nn 次正方行列 ARn×nA \in \mathbb{R}^{n \times n} に対して, ta{}^t a , tA{}^t A をそれぞれ aa の転置ベクトル, AA の転置行列とする。次の3つの3次元実ベクトル

p1=(101),p2=(110),p3=(010)p_1 = \begin{pmatrix} 1 \\ 0 \\ 1 \end{pmatrix}, p_2 = \begin{pmatrix} -1 \\ 1 \\ 0 \end{pmatrix}, p_3 = \begin{pmatrix} 0 \\ -1 \\ 0 \end{pmatrix}

に対して,3つの3次正方行列

P1=p1tp1=(101000101),P2=p2tp2=(110110000),P3=p3tp3=(000010000)P_1 = p_1 {}^t p_1 = \begin{pmatrix} 1 & 0 & 1 \\ 0 & 0 & 0 \\ 1 & 0 & 1 \end{pmatrix}, P_2 = p_2 {}^t p_2 = \begin{pmatrix} 1 & -1 & 0 \\ -1 & 1 & 0 \\ 0 & 0 & 0 \end{pmatrix}, P_3 = p_3 {}^t p_3 = \begin{pmatrix} 0 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 0 \end{pmatrix}

を定義し,さらに3つの実数 α1,α2,α3R\alpha_1, \alpha_2, \alpha_3 \in \mathbb{R} に対して3次正方行列

A(α1,α2,α3)=α1P1+α2P2+α3P3A(\alpha_1, \alpha_2, \alpha_3) = \alpha_1 P_1 + \alpha_2 P_2 + \alpha_3 P_3

を定義する。以下の(1)-(5)に答えよ。

(1) 以下の3次正方行列の各要素を計算せよ。

(tp1p1tp1p2tp1p3tp2p1tp2p2tp2p3tp3p1tp3p2tp3p3)\begin{pmatrix} {}^t p_1 p_1 & {}^t p_1 p_2 & {}^t p_1 p_3 \\ {}^t p_2 p_1 & {}^t p_2 p_2 & {}^t p_2 p_3 \\ {}^t p_3 p_1 & {}^t p_3 p_2 & {}^t p_3 p_3 \end{pmatrix}

(2) 3つの3次元実ベクトル A(α1,α2,α3)p1A(\alpha_1, \alpha_2, \alpha_3) p_1 , A(α1,α2,α3)p2A(\alpha_1, \alpha_2, \alpha_3) p_2 , A(α1,α2,α3)p3A(\alpha_1, \alpha_2, \alpha_3) p_3 それぞれを, α1,α2,α3\alpha_1, \alpha_2, \alpha_3 と具体的な実数で表される x1,x2,x3x_1, x_2, x_3 を用いて, x1p1+x2p2+x3p3x_1 p_1 + x_2 p_2 + x_3 p_3 の形で表せ。

(3) 3次正方行列 A(1,2,3)A(1, 2, 3) の固有値をすべて求めよ。

(4) nn 次正方行列 AA について,任意の nn 次元実ベクトル dRnd \in \mathbb{R}^n に対して tdAd>0{}^t d A d > 0 が成り立つとき, AA は正定値行列であるという。3次正方行列 A(1,2,3)A(1, 2, 3) が正定値行列であるか否か,理由を付して述べよ。

(5) 3次正方行列 A(α1,α2,α3)A(\alpha_1, \alpha_2, \alpha_3) が半正定値行列であるための, 3つの実数 α1,α2,α3R\alpha_1, \alpha_2, \alpha_3 \in \mathbb{R} に関する必要十分条件について,理由を付して述べよ。

题目描述

原题开头要求问题 I、II 均作答并分别使用答题纸,但当前文件的 Description 与 Kai 只包含问题 I;问题 II 的题面缺失,无法由本文件唯一恢复。以下为现有问题 I。

nn 维实向量 aRna\in\mathbb R^nnn 阶实方阵 ARn×nA\in\mathbb R^{n\times n},分别以 ta{}^tatA{}^tA 表示其转置。给定

p1=(101),p2=(110),p3=(010),p_1=\begin{pmatrix}1\\0\\1\end{pmatrix},\qquad p_2=\begin{pmatrix}-1\\1\\0\end{pmatrix},\qquad p_3=\begin{pmatrix}0\\-1\\0\end{pmatrix},

并定义三个外积矩阵

P1=p1tp1=(101000101),P_1=p_1{}^tp_1= \begin{pmatrix}1&0&1\\0&0&0\\1&0&1\end{pmatrix},
P2=p2tp2=(110110000),P3=p3tp3=(000010000).P_2=p_2{}^tp_2= \begin{pmatrix}1&-1&0\\-1&1&0\\0&0&0\end{pmatrix},\qquad P_3=p_3{}^tp_3= \begin{pmatrix}0&0&0\\0&1&0\\0&0&0\end{pmatrix}.

α1,α2,α3R\alpha_1,\alpha_2,\alpha_3\in\mathbb R,再令

A(α1,α2,α3)=α1P1+α2P2+α3P3.A(\alpha_1,\alpha_2,\alpha_3) =\alpha_1P_1+\alpha_2P_2+\alpha_3P_3.

回答下列问题:

  1. 计算 Gram 矩阵

    (tp1p1tp1p2tp1p3tp2p1tp2p2tp2p3tp3p1tp3p2tp3p3)\begin{pmatrix} {}^tp_1p_1&{}^tp_1p_2&{}^tp_1p_3\\ {}^tp_2p_1&{}^tp_2p_2&{}^tp_2p_3\\ {}^tp_3p_1&{}^tp_3p_2&{}^tp_3p_3 \end{pmatrix}

    的每个元素。

  2. 分别将

    A(α1,α2,α3)p1,A(α1,α2,α3)p2,A(α1,α2,α3)p3A(\alpha_1,\alpha_2,\alpha_3)p_1,\quad A(\alpha_1,\alpha_2,\alpha_3)p_2,\quad A(\alpha_1,\alpha_2,\alpha_3)p_3

    表示成 x1p1+x2p2+x3p3x_1p_1+x_2p_2+x_3p_3 的形式,其中 x1,x2,x3x_1,x_2,x_3 要用 α1,α2,α3\alpha_1,\alpha_2,\alpha_3 和具体实数写出。

  3. A(1,2,3)A(1,2,3) 的全部特征值。

  4. 判断 A(1,2,3)A(1,2,3) 是否为正定矩阵,并说明理由。题面以 tdAd>0{}^tdAd>0 定义正定性;按 Kai 可唯一确认的标准含义,此条件针对任意非零 dRnd\in\mathbb R^n(若含 d=0d=0,严格不等式不可能成立)。

  5. 给出 A(α1,α2,α3)A(\alpha_1,\alpha_2,\alpha_3) 为半正定矩阵时,α1,α2,α3\alpha_1,\alpha_2,\alpha_3 所满足的必要充分条件,并说明理由。

Kai

(1)

(tp1p1tp1p2tp1p3tp2p1tp2p2tp2p3tp3p1tp3p2tp3p3)=(210121011).\begin{pmatrix} {}^t p_1p_1 & {}^t p_1p_2 & {}^t p_1p_3\\ {}^t p_2p_1 & {}^t p_2p_2 & {}^t p_2p_3\\ {}^t p_3p_1 & {}^t p_3p_2 & {}^t p_3p_3 \end{pmatrix} = \begin{pmatrix} 2 & -1 & 0\\ -1 & 2 & -1\\ 0 & -1 & 1 \end{pmatrix}.

(2)

A(α1,α2,α3)p1=2α1p1α2p2,A(α1,α2,α3)p2=α1p1+2α2p2α3p3,A(α1,α2,α3)p3=α2p2+α3p3.\begin{aligned} A(\alpha_1,\alpha_2,\alpha_3)p_1 &=2\alpha_1p_1-\alpha_2p_2,\\ A(\alpha_1,\alpha_2,\alpha_3)p_2 &=-\alpha_1p_1+2\alpha_2p_2-\alpha_3p_3,\\ A(\alpha_1,\alpha_2,\alpha_3)p_3 &=-\alpha_2p_2+\alpha_3p_3. \end{aligned}

実際, P=(p1,p2,p3)P=(p_1,p_2,p_3) とおけば A=Pdiag(α1,α2,α3)tPA=P\operatorname{diag}(\alpha_1,\alpha_2,\alpha_3){}^tP であり, 上の係数は (1) の内積から得られる。

(3)

α1=1,α2=2,α3=3\alpha_1=1,\alpha_2=2,\alpha_3=3 を代入すると

A(1,2,3)=(321250101).A(1,2,3) = \begin{pmatrix} 3 & -2 & 1\\ -2 & 5 & 0\\ 1 & 0 & 1 \end{pmatrix} .

その固有方程式は

λ39λ2+18λ6=0\lambda^3-9\lambda^2+18\lambda-6=0

θ=13arccos(1/3)\theta=\frac13\arccos(1/\sqrt3) とおけば,3つの固有値は

λk=3+23cos(θ+2kπ3)(k=0,1,2)\lambda_k=3+2\sqrt3\cos\left(\theta+\frac{2k\pi}{3}\right) \qquad(k=0,1,2)

である(近似値は 0.415775,2.294280,6.2899450.415775,\,2.294280,\,6.289945 )。

(4)

対称行列 A(1,2,3)A(1,2,3) の首座主小行列式は

3,3225=11,detA(1,2,3)=63,\qquad \begin{vmatrix}3&-2\\-2&5\end{vmatrix}=11,\qquad \det A(1,2,3)=6

であり,すべて正である。したがって Sylvester の判定法により A(1,2,3)A(1,2,3) は正定値行列である。

(5)

P=(p1,p2,p3)P=(p_1,p_2,p_3) とおくと detP=1\det P=1 なので PP は正則であり,

A(α1,α2,α3)=Pdiag(α1,α2,α3)tPA(\alpha_1,\alpha_2,\alpha_3) =P\operatorname{diag}(\alpha_1,\alpha_2,\alpha_3){}^tP

と書ける。任意の dR3d\in\mathbb R^3 に対して

tdA(α1,α2,α3)d=α1(tp1d)2+α2(tp2d)2+α3(tp3d)2.{}^tdA(\alpha_1,\alpha_2,\alpha_3)d =\alpha_1({}^tp_1d)^2+\alpha_2({}^tp_2d)^2+\alpha_3({}^tp_3d)^2.

PP が正則であるから, (tp1d,tp2d,tp3d)({}^tp_1d,{}^tp_2d,{}^tp_3d)R3\mathbb R^3 のすべての値をとる。よって半正定値であるための必要十分条件は

α10,α20,α30\alpha_1\ge 0,\quad \alpha_2\ge 0,\quad \alpha_3\ge 0

である。