Ax=b
Ax=b의 의미
m×n 행렬 A에 대해서 Ax는 아래와 같이 표현할 수 있다.
Ax=[a1 a2 ⋯ an]⎣⎢⎢⎢⎢⎡x1x2⋮xn⎦⎥⎥⎥⎥⎤=x1a1+x2a2+⋯+xnan
- n개의 열벡터(ai)로 이루어진 행렬 A와 n개의 미지수로 이루어진 벡터 x의 곱
- 각 열벡터(ai)와 벡터 x의 요소들(xi)의 선형결합 → b
= A의 column들의 liner combination이다.
따라서 Ax=b 방정식은 벡터 b가 A의 Column space 안에 들어있다면 해를 가진다. (if and only if)
x is exist, if and only if b∈C(A)
Some Terminology : 용어
elimination이 진행된 행렬 A에 대해, Ax=0 에서
- Pivot column: 행렬 A에서 pivot을 가지는 column
- Free column: 행렬 A에서 pivot을 가지지 않는 column
- Pivot variable: 벡터 x에서 행렬 A의 pivot column에 곱해지는 변수 (요소)
- Free variable: 벡터 x에서 행렬 A의 free column에 곱해지는 변수 (요소)
예시 : 아래 행렬에서 pivot/free column과 pivot/free variable은?
⎣⎢⎡1002101400−10⎦⎥⎤⎣⎢⎢⎢⎡x1x2x3x4⎦⎥⎥⎥⎤=⎣⎢⎡000⎦⎥⎤
- pivot column : [1 0 0]T, [2 1 0]T (pivot : 1,1)
- free column : [1 4 0]T, [0 −1 0]T
- pivot variable : x1,x2
- free variable : x3,x4
NullSpace 구하기 : special solution
method 1
각각의 pivot variable을 free variable를 이용해서 표현하고 벡터 x에 대입한다.
-
Ax=0 꼴로 행렬을 표현한다.
-
A에 대하여 Gauss Elimination을 진행한다.
-
이후 x의 pivot variable을 free variable로 표현한다.
→ fv의 선형결합으로 표현된 pc를 각각의 표현으로 분리한다.
-
그렇게 free variable로 표현된 변수들을 벡터 x에 대입한다.
-
벡터 x를 각각의 free variable과 벡터*의 linear combination 으로 표현한다. ⇒ solution to Ax=0, N(A)
이때, 각각의 벡터를 special solution이라고 부른다.
- 표현 예시
-
variable에 대해 다음 두 방정식이 주어졌다. (fv : x2,x4,x5 일 때)
- x1=−2x2+2x5
- x3=−x4−2x5
-
식을 활용하여 벡터 x에 대입한다.
x=⎣⎢⎢⎢⎢⎢⎡x1x2x3x4x5⎦⎥⎥⎥⎥⎥⎤=⎣⎢⎢⎢⎢⎢⎡−2x2+2x5x2−x4−2x5x4x5⎦⎥⎥⎥⎥⎥⎤
-
각각의 free variable에 대한 선형결합으로 표현한다.
x=x2⎣⎢⎢⎢⎢⎢⎡−21000⎦⎥⎥⎥⎥⎥⎤+x4⎣⎢⎢⎢⎢⎢⎡00−110⎦⎥⎥⎥⎥⎥⎤+x5⎣⎢⎢⎢⎢⎢⎡20−201⎦⎥⎥⎥⎥⎥⎤, x2,x4,x5∈R
method 2
free variable중 하나의 변수에는 1을 대입하고 나머지에는 0을 대입하는 과정을 반복한다.
-
Ax=0 꼴로 행렬을 표현한다.
-
A에 대하여 Gauss Elimination을 진행한다.
-
여러개의 free variable중 하나를 1로 하고 다른 free variale들은 0을 대입하여 표현한다.
-
3번 과정을 거쳐 이제 pivot variable의 값만 알아내면 되는 벡터 x에 대해 Ax=0 을 풀어 벡터: special solution를 구한다.
-
각각의 모든 free variable에 해당 과정을 반복해서(1인 variable을 변경해줌) 적용하여 모든 special solution들을 구한다.
-
이렇게 구한 special colution들의 linear combination ⇒ solution to Ax=0, N(A)
- 표현 예시
fv : x2,x4,x5 일 때
- variable에 대해 다음 두 방정식이 주어졌다. (fv : x2,x4,x5 일 때)
- x1=−2x2+2x5
- x3=−x4−2x5
- 각 case별 Special solution
- x2=1,x4=0,x5=0
x=⎣⎢⎢⎢⎢⎢⎡−2x2+2x51−x4−2x500⎦⎥⎥⎥⎥⎥⎤=⎣⎢⎢⎢⎢⎢⎡−21000⎦⎥⎥⎥⎥⎥⎤
- x2=0,x4=1,x5=0
x=⎣⎢⎢⎢⎢⎢⎡−2x2+2x50−x4−2x510⎦⎥⎥⎥⎥⎥⎤=⎣⎢⎢⎢⎢⎢⎡00−110⎦⎥⎥⎥⎥⎥⎤
- x2=0,x4=0,x5=1
x=⎣⎢⎢⎢⎢⎢⎡−2x2+2x50−x4−2x501⎦⎥⎥⎥⎥⎥⎤=⎣⎢⎢⎢⎢⎢⎡20−201⎦⎥⎥⎥⎥⎥⎤
- Special solution들의 linear combination
x=x2⎣⎢⎢⎢⎢⎢⎡−21000⎦⎥⎥⎥⎥⎥⎤+x4⎣⎢⎢⎢⎢⎢⎡00−110⎦⎥⎥⎥⎥⎥⎤+x5⎣⎢⎢⎢⎢⎢⎡20−201⎦⎥⎥⎥⎥⎥⎤, x2,x4,x5∈R
special solution?
- NullSpace를 이루는 vector
- special solution은 free variable(=free column)의 개수와 동일하다.
- free variable의 개수 =
전체 column의 갯수 - pivot의 개수
(n−r)
Complete Solution
Theorem
벡터 w가 Ax=b의 임의의 solution이라고 가정하자. Then,
- Ax=b 의 solution이 y이면 y=w+z 이다. (z∈N(A))
pf)
- if Ay=b 이라면,
Ay−Aw=b−b=0 A(y−w)=0
- 따라서 벡터 y−w는 NullSpace에 속하는 vector이다.
y−w∈N(A)
- 벡터 z를 NullSpace에 속하는 또다른 vector라고 하면, linear combination 성질에 의해
y−w=z (z∈N(A)) y=w+z
- ↔ y=w+z (z∈N(A)) 이면 y는 Ax=b의 solution이다.
pf)
- if y=w+z 이라면,
Ay=A(w+z)=Aw+Az
- 이때, 처음 Theorem의 가정에 따라 Aw=b이고 z는 NullSpace에 속하는 벡터이므로 Az=0
b+0=b
- 따라서 y는 Ax=b의 solution이다.
Ay=b
y=w+z 를 이용하여 모든 solution을 표현할 수 있다.
→ 즉, NullSpace를 구하여 그 벡터 (z∈N(A))에 임의의 solution w를 더하면 complete solution을 표현할 수 있다.
기하학적 의미
- 임의의 solution과 NullSpace의 vector(special solution)을 더한 vector의 종점은 여전히 Ax=b 안에 들어있다.
- 따라서 NullSpace의 어떤 solution과 임의의 solution의 합으로 Ax=b의 모든 solution을 표현할 수 있다.
과정
-
Ax=b의 아무 solution을 찾는다. (=xp)
- 임의의 어떤 solution이던간에 관계없음
- 따라서 가장 쉽게 solution을 찾을 수 있도록 fv를 0으로 가정하고 Ax=b를 풀어서 particular solution을 구한다.
-
Ax=0의 solution을 찾는다. → Null Space를 찾는다. (N(A)=xn)
-
xp+xn이 Ax=b의 complete solution이 된다.
(NullSpace의 모든 vector에 대해 진행하여 새로운 벡터 공간을 정의함)
완전해(complete solution)은 특수해들의 선형조합이다.
- 특수헤?
- Ax=0 의 해 (=Special solution) : NullSpace의 벡터
- Ax=b의 임의의 solution (=particular solution)