베셀 함수 계산기 쉽게 이해하기 — 개념과 원리
이 계산기는 0부터 100까지의 정수 또는 소수 차수 ν와 인수 x에 대해 제1종 베셀 함수 Jν(x)와 제2종 베셀 함수 Yν(x)의 값을 구합니다. 베셀 함수는 진동하는 북의 막, 원기둥 속 열 흐름, 원형 도파관의 전자기파, 신호 처리 등에 나옵니다.
계산 방법
베셀 방정식: x² w″ + x w′ + (x² − ν²) w = 0
급수: Jν(x) = (x/2)^ν Σ (−x²/4)^k ÷ (k! Γ(ν + k + 1))
정수 차수: Jₙ(x) = (1/π) ∫₀^π cos(x sin θ − nθ) dθ
Yν(x) = (1/π) ∫₀^π sin(x sin θ − νθ) dθ − (1/π) ∫₀^∞ (e^(νt) + e^(−νt) cos νπ) e^(−x sinh t) dt (x > 0)
- 작은 x에는 급수를, 그 밖에는 적분 표현을 쓰고, 높은 차수는 세 항 점화식으로 구하며 유효숫자 6자리로 표시합니다.
- x → 0이면 Yν(x) → −∞이므로 x ≤ 0에서는 —로 표시합니다. 정수 n이면 Jₙ(−x) = (−1)ⁿ Jₙ(x)입니다.
원리
베셀 방정식은 파동이나 확산 문제를 원기둥 좌표로 쓰면 나타납니다. Jν는 중심에서 유한한 해이고, Yν는 중심에서 발산하는 또 하나의 독립된 해입니다. x가 크면 두 함수 모두 진폭이 약 √(2/(πx))인 천천히 줄어드는 코사인 파처럼 행동하는데, 그래서 북이나 연못의 물결이 바깥으로 갈수록 약해집니다. 이웃한 차수는 점화식 J_(ν−1)(x) + J_(ν+1)(x) = (2ν/x) Jν(x)로 이어집니다.
예제 풀이
기본값: ν = 2, x = 5.
- 알려진 값 J₀(5) = −0.177597, J₁(5) = −0.327579입니다.
- 점화식: J₂(5) = (2 × 1 ÷ 5) J₁(5) − J₀(5) = 0.4 × (−0.327579) + 0.177597 = 0.0465651입니다.
- 같은 방식으로 Y₀(5) = −0.308518, Y₁(5) = 0.147863이므로 Y₂(5) = 0.4 × 0.147863 + 0.308518 = 0.367663입니다.
- 차수 2, 인수 5이며, 정수 차수(Jₙ, Yₙ로 표기)입니다.
| 차수 ν (0~100, 소수 가능) | Jν(x) — 제1종 베셀 함수 |
|---|---|
| 0.5 | -0.342168 |
| 1 | -0.327579 |
| 1.5 | -0.169651 |
| 2 | 0.0465651 |
| 2.5 | 0.240377 |
| 3 | 0.364831 |
| 4 | 0.391232 |
결과 해석
Jν(x)는 진동하며 부호가 바뀝니다. 그 영점은 예를 들어 원형 막의 진동수를 정합니다. J₂(5) ≈ 0.047처럼 작은 값은 x = 5가 J₂의 영점 근처라는 뜻으로, J₂의 첫 영점은 약 5.14입니다. 차수가 크고 x가 작으면 Jν는 매우 작고 Yν는 매우 큰 음수가 됩니다.
흔한 실수
- Yν(노이만 함수, Nν로도 씀)를 변형 베셀 함수 Iν, Kν와 혼동하는 경우
- x = 0에서 Yν 값이 있다고 기대하는 경우
- 차수와 인수를 바꿔 넣는 경우
- 영점 근처에서 자릿수를 너무 적게 쓰는 경우
자주 묻는 질문
베셀 함수는 어디에 쓰이나요?
북, 관, 안테나처럼 둥근 구조에서 파동, 열, 확산 문제를 풀 때 쓰입니다.
J와 Y는 무엇이 다른가요?
J는 x = 0에서 유한하고, Y는 그곳에서 −∞로 가는 또 하나의 독립된 해입니다.
차수가 소수여도 되나요?
됩니다. 반정수 차수는 사인과 코사인으로 쓸 수 있는 구면 베셀 함수와 연결됩니다.
출처: NIST Digital Library of Mathematical Functions (DLMF), §10.2 Definitions and §10.9 Integral Representations