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

๊ณ„์‚ฐ ์ž…๋ ฅ

๊ณต์‹

๊ด‘๊ณ 

๊ฒฐ๊ณผ

์ œ1์ข… ๋ณ€ํ˜• ๋ฒ ์…€ ํ•จ์ˆ˜ I_v(x)
1.2660658778
๋ฌด์ฐจ์›
K_v(x) โ€” ์ œ2์ข… 0.4210244382
I'_v(x) โ€” I์˜ ๋„ํ•จ์ˆ˜ 0.565159104
K'_v(x) โ€” K์˜ ๋„ํ•จ์ˆ˜ -0.6019072302

์ด ๊ณ„์‚ฐ๊ธฐ์˜ ๊ธฐ๋Šฅ

์ด ๋„๊ตฌ๋Š” ์ œ1์ข… ๋ณ€ํ˜• ๋ฒ ์…€ ํ•จ์ˆ˜ \(I_v(x)\)์™€ ์ œ2์ข… ๋ณ€ํ˜• ๋ฒ ์…€ ํ•จ์ˆ˜ \(K_v(x)\), ๊ทธ๋ฆฌ๊ณ  ์ด๋“ค์˜ 1์ฐจ ๋„ํ•จ์ˆ˜ \(I'_v(x)\)์™€ \(K'_v(x)\)๋ฅผ ๊ณ„์‚ฐํ•ฉ๋‹ˆ๋‹ค. ๋‘ ํ•จ์ˆ˜๋Š” ๋ณ€ํ˜• ๋ฒ ์…€ ๋ฐฉ์ •์‹ \(x^2 y'' + x y' - (x^2 + v^2)y = 0\)์˜ ์„œ๋กœ ๋…๋ฆฝ์ ์ธ ๋‘ ํ•ด์ž…๋‹ˆ๋‹ค. ์›ํ†ต์—์„œ์˜ ์—ด์ „๋„, ํ™•์‚ฐ, ์ „์†ก์„ ยท๋„ํŒŒ๊ด€ ์ด๋ก , ํ†ต๊ณ„ํ•™ ๋“ฑ ๋ฌผ๋ฆฌํ•™๊ณผ ๊ณตํ•™ ์ „๋ฐ˜์— ๊ฑธ์ณ ๋“ฑ์žฅํ•ฉ๋‹ˆ๋‹ค. ์ˆœ์ˆ˜ ์ˆ˜ํ•™์ด๋ฏ€๋กœ ๊ฒฐ๊ณผ๋Š” ์–ด๋””์„œ๋‚˜ ๋™์ผํ•˜๋ฉฐ, ์ง€์—ญ๋ณ„ ๊ทœ์น™์€ ์ ์šฉ๋˜์ง€ ์•Š์Šต๋‹ˆ๋‹ค.

x์— ๋Œ€ํ•œ ์ œ1์ข… ๋ฐ ์ œ2์ข… ๋ณ€ํ˜• ๋ฒ ์…€ ํ•จ์ˆ˜๋ฅผ ๋ณด์—ฌ์ฃผ๋Š” ๋‘ ๊ณก์„ 
x๊ฐ€ ์ฆ๊ฐ€ํ•˜๋ฉด \(I_v(x)\)๋Š” ์ง€์ˆ˜์ ์œผ๋กœ ์ฆ๊ฐ€ํ•˜๊ณ  \(K_v(x)\)๋Š” 0์œผ๋กœ ๊ฐ์†Œํ•ฉ๋‹ˆ๋‹ค.

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

์ฐจ์ˆ˜ v(์ž„์˜์˜ ์‹ค์ˆ˜)์™€ ์ธ์ˆ˜ x๋ฅผ ์ž…๋ ฅํ•˜์„ธ์š”. \(I_v(x)\)๋Š” v๊ฐ€ ์ •์ˆ˜์ผ ๋•Œ ๋ชจ๋“  ์‹ค์ˆ˜ x์— ๋Œ€ํ•ด ๊ณ„์‚ฐ๋˜๋ฉฐ, ๊ทธ ์™ธ์˜ ๊ฒฝ์šฐ์—๋Š” \(x \ge 0\)์—์„œ ๊ณ„์‚ฐ๋ฉ๋‹ˆ๋‹ค. \(K_v(x)\)๋Š” \(x \to 0^+\)์—์„œ ๋ฐœ์‚ฐํ•˜๊ธฐ ๋•Œ๋ฌธ์— \(x > 0\)์ด์–ด์•ผ ํ•˜๋ฉฐ, \(x \le 0\)์ธ ๊ฒฝ์šฐ์—๋Š” ์ •์˜๋˜์ง€ ์•Š์Œ(NaN)์œผ๋กœ ํ‘œ์‹œ๋ฉ๋‹ˆ๋‹ค.

๊ณต์‹ ์„ค๋ช…

\(I_v(x)\)๋Š” ๊ฑฐ๋“ญ์ œ๊ณฑ ๊ธ‰์ˆ˜๋กœ๋ถ€ํ„ฐ ํ•ฉ์‚ฐํ•˜๋ฉฐ, ๊ฐ๋งˆ ํ•จ์ˆ˜์—๋Š” ๋ž€์ดˆ์Šค(Lanczos) ๊ทผ์‚ฌ๋ฅผ ์‚ฌ์šฉํ•ฉ๋‹ˆ๋‹ค.

$$I_v(x) = \sum_{k=0}^{\infty} \frac{1}{k!\,\Gamma(k + v + 1)} \left(\frac{x}{2}\right)^{2k + v}$$

\(K_v(x)\)๋Š” $$K_v(x) = \int_{0}^{\infty} e^{-x\cdot\cosh t}\cdot\cosh(vt)\,dt$$๋ฅผ ์ˆ˜์น˜ ์ ๋ถ„ํ•˜์—ฌ ๊ตฌํ•˜๋Š”๋ฐ, ์ด ๋ฐฉ๋ฒ•์€ ์ •์ˆ˜ ์ฐจ์ˆ˜์™€ ๋น„์ •์ˆ˜ ์ฐจ์ˆ˜ ๋ชจ๋‘์—์„œ ์•ˆ์ •์ ์ž…๋‹ˆ๋‹ค. ๋„ํ•จ์ˆ˜๋Š” ๋Œ€์นญ ์ ํ™”์‹ $$I'_v(x) = \tfrac{1}{2}\left(I_{v-1}(x) + I_{v+1}(x)\right)$$์™€ $$K'_v(x) = -\tfrac{1}{2}\left(K_{v-1}(x) + K_{v+1}(x)\right)$$๋ฅผ ์‚ฌ์šฉํ•˜๋ฏ€๋กœ x๋กœ ๋‚˜๋ˆ„๋Š” ์—ฐ์‚ฐ์„ ํ”ผํ•  ์ˆ˜ ์žˆ์Šต๋‹ˆ๋‹ค.

๋ฌดํ•œ ๊ธ‰์ˆ˜์˜ ํ•ญ๋“ค์ด ํ•ฉ์ณ์ ธ ๋ฒ ์…€ ํ•จ์ˆ˜ ๊ฐ’์„ ์ด๋ฃจ๋Š” ๋„ํ•ด
๊ธ‰์ˆ˜์˜ ๊ฐ ํ•ญ์€ ์ ์  ์ž‘์€ ๊ธฐ์—ฌ๋ฅผ ๋”ํ•˜๋ฉฐ \(I_v(x)\)๋กœ ์ˆ˜๋ ดํ•ฉ๋‹ˆ๋‹ค.

๊ณ„์‚ฐ ์˜ˆ์‹œ (v = 0, x = 1)

๊ธ‰์ˆ˜๋กœ๋ถ€ํ„ฐ $$I_0(1) = 1 + 0.25 + 0.015625 + \ldots \approx 1.26606588$$์„ ์–ป์Šต๋‹ˆ๋‹ค. ์ ๋ถ„์œผ๋กœ๋Š” \(K_0(1) \approx 0.42102444\)๋ฅผ ์–ป์Šต๋‹ˆ๋‹ค. \(I_{-1} = I_1\)์ด๋ฏ€๋กœ ๋Œ€์นญ ํ˜•ํƒœ์—์„œ \(I'_0(1) = I_1(1) \approx 0.56515910\), \(K'_0(1) = -K_1(1) \approx -0.60190723\)์ด ๋ฉ๋‹ˆ๋‹ค.

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

\(K_v(x)\)๊ฐ€ ์™œ ์ •์˜๋˜์ง€ ์•Š์Œ์œผ๋กœ ํ‘œ์‹œ๋˜๋‚˜์š”? \(K_v(x)\)๋Š” \(x > 0\)์—์„œ๋งŒ ์ •์˜๋ฉ๋‹ˆ๋‹ค. 0 ์ดํ•˜์—์„œ๋Š” ๋ฐœ์‚ฐํ•˜๊ธฐ ๋•Œ๋ฌธ์ž…๋‹ˆ๋‹ค.

๋ถ„์ˆ˜ ์ฐจ์ˆ˜๋ฅผ ์‚ฌ์šฉํ•  ์ˆ˜ ์žˆ๋‚˜์š”? ๋„ค. ๋‘ ํ•จ์ˆ˜ ๋ชจ๋‘ ๋น„์ •์ˆ˜์™€ ์Œ์ˆ˜๋ฅผ ํฌํ•จํ•œ ์ž„์˜์˜ ์‹ค์ˆ˜ ์ฐจ์ˆ˜๋ฅผ ๋ฐ›์Šต๋‹ˆ๋‹ค.

์ •ํ™•๋„๋Š” ์–ด๋А ์ •๋„์ธ๊ฐ€์š”? ๊ฒฐ๊ณผ๋Š” ๋ฐฐ์ •๋ฐ€๋„(์œ ํšจ์ˆซ์ž ์•ฝ 12~15์ž๋ฆฌ)๋ฅผ ์‚ฌ์šฉํ•˜๋ฉฐ, ์ ๋‹นํ•œ ๋ฒ”์œ„์˜ x์— ๋Œ€ํ•ด ํ‘œ์ค€ ์ฐธ์กฐํ‘œ์™€ ์ผ์น˜ํ•ฉ๋‹ˆ๋‹ค.

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