이 계산기의 기능
이 도구는 제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)$$ 가 더 적합합니다.
