한 변수의 정적분은 곡선 아래의 넓이였습니다. 두 변수 함수 f(x,y)를 평면의 영역에서 적분하면 곡면 아래의 부피가 되고, 밀도를 적분하면 질량이, 결합 확률밀도를 적분하면 확률이 됩니다. 계산은 놀랍도록 간단합니다. 한 변수씩 차례로 적분하면 됩니다(반복적분). 이 장의 어려움은 계산보다 “적분 범위를 정하는 것”에 있습니다.
7.1 이중적분
정의 7.1 (이중적분). 직사각형 R=[a,b]×[c,d]를 작은 직사각형으로 나누고, 각 조각의 넓이 ΔA에 그 조각 안의 함숫값을 곱해 더한 리만 합의 극한을
∬Rf(x,y)dA
로 쓰고 이중적분(double integral)이라 합니다.
미적분학 7장의 리만 합에서 폭 Δx 대신 넓이 ΔA를 쓴 것입니다. f≥0이면 곡면 z=f(x,y) 아래, R 위의 부피입니다. f=1이면 R의 넓이입니다.
7.2 반복적분과 푸비니 정리
정리 7.2 (푸비니). f가 R에서 연속이면
∬RfdA=∫ab(∫cdf(x,y)dy)dx=∫cd(∫abf(x,y)dx)dy.
안쪽 적분에서는 바깥 변수를 상수로 봅니다(편미분의 거꾸로). 부피를 얇은 조각으로 썰어, 조각의 넓이를 쌓는 것과 같습니다. 어느 방향으로 썰어도 부피는 같습니다.
예제 7.3. ∬[0,1]×[0,2](x2+y)dA=∫01[x2y+2y2]02dx=∫01(2x2+2)dx=32+2=38.
곱으로 나뉘는 경우. f(x,y)=g(x)h(y)이면 ∬RfdA=(∫abgdx)(∫cdhdy)입니다. 독립인 확률변수의 결합밀도가 이 꼴입니다.
7.3 일반 영역
영역이 직사각형이 아니면, 한 변수의 범위가 다른 변수에 따라 바뀝니다.
정리 7.4. 영역 D={(x,y)∣a≤x≤b, g1(x)≤y≤g2(x)}(위아래가 곡선)이면
∬DfdA=∫ab∫g1(x)g2(x)f(x,y)dydx.
예제 7.5. D를 y=x2과 y=x 사이(0≤x≤1)라 하면, D의 넓이는 ∫01∫x2x1dydx=∫01(x−x2)dx=61입니다(미적분학 10장과 같은 답).
적분 순서 바꾸기. 같은 영역을 “가로로” 읽으면 순서를 바꿀 수 있고, 때로 안 되던 적분이 됩니다.
예제 7.6. ∫01∫x1ey2dydx는 안쪽 적분 ∫ey2dy를 초등함수로 할 수 없습니다. 영역 {0≤x≤1, x≤y≤1}은 {0≤y≤1, 0≤x≤y}와 같으므로
∫01∫0yey2dxdy=∫01yey2dy=21(e−1).
영역을 반드시 그림으로 그려 보는 것이 요령입니다.
7.4 삼중적분
∭EfdV도 같은 방식으로 세 번 반복합니다. f=1이면 부피, f가 밀도이면 질량입니다.
예제 7.7. 사면체 x,y,z≥0, x+y+z≤1의 부피: ∫01∫01−x∫01−x−ydzdydx=∫01∫01−x(1−x−y)dydx=∫012(1−x)2dx=61.
7.5 결합 확률밀도
정의 7.8. f(x,y)≥0이고 ∬R2fdA=1이면 (X,Y)의 결합 확률밀도(joint density)이고, P((X,Y)∈D)=∬DfdA입니다.
- 주변밀도: fX(x)=∫−∞∞f(x,y)dy (다른 변수를 적분해 없앰).
- 기댓값: E[g(X,Y)]=∬g(x,y)f(x,y)dA.
- 독립: f(x,y)=fX(x)fY(y).
예제 7.9. X, Y가 각각 [0,1]에서 균등하고 독립이면 f=1([0,1]2 위). P(X+Y≤1)은 삼각형의 넓이 21, E[XY]=∫01∫01xydxdy=41입니다.
7.6 이 장의 요약
- 이중적분 = 곡면 아래 부피 = 리만 합의 극한. 계산은 반복적분.
- 일반 영역은 범위가 곡선. 순서를 바꾸면 쉬워질 수 있습니다(영역을 그립니다).
- 결합밀도의 적분이 확률, g⋅f의 적분이 기댓값입니다.
쓰이는 곳. 머신러닝 논문의 기댓값 E(x,y)∼p[ℓ(f(x),y)](기대 위험), 주변화(marginalization) p(x)=∫p(x,z)dz(잠재변수 모델, VAE)가 다중적분입니다. 물리에서는 질량, 질량중심, 관성 모멘트, 전하량, 전기장의 계산에 쓰입니다.
문제
7.1 [계산] · 확인 문제
∬[0,2]×[0,3]xydA는?
정답
풀이
곱으로 나뉘므로 ∫02xdx∫03ydy=2⋅29=9입니다.
7.2 [계산] · 확인 문제
삼각형 D={0≤x≤1, 0≤y≤x}에서 ∬D(x+y)dA는?
정답
풀이
∫01∫0x(x+y)dydx=∫01(x2+2x2)dx=23⋅31=21입니다.
7.3 [계산] · 확인 문제
적분 순서를 바꿔 ∫01∫y1xsinxdxdy를 구합니다. (소수 넷째 자리까지)
힌트 1
영역은 0≤y≤x≤1입니다. x를 바깥으로 보내면 안쪽 적분이 ∫0xdy=x가 됩니다.
정답
풀이
∫01∫0xxsinxdydx=∫01sinxdx=1−cos1≈0.4597입니다.
7.4 [계산]
X, Y가 [0,1]에서 독립인 균등분포일 때 P(Y≤X2)은?
정답
풀이
∫01∫0x21dydx=∫01x2dx=31입니다.
7.5 [증명]
결합밀도 f(x,y)=g(x)h(y)(∫g=∫h=1)이면 E[XY]=E[X]E[Y]임을 보입니다.
풀이
E[XY]=∬xyg(x)h(y)dA이고, 함수가 x의 부분과 y의 부분의 곱이므로 반복적분이 나뉩니다: ∫xg(x)dx∫yh(y)dy. 결합밀도가 곱이면 주변밀도는 fX=g, fY=h이므로 이 값은 E[X]E[Y]입니다.