๋ฒ ์ Y ํจ์ ํ ๊ณ์ฐ๊ธฐ๋?
์ด ๋๊ตฌ๋ ์ 2์ข ๋ฒ ์ ํจ์๋ฅผ ํ๋ก ๊ณ์ฐํฉ๋๋ค. ๋ฒ ๋ฒ(Weber) ํจ์ ๋๋ ๋ ธ์ด๋ง(Neumann) ํจ์๋ผ๊ณ ๋ ๋ถ๋ฆฌ๋ฉฐ \(Y_{v}(x)\)๋ก ํ๊ธฐํฉ๋๋ค. ์ด ํจ์๋ ๋ฒ ์ ๋ฏธ๋ถ๋ฐฉ์ ์์ ๋ ๋ฒ์งธ ์ ํ ๋ ๋ฆฝ ํด์ ๋๋ค. ์ค์ ์ฐจ์ \(v\)๋ฅผ ๊ณ ์ ํ ๋ค, ์์๊ฐยท์ฆ๋ถยท์ ๊ฐ์๋ก ์ ์๋ ์ผ๋ จ์ \(x\) ๊ฐ์์ \(Y_{v}(x)\)๋ฅผ ๊ณ์ฐํ์ฌ ์์ ํ ์์น ํ๋ฅผ ๋ง๋ค์ด ์ค๋๋ค.
์ฌ์ฉ ๋ฐฉ๋ฒ
์ฐจ์ \(v\)(์ ์๊ฐ ์๋์ด๋ ๋๊ณ ์์๋ ๊ฐ๋ฅ), \(x\)์ ์์๊ฐ, ์ ๊ณผ ์ ์ฌ์ด์ ์ฆ๋ถ(์คํ ), ๊ทธ๋ฆฌ๊ณ ๋ฐ๋ณต ํ์(ํ ์)๋ฅผ ์ ๋ ฅํ์ธ์. ๊ณ์ฐ๊ธฐ๋ \(i = 0\)๋ถํฐ \(\text{pointCount}-1\)๊น์ง \(x_i = \text{startX} + i \cdot \text{stepX}\) ๋ฅผ ๋ง๋ค๊ณ ๊ฐ ๊ฐ์ ๋ํ \(Y_{v}(x)\)๋ฅผ ๋์ดํฉ๋๋ค. ๋จ, \(Y_{v}(x)\)๋ \(x = 0\)์์ ์์ ๋ฌดํ๋๋ก ๋ฐ์ฐํ๋ฉฐ \(x > 0\) ์์๋ง ์ค์๊ฐ์ ๊ฐ์ง๋ฏ๋ก, \(x \le 0\) ์ธ ํ์ ์ ์๋์ง ์์(undefined)์ผ๋ก ํ์๋ฉ๋๋ค.
๊ณต์
์ ์๊ฐ ์๋ ์ฐจ์์ ๊ฒฝ์ฐ:
$$Y_{\nu}(x) = \frac{J_{\nu}(x)\cos(\nu\pi) - J_{-\nu}(x)}{\sin(\nu\pi)}$$์ ์ ์ฐจ์ \(n\)์ ๊ฒฝ์ฐ์๋ ๊ทนํ์ ์ทจํ๋ฉด \(J_{n}(x)\cdot\ln(x/2)\) ํํ์ ๋ก๊ทธ ํญ, ์ ํํ ๊ฑฐ๋ญ์ ๊ณฑ ๊ธ์ ๋ณด์ ํญ, ๋๊ฐ๋ง(digamma) ๊ธ์๋ฅผ ํฌํจํ ๋ซํ ํํ์ ์์ด ์ป์ด์ง๋๋ค. ์ 1์ข ํจ์ \(J_{v}(x)\)๋ ๊ฑฐ๋ญ์ ๊ณฑ ๊ธ์๋ก๋ถํฐ ํฉ์ฐํ๋ฉฐ, ๊ฐ๋ง ํจ์๋ ๋์ด์ค(Lanczos) ๊ทผ์ฌ๋ก ๊ณ์ฐํฉ๋๋ค.
๊ณ์ฐ ์์
\(v = 0\), \(\text{startX} = 0\), \(\text{stepX} = 0.2\), \(\text{pointCount} = 51\)๋ก ์ค์ ํ๋ฉด ํ์ \(x = 0.0\)๋ถํฐ \(10.0\)๊น์ง ์งํ๋ฉ๋๋ค. \(Y_{0}(0)\)์ ์ ์๋์ง ์์(\(-\infty\)), \(Y_{0}(0.2) \approx -1.0811\), \(Y_{0}(1.0) \approx 0.0883\), \(Y_{0}(2.0) \approx 0.5104\), \(Y_{0}(10.0) \approx 0.0557\) ์ ๋๋ค. ์๋จ์ "์ฒซ ์ ํ๊ฐ"์๋ \(-1.0811\)์ด ํ์๋ฉ๋๋ค.
์ ์ ๋ฐ ์ฉ์ด์ง
- ์ฐจ์ \(\nu\)
- ๋ฒ ์ ํจ์ ์กฑ์ ์ธ๋ฑ์ฑํ๋ ๋งค๊ฐ๋ณ์(order ํ๋)์ ๋๋ค. ์์์ ์ค์๊ฐ ๋ ์ ์์ต๋๋ค. ์ ์ ์ฐจ์(0, 1, 2, โฆ)๋ ์ํต ๋์นญ์ ๊ฐ์ง ๋ฌผ๋ฆฌ ๋ฌธ์ ์์ ๊ฐ์ฅ ์ผ๋ฐ์ ์ ๋๋ค. ๋ฐ์ ์ ์ฐจ์๋ ๊ตฌ ๋ฒ ์ ํจ์๋ฅผ ์ ๊ณตํฉ๋๋ค.
- ์ 2์ข ๋ฒ ์ ํจ์ \(Y_\nu(x)\)
- ์จ๋ฒ ํจ์ ๋๋ ๋ ธ์ด๋ง ํจ์๋ผ๊ณ ๋ ํฉ๋๋ค(๋๋๋ก \(N_\nu\)๋ก ํ๊ธฐ). ์ด๊ฒ์ ์์ ์์ ๋ฌดํ(ํน์ด)์ธ ๋ฒ ์ ๋ฐฉ์ ์์ ํด์ ๋๋ค. ์ ์๊ฐ ์๋ \(\nu\)์ ๋ํด \(Y_\nu(x) = \dfrac{J_\nu(x)\cos(\nu\pi) - J_{-\nu}(x)}{\sin(\nu\pi)}\)๋ก ์ ์๋๋ฉฐ, ์ ์ ๊ฒฝ์ฐ๋ ๊ทนํ์ผ๋ก ์ป์ด์ง๋๋ค.
- \(J_\nu\) ๋ \(Y_\nu\)
- \(J_\nu(x)\)(์ 1์ข )๋ \(x=0\)์์ ์ ํํฉ๋๋ค. \(Y_\nu(x)\)(์ 2์ข )๋ \(x\to 0^+\)์ผ ๋ \(-\infty\)๋ก ๋ฐ์ฐํฉ๋๋ค. ํจ๊ป ๊ทธ๋ค์ ๋ฒ ์ ๋ฐฉ์ ์์ ์์ ํ ๋ ๋ฆฝ ํด์ ์์ ํ์ฑํฉ๋๋ค.
- ๋ฒ ์ ๋ฏธ๋ถ ๋ฐฉ์ ์
- ์ ํ ์๋ฏธ๋ถ ๋ฐฉ์ ์ \(x^2 y'' + x y' + (x^2 - \nu^2) y = 0\)์ ๋๋ค. ์ผ๋ฐํด๋ \(y = c_1 J_\nu(x) + c_2 Y_\nu(x)\)์ ๋๋ค.
- ๊ฐ๋ง ํจ์ \(\Gamma(z)\)
- ๊ณ์น์ ์ฐ์ ํ์ฅ์ผ๋ก, \(\Gamma(n+1) = n!\)์ด๋ฉฐ, \(J_\nu\)์ \(Y_\nu\)์ ๊ธ์ ๊ณ์์ ๋ํ๋ฉ๋๋ค.
- ๋๊ฐ๋ง ํจ์ \(\psi(z)\)
- ๋ก๊ทธ ๋ฏธ๋ถ \(\psi(z) = \Gamma'(z)/\Gamma(z)\)์ ๋๋ค. ์ ์ ์ฐจ์ \(Y_n(x)\)์ ๊ธ์์ ๋ช ์์ ์ผ๋ก ๋ํ๋๋ฉฐ, ๋ก๊ทธ ํญ \(\tfrac{2}{\pi}\ln(x/2)J_n(x)\)๊ณผ ๋๊ฐ๋ง ๊ฐ์ค ๊ณ์๋ฅผ ํฌํจํฉ๋๋ค.
- ๋์ด์ค ๊ทผ์ฌ
- ๋ณต์์ ๋๋ ์ค์ ์ธ์์ ๋ํด ๊ฐ๋ง ํจ์ \(\Gamma(z)\)๋ฅผ ๊ณ์ฐํ๋ ๋งค์ฐ ์ ํํ ์์น ๋ฐฉ๋ฒ์ผ๋ก, ๋ฒ ์ ํจ์ ๋ฃจํด ๋ด์์ ๊ธ์ ๊ณ์๋ฅผ ๊ณ์ฐํ๋ ๋ฐ ์ผ๋ฐ์ ์ผ๋ก ์ฌ์ฉ๋ฉ๋๋ค.
- ์ ํ๋ ๋ฆฝ ํด
- ์ฒซ ๋ฒ์งธ ํด์ ์์๋ฐฐ๋ก ํํํ ์ ์๋ ๋ ๋ฒ์งธ ํด์ ๋๋ค. \(J_\nu\)๋ง์ผ๋ก๋ ์์ ์์ ํน์ด์ธ ํด๋ฅผ ๋ํ๋ผ ์ ์๊ธฐ ๋๋ฌธ์, \(Y_\nu\)๋ ์ผ๋ฐํด์ ํ์ํ ๋ ๋ฆฝ์ ์ธ ๋๋ฐ์๋ฅผ ์ ๊ณตํฉ๋๋ค.
์์ฃผ ๋ฌป๋ ์ง๋ฌธ
์ ์ฒซ ๋ฒ์งธ ํ์ด ์ ์๋์ง ์์์ผ๋ก ๋์ค๋์? \(Y_{v}(x)\)๋ \(x = 0\)์์ ํน์ด์ ์ ๊ฐ์ง๋ฉฐ \(-\infty\)๋ก ๋ฐ์ฐํ๊ธฐ ๋๋ฌธ์ ๊ทธ ์ง์ ์์๋ ์ ํํ ๊ฐ์ด ์กด์ฌํ์ง ์์ต๋๋ค.
์ฐจ์๋ฅผ ์์๋ก ๋ ์ ์๋์? ๋ค. ์์ ์ ์ ์ฐจ์์๋ ๋์นญ์ฑ \(Y_{-n}(x) = (-1)^{n}Y_{n}(x)\) ๊ฐ ์ ์ฉ๋๋ฉฐ, ์์ ๋น์ ์ ์ฐจ์์๋ ์ผ๋ฐ ๊ณต์์ ๊ทธ๋๋ก ์ฌ์ฉํฉ๋๋ค.
์ ํ๋๋ ์ด๋ ์ ๋์ธ๊ฐ์? ๊ธ์๋ ๊ฐ ํญ์ด ๊ธฐ๊ณ ํ์ฉ์ค์ฐจ ์๋๋ก ๋จ์ด์ง ๋๊น์ง ํฉ์ฐ๋๋ฉฐ, ์ ๋นํ ํฌ๊ธฐ์ \(x\)์ ๋ํด ์ฝ 6~7์๋ฆฌ์ ์ ํจ์ซ์๋ฅผ ์ ๊ณตํฉ๋๋ค.