행렬의 스미스 노멀 폼
📂위상데이터분석 행렬의 스미스 노멀 폼 정의 PID, 즉 주아이디얼정역 R R R 이 주어져 있다고 하자. 행렬 A ∈ R m × n A \in R^{m \times n} A ∈ R m × n 에 대해 다음을 만족시키는 d 1 , ⋯ , d r ∈ R d_{1} , \cdots , d_{r} \in R d 1 , ⋯ , d r ∈ R 과 가역행렬 P ∈ R m × m P \in R^{m \times m} P ∈ R m × m , Q ∈ R n × n Q \in R^{n \times n} Q ∈ R n × n 이 존재하면 P A Q ∈ R m × n PAQ \in R^{m \times n} P A Q ∈ R m × n 을 A A A 의 스미스 노멀 폼 smith Normal form 이라 한다.
P A Q = [ d 1 0 0 0 ⋯ 0 0 ⋱ 0 0 ⋯ 0 0 0 d r 0 ⋯ 0 0 0 0 0 ⋯ 0 ⋮ ⋮ ⋮ ⋮ ⋱ ⋮ 0 0 0 0 ⋯ 0 ] ∈ R m × n
PAQ = \begin{bmatrix}
d_{1} & 0 & 0 & 0 & \cdots & 0
\\ 0 & \ddots & 0 & 0 & \cdots & 0
\\ 0 & 0 & d_{r} & 0 & \cdots & 0
\\ 0 & 0 & 0 & 0 & \cdots & 0
\\ \vdots & \vdots & \vdots & \vdots & \ddots & \vdots
\\ 0 & 0 & 0 & 0 & \cdots & 0
\end{bmatrix} \in R^{m \times n}
P A Q = d 1 0 0 0 ⋮ 0 0 ⋱ 0 0 ⋮ 0 0 0 d r 0 ⋮ 0 0 0 0 0 ⋮ 0 ⋯ ⋯ ⋯ ⋯ ⋱ ⋯ 0 0 0 0 ⋮ 0 ∈ R m × n
여기서 d 1 , ⋯ , d r d_{1} , \cdots , d_{r} d 1 , ⋯ , d r 은 R R R 에서 0 0 0 이 아니어야 하며, d 1 ∣ ⋯ ∣ d r d_{1} \mid \cdots \mid d_{r} d 1 ∣ ⋯ ∣ d r 이어야한다. 다시 말해 d k ≠ 0 d_{k} \ne 0 d k = 0 는 d k + 1 d_{k+1} d k + 1 의 인자divisor 여야한다.
예시 역행렬이 나와서 복잡해보이지만 실상은 가우스 소거법 을 통해서 구할 수 있다. 공식적으로 쓰는 표현은 절대 아니지만, 예시를 읽기 쉽도록 행렬 우측에 어떤 연산을 취했는지 표기하도록 하겠다. r − r- r − 는 행row , c − c- c − 는 열column 이고 − x -x − x 는 교환exchange , − e -e − e 는 소거elimination , − + -+ − + 는 해당 계수를 곱해서 더하는 연산이다.
물론 가우스 소거법을 통해 계산하는 것은 계산 측면에서 별로 현명하지 않을지도 모른다. 그냥 손계산 할 때는 이렇게도 구할 수 있구나 하는 정도로만 참고하자.
정수환 Z \mathbb{Z} Z [ 2 3 1 4 − 1 − 1 ] ∼ [ 1 3 2 − 1 − 1 4 ] c x ( 1 , 3 ) ∼ [ 1 0 0 − 1 2 6 ] c e ( 1 → 2 , 3 ) ∼ [ 1 0 0 0 2 6 ] r e ( 1 → 2 ) ∼ [ 1 0 0 0 2 0 ] c e ( 2 → 3 )
\begin{align*}
\begin{bmatrix}
2 & 3 & 1
\\ 4 &-1 &-1
\end{bmatrix} \sim & & \begin{bmatrix}
1 & 3 & 2
\\ -1 &-1 & 4
\end{bmatrix} & \qquad cx(1,3)
\\ \sim & & \begin{bmatrix}
1 & 0 & 0
\\ -1 & 2 & 6
\end{bmatrix} & \qquad ce(1 \to 2,3)
\\ \sim & & \begin{bmatrix}
1 & 0 & 0
\\ 0 & 2 & 6
\end{bmatrix} & \qquad re(1 \to 2)
\\ \sim & & \begin{bmatrix}
1 & 0 & 0
\\ 0 & 2 & 0
\end{bmatrix} & \qquad ce(2 \to 3)
\end{align*}
[ 2 4 3 − 1 1 − 1 ] ∼ ∼ ∼ ∼ [ 1 − 1 3 − 1 2 4 ] [ 1 − 1 0 2 0 6 ] [ 1 0 0 2 0 6 ] [ 1 0 0 2 0 0 ] c x ( 1 , 3 ) ce ( 1 → 2 , 3 ) re ( 1 → 2 ) ce ( 2 → 3 )
여기서 1 ∣ 2 1 \mid 2 1 ∣ 2 임에 주의하라.
다항식환 Q [ x ] \mathbb{Q} [x] Q [ x ] [ x 2 0 0 1 2 x − 1 0 x + 1 ] ∼ [ x 2 x 0 1 2 x − 1 0 x ] c + ( 1 → 3 , 1 ) ∼ [ x + 1 2 0 2 1 0 − 1 0 x ] r e ( 3 → 1 , 2 ) ∼ [ 2 1 0 x + 1 2 0 − 1 0 x ] r x ( 1 , 2 ) ∼ [ 1 2 0 2 x + 1 0 0 − 1 x ] c x ( 1 , 2 ) ∼ [ 1 2 0 0 x − 3 0 0 − 1 x ] r e ( 1 → 2 ) ∼ [ 1 0 0 0 x − 3 0 0 − 1 x ] c e ( 1 → 2 ) ∼ [ 1 0 0 0 1 − x 0 x − 3 0 ] r x ( 2 , 3 ) ∼ [ 1 0 0 0 1 0 0 x − 3 x ( x − 3 ) ] c + ( 2 → 3 , x ) ∼ [ 1 0 0 0 1 0 0 0 x ( x − 3 ) ] r + ( 2 → 3 , − ( x − 3 ) )
\begin{align*}
\begin{bmatrix}
x & 2 & 0
\\ 0 & 1 & 2x
\\ -1 & 0 & x+1
\end{bmatrix} \sim & \begin{bmatrix}
x & 2 & x
\\ 0 & 1 & 2x
\\ -1 & 0 & x
\end{bmatrix} & \qquad c+(1 \to 3, 1)
\\ \sim & \begin{bmatrix}
x+1 & 2 & 0
\\ 2 & 1 & 0
\\ -1 & 0 & x
\end{bmatrix} & \qquad re(3 \to 1, 2)
\\ \sim & \begin{bmatrix}
2 & 1 & 0
\\ x+1 & 2 & 0
\\ -1 & 0 & x
\end{bmatrix} & \qquad rx(1,2)
\\ \sim & \begin{bmatrix}
1 & 2 & 0
\\ 2 & x+1 & 0
\\ 0 & -1 & x
\end{bmatrix} & \qquad cx(1,2)
\\ \sim & \begin{bmatrix}
1 & 2 & 0
\\ 0 & x-3 & 0
\\ 0 & -1 & x
\end{bmatrix} & \qquad re(1 \to 2)
\\ \sim & \begin{bmatrix}
1 & 0 & 0
\\ 0 & x-3 & 0
\\ 0 & -1 & x
\end{bmatrix} & \qquad ce(1 \to 2)
\\ \sim & \begin{bmatrix}
1 & 0 & 0
\\ 0 & 1 & -x
\\ 0 & x-3 & 0
\end{bmatrix} & \qquad rx(2,3)
\\ \sim & \begin{bmatrix}
1 & 0 & 0
\\ 0 & 1 & 0
\\ 0 & x-3 & x(x-3)
\end{bmatrix} & \qquad c+(2 \to 3, x)
\\ \sim & \begin{bmatrix}
1 & 0 & 0
\\ 0 & 1 & 0
\\ 0 & 0 & x(x-3)
\end{bmatrix} & \qquad r+(2 \to 3, -(x-3))
\end{align*}
x 0 − 1 2 1 0 0 2 x x + 1 ∼ ∼ ∼ ∼ ∼ ∼ ∼ ∼ ∼ x 0 − 1 2 1 0 x 2 x x x + 1 2 − 1 2 1 0 0 0 x 2 x + 1 − 1 1 2 0 0 0 x 1 2 0 2 x + 1 − 1 0 0 x 1 0 0 2 x − 3 − 1 0 0 x 1 0 0 0 x − 3 − 1 0 0 x 1 0 0 0 1 x − 3 0 − x 0 1 0 0 0 1 x − 3 0 0 x ( x − 3 ) 1 0 0 0 1 0 0 0 x ( x − 3 ) c + ( 1 → 3 , 1 ) re ( 3 → 1 , 2 ) r x ( 1 , 2 ) c x ( 1 , 2 ) re ( 1 → 2 ) ce ( 1 → 2 ) r x ( 2 , 3 ) c + ( 2 → 3 , x ) r + ( 2 → 3 , − ( x − 3 ))
설명 주의해야할 것은 주어진 주아이디얼정역 R R R 은 체 라는 보장이 없기 때문에 일반적으로 나눗셈이 허용되지 않고, PID 가 베주 정역 임을 이용해서 소거를 해야한다는 것이다.
정리 존재성과 유일성 R R R 이 주아이디얼정역 일 때, 모든 행렬 A ∈ R m × n A \in R^{m \times n} A ∈ R m × n 에 대해 스미스 노멀 폼이 유일하게 존재한다. 이 때 유일성은 d k d_{k} d k 들의 단원 배를 무시한다.
증명 구체적으로 스미스 노멀 폼을 구하는 알고리즘을 통해 보인다.
■
같이보기