ํผ๋ฆฌ ๋ฐฉ๋ฒ์ด๋?
ํผ๋ฆฌ ๋ฐฉ๋ฒ(Halley's method)์ f(x) = 0 ํํ์ ๋ฐฉ์ ์์ ํธ๋ ๋ฐ๋ณต์ ์์นํด์ ๊ธฐ๋ฒ์ ๋๋ค. ๋ดํด-๋ฉ์จ(Newton-Raphson) ๋ฐฉ๋ฒ๊ณผ ๊ฐ๊น์ด ์น์ฒ ๊ฒฉ์ด์ง๋ง, ์๋ ด ์ฐจ์๊ฐ 3์ฐจ(3์ฐจ ์๋ ด)๋ผ๋ ์ ์ด ๋ค๋ฆ ๋๋ค. ๋ดํด๋ฒ์ด ํจ์์ 1์ฐจ ๋ํจ์๋ง ์ฌ์ฉํ๋ ๋ฐ๋ฉด, ํผ๋ฆฌ ๋ฐฉ๋ฒ์ ์ฌ๊ธฐ์ 2์ฐจ ๋ํจ์๊น์ง ํ์ฉํ๊ธฐ ๋๋ฌธ์ ๊ฐ์ ์ ํ๋์ ๋๋ฌํ๋ ๋ฐ ์ผ๋ฐ์ ์ผ๋ก ๋ ์ ์ ๋ฐ๋ณต ํ์๊ฐ ํ์ํฉ๋๋ค. ์ด ๋๊ตฌ๋ ํน์ ๊ตญ๊ฐ์ ๊ตญํ๋์ง ์๋ ๋ณดํธ์ ์ธ ์ํยท์์นํด์ ๊ณ์ฐ๊ธฐ๋ก ์ด๋์๋ ๋์ผํ๊ฒ ์ฌ์ฉํ ์ ์์ผ๋ฉฐ, ์ผ๊ฐํจ์ ์์ ๊ฐ๋๋ ๋ผ๋์(radian)์ผ๋ก ํด์ํฉ๋๋ค.
์ฌ์ฉ ๋ฐฉ๋ฒ
๊ทผ์ ๊ตฌํ๊ณ ์ถ์ ํจ์๋ฅผ \(f(x)\)์ ์ ๋ ฅํ ๋ค, ๊ทธ 1์ฐจ ๋ํจ์ \(f'(x)\)์ 2์ฐจ ๋ํจ์ \(f''(x)\)๋ฅผ \(x\)์ ๋ํ ์์ผ๋ก ๊ฐ๊ฐ ๋ฃ์ด ์ฃผ์ธ์. ์ง์ํ๋ ํ๊ธฐ๋ฒ์ผ๋ก๋ ๊ฑฐ๋ญ์ ๊ณฑ์ ์ํ +, -, *, /, ^, ๊ดํธ, ๊ทธ๋ฆฌ๊ณ sin, cos, tan, asin, acos, atan, sinh, cosh, tanh, exp, log/ln, log10, sqrt, abs ๊ฐ์ ํจ์์ ์์ pi, e๊ฐ ์์ต๋๋ค. ์ฐพ์ผ๋ ค๋ ๊ทผ์ ๊ฐ๊น์ด ์ด๊ธฐ ์ถ์ ๊ฐ \(x_0\)๋ฅผ ์ ํ๊ณ , ์ต๋ ๋ฐ๋ณต ํ์ \(n\)๊ณผ ํ์ํ ์ ํจ์ซ์ ์๋ฆฟ์๋ฅผ ์ ํํ๋ฉด ๋ฉ๋๋ค.
๊ณต์ ํ์ด
๋งค ๋จ๊ณ๋ง๋ค $$x_{n+1} = x_n - \frac{2\,f(x_n)\,f'(x_n)}{2\,[f'(x_n)]^2 - f(x_n)\,f''(x_n)}$$ ๋ฅผ ๊ณ์ฐํฉ๋๋ค. ๋ถ์๋ ๊ธฐ์กด ๋ดํด ๋ณด์ ํญ(์ ๋ฐฐ์จ ์กฐ์ ํ ํํ)์ด๊ณ , ๋ถ๋ชจ์ ์ถ๊ฐ๋ \(- f(x_n)\,f''(x_n)\) ํญ์ด ํจ์ \(f\)์ ๊ณก๋ฅ ์ ๋ณด์ ํด ์ฃผ์ด ์๋ ด์ ๊ฐ์ํฉ๋๋ค. ๋ฐ๋ณต์ \(x\)์ ๋ณํ๋์ด๋ ์์ฐจ \(f(x)\)๊ฐ ์์ฃผ ์์ ํ์ฉ์ค์ฐจ ์๋๋ก ๋ด๋ ค๊ฐ๊ฑฐ๋, ๋ฐ๋ณต ํ๋์ ๋๋ฌํ๋ฉด ๋ฉ์ถฅ๋๋ค. ๋ง์ฝ ๋ถ๋ชจ๊ฐ 0์ด ๋๋ฉด ๋ฐฉ๋ฒ์ด ์๋ํ์ง ์์ผ๋ฏ๋ก, ๋ค๋ฅธ \(x_0\) ๊ฐ์ผ๋ก ๋ค์ ์๋ํด์ผ ํฉ๋๋ค.
๊ณ์ฐ ์์
\(f(x) = x - \cos(x)\), \(f'(x) = 1 + \sin(x)\), \(f''(x) = \cos(x)\)์์ \(x_0 = 1\)๋ก ์์ํด ๋ณด๊ฒ ์ต๋๋ค. ์ฒซ ๋จ๊ณ๋ $$x_1 = 1 - \frac{2 \times 0.4596977 \times 1.8414710}{2 \times 1.8414710^2 - 0.4596977 \times 0.5403023} = 1 - \frac{1.6930504}{6.5336550} = 0.7408769$$ ๊ฐ ๋ฉ๋๋ค. ์ดํ ๋ฐ๋ณต์ ๋น ๋ฅด๊ฒ ๋ํฐ ์(Dottie number) \(x = 0.7390851332151607\)๋ก ์๋ ดํ๋๋ฐ, ์ด๋ \(x = \cos(x)\)์ ์ ์ผํ ํด์ ๋๋ค.
์์ฃผ ๋ฌป๋ ์ง๋ฌธ
ํผ๋ฆฌ ๋ฐฉ๋ฒ์ ๋ดํด๋ฒ๊ณผ ์ด๋ป๊ฒ ๋ค๋ฅธ๊ฐ์? ๋ดํด๋ฒ์ ๊ณก๋ฅ ์ ๊ณ ๋ คํ์ง ์์ง๋ง, ํผ๋ฆฌ ๋ฐฉ๋ฒ์ 2์ฐจ ๋ํจ์ ํญ์ ๋ํด 2์ฐจ ์๋ ด์ด ์๋ 3์ฐจ ์๋ ด์ ๋ฌ์ฑํฉ๋๋ค. ๊ทธ ๊ฒฐ๊ณผ ์ ํํ ์๋ฆฟ์ ํ๋๋ฅผ ์ป๋ ๋ฐ ๋ณดํต ๋ ์ ์ ๋ฐ๋ณต์ด ํ์ํฉ๋๋ค.
์ ๋ํจ์๋ฅผ ์ง์ ์ ๋ ฅํด์ผ ํ๋์? ์ด ๊ณ์ฐ๊ธฐ๋ ์ฌ์ฉ์๊ฐ ์ ๋ ฅํ ๋ํจ์๋ฅผ ๊ทธ๋๋ก ์ ๋ขฐํฉ๋๋ค. ๋ํจ์๊ฐ ํ๋ฆฌ๋ฉด ์๋ ด์ด ๋๋น ์ง๊ฑฐ๋ ์คํจํฉ๋๋ค. \(f(x)\)๋ฅผ ์ ์คํ๊ฒ ๋ฏธ๋ถํด \(f'(x)\)์ \(f''(x)\)๋ฅผ ์ ํํ ๊ตฌํด ์ ๋ ฅํ์ธ์.
์๋ ดํ์ง ์์ผ๋ฉด ์ด๋ป๊ฒ ํ๋์? ์ด ๋ฐฉ๋ฒ์ \(x_0\)์ ๊ฐ์ฅ ๊ฐ๊น์ด ๊ทผ์ ์ฐพ์ต๋๋ค. ์์์ ์ด ์ข์ง ์์ผ๋ฉด ๋ฐ์ฐํ๊ฑฐ๋ ๋ค๋ฅธ ๊ทผ์ผ๋ก ๋น ์ง ์ ์๊ณ , ๋ถ๋ชจ๊ฐ 0์ด ๋๋ฉด ๊ณ์ฐ์ด ์์ ํ ๋ฉ์ถฅ๋๋ค. ์ด๋ด ๋๋ \(x_0\)๋ฅผ ๋ฐ๊พธ๊ฑฐ๋ ๋ํจ์ ์์ ๋ค์ ํ์ธํด ๋ณด์ธ์.