์ด ๊ณ์ฐ๊ธฐ์ ๊ธฐ๋ฅ
์ด ๋๊ตฌ๋ ์ 1์ข ๊ตฌ๋ฉด ๋ฒ ์ ํจ์ \(j_{\nu}(x)\)์ ๊ฐ์ ์ผ๋ จ์ x ๊ฐ์ ๋ํด ํ๋ก ์ ๋ฆฌํด ์ค๋๋ค. ์ฐจ์ \(\nu\), ์์ x ๊ฐ, ์ฆ๊ฐ๋(์คํ ), ์์ฑํ ํ ์๋ฅผ ์ ๋ ฅํ๋ฉด \((x, j_{\nu}(x))\) ์์ผ๋ก ์ด๋ฃจ์ด์ง ๋ ์ด์ง๋ฆฌ ํ๋ฅผ ๋๋ ค์ค๋๋ค. ์์ํ ์ํ ๊ณ์ฐ์ด๋ฏ๋ก ์ด๋ ๋๋ผ์์๋ ๋์ผํ๊ฒ ์ ์ฉ๋๋ฉฐ, ๊ตญ๊ฐ๋ ๋จ์์ ๋ํ ๊ฐ์ ์ ์ ํ ๋ค์ด๊ฐ์ง ์์ต๋๋ค.
์ฌ์ฉ ๋ฐฉ๋ฒ
์ฐจ์ \(\nu\)(0, 1, 2์ฒ๋ผ ์ ์๋ 1.5 ๊ฐ์ ์ค์๋ ์ด๋ค ์ค์ ๊ฐ์ด๋ ๊ฐ๋ฅํฉ๋๋ค), x์ ์ด๊ธฐ๊ฐ, ๊ฐ ํ๋ง๋ค x์ ๋ํด์ง๋ ์ฆ๊ฐ๋, ๊ทธ๋ฆฌ๊ณ ํ ์๋ฅผ ์ ๋ ฅํ์ธ์. ๊ฐ ํ \(k\)์ x๋ $$x_k = \text{์ด๊ธฐ๊ฐ} + k\cdot\text{์ฆ๊ฐ๋}$$ ์ผ๋ก ๊ณ์ฐ๋ฉ๋๋ค. ๋งจ ์์ ํฌ๊ฒ ํ์๋๋ ๋ํ ์ซ์๋ ์ฒซ ๋ฒ์งธ x์์์ \(j_{\nu}\) ๊ฐ์ด๋ฉฐ, ํ์๋ ์์ฑ๋ ๋ชจ๋ ๊ฐ์ด ๋์ด๋ฉ๋๋ค.
๊ณต์ ํ์ด
์ผ๋ฐ์ ์ธ ์ค์ ์ฐจ์์ ๊ฒฝ์ฐ $$j_{\nu}(x) = \sqrt{\frac{\pi}{2x}}\; J_{\nu+\frac{1}{2}}(x)$$ ๋ก ๊ณ์ฐํฉ๋๋ค. ์ฌ๊ธฐ์ \(J\)๋ ๋ณดํต์ ์ 1์ข ๋ฒ ์ ํจ์๋ก, ๊ฑฐ๋ญ์ ๊ณฑ ๊ธ์์ Lanczos ๊ฐ๋ง ํจ์๋ฅผ ํตํด ๊ณ์ฐ๋ฉ๋๋ค. ์ ์ ์ฐจ์์์๋ ์์น์ ์ผ๋ก ์์ ์ ์ธ ๋ซํ ํํ์ธ \(j_0(x) = \sin(x)/x\) ์ \(j_1(x) = \sin(x)/x^2 - \cos(x)/x\) ๋ฅผ ์ฌ์ฉํ ๋ค, ์ํฅ ์ ํ์ $$j_{n+1}(x) = \frac{2n+1}{x}\cdot j_n(x) - j_{n-1}(x)$$ ๋ก ์ฐจ์๋ฅผ ์ฌ๋ ค ๊ฐ๋๋ค. \(x = 0\)์์๋ ๊ทนํ๊ฐ์ ์ ์ฉํด \(j_0(0) = 1\), \(\nu > 0\)์์๋ \(j_n(0) = 0\) ์ผ๋ก ์ฒ๋ฆฌํ์ฌ 0์ผ๋ก ๋๋๋ ๋ฌธ์ ๋ฅผ ํผํฉ๋๋ค.
๊ณ์ฐ ์์
์ฐจ์ \(\nu = 0\), ์ด๊ธฐ๊ฐ 0, ์ฆ๊ฐ๋ 0.2, ํ ์ 6์ผ๋ก ์ค์ ํ๋ฉด \(x = 0, 0.2, 0.4, 0.6, 0.8, 1.0\) ์ด ๋ฉ๋๋ค. \(j_0(x) = \sin(x)/x\) ๋ฅผ ์ฌ์ฉํ๋ฉด $$j_0(0)=1,\quad j_0(0.2)=0.993347,\quad j_0(0.4)=0.973546,$$ $$j_0(0.6)=0.941071,\quad j_0(0.8)=0.896695,\quad j_0(1.0)=0.841471$$ ๋ก, ์ฐ๋ฆฌ๊ฐ ์ตํ ์๋ ๊ฐ์ ํ sinc ๊ณก์ ๋ชจ์์ด ๋ํ๋ฉ๋๋ค.
์์ฃผ ๋ฌป๋ ์ง๋ฌธ
์ฐจ์๋ฅผ ์์๋ก ๋ฃ์ด๋ ๋๋์? ๋ค. ์ ์๊ฐ ์๋ \(\nu\)์ ๋ํด์๋ \(\sqrt{\pi/2x}\cdot J\) ๊ธ์ ํํ๋ฅผ ์ฌ์ฉํฉ๋๋ค.
x๊ฐ 0์์ ์์ํ๋๋ฐ ์ฒซ ๊ฐ์ด ์ ํํ 1์ธ ์ด์ ๋? ๊ทนํ์ ์ํด \(j_0(0) = 1\) ์ด๊ธฐ ๋๋ฌธ์ ๋๋ค. \(\nu > 0\)์ธ ๊ฒฝ์ฐ์๋ ๊ทนํ๊ฐ์ด 0์ ๋๋ค.
์ํฅ ์ ํ์์ ์ธ์ ๋ ์์ ํ๊ฐ์? ํ๋ฅผ ๋ณด๊ธฐ ์ํ ์ ๋นํ ์ฐจ์์ ์ผ๋ฐ์ ์ธ x ๋ฒ์์์๋ ์์ ํฉ๋๋ค. x์ ๋นํด ์ฐจ์๊ฐ ๋งค์ฐ ํฐ ๊ฒฝ์ฐ์๋ ํํฅ ์ ํ์์ด ๋ ์์ ์ ์ด์ง๋ง, ์ด ๊ณ์ฐ๊ธฐ์์๋ ๊ทธ๋ฐ ์ํฉ์ด ๊ฑฐ์ ์์ต๋๋ค.