이 계산기의 기능
이 도구는 제2종 변형 베셀 함수, 즉 \(K_{\nu}(x)\)를 표로 정리하고 그래프로 그려 줍니다. 고정된 실수 차수 \(\nu\)와 훑어볼 \(x\) 값들을 입력하면 \(x\)와 \(K_{\nu}(x)\)를 짝지은 2열 표를 돌려주고, 함수가 어떻게 감쇠하는지 보여 주는 선 그래프도 함께 표시합니다. 순수 수학이므로 지역이나 국가에 관계없이 어디서나 그대로 적용됩니다.
배경 지식
변형 베셀 함수는 변형 베셀 방정식 $$x^{2}y'' + xy' - (x^{2} + \nu^{2})y = 0$$의 두 독립적인 해입니다. 제1종 \(I_{\nu}(x)\)는 지수적으로 증가하고, 제2종 \(K_{\nu}(x)\)는 지수적으로 감쇠합니다. 한 가지 알아둘 점은 \(K_{\nu}(x) = K_{-\nu}(x)\)라는 사실입니다. 즉 차수의 부호는 결과에 영향을 주지 않으며, 계산기는 내부적으로 \(|\nu|\)를 사용합니다.
사용 방법
차수 \(\nu\)(임의의 실수), \(x\)의 시작값, 각 행마다 \(x\)에 더해지는 증가폭, 그리고 반복 횟수(\(x\) 표본점 개수, 곧 표의 행 수)를 입력하세요. \(i\)번째 행은 $$x = \text{시작값} + i\cdot\text{증가폭}$$으로 계산됩니다. \(K_{\nu}(x)\)는 \(x > 0\)에서만 정의되며 \(x \to 0^{+}\)일 때 \(+\infty\)로 발산하므로, \(x = 0\)인 행(또는 음수 \(x\))은 무한대(Infinity)로 표시됩니다.
공식 풀이
이 계산기는 적분 표현 $$K_{\nu}(x) = \int_{0}^{\infty} e^{-x \cosh t}\,\cosh\!\left(\nu\, t\right)\,dt$$를 합성 심프슨(Simpson) 공식으로 수치 적분하여 계산합니다. 피적분 함수가 무시할 만큼 작아지는 지점(\(x\cdot\cosh t\)가 약 45를 넘는 곳)에서 적분 상한을 잘라냅니다. 이 표현은 정수 차수에서 특이점이 없으므로 \(K_{0}\), \(K_{1}\), \(\ldots\) 같은 값도 \(I_{\nu}\)를 포함한 닫힌 형식의 0/0 문제 없이 곧바로 계산됩니다.
계산 예시
\(\nu = 0\), 시작값 = 0.1, 증가폭 = 0.1, 반복 횟수 = 3으로 두면 표는 다음과 같습니다. $$x = 0.1 \rightarrow 2.427069, \quad x = 0.2 \rightarrow 1.752704, \quad x = 0.3 \rightarrow 1.372460.$$ 이 값들은 표준 \(K_{0}\) 수표 값과 일치합니다.
자주 묻는 질문
첫 번째 값이 가끔 \(\infty\)로 나오는 이유는? \(K_{\nu}(0)\)이 발산하기 때문입니다. 0.1처럼 작은 양수를 시작값으로 잡으세요.
음수 차수도 쓸 수 있나요? 네. \(K_{\nu}(x) = K_{-\nu}(x)\)이므로 \(-\nu\)의 결과는 \(\nu\)의 결과와 똑같습니다.
\(x\)가 클 때 값이 0이 되는 이유는? \(K_{\nu}(x)\)는 \(\sqrt{\pi/2x}\cdot e^{-x}\)처럼 감쇠하므로 \(x\)가 커지면 0으로 언더플로됩니다. 이는 올바른 결과입니다.
