์ด ๊ณ์ฐ๊ธฐ์ ๊ธฐ๋ฅ
์ด ๋๊ตฌ๋ ๋ฒ ์ ๋ฏธ๋ถ๋ฐฉ์ ์ \(x^2y'' + xy' + (x^2 - v^2)y = 0\)์ ์ ํ ๋ ๋ฆฝ์ธ ๋ ํด, ์ฆ ์ 1์ข ๋ฒ ์ ํจ์ \(J_v(x)\), ์ 2์ข ๋ฒ ์ ํจ์ \(Y_v(x)\)(๋ ธ์ด๋ง ํจ์๋ผ๊ณ ๋ ํจ), ๊ทธ๋ฆฌ๊ณ ์ด๋ค์ 1์ฐจ ๋ํจ์ \(J'_v(x)\)์ \(Y'_v(x)\)๋ฅผ ๊ณ์ฐํฉ๋๋ค. ์ฐจ์ \(v\)๋ ์ ์, ๋ถ์, ์์ ๋ฑ ์ด๋ค ์ค์๋ ์ ๋ ฅํ ์ ์์ผ๋ฉฐ ์ธ์ \(x\)๋ ์ค์๋ก ๋ฐ์ต๋๋ค. ๋ฒ ์ ํจ์๋ ๋ฌผ๋ฆฌํ๊ณผ ๊ณตํ ์ ๋ฐ์ ๋ฑ์ฅํฉ๋๋ค. ์ํ ๋ง์ ์ง๋, ์ํต์์์ ์ด์ ๋, ๋ํ๊ด ๋ด ์ ์๊ธฐํ, ์ ํธ ์ฒ๋ฆฌ ๋ฑ์ด ๋ํ์ ์ธ ์์ ๋๋ค.
์ฌ์ฉ ๋ฐฉ๋ฒ
์ฐจ์ \(v\)(์: 0, 1, 0.5), ์ธ์ \(x\)๋ฅผ ์ ๋ ฅํ๊ณ ํ์ํ ์๋ฆฟ์๋ฅผ ์ ํํ์ธ์. ๊ณ์ฐ ๋ฒํผ์ ๋๋ฅด๋ฉด ๋ค ๊ฐ์ง ๊ฐ์ด ๋ชจ๋ ๋ํ๋ฉ๋๋ค. ์ฐจ์๊ฐ ์ ์๊ฐ ์๋ ๊ฒฝ์ฐ์๋ \(x \geq 0\)์ผ๋ก ์ ๋ ฅํ์ธ์. \(x\)๊ฐ ์์์ด๋ฉด \((x/2)^v\)๊ฐ ๋ณต์์๊ฐ ๋๊ธฐ ๋๋ฌธ์ ๋๋ค. \(x = 0\)์์๋ ์ 2์ข ํจ์๊ฐ ํน์ด์ ์ ๊ฐ์ง๋ฏ๋ก ์ ์๋์ง ์์์ผ๋ก ํ์๋ฉ๋๋ค.
๊ณต์ ์ค๋ช
\(J_v(x)\)๋ ๊ฐ๋ง ํจ์๋ฅผ ์ด์ฉํ ๋ฉฑ๊ธ์๋ก ๊ณ์ฐํ๋ฉฐ, ์์ ์ ์ธ ์ ํ์์ ์ฌ์ฉํด ๊ฐ ํญ์ด ํ์ฉ ์ค์ฐจ ์ดํ๋ก ์์์ง ๋๊น์ง ํญ๋ณ๋ก ํฉ์ฐํฉ๋๋ค.
$$J_{\nu}(x) = \sum_{k=0}^{\infty} \frac{(-1)^{k}}{k!\,\Gamma(\nu+k+1)} \left(\frac{x}{2}\right)^{2k+\nu}, \qquad \nu = \text{Order }\nu,\; x = \text{Argument }x$$์ ์๊ฐ ์๋ ์ฐจ์์ \(Y_v(x)\)๋ ๋ค์ ๊ณต์์ ์ฌ์ฉํฉ๋๋ค.
$$Y_{\nu}(x) = \frac{J_{\nu}(x)\cos(\nu\pi) - J_{-\nu}(x)}{\sin(\nu\pi)}, \qquad \nu = \text{Order }\nu,\; x = \text{Argument }x$$์ ์ ์ฐจ์์ผ ๋๋ \(\sin(v\pi)=0\)์ผ๋ก ์ธํ 0์ผ๋ก ๋๋๊ธฐ๋ฅผ ํผํ๊ธฐ ์ํด \(v\)๋ฅผ ์์ฃผ ์ฝ๊ฐ(\(10^{-7}\)๋งํผ) ๋ฏธ์ธํ๊ฒ ๋ณํํฉ๋๋ค. ๋ํจ์๋ ์ ํ์์ ์ฌ์ฉํ๋ฉฐ, ํน์ํ ๊ฒฝ์ฐ๋ก \(C'_0(x) = -C_1(x)\)๋ฅผ ์ ์ฉํฉ๋๋ค.
$$C_{\nu}^{\prime}(x) = \tfrac{1}{2}\bigl(C_{\nu-1}(x) - C_{\nu+1}(x)\bigr), \quad C \in \{J,\,Y\}, \quad \nu = \text{Order }\nu,\; x = \text{Argument }x$$
๊ณ์ฐ ์์
\(v = 0\), \(x = 1\)์ธ ๊ฒฝ์ฐ: \(J_0(1)\)์ ๊ธ์๋ ๋ค์๊ณผ ๊ฐ์ต๋๋ค.
$$J_0(1) = 1 - 0.25 + 0.015625 - 0.000434 + \ldots \approx 0.7651977$$์๋ ค์ง ๊ฐ๋ค์ \(Y_0(1) \approx 0.0882570\), \(J'_0(1) = -J_1(1) \approx -0.4400506\), \(Y'_0(1) = -Y_1(1) \approx 0.7812128\)์ ๋๋ค.
์ฃผ์ ์ฉ์ด ๋ฐ ๋ณ์
- ์ฐจ์ \(\nu\)
- ๋ฏธ๋ถ๋ฐฉ์ ์ \(x^2 y'' + x y' + (x^2-\nu^2)y = 0\)์ ํํ๋ฅผ ๊ฒฐ์ ํ๋ \(J_\nu(x)\) ๋ฐ \(Y_\nu(x)\)์ ๋งค๊ฐ๋ณ์ \(\nu\)์ ๋๋ค. ์์์ ์ค์์ผ ์ ์์ผ๋ฉฐ, ์ ์ ์ฐจ์ (\(\nu = 0,1,2,\dots\))๋ ์ํต ์ขํ์์์ ๊ฐ๋ ๋ถ๋ฆฌ๋ก๋ถํฐ ๋ํ๋๊ณ , ๋ฐ์ ์ ์ฐจ์๋ ์ด๋ฑํจ์๋ก ํํ ๊ฐ๋ฅํ ๊ตฌ๋ฉด ๋ฒ ์ ํจ์๋ฅผ ์ ๊ณตํฉ๋๋ค.
- ๋ ๋ฆฝ๋ณ์ \(x\)
- ํจ์๊ฐ ํ๊ฐ๋๋ ๋ ๋ฆฝ๋ณ์๋ก, ์ผ๋ฐ์ ์ผ๋ก ์ค์ผ์ผ๋ ๋ฐ์ง๋ฆ ๊ฑฐ๋ฆฌ \(x = kr\)์ ๋๋ค. ์ค์ \(x\)์ ๋ํด, \(J_\nu\)๋ ์ ์ \(\nu\)์ผ ๋ ์ค์๊ฐ์ด๊ณ , \(Y_\nu\)๋ \(x>0\)์์๋ง ์ ์๋ฉ๋๋ค.
- ์ 1์ข ๋ฒ ์ ํจ์ \(J_\nu(x)\)
- ์์ ์์ ์ ํํ ํด (\(\nu\ge 0\)์ ๊ฒฝ์ฐ)๋ก, ๊ธ์ \(J_\nu(x)=\sum_{k=0}^{\infty}\frac{(-1)^k}{k!\,\Gamma(\nu+k+1)}\left(\tfrac{x}{2}\right)^{2k+\nu}\)๋ก ์ ์๋ฉ๋๋ค. \(x\)๊ฐ ์ฆ๊ฐํจ์ ๋ฐ๋ผ ์งํญ์ด ์ฒ์ฒํ ๊ฐ์ํ๋ฉฐ ์ง๋ํฉ๋๋ค.
- ์ 2์ข ๋ฒ ์ ํจ์ \(Y_\nu(x)\)
- ๋ ธ์ด๋ง(๋๋ ๋ฒ ๋ฒ) ํจ์๋ผ๊ณ ๋ ๋ถ๋ฆฌ๋ฉฐ, ์ด๋ ๋ ๋ฒ์งธ ์ ํ๋ ๋ฆฝ ํด์ ๋๋ค. \(Y_\nu(x)=\dfrac{J_\nu(x)\cos(\nu\pi)-J_{-\nu}(x)}{\sin(\nu\pi)}\)๋ฅผ ํตํด ์ ์๋๊ณ (์ ์ \(\nu\)์ ๊ฒฝ์ฐ ๊ทนํํ), ์์ ์์ ๋ก๊ทธ์ ์ผ๋ก ๋๋ \(x\)์ ๊ฑฐ๋ญ์ ๊ณฑ์ผ๋ก ๋ฐ์ฐํฉ๋๋ค.
- ๋ํจ์ \(J'_\nu(x)\), \(Y'_\nu(x)\)
- \(x\)์ ๋ํ ๋ํจ์์ ๋๋ค. ์ด๋ค์ ์ ํ์ \(C'_\nu(x)=\tfrac{1}{2}\bigl(C_{\nu-1}(x)-C_{\nu+1}(x)\bigr)\) ๋ฐ \(C'_\nu(x)=C_{\nu-1}(x)-\tfrac{\nu}{x}C_\nu(x)\)๋ฅผ ๋ง์กฑํฉ๋๋ค. ์ฌ๊ธฐ์ \(C\)๋ \(J\) ๋๋ \(Y\)๋ฅผ ๋ํ๋ ๋๋ค. ํนํ \(J'_0(x)=-J_1(x)\)์ ๋๋ค.
- ๊ฐ๋ง ํจ์ \(\Gamma(z)\)
- ๊ณ์น์ ์ฐ์ ํ์ฅ์ผ๋ก, ์์ด ์๋ ์ ์์ ๋ํด \(\Gamma(n+1)=n!\)์ด๋ฉฐ, ๋น์ ์ ์ฐจ์๋ฅผ ํ์ฉํ๊ธฐ ์ํด \(J_\nu\) ๊ธ์์ ๋ถ๋ชจ์ ๋ํ๋ฉ๋๋ค. ๊ฐ๋ง ํจ์ ๊ณ์ฐ๊ธฐ์์ ๊ฐ๋ณ ๊ฐ์ ํ์ธํ์ธ์.
- ์์ (๊ทผ)
- \(J_\nu(x)=0\) ๋๋ \(Y_\nu(x)=0\)์ธ ๊ฐ \(j_{\nu,m}\) ๋ฐ \(y_{\nu,m}\)์ ๋๋ค. ์ด๋ค์ ๊ฒฝ๊ณ๊ฐ ๋ฌธ์ ์์ ๊ณ ์ ๊ฐ์ผ๋ก ์์ฉํ๋ฉฐ, ์๋ฅผ ๋ค์ด ๊ณ ์ ๋ ๊ฐ์ฅ์๋ฆฌ ์ํ ๋ง์ ์ง๋ ์ฃผํ์๋ ์์ \(j_{\nu,m}\)์ ๋น๋กํฉ๋๋ค.
์์ฃผ ๋ฌป๋ ์ง๋ฌธ
์ฐจ์๊ฐ ์์๋ ๋ถ์์ฌ๋ ๋๋์? ๋ค. ๊ธ์์ ๊ฐ๋ง ํจ์๊ฐ ๋ชจ๋ ์ค์ \(v\)๋ฅผ ์ฒ๋ฆฌํ ์ ์์ต๋๋ค. ๋ค๋ง ์ ์๊ฐ ์๋ \(v\)์์๋ \(x \geq 0\)์ ์ ์งํ์ธ์.
\(x = 0\)์์ Y๊ฐ ์ ์๋์ง ์๋ ์ด์ ๋? ๋ชจ๋ ์ 2์ข ๋ฒ ์ ํจ์๋ \(x \to 0\)์ผ ๋ ์์ ๋ฌดํ๋๋ก ๋ฐ์ฐํ๋ฏ๋ก ์ ํํ ๊ฐ์ด ์กด์ฌํ์ง ์์ต๋๋ค.
์ ํ๋๋ ์ด๋ ์ ๋์ธ๊ฐ์? ๊ณ์ฐ์ ๋ฐฐ์ ๋ฐ๋(์ฝ 15์๋ฆฌ์ ์ ํจ ์ซ์)๋ก ์ํ๋ฉ๋๋ค. ํ์ ์๋ฆฟ์ ์ต์ ์ ํ์ ํ์๋ง ์กฐ์ ํ ๋ฟ, ๋ด๋ถ ๊ณ์ฐ์๋ ์ํฅ์ ์ฃผ์ง ์์ต๋๋ค.