MCP๋กœ ์—ฐ๊ฒฐ โ†’

๊ณ„์‚ฐ ์ž…๋ ฅ

๊ณต์‹

๊ด‘๊ณ 

๊ฒฐ๊ณผ

์ ๋ถ„๊ฐ’
0.7468241328
๊ตฌ๊ฐ„ [a, b]์—์„œ ์ •์ ๋ถ„์˜ ๊ทผ์‚ฟ๊ฐ’
๊ณ„์‚ฐ ๋ฐฉ๋ฒ• Tanh-Sinh(์ด์ค‘ ์ง€์ˆ˜)
๋ถ„ํ•  ์ˆ˜ n 100
์‚ฌ์šฉ๋œ ํ•จ์ˆ˜ ํ‰๊ฐ€ ํšŸ์ˆ˜ 201

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\)์—์„œ ๋ฐœ์‚ฐํ•˜๋Š” ํ•จ์ˆ˜์กฐ์ฐจ ์ •ํ™•ํ•˜๊ฒŒ ์ ๋ถ„ํ•  ์ˆ˜ ์žˆ์Šต๋‹ˆ๋‹ค. ์ด ๋„๊ตฌ๋Š” ๋ณดํŽธ์ ์ธ ์ˆ˜ํ•™ ์›๋ฆฌ์— ๊ธฐ๋ฐ˜ํ•˜๋ฏ€๋กœ ์–ด๋””์„œ๋‚˜ ๋™์ผํ•˜๊ฒŒ ์ ์šฉ๋ฉ๋‹ˆ๋‹ค.

๋งˆ์ด๋„ˆ์Šค 1๊ณผ ํ”Œ๋Ÿฌ์Šค 1์˜ ์ ๊ทผ์„ ์„ ํ–ฅํ•ด ํ‰ํƒ„ํ•ด์ง€๋Š” S์žํ˜• tanh-sinh ๋ณ€ํ™˜ ๊ณก์„ 
์ด์ค‘ ์ง€์ˆ˜ ๋ณ€์ˆ˜ ๋ณ€ํ™˜์€ ์‹ค์ˆ˜ ์ง์„ ์„ (-1, 1)๋กœ ์‚ฌ์ƒํ•˜์—ฌ ๋…ธ๋“œ๋ฅผ ์–‘ ๋์  ์ชฝ์œผ๋กœ ๋ชจ์๋‹ˆ๋‹ค.

์‚ฌ์šฉ ๋ฐฉ๋ฒ•

ํ‘œ์ค€ ํ‘œ๊ธฐ๋ฒ•์œผ๋กœ ํ•จ์ˆ˜ \(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์œผ๋กœ ์–ธ๋”ํ”Œ๋กœ๋˜๋Š” ๋…ธ๋“œ(ํฌํ™”๋œ ๋์ )๋Š” ๊ฑด๋„ˆ๋›ฐ๋ฏ€๋กœ, ํŠน์ด์ ์ด ์žˆ์„ ์ˆ˜ ์žˆ๋Š” ๊ฒฝ๊ณ„์—์„œ์˜ ํ•จ์ˆ˜ ํ‰๊ฐ€๋ฅผ ํ”ผํ•  ์ˆ˜ ์žˆ์Šต๋‹ˆ๋‹ค.

a์—์„œ b๊นŒ์ง€์˜ ๊ตฌ๊ฐ„์—์„œ ์Œ์˜ ์˜์—ญ ์•„๋ž˜ ์–‘ ๋์  ๊ทผ์ฒ˜์— ๋นฝ๋นฝํ•˜๊ฒŒ ๋ชจ์ธ ๊ตฌ์  ๋…ธ๋“œ
๋…ธ๋“œ๊ฐ€ ๋์  a์™€ b ๊ทผ์ฒ˜์— ๋ชจ์ด๋ฏ€๋กœ ๋์  ํŠน์ด์ ์„ ์ •ํ™•ํ•˜๊ฒŒ ์ฒ˜๋ฆฌํ•ฉ๋‹ˆ๋‹ค.

๊ณ„์‚ฐ ์˜ˆ์‹œ

๊ตฌ๊ฐ„ [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}\) ์ธ์ˆ˜๊ฐ€ ๋ถ€ํ˜ธ๋ฅผ ๊ทธ๋Œ€๋กœ ๋ฐ˜์˜ํ•˜๋ฏ€๋กœ ๊ฒฐ๊ณผ์˜ ๋ถ€ํ˜ธ๊ฐ€ ์˜ฌ๋ฐ”๋ฅด๊ฒŒ ๋ฐ˜์ „๋ฉ๋‹ˆ๋‹ค.

์ตœ์ข… ์—…๋ฐ์ดํŠธ: