How To Compute The Inverse of a Matrix

2023. 9. 10. 13:53·Mathematics/Linear Algebra

Inverse of a Matrix

Definition 1. Let $A \in M_{n \times n}(F)$. Then $A$ is invertible if $\exists B \in M_{n \times n}(F)$ such that $AB = BA = I$. The matrix $B$ is called the inverse of $A$, denoted $A^{-1}$.

    위와 같이 정의되는 역행렬을 구하는 방법은 여러가지가 있으나, 여기서는 elementary operation을 이용하여 구하는 방법만을 다룬다.

Augmented Matrix

Definition 2. Let $A \in M_{m \times n}(F)$ and $B \in M_{m \times p}(F)$. Then the augmented matrix $(A \vert B)$ is the $m \times (n + p)$ matrix $(A \,\, B)$.

    위와 같이 두 개의 행렬을 붙여서 만든 행렬을 augmented matrix, 첨가 행렬이라고 부른다. 역행렬을 구하고자 하는 행렬과 동일한 크기의 단위행렬로 만든 첨가 행렬을 elementary row operation만을 이용해 변환하면 역행렬을 구할 수 있다. 아래의 정리가 그것을 보장해준다. 

Lemma

Lemma. Suppose that $A$ and $B$ are matrices having $n$ rows. Then $M(A \vert B) = (MA \vert MB), \forall M \in M_{m \times n}(F)$.
Proof. Note that $(A \vert B)_{ij} = \begin{cases} A_{ij} & \text{if } 1 \leq j \leq n \\ B_{ij} & \text{if } n+1 \leq j \leq p \end{cases}$
For $1 \leq j \leq n, (M(A \vert B))_{ij} = M_{ik}(A \vert B)_{kj} = M_{ik}A_{kj} = (MA)_{ij}.$
For $n+1 \leq j \leq p, (M(A \vert B))_{ij} = M_{ik}(A \vert B)_{kj} = M_{ik}B_{kj} = (MB)_{ij}.$
$\Longrightarrow (M(A \vert B))_{ij} = \begin{cases} (MA)_{ij} & \text{if } 1 \leq j \leq n \\ (MB)_{ij} & \text{if } n+1 \leq j \leq p \end{cases} = (MA \vert MB)_{ij}$. $\blacksquare$

Theorem 1

Theorem 1. Let $A \in M_{n \times n}(F)$. Then
(a) If $A$ is invertible, then $(A \vert I_n)$ can be transformed into $(I_n \vert A^{-1})$ by means of a finite number of elementary row operations,
(b) If $A$ is invertible and $(A \vert I_n)$ is transformed into $(I_n \vert B)$ by means of a finite number of elementary row operations, then $B = A^{-1}$,
(c) If $A$ is not invertible, then any attempt to transform $(A \vert I_n)$ into $(I_n \vert B)$ produces a row whose first $n$ entries are zeros.
Proof.
(a) Note that $A^{-1}$ is invertible so that $A^{-1}$ can be written as $A^{-1} = R_p(I_n) \cdots R_1(I_n)$. Then $R_p \cdots R_1(A \vert I_n) = R_p(I_n) \cdots R_1(I_n)(A \vert I_n)$ = $A^{-1}(A \vert I_n) = (A^{-1}A \vert A^{-1}I_n) = (I_n \vert A^{-1})$. Thus we can obtain $(I_n \vert A^{-1})$ from $(A \vert I_n)$ by means of a finite number of elementary row operations.
(b) Suppose that $R_q \cdots R_1(A \vert I_n) = (I_n \vert B)$. Then $R_q \cdots R_1(A \vert I_n) = R_p(I_n) \cdots R_1(I_n)(A \vert I_n)$. Define $B := R_p(I_n) \cdots R_1(I_n)$. Then $B(A \vert I_n) = (BA \vert BI_n) = (BA \vert B) = (I_n \vert B)$. Thus $BA = I_n$, so $B = A^{-1}$.
(c) Since $A$ is not invertible, rank($A$) < $n$. Then $(A \vert I_n)$ cannot be transformed into $(I_n \vert B)$ because elementary operations are rank-preserving (rank($A$) = rank($I_n$) = $n < n \bigotimes$). $\blacksquare$

Example

역행렬은 다양한 상황에서 계산되어 지는데, 그 중 하나로 선형 변환의 역함수, 즉 역변환을 구하는 방법을 알아보자. $T(f(x)) = f(x) + f'(x) + f''(x)$로 정의되는 선형 변환 $T \in \mathcal{L}(P_2(\mathbb{R}))$을 생각하자. 이때 역변환 $T^{-1}$을 구해보자. $P_2(\mathbb{R})$의 standard ordered basis $\beta = \{1, x, x^2\}$에 대한 $T$의 matrix representation은 다음과 같다. $$[T]_{\beta} = \begin{pmatrix} 1 & 1 & 2 \\ 0 & 1 & 2 \\ 0 & 0 & 1 \end{pmatrix}$$ Theorem 1을 사용하여 augmented matrix $([T]_{\beta} \,\, | \,\, I_3)$를 $(I_3 \,\, | \,\, B)$로 변환하면 $B = ([T]_{\beta})^{-1} = [T^{-1}]_{\beta}$이고 $$[T^{-1}]_{\beta} = \begin{pmatrix} 1 & -1 & 0 \\ 0 & 1 & -2 \\ 0 & 0 & 1 \end{pmatrix}$$이다. 따라서 $$\begin{align*} &T^{-1}(1) = 1 \\ &T^{-1}(x) = x - 1 \\ &T^{-1}(x^2) = x^2 - 2x \end{align*}$$이다. 때문에 역변환 $T^{-1}$은 $T^{-1}(a_0 + a_1 x + a_2 x^2) = (a_0 - a_1) + (a_1 - 2a_2)x + a_2 x^2$으로 계산된다. 이와 같은 방법으로 invertible linear transformation의 역변환을 찾을 수 있다.

Components of The Inverse of a Matrix

Theorem 2. Let $A \in M_{n}(F)$ be invertible. Then $$(A^{-1})_{ij} = \frac{(-1)^{i+j} \det(\widetilde{A_{ji}})}{\det(A)}.$$
저작자표시 (새창열림)
'Mathematics/Linear Algebra' 카테고리의 다른 글
  • A Homogeneous System of Linear Equations
  • A System of Linear Equations
  • How To Calculate The Rank of a Matrix
  • rank($AB$) $\leq$ rank($A$), rank($B$)
Erdos
Erdos
수학과, 물리학과 학부생들이 운영하는 팀블로그입니다.
  • Erdos
    SAMICO
    Erdos
  • 전체
    오늘
    어제
    • 분류 전체보기 (273) N
      • Mathematics (183)
        • Real analysis (30)
        • Linear Algebra (64)
        • Number Thoery (11)
        • Calculus (55)
        • Probability (6)
        • Set Theory (13)
        • Writing (2)
        • Problems (1)
        • Abstract Algebra (1)
      • Physics (71) N
        • 일반물리 (2)
        • 상대성이론과 양자역학 입문 (35)
        • 열물리 (15) N
        • 수리물리 (13)
        • 고전역학 (6)
      • Computer (7)
      • 독서 (12)
        • 과학 (5)
        • 문학 (2)
        • 자기계발서 (4)
  • 공지사항

    • 참고서적
  • 최근 댓글

  • 최근 글

  • hELLO· Designed By정상우.v4.10.2
Erdos
How To Compute The Inverse of a Matrix
상단으로

티스토리툴바