베셀 함수 영점 계산기

계산 입력

차수 v

영점 번호 s (s번째)

결과

J_v(x)의 s번째 양의 영점

2.4048255577

j_v,s (제1종 베셀 함수)

Y_v(x)의 s번째 양의 영점	0.8935769664
표기법	y_v,s (제2종 베셀 함수 / 노이만)

이 계산기의 기능

이 도구는 실수 차수 v에 대해 제1종 베셀 함수 \(J_{v}(x)\)와 제2종 베셀 함수 \(Y_{v}(x)\)(노이만 함수 또는 베버 함수라고도 함)의 s번째 양의 영점을 찾아 줍니다. \(j_{v,s}\), \(y_{v,s}\)로 표기되는 이 영점들은 물리학과 공학 곳곳에서 등장합니다. 원형 막의 진동(북 가죽), 원통 내부의 열전도, 전자기 도파관, 그리고 푸리에–베셀 급수 등이 대표적입니다. 순수 수학의 특수함수 도구로서 지역이나 단위에 전혀 의존하지 않는 범용 계산기입니다.

영점은 \(J_v(x)\)와 \(Y_v(x)\)가 가로축과 만나는 x 값입니다.

사용 방법

차수 v(보통 0에서 200 사이의 실수)와 영점 번호 s(1, 2, 3, …의 양의 정수)를 입력하세요. 계산기는 \(J_{v}(x)\)의 s번째 양의 근인 \(j_{v,s}\)와 \(Y_{v}(x)\)의 s번째 양의 근인 \(y_{v,s}\)를 반환합니다. 예를 들어 v = 0, s = 1은 북 가죽의 기본 진동 모드에 해당합니다.

공식과 계산 방법

두 함수 모두 베셀 방정식 \(x^2 y'' + x y' + (x^2 - v^2)y = 0\)의 해입니다. \(J_{v}\)는 수렴하는 거듭제곱 급수로 계산하고, \(Y_{v}\)는 다음 공식을 사용하며,

$$Y_{v} = \frac{J_{v}\cos(v\pi) - J_{-v}}{\sin(v\pi)}$$

정수 차수인 경우는 수치적 극한으로 처리합니다. s번째 영점을 찾기 위해 먼저 맥마흔(McMahon) 점근 근사로 초기값을 잡고, 근을 구간으로 가둔 다음, 수렴할 때까지 이분법으로 정밀하게 다듬습니다. \(Y_{v}(x)\)는 x = 0에서 −∞로 발산하므로, 로그 특이점을 피하기 위해 탐색은 0보다 약간 큰 양수에서 시작합니다.

연속된 양의 영점 \(x_1\), \(x_2\), \(x_3\)은 곡선을 따라 s로 인덱싱됩니다.

계산 예시

v = 0, s = 1일 때: \(J_{0}(x)\)의 첫 번째 양의 영점은 \(2.4048255577\)이고, \(Y_{0}(x)\)의 첫 번째 양의 영점은 \(0.8935769663\)입니다. v = 1, s = 1일 때: \(j_{1,1} = 3.8317059702\), \(y_{1,1} = 2.1971413260\)입니다.

응용에서 영점들의 의미

베셀 함수 영점은 추상적 호기심이 아니다 — 이들은 파동, 열, 또는 퍼텐셜 문제가 원형 또는 원통형 영역에서 제기될 때마다 나타나는 이산 고유값이다. 경계 조건은 해의 방사형 부분이 경계에서 소멸하도록(또는 그 도함수가 소멸하도록) 강제하며, 이 조건은 영점 \(j_{v,s}\)에서만 만족된다.

진동하는 원형 막(드럼헤드)

가장자리가 고정된 반지름 \(a\)의 이상적인 원형 드럼의 경우, 변위는 \((v,s)\)로 표시된 모드로 분리되며, \(v\)는 각도 마디 지름의 개수를 세고 \(s\)는 방사형 마디 원의 개수를 센다. 허용되는 고유 진동수는 \(f_{v,s}=\frac{c}{2\pi a}\,j_{v,s}\)이며, 여기서 \(c\)는 파동 속도이다. 기본 음은 \(j_{0,1}=2.4048256\)을 사용하고; 더 높은 \(s\)와 더 높은 \(v\) 모두 음정을 높이며, \(j_{v,s}\)들이 서로의 정수 배수가 아니기 때문에 드럼의 배음은 비조화적이다.

원통형 열 전도

고정된 온도 표면을 가진 긴 원통에서 열 방정식을 풀 때, 방사형 고유 함수는 \(J_0(\lambda_s r/a)\)이며 \(\lambda_s=j_{0,s}\)이다. 각 모드는 시간에서 \(\exp\!\left(-\alpha (j_{0,s}/a)^2 t\right)\)로 감소하므로, 가장 작은 영점 \(j_{0,1}\)은 가장 천천히 감소하는, 가장 오래 지속되는 온도 분포를 지배한다. 더 큰 \(s\)는 더 큰 고유값을 제공하고 따라서 더 빠른 감소를 제공한다.

도파관 차단 진동수

반지름 \(a\)의 속이 빈 원형 금속 도파관에서, 횡자기(TM) 모드는 \(j_{v,s}\)로 설정된 진동수에서 차단되고 횡전기(TE) 모드는 도함수 \(J_v'\)의 영점으로 차단된다. TM 모드의 경우 차단 진동수는 \(f_{c}=\frac{c\,j_{v,s}}{2\pi a}\)이며; 이 진동수 이상에서만 모드가 전파된다. 가장 낮은 TM 모드(TM\(_{01}\))는 다시 \(j_{0,1}\)을 사용한다.

푸리에–베셀 급수

디스크 위의 모든 합리적 함수는 \(f(r)=\sum_{s=1}^{\infty} c_s\,J_v\!\left(j_{v,s}\,r/a\right)\)로 전개될 수 있다. 척도화된 영점 \(j_{v,s}/a\)은 정확히 일반 푸리에 사인 급수의 파수 \(n\pi/L\)처럼 작용하고, \(J_v(j_{v,s}r/a)\)의 직교성은 디스크 위에서(가중치 \(r\) 포함) 계수가 적분으로 계산되도록 한다, \(c_s=\frac{2}{a^2 J_{v+1}^2(j_{v,s})}\int_0^a f(r)\,J_v\!\left(j_{v,s}r/a\right) r\,dr\).

영점들이 어떻게 이동하는가

고정된 차수에 대해, 지수 \(s\)를 증가시키면 \(j_{v,s}\)가 \(\pi\)에 가까워지는 단계로 증가한다(진동은 원점에서 멀리 떨어진 곳에서 거의 주기적이 된다). 고정된 지수에 대해, 차수 \(v\)를 증가시키면 첫 번째 영점을 바깥쪽으로 대략 \(j_{v,1}\approx v + 1.8557\,v^{1/3}+\cdots\) 형태로 밀어내며, 이는 큰 \(v\)에 대해 베셀 방정식의 원심력 \(v^2/x^2\) 항이 진동의 시작을 얼마나 지연시키는지 반영한다. 실질적인 교훈: 더 높은 각도 복잡성(\(v\))과 더 많은 방사형 마디(\(s\)) 모두 더 큰 고유값 및 따라서 더 높은 진동수 또는 더 빠른 감소에 대응한다.

\(c=100\ \text{m/s}\)이고 \(a=0.30\ \text{m}\)인 구체적인 드럼 예제의 경우, 기본 진동수는 \(f_{0,1}=\frac{100}{2\pi(0.30)}\,(2.4048256)\approx 127.6\ \text{Hz}\)이다. 첫 번째 배음은 \(j_{1,1}=3.8317060\)을 사용하여 \(\approx 203\ \text{Hz}\)를 제공하므로, 배음 대 기본음 비율은 \(j_{1,1}/j_{0,1}\approx 1.593\)이다 — 청각적으로 1옥타브나 5도음정이 아니므로, 드럼이 현악기와 비교하여 음정이 없는 소리가 나는 이유이다.

자주 묻는 질문

\(Y_{v}\)의 첫 번째 영점은 왜 1보다 작은가요? \(Y_{0}\)는 원점에서 −∞로 발산했다가 첫 번째 최댓값에 도달하기 전인 \(x \approx 0.894\) 부근에서 0을 가로지릅니다. 그래서 첫 번째 영점이 \(J_{0}\)의 영점보다 훨씬 작습니다.

v가 정수가 아니어도 되나요? 됩니다. \(Y_{v}\)의 정의 공식은 정수가 아닌 모든 v에 대해 그대로 유효하며, 정수 차수는 매끄러운 극한으로 처리합니다.

결과는 얼마나 정확한가요? 계산은 배정밀도(double-precision) 연산을 사용하므로, 적당한 크기의 v와 s에 대해 약 10자리의 유효 숫자를 보장합니다.