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

๊ณ„์‚ฐ ์ž…๋ ฅ

๊ณต์‹

๊ด‘๊ณ 

๊ฒฐ๊ณผ

ํ”„๋ ˆ๋„ฌ ์ฝ”์‚ฌ์ธ ์ ๋ถ„ C(x)
0.7798934004
๋ฌด์ฐจ์›(๋‹จ์œ„ ์—†์Œ)
์ •์˜ C(x) = โˆซโ‚€หฃ cos(ฯ€tยฒ/2) dt
๊ณ„์‚ฐ ๋ฐฉ๋ฒ• Composite Simpson's rule (asymptotic for |x| > 100)
๊ทนํ•œ x๊ฐ€ ยฑ๋ฌดํ•œ๋Œ€๋กœ ๊ฐˆ ๋•Œ C(x)๋Š” ยฑ0.5์— ์ˆ˜๋ ด

ํ”„๋ ˆ๋„ฌ ์ฝ”์‚ฌ์ธ ์ ๋ถ„์ด๋ž€?

ํ”„๋ ˆ๋„ฌ ์ฝ”์‚ฌ์ธ ์ ๋ถ„ \(C(x)\)๋Š” 0๋ถ€ํ„ฐ x๊นŒ์ง€ \(\cos(\pi t^{2}/2)\)๋ฅผ ์ ๋ถ„ํ•œ ๊ฐ’์œผ๋กœ ์ •์˜๋˜๋Š” ํŠน์ˆ˜ํ•จ์ˆ˜์ž…๋‹ˆ๋‹ค. ์ง์„  ๋ชจ์„œ๋ฆฌ์—์„œ ์ผ์–ด๋‚˜๋Š” ๊ทผ๊ฑฐ๋ฆฌ์žฅ ํšŒ์ ˆ์˜ ๊ฐ•๋„ ๋ถ„ํฌ๋ฅผ ๋น„๋กฏํ•ด ๊ด‘ํ•™ ์ „๋ฐ˜, ํŒŒ๋™ ๋ฌผ๋ฆฌํ•™, ๊ทธ๋ฆฌ๊ณ  ํ† ๋ชฉ๊ณตํ•™์—์„œ ํญ๋„“๊ฒŒ ๋“ฑ์žฅํ•ฉ๋‹ˆ๋‹ค. ํŠนํžˆ ํ† ๋ชฉ ๋ถ„์•ผ์—์„œ๋Š” ์ด์™€ ๋ฐ€์ ‘ํ•œ ํด๋กœ์†Œ์ด๋“œ(์˜ค์ผ๋Ÿฌ ๋‚˜์„ )๊ฐ€ ํ™œ์šฉ๋˜๋Š”๋ฐ, ๊ณก๋ฅ ์ด ํ˜ธ์˜ ๊ธธ์ด์— ๋”ฐ๋ผ ์„ ํ˜•์œผ๋กœ ์ฆ๊ฐ€ํ•˜๋Š” ์ด ๊ณก์„ ์€ ๊ณ ์†๋„๋กœ๋‚˜ ์ฒ ๋„์˜ ์™„๋งŒํ•œ ์™„ํ™”๊ณก์„ (์ „์ด๊ณก์„ )์„ ์„ค๊ณ„ํ•˜๋Š” ๋ฐ ์“ฐ์ž…๋‹ˆ๋‹ค.

x์— ๋Œ€ํ•œ ํ”„๋ ˆ๋„ฌ ์ฝ”์‚ฌ์ธ ์ ๋ถ„ C(x)์˜ ๊ทธ๋ž˜ํ”„
ํ”„๋ ˆ๋„ฌ ์ฝ”์‚ฌ์ธ ์ ๋ถ„ C(x)๋Š” ์ง„๋™ํ•˜๋ฉฐ x๊ฐ€ ์ปค์งˆ์ˆ˜๋ก 1/2๋กœ ์ˆ˜๋ ดํ•ฉ๋‹ˆ๋‹ค.

๊ณ„์‚ฐ๊ธฐ ์‚ฌ์šฉ ๋ฐฉ๋ฒ•

์ ๋ถ„ ์ƒํ•œ x์— ์ž„์˜์˜ ์‹ค์ˆ˜(์–‘์ˆ˜, ์Œ์ˆ˜, 0 ๋ชจ๋‘ ๊ฐ€๋Šฅ)๋ฅผ ์ž…๋ ฅํ•˜๋ฉด \(C(x)\) ๊ฐ’์ด ์ถœ๋ ฅ๋ฉ๋‹ˆ๋‹ค. x ์ž์ฒด๊ฐ€ ์ˆœ์ˆ˜ํ•œ ์ˆ˜์ด๋ฏ€๋กœ ๊ฒฐ๊ณผ๊ฐ’์—๋Š” ๋‹จ์œ„๊ฐ€ ์—†์Šต๋‹ˆ๋‹ค. \(|x|\)๊ฐ€ ์ปค์งˆ์ˆ˜๋ก \(C(x)\)๋Š” ์ง„๋™ํ•˜๋ฉด์„œ \(+0.5\)(x๊ฐ€ ์–‘์˜ ๋ฌดํ•œ๋Œ€๋กœ ๊ฐˆ ๋•Œ) ๋˜๋Š” \(-0.5\)(x๊ฐ€ ์Œ์˜ ๋ฌดํ•œ๋Œ€๋กœ ๊ฐˆ ๋•Œ)๋กœ ์ˆ˜๋ ดํ•ฉ๋‹ˆ๋‹ค.

๊ณต์‹๊ณผ ๊ทœ์•ฝ

์ด ๋„๊ตฌ๋Š” ์ฝ”์‚ฌ์ธ ์•ˆ์— \(\pi/2\) ์ธ์ž๋ฅผ ๋‘๋Š” ์ •๊ทœํ™”(normalized) ๊ทœ์•ฝ์„ ์‚ฌ์šฉํ•ฉ๋‹ˆ๋‹ค:

$$C(x) = \int_{0}^{x} \cos\!\left(\frac{\pi}{2}\,t^{2}\right) dt$$

์ด๋Š” ๋น„์ •๊ทœํ™” ํ˜•ํƒœ์ธ \(\int \cos(t^{2})\)์™€๋Š” ๋‹ค๋ฆ…๋‹ˆ๋‹ค. ๋‹ซํžŒ ํ˜•ํƒœ์˜ ํ•ด๊ฐ€ ์กด์žฌํ•˜์ง€ ์•Š๊ธฐ ๋•Œ๋ฌธ์—, ๊ฐ’์€ x์— ๋”ฐ๋ผ ์กฐ๋ฐ€ํ•ด์ง€๋Š” ๊ฒฉ์ž(๋ถ€๋ถ„๊ตฌ๊ฐ„ ์ˆ˜ \(n = \max(1000, \lceil 200\cdot|x| \rceil)\))๋ฅผ ์‚ฌ์šฉํ•œ ๋ณตํ•ฉ ์‹ฌํ”„์Šจ ๊ณต์‹์œผ๋กœ ๊ตฌํ•ฉ๋‹ˆ๋‹ค. \(|x|\)๊ฐ€ ๋งค์šฐ ํฐ ๊ฒฝ์šฐ์—๋Š” ๋ง‰๋Œ€ํ•œ ์ˆ˜์˜ ์ง„๋™์„ ์ผ์ผ์ด ์ ๋ถ„ํ•˜์ง€ ์•Š๋„๋ก ์ ๊ทผ ์ „๊ฐœ(asymptotic expansion)๋ฅผ ์ ์šฉํ•ฉ๋‹ˆ๋‹ค.

0๋ถ€ํ„ฐ x๊นŒ์ง€ cos(ฯ€tยฒ/2) ์•„๋ž˜์˜ ์Œ์˜ ์˜์—ญ
C(x)๋Š” 0๋ถ€ํ„ฐ x๊นŒ์ง€ cos(ฯ€tยฒ/2) ์•„๋ž˜์˜ ๋ถ€ํ˜ธ ์žˆ๋Š” ๋„“์ด์ž…๋‹ˆ๋‹ค.

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

x = 1์ผ ๋•Œ \(C(1) = \int_{0}^{1} \cos\!\left(\frac{\pi}{2}\,t^{2}\right) dt\)์ด๋ฉฐ, ์ˆ˜์น˜ ์ ๋ถ„ ๊ฒฐ๊ณผ ํ‘œ์ค€๊ฐ’ \(C(1) \approx 0.7798934004\)๋ฅผ ์–ป์Šต๋‹ˆ๋‹ค. x = 0.5์ผ ๋•Œ๋Š” \(C(0.5) \approx 0.4923442275\), x = 0์ผ ๋•Œ๋Š” \(C(0) = 0\)(์ •ํ™•๊ฐ’)์ž…๋‹ˆ๋‹ค.

์ž์ฃผ ๋ฌป๋Š” ์งˆ๋ฌธ

\(C(x)\)๋Š” ๊ธฐํ•จ์ˆ˜์ธ๊ฐ€์š”, ์šฐํ•จ์ˆ˜์ธ๊ฐ€์š”? ๊ธฐํ•จ์ˆ˜(odd function)์ž…๋‹ˆ๋‹ค. \(C(-x) = -C(x)\)๊ฐ€ ์„ฑ๋ฆฝํ•˜๋ฏ€๋กœ, ์Œ์ˆ˜๋ฅผ ์ž…๋ ฅํ•˜๋ฉด \(C(|x|)\)์˜ ๋ถ€ํ˜ธ๋ฅผ ๋’ค์ง‘์€ ๊ฐ’์ด ๋ฐ˜ํ™˜๋ฉ๋‹ˆ๋‹ค.

๋ฌดํ•œ๋Œ€์—์„œ์˜ ๊ทนํ•œ์€ ์–ผ๋งˆ์ธ๊ฐ€์š”? x๊ฐ€ ์–‘์˜ ๋ฌดํ•œ๋Œ€๋กœ ๊ฐˆ ๋•Œ \(C(x)\)๋Š” \(+1/2\)์—, ์Œ์˜ ๋ฌดํ•œ๋Œ€๋กœ ๊ฐˆ ๋•Œ \(-1/2\)์— ์ˆ˜๋ ดํ•ฉ๋‹ˆ๋‹ค.

๊ฒฐ๊ณผ๋Š” ์–ผ๋งˆ๋‚˜ ์ •ํ™•ํ•œ๊ฐ€์š”? ๋ฐฐ์ •๋ฐ€๋„(double-precision) ์‹ฌํ”„์Šจ ๋ฐฉ์‹์€ ์ผ๋ฐ˜์ ์ธ ์ž…๋ ฅ๊ฐ’์— ๋Œ€ํ•ด ์‹ ๋ขฐํ•  ์ˆ˜ ์žˆ๋Š” ์œ ํšจ์ˆซ์ž๋ฅผ ์•ฝ 10์ž๋ฆฌ๊นŒ์ง€ ์ œ๊ณตํ•ฉ๋‹ˆ๋‹ค. ์ง„์ •ํ•œ 50์ž๋ฆฌ ์ถœ๋ ฅ์„ ์–ป์œผ๋ ค๋ฉด ์ž„์˜ ์ •๋ฐ€๋„ ์—ฐ์‚ฐ์ด ํ•„์š”ํ•ฉ๋‹ˆ๋‹ค.

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