logo

평균이 0인 정규분포를 따르는 확률변수 거듭제곱의 기대값 📂확률분포론

평균이 0인 정규분포를 따르는 확률변수 거듭제곱의 기대값

공식

확률변수 XX정규분포 N(0,σ2)N \left( 0 , \sigma^{2} \right) 를 따른다고 하면 그 거듭제곱 XnX^{n}기대값은 다음과 같이 재귀적인 공식으로 나타난다1. E(Xn)=(n1)σ2E(Xn2) E \left( X^{n} \right) = (n - 1) \sigma^{2} E \left( X^{n-2} \right) E(Xn)E \left( X^{n} \right)nn 이 홀수일 때 00 이고, 짝수일 때 다음과 같다2. E(X2n)=(2n1)!!σ2n E \left( X^{2n} \right) = \left( 2n - 1 \right)!! \sigma^{2n} 여기서 느낌표가 두 개 들어간 k!!=k(k2)k!! = k \cdot \left( k - 2 \right) \cdots더블 팩토리얼을 나타낸다.

설명

널리 알려진 따름정리로써 E(X4)=3σ4E \left( X^{4} \right) = 3 \sigma^{4} 이 성립하고, 이는 이토 테이블의 유도 등에서 쓰일 수 있다.

유도

결과를 얻는 두가지 방법이 있다. 가우스 적분을 통한 방식은 일반화된 공식을 통한 숏컷이라 유도가 쉽고 빠르다. 부분 적분을 통한 방식은 다소 까다롭지만 그 과정에서 재귀공식도 얻을 수 있다. 두 방법 모두 다음과 같이 E(Xn)E \left( X^{n} \right) 을 적분으로 풀어내는 것에서 시작한다. E(Xn)=12πσxnex2/2σ2dx E \left( X^{n} \right) = \int_{-\infty}^{\infty} {\frac{ 1 }{ \sqrt{2 \pi} \sigma }} x^{n} e^{-x^{2} / 2 \sigma^{2}} dx

가우스 적분을 이용한 방식

가우스 적분의 일반화: nn자연수라 하자. x2neαx2dx=(2n)!n!22nπα2n+1x2n+1eαx2dx=0 \begin{align*} \int_{-\infty}^{\infty} x^{2n} e^{-\alpha x^{2}}dx =& \dfrac{(2n)!}{n! 2^{2n}}\sqrt{\dfrac{\pi}{\alpha^{2n+1}}} \\ \int_{-\infty}^{\infty} x^{2n+1} e^{-\alpha x^{2}}dx =& 0 \end{align*}

연속된 홀수의 곱: 정수 n0n \ge 0에 대해서 다음이 성립한다. (2n1)(2n3)531=(2n)!2n(n!)=(2n1)!! (2n-1) \cdot (2n-3) \cdots 5 \cdot 3 \cdot 1 = \dfrac{(2n)!}{2^{n} (n!)} = (2n-1)!!

nn 이 홀수면 공식에 의해 더 볼 것도 없이 E(Xn)=0E \left( X^{n} \right) = 0 이다. nn 이 짝수면 α=1/2σ2\alpha = 1/2 \sigma^{2} 을 대입해서 다음을 얻는다. E(X2n)=12πσx2nex2/2σ2dx=12πσ(2n)!n!22nπ22n+1σ4n+2=(2n)!n!2n2n22nσ4n=(2n)!n!2n2n2nσ2n=(2n)!n!2nσ2n=(2n1)!!σ2n \begin{align*} E \left( X^{2n} \right) =& {\frac{ 1 }{ \sqrt{2 \pi} \sigma }} \int_{-\infty}^{\infty} x^{2n} e^{-x^{2} / 2 \sigma^{2}} dx \\ =& {\frac{ 1 }{ \sqrt{2 \pi} \sigma }} \frac{(2n)!}{n! 2^{2n}}\sqrt{\pi 2^{2n+1} \sigma^{4n+2}} \\ =& \frac{(2n)!}{n! 2^{n} 2^{n}}\sqrt{2^{2n} \sigma^{4n}} \\ =& \frac{(2n)!}{n! 2^{n} 2^{n}} 2^{n} \sigma^{2n} \\ =& \frac{(2n)!}{n! 2^{n}} \sigma^{2n} \\ =& (2n-1) !! \sigma^{2n} \end{align*}

부분 적분을 이용한 방식

In:=tnet2/2dt I_{n} := \int_{-\infty}^{\infty} t^{n} e^{-t^{2} / 2} dt InI_{n} 을 위와 같이 두고 부분적분법을 사용하자. 이 부분이 다소 까다롭다. In=tnet2/2dt=tn12t2et2/2dt=[tn1et2/2](n1)tn2et2/2dt=0(n1)In2 \begin{align*} - I_{n} =& - \int_{-\infty}^{\infty} t^{n} e^{-t^{2} / 2} dt \\ =& \int_{-\infty}^{\infty} t^{n-1} \cdot {\frac{ - 2 t }{ 2 }} e^{-t^{2} / 2} dt \\ =& \left[ t^{n-1} \cdot e^{-t^{2} / 2} \right]_{-\infty}^{\infty} - \int_{-\infty}^{\infty} (n-1) t^{n-2} e^{-t^{2} / 2} dt \\ =& 0 - (n-1) I_{n-2} \end{align*} 정리하면 In=(n1)In2I_{n} = (n-1) I_{n-2} 이고, 이를 E(Xn)E \left( X^{n} \right) 을 계산하는 과정에 적용한다. t=x/σt = x / \sigma 이라 치환하면 dx=σdtdx = \sigma dt 이므로 E(Xn)=12πσxnex2/2σ2dx=12πσ(σt)net2/2σdt=σn2πtnet2/2dt=σn2πIn=σn2π(n1)In2=(n1)σ2σn22πtn2et2/2dt=(n1)σ2σn22πtn2et2/2dt=(n1)σ2E(Xn2) \begin{align*} E \left( X^{n} \right) =& \int_{-\infty}^{\infty} {\frac{ 1 }{ \sqrt{2 \pi} \sigma }} x^{n} e^{-x^{2} / 2 \sigma^{2}} dx \\ =& \int_{-\infty}^{\infty} {\frac{ 1 }{ \sqrt{2 \pi} \sigma }} \left( \sigma t \right)^{n} e^{-t^{2} / 2 } \cdot \sigma dt \\ =& {\frac{ \sigma^{n} }{ \sqrt{2 \pi} }} \int_{-\infty}^{\infty} t^{n} e^{-t^{2} / 2 } dt \\ =& {\frac{ \sigma^{n} }{ \sqrt{2 \pi} }} I_{n} \\ =& {\frac{ \sigma^{n} }{ \sqrt{2 \pi} }} (n-1) I_{n-2} \\ =& (n-1) {\frac{ \sigma^{2} \cdot \sigma^{n-2} }{ \sqrt{2 \pi} }} \int_{-\infty}^{\infty} t^{n-2} e^{-t^{2} / 2} dt \\ =& (n-1) \sigma^{2} \int_{-\infty}^{\infty} {\frac{ \sigma^{n-2} }{ \sqrt{2 \pi} }} t^{n-2} e^{-t^{2} / 2} dt \\ =& (n-1) \sigma^{2} E \left( X^{n-2} \right) \end{align*} 을 얻는다. XX 는 평균이 00 인 정규분포로 가정했으므로 E(X1)=0E \left( X^{1} \right) = 0 이고, 모든 홀수 nn 에 대해 E(Xn)=0E \left( X^{n} \right) = 0 이다. 짝수에 대해서는 재귀공식을 전개함으로써 다음의 결과를 얻는다. E(X2n)=(2n1)σ2E(X2n2)=(2n1)σ2(2n3)σ2E(X2n4)=[(2n1)(2n3)1]σ2(n1)E(X2)=(2n1)!!σ2n \begin{align*} E \left( X^{2n} \right) =& (2n-1) \sigma^{2} E \left( X^{2n-2} \right) \\ =& (2n-1) \sigma^{2} \cdot (2n-3) \sigma^{2} E \left( X^{2n-4} \right) \\ \vdots& \\ =& [ (2n-1) (2n-3) \cdots 1 ] \sigma^{2(n-1)} E \left( X^{2} \right) \\ =& (2n-1)!! \sigma^{2n} \end{align*}


  1. grand_chat, Expected value of XnX^n for normal distribution, URL (version: 2020-11-10): https://math.stackexchange.com/q/2752327 ↩︎

  2. user65203, Proving E[X4]=3σ4E[X^4]=3σ^4, URL (version: 2018-11-26): https://math.stackexchange.com/q/1917666 ↩︎