跳到主要内容

東京大学 新領域創成科学研究科 人間環境学専攻 2022年8月実施 线性代数

Author

思齐塾, 祭音Myyura

Description

行列 AA について以下の全ての設問に解答せよ。

A=[201020112]A = \begin{bmatrix} 2 & 0 & -1 \\ 0 & 2 & 0 \\ -1 & 1 & 2 \end{bmatrix}

(1) 行列 AA の固有値および固有ベクトルを全て求めよ。 (2) 行列 AA に対し P1APP^{-1}AP が対角行列となる正則行列 PP 、および、 P1APP^{-1}AP を求めよ。 (3) 行列 AnA^n を求めよ。ただし、 nn は自然数とする。

题目描述

给定矩阵

A=[201020112].A=\begin{bmatrix} 2&0&-1\\ 0&2&0\\ -1&1&2 \end{bmatrix}.

回答以下全部问题:

  1. 求矩阵 AA 的所有特征值及各自对应的特征向量。
  2. 求一个可逆矩阵 PP,使 P1APP^{-1}AP 为对角矩阵,并明确写出 P1APP^{-1}AP
  3. nn 为自然数,求 AnA^n 的显式表达式。

Kai

(1) 固有値の計算 Characteristic polynomial AλI=0|A - \lambda I| = 0 より、 2λ0102λ0112λ=(2λ)2λ112λ=(2λ)((2λ)21)=(2λ)(44λ+λ21)=(2λ)(λ24λ+3)=(2λ)(λ1)(λ3)=0\begin{vmatrix} 2-\lambda & 0 & -1 \\ 0 & 2-\lambda & 0 \\ -1 & 1 & 2-\lambda \end{vmatrix} = (2-\lambda)\begin{vmatrix} 2-\lambda & -1 \\ -1 & 2-\lambda \end{vmatrix} = (2-\lambda)((2-\lambda)^2 - 1) = (2-\lambda)(4 - 4\lambda + \lambda^2 - 1) = (2-\lambda)(\lambda^2 - 4\lambda + 3) = (2-\lambda)(\lambda-1)(\lambda-3) = 0 固有値は λ=1,2,3\lambda = 1, 2, 3

固有ベクトルの計算 λ=1\lambda = 1 の時、 (AI)v=0(A - I)v = 0 を解く。 [101010111][xyz]=[000]\begin{bmatrix} 1 & 0 & -1 \\ 0 & 1 & 0 \\ -1 & 1 & 1 \end{bmatrix} \begin{bmatrix} x \\ y \\ z \end{bmatrix} = \begin{bmatrix} 0 \\ 0 \\ 0 \end{bmatrix} より、 x=z,y=0x = z, y = 0 。固有ベクトルは v1=[101]v_1 = \begin{bmatrix} 1 \\ 0 \\ 1 \end{bmatrix} (またはその定数倍). λ=2\lambda = 2 の時、 (A2I)v=0(A - 2I)v = 0 を解く。 [001000110][xyz]=[000]\begin{bmatrix} 0 & 0 & -1 \\ 0 & 0 & 0 \\ -1 & 1 & 0 \end{bmatrix} \begin{bmatrix} x \\ y \\ z \end{bmatrix} = \begin{bmatrix} 0 \\ 0 \\ 0 \end{bmatrix} より、 z=0,x=yz = 0, x = y 。固有ベクトルは v2=[110]v_2 = \begin{bmatrix} 1 \\ 1 \\ 0 \end{bmatrix} (またはその定数倍). λ=3\lambda = 3 の時、 (A3I)v=0(A - 3I)v = 0 を解く。 [101010111][xyz]=[000]\begin{bmatrix} -1 & 0 & -1 \\ 0 & -1 & 0 \\ -1 & 1 & -1 \end{bmatrix} \begin{bmatrix} x \\ y \\ z \end{bmatrix} = \begin{bmatrix} 0 \\ 0 \\ 0 \end{bmatrix} より、 x=z,y=0x = -z, y = 0 。固有ベクトルは v3=[101]v_3 = \begin{bmatrix} 1 \\ 0 \\ -1 \end{bmatrix} (またはその定数倍).

(2) P=[111010101]P = \begin{bmatrix} 1 & 1 & 1 \\ 0 & 1 & 0 \\ 1 & 0 & -1 \end{bmatrix} とすると、 P1AP=[100020003]P^{-1}AP = \begin{bmatrix} 1 & 0 & 0 \\ 0 & 2 & 0 \\ 0 & 0 & 3 \end{bmatrix} (対角行列).

P1P^{-1} を求める。 [111100010010101001][111100010010012101][101110010010002111][1001/21/21/20100100011/21/21/2]\begin{bmatrix} 1 & 1 & 1 & | & 1 & 0 & 0 \\ 0 & 1 & 0 & | & 0 & 1 & 0 \\ 1 & 0 & -1 & | & 0 & 0 & 1 \end{bmatrix} \to \begin{bmatrix} 1 & 1 & 1 & | & 1 & 0 & 0 \\ 0 & 1 & 0 & | & 0 & 1 & 0 \\ 0 & -1 & -2 & | & -1 & 0 & 1 \end{bmatrix} \to \begin{bmatrix} 1 & 0 & 1 & | & 1 & -1 & 0 \\ 0 & 1 & 0 & | & 0 & 1 & 0 \\ 0 & 0 & -2 & | & -1 & 1 & 1 \end{bmatrix} \to \begin{bmatrix} 1 & 0 & 0 & | & 1/2 & -1/2 & 1/2 \\ 0 & 1 & 0 & | & 0 & 1 & 0 \\ 0 & 0 & 1 & | & 1/2 & -1/2 & -1/2 \end{bmatrix} P1=[1/21/21/20101/21/21/2]P^{-1} = \begin{bmatrix} 1/2 & -1/2 & 1/2 \\ 0 & 1 & 0 \\ 1/2 & -1/2 & -1/2 \end{bmatrix}

直接計算すると

P1AP=[100020003].P^{-1}AP= \begin{bmatrix} 1&0&0\\ 0&2&0\\ 0&0&3 \end{bmatrix}.

(3) An=(PDP1)n=PDnP1A^n = (PDP^{-1})^n = PD^nP^{-1} (where D=[100020003]D = \begin{bmatrix} 1 & 0 & 0 \\ 0 & 2 & 0 \\ 0 & 0 & 3 \end{bmatrix} )

Dn=[1n0002n0003n]=[10002n0003n]D^n = \begin{bmatrix} 1^n & 0 & 0 \\ 0 & 2^n & 0 \\ 0 & 0 & 3^n \end{bmatrix} = \begin{bmatrix} 1 & 0 & 0 \\ 0 & 2^n & 0 \\ 0 & 0 & 3^n \end{bmatrix}

An=PDnP1=[111010101][10002n0003n][1/21/21/20101/21/21/2]=[12n3n02n0103n][1/21/21/20101/21/21/2]=[(1+3n)/2(1+22n3n)/2(13n)/202n0(13n)/2(3n1)/2(1+3n)/2]A^n = PD^nP^{-1} = \begin{bmatrix} 1 & 1 & 1 \\ 0 & 1 & 0 \\ 1 & 0 & -1 \end{bmatrix} \begin{bmatrix} 1 & 0 & 0 \\ 0 & 2^n & 0 \\ 0 & 0 & 3^n \end{bmatrix} \begin{bmatrix} 1/2 & -1/2 & 1/2 \\ 0 & 1 & 0 \\ 1/2 & -1/2 & -1/2 \end{bmatrix} = \begin{bmatrix} 1 & 2^n & 3^n \\ 0 & 2^n & 0 \\ 1 & 0 & -3^n \end{bmatrix} \begin{bmatrix} 1/2 & -1/2 & 1/2 \\ 0 & 1 & 0 \\ 1/2 & -1/2 & -1/2 \end{bmatrix} = \begin{bmatrix} (1+3^n)/2 & (-1 + 2*2^n - 3^n)/2 & (1-3^n)/2 \\ 0 & 2^n & 0 \\ (1-3^n)/2 & (3^n-1)/2 & (1+3^n)/2 \end{bmatrix}