Tanh-Sinh ๊ตฌ์ ๋ฒ์ด๋?
Tanh-Sinh ๊ตฌ์ ๋ฒ์ ์ด์ค ์ง์(Double Exponential, DE) ๊ณต์์ด๋ผ๊ณ ๋ ๋ถ๋ฆฌ๋ฉฐ, ์ ํ ๊ตฌ๊ฐ [a, b]์์์ ์ ์ ๋ถ์ ๊ณ์ฐํ๋ ๋ฐ ํ์ํ ์์น ์ ๋ถ ๊ธฐ๋ฒ์ ๋๋ค. ํนํ ํผ์ ๋ถํจ์๊ฐ ๊ตฌ๊ฐ์ ์ ๋์ ์์ ํน์ด์ ์ ๊ฐ์ง ๋ ์ง๊ฐ๋ฅผ ๋ฐํํฉ๋๋ค. ์ด ๋ฐฉ๋ฒ์ \(u = \tanh\!\left(\tfrac{\pi}{2}\sinh(t)\right)\)๋ผ๋ ๋ณ์ ๋ณํ์ ์ ์ฉํ์ฌ ๋์ ์ \(t = \pm\infty\)๋ก ๋ณด๋ ๋๋ค. ๊ทธ๋ฌ๋ฉด ๋์ ๊ทผ์ฒ์์ ํผ์ ๋ถํจ์์ ๊ธฐ์ฌ๋๊ฐ ์ด์ค ์ง์์ ์ผ๋ก ๊ฐ์ํ๊ธฐ ๋๋ฌธ์, \(\frac{1}{\sqrt{1-x^2}}\)์ฒ๋ผ \(a\)๋ \(b\)์์ ๋ฐ์ฐํ๋ ํจ์์กฐ์ฐจ ์ ํํ๊ฒ ์ ๋ถํ ์ ์์ต๋๋ค. ์ด ๋๊ตฌ๋ ๋ณดํธ์ ์ธ ์ํ ์๋ฆฌ์ ๊ธฐ๋ฐํ๋ฏ๋ก ์ด๋์๋ ๋์ผํ๊ฒ ์ ์ฉ๋ฉ๋๋ค.
์ฌ์ฉ ๋ฐฉ๋ฒ
ํ์ค ํ๊ธฐ๋ฒ์ผ๋ก ํจ์ \(f(x)\)๋ฅผ ์ ๋ ฅํ์ธ์(์ฐ์ฐ์ + - * / ^, ๊ดํธ, ๊ทธ๋ฆฌ๊ณ sin, cos, exp, log, sqrt, abs ๊ฐ์ ํจ์์ ์์ pi, e๋ฅผ ์ฌ์ฉํ ์ ์์ต๋๋ค). ์๋์ชฝ ์ ๋ถ ํ๊ณ \(a\), ์์ชฝ ์ ๋ถ ํ๊ณ \(b\), ๊ทธ๋ฆฌ๊ณ ๋ ธ๋ ๋ฐ๋๋ฅผ ์กฐ์ ํ๋ ๋ถํ ์ \(n\)์ ์ ๋ ฅํฉ๋๋ค. \(n\)์ด ํด์๋ก ๊ณ์ฐ ๋น์ฉ์ ๋์ด๋์ง๋ง ์ ํ๋๊ฐ ๋์์ง๋ฉฐ, ์ค๋ฌด์์๋ ๋ณดํต 50โ400 ์ ๋๋ฅผ ์ฌ์ฉํฉ๋๋ค. ํผ์ ๋ถํจ์๋ ์ด๋ฆฐ๊ตฌ๊ฐ์์ ํด์์ ์ด์ด์ผ ํ๊ณ (๋์ ํน์ด์ ์ ๊ด์ฐฎ์ต๋๋ค), ์ฃผ๊ธฐํจ์์ฌ์๋ ์ ๋ฉ๋๋ค.
๊ณต์ ์ค๋ช
๋จผ์ ๊ตฌ๊ฐ์ \(x = \frac{b-a}{2}u + \frac{a+b}{2}\), \(dx = \frac{b-a}{2}\,du\)๋ฅผ ํตํด [-1, 1]๋ก ์ ๊ทํํฉ๋๋ค. ์ด์ด์ DE ๊ณต์์ ๋ ธ๋ \(t_k = k\,h\), ๊ฐ๋ก ์ขํ \(u_k = \tanh\!\left(\tfrac{\pi}{2}\sinh(t_k)\right)\), ๊ฐ์ค์น \(w_k = \dfrac{\tfrac{\pi}{2}\cosh(t_k)}{\cosh^{2}\!\left(\tfrac{\pi}{2}\sinh(t_k)\right)}\)๋ฅผ ์ฌ์ฉํฉ๋๋ค. ์ ๋ถ๊ฐ์ ๋ค์๊ณผ ๊ฐ์ด ๊ทผ์ฌ๋ฉ๋๋ค.
$$\int_{a}^{b} f(x)\,dx \approx \frac{b-a}{2}\, h \sum_{k=-N}^{N} w_k\, f\!\left(\frac{a+b}{2} + \frac{b-a}{2}\,x_k\right)$$๊ฐ์ค์น๊ฐ 0์ผ๋ก ์ธ๋ํ๋ก๋๋ ๋ ธ๋(ํฌํ๋ ๋์ )๋ ๊ฑด๋๋ฐ๋ฏ๋ก, ํน์ด์ ์ด ์์ ์ ์๋ ๊ฒฝ๊ณ์์์ ํจ์ ํ๊ฐ๋ฅผ ํผํ ์ ์์ต๋๋ค.
๊ณ์ฐ ์์
๊ตฌ๊ฐ [0, 1]์์ \(f(x) = \exp(-x^2)\)์ ์ ๋ถํด ๋ด ์๋ค. ์ ํํ ๊ฐ์ \(\frac{\sqrt{\pi}}{2}\cdot\operatorname{erf}(1) \approx 0.7468241\)์ ๋๋ค. ๊ฑฐ์น ๊ฐ๊ฒฉ(\(h = 0.5\), \(N = 4\))๋ง์ผ๋ก๋ ์ฝ 0.7467์ด๋ผ๋ ๊ฐ์ ์ป์ ์ ์์ผ๋ฉฐ, ๊ธฐ๋ณธ๊ฐ \(n = 100\)์ ์ฌ์ฉํ๋ฉด ์ฝ 12์๋ฆฌ๊น์ง ์ ํํ๊ฒ ์ผ์นํฉ๋๋ค.
์์ฃผ ๋ฌป๋ ์ง๋ฌธ
๋์ ์ ํน์ด์ ์ด ์๋ ํจ์๋ ์ ๋ถํ ์ ์๋์? ๋ค โ ๋ฐ๋ก ์ด๊ฒ์ด ์ด ๋ฐฉ๋ฒ์ ํต์ฌ ๊ฐ์ ์ ๋๋ค. \(a\)๋ \(b\)์ ์๋ ์ ๋ถ ๊ฐ๋ฅํ ํน์ด์ ์ ๋งค๋๋ฝ๊ฒ ์ฒ๋ฆฌํฉ๋๋ค.
์ ์ฃผ๊ธฐ์ฑ์ด ์ค์ํ๊ฐ์? ์ด์ค ์ง์ ๊ณต์์ ๋น์ฃผ๊ธฐ ํผ์ ๋ถํจ์์ ์ต์ ํ๋์ด ์์ต๋๋ค. ์ฃผ๊ธฐํจ์์ ๊ฒฝ์ฐ์๋ ์ฌ๋ค๋ฆฌ๊ผด ๊ณต์์ด ๋ ๋น ๋ฅด๊ฒ ์๋ ดํ๋ฉฐ, DE ๊ณต์์ ๋ถ์ ํํ ์ ์์ต๋๋ค.
a์ b๊ฐ ๊ฐ์ผ๋ฉด ์ด๋ป๊ฒ ๋๋์? ์ ๋ถ๊ฐ์ 0์ด ๋ฉ๋๋ค. \(a > b\)์ธ ๊ฒฝ์ฐ์๋ \(\frac{b-a}{2}\) ์ธ์๊ฐ ๋ถํธ๋ฅผ ๊ทธ๋๋ก ๋ฐ์ํ๋ฏ๋ก ๊ฒฐ๊ณผ์ ๋ถํธ๊ฐ ์ฌ๋ฐ๋ฅด๊ฒ ๋ฐ์ ๋ฉ๋๋ค.