Processing math: 0%

반응형

22. 베셀방정식



다음의 미분방정식을 베셀방정식(Bessel's equation)이라고 한다.x2y베셀방정식을 \displaystyle y''+\frac{1}{x}y'+\frac{x^{2}-\nu^{2}}{x^{2}}y=0으로 나타낼 수 있고, p(x)=1, q(x)=x^{2}-\nu^{2}는 해석적이므로 프로베니우스 해법으로 이 방정식을 풀 수 있다.

\displaystyle y=\sum_{n=0}^{\infty}{a_{n}x^{n+r}}을 미분하여 베셀방정식에 대입하면\sum_{n=0}^{\infty}{(n+r)(n+r-1)a_{n}x^{n+r}}+\sum_{n=0}^{\infty}{(n+r)a_{n}x^{n+r}}+\sum_{n=0}^{\infty}{a_{n}x^{n+r+2}}-\nu^{2}\sum_{n=0}^{\infty}{a_{n}x^{n+r}}=0이므로\sum_{n=0}^{\infty}{(n+r+\nu)(n+r-\nu)a_{n}x^{n+r}}+\sum_{n=0}^{\infty}{x^{n+r+2}}=0이고 지수를 같게 만들면\sum_{m=0}^{\infty}{(m+r+\nu)(m+r-\nu)a_{m}x^{m+r}}+\sum_{m=2}^{\infty}{a_{m-2}x^{m+r}}=0이므로(r+\nu)(r-\nu)a_{0}x^{r}+(r+1+\nu)(r+1-\nu)a_{1}x^{r+1}+\sum_{m=2}^{\infty}{\{(m+r+\nu)(m+r-\nu)a_{m}+a_{m-2}\}x^{m+r}}=0이다. x^{r}의 계수 (r+\nu)(r-\nu)a_{0}는 0이어야 하므로 r=\nu, r=-\nu이다. r=\nu일 때(2\nu+1)a_{1}=0,\,(m+2\nu)ma_{m}+a_{m-2}=0\,(m=2,\,3,\,...)이므로a_{1}=0,\,a_{m}=-\frac{a_{m-2}}{m(m+2\nu)}\,(m=2,\,3,\,...)이고 a_{1}=a_{3}=a_{5}=\cdots=0\begin{align*}a_{2}&=-\frac{a_{0}}{2(2+2\nu)}=-\frac{a_{0}}{2^{2}(\nu+1)},\,a_{4}=-\frac{a_{2}}{4(4+2\nu)}=\frac{a_{0}}{2^{4}2!(\nu+1)(\nu+2)},\,...\\a_{2m}&=\frac{(-1)^{m}a_{0}}{2^{2m}m!(\nu+1)(\nu+2)\cdots(\nu+m)}\,(m=1,\,2,\,...)\end{align*}이므로 y_{1}은 다음과 같다.y_{1}=a_{0}x^{\nu}\sum_{n=0}^{\infty}{\frac{(-1)^{n}x^{2n}}{2^{2n}n!(\nu+1)(\nu+2)\cdots(\nu+n)}}감마함수 \displaystyle\Gamma(\nu)=\int_{0}^{\infty}{e^{-t}t^{\nu-1}dt}에 대해 다음의 성질들이 성립한다.\Gamma(1)=1,\,\Gamma(n+1)=n!\,(n\in\mathbb{N}),\,\Gamma\left(\frac{1}{2}\right)=\sqrt{\pi}\displaystyle a_{0}=\frac{1}{2^{\nu}\Gamma(\nu+1)}라고 하면a_{2m}=\frac{(-1)^{m}}{2^{2m}m!(\nu+1)(\nu+2)\cdots(\nu+m)2^{\nu}\Gamma(\nu+1)}=\frac{(-1)^{m}}{2^{2m+\nu}m!\Gamma(\nu+m+1)}이고 다음과 같이 y_{1}을 나타낼 수 있는데 이 함수를 제 1종 베셀함수(Bessel function of the 1st kind)라고 한다.J_{\nu}(x)=x^{\nu}\sum_{n=0}^{\infty}{\frac{(-1)^{n}x^{2n}}{2^{2n+\nu}n!\Gamma(\nu+n+1)}}\nu가 정수가 아니면 J_{\nu}, J_{-\nu}는 독립인 반면 \nu가 정수이면 J_{\nu}J_{-\nu}는 종속이다. 따라서 \nu가 정수가 아니면 베셀방정식의 일반해는 다음과 같다.y=c_{1}J_{\nu}(x)+c_{2}J_{-\nu}(x)\nu가 0 또는 자연수이면, J_{-m}(x)=(-1)^{m}J_{m}(x)이다.

m=0,\,1,\,2,\,...라고 하면 \displaystyle J_{-m}(x)=\sum_{m=0}^{\infty}{\frac{(-1)^{n}x^{2n-m}}{2^{2n-m}n!\Gamma(-m+n+1)}}이고 0 또는 음의 정수값에서 감마함수값은 무한대(\infty)이므로J_{-m}=\sum_{n=m}^{\infty}{\frac{(-1)^{n}x^{2n-m}}{2^{2n-m}n!(n-m)!}}s=n-m이라 하면 다음의 식으로부터 성립한다.\begin{align*}J_{-m}(x)&=\sum_{s=0}^{\infty}{\frac{(-1)^{s+m}x^{2s+m}}{2^{2s+m}(s+m)!s!}}=(-1)^{m}x^{m}\sum_{s=0}^{\infty}{\frac{(-1)^{s}x^{2s}}{2^{2s+m}s!\Gamma(m+s+1)}}\\&=(-1)^{m}J_{m}(x)\end{align*}\nu가 정수이면 계수축소법으로 J_{\nu}와 일차독립인 해를 구해야 한다.


\nu가 정수이면 베셀방정식의 일반해는 다음과 같다.y=c_{1}J_{\nu}(x)+c_{2}Y_{\nu}(x)\,\left(Y_{\nu}(x)=\frac{1}{\sin\nu x}\left\{J_{\nu}(x)\cos\nu\pi-J_{-\nu}(x)\right\}\right)여기서 Y_{\nu}(x)를 제 2종 베셀함수(Bessel function of 2nd kind)라고 한다. Y_{\nu}(x)J_{\nu}에 계수축소법을 적용하여 구할 수 있으나 과정이 너무 복잡하기 때문에 식 유도는 생략하겠다. 


베셀함수 

1. (x^{\nu}J_{\nu}(x))'=x^{\nu}J_{\nu-1}(x) 

2. (x^{-\nu}J_{\nu}(x))'=-x^{-\nu}J_{\nu+1}(x) 

증명: 

1. \displaystyle x^{\nu}J_{\nu}(x)=\sum_{n=0}^{\infty}{\frac{(-1)^{n}x^{2n+2\nu}}{2^{2n+\nu}n!\Gamma(\nu+n+1)}}x에 대해 미분하면 다음의 결과를 얻는다.(x^{\nu}J_{\nu}(x))'=\sum_{n=0}^{\infty}{\frac{(-1)^{n}2(n+\nu)x^{2n+2\nu-1}}{2^{2n+\nu}n!\Gamma(\nu+n+1)}}=x^{\nu}x^{\nu-1}\sum_{n=0}^{\infty}{\frac{(-1)^{n}x^{2n}}{2^{2n+\nu-1}n!\Gamma(\nu+n)}}=x^{\nu}J_{\nu-1}(x) 

2. \displaystyle(x^{-\nu}J_{\nu}(x))'=\sum_{n=1}^{\infty}{\frac{(-1)^{n}2nx^{2n-1}}{2^{2n+\nu}n!\Gamma(\nu+n+1)}}이고 n=s+1이라고 하면\sum_{n=1}^{\infty}{\frac{(-1)^{n}2nx^{2n+1}}{2^{2n+\nu}n!\Gamma(\nu+n+1)}}=\sum_{s=0}^{\infty}{\frac{(-1)^{s+1}2(s+1)x^{2s+1}}{2^{2s+\nu+2}(s+1)!\Gamma(\nu+s+2)}}=-x^{-\nu}\sum_{s=0}^{\infty}{\frac{(-1)^{s}x^{2s+(\nu+1)}}{2^{2s+(\nu+1)}s!\Gamma((\nu+1)+s+1)}}=-x^{-\nu}J_{\nu+1}(x)이므로 (x^{-\nu}J_{\nu}(x))'=-x^{-\nu}J_{\nu+1}(x)이다.      


베셀함수 J_{\nu}(x)에 대하여 다음의 두 식 들이 성립한다.  

1. J_{\nu-1}(x)+J_{\nu+1}(x)=\frac{2\nu}{x}J_{\nu}(x) 

2. J_{\nu-1}(x)-J_{\nu+1}(x)=2J_{\nu}'(x) 

증명: 다음의 두 식x^{\nu}J_{\nu-1}(x)=\nu x^{\nu-1}J_{\nu}(x)+x^{\nu}J_{\nu}'(x),\,-x^{-\nu}J_{\nu+1}(x)=-\nu x^{-\nu-1}J_{\nu}(x)+x^{-\nu}J_{\nu}'(x)을 다음과 같이 나타낼 수 있고,J_{\nu-1}(x)=\nu x^{-1}J_{\nu}(x)+J_{\nu}'(x),\,-J_{\nu+1}(x)=-\nu x^{-1}J_{\nu}(x)+J_{\nu}'(x)이 두 식들을 서로 더하고 뺌으로써 위의 결과를 얻는다. 


베셀함수 J_{\nu}(x)에 대해 다음의 두 식들이 성립한다.J_{\frac{1}{2}}(x)=\sqrt{\frac{2}{\pi x}}\sin x,\,J_{-\frac{1}{2}}(x)=\sqrt{\frac{2}{\pi x}}\cos x그 이유는J_{\frac{1}{2}}(x)=\sqrt{x}\sum_{n=0}^{\infty}{\frac{(-1)^{n}x^{2n}}{2^{2n+\frac{1}{2}}n!\Gamma\left(n+\frac{3}{2}\right)}}=\sqrt{\frac{2}{x}}\sum_{n=0}^{\infty}{\frac{(-1)^{n}x^{2n+1}}{2^{2n+1}n!\Gamma\left(n+\frac{3}{2}\right)}}이고2^{n}n!=2^{n}n(n-1)\cdots2\cdot1=2n(2n-2)\cdots4\cdot2,\,\Gamma\left(\frac{1}{2}\right)=\sqrt{\pi}이므로\begin{align*}2^{n+1}\Gamma\left(n+\frac{3}{2}\right)&=2^{n+1}\left(n+\frac{1}{2}\right)\left(n-\frac{1}{2}\right)\cdots\frac{3}{2}\cdot\frac{1}{2}\Gamma\left(\frac{1}{2}\right)\\&=(2n+1)(2n-1)\cdots3\cdot1\cdot\sqrt{\pi}\end{align*}이고 J_{\frac{1}{2}}(x)는 다음과 같다.J_{\frac{1}{2}}(x)=\sqrt{\frac{2}{\pi x}}\sum_{n=0}^{\infty}{\frac{(-1)^{n}x^{2n+1}}{(2n+1)!}}=\sqrt{\frac{2}{\pi x}}\sin x또한 \displaystyle\sqrt{x}J_{\frac{1}{2}}(x)=\sqrt{\frac{2}{\pi}}\sin x이므로 \displaystyle\left(\sqrt{x}J_{\frac{1}{2}}(x)\right)'=\sqrt{\frac{2}{\pi}}\cos x이고 \displaystyle\left(\sqrt{x}J_{\frac{1}{2}}(x)\right)'=\sqrt{x}J_{-\frac{1}{2}}(x)이므로 \displaystyle J_{-\frac{1}{2}}(x)=\sqrt{\frac{2}{\pi x}}\cos x이다.


참고자료:

미분방정식 입문, 김병수, 양영균, 북스힐   

반응형
Posted by skywalker222