์ด๋ถ๋ฒ์ด๋?
์ด๋ถ๋ฒ์ ์ฐ์ํจ์ f(x)์ ๊ทผ, ์ฆ f(x) = 0์ด ๋๋ x ๊ฐ์ ์ฐพ๋ ๊ฐ์ฅ ์ค๋๋๊ณ ์์ ์ ์ธ ๋ฐฉ๋ฒ ์ค ํ๋์ ๋๋ค. ํจ์์ ๋ถํธ๊ฐ ๋ฐ๋๋ ๊ตฌ๊ฐ [a, b]์์ ์ถ๋ฐํด, ๊ทธ ๊ตฌ๊ฐ์ ์ ๋ฐ์ผ๋ก ๊ณ์ ๋๋๋ฉด์ ์ฌ์ ํ ๊ทผ์ ํ๊ณ ์๋ ์ชฝ ์ ๋ฐ๋ง ๋จ๊ธฐ๋ ๋ฐฉ์์ผ๋ก ์๋ํฉ๋๋ค. ๋งค ๋จ๊ณ๋ง๋ค ๋ถํธ ๋ณํ๊ฐ ์ ์ง๋๊ธฐ ๋๋ฌธ์, f๊ฐ ์ฐ์์ด๊ณ f(a)์ f(b)์ ๋ถํธ๊ฐ ์๋ก ๋ฐ๋์ด๊ธฐ๋ง ํ๋ฉด ๋ฐ๋์ ์๋ ดํ๋ ๊ฒ์ด ๋ณด์ฅ๋ฉ๋๋ค. ์ด ๊ณ์ฐ๊ธฐ๋ ๋ฌด์ฐจ์์ ์คํจ์๋ผ๋ฉด ๋ชจ๋ ์ฒ๋ฆฌํ ์ ์์ผ๋ฉฐ, ์ผ๊ฐํจ์๋ ๋ผ๋์์ ๊ธฐ์ค์ผ๋ก ๊ณ์ฐํฉ๋๋ค.
๊ณ์ฐ๊ธฐ ์ฌ์ฉ๋ฒ
f(x) ์นธ์ x๋ฅผ ๋ณ์๋ก ์ฌ์ฉํด ํจ์๋ฅผ ์
๋ ฅํ์ธ์(์: x-cos(x), x^2-2, exp(x)-3). ๊ทผ์ด ๋ ๋์ ์ฌ์ด์ ์ค๋๋ก ํํ a์ ์ํ b๋ฅผ ์ค์ ํ๋ฉด ๋ฉ๋๋ค. ์ด๋ \(f(a)\cdot f(b)\)์ ๊ฐ์ด 0 ์ดํ์ฌ์ผ ํฉ๋๋ค. ์ต๋ ๋ฐ๋ณต ํ์ n๊ณผ ํ์ํ ์๋ฆฟ์๋ฅผ ๊ณ ๋ฅธ ๋ค ๊ณ์ฐํ๋ฉด, ๊ทผ์ ๊ทผ์ฟ๊ฐ๊ณผ ๊ทธ ์ง์ ์์์ ํจ์ซ๊ฐ(0์ ๊ฐ๊น์์ผ ํฉ๋๋ค), ๊ทธ๋ฆฌ๊ณ ์ค์ ๋ก ํ์ํ๋ ๋ฐ๋ณต ํ์๋ฅผ ์๋ ค์ค๋๋ค.
๊ณ์ฐ ๊ณต์
๊ฐ ๋จ๊ณ์์ ๊ทผ์ ์ถ์ ๊ฐ์ ๊ตฌ๊ฐ์ ์ค์ \(x_n = (a_n + b_n) / 2\)์ ๋๋ค.
$$c = \frac{a + b}{2} \\[1.5em] \text{where}\quad \left\{ \begin{aligned} &\text{if } f(a)\cdot f(c) < 0 \Rightarrow b \leftarrow c \\ &\text{else} \Rightarrow a \leftarrow c \\ &\text{repeat up to } n \text{ times} \end{aligned} \right.$$\(|f(x_n)|\)์ด ํ์ฉ ์ค์ฐจ๋ณด๋ค ์์์ง๋ฉด \(x_n\)์ ๊ทผ์ผ๋ก ๋ฐ์๋ค์ ๋๋ค. ๊ทธ๋ ์ง ์์ผ๋ฉด ๋ถํธ ๋ณํ๊ฐ ๋จ์ ์๋ ์ชฝ ์ ๋ฐ์ ํํฉ๋๋ค. ์ฆ, \(f(a_n)\cdot f(x_n) > 0\)์ด๋ฉด ๊ทผ์ ์ค๋ฅธ์ชฝ์ ์์ผ๋ฏ๋ก \(a = x_n\)์ผ๋ก, ๊ทธ๋ ์ง ์์ผ๋ฉด ์ผ์ชฝ์ ์์ผ๋ฏ๋ก \(b = x_n\)์ผ๋ก ๊ฐฑ์ ํฉ๋๋ค. ๊ตฌ๊ฐ ํญ์ \((b-a)/2^n\)์ผ๋ก ์ค์ด๋ค๊ธฐ ๋๋ฌธ์ ์๋ ด์ ์ ํ์ด๋ฉฐ, ํ ๋จ๊ณ๋ง๋ค ๋๋ต ์ ํํ ์ด์ง ์๋ฆฟ์ ํ๋์ฉ์ ์ป๋ ์ ์ ๋๋ค.
์์ ํ์ด
๊ตฌ๊ฐ \([-10, 10]\)์์ \(f(x) = x - \cos(x)\)๋ฅผ ์๊ฐํด ๋ด ์๋ค. \(f(-10) \approx -10.84\)(์์), \(f(10) \approx 10.84\)(์์)์ด๋ฏ๋ก ๊ทผ์ด ๊ตฌ๊ฐ ์์ ๊ฐํ ์์ต๋๋ค. ๊ตฌ๊ฐ์ ๋ฐ๋ณตํด์ ์ ๋ฐ์ผ๋ก ๋๋๋ฉด \(x \approx 0.7390851332151607\)๋ก ์๋ ดํ๋๋ฐ, ์ด๋ \(x = \cos x\)๋ฅผ ๋ง์กฑํ๋ ์ ๋ช ํ ๋ํฐ ์(Dottie number)์ด๋ฉฐ ์ด๋ \(f(x)\)๋ ์ฌ์ค์ 0์ด ๋ฉ๋๋ค.
์ด๋ถ๋ฒ์ ์์ผ๋ก ํ๋ ๋ฐฉ๋ฒ
์ด๋ถ๋ฒ์ ๊ทผ์ด ์๋ ๊ฒ์ผ๋ก ์๋ ค์ง ๊ตฌ๊ฐ์ ๋ฐ๋ณต์ ์ผ๋ก ์ด๋ฑ๋ถํ์ฌ \(f(x)=0\)์ ๊ทผ์ ์ฐพ์ต๋๋ค. ์ด๊ฒ์ ์ค๊ฐ๊ฐ ์ ๋ฆฌ์ ์์กดํฉ๋๋ค: \(f\)๊ฐ \([a,b]\)์์ ์ฐ์์ด๊ณ \(f(a)\)์ \(f(b)\)๊ฐ ๋ฐ๋ ๋ถํธ๋ฅผ ๊ฐ์ง๋ฉด, ๊ทผ์ด ๋ฐ๋์ ๊ทธ ์ฌ์ด์ ์์ด์ผ ํฉ๋๋ค.
- ๊ดํธ๋ฅผ ํ์ธํ์ญ์์ค. \(f(a)\cdot f(b)<0\)์์ ํ์ธํฉ๋๋ค. ๊ณฑ์ด ์์์ด๋ฉด ๊ตฌ๊ฐ์ด ๊ทผ์ ํฌํจํ๋ค๊ณ ๋ณด์ฅ๋์ง ์์ผ๋ฏ๋ก ๋ค๋ฅธ \([a,b]\)๋ฅผ ์ ํํ์ญ์์ค.
- ์ค์ ์ ๊ณ์ฐํ์ญ์์ค. \(m=\dfrac{a+b}{2}\).
- ํจ์๋ฅผ ํ๊ฐํ์ญ์์ค. \(f(m)\)์ ๊ตฌํฉ๋๋ค. \(f(m)=0\)์ด๊ฑฐ๋ ํ์ฉ ๋ฒ์ ๋ด์ ์์ผ๋ฉด \(m\)์ด ๊ทผ์ด๊ณ ๋ฉ์ถฅ๋๋ค.
- ๋ถํธ์ ๋ฐ๋ผ ๋์ ์ ๋ฐ๊พธ์ญ์์ค. \(f(a)\cdot f(m)<0\)์ด๋ฉด ๊ทผ์ \([a,m]\)์ ์์ผ๋ฏ๋ก \(b\leftarrow m\)์ ์ค์ ํฉ๋๋ค. ๊ทธ๋ ์ง ์์ผ๋ฉด ๊ทผ์ \([m,b]\)์ ์์ผ๋ฏ๋ก \(a\leftarrow m\)์ ์ค์ ํฉ๋๋ค.
- ๋ฐ๋ณต \(|b-a|<\text{ํ์ฉ๊ฐ}\) ๋๋ \(|f(m)|<\text{ํ์ฉ๊ฐ}\)์ด ๋ ๋๊น์ง, ๋๋ ์ต๋ ๋ฐ๋ณต ํ์์ ๋๋ฌํ ๋๊น์ง 2๋จ๊ณ๋ถํฐ 4๋จ๊ณ๊น์ง ๋ฐ๋ณตํฉ๋๋ค.
์: \(f(x)=x^{3}-x-2\)๋ฅผ \([1,2]\)์์ ๊ณ์ฐํฉ๋๋ค. ํ์ธ: \(f(1)=-2\), \(f(2)=4\), ๊ณฑ \(<0\) โ ๊ดํธ๊ฐ ์ ํจํฉ๋๋ค.
| ๋ฐ๋ณต | a | b | m=(a+b)/2 | f(m) | ์ ๊ตฌ๊ฐ |
|---|---|---|---|---|---|
| 1 | 1.0000 | 2.0000 | 1.5000 | โ0.125 | [1.5, 2] |
| 2 | 1.5000 | 2.0000 | 1.7500 | 1.6094 | [1.5, 1.75] |
| 3 | 1.5000 | 1.7500 | 1.6250 | 0.6660 | [1.5, 1.625] |
| 4 | 1.5000 | 1.6250 | 1.5625 | 0.2522 | [1.5, 1.5625] |
| 5 | 1.5000 | 1.5625 | 1.5313 | 0.0591 | [1.5, 1.5313] |
| 6 | 1.5000 | 1.5313 | 1.5156 | โ0.0340 | [1.5156, 1.5313] |
๋ ๋ง์ ๋ฐ๋ณต ํ ๊ตฌ๊ฐ์ด ์ฐธ ๊ทผ์ธ 1.521380์ผ๋ก ์๋ ดํฉ๋๋ค. ์ด๋ 3์ฐจ ๋ฐฉ์ ์ \(x^{3}-x-2=0\)์ ์ ์ผํ ์ค์ ๊ทผ์ด๊ธฐ๋ ํ๋ฉฐ, ์ง์ ํ์ด๊ธฐ๋ก ์ฐพ์ 1.521380์ ๋๋ค.
๋ฐ๋ณต ํ์ ๋ ํ์ฉ๊ฐ ๋ฐ ๊ตฌ๊ฐ ๋๋น
๊ฐ ์ด๋ถ ๋จ๊ณ์์ ๊ตฌ๊ฐ์ ์ด๋ฑ๋ถํ๋ฏ๋ก \(n\)๋ฒ ๋ฐ๋ณต ํ ๊ดํธ ๋๋น๋ \((b-a)/2^{\,n}\)์ ๋๋ค. ๊ทผ์ ์์น์ ๋ํด ํ์ฉ๊ฐ \(\text{ํ์ฉ๊ฐ}\)์ ๋๋ฌํ๋ ค๋ฉด ๋๋ต
$$n \approx \log_2\!\left(\frac{b-a}{\text{ํ์ฉ๊ฐ}}\right).$$ํ์๋ ์ ๋ฐ๋์ ๋ก๊ทธ์๋ง ์ฆ๊ฐํ๋ฏ๋ก ๋งค์ฐ ์๊ฒฉํ ํ์ฉ๊ฐ๋ ์๋์ ์ผ๋ก ์ ์ ๋จ๊ณ๋ง ํ์ํฉ๋๋ค. ํ๋ ์ฌ๋ฌ ์กฐํฉ์ ๋ํด ๋ฐ์ฌ๋ฆผ๋ ๋ฐ๋ณต ํ์๋ฅผ ๋ณด์ฌ์ค๋๋ค.
| ์ด๊ธฐ ๋๋น \(b-a\) | ๋ชฉํ ํ์ฉ๊ฐ | \(\log_2((b-a)/\text{ํ์ฉ๊ฐ})\) | ํ์ํ ๋ฐ๋ณต ํ์ |
|---|---|---|---|
| 1 | \(10^{-3}\) | 9.97 | 10 |
| 1 | \(10^{-6}\) | 19.93 | 20 |
| 1 | \(10^{-10}\) | 33.22 | 34 |
| 10 | \(10^{-6}\) | 23.25 | 24 |
| 20 | \(10^{-6}\) | 24.25 | 25 |
| 100 | \(10^{-8}\) | 33.22 | 34 |
| 0.5 | \(10^{-12}\) | 38.86 | 39 |
์๋ฅผ ๋ค์ด, ๋๋น 20์ธ ๊ดํธ๋ฅผ \(10^{-6}\)๋ก ์ ์ ํ๋ ค๋ฉด \(\lceil\log_2(20/10^{-6})\rceil=\lceil 24.25\rceil=\) 25๋ฒ์ ๋ฐ๋ณต์ด ํ์ํ๊ณ , ๋๋น 1์ธ ๊ดํธ๋ฅผ \(10^{-10}\)๋ก ์ ์ ํ๋ ค๋ฉด \(\lceil 33.22\rceil=\) 34๋ฒ์ด ํ์ํฉ๋๋ค. ์์ ๋๋น๋ฅผ ์ด๋ฑ๋ถํ๋ฉด ์ ํํ ํ ๋ฒ์ ๋ฐ๋ณต์ ์ ์ฝํฉ๋๋ค. ์ ๋ฐ๋๋ฅผ ์ ๊ณฑํ๋ฉด(์์์ ํ ์๋ฆฌ ์ถ๊ฐ) ์ฝ 3.3๋ฒ์ ๋ฐ๋ณต์ด ์์๋ฉ๋๋ค.
์ฃผ์ ์ฉ์ด
- ๊ทผ. ํจ์๊ฐ 0์ธ ๊ฐ \(x^{*}\), ์ฆ \(f(x^{*})=0\). ์์ ๋๋ ๋ฐฉ์ ์์ ํด๋ผ๊ณ ๋ ํฉ๋๋ค.
- ๊ดํธ / ๊ตฌ๊ฐ \([a,b]\). ๊ทผ์ ํฌํจํ๋ค๊ณ ์๊ฐ๋๋ ๋์ ์์ ๋๋ค. ์ด๋ถ๋ฒ์ ๊ฒฝ์ฐ ๋ถํธ ๋ณํ ์กฐ๊ฑด์ ๋ง์กฑํด์ผ ํฉ๋๋ค.
- ๋ถํธ ๋ณํ. ์กฐ๊ฑด \(f(a)\cdot f(b)<0\), ์ฆ \(f\)๊ฐ ๋์ ์์ ๋ฐ๋ ๋ถํธ๋ฅผ ์ทจํ๋ ๊ฒ์ ์๋ฏธํฉ๋๋ค. ์ฐ์ \(f\)์ ๊ฒฝ์ฐ ์ด๊ฒ์ ๊ทธ ์ฌ์ด์ ์ต์ํ ํ๋์ ๊ทผ์ด ์์์ ๋ณด์ฅํฉ๋๋ค(์ค๊ฐ๊ฐ ์ ๋ฆฌ).
- ์ค์ . ํ์ฌ ๊ตฌ๊ฐ์ ์ค์ฌ, \(m=(a+b)/2\). ๊ฐ ๋จ๊ณ์์ ์ด ์ ์ ๊ฒ์ฌํ๊ณ ๊ทผ์ ํฌํจํ ์ ์๋ ์ ๋ฐ์ ๋ฒ๋ฆฝ๋๋ค.
- ํ์ฉ๊ฐ. ๋ฐ๋ณต์ ๋ฉ์ถ๋ ์ ํ๋ ๋ชฉํ๋ก, ๊ตฌ๊ฐ ๋๋น \(|b-a|\) ๋๋ ์์ฐจ \(|f(m)|\)์ ์ ์ฉ๋ฉ๋๋ค.
- ์๋ ด(์ ํ). ์ด๋ถ๋ฒ์ ์ ํ์ผ๋ก ์๋ ดํฉ๋๋ค: ์ค์ฐจ๋ ๊ฐ ๋จ๊ณ์์ ๋๋ต ์ด๋ฑ๋ถ๋ฉ๋๋ค(์ค์ฐจ \(\le (b-a)/2^{n}\)), ์ผ์ ํ์ง๋ง ๊ฐ์ํ์ง ์๋ ์งํ์ ์ ๊ณตํฉ๋๋ค โ ๋ฐ๋ณต๋น ์ฝ 1๊ฐ์ ์ถ๊ฐ ์ฌ๋ฐ๋ฅธ ์ด์ง ์๋ฆฟ์์ ๋๋ค.
- ๋ฐ๋ณต. ์ค์ ๊ณ์ฐ, ํจ์ ํ๊ฐ, ๋์ ์ ๋ฐ์ดํธ์ ํ ์ ์ฒด ์ฃผ๊ธฐ์ ๋๋ค. ๋ฐ๋ณต ํ์๋ ์ต๋ ๋ฐ๋ณต ์ค์ ์ผ๋ก ์ ํ๋ฉ๋๋ค.
์์ฃผ ๋ฌป๋ ์ง๋ฌธ
"๋ถํธ ๋ณํ ์์" ์ค๋ฅ๊ฐ ์ ๋์ค๋์? ์ด๋ถ๋ฒ์ f(a)์ f(b)์ ๋ถํธ๊ฐ ์๋ก ๋ฐ๋์ฌ์ผ ํฉ๋๋ค. ๊ทผ์ ์ฌ์ด์ ๋๋๋ก a์ b๋ฅผ ์กฐ์ ํด ๋ณด์ธ์.
๊ทผ์ ์ฌ๋ฌ ๊ฐ ์ฐพ์ ์ ์๋์? ์๋์. ๊ตฌ๊ฐ ์์ ๊ทผ ํ๋๋ง ๋ฐํํ๋ฉฐ, ๊ณก์ ์ด ์ถ์ ๋์ง ์๊ณ ๋ฟ๊ธฐ๋ง ํ๋ ๊ทผ์ ์ฐพ์๋ด์ง ๋ชปํฉ๋๋ค.
๋ดํด๋ฒ๋ณด๋ค ์ ๋๋ฆฐ๊ฐ์? ์ด๋ถ๋ฒ์ ํ ๋จ๊ณ๋ง๋ค ์ด์ง ์๋ฆฟ์ ํ๋ ์ ๋๋ฅผ ์ป๋ ์ ํ ์๋ ด์ธ ๋ฐ๋ฉด, ๋ดํด๋ฒ์ 2์ฐจ ์๋ ด์ด๊ธฐ ๋๋ฌธ์ ๋๋ค. ๊ทธ๋์ ์ด๋ถ๋ฒ์ ์์ ํ ์ถ๋ฐ์ ์ ๋จผ์ ์ฐพ์ ๋๊ณ , ๋ ๋น ๋ฅธ ๋ฐฉ๋ฒ์ผ๋ก ์ ๋ฐํ๊ฒ ๋ค๋ฌ๋ ๋ฐ ์์ฃผ ์ฐ์ ๋๋ค.