Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

練習問題メモ 21(対角化)

問1

2 次の正方行列 A=(1322)A=\left(\begin{array}{ll}1 & 3 \\ 2 & 2\end{array}\right) を考える。

  1. AA は対角化可能であることを示せ。

  2. P1APP^{-1} A P が対角行列となるような正則行列 PP を 1 つ求めよ。

  1. AA は対角化可能であることを示せ。

行列AMn(R)A\in M_n(\mathbb{R})が対角化可能であるとは、AAの固有ベクトルがRn\mathbb{R}^nの基底であることと等しい。

(1) AAの固有ベクトルを求める

AAの固有値λ\lambdaは、

det(1λ322λ)=(1λ)(2λ)6=0    λ23λ+26=0    λ23λ4=0    (λ+1)(λ4)=0λ=1,4\det \begin{pmatrix} 1- \lambda & 3\\ 2 & 2 - \lambda \end{pmatrix} = (1- \lambda) (2 - \lambda) - 6 = 0\\ \iff \lambda^2 - 3\lambda + 2 - 6 = 0\\ \iff \lambda^2 - 3\lambda - 4 = 0\\ \iff (\lambda + 1) (\lambda - 4) = 0\\ \therefore \lambda = -1, 4

相異なる固有値とそれに属する固有ベクトルが存在するため、これらはRn\mathbb{R}^nの基底となる。(固有ベクトルが基底となっていれば対角化可能であるという)定理により、AAは対角化可能である。

  1. P1APP^{-1} A P が対角行列となるような正則行列 PP を 1 つ求めよ。

λ=1\lambda = -1のとき、

2x1+3x2=02 x_1 + 3 x_2 = 0

より、x1=32x2x_1 = -\frac{3}{2} x_2。よって固有空間は

V(1)={c(321)cR}V(-1) = \left . \left\{ c \begin{pmatrix} -\frac{3}{2} \\ 1 \end{pmatrix} \right | c\in \mathbb{R} \right\}

λ=4\lambda = 4のとき、

3x1=3x2=0-3x_1 = 3x_2 = 0

よってx1=x2x_1 = x_2。固有空間は

V(4)={c(11)cR}V(4) = \left . \left\{ c \begin{pmatrix} 1 \\ 1 \end{pmatrix} \right | c\in \mathbb{R} \right\}

各固有値に属する固有ベクトルを列ベクトルとして行列を作り

P=(32111)P = \begin{pmatrix} -\frac{3}{2} & 1\\ 1 & 1 \end{pmatrix}

とする。

P1AP=(1004)P^{-1}AP = \begin{pmatrix} -1 & 0\\ 0 & 4 \end{pmatrix}

であるため、PPP1APP^{-1}APが対角行列となる正則行列である。

Source
import sympy as sp
A = sp.Matrix([
    [1, 3],
    [2, 2]
])
P = sp.Matrix([
    [-3/2, 1],
    [1, 1]
])

P.inv() @ A @ P
Output
Loading...
Source
import sympy as sp
A = sp.Matrix([
    [1, 3],
    [2, 2]
])
print("λ=", A.eigenvals().keys())
print("x=", A.eigenvects())
Output
λ= dict_keys([4, -1])
x= [(-1, 1, [Matrix([
[-3/2],
[   1]])]), (4, 1, [Matrix([
[1],
[1]])])]

問2

次の文章は、「対角化可能なべき零行列は零行列である」ことの証明である。証明文中の \square を埋めよ。

[1]=[1] =固有値

[2]=[2] =正則

[3]=0[3] = 0

[4]=0[4] = 0

[5]=[5] =零行列

問3

3 次の正方行列 A=(111021001)A=\left(\begin{array}{lll}1 & 1 & 1 \\ 0 & 2 & 1 \\ 0 & 0 & 1\end{array}\right) は上三角行列なので、 AA の固有値 λ\lambdaλ=1\lambda=1 (重解),2の2個であることがわかる。

  1. AA は対角化可能であることを示せ。

  2. P1APP^{-1} A P が対角行列となるような正則行列 PP を 1 つ求めよ。

  1. AA は対角化可能であることを示せ。

λ=1\lambda=1のとき

(AλE)x=0(011011000)(x1x2x3)=(000)x2+x3=0x2=x3(A-\lambda E) x = 0\\ \to \begin{pmatrix} 0 & 1 & 1 \\ 0 & 1 & 1 \\ 0 & 0 & 0 \end{pmatrix} \begin{pmatrix} x_1 \\ x_2 \\ x_3 \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \\ 0 \end{pmatrix} \\ \to x_2 + x_3 = 0\\ \to x_2 = -x_3

より、固有空間V(1)V(1)

V(1)={c1(100)+c2(011)c1,c2R}V(1)=\left\{\left. c_1 \begin{pmatrix} 1 \\ 0 \\ 0 \end{pmatrix} + c_2 \begin{pmatrix} 0 \\ -1 \\ 1 \end{pmatrix} \right\rvert\, c_1,c_2 \in \mathbb{R}\right\}

λ=2\lambda=2のとき

(AλE)x=0(111001001)(x1x2x3)=(000){x1+x2+x3=0x3=0x3=0(A-\lambda E) x = 0\\ \to \begin{pmatrix} -1 & 1 & 1 \\ 0 & 0 & 1 \\ 0 & 0 & -1 \end{pmatrix} \begin{pmatrix} x_1 \\ x_2 \\ x_3 \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \\ 0 \end{pmatrix} \\ \to \begin{cases} -x_1 + x_2 + x_3 = 0\\ x_3 = 0\\ -x_3 = 0 \end{cases} \\

より、

x1=x2+x30=x2x_1 = x_2 + \underbrace{ x_3 }_{0} = x_2

よって固有空間V(2)V(2)

V(2)={c(110)cR}V(2)=\left\{\left. c \begin{pmatrix} 1 \\ 1 \\ 0 \end{pmatrix} \right\rvert\, c \in \mathbb{R}\right\}

したがって、対角化可能であることの必要十分条件のひとつである

「固有空間の次元数の総和と行列AAの次元数が一致する」

あるいは

AAnn個の1次独立な固有ベクトルが存在する」

を満たすため、AAは対角化可能である。

  1. P1APP^{-1} A P が対角行列となるような正則行列 PP を 1 つ求めよ。

P=(101011010)P= \begin{pmatrix} 1 & 0 & 1\\ 0 & -1 & 1\\ 0 & 1 & 0 \end{pmatrix}
Source
import sympy as sp
A = sp.Matrix([
    [1, 1, 1],
    [0, 2, 1],
    [0, 0, 1]
])

P = sp.Matrix([
    [1, 0, 1],
    [0, -1, 1],
    [0, 1, 0]
])

P.inv() @ A @ P
Output
Loading...
Source
e = A.eigenvects()
print(e)
P = sp.Matrix(e[0][2] + e[1][2]).reshape(3,3).T
P.inv() @ A @ P
Output
[(1, 2, [Matrix([
[1],
[0],
[0]]), Matrix([
[ 0],
[-1],
[ 1]])]), (2, 1, [Matrix([
[1],
[1],
[0]])])]
Loading...

問4

2 次の正方行列 (ab0a)\left(\begin{array}{ll}a & b \\ 0 & a\end{array}\right) が対角化可能であるための必要十分条件は b=0b=0 であることを示せ。

A=(ab0a)A = \begin{pmatrix} a & b \\ 0 & a \end{pmatrix}

とおく。

b=0    b=0 \implies AAは対角化可能 を示す。 もしb=0b=0であれば

A=(a00a)A = \begin{pmatrix} a & 0 \\ 0 & a \end{pmatrix}

であり、AAはすでに対角行列であるため、2次の単位行列をPPとおけばP1APP^{-1} APは対角行列になり、AAは対角化可能である。

b0    b \neq 0 \implies AAは対角化不可能 を示す。

AAの固有値を求めると、

det(AλE)=det(aλb0aλ)=0    (aλ)2=0\det(A-\lambda E) = \det \begin{pmatrix} a-\lambda & b \\ 0 & a-\lambda \end{pmatrix} = 0\\ \iff (a-\lambda)^2 = 0

より、λ=a\lambda = aである。

AAが対角化可能であることと同値の条件のひとつは「AAの各固有値に対応する固有空間の次元の総和がAAの次数2と等しい」である。 AAの固有値λ\lambdaは1つ(λ=a\lambda=a)のみであるため、AAが対角化可能であれば、対応する固有空間V(a)V(a)の次元は2である。

λ=a\lambda = aのとき、

(AλE)=(0b00)(x1x2)=(00)(A -\lambda E) = \begin{pmatrix} 0 & b \\ 0 & 0 \end{pmatrix} \begin{pmatrix} x_1 \\ x_2 \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \end{pmatrix}

より、

bx2=0bx_2 = 0

となる。 もしb0b\neq 0であればx2=0x_2 = 0となり、固有空間V(λ)V(\lambda)の次元は

dimV(a)=dim{c(10)cR}=1\operatorname{\dim} V(a) = \operatorname{\dim} \left\{\left. c \begin{pmatrix} 1 \\ 0 \end{pmatrix} \right\rvert\, c \in \mathbb{R}\right\} = 1

となり、AAの次元数とは異なるため対角化可能にならない。

間違いではないが冗長なのでボツ

AAが対角化可能     b=0\implies b = 0を示す。

AAが対角化可能であることと同値の条件のひとつは「AAの各固有値に対応する固有空間の次元の総和がAAの次数2と等しい」である。

AAの固有値λ\lambdaは1つ(λ=a\lambda=a)のみであるため、AAが対角化可能であれば、対応する固有空間V(a)V(a)の次元は2である。

λ=a\lambda = aのとき、

(AλE)=(0b00)(x1x2)=(00)(A -\lambda E) = \begin{pmatrix} 0 & b \\ 0 & 0 \end{pmatrix} \begin{pmatrix} x_1 \\ x_2 \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \end{pmatrix}

より、

bx2=0bx_2 = 0

となる。

もしb0b\neq 0であればx2=0x_2 = 0となり、固有空間V(λ)V(\lambda)の次元は

dimV(a)=dim{c(10)cR}=1\operatorname{\dim} V(a) = \operatorname{\dim} \left\{\left. c \begin{pmatrix} 1 \\ 0 \end{pmatrix} \right\rvert\, c \in \mathbb{R}\right\} = 1

もしb=0b = 0であればx2x_2が0かどうかは未知であり、固有空間V(λ)V(\lambda)の次元は

dimV(a)=dim{c1(10)+c2(01)c1,c2R}=2\operatorname{\dim} V(a) = \operatorname{\dim} \left\{\left. c_1 \begin{pmatrix} 1 \\ 0 \end{pmatrix} + c_2 \begin{pmatrix} 0 \\ 1 \end{pmatrix} \right\rvert\, c_1,c_2 \in \mathbb{R}\right\} = 2

となる。よって

AAが対角化可能     b=0\implies b = 0

である。

以上により、

b=0    b=0 \implies AAは対角化可能

かつ

AAが対角化可能     b=0\implies b = 0

であるため、b=0b=0AAが対角化可能であるための必要十分条件(同値の条件)である

Source
import sympy as sp
A = sp.Matrix([
    [1, 0],
    [0, 1],
])

P = A

P.inv() @ A @ P
Output
Loading...