이번 포스트에서는 Cramer's rule에 대해서 알아보겠습니다.
1) Cramer's Rule
Cramer's rule을 정의하기 위해서는 하나의 notation 정의가 필요합니다.
for any n×n matrix A and any b in Rn, let Ai(b) be the matrix obtained from A by replaceing column i by the vector b
Ai(b)=[a1a2...ai−1bai+1...an]
즉 Ai(b) 는 matrix A의 i 번째 column 대신 b를 넣은 새로운 matrix입니다.
Let A be an invertible n×n matrix. For any b in Rn, the unique solution x of Ax=b has entries given by
xi=detAdetAi(b), i=1,2,...,n
Cramer's rule을 이용하면, determinant를 이용하여 linear system의 solution을 구할 수 있습니다.
(증명은 appendix 참고)
example
3x1−2x2−5x1+4x2=6=8
위 linear system을 matrix equation으로 바꾸면
[3−5−24][x1x2]=[68]
Cramer's rule을 적용하기 위해 detAi(b)와 detA를 게산하면
detA1(b)=∣∣∣∣∣68−24∣∣∣∣∣=40 detA2(b)=∣∣∣∣∣3−568∣∣∣∣∣=54 detA=∣∣∣∣∣3−5−24∣∣∣∣∣=2
Cramer's rule을 적용하면
x1=detAdetA1(b)=20 x2=detAdetA2(b)=27
따라서 위 linear system의 solution은
x=[2027]
입니다.
Cramer's rule을 적용하여 A−1를 찾을 수 있습니다.
1) FInding A−1
Let A be an invertible n×n matrix. Then, the jth column of A−1 is a vector x that satisfies
Ax=ej
where ej is the jth column of identity matrix.
Inverse의 정의에 의해
를 만족합니다. 따라서 위 matrix 식의 결과를 각 column별로 살펴보면
Ax=ej
을 만족하는 x가 A−1의 jth column이 되는 것을 알 수 있습니다.
위 방정식의 solution을 구할 때, Cramer's rule을 이용하면
xij=detAdetAi(ej)
가 되고, xij는 A−1의 (i,j) entry가 됩니다.
Ai(ej)를 살펴보면
Ai(ej)=[a1...ai−1ejai+1...an]=⎣⎢⎢⎢⎢⎢⎢⎢⎢⎢⎡a1...ai−100⋮1⋮0ai+1...an⎦⎥⎥⎥⎥⎥⎥⎥⎥⎥⎤
입니다. 따라서 detAi(ej)를 co-factor expansion을 이용하여 구할 때 i 번째 column을 기준으로 구하게 됩니다. Ai(ej)의 (j,i) 위치에서의 co-factor는 Ai(ej)에서 j 번째 row, i 번째 column을 제외한 matrix의 determinant입니다.
해당 determinant는 A에서 j 번째 row, i 번째 column을 제외한 matrix의 determinant와 같습니다. 즉 Ai(ej)의 (i,j) 위치에 해당하는 cofactor는 A의 (i,j) 위치에 해당하는 cofactor와 일치합니다. 따라서
detAi(ej)=Cji
입니다.
A−1의 (i,j) entry에 해당하는 값을 알았으니, A−1를 표현하면
A−1=detA1⎣⎢⎢⎢⎢⎡C11C12⋮C1nC21C22⋮C2n......⋱...Cn1Cn2⋮Cnn⎦⎥⎥⎥⎥⎤
이 됩니다. 여기서 detA1를 제외한 matrix 부분
⎣⎢⎢⎢⎢⎡C11C12⋮C1nC21C22⋮C2n......⋱...Cn1Cn2⋮Cnn⎦⎥⎥⎥⎥⎤
를 adjugate of A라고 하고, adjA로 표시합니다.
A−1=detA1adjA
example
A=⎣⎢⎡2111−14312⎦⎥⎤
A의 determinant를 첫 번째 행을 기준으로 co-factor expansion을 이용하여 구해보면
detA=2⋅(−1)1+1∣∣∣∣∣−1412∣∣∣∣∣+(−1)1+2∣∣∣∣∣1112∣∣∣∣∣+3⋅(−1)1+3∣∣∣∣∣11−14∣∣∣∣∣=2
determinant가 0이 아니므로, A는 invertible matrix입니다. 각 entry에 해당하는 cofactor를 구해보면
C11=3,C21=−10,C31=4C12=1,C22=1,C32=−1C13=5,C23=7,C33=−3
이를 이용하여 A−1를 구하면
A−1=detA1adjA=21⎣⎢⎡315−10174−1−3⎦⎥⎤
가 됩니다.
지금까지 Cramer's rule에 대해서 알아보았습니다. 다음 포스트에서는 Vector space와 subspace에 대해서 알아보겠습니다. 질문이나 오류 있으면 댓글 남겨주세요! 감사합니다!
Appendix : Proof of Theorem
Theorem : Cramer's Rule
Let A be an invertible n×n matrix. For any b in Rn, the unique solution x of Ax=b has entries given by
xi=detAdetAi(b), i=1,2,...,n
Let
A=⎣⎢⎢⎢⎢⎡a11a21⋮an1a12a22⋮an2......⋱⋯a1na2n⋮ann⎦⎥⎥⎥⎥⎤=[a1a2...an], b=⎣⎢⎢⎢⎢⎡b1b2⋮bn⎦⎥⎥⎥⎥⎤
Ax=b를 linear system으로 표현하면 다음과 같이 표현됩니다.
a11x1+a12x2+⋯+a1nxna21x1+a22x2+⋯+a2nxn⋮an1x1+an2x2+⋯+annxn=b1⋯1.=b2⋯2.=bn⋯n.
첫 번째 식부터 n 번째 식까지 오른쪽에 번호를 통하여 나타내었습니다.
여기서 i번 째 식의 양변에 Cij를 곱해주겠습니다. (i=1,2,...,n)
a11C1jx1+a12C1jx2+⋯+a1nC1jxna21C2jx1+a22C2jx2+⋯+a2nC2jxn⋮an1Cnjx1+an2Cnjx2+⋯+annCnjxn=C1jb1⋯1.=C2jb2⋯2.=Cnjbn⋯n.
이 후, 1번 식부터 n번 식 모두를 더하고 x1,x2,...,xn에 대해서 정리를 하면
p1x1+p2x2+⋯+pnxn=b′
where
pi=a1iC1j+a2iC2j+⋯+aniCnj (i=1,2,...,n)b′=b1C1j+b2C2j+...bnCnj
가 됩니다.
pi에 대해서 살펴보겠습니다.
만약 i=j라면
pi=a1iC1i+a2iC2i+⋯+aniCni=detA
가 됩니다. i 번째 column을 기준으로 co-factor expansion을 한 determinant입니다.
만약 i=j라면
pi=a1iC1j+a2iC2j+⋯+aniCnj=detAj(ai)
가 됩니다. A의 j 번째 column을 ai로 바꾼 matrix의 j 번째 column을 기준으로 co-factor expansion을 한 결과가 됩니다.
여기서, j=i이므로, Aj(ai)에서 i 번째, j 번째 column이 ai입니다.
Aj(ai)에서 똑같은 column이 존재하기 때문에, Aj(ai)의 determinant가 0입니다.
detAj(ai)=0
마지막으로 b′을 살펴보면
b′=b1C1j+b2C2j+...bnCnj=detAj(b)
인 것을 알 수 있습니다.
이를 이용하여 식을 정리하면
p1x1+p2x2+⋯+pnxn=b′⇒pjxj=b′⇒(detA)xj=detAj(b)
이 성립하고, A가 invertible하므로
xj=detAdetAj(b)
결과를 얻을 수 있습니다.