7.33 이중 적분의 정의와 계산 순서

이중 적분(double integral)은 2차원 영역 위에서 정의된 함수의 총량을 계산하는 수학적 도구이다. 로봇공학에서는 평면 링크의 질량 분포, 관성 모멘트, 작업 공간 내 확률 밀도의 적분 등 다양한 문제에 이중 적분이 적용된다. 이 절에서는 이중 적분의 엄밀한 정의, 반복 적분(iterated integral)으로의 환원, 그리고 적분 순서 교환의 조건과 방법을 다룬다.

1. 이중 적분의 정의

1.1 리만 합에 의한 정의

$xy$ -평면의 유계 닫힌 영역 $D$ 에서 정의된 함수 $f(x, y)$ 에 대해, $D$ 를 $n$ 개의 소영역 $\Delta A_1, \Delta A_2, \ldots, \Delta A_n$ 으로 분할한다. 각 소영역 $\Delta A_k$ 내에서 임의의 점 $(x_k^*, y_k^*)$ 를 선택하면, 리만 합은 다음과 같다.

$S_n = \sum_{k=1}^{n} f(x_k^*, y_k^*) \Delta A_k$

분할을 점점 더 세밀하게 하여, 모든 소영역의 지름이 $0$ 으로 수렴할 때 리만 합의 극한이 존재하면, 이를 $f$ 의 $D$ 위에서의 이중 적분이라 정의한다.

$\iint_D f(x, y) \, dA = \lim_{\| P \| \to 0} \sum_{k=1}^{n} f(x_k^*, y_k^*) \Delta A_k$

여기서 $\| P \|$ 는 분할의 노름(norm)으로, 가장 큰 소영역의 지름을 의미한다.

1.2 기하학적 해석

$f(x, y) \geq 0$ 인 경우, 이중 적분 $\iint_D f(x, y) \, dA$ 는 영역 $D$ 위에 곡면 $z = f(x, y)$ 아래에 놓인 입체의 부피를 나타낸다. $f(x, y) = 1$ 이면 이중 적분은 영역 $D$ 의 면적 $A(D)$ 를 반환한다.

$A(D) = \iint_D 1 \, dA$

2. 반복 적분과 푸비니 정리

2.1 푸비니 정리 (Fubini’s Theorem)

이중 적분을 실제로 계산하기 위해서는 단일 변수 적분의 반복으로 변환해야 한다. 푸비니 정리는 이 변환의 조건과 방법을 제시한다.

정리 (푸비니). 함수 $f(x, y)$ 가 직사각형 영역 $R = [a, b] \times [c, d]$ 위에서 연속이면, 다음이 성립한다.

$\iint_R f(x, y) \, dA = \int_a^b \left( \int_c^d f(x, y) \, dy \right) dx = \int_c^d \left( \int_a^b f(x, y) \, dx \right) dy$

즉, 적분 순서를 교환하더라도 결과가 동일하다.

2.2 일반 영역에서의 반복 적분

직사각형이 아닌 일반 영역에서는 영역의 형태에 따라 적분 한계가 함수로 표현된다.

제1형 영역 (Type I). $x$ 에 대해 먼저 한계를 고정하고, $y$ 의 한계가 $x$ 의 함수인 경우이다.

$D = \{(x, y) \mid a \leq x \leq b, \; g_1(x) \leq y \leq g_2(x)\}$

이때 이중 적분은 다음과 같이 계산된다.

$\iint_D f(x, y) \, dA = \int_a^b \int_{g_1(x)}^{g_2(x)} f(x, y) \, dy \, dx$

제2형 영역 (Type II). $y$ 에 대해 먼저 한계를 고정하고, $x$ 의 한계가 $y$ 의 함수인 경우이다.

$D = \{(x, y) \mid c \leq y \leq d, \; h_1(y) \leq x \leq h_2(y)\}$

이때 이중 적분은 다음과 같이 계산된다.

$\iint_D f(x, y) \, dA = \int_c^d \int_{h_1(y)}^{h_2(y)} f(x, y) \, dx \, dy$

3. 적분 순서의 결정

3.1 영역 분석에 의한 순서 선택

적분 순서는 영역 $D$ 의 경계가 어떤 변수의 함수로 더 자연스럽게 표현되는가에 따라 결정된다. 예를 들어 삼각형 영역

$D = \{(x, y) \mid 0 \leq x \leq 1, \; 0 \leq y \leq x\}$

에 대해 $\iint_D f(x, y) \, dA$ 를 계산하면, 제1형으로는 다음과 같다.

$\int_0^1 \int_0^x f(x, y) \, dy \, dx$

동일한 영역을 제2형으로 표현하면 $D = \{(x, y) \mid 0 \leq y \leq 1, \; y \leq x \leq 1\}$ 이므로 다음과 같다.

$\int_0^1 \int_y^1 f(x, y) \, dx \, dy$

3.2 적분 순서 교환

때로는 한 순서에서는 적분이 불가능하지만 순서를 교환하면 계산이 가능해지는 경우가 있다. 다음 적분을 고려하자.

$\int_0^1 \int_x^1 e^{y^2} \, dy \, dx$

내부 적분 $\int_x^1 e^{y^2} \, dy$ 는 초등 함수로 계산할 수 없다. 그러나 적분 순서를 교환하면 다음과 같다.

$\int_0^1 \int_0^y e^{y^2} \, dx \, dy = \int_0^1 y e^{y^2} \, dy = \frac{1}{2}\left(e - 1\right)$

순서를 교환하기 위해서는 적분 영역을 정확하게 파악해야 한다. 원래 영역은 $D = \{(x, y) \mid 0 \leq x \leq 1, \; x \leq y \leq 1\}$ 이며, 이를 $y$ 기준으로 다시 기술하면 $D = \{(x, y) \mid 0 \leq y \leq 1, \; 0 \leq x \leq y\}$ 이 된다.

4. 적분 순서 교환의 체계적 절차

적분 순서를 교환하는 체계적 절차는 다음과 같다.

원래 반복 적분에서 적분 영역 $D$ 를 $xy$ -평면에 그린다.
영역 $D$ 의 경계 곡선을 확인한다.
새로운 순서에 맞게 외부 적분 변수의 범위를 결정한다.
각 외부 변수 값에 대해 내부 변수의 범위를 경계 곡선으로부터 구한다.
영역이 단순하지 않으면 여러 부분으로 분할하여 각각에 대해 반복 적분을 설정한다.

5. 로봇공학에서의 응용 예제

5.1 평면 링크의 면적 계산

반원형 단면을 가진 로봇 링크의 면적을 계산한다. 영역이 $D = \{(x, y) \mid -a \leq x \leq a, \; 0 \leq y \leq \sqrt{a^2 - x^2}\}$ 인 경우 다음과 같다.

$A = \iint_D dA = \int_{-a}^{a} \int_0^{\sqrt{a^2 - x^2}} dy \, dx = \int_{-a}^{a} \sqrt{a^2 - x^2} \, dx = \frac{\pi a^2}{2}$

5.2 작업 공간 내 확률 질량 계산

로봇 말단 장치의 위치 불확실성이 확률 밀도 함수 $p(x, y)$ 로 주어질 때, 특정 영역 $D$ 내에 말단 장치가 위치할 확률은 다음과 같다.

$P = \iint_D p(x, y) \, dA$

2차원 가우시안 분포의 경우, 공분산 행렬 $\boldsymbol{\Sigma}$ 와 평균 $\boldsymbol{\mu}$ 에 대해 다음과 같다.

$p(x, y) = \frac{1}{2\pi \sqrt{\det \boldsymbol{\Sigma}}} \exp\left(-\frac{1}{2}(\mathbf{x} - \boldsymbol{\mu})^T \boldsymbol{\Sigma}^{-1} (\mathbf{x} - \boldsymbol{\mu})\right)$

여기서 $\mathbf{x} = (x, y)^T$ 이다.

6. 이중 적분의 성질

이중 적분은 다음의 기본 성질을 만족한다.

$\iint_D [f(x, y) + g(x, y)] \, dA = \iint_D f(x, y) \, dA + \iint_D g(x, y) \, dA$

$\iint_D c f(x, y) \, dA = c \iint_D f(x, y) \, dA \quad (c \text{는 상수})$

$D = D_1 \cup D_2, \; D_1 \cap D_2 = \emptyset \implies \iint_D f \, dA = \iint_{D_1} f \, dA + \iint_{D_2} f \, dA$

$f(x, y) \leq g(x, y), \; \forall (x, y) \in D \implies \iint_D f \, dA \leq \iint_D g \, dA$

참고 문헌

Stewart, J. (2015). Calculus: Early Transcendentals. 8th ed., Cengage Learning.
Kreyszig, E. (2011). Advanced Engineering Mathematics. 10th ed., Wiley.
Siciliano, B., Sciavicco, L., Villani, L., & Oriolo, G. (2009). Robotics: Modelling, Planning and Control. Springer.
Apostol, T. M. (1967). Calculus, Vol. 2: Multi-Variable Calculus and Linear Algebra. 2nd ed., Wiley.

v 0.1