[선형대수] Lecture 8: Solving Ax = b: row reduced form R

이재호·2025년 3월 4일

선형대수

목록 보기
8/31

https://ocw.mit.edu/courses/18-06-linear-algebra-spring-2010/video_galleries/video-lectures/

"어떻게 Ax=bAx=b에 대한 값을 찾을 수 있을까."
해당 강의는 위 내용에 대해서 다룬다.

먼저 예시와 함께 시작하자.

x1+2x2+2x3+2x4=b1x_1+ 2x_2 + 2x_3 + 2x_4 = b_1
2x1+4x2+6x3+8x4=b22x_1+ 4x_2 + 6x_3 + 8x_4 = b_2
3x1+6x2+8x3+10x4=b33x_1+ 6x_2 + 8x_3 + 10x_4 = b_3

위 방정식들을 Ax=bAx=b로 표현하면 다음과 같다.

[1222246836810][x1x2x3x4]=[b1b2b3]\begin{bmatrix} 1 & 2 & 2 & 2 \\ 2 & 4 & 6 & 8 \\ 3 & 6 & 8 & 10 \\ \end{bmatrix} \begin{bmatrix} x_1 \\ x_2 \\ x_3 \\ x_4 \end{bmatrix} = \begin{bmatrix} b_1 \\ b_2 \\ b_3 \\ \end{bmatrix}

그리고 AAbb를 합쳐 보자.

[1222b12468b236810b3]=Augmented matrix [A b]\begin{bmatrix} \begin{array}{c c c c | c} 1 & 2 & 2 & 2 & b_1 \\ 2 & 4 & 6 & 8 & b_2 \\ 3 & 6 & 8 & 10 & b_3 \\ \end{array} \end{bmatrix} = \text{Augmented matrix [A b]}

그리고 해당 행렬을 Elimination을 적용해 보자.

[1222b12468b236810b3][1222b10024b22b10024b33b1][1222b10024b22b10000b3b2b1]\begin{bmatrix} \begin{array}{c c c c | c} 1 & 2 & 2 & 2 & b_1 \\ 2 & 4 & 6 & 8 & b_2 \\ 3 & 6 & 8 & 10 & b_3 \\ \end{array} \end{bmatrix} \rightarrow \begin{bmatrix} 1 & 2 & 2 & 2 & b_1 \\ 0 & 0 & 2 & 4 & b_2-2b_1 \\ 0 & 0 & 2 & 4 & b_3-3b_1 \\ \end{bmatrix} \rightarrow \begin{bmatrix} 1 & 2 & 2 & 2 & b_1 \\ 0 & 0 & 2 & 4 & b_2-2b_1 \\ 0 & 0 & 0 & 0 & b_3-b_2-b_1 \\ \end{bmatrix}

그러면 위 행렬에서 0=b3b2b10=b_3-b_2-b_1라는 정보를 얻을 수 있다.
그리고 위 정보를 통해 임의로 b=[156]b=\begin{bmatrix}1 \\ 5\\ 6\end{bmatrix}을 대입해보자. 그리고 위 행렬에 대입해보자.

[122210024300000]\begin{bmatrix} \begin{array}{c c c c | c} 1 & 2 & 2 & 2 & 1 \\ 0 & 0 & 2 & 4 & 3 \\ 0 & 0 & 0 & 0 & 0 \\ \end{array} \end{bmatrix}

문제없이 잘 나오는 것을 확인할 수 있다.


그렇다면 "Ax=bAx=b"를 풀이할 수 있는 조건은 무엇일까? 강의에서 제시하는 내용은 다음과 같다.

  • Ax=bAx=b is solvable when bb is in C(A)C(A)(=AA의 칼럼 공간)
  • If a comb. of rows of AA gives zero row, then the same comb. of entries of BB must give 00. (ex. 위에서 [0 0 0 0  0][0 \ 0 \ 0 \ 0 \ | \ 0]과 같은 경우)
[122210024300000]\begin{bmatrix} \begin{array}{c c c c | c} 1 & 2 & 2 & 2 & 1 \\ 0 & 0 & 2 & 4 & 3 \\ 0 & 0 & 0 & 0 & 0 \\ \end{array} \end{bmatrix}

이어서 이를 일반화하면 다음과 같다. (위 행렬을 예시로 설명한다)

  1. xparticularx_{particular}
    • set all free variables to zero. (x2=0,x4=0x_2=0, x_4=0)
    • sovle Axp=bAx_p=b for the pivot variables.
      x1+2x3=12x3=3x3=32,x1=2x_1+2x_3=1 \\ 2x_3=3 \\ \rightarrow x_3=\frac{3}{2}, x_1=-2
    • 따라서 위 예시에서 xpx_pxp=[20230]x_p=\begin{bmatrix}-2 \\ 0 \\ \frac{2}{3} \\ 0 \end{bmatrix}이 된다.
  2. xnullspacex_{nullspace}
    • free variable 중 하나는 1로, 나머진 0으로. (x2=1,x4=0x_2=1, x_4=0 or x2=0,x4=1x_2=0, x_4=1)
    • solve Axn=0Ax_n=0.
    • xn=c1[2100]+c2[2021]x_n=c_1\begin{bmatrix}-2 \\ 1 \\ 0 \\ 0\end{bmatrix}+c_2\begin{bmatrix}2 \\ 0 \\ -2 \\ 1\end{bmatrix}
  3. xx(정답 벡터)
    • x=xp+xnx=x_p+x_n
    • xnsubspace of xpx_n \in \text{subspace of $x_p $}
    • A(xp+xn)=bA(x_p + x_n) = b

이어서 rank(pivot 칼럼의 수)에 대해서 좀더 알아보자.

  • Am×nA_{m\times n} 행렬의 rank rr은 항상 다음과 같은 조건을 만족한다.
    • rm,rnr \le m, r \le n : pivot의 개수는 항상 행 혹은 열의 개수보다 많을 수 없다.

그렇다면 다음 상황들에 대해서 생각해보자.

  1. r=n<mr=n<m

    • ex. A=[13216151][10010000]A=\begin{bmatrix}1 & 3 \\ 2 & 1 \\ 6 & 1 \\ 5 & 1\end{bmatrix} \rightarrow \begin{bmatrix}1 & 0 \\ 0 & 1 \\ 0 & 0 \\ 0 & 0\end{bmatrix}
    • no free variables.
    • N(A)=0N(A)=0 only.
    • x=xpx=x_p (unique solution if it exists)
    • R=[I0]R=\begin{bmatrix}I \\ 0\end{bmatrix}
    • Ax=b0 or 1 solution.Ax=b \rightarrow \text{0 or 1 solution.}
  2. r=m<nr=m<n

    • ex. A=[12653111][10f1f201f3f4]A=\begin{bmatrix}1 & 2 & 6 & 5 \\ 3 & 1 & 1 & 1\end{bmatrix} \rightarrow \begin{bmatrix}1 & 0 & f_1 & f_2 \\ 0 & 1 & f_3 & f_4\end{bmatrix}
    • nr(=nm)n-r(=n-m) free variables.
    • R=[IF]R=\begin{bmatrix}I & F\end{bmatrix}
    • Ax=b solutions.Ax=b \rightarrow \text{$\infty$ solutions.}
  3. r=m=nr=m=n

    • ex. A=[1231][1001]A=\begin{bmatrix}1 & 2 \\ 3 & 1\end{bmatrix} \rightarrow \begin{bmatrix}1 & 0 \\ 0 & 1\end{bmatrix}
    • R=IR=I
    • AA is a "invertible matrix".
    • Ax=b1 solution only.Ax=b \rightarrow \text{1 solution only.}
  4. r<m,r<nr<m,r<n

    • R=[IF00]R=\begin{bmatrix}I & F \\ 0 & 0\end{bmatrix}
    • Ax=b0 or  solution(s).Ax=b \rightarrow \text{0 or $\infty$ solution(s).}
profile
천천히, 그리고 꾸준히.

0개의 댓글