์ด ๊ณ์ฐ๊ธฐ์ ๊ธฐ๋ฅ
์ด ๋๊ตฌ๋ ์ 1์ข ๋ฒ ์ ํจ์ \(J_{v}(x)\)๋ฅผ ๊ณ ์ ๋ ์ฐจ์ \(v\)์ ๋ํด ๋ณ์ \(x\)๋ฅผ ๋ณํ์ํค๋ฉฐ ํ๋ก ์ ๋ฆฌํฉ๋๋ค. \(x\)์ ์์๊ฐ, ์ฆ๊ฐ๋, ์์ฑํ ํ์ ๊ฐ์๋ฅผ ์ ๋ ฅํ๋ฉด \(x\)์ \(J_{v}(x)\)๋ฅผ ๋๋ํ ๋ณด์ฌ ์ฃผ๋ ๊น๋ํ 2์ด ํ๋ฅผ ์ถ๋ ฅํฉ๋๋ค. ์ 1์ข ๋ฒ ์ ํจ์๋ ๋ฌผ๋ฆฌํ๊ณผ ๊ณตํ ์ ๋ฐ์ ๋ฑ์ฅํฉ๋๋ค. ์ํ ๋ถ์ ์ง๋, ์ํต ๋ด๋ถ์ ์ด์ ๋, ๋ํ๊ด ์ ์ ์๊ธฐํ, ์ ํธ ์ฒ๋ฆฌ(FM ๋ณ์กฐ์ ์ธกํ๋) ๋ฑ์์ ํญ๋๊ฒ ์ฐ์ ๋๋ค.
์ฌ์ฉ ๋ฐฉ๋ฒ
์ฐจ์ \(v\)๋ฅผ ์ ๋ ฅํฉ๋๋ค(0, 1, 2 ๊ฐ์ ์ ์๋ ๋ฌผ๋ก 0.5 ๊ฐ์ ๋ถ์๋ ์์ ๋ฑ ์์์ ์ค์๊ฐ ๊ฐ๋ฅํฉ๋๋ค). x์ ์ด๊ธฐ๊ฐ, ์ฆ๊ฐ๋(์ฐ์๋ \(x\) ๊ฐ ์ฌ์ด์ ๊ฐ๊ฒฉ์ผ๋ก, ๋ด๋ฆผ์ฐจ์ ๊ณ์ฐ์ ์ํ๋ฉด ์์, ํ ์ ์ ๋ฐ๋ณตํ๋ ค๋ฉด 0๋ ๊ฐ๋ฅ), ๊ทธ๋ฆฌ๊ณ ๋ฐ๋ณต ํ์(์์ฑํ ํ ์๋ก 1๋ถํฐ 10000๊น์ง)๋ฅผ ์ง์ ํ์ธ์. \(i\)๋ฒ์งธ ํ์ \(x = \text{startX} + i \times \text{stepX}\) ๋ก ๊ณ์ฐ๋ฉ๋๋ค.
๊ณต์ ์ค๋ช
์ด ํจ์๋ ๋ฉฑ๊ธ์ $$J_{v}(x) = \sum_{k=0}^{\infty} \frac{(-1)^{k}}{k!\,\Gamma(k+v+1)} \left(\frac{x}{2}\right)^{v+2k}$$ ๋ก ์ ์๋๋ฉฐ, ์ฌ๊ธฐ์ \(\Gamma\)๋ ๊ฐ๋ง ํจ์์ ๋๋ค. ๊ณ์ฐ๊ธฐ๋ ์ด ๊ธ์๋ฅผ ํญ๋ณ๋ก ๊ณ์ฐํ๋ ์์ ์ ์ธ ์ ํ์์ ์ฌ์ฉํฉ๋๋ค. ์ฆ, ๊ฐ ํญ์ ์ด์ ํญ์ \(-\dfrac{x^{2}/4}{(k+1)(k+v+1)}\) ์ ๊ณฑํด ์ป์ผ๋ฏ๋ก ํฉํ ๋ฆฌ์ผ ์ค๋ฒํ๋ก๋ฅผ ํผํ ์ ์์ต๋๋ค. ๊ฐ๋ง ํจ์๋ ๋์ด์ค(Lanczos) ๊ทผ์ฌ๋ก ๊ณ์ฐํ์ฌ ์ ์๊ฐ ์๋ ์ฐจ์๋ ์์ ์ฐจ์๋ ์ฒ๋ฆฌํ ์ ์์ต๋๋ค. ์์ ์ ์ ์ฐจ์์ผ ๋๋ ํญ๋ฑ์ \(J_{-n}(x) = (-1)^{n} J_{n}(x)\) ๋ฅผ ์ฌ์ฉํฉ๋๋ค.
๊ณ์ฐ ์์
\(v = 0\), \(\text{startX} = 0\), \(\text{stepX} = 0.2\), \(\text{loopCount} = 6\) ์ผ๋ก ์ค์ ํ๋ฉด ํ๋ \(J_{0}(0) = 1\), \(J_{0}(0.2) \approx 0.990025\), \(J_{0}(0.4) \approx 0.960398\), \(J_{0}(0.6) \approx 0.912005\), \(J_{0}(0.8) \approx 0.846287\), \(J_{0}(1.0) \approx 0.765198\) ์ ์ถ๋ ฅํฉ๋๋ค. ์ด๋ ํ์ค ์์นํ์ ๊ฐ \(J_{0}(1) = 0.7651976866\) ๊ณผ ์ผ์นํฉ๋๋ค.
J_v(x)์ ์ฐธ์กฐ๊ฐ
๋ค์ ํ๋ ์ฌ๋ฌ ํ์ค ์ธ์์์ 1์ข ๋ฒ ์ ํจ์ \(J_v(x)\)๋ฅผ ์ฐจ์ \(v=0,1,2\)์ ๋ํด ๋์ดํฉ๋๋ค. ๊ฐ์ ์์์ ์ดํ 6์๋ฆฌ๋ก ๋ฐ์ฌ๋ฆผ๋์์ผ๋ฉฐ ๊ธ์ \(J_{v}(x)=\sum_{k=0}^{\infty}\frac{(-1)^{k}}{k!\,\Gamma(v+k+1)}\left(\frac{x}{2}\right)^{2k+v}\)์์ ๋์ต๋๋ค.
| \(x\) | \(J_0(x)\) | \(J_1(x)\) | \(J_2(x)\) |
|---|---|---|---|
| 0 | 1.000000 | 0.000000 | 0.000000 |
| 0.5 | 0.938470 | 0.242268 | 0.030604 |
| 1 | 0.765198 | 0.440051 | 0.114903 |
| 2 | 0.223891 | 0.576725 | 0.352834 |
| 3 | โ0.260052 | 0.339059 | 0.486091 |
| 5 | โ0.177597 | โ0.327579 | 0.046565 |
| 10 | โ0.245936 | 0.043473 | 0.254630 |
ํ์ธ ์์ ๋ก์ \(J_0(1)\)์ ๊ณ์ฐํ๋ฉด: ์ฒ์ ๋ช ํญ์ ๋ํ๋ฉด \(1-\tfrac{(0.5)^2}{1}+\tfrac{(0.5)^4}{4}-\tfrac{(0.5)^6}{36}+\dots = 1-0.25+0.015625-0.000434+\dots\approx\) 0.765198.
์ฃผ์ ์์ (๊ทผ)
์์ ์์ ์ \(J_v(x)=0\)์ธ \(x\)์ ๊ฐ์ ๋๋ค. ์ด๋ค์ ๋๋ผ์ ์ง๋ ๋ชจ๋, ํ์ดํ ์ฐจ๋จ ์ฃผํ์ ๋ฐ ์ ์ฌํ ๊ฒฝ๊ณ ์กฐ๊ฑด์ ๊ฒฐ์ ํฉ๋๋ค.
| ๊ทผ ์ง์ \(s\) | \(J_0\)์ \(s\)๋ฒ์งธ ์์ | \(J_1\)์ \(s\)๋ฒ์งธ ์์ |
|---|---|---|
| 1 | 2.404826 | 3.831706 |
| 2 | 5.520078 | 7.015587 |
| 3 | 8.653728 | 10.173468 |
| 4 | 11.791534 | 13.323692 |
\(x=0\)์ ๋ชจ๋ ์ฐจ์ \(v>0\)์์ \(J_v\)์ ์์ ์ด์ง๋ง, ์์ ์์ ๊ทผ์๋ ํฌํจ๋์ง ์์ต๋๋ค.
์ ์ & ์ฉ์ด์ง
- ์ฐจ์ \(v\)
- ๋ฒ ์ ํจ์์กฑ์ ์ด๋ ๋ฉค๋ฒ๋ฅผ ๊ณ์ฐํ ์ง ์ ํํ๋ ๋งค๊ฐ๋ณ์์ ๋๋ค(์ฌ๊ธฐ์ ์์ ํ๋ order). ์์์ ์ค์์ผ ์ ์์ต๋๋ค โ ์ ์ ์ฐจ์๋ ์ํ ๋ฌธ์ ์์ ๋ํ๋๊ณ , ๋ฐ์ ์ ์ฐจ์ \(v=n+\tfrac12\)๋ ๊ตฌ๋ฉด ๋ฒ ์ ํจ์๋ฅผ ์ ๊ณตํฉ๋๋ค.
- ์ธ์ \(x\)
- \(J_v\)๊ฐ ๊ณ์ฐ๋๋ ๋ ๋ฆฝ ๋ณ์์ ๋๋ค. ์ด ํ์์๋ startX์์ ์์ํ์ฌ loopCountํ์ ๊ฑธ์ณ stepX์ฉ ์ฆ๊ฐํฉ๋๋ค.
- ๊ฐ๋ง ํจ์ \(\Gamma\)
- ํฉํ ๋ฆฌ์ผ์ ์ฐ์ ํ์ฅ์ผ๋ก, ์์ด ์๋ ์ ์์ ๋ํด \(\Gamma(n+1)=n!\)์ ๋๋ค. ์ด๋ ๊ธ์์ ๋ถ๋ชจ \(\Gamma(v+k+1)\)์ ๋ํ๋๋ฏ๋ก ๋น์ ์ ์ฐจ์๋ ์ ์ ์๋ฉ๋๋ค.
- 1์ข ๋ฒ ์ ํจ์ \(J_v(x)\)
- ๋ฒ ์ ์ ๋ฏธ๋ถ๋ฐฉ์ ์ \(x^2 y''+x y'+(x^2-v^2)y=0\)์ ํด๋ก์, ์์ ์์ ์ ํํ๊ฒ ์ ์ง๋ฉ๋๋ค(\(v\ge 0\)์ธ ๊ฒฝ์ฐ). ์ด๋ ์์ ๊ณต์์ ์๋ ๋ฉฑ๊ธ์๋ก ์ฃผ์ด์ง๋๋ค.
- ์์ /๊ทผ
- \(J_v(x)=0\)์ธ \(x\)์ ๊ฐ์ ๋๋ค. ๊ฐ ์ฐจ์๋ ๋ฌดํํ ๋ง์ ์์ ์์ ์ ๊ฐ์ง๋ฉฐ, ์ ์ ๋ ๊ท ๋ฑํ ๊ฐ๊ฒฉ์ผ๋ก ๋ฐฐ์น๋๊ณ ์ ๊ทผ์ ์ผ๋ก \(\pi\)์ฉ ๋ถ๋ฆฌ๋ฉ๋๋ค.
- ๋ฐ์ ์(๊ตฌ๋ฉด) ์ฐจ์
- \(v=n+\tfrac12\)์ผ ๋, \(J_v\)๋ ๊ตฌ๋ฉด ๋ฒ ์ ํจ์ \(j_n(x)=\sqrt{\tfrac{\pi}{2x}}\,J_{n+1/2}(x)\)์ ๊ด๋ จ์ด ์์ผ๋ฉฐ, ์ด๋ ๊ตฌ๋ฉด ์ขํ์์ ํ๋ ๋ฐฉ์ ์์ ๋ฐ์ง๋ฆ ๋ถ๋ถ์ ์ค๋ช ํฉ๋๋ค.
- ์ฌ๊ท ํญ์ ๋น
- ๊ธ์์ ์ฐ์ ํญ์ \(\frac{a_{k+1}}{a_k}=\frac{-(x/2)^2}{(k+1)(v+k+1)}\)๋ฅผ ๋ง์กฑํฉ๋๋ค. ์ด๋ ์ด์ ํญ์์ ๊ฐ ํญ์ ์์ฑํ๊ณ ์๋ ด์ฑ์ ํ๊ฐํ๊ธฐ ์ํด ๋ด๋ถ์ ์ผ๋ก ์ฌ์ฉ๋ฉ๋๋ค.
ํ ํด์ํ๊ธฐ
์์ฑ๋ ํ์ ์ด์ ์ฝ๋ ๋ฐ ๋์์ด ๋๋ ๋ช ๊ฐ์ง ์ฌ์ค์ด ์์ต๋๋ค:
- ์์๊ฐ. \(J_0(0)=1\)์ด์ง๋ง, ๋ชจ๋ ์ฐจ์ \(v>0\)์ ๋ํด \(J_v(0)=0\)์ ๋๋ค. ๋ฐ๋ผ์ \(x=0\)์์ ์์ํ๋ ํ๋ 0์ฐจ ๊ฒฝ์ฐ์๋ง 1์์ ์์ํฉ๋๋ค.
- ๊ฐ์ ๋ฅผ ์๋ฐํ ์ง๋. ํฐ \(x\)์ ๋ํด, \(J_v(x)\approx\sqrt{\tfrac{2}{\pi x}}\cos\!\left(x-\tfrac{v\pi}{2}-\tfrac{\pi}{4}\right)\)์ ๋๋ค. ํจ์๋ ์์์ด ์ด๋ํ ์ฝ์ฌ์ธ์ฒ๋ผ ์ง๋ํ๋ฉด์ ์งํญ์ \(1/\sqrt{x}\)๋ก ๊ฐ์ ํฉ๋๋ค. ๋ฐ๋ผ์ ์ฐ์ ๊ทน๋๊ฐ์ \(x\)๊ฐ ์ฆ๊ฐํจ์ ๋ฐ๋ผ ์ฒ์ฒํ ์ถ์๋ฉ๋๋ค.
- ๋ถํธ ๋ณํ๋ ์์ ์ ํ์ํฉ๋๋ค. ์ด์ด ๋ ํ ์ฌ์ด์์ ๋ถํธ๋ฅผ ๋ฐ๊พธ๋ ๊ณณ๋ง๋ค, \(J_v\)์ ๊ทผ์ด ๊ทธ ๊ตฌ๊ฐ์ ์์นํฉ๋๋ค(์: \(J_0\)์ \(x=2\)์ \(x=3\) ์ฌ์ด์์ ๋ถํธ๊ฐ ๋ฐ๋๋ฉฐ, ์ฒซ ๋ฒ์งธ ์์ \(\approx 2.4048\)์ ํฌํจํฉ๋๋ค). ํฐ ์ธ์์ ๊ฒฝ์ฐ ์ฐ์๋ ์์ ๋ค์ ์ฝ \(\pi\)๋งํผ ๋จ์ด์ ธ ์์ต๋๋ค.
- ๋ฌผ๋ฆฌ์ ๋ ธ๋. ์ด๋ฌํ ์์ ์ ๋ฌผ๋ฆฌ์ ๊ฒฝ๊ณ ์กฐ๊ฑด์ ํด๋นํฉ๋๋ค: ์ง๋ํ๋ ์ํ ๋๋ผํค๋์ ๋ฐ์ง๋ฆ ๋ชจ๋, ์ํ ํ์ดํ์ ์ฐจ๋จ ์ฃผํ์, ๊ดํ ์ฌ์ ์ ์ฅ ํจํด์ ๋ชจ๋ \(J_v\)์ ์์ ์ผ๋ก ์ธ๋ฑ์ฑ๋ฉ๋๋ค.
- ํฌ๊ธฐ. ๊ณ ์ ๋ \(x\)์ ๋ํด, ๋ ๋์ ์ฐจ์ \(v\)๋ ์์ ์ ๊ฐ๊น์ด ์์ํ์ฌ ๋ ์ฒ์ฒํ ์์นํฉ๋๋ค. ์์ \(x\)์ ๋ํด ์ ํ ๋์์ \(J_v(x)\sim \frac{1}{\Gamma(v+1)}\left(\frac{x}{2}\right)^{v}\)์ด๋ฏ๋ก, ๋ ํฐ \(v\)๋ \(x\)๊ฐ \(v\)์ ๋น๊ต ๊ฐ๋ฅํด์ง ๋๊น์ง ๋ ์๊ฒ ์ ์ง๋ฉ๋๋ค.
์ด๋ฌํ ๊ด์ฐฐ์ ์์ ํ๋ฆฝ๋ ๊ธ์ ๋ฐ ์ ๊ทผ ํํ์์ ๋ฐ๋ฅด๋ฉฐ ์ ๋ ฅํ ๋ชจ๋ ์ฐจ์์ ์ ์ฉ๋ฉ๋๋ค.
์์ฃผ ๋ฌป๋ ์ง๋ฌธ
์ฐจ์๋ฅผ ๋ถ์๋ ์์๋ก ์ง์ ํ ์ ์๋์? ๊ฐ๋ฅํฉ๋๋ค. ๊ฐ๋ง ํจ์ ๊ธฐ๋ฐ์ ๊ธ์์ด๋ฏ๋ก ๋ฐ์ ์(๊ตฌ๋ฉด ๋ฒ ์ ํจ์ ํํ๊ฐ ๋ฉ๋๋ค)๋ ์์๋ฅผ ํฌํจํ ์์์ ์ค์ ์ฐจ์๋ฅผ ์ง์ํฉ๋๋ค.
x = 0์ผ ๋๋ ์ด๋ป๊ฒ ๋๋์? \(J_{0}(0) = 1\) ์ด๊ณ , \(v > 0\) ์ธ ๊ฒฝ์ฐ \(J_{v}(0) = 0\) ์ ๋๋ค. ์ ๋์ \((x/2)^{v}\) ์ธ์๊ฐ 0์ด ๋๊ธฐ ๋๋ฌธ์ ๋๋ค.
x๊ฐ ํด ๋ ์ ํ๋๋ ์ด๋ค๊ฐ์? ๋ฐฐ์ ๋ฐ๋ ๊ธ์๋ ์ผ๋ฐ์ ์ธ ๋ฒ์(๋๋ต \(x\)๊ฐ 20โ30 ์ดํ)์์ ์ ํํฉ๋๋ค. \(x\)๊ฐ ์์ฃผ ํด ๋๋ ์๋ฆฟ์ ์๊ฑฐ ์ค์ฐจ๋ก ์ ๋ฐ๋๊ฐ ๋จ์ด์ง ์ ์๋๋ฐ, ์ด๋ฐ ๊ฒฝ์ฐ์๋ ์ ๊ทผ ํํ $$J_{v}(x) \approx \sqrt{\frac{2}{\pi x}}\,\cos\!\left(x - \frac{v\pi}{2} - \frac{\pi}{4}\right)$$ ๊ฐ ๋ ์ ํฉํฉ๋๋ค.