logo

버금 르장드르 다항식의 직교성 📂함수

버금 르장드르 다항식의 직교성

정리

구간 [1,1][-1,1]에서 고정된 mm에 대한 버금 르장드르 다항식직교 집합을 이룬다.

11Plm(x)Pkm(x)dx=22l+1(l+m)!(lm)!δlk \int_{-1}^{1} P_{l}^{m}(x)P_{k}^{m}(x)dx =\frac{ 2}{ 2l+1 }\frac{(l+m)!}{(l-m)!}\delta_{lk}

x=cosθx=\cos \theta일 경우에는,

0πPlm(cosθ)Pkm(cosθ)sinθdθ=22l+1(l+m)!(lm)!δlk \int_{0}^{\pi} P_{l}^{m}(\cos \theta)P_ {k}^{m}(\cos\theta)\sin \theta d\theta =\frac{ 2}{ 2l+1 }\frac{(l+m)!}{(l-m)!}\delta_{lk}

버금 르장드르 다항식 Plm(x)=(1x2)m212ll!dl+mdxl+m(x21)l P_{l}^{m}(x) = (1-x ^{2})^{\frac{m}{2}} \dfrac{1}{2^l l!} \dfrac{d^{l+m}}{dx^{l+m}}(x^2-1)^l

증명

우선 편의를 위해 간단히 Plm=Plm(x)P_{lm} = P_{l}^{m}(x)과 같이 표기하자. 버금 르장드르 미분방정식은 아래와 같다.

ddx[(1x2)Plm]+[l(l+1)m21x2]Plm=0 \begin{equation} \frac{d}{dx} \left[ (1-x^{2})P_{lm}^{\prime} \right] +\left[ l(l+1)-\frac{m^{2}}{1-x^{2}} \right]P_{lm} = 0 \end{equation}

경우 1: lkl \ne k

증명 방식은 르장드르 다항식의 직교성을 보이는 것과 같다. (1)(1)ll, kk에 대해서 쓰면,

ddx[(1x2)Plm]+[l(l+1)m21x2]Plm=0ddx[(1x2)Pkm]+[k(k+1)m21x2]Pkm=0 \frac{d}{dx} \left[(1-x^{2}) P_{lm}^{\prime}\right] + \left[l(l+1) - \frac{m^{2}}{1-x^{2}}\right] P_{lm} = 0 \\ \frac{d}{dx} \left[(1-x^{2}) P_{km}^{\prime}\right] + \left[k(k+1) - \frac{m^{2}}{1-x^{2}}\right] P_{km} = 0

위의 식(ll에 대한 식)에 PkmP_{km}을 곱하고, 아래의 식(kk에 대한 식)에 PlmP_{lm}을 곱한 뒤 서로 빼주면 다음과 같다.

Pkmddx[(1x2)Plm]Plmddx[(1x2)Pkm]+[l(l+1)k(k+1)]PlmPkm=0 \begin{equation} P_{km} \frac{d}{dx} \left[ (1-x^{2})P_{lm}^{\prime} \right]-P_{lm}\frac{d}{dx} \left[ (1-x^{2})P_{km}^{\prime} \right]+\left[l(l+1)- k(k+1) \right]P_{lm}P_{km} = 0 \end{equation}

여기서 첫째항, 둘째항을 다음과 같이 정리할 수 있다.

 Pkmddx[(1x2)Plm]Plmddx[(1x2)Pkm]=Pkm(1x2)Plm+Pkm(1x2)PlmPlm(1x2)PkmPlm(1x2)Pkm=Pkm(1x2)Plm+Pkm(1x2)PlmPlm(1x2)PkmPlm(1x2)Pkm+Pkm(1x2)PlmPkm(1x2)Plm=ddx[(1x2)PlmPkm(1x2)PlmPkm]=ddx[(1x2)(PlmPkmPlmPkm)] \begin{align*} &\quad \ P_{km} \frac{d}{dx} \left[ (1-x^{2})P_{lm}^{\prime} \right]-P_{lm}\frac{d}{dx} \left[ (1-x^{2})P_{km}^{\prime} \right] \\ &= P_{km}(1-x^{2})^{\prime}P_{lm}^{\prime}+P_{km}(1-x^{2})P_{lm}^{\prime \prime}-P_{lm}(1-x^{2})^{\prime}P_{km}^{\prime}-P_{lm}(1-x^{2})P_{km}^{\prime \prime} \\ &= {\color{blue}P_{km}(1-x^{2})^{\prime}P_{lm}^{\prime}+P_{km}(1-x^{2})P_{lm}^{\prime \prime}}-{\color{orange}P_{lm}(1-x^{2})^{\prime}P_{km}^{\prime}-P_{lm}(1-x^{2})P_{km}^{\prime \prime} } \\ &\quad + {\color{blue}P_{km}^{\prime}(1-x^{2})P_{lm}^{\prime}}-{\color{orange}P_{km}^{\prime}(1-x^{2})P_{lm}^{\prime}} \\ &= \frac{d}{dx}\left[{\color{blue}(1-x^{2})P_{lm}^{\prime}P_{km}}-{\color{orange}(1-x^{2})P_{lm}P_{km}^{\prime} } \right] \\ &= \frac{d}{dx}\left[(1-x^{2})(P_{lm}^{\prime}P_{km}-P_{lm}P_{km}^{\prime}) \right] \end{align*}

이를 (2)(2)에 대입하면,

ddx[(1x2)(PlmPkmPlmPkm)]+[l(l+1)k(k+1)]PlmPkm=0 \frac{d}{dx}\left[(1-x^{2})(P_{lm}^{\prime}P_{km}-P_{lm}P_{km}^{\prime}) \right]+\left[l(l+1)- k(k+1) \right]P_{lm}P_{km}=0

양변을 구간 [1,1][-1,1]에서 적분하면 다음을 얻는다.

[(1x2)(PlmPkmPlmPkm)]11+[l(l+1)k(k+1)]11PlmPkmdx=0 \left[(1-x^{2})(P_{lm}^{\prime}P_{km}-P_{lm}P_{km}^{\prime}) \right]_{-1}^{1}+\left[l(l+1)- k(k+1) \right]\int_{-1}^{1}P_{lm}P_{km}dx=0

첫번째 항은 00이므로 아래와 같다.

[l(l+1)k(k+1)]11PlmPkmdx=0 \left[l(l+1)- k(k+1) \right]\int_{-1}^{1}P_{lm}P_{km}dx=0

이때 lkl \ne k이므로 l(l+1)k(k+1)0l(l+1)-k(k+1)\ne 0이다. 그러므로 다음과 같다.

11Plm(x)Pkm(x)dx=0 \int_{-1}^{1}P_{l}^{m}(x)P_{k}^{m}(x)dx=0

경우 2: l=kl=k

보조 정리

다음이 성립한다.

dl+mdxl+m(x21)l=(l+m)!(lm)!(x21)mdlmdxlm(x21)l(3) \frac{d^{l+m}}{dx^{l+m}}(x^{2}-1)^{l}=\frac{(l+m)!}{(l-m)!}(x^{2}-1)^{-m}\frac{ d ^{l-m}}{ dx^{l-m} }(x^{2}-1)^{l} \tag {3}

위의 공식을 버금 르장드르 다항식에 대입하면,

Plm(x)=(1x2)m212ll!dl+mdxl+m(x21)l=(1x2)m212ll!(l+m)!(lm)!(x21)mdlmdxlm(x21)l=(1)m2ll!(l+m)!(lm)!(1x2)m2dlmdxlm(x21)l \begin{align*} P_{l}^{m}(x) &= (1-x ^{2})^{\frac{m}{2}} \dfrac{1}{2^l l!} \dfrac{d^{l+m}}{dx^{l+m}}(x^2-1)^l \\ &= (1-x ^{2})^{\frac{m}{2}}\dfrac{1}{2^l l!} \frac{(l+m)!}{(l-m)!}(x^{2}-1)^{-m}\frac{ d ^{l-m}}{ dx^{l-m} }(x^{2}-1)^{l} \\ &= \dfrac{(-1)^{m}}{2^l l!} \frac{(l+m)!}{(l-m)!}(1-x^{2})^{-\frac{m}{2}}\frac{ d ^{l-m}}{ dx^{l-m} }(x^{2}-1)^{l} \end{align*}

위 식의 양변을 제곱하면 다음을 얻는다.

[Plm(x)]2=122l(l!)2[(l+m)!(lm)!]2(1x2)mdlmdxlm(x21)ldlmdxlm(x21)l [P_{l}^{m}(x)]^{2} =\dfrac{1}{2^{2l}(l!)^{2}} \left[ \frac{(l+m)!}{(l-m)!} \right]^{2}(1-x^{2})^{-m}\frac{ d ^{l-m}}{ dx^{l-m} }(x^{2}-1)^{l}\frac{ d ^{l-m}}{ dx^{l-m} }(x^{2}-1)^{l}

위 식에 다시 (3)(3)을 대입하면,

 [Plm(x)]2=122l(l!)2[(l+m)!(lm)!]2(1x2)m[(lm)!(l+m)!(x21)mdl+mdxl+m(x21)l]dlmdxlm(x21)l=(1)m22l(l!)2(l+m)!(lm)!dl+mdxl+m(x21)ldlmdxlm(x21)l \begin{align*} &\quad \ [P_{l}^{m}(x)]^{2} \\ &= \dfrac{1}{2^{2l}(l!)^{2}} \left[ \frac{(l+m)!}{(l-m)!} \right]^{2}(1-x^{2})^{-m}\left[ \frac{(l-m)!}{(l+m)!}(x^{2}-1)^{m}\frac{ d ^{l+m}}{ dx^{l+m} }(x^{2}-1)^{l} \right]\frac{ d ^{l-m}}{ dx^{l-m} }(x^{2}-1)^{l} \\ &= \dfrac{(-1)^{m}}{2^{2l}(l!)^{2}} \frac{(l+m)!}{(l-m)!} \frac{ d ^{l+m}}{ dx^{l+m} }(x^{2}-1)^{l}\frac{ d ^{l-m}}{ dx^{l-m} }(x^{2}-1)^{l} \end{align*}

이제 양변을 구간 [1,1][-1,1]에서 적분하면 아래와 같다.

 11[Plm(x)]2dx=(1)m22l(l!)2(l+m)!(lm)!11[dl+mdxl+m(x21)ldlmdxlm(x21)l]dx \begin{align*} &\quad \ \int_{-1}^{1}[P_{l}^{m}(x)]^{2}dx \\ &= \dfrac{(-1)^{m}}{2^{2l}(l!)^{2}} \frac{(l+m)!}{(l-m)!}\int_{-1}^{1} \left[ \frac{ d ^{l+m}}{ dx^{l+m} }(x^{2}-1)^{l}\frac{ d ^{l-m}}{ dx^{l-m} }(x^{2}-1)^{l} \right]dx \tag{4} \end{align*}

우변의 적분 부분만 살펴보자. 부분적분으로 풀어내면 다음과 같다.

 11[dl+mdxl+m(x21)ldlmdxlm(x21)l]dx=11[dl+m1dxl+m1(x21)l]dlmdxlm(x21)ldx=[dl+m1dxl+m1(x21)ldlmdxlm(x21)l]1111dl+m1dxl+m1(x21)ldlm+1dxlm+1(x21)ldx \begin{align*} &\quad \ \int_{-1}^{1} \left[ \frac{ d ^{l+m}}{ dx^{l+m} }(x^{2}-1)^{l}\frac{ d ^{l-m}}{ dx^{l-m} }(x^{2}-1)^{l} \right]dx \\ &= \int_{-1}^{1} \left[ \frac{ d ^{l+m-1}}{ dx^{l+m-1} }(x^{2}-1)^{l}\right]^{\prime}\frac{ d ^{l-m}}{ dx^{l-m} }(x^{2}-1)^{l} dx \\ &= \left[ \frac{ d ^{l+m-1}}{ dx^{l+m-1} }(x^{2}-1)^{l}\frac{ d ^{l-m}}{ dx^{l-m} }(x^{2}-1)^{l}\right]_{-1}^{1}-\int_{-1}^{1}\frac{ d ^{l+m-1}}{ dx^{l+m-1} }(x^{2}-1)^{l}\frac{ d ^{l-m+1}}{ dx^{l-m+1} }(x^{2}-1)^{l}dx \end{align*}

여기서 첫째항은 00이다. (x21)l(x^{2}-1)^{l}2l2l차 다항식이고 m<l|m| \lt l이므로 l+m1l+m-1lml-m모두 ll보다 작아서 적어도 (x21)(x^{2}-1)가 미분되지 않고 남아있기 때문이다. 여기에 ±1\pm 1을 대입하면 00이다. 남은 항을 다시 부분 적분으로 풀업보면,

 11dl+m1dxl+m1(x21)ldlm+1dxlm+1(x21)ldx=[dl+m2dxl+m2(x21)ldlm+1dxlm+1(x21)l]+11dl+m2dxl+m2(x21)ldlm+2dxlm+2(x21)ldx \begin{align*} &\quad \ -\int_{-1}^{1}\frac{ d ^{l+m-1}}{ dx^{l+m-1} }(x^{2}-1)^{l}\frac{ d ^{l-m+1}}{ dx^{l-m+1} }(x^{2}-1)^{l}dx \\ &= \left[- \frac{ d ^{l+m-2}}{ dx^{l+m-2} }(x^{2}-1)^{l}\frac{ d ^{l-m+1}}{ dx^{l-m+1} }(x^{2}-1)^{l} \right]+\int_{-1}^{1} \frac{ d ^{l+m-2}}{ dx^{l+m-2} }(x^{2}-1)^{l}\frac{ d ^{l-m+2}}{ dx^{l-m+2} }(x^{2}-1)^{l}dx \end{align*}

여기서 첫번째 항은 위와 같은 이유로 00이다. 이런 식으로 부분 적분을 mm번 반복하면 아래와 같은 식을 얻는다.

11dl+mdxl+m(x21)ldlmdxlm(x21)ldx=(1)m11dldxl(x21)ldldxl(x21)ldx \int_{-1}^{1} \frac{ d ^{l+m}}{ dx^{l+m} }(x^{2}-1)^{l}\frac{ d ^{l-m}}{ dx^{l-m} }(x^{2}-1)^{l} dx=(-1)^{m}\int_{-1}^{1}\frac{ d ^{l}}{ dx^{l} }(x^{2}-1)^{l}\frac{d^{l}}{dx^{l}}(x^{2}-1)^{l}dx

따라서 (4)(4)는 다음과 같다.

11[Plm(x)]2dx=122l(l!)2(l+m)!(lm)!11[dldxl(x21)l]2dx \int_{-1}^{1}[P_{l}^{m}(x)]^{2}dx= \dfrac{1}{2^{2l}(l!)^{2}} \frac{(l+m)!}{(l-m)!}\int_{-1}^{1}\left[ \frac{ d ^{l}}{ dx^{l} }(x^{2}-1)^{l}\right]^{2}dx

로드리게스 공식

Pl(x)=12ll!dldxl(x21)l P_{l}(x)=\dfrac{1}{2^l l!} \dfrac{d^l}{dx^l}(x^2-1)^l

그러면 로드리게스 공식에 의해 다음과 같다.

11[Plm(x)]2dx=122l(l!)2(l+m)!(lm)!22l(l!)211[Pl(x)]2dx=(l+m)!(lm)!11[Pl(x)]2dx \begin{align*} \int_{-1}^{1}[P_{l}^{m}(x)]^{2}dx &= \dfrac{1}{2^{2l}(l!)^{2}} \frac{(l+m)!}{(l-m)!}2^{2l}(l!)^{2}\int_{-1}^{1}\left[ P_{l}(x)\right]^{2}dx \\ &= \frac{(l+m)!}{(l-m)!}\int_{-1}^{1}\left[ P_{l}(x)\right]^{2}dx \end{align*}

그러면 르장드르 다항식의 직교성에 의해 최종적으로 다음을 얻는다.

11[Plm(x)]2dx=22l+1(l+m)!(lm)! \int_{-1}^{1}[P_{l}^{m}(x)]^{2}dx = \frac{2}{2l+1}\frac{(l+m)!}{(l-m)!}