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

๊ณ„์‚ฐ ์ž…๋ ฅ

๊ณต์‹

๊ด‘๊ณ 

๊ฒฐ๊ณผ

ํ‘œ์˜ ์ฒซ ๋ฒˆ์งธ ๊ฐ’ โ€” x = ์ดˆ๊ธฐ๊ฐ’์—์„œ์˜ I_v(x)
1
51 rows generated
x I_v(x)
0 1
0.1 1.0025015629
0.2 1.0100250278
0.3 1.0226268794
0.4 1.0404017822
0.5 1.0634833707
0.6 1.0920453643
0.7 1.1263030183
0.8 1.1665149229
0.9 1.2129851657
1 1.2660658778
1.1 1.3261601837
1.2 1.3937255841
1.3 1.4692777979
1.4 1.5533950997
1.5 1.6467231898
1.6 1.7499806397
1.7 1.8639649621
1.8 1.9895593566
1.9 2.1277401941
2 2.2795853023
2.1 2.4462831294
2.2 2.6291428636
2.3 2.8296056006
2.4 3.049256658
2.5 3.2898391441
2.6 3.5532689042
2.7 3.8416509766
2.8 4.1572977035
2.9 4.5027486613
3 4.8807925859
3.1 5.2944914897
3.2 5.7472071872
3.3 6.2426304652
3.4 6.7848131604
3.5 7.3782034322
3.6 8.0276845471
3.7 8.7386175242
3.8 9.5168880261
3.9 10.3689579167
4 11.3019219521
4.1 12.323570116
4.2 13.4424561633
4.3 14.6679729918
4.4 16.0104355249
4.5 17.4811718556
4.6 19.0926234795
4.7 20.8584555266
4.8 22.7936779931
4.9 24.9147790758
5 27.2398718236

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

์ด ๋„๊ตฌ๋Š” ๊ณ ์ •๋œ ์‹ค์ˆ˜ ์ฐจ์ˆ˜ v์— ๋Œ€ํ•ด ์ œ1์ข… ๋ณ€ํ˜• ๋ฒ ์…€ ํ•จ์ˆ˜ \(I_{v}(x)\)๋ฅผ ์ผ๋ จ์˜ x ๊ฐ’ ๊ตฌ๊ฐ„์—์„œ ํ‘œ๋กœ ๊ณ„์‚ฐํ•ฉ๋‹ˆ๋‹ค. ์ฐจ์ˆ˜, ์‹œ์ž‘ x ๊ฐ’, ์ฆ๊ฐ€๋Ÿ‰(๊ฐ„๊ฒฉ), ์ƒ์„ฑํ•  ํ–‰ ์ˆ˜๋ฅผ ์ž…๋ ฅํ•˜๋ฉด \(x_{i} = \text{์‹œ์ž‘๊ฐ’} + i\cdot\text{๊ฐ„๊ฒฉ}\) ํ˜•ํƒœ์˜ ๋ชฉ๋ก์„ ๋งŒ๋“ค๊ณ  ๊ฐ ์ง€์ ์—์„œ \(I_{v}(x_{i})\)๋ฅผ ๊ณ„์‚ฐํ•˜์—ฌ ํ‘œ์™€ ๊ทธ๋ž˜ํ”„๋ฅผ ํ•จ๊ป˜ ์ œ๊ณตํ•ฉ๋‹ˆ๋‹ค. ์ˆœ์ˆ˜ ์ˆ˜ํ•™์˜ ํŠน์ˆ˜ํ•จ์ˆ˜ ๊ณ„์‚ฐ ๋„๊ตฌ์ด๋ฏ€๋กœ ์ง€์—ญ๋ณ„ ๊ทœ์น™์ด๋‚˜ ๋‹จ์œ„์— ๊ตฌ์• ๋ฐ›์ง€ ์•Š๊ณ  ์–ด๋””์„œ๋‚˜ ๋™์ผํ•˜๊ฒŒ ์ ์šฉ๋ฉ๋‹ˆ๋‹ค.

์—ฌ๋Ÿฌ ์ฐจ์ˆ˜์— ๋Œ€ํ•œ ์ œ1์ข… ๋ณ€ํ˜• ๋ฒ ์…€ ํ•จ์ˆ˜์˜ ๋‹จ์กฐ ์ฆ๊ฐ€ํ•˜๋Š” ์ง€์ˆ˜ ํ•จ์ˆ˜ํ˜• ๊ณก์„ ๊ตฐ
์ฐจ์ˆ˜ v = 0, 1, 2, 3์— ๋Œ€ํ•œ ๋ณ€ํ˜• ๋ฒ ์…€ ํ•จ์ˆ˜ I_v(x)๋Š” x์— ๋”ฐ๋ผ ๊ธ‰๊ฒฉํžˆ ์ฆ๊ฐ€ํ•ฉ๋‹ˆ๋‹ค.

๊ณต์‹

๋ณ€ํ˜• ๋ฒ ์…€ ํ•จ์ˆ˜ \(I_{v}(x)\)๋Š” ๋ณ€ํ˜• ๋ฒ ์…€ ๋ฐฉ์ •์‹ \(x^{2}y'' + xy' - (x^{2} + v^{2})y = 0\)์˜ ํ•ด์ž…๋‹ˆ๋‹ค. ๋ณธ ๊ณ„์‚ฐ๊ธฐ์—์„œ๋Š” ๋‹ค์Œ ๋ฉฑ๊ธ‰์ˆ˜๋กœ ๊ฐ’์„ ๊ตฌํ•ฉ๋‹ˆ๋‹ค.

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

๊ณ„์Šน(ํŒฉํ† ๋ฆฌ์–ผ)๊ณผ ๊ฐ๋งˆ ํ•จ์ˆ˜๋ฅผ ์‚ฌ์šฉํ•˜๋ฏ€๋กœ v์—๋Š” ์ž„์˜์˜ ์‹ค์ˆ˜๋ฅผ ๋„ฃ์„ ์ˆ˜ ์žˆ์Šต๋‹ˆ๋‹ค. ์ˆ˜์น˜์  ์•ˆ์ •์„ฑ์„ ์œ„ํ•ด ๊ฐ ํ•ญ์€ \(\ln\Gamma\)์˜ ๋ž€์ดˆ์Šค(Lanczos) ๊ทผ์‚ฌ๋ฅผ ์ด์šฉํ•ด ๋กœ๊ทธ ๊ณต๊ฐ„์—์„œ ๊ณ„์‚ฐํ•œ ๋’ค, ํ•ญ์˜ ํฌ๊ธฐ๊ฐ€ ๋ฌด์‹œํ•  ๋งŒํผ ์ž‘์•„์งˆ ๋•Œ๊นŒ์ง€ ํ•ฉ์‚ฐํ•ฉ๋‹ˆ๋‹ค.

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

์ฐจ์ˆ˜ v(์˜ˆ: 0, 1, 2.5), x์˜ ์ดˆ๊ธฐ๊ฐ’, ๊ฐ ํ–‰๋งˆ๋‹ค x์— ๋”ํ•ด์ง€๋Š” ์ฆ๊ฐ€๋Ÿ‰, ๊ทธ๋ฆฌ๊ณ  ๋ฐ˜๋ณต ํšŸ์ˆ˜(ํ–‰ ์ˆ˜)๋ฅผ ์ž…๋ ฅํ•˜์„ธ์š”. ๊ณ„์‚ฐ ๋ฒ„ํŠผ์„ ๋ˆ„๋ฅด๋ฉด x์™€ \(I_{v}(x)\)๋กœ ๊ตฌ์„ฑ๋œ ๋‘ ์—ด์งœ๋ฆฌ ํ‘œ์™€, ๊ฐ™์€ ๊ตฌ๊ฐ„์— ๋Œ€ํ•œ ๊ทธ๋ž˜ํ”„๊ฐ€ ๋‚˜ํƒ€๋‚ฉ๋‹ˆ๋‹ค.

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

v = 0, ์‹œ์ž‘๊ฐ’ = 0, ๊ฐ„๊ฒฉ = 0.5, ํšŸ์ˆ˜ = 5๋กœ ์„ค์ •ํ•˜๋ฉด \(x = 0, 0.5, 1, 1.5, 2\)๊ฐ€ ๋˜๊ณ  ๊ฒฐ๊ณผ๋Š” ๋‹ค์Œ๊ณผ ๊ฐ™์Šต๋‹ˆ๋‹ค.

\(I_{0}(0) = 1\), \(I_{0}(0.5) \approx 1.0634834\), \(I_{0}(1) \approx 1.2660658\), \(I_{0}(1.5) \approx 1.6467232\), \(I_{0}(2) \approx 2.2795853\). ์ด ๊ฐ’๋“ค์€ ํ‘œ์ค€ ์ฐธ์กฐํ‘œ์™€ ์ผ์น˜ํ•ฉ๋‹ˆ๋‹ค.

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

์ฐจ์ˆ˜๊ฐ€ ์Œ์ˆ˜์ด๊ฑฐ๋‚˜ ์ •์ˆ˜๊ฐ€ ์•„๋‹ˆ์–ด๋„ ๋˜๋‚˜์š”? ๋ฉ๋‹ˆ๋‹ค. ์Œ์˜ ์ •์ˆ˜ ์ฐจ์ˆ˜์—๋Š” ํ•ญ๋“ฑ์‹ \(I_{-n}(x) = I_{n}(x)\)๊ฐ€ ์ ์šฉ๋ฉ๋‹ˆ๋‹ค. ์ •์ˆ˜๊ฐ€ ์•„๋‹Œ v๋Š” \(x \geq 0\)์—์„œ ์ง€์›๋˜๋ฉฐ, \(x < 0\)์ด๋ฉด์„œ v๊ฐ€ ์ •์ˆ˜๊ฐ€ ์•„๋‹Œ ๊ฒฝ์šฐ์—๋Š” ๊ฐ’์ด ๋ณต์†Œ์ˆ˜๊ฐ€ ๋˜๋ฏ€๋กœ NaN์„ ๋ฐ˜ํ™˜ํ•ฉ๋‹ˆ๋‹ค.

\(I_{v}(x)\)๋Š” ์™œ ์ด๋ ‡๊ฒŒ ๋น ๋ฅด๊ฒŒ ์ปค์ง€๋‚˜์š”? ์ง„๋™ํ•˜๋Š” ์ผ๋ฐ˜ ๋ฒ ์…€ ํ•จ์ˆ˜ \(J_{v}\)์™€ ๋‹ฌ๋ฆฌ, ๋ณ€ํ˜• ๋ฒ ์…€ ํ•จ์ˆ˜๋Š” x๊ฐ€ ์ปค์งˆ์ˆ˜๋ก ๋Œ€๋žต \(e^{x}/\sqrt{2\pi x}\)์— ๋น„๋ก€ํ•ด ์ฆ๊ฐ€ํ•ฉ๋‹ˆ๋‹ค. ๊ทธ๋ž˜์„œ x๊ฐ€ ์•„์ฃผ ํฌ๋ฉด ๊ฐ’์ด ๋ฌดํ•œ๋Œ€๋กœ ๋ฐœ์‚ฐ(overflow)ํ•  ์ˆ˜ ์žˆ์Šต๋‹ˆ๋‹ค.

\(I_{v}(0)\)์˜ ๊ฐ’์€ ๋ฌด์—‡์ธ๊ฐ€์š”? \(I_{0}(0) = 1\)์ด๊ณ , \(v > 0\)์ธ ๊ฒฝ์šฐ \(I_{v}(0) = 0\)์ž…๋‹ˆ๋‹ค.

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