์ด ๊ณ์ฐ๊ธฐ์ ๊ธฐ๋ฅ
์ด ๋๊ตฌ๋ ๋ฏธ์ง์๊ฐ n๊ฐ์ธ n์ ์ฐ๋ฆฝ์ผ์ฐจ๋ฐฉ์ ์์ ํ๋๋ค. ์ด๋ฅผ ๊ฐ๋จํ \(\mathbf{A}\,\mathbf{x} = \mathbf{b}\)๋ก ํ๊ธฐํ๋๋ฐ, ์ฌ๊ธฐ์ A๋ nรn ๊ณ์ ํ๋ ฌ, x๋ ๋ฏธ์ง์ ๋ฒกํฐ, b๋ ์์ ๋ฒกํฐ์ ๋๋ค. ๊ณ์ฐ๊ธฐ๋ ์ ์ผํด ๋ฒกํฐ x์ ํจ๊ป ํ๋ ฌ A์ ํ๋ ฌ์์ ํจ๊ป ๋ฐํํฉ๋๋ค. ์ฌ์ฉํ๋ ๋ฐฉ๋ฒ์ ์์ํ ์ ํ๋์์ด๋ฏ๋ก ์ด๋ ๋๋ผ์์๋ ๋์ผํ๊ฒ ์๋ํ๋ฉฐ, ๊ตญ๊ฐ๋ ๋จ์์ ๋ฐ๋ฅธ ์ฐจ์ด๊ฐ ์ ํ ์์ต๋๋ค. ์ ๋ ฅํ๋ ๋ชจ๋ ๊ฐ์ ๊ทธ์ ์ค์์ผ ๋ฟ์ ๋๋ค.
์ฌ์ฉ ๋ฐฉ๋ฒ
๋จผ์ n(๋ฐฉ์ ์๊ณผ ๋ฏธ์ง์์ ๊ฐ์)์ ์ ํฉ๋๋ค. ๊ณ์ ํ๋ ฌ A๋ฅผ ํ ์ค์ ํ ํ์ฉ ์ ๋ ฅํ๋ ์ซ์๋ ๊ณต๋ฐฑ์ด๋ ์ผํ๋ก ๊ตฌ๋ถํ๊ณ , ์ด์ด์ ์์ ๋ฒกํฐ b๋ฅผ ๊ธธ์ด n์ ๋ชฉ๋ก์ผ๋ก ์ ๋ ฅํฉ๋๋ค. ํ์ํ ์ ๋ฐ๋๋ฅผ ์ ํํ ๋ค ํ์ด๋ฅผ ์คํํ๋ฉด ๋ฉ๋๋ค. ์์, ์์, ์์ ํํ์ ๊ฐ ๋ชจ๋ ์ ๋ ฅํ ์ ์์ต๋๋ค. A์ ํ๊ณผ ์ด ๊ฐ์๊ฐ n๊ณผ ๊ฐ๊ณ b์ ๊ธธ์ด๋ n๊ณผ ์ผ์นํ๋ฉด ํด๊ฐ ๊ณ์ฐ๋๋ฉฐ, ๊ทธ๋ ์ง ์์ผ๋ฉด ์ฐจ์์ด ๋ง์ง ์๋๋ค๋ ์๋ด๊ฐ ํ์๋ฉ๋๋ค.
๊ณ์ฐ ์๋ฆฌ
์ด ๊ณ์ฐ๊ธฐ๋ ๋ถ๋ถ ํผ๋ฒํ (partial pivoting)์ ์ ์ฉํ ๊ฐ์ฐ์ค ์๊ฑฐ๋ฒ์ ์ํํฉ๋๋ค. ์ด๋ ์ํ์ ์ผ๋ก LU ๋ถํด \(\mathbf{P}\mathbf{A} = \mathbf{L}\mathbf{U}\)์ ๋์ผํฉ๋๋ค. ๊ฐ ์ด๋ง๋ค ์ฌ์ฉํ ์ ์๋ ํผ๋ฒ ์ค ์ ๋๊ฐ์ด ๊ฐ์ฅ ํฐ ๊ฐ์ ์ ํํด ์์น์ ์์ ์ฑ์ ํ๋ณดํ๊ณ , ํผ๋ฒ ์๋์ ์ฑ๋ถ์ ์๊ฑฐํ ๋ค, ๋ง์ง๋ง ๋ฏธ์ง์๋ถํฐ ์๋ก ์ฌ๋ผ๊ฐ๋ฉฐ ํ์ง ๋์ ์ ์งํํฉ๋๋ค. ๋ง์ฝ ํผ๋ฒ์ด ์ฌ์ค์ 0์ด๋ผ๋ฉด ํ๋ ฌ์์ด 0์ด ๋์ด ํ๋ ฌ์ด ํน์ด ํ๋ ฌ(singular)์ด ๋๊ณ , ๋ฐ๋ผ์ ์ ์ผํด๊ฐ ์กด์ฌํ์ง ์์ต๋๋ค. ์ด ๊ฒฝ์ฐ ๊ณ์ฐ๊ธฐ๋ 0์ผ๋ก ๋๋๋ ๋์ ์ด๋ฅผ ์๋ ค์ค๋๋ค.
$$\begin{gathered} \mathbf{A}\,\mathbf{x} = \mathbf{b}, \qquad \mathbf{A} = \mathbf{L}\mathbf{U} \\[1.5em] \text{solve}\quad \left\{ \begin{aligned} \mathbf{L}\mathbf{y} &= \mathbf{b} \quad(\text{forward}) \\ \mathbf{U}\mathbf{x} &= \mathbf{y} \quad(\text{back substitution}) \end{aligned} \right. \end{gathered}$$ $$\det(\mathbf{A}) = (-1)^{s}\prod_{k=1}^{n} u_{kk}$$
์์ ํ์ด
๋ค์ ์์ ์ดํด๋ด ์๋ค. \(2x + y - z = 8\), \(-3x - y + 2z = -11\), \(-2x + y + 2z = -3\). ์ฆ \(\mathbf{A} = [[2,1,-1],[-3,-1,2],[-2,1,2]]\), \(\mathbf{b} = [8,-11,-3]\)์ ๋๋ค. ์๊ฑฐ ๊ณผ์ ์ ๊ฑฐ์น๋ฉด \(x = 2\), \(y = 3\), \(z = -1\)์ด ๋์ต๋๋ค. ์ฒซ ๋ฒ์งธ ์์ ๋์ ํด ํ์ธํ๋ฉด $$2(2) + 3 - (-1) = 4 + 3 + 1 = 8$$๋ก ์ ํํ ๋ง์ต๋๋ค.
์์ฃผ ๋ฌป๋ ์ง๋ฌธ
ํ๋ ฌ์์ด 0์ด๋ฉด ์ด๋ป๊ฒ ๋๋์? ํ๋ ฌ์ด ํน์ด ํ๋ ฌ์ด๋ผ๋ ๋ป์ผ๋ก, ๋ฐฉ์ ์์ด ์๋ก ์ข ์์ด๊ฑฐ๋ ๋ชจ์์ด ๋ฉ๋๋ค. ์ด ๊ฒฝ์ฐ ์ ์ผํ ํด๊ฐ ์กด์ฌํ์ง ์์ผ๋ฏ๋ก ๊ณ์ฐ๊ธฐ๋ ํน์ด ํ๋ ฌ์์ ์๋ ค์ค๋๋ค.
์ ๋ถ๋ถ ํผ๋ฒํ ์ ์ฌ์ฉํ๋์? ์ ๋๊ฐ์ด ๊ฐ์ฅ ํฐ ๊ฐ์ ํผ๋ฒ์ผ๋ก ์ ํํ๋ฉด ๋ฐ์ฌ๋ฆผ ์ค์ฐจ๊ฐ ์ฆํญ๋๋ ๊ฒ์ ๋ง์ ์ ์์ด, ๊น๋ค๋ก์ด ํ๋ ฌ์์๋ ์ ํํ ๊ฒฐ๊ณผ๋ฅผ ์ป์ ์ ์์ต๋๋ค.
ํด๊ฐ ์ ์๊ฐ ์๋ ์๋ ์๋์? ๊ทธ๋ ์ต๋๋ค. ํด๋ ๋ถ๋์์์ ์ผ๋ก ๊ณ์ฐ๋์ด ์์๊ฐ ๋ ์ ์์ผ๋ฉฐ, ํ์ ์ ๋ฐ๋ ์ค์ ์ผ๋ก ์ ํจ์ซ์๋ฅผ ๋ช ์๋ฆฌ๊น์ง ๋ณด์ฌ์ค์ง ์กฐ์ ํ ์ ์์ต๋๋ค.