logo

베셀 방정식의 두번째 급수해: 제2 종 베셀 함수, 노이만 함수, 베버 함수 📂함수

베셀 방정식의 두번째 급수해: 제2 종 베셀 함수, 노이만 함수, 베버 함수

정의1

베셀 방정식의 두번째 해second solution를 노이만 함수라 부르고 Nν(x)N_{\nu}(x) 혹은 Yν(x)Y_{\nu}(x)로 표기한다. 정수가 아닌 ν\nu에 대해서

Nν(x)=Yν(x)=cos(νπ)Jν(x)Jν(x)sin(νπ) N_{\nu}(x)=Y_{\nu}(x)=\frac{\cos (\nu \pi)J_{\nu}(x)-J_{-\nu}(x)}{\sin (\nu\pi)}

ν\nu가 정수일 경우 극한으로 정의한다. nZn\in \mathbb{Z}, νRZ\nu \in \mathbb{R}\setminus\mathbb{Z}에 대해서

Nn(x)=limνnNν(x) N_{n}(x)=\lim \limits_{\nu \rightarrow n}N_{\nu}(x)

이때 J±ν(x)J_{\pm \nu}(x)제1 종 베셀함수이다. 따라서 베셀 방정식의 일반해는 다음과 같다.

y(x)=AJν(x)+BNν(x) y(x)=AJ_{\nu}(x)+BN_{\nu}(x)

여기서 AA, BB는 임의의 상수

설명

x2y+xy+(x2ν2)y=0 x^{2}y^{\prime \prime} + xy^{\prime} +(x^{2}-\nu^{2})y=0

베셀 방정식의 급수 해J±ν(x)J_{\pm\nu}(x)라고 표기하고 ν\nu차 제1 종 베셀함수라 한다.

Jν(x)=n=0(1)nΓ(n+1)Γ(n+ν+1)(x2)2n+ν J_{\nu}(x)=\sum \limits_{n=0}^{\infty} \frac{(-1)^{n} }{\Gamma (n+1) \Gamma (n+\nu+1)} \left(\frac{x}{2} \right)^{2n+\nu}

Jν(x)=n=0(1)nΓ(n+1)Γ(nν+1)(x2)2nν J_{-\nu}(x)=\sum \limits_{n=0}^{\infty}\frac{(-1)^{n}}{\Gamma (n+1)\Gamma (n-\nu+1)} \left( \frac{x}{2} \right)^{2n-\nu}

보다시피 두 해는 독립이기 때문에 일반해는 아래와 같다.

y(x)=AJν(x)+BJν(x) y(x)=AJ_{\nu}(x)+BJ_{-\nu}(x)

그런데 ν\nu가 정수일 경우 두 해는 선형 독립이 아니다. 따라서 ν\nu가 정수일 때도 Jν(x)J_{\nu}(x)와 독립인 두번째 해를 구해야 한다.

잠시 sinx\sin xcosx\cos x를 생각해보자. 두 함수는 선형 독립이다. 그런데 sinx\sin xcosx\cos x의 어떤 선형 결합2sinx5cosx2\sin x -5\cos xsinx\sin x와 선형 독립이다. 이와 같은 아이디어로 Jν(x)J_{\nu}(x)Jν(x)J_{-\nu}(x)의 임의의 선형 결합을 베셀 방정식의 두번째 해라고 둔다.

Nν(x)=cos(νπ)Jν(x)Jν(x)sin(νπ) \begin{equation} N_{\nu}(x)=\frac{\cos (\nu \pi)J_{\nu}(x)-J_{-\nu}(x)}{\sin (\nu\pi)} \label{eq1} \end{equation}

Nν(x)N_{\nu}(x)ν\nu의 조건과 무관하게 Jν(x)J_{\nu}(x)와 독립이라는 것을 알 수 있다. 그런데 여기서 다시 문제가 생기는데 ν\nu가 정수이면

Nν(x)=cos(νπ)Jν(x)Jν(x)sin(νπ)=(1)νJν(x)(1)νJν(x)0=00 N_{\nu}(x)=\frac{\cos (\nu \pi)J_{\nu}(x)-J_{-\nu}(x)}{\sin (\nu\pi)}=\frac{(-1)^{\nu}J_{\nu}(x)-(-1)^{\nu}J_{\nu}(x)}{0}=\frac{0}{0}

이므로 정의되지 않는다. 따라서 ν\nu가 정수일 때는 다음과 같이 극한을 사용해서 정의한다.

Nn(x)=limνnNν(x)for nZ, νRZ N_{n}(x)=\lim \limits_{\nu \rightarrow n}N_{\nu}(x)\quad \text{for }n\in \mathbb{Z},\ \nu \in \mathbb{R}\setminus \mathbb{Z}

이때 임의의 x0x \ne 0에 대해서 위 극한이 존재한다.

정리

정수인 ν\nu에 대해서 베셀 함수 J±ν(x)J_{\pm \nu}(x)는 아래의 식을 만족한다. 다시말해 독립이 아니다.

Jν(x)=(1)νJν(x) J_{-\nu}(x)=(-1)^{\nu}J_{\nu}(x)

증명

Jν(x)=n=0(1)nΓ(n+1)Γ(nν+1)(x2)2nν J_{-\nu}(x)=\sum \limits_{n=0}^{\infty}\frac{(-1)^{n}}{\Gamma (n+1)\Gamma (n-\nu+1)} \left( \frac{x}{2} \right)^{2n-\nu}

에서 n=k+νn=k+\nu로 치환하자. 그러면

Jν(x)=k=ν(1)k+νΓ(k+ν+1)Γ(k+1)(x2)2k+ν J_{-\nu}(x)=\sum \limits_{k=-\nu}^{\infty}\frac{(-1)^{k+\nu}}{\Gamma (k+\nu+1)\Gamma (k+1)} \left( \frac{x}{2} \right)^{2k+\nu}

감마함수는 00과 음의 정수에서 발산하므로 k=ν,ν+1,,1k=-\nu,-\nu+1,\cdots,-1일 때 분모의 Γ(k+1)\Gamma (k+1)은 발산하여 Jν(x)=0J_{-\nu}(x)=0이다. 따라서

Jν(x)=k=ν(1)k+νΓ(k+ν+1)Γ(k+1)(x2)2k+ν=k=0(1)k+νΓ(k+ν+1)Γ(k+1)(x2)2k+ν=(1)νk=0(1)kΓ(k+ν+1)Γ(k+1)(x2)2k+ν=(1)νJν(x) \begin{align*} J_{-\nu}(x)&=\sum \limits_{k=-\nu}^{\infty}\frac{(-1)^{k+\nu}}{\Gamma (k+\nu+1)\Gamma (k+1)} \left( \frac{x}{2} \right)^{2k+\nu} \\ &=\sum \limits_{k=0}^{\infty}\frac{(-1)^{k+\nu}}{\Gamma (k+\nu+1)\Gamma (k+1)} \left( \frac{x}{2} \right)^{2k+\nu} \\ &=(-1)^{\nu}\sum \limits_{k=0}^{\infty}\frac{(-1)^{k}}{\Gamma (k+\nu+1)\Gamma (k+1)} \left( \frac{x}{2} \right)^{2k+\nu} \\ &=(-1)^{\nu}J_{\nu}(x) \end{align*}


  1. Mary L. Boas, 수리물리학(Mathematical Methods in the Physical Sciences, 최준곤 역) (3rd Edition, 2008), p604-606 ↩︎