3.3 Cramer's Rule

Jaehyun_onelion·2023년 3월 17일
0

선형대수학

목록 보기
16/42

이번 포스트에서는 Cramer's rule에 대해서 알아보겠습니다.


1) Cramer's Rule


Cramer's rule을 정의하기 위해서는 하나의 notation 정의가 필요합니다.


  • Definition

for any n×nn \times n matrix AA and any b\boldsymbol{b} in Rn\mathbb R^n, let Ai(b)A_i(\boldsymbol b) be the matrix obtained from AA by replaceing column ii by the vector b\boldsymbol{b}

Ai(b)=[a1a2...ai1bai+1...an]A_i(\boldsymbol{b})=\begin{bmatrix}\boldsymbol{a_1} & \boldsymbol{a_2} & ...& \boldsymbol{a_{i-1}} & \boldsymbol{b} & \boldsymbol{a_{i+1}} & ... & \boldsymbol{a_n} \end{bmatrix}

Ai(b)A_i(\boldsymbol b) 는 matrix AAii 번째 column 대신 b\boldsymbol{b}를 넣은 새로운 matrix입니다.


  • Theorem : Cramer's Rule

Let AA be an invertible n×nn \times n matrix. For any b\boldsymbol{b} in Rn\mathbb R ^n, the unique solution x\boldsymbol{x} of Ax=bA\boldsymbol{x}=\boldsymbol{b} has entries given by

xi=detAi(b)detA,  i=1,2,...,nx_i = \frac{detA_i(\boldsymbol{b})}{detA}, \ \ i=1, 2, ..., n

Cramer's rule을 이용하면, determinant를 이용하여 linear system의 solution을 구할 수 있습니다.

(증명은 appendix 참고)


example

3x12x2=65x1+4x2=8\begin{aligned} 3x_1-2x_2 &= 6 \\ -5x_1+4x_2 &= 8 \end{aligned}

위 linear system을 matrix equation으로 바꾸면

[3254][x1x2]=[68]\begin{bmatrix} 3 & -2 \\ -5 & 4 \end{bmatrix} \begin{bmatrix} x_1 \\ x_2 \end{bmatrix} = \begin{bmatrix} 6 \\ 8 \end{bmatrix}

Cramer's rule을 적용하기 위해 detAi(b)detA_i(\boldsymbol{b})detAdetA를 게산하면

detA1(b)=6284=40 detA2(b)=3658=54 detA=3254=2detA_1(\boldsymbol{b}) =\begin{vmatrix}6 & -2 \\ 8 & 4 \end{vmatrix} =40 \\ \ \\ detA_2(\boldsymbol{b}) =\begin{vmatrix}3 & 6 \\ -5 & 8 \end{vmatrix} = 54 \\ \ \\ detA = \begin{vmatrix}3 & -2 \\ -5 & 4 \end{vmatrix} = 2

Cramer's rule을 적용하면

x1=detA1(b)detA=20 x2=detA2(b)detA=27x_1 = \frac{detA_1(\boldsymbol{b})}{detA} = 20 \\ \ \\ x_2 = \frac{detA_2(\boldsymbol{b})}{detA}=27

따라서 위 linear system의 solution은

x=[2027]\boldsymbol{x} = \begin{bmatrix} 20 \\ 27\end{bmatrix}

입니다.


2) A formula for A1A^{-1}


Cramer's rule을 적용하여 A1A^{-1}를 찾을 수 있습니다.


1) FInding A1A^{-1}


Let AA be an invertible n×nn\times n matrix. Then, the jjth column of A1A^{-1} is a vector x\boldsymbol{x} that satisfies

Ax=ejA\boldsymbol{x} = \boldsymbol{e_j}

where ej\boldsymbol{e_j} is the jjth column of identity matrix.

Inverse의 정의에 의해

AA1=IAA^{-1} = I

를 만족합니다. 따라서 위 matrix 식의 결과를 각 column별로 살펴보면

Ax=ejA\boldsymbol{x} = \boldsymbol{e_j}

을 만족하는 x\boldsymbol{x}A1A^{-1}jjth column이 되는 것을 알 수 있습니다.

위 방정식의 solution을 구할 때, Cramer's rule을 이용하면

xij=detAi(ej)detAx_{ij} = \frac{detA_i({\boldsymbol{e_j})}}{detA}

가 되고, xijx_{ij}A1A^{-1}(i,j)(i, j) entry가 됩니다.

Ai(ej)A_i(\boldsymbol{e_j})를 살펴보면

Ai(ej)=[a1...ai1ejai+1...an]=[00a1...ai11ai+1...an0]\begin{aligned} A_i(\boldsymbol{e_j}) &= \begin{bmatrix}\boldsymbol{a_1} & ... & \boldsymbol{a_{i-1}} & \boldsymbol{e_j} & \boldsymbol{a_{i+1}} & ... & \boldsymbol{a_n} \end{bmatrix} \\ \\ &=\begin{bmatrix} & & & 0 & & & \\ & & & 0 & & & \\ & & & \vdots & & & \\\boldsymbol{a_1} & ... & \boldsymbol{a_{i-1}} & 1 & \boldsymbol{a_{i+1}} & ... & \boldsymbol{a_n} \\ & & & \vdots & & & \\ & & & 0 & & & \end{bmatrix} \end{aligned}

입니다. 따라서 detAi(ej)detA_i(\boldsymbol{e_j})를 co-factor expansion을 이용하여 구할 때 ii 번째 column을 기준으로 구하게 됩니다. Ai(ej)A_i(\boldsymbol{e_j})(j,i)(j, i) 위치에서의 co-factor는 Ai(ej)A_i(\boldsymbol{e_j})에서 jj 번째 row, ii 번째 column을 제외한 matrix의 determinant입니다.

해당 determinant는 AA에서 jj 번째 row, ii 번째 column을 제외한 matrix의 determinant와 같습니다. 즉 Ai(ej)A_i(\boldsymbol{e_j})(i,j)(i, j) 위치에 해당하는 cofactor는 AA(i,j)(i, j) 위치에 해당하는 cofactor와 일치합니다. 따라서

detAi(ej)=CjidetA_i(\boldsymbol{e_j}) = C_{ji}

입니다.

A1A^{-1}(i,j)(i, j) entry에 해당하는 값을 알았으니, A1A^{-1}를 표현하면

A1=1detA[C11C21...Cn1C12C22...Cn2C1nC2n...Cnn]A^{-1}= \frac{1}{detA}\begin{bmatrix}C_{11} & C_{21} & ... & C_{n1} \\ C_{12} & C_{22} & ... & C_{n2} \\ \vdots & \vdots & \ddots & \vdots \\ C_{1n} & C_{2n} & ... & C_{nn}\end{bmatrix}

이 됩니다. 여기서 1detA\frac{1}{detA}를 제외한 matrix 부분

[C11C21...Cn1C12C22...Cn2C1nC2n...Cnn]\begin{bmatrix}C_{11} & C_{21} & ... & C_{n1} \\ C_{12} & C_{22} & ... & C_{n2} \\ \vdots & \vdots & \ddots & \vdots \\ C_{1n} & C_{2n} & ... & C_{nn}\end{bmatrix}

adjugate of AA라고 하고, adjAadjA로 표시합니다.

A1=1detAadjAA^{-1} = \frac{1}{detA}adjA

example

A=[213111142]A = \begin{bmatrix} 2 & 1 & 3 \\ 1 & -1 & 1 \\ 1 & 4 & 2\end{bmatrix}

AA의 determinant를 첫 번째 행을 기준으로 co-factor expansion을 이용하여 구해보면

detA=2(1)1+11142+(1)1+21112+3(1)1+31114=2detA = 2 \cdot(-1)^{1+1}\begin{vmatrix}-1 & 1 \\ 4 & 2\end{vmatrix} +(-1)^{1+2}\begin{vmatrix}1 & 1 \\ 1 & 2\end{vmatrix} + 3\cdot(-1)^{1+3} \begin{vmatrix}1 & -1 \\ 1 & 4\end{vmatrix} = 2

determinant가 0이 아니므로, AA는 invertible matrix입니다. 각 entry에 해당하는 cofactor를 구해보면

C11=3,C21=10,C31=4C12=1,C22=1,C32=1C13=5,C23=7,C33=3C_{11}=3, C_{21}=-10, C_{31}=4 \\C_{12}=1, C_{22}=1, C_{32}=-1 \\C_{13}=5, C_{23}=7, C_{33}=-3

이를 이용하여 A1A^{-1}를 구하면

A1=1detAadjA=12[3104111573]A^{-1}=\frac{1}{detA}adjA = \frac{1}{2}\begin{bmatrix}3 & -10 & 4 \\ 1 & 1& -1 \\ 5 & 7 & -3 \end{bmatrix}

가 됩니다.


지금까지 Cramer's rule에 대해서 알아보았습니다. 다음 포스트에서는 Vector space와 subspace에 대해서 알아보겠습니다. 질문이나 오류 있으면 댓글 남겨주세요! 감사합니다!


Appendix : Proof of Theorem


Theorem : Cramer's Rule

Let AA be an invertible n×nn \times n matrix. For any b\boldsymbol{b} in Rn\mathbb R ^n, the unique solution x\boldsymbol{x} of Ax=bA\boldsymbol{x}=\boldsymbol{b} has entries given by

xi=detAi(b)detA,  i=1,2,...,nx_i = \frac{detA_i(\boldsymbol{b})}{detA}, \ \ i=1, 2, ..., n
  • Proof

Let

A=[a11a12...a1na21a22...a2nan1an2ann]=[a1a2...an],  b=[b1b2bn]A = \begin{bmatrix} a_{11} & a_{12} & ... & a_{1n} \\ a_{21} & a_{22} & ... & a_{2n} \\ \vdots & \vdots & \ddots & \vdots \\ a_{n1} & a_{n2} & \cdots & a_{nn} \end{bmatrix} = \begin{bmatrix}\boldsymbol{a_1} & \boldsymbol{a_2} & ... & \boldsymbol{a_n} \end{bmatrix}, \ \ \boldsymbol{b}=\begin{bmatrix}b_1 \\ b_2 \\ \vdots \\ b_n\end{bmatrix}

Ax=bA\boldsymbol{x}=\boldsymbol{b}를 linear system으로 표현하면 다음과 같이 표현됩니다.

a11x1+a12x2++a1nxn=b11.a21x1+a22x2++a2nxn=b22.an1x1+an2x2++annxn=bnn.\begin{aligned} a_{11}x_1 + a_{12}x_2 + \cdots + a_{1n}x_n&=b_1 \cdots 1. \\ a_{21}x_1 + a_{22}x_2 + \cdots + a_{2n}x_n&=b_2 \cdots 2. \\ \vdots \\ a_{n1}x_1 + a_{n2}x_2 + \cdots + a_{nn}x_n&=b_n \cdots n. \end{aligned}

첫 번째 식부터 n 번째 식까지 오른쪽에 번호를 통하여 나타내었습니다.

여기서 ii번 째 식의 양변에 CijC_{ij}를 곱해주겠습니다. (i=1,2,...,ni=1, 2, ... ,n)

a11C1jx1+a12C1jx2++a1nC1jxn=C1jb11.a21C2jx1+a22C2jx2++a2nC2jxn=C2jb22.an1Cnjx1+an2Cnjx2++annCnjxn=Cnjbnn.\begin{aligned} a_{11}C_{1j}x_1 + a_{12}C_{1j}x_2 + \cdots + a_{1n}C_{1j}x_n&=C_{1j}b_1 \cdots 1. \\ a_{21}C_{2j}x_1 + a_{22}C_{2j}x_2 + \cdots + a_{2n}C_{2j}x_n&=C_{2j}b_2 \cdots 2. \\ \vdots \\ a_{n1}C_{nj}x_1 + a_{n2}C_{nj}x_2 + \cdots + a_{nn}C_{nj}x_n&=C_{nj}b_n \cdots n. \end{aligned}

이 후, 1번 식부터 n번 식 모두를 더하고 x1,x2,...,xnx_1, x_2, ..., x_n에 대해서 정리를 하면

p1x1+p2x2++pnxn=bp_1x_1+p_2x_2+\cdots+p_nx_n = b'

where

pi=a1iC1j+a2iC2j++aniCnj  (i=1,2,...,n)b=b1C1j+b2C2j+...bnCnjp_i = a_{1i}C_{1j}+a_{2i}C_{2j}+\cdots+a_{ni}C_{nj} \ \ (i=1, 2, ..., n) \\ b' =b_1C_{1j}+b_2C_{2j}+...b_nC_{nj}

가 됩니다.

pip_i에 대해서 살펴보겠습니다.

만약 i=ji= j라면

pi=a1iC1i+a2iC2i++aniCni=detAp_i = a_{1i}C_{1i}+a_{2i}C_{2i}+\cdots+a_{ni}C_{ni} = detA

가 됩니다. ii 번째 column을 기준으로 co-factor expansion을 한 determinant입니다.

만약 iji\neq j라면

pi=a1iC1j+a2iC2j++aniCnj=detAj(ai)p_i = a_{1i}C_{1j}+a_{2i}C_{2j}+\cdots+a_{ni}C_{nj} = detA_{j}(\boldsymbol{a_i})

가 됩니다. A의 jj 번째 column을 ai\boldsymbol{a_i}로 바꾼 matrix의 jj 번째 column을 기준으로 co-factor expansion을 한 결과가 됩니다.

여기서, jij\neq i이므로, Aj(ai)A_{j}(\boldsymbol{a_i})에서 ii 번째, jj 번째 column이 ai\boldsymbol{a_i}입니다.

Aj(ai)A_{j}(\boldsymbol{a_i})에서 똑같은 column이 존재하기 때문에, Aj(ai)A_{j}(\boldsymbol{a_i})의 determinant가 0입니다.

detAj(ai)=0detA_{j}(\boldsymbol{a_i})=0

마지막으로 b\boldsymbol{b}'을 살펴보면

b=b1C1j+b2C2j+...bnCnj=detAj(b)b' =b_1C_{1j}+b_2C_{2j}+...b_nC_{nj} = detA_j(\boldsymbol{b})

인 것을 알 수 있습니다.

이를 이용하여 식을 정리하면

p1x1+p2x2++pnxn=bpjxj=b(detA)xj=detAj(b)p_1x_1+p_2x_2+\cdots+p_nx_n = b' \\ \Rightarrow p_jx_j= b' \\ \Rightarrow (detA) x_j = detA_j(\boldsymbol{b})

이 성립하고, AA가 invertible하므로

xj=detAj(b)detAx_j = \frac{ detA_j(\boldsymbol{b})}{detA}

결과를 얻을 수 있습니다.

profile
데이터 분석가 새싹

0개의 댓글