본문 바로가기
@ 필수과목/손으로 푸는 통계

[손으로 푸는 통계 ver1.0] 62. 표본분산의 분포 유도 (27) 감마함수 적분형 유도

by bigpicture 2021. 7. 10.
반응형

 

 

지난시간에 배운 내용을 잠깐 리뷰해봅시다. 

 

오일러는 12!=π212!=π2 라는 결과에 영감을 받아 팩토리얼이 적분과 관련이 있을 것이라고 생각하게 되고, 아래 적분을 떠올립니다. 

 

10xe(1x)ndx10xe(1x)ndx

 

이 적분식을 지난시간에 아래와 같이 변형했습니다. 

 

10xe(1x)ndx=n!(e+1)(e+2)(e+n)(e+n+1)10xe(1x)ndx=n!(e+1)(e+2)(e+n)(e+n+1)

 

팩토리얼과 적분을 연결한 식이 유도되었습니다. 오일러는 위 수식을 변형해서 감마함수 적분형을 유도합니다. 유도해보겠습니다. 

 

먼저 eefgfg 로 치환합시다. 

 

10xfg(1x)ndx=n!(fg+1)(fg+2)(fg+n)(fg+n+1)10xfg(1x)ndx=n!(fg+1)(fg+2)(fg+n)(fg+n+1)

 

아래와 같이 우변 분모의 각 항을 통분해줍니다. 

 

10xfg(1x)ndx=n!(f+gg)(f+2gg)(f+ngg)(f+(n+1)gg)10xfg(1x)ndx=n!(f+gg)(f+2gg)(f+ngg)(f+(n+1)gg)

 

우변을 아래와 같이 변형합시다. 

 

10xfg(1x)ndx=n!(f+g)(f+2g)(f+ng)(f+(n+1)g)gn+110xfg(1x)ndx=n!(f+g)(f+2g)(f+ng)(f+(n+1)g)gn+1

 

아래와 같이 우변을 둘로 분리해줍니다. 

 

10xfg(1x)ndx=gn+1f+(n+1)gn!(f+g)(f+2g)(f+ng)10xfg(1x)ndx=gn+1f+(n+1)gn!(f+g)(f+2g)(f+ng)

 

아래와 같이 우변의 첫항을 좌변으로 옮겨줍니다. 

 

f+(n+1)ggn+110xfg(1x)ndx=n!(f+g)(f+2g)(f+ng)f+(n+1)ggn+110xfg(1x)ndx=n!(f+g)(f+2g)(f+ng)

 

좌우 변을 바꿔줍니다. 

 

n!(f+g)(f+2g)(f+ng)=f+(n+1)ggn+110xfg(1x)ndxn!(f+g)(f+2g)(f+ng)=f+(n+1)ggn+110xfg(1x)ndx

 

치환을 할겁니다. 아래와 같이 치환합니다. 

 

x=sgf+gx=sgf+g 

 

dx=gf+gsff+gdsdx=gf+gsff+gds

 

위 식에 대입합시다. 

 

n!(f+g)(f+2g)(f+ng)=f+(n+1)ggn+110sff+g(1sgf+g)ngf+gsff+gdsn!(f+g)(f+2g)(f+ng)=f+(n+1)ggn+110sff+g(1sgf+g)ngf+gsff+gds

 

아래와 같이 소거합니다. 

 

n!(f+g)(f+2g)(f+ng)=f+(n+1)ggn+110gf+g(1sgf+g)ndsn!(f+g)(f+2g)(f+ng)=f+(n+1)ggn+110gf+g(1sgf+g)nds

 

아래와 같이 변형합니다. 

 

n!(f+g)(f+2g)(f+ng)=f+(n+1)ggn+110(gf+g)n+1(gf+g)n(1sgf+g)ndsn!(f+g)(f+2g)(f+ng)=f+(n+1)ggn+110(gf+g)n+1(gf+g)n(1sgf+g)nds

 

아래와 같이 변형합니다. 적분식 안의 첫 항의 분자는 적분식 바깥의 식과 곱해주고, 분모는 적분식 안의 두번째 항과 합쳐준 것입니다. 

 

n!(f+g)(f+2g)(f+ng)=f+(n+1)g(f+g)n+110(1sgf+gg/f+g)ndsn!(f+g)(f+2g)(f+ng)=f+(n+1)g(f+g)n+110(1sgf+gg/f+g)nds

 

이제 f를 1로 보내고, g를 0으로 보낼 것입니다. 

 

limf1,g0n!(f+g)(f+2g)(f+ng)=limf1,g0f+(n+1)g(f+g)n+110(1sgf+gg/f+g)ndslimf1,g0n!(f+g)(f+2g)(f+ng)=limf1,g0f+(n+1)g(f+g)n+110(1sgf+gg/f+g)nds

 

적분안의 항만 따로 계산해봅시다. 아래와 같습니다. 

 

limf1,g01sgf+gg/f+glimf1,g01sgf+gg/f+g

 

계산을 편하게 하기 위해 gf+ggf+g 를 z로 치환합시다. 

 

limz01szzlimz01szz

 

로피탈 정리를 이용하여 극한값을 구하겠습니다. 분자와 분모를 z로 미분합니다. (ax)=lnaax 를 이용합니다. 

 

limz01szz=limz0lnssz1=lns

 

우리가 계산하던 식을 다시 가져옵시다. 

 

limf1,g0n!(f+g)(f+2g)(f+ng)=limf1,g0f+(n+1)g(f+g)n+110(1sgf+gg/f+g)nds

 

극한값은 아래와 같습니다. 

 

n!=10(lns)ndt

 

치환을 하겠습니다. lns를 t로 치환합시다. 

 

lns=ts=etds=etdt

 

적분범위를 구하기 위해 그래프를 그리면 아래와 같습니다. 

 

 

s가 0에서 1로 변할 때, t는 무한대에서 0으로 변합니다. 

 

n!=0tn(et)dt

 

아래와 같이 변형합시다. 적분구간을 뒤집으면 부호가 바뀝니다. 

 

n!=0tnetdt

 

감마함수는 팩토리얼함수에서 정의역을 확장한 것이므로 n대신 n-1을 넣어서 팩토리얼 함수로 만들어줍니다. 

 

(n1)!=0tn1etdt

 

n을 실수 영역으로 확장하면 감마함수가 됩니다. 

 

Γ(z)=0tz1etdt

반응형

댓글

bigpicture님의
글이 좋았다면 응원을 보내주세요!