[선형대수학] Solvability of Ax = b

Vaughan·2022년 8월 10일
0

선형대수학

목록 보기
8/14
post-thumbnail

Ax=bA \bold x=\bold b

Ax=bA\bold x= \bold b의 의미

m×nm\times n 행렬 AA에 대해서 AxA\bold x는 아래와 같이 표현할 수 있다.

Ax=[a1  a2    an][x1x2xn]=x1a1+x2a2++xnanA\bold x = \begin{bmatrix}\bold{a_1}\ \ \bold{a_2}\ \ \cdots \ \ \bold{a_n}\end{bmatrix} \begin{bmatrix} x_1 \\ x_2 \\ \vdots \\ x_n \end{bmatrix} = x_1\bold{a_1} + x_2\bold{a_2} + \cdots + x_n\bold{a_n}
  • nn개의 열벡터(ai\bold{a_i})로 이루어진 행렬 AAnn개의 미지수로 이루어진 벡터 x\bold x의 곱
  • 각 열벡터(ai\bold{a_i})와 벡터 x\bold x의 요소들(xix_i)의 선형결합b\bold b
    = AA의 column들의 liner combination이다.

따라서 Ax=bA\bold x=\bold b 방정식은 벡터 b\bold bAA의 Column space 안에 들어있다면 해를 가진다. (if and only if)

x is exist,  if and only if bC(A)\bold x \text{ is exist, }\ \text{if and only if }\bold b ∈ C(A)

Some Terminology : 용어

elimination이 진행된 행렬 AA에 대해, Ax=0A\bold x = 0 에서

  • Pivot column: 행렬 AA에서 pivot을 가지는 column
  • Free column: 행렬 AA에서 pivot을 가지지 않는 column
  • Pivot variable: 벡터 x\bold x에서 행렬 AA의 pivot column에 곱해지는 변수 (요소)
  • Free variable: 벡터 x\bold x에서 행렬 AA의 free column에 곱해지는 변수 (요소)

예시 : 아래 행렬에서 pivot/free column과 pivot/free variable은?

[121001410000][x1x2x3x4]=[000]\begin{bmatrix} 1 & 2 & 1 & 0 \\ 0 & 1 & 4 & -1 \\ 0 & 0 & 0 & 0 \end{bmatrix} \begin{bmatrix} x_1 \\ x_2 \\ x_3 \\ x_4 \end{bmatrix} = \begin{bmatrix} 0 \\ 0 \\ 0\\ \end{bmatrix}
  • pivot column : [1  0  0]T[1\ \ 0 \ \ 0]^T, [2  1  0]T[2\ \ 1 \ \ 0]^T (pivot : 1,11, 1)
  • free column : [1  4  0]T[1\ \ 4 \ \ 0]^T, [0  1  0]T[0\ \ -1 \ \ 0]^T
  • pivot variable : x1,x2x_1, x_2
  • free variable : x3,x4x_3, x_4

NullSpace 구하기 : special solution

method 1

각각의 pivot variable을 free variable를 이용해서 표현하고 벡터 x\bold x에 대입한다.

  1. Ax=0A\bold x = \bold 0 꼴로 행렬을 표현한다.

  2. AA에 대하여 Gauss Elimination을 진행한다.

  3. 이후 x\bold x의 pivot variable을 free variable로 표현한다.

    → fv의 선형결합으로 표현된 pc를 각각의 표현으로 분리한다.

  4. 그렇게 free variable로 표현된 변수들을 벡터 x\bold x에 대입한다.

  5. 벡터 x\bold x를 각각의 free variable과 벡터*의 linear combination 으로 표현한다. ⇒ solution to Ax=0A\bold x= \bold 0, N(A)N(A)

    이때, 각각의 벡터를 special solution이라고 부른다.


  • 표현 예시
    1. variable에 대해 다음 두 방정식이 주어졌다. (fv : x2,x4,x5x_2, x_4, x_5 일 때)

      • x1=2x2+2x5x_1 = -2x_2 + 2x_5
      • x3=x42x5x_3 = -x_4-2x_5
    2. 식을 활용하여 벡터 x\bold x에 대입한다.

      x=[x1x2x3x4x5]=[2x2+2x5x2x42x5x4x5]\bold x = \begin{bmatrix} x_1\\x_2\\x_3\\x_4\\x_5 \end{bmatrix} = \begin{bmatrix} -2x_2 + 2x_5\\x_2\\-x_4-2x_5\\x_4\\x_5 \end{bmatrix}
    3. 각각의 free variable에 대한 선형결합으로 표현한다.

      x=x2[21000]+x4[00110]+x5[20201],   x2,x4,x5R\bold x = x_2 \begin{bmatrix} -2\\1\\0\\0\\0 \end{bmatrix}+x_4 \begin{bmatrix} 0\\0\\-1\\1\\0 \end{bmatrix}+x_5 \begin{bmatrix} 2\\0\\-2\\0\\1 \end{bmatrix},\ \ \ x_2,x_4,x_5\in \mathbb R

method 2

free variable중 하나의 변수에는 1을 대입하고 나머지에는 0을 대입하는 과정을 반복한다.

  1. Ax=0A\bold x = \bold 0 꼴로 행렬을 표현한다.

  2. AA에 대하여 Gauss Elimination을 진행한다.

  3. 여러개의 free variable중 하나를 1로 하고 다른 free variale들은 0을 대입하여 표현한다.

  4. 3번 과정을 거쳐 이제 pivot variable의 값만 알아내면 되는 벡터 x\bold x에 대해 Ax=0A\bold x = \bold 0 을 풀어 벡터: special solution를 구한다.

  5. 각각의 모든 free variable에 해당 과정을 반복해서(1인 variable을 변경해줌) 적용하여 모든 special solution들을 구한다.

  6. 이렇게 구한 special colution들의 linear combination ⇒ solution to Ax=0A\bold x= \bold 0, N(A)N(A)


  • 표현 예시
    fv : x2,x4,x5x_2, x_4, x_5 일 때
    1. variable에 대해 다음 두 방정식이 주어졌다. (fv : x2,x4,x5x_2, x_4, x_5 일 때)
      • x1=2x2+2x5x_1 = -2x_2 + 2x_5
      • x3=x42x5x_3 = -x_4-2x_5
    2. 각 case별 Special solution
      • x2=1,x4=0,x5=0x_2 = 1, x_4 = 0, x_5 = 0
        x=[2x2+2x51x42x500]=[21000]\bold x = \begin{bmatrix} -2x_2 + 2x_5\\1\\-x_4-2x_5\\0\\0 \end{bmatrix} = \begin{bmatrix} -2 \\1\\0\\0\\0 \end{bmatrix}
      • x2=0,x4=1,x5=0x_2 = 0, x_4 = 1, x_5 = 0
        x=[2x2+2x50x42x510]=[00110]\bold x = \begin{bmatrix} -2x_2 + 2x_5\\0\\-x_4-2x_5\\1\\0 \end{bmatrix} = \begin{bmatrix} 0 \\0\\-1\\1\\0 \end{bmatrix}
      • x2=0,x4=0,x5=1x_2 = 0, x_4 = 0, x_5 = 1
        x=[2x2+2x50x42x501]=[20201]\bold x = \begin{bmatrix} -2x_2 + 2x_5\\0\\-x_4-2x_5\\0\\1 \end{bmatrix} = \begin{bmatrix} 2 \\0\\-2\\0\\1 \end{bmatrix}
    3. Special solution들의 linear combination
      x=x2[21000]+x4[00110]+x5[20201],   x2,x4,x5R\bold x = x_2 \begin{bmatrix} -2\\1\\0\\0\\0 \end{bmatrix}+x_4 \begin{bmatrix} 0\\0\\-1\\1\\0 \end{bmatrix}+x_5 \begin{bmatrix} 2\\0\\-2\\0\\1 \end{bmatrix},\ \ \ x_2,x_4,x_5\in \mathbb R

special solution?

  • NullSpace를 이루는 vector
  • special solution은 free variable(=free column)의 개수와 동일하다.
  • free variable의 개수 = 전체 column의 갯수 - pivot의 개수 (nrn-r)

Complete Solution

Theorem

벡터 w\bold wAx=bA\bold x= \bold b의 임의의 solution이라고 가정하자. Then,

  • Ax=bA\bold x= \bold b 의 solution이 y\bold y이면 y=w+z\bold y = \bold w+\bold z 이다. (zN(A)\bold z \in N(A))
    pf)
    1. if  Ay=b\text{if }\ A\bold y = \bold b 이라면,
      AyAw=bb=0 A(yw)=0A\bold y - A\bold w = \bold b - \bold b =0\\\ \\A(\bold y - \bold w) = 0
    2. 따라서 벡터 yw\bold y - \bold w는 NullSpace에 속하는 vector이다.
      ywN(A)\bold y - \bold w \in N(A)
    3. 벡터 z\bold z를 NullSpace에 속하는 또다른 vector라고 하면, linear combination 성질에 의해
      yw=z  (zN(A)) y=w+z\bold y-\bold w = \bold z\ \ (\bold z \in N(A))\\\ \\ \bold y = \bold w + \bold z
  • ↔  y=w+z\bold y = \bold w+\bold z (zN(A)\bold z \in N(A)) 이면 y\bold yAx=bA\bold x= \bold b의 solution이다.
    pf)
    1. if  y=w+z\text{if }\ \bold y = \bold w+\bold z 이라면,
      Ay=A(w+z)=Aw+AzA\bold y = A(\bold w + \bold z)= A\bold w + A\bold z
    2. 이때, 처음 Theorem의 가정에 따라 Aw=bA\bold w = \bold b이고 z\bold z는 NullSpace에 속하는 벡터이므로 Az=0A\bold z = \bold 0
      b+0=b\bold b + \bold 0 = \bold b
    3. 따라서 y\bold yAx=bA\bold x= \bold b의 solution이다.
      Ay=bA\bold y = \bold b

y=w+z\bold y = \bold w + \bold z 를 이용하여 모든 solution을 표현할 수 있다.

→ 즉, NullSpace를 구하여 그 벡터 (zN(A)\bold z \in N(A))에 임의의 solution w\bold w를 더하면 complete solution을 표현할 수 있다.


기하학적 의미

  • 임의의 solution과 NullSpace의 vector(special solution)을 더한 vector의 종점은 여전히 Ax=bA\bold x = \bold b 안에 들어있다.
  • 따라서 NullSpace의 어떤 solution과 임의의 solution의 합으로 Ax=bA\bold x = \bold b의 모든 solution을 표현할 수 있다.


과정

  1. Ax=bA\bold x = \bold b의 아무 solution을 찾는다. (=xp)(=x_p)

    • 임의의 어떤 solution이던간에 관계없음
    • 따라서 가장 쉽게 solution을 찾을 수 있도록 fv를 0으로 가정하고 Ax=bA\bold x = \bold b를 풀어서 particular solution을 구한다.
  2. Ax=0A\bold x = \bold 0의 solution을 찾는다. → Null Space를 찾는다. (N(A)=xn)(N(A) = x_n)

  3. xp+xnx_p + x_nAx=bA\bold x = \bold b의 complete solution이 된다.

    (NullSpace의 모든 vector에 대해 진행하여 새로운 벡터 공간을 정의함)

완전해(complete solution)은 특수해들의 선형조합이다.

  • 특수헤?
    • Ax=0A\bold x = \bold 0 의 해 (=Special solution) : NullSpace의 벡터
    • Ax=bA\bold x = \bold b의 임의의 solution (=particular solution)
profile
우주의 아름다움도 다양한 지식을 접하며 스스로의 생각이 짜여나갈 때 불현듯 나를 덮쳐오리라.

0개의 댓글