2020년 7월 20일 월요일

벡터 미적분학(Vector Calculus)

[경고] 아래 글을 읽지 않고 "벡터 미적분학"를 보면 바보로 느껴질 수 있습니다.


3차원 공간에서 움직이는 물체의 자취를 추적하는 가장 효과적인 방법은 좌표계 기반 벡터(vector in coordinate system)이다. [그림 1]과 같은 데카르트 좌표계(Cartesian coordinate system)에서, 물체가 움직여 정지해 있는 현재 위치는 원점에서 $x, y, z$축을 따라 점 $P$까지 가기 위한 변화량으로 정한다.

[그림 1] 데카르트 좌표계의 위치 벡터(출처: wikipedia.org)

이 경우 물체의 위치는 간단히 $P$ = $(x, y, z)$라 표현한다. 하지만 [그림 1]처럼 좌표점을 표현하면 유클리드 기하학(Euclidean geometry)에 나오는 점일 뿐이다. 그래서 점 $P$에 변화되는 특성을 실제로 부여하려면 벡터(vector) 개념을 도입해야 한다. 벡터는 시작점과 끝점이 필요하므로, 시작점을 원점 $(0, 0, 0)$로 고정한다. 이렇게 현재 위치를 원점 대비해 변화한 양으로 표현하면, 점 $P$는 위치 벡터(position vector) $\bar P$ = $(x, y, z) - (0, 0, 0)$ = $(x, y, z)$가 된다. 점 $P$와 위치 벡터 $\bar P$는 동일한 값을 가지지만 내재적 의미는 큰 차이가 있다. 점은 단순한 숫자의 나열일 뿐이지만, 위치 벡터는 직교하는 좌표축으로 변화한 양을 명확히 의미한다.

[그림 2] 매개변수 $t$에 대한 위치 벡터 $\bar r(t)$의 변화(출처: wikipedia.org)

원래 문제로 돌아가서 3차원 공간에서 움직이는 물체를 위치 벡터로 표현한다. 보통 시간에 따라 변화되는 특성을 추적하므로 매개변수는 $t$로 정한다. 그러면 매개변수 $t$의 변화에 대한 임의의 위치 벡터 $\bar r(t)$는 다음처럼 표현할 수 있다.

                  (1)

매개변수 $t$에 대한 $\bar r(t)$의 미분(differentiation)을 구한다.

                  (2)

식 (2)의 의미를 이해하기 위해 [그림 2, 3]을 본다.

[그림 3] 단위 접선과 법선 벡터(출처: wikipedia.org)

[그림 2]에서 $\bar r(t)$는 $t$에 따라 변하는 곡선(curve)을 표현한다. 식 (2)의 우변처럼 벡터의 차를 정하면, $d \bar t(t)$는 $\bar r(t)$에서 시작해 $\bar r(t + \Delta t)$로 향하는 벡터가 된다. 차분 $\Delta t$가 $0$으로 가면, 이 벡터는 [그림 3]처럼 곡선에 접하는 접선 벡터(tangent vector)이다. 따라서 단위 접선 벡터(unit tangent vector) $\hat T$를 다음처럼 공식화할 수 있다.

                  (3)

여기서 $|\bar A|$ = $\sqrt{\bar A \cdot \bar A}$, $\bar A \cdot \bar A$는 벡터의 내적(inner product), $\hat T$의 크기인 $|\hat T|$는 $t$에 관계없이 항상 $1$이다. 식 (3)에 의해 $\hat T$의 크기를 미분하면 다음처럼 항상 $0$이 나온다.

                  (4)

벡터 $\hat T$와 $d\hat T/dt$는 서로 수직이므로, [그림 3]처럼 단위 접선 벡터 $\hat T$에 수직한 단위 법선 벡터(unit normal vector) $\hat N$도 정의한다.

                  (5)

법선 벡터는 접선 벡터의 미분이며 위치 벡터의 2차 미분이기도 하다. 그래서 매개변수가 시간을 의미한다면, 접선 벡터는 속도(velocity), 법선 벡터는 가속도(acceleration)를 뜻한다. 벡터의 외적(outer product)을 이용하면 단위 접선과 법선 벡터에 각각 수직인 단위 종법선 벡터(單位縱法線, unit binormal vector)도 다음처럼 얻을 수 있다.

                  (6)

[그림 4] 3차원에 위치한 단위 접선, 법선, 종법선 벡터(출처: wikipedia.org)

단위 벡터(unit vector) TNB를 3차원 공간의 기저(basis)로 사용하면, 위치 벡터의 궤적인 곡선 $\bar r(t)$에 국소 좌표계(局所座標系, local coordinate system)를 설정할 수 있다. 물론 우리는 3차원 공간의 위치를 표현하기 위해 전역 좌표계(global coordinate system)의 기저인 $\hat x, \hat y, \hat z$를 사용할 수도 있다. 하지만 주어진 문제에 따라 국소 좌표계의 표현식이 해를 더 간략화할 수도 있다. 예를 들면, 식 (29)와 (32)에 정의한 선 적분(line integral)이나 표면 적분(surface integral)을 정의할 때, 국소 좌표계를 도입하면 전체 적분식을 더 명시적이고 간단하게 기술할 수 있다.
현재까지 유도한 결과에서 매개변수를 $t$에서 호의 길이(arc length)인 $s$로 바꾸어본다.

                  (7)

여기서 $t = 0$일 때 $s = 0$이라 가정한다.

[그림 5] 매개변수 $\theta$에 따라 변하는 호의 길이 $s$(출처: wikipedia.org)

매개변수를 $t$에서 $s$로 바꾼 이유는 식 (2)를 $s$에 대한 미분으로 바꾸면 쉽게 알 수 있다.

                  (8)

여기서 $ds/dt$ = $\sqrt{(dx/dt)^2 + (dy/dt)^2 + (dz/dt)^2}$ = $|d \bar r(t)/dt|$이다. 즉 위치 벡터 $\bar r(t)$의 $s$에 대한 미분은 바로 단위 접선 벡터 $\hat T$가 된다. 추가적으로 단위 벡터 NB를 $s$에 대해 미분하면, 식 (8)과 같은 재미있는 결과가 나오는지 차례대로 확인해본다. 먼저 식 (5)를 변형해서 다음처럼 써본다.

                  (9)

식 (9)에서 새롭게 정의한 $\kappa$의 의미를 찾기 위해 위치 벡터를 2차원으로 축소해서 $\bar r(t) = (x(t), y(t))$라 한다. 그러면 다음 관계가 성립한다.

                  (10)

                  (11)

여기서 $d^2 \bar r(t)/ds^2$ = $d/ds [d \bar r(t)/ds]$, $(x', y')$는 2차원 좌표점이면서 벡터, $(x', y')$의 크기는 $\sqrt{x'^2 + y'^2}$이다. 식 (11)을 식 (10)에 대입해 정리하면 $\kappa$를 $\bar r(t)$의 성분으로 표현할 수 있다.

                  (12)

여기서 $(\cdot)'$과 $(\cdot)''$은 각각 $t$에 대한 1차 및 2차 미분이다. 복잡해 보이는 식 (12)에 반지름이 $R$인 원의 매개변수 방정식 $r(t)$ = $(R\cos t, R\sin t)$를 대입해본다.

                  (13)

식 (13)에 의해 $\kappa$는 곡선 $r(t)$에 접하는 원의 반지름에 대한 역수, 즉 곡률(曲率, curvature)임을 알 수 있다. 곡선 $r(t)$와 곡률 $\kappa$의 기하학적 관계는 [그림 6]을 보면 더 쉽게 이해할 수 있다.

[그림 6] 곡선과 곡률의 기하학적 관계(출처: wikipedia.org)

식 (9), (13)과 [그림 6]을 보면, 곡선 $r(t)$를 호의 길이 $s$에 대해 두 번 미분한 벡터의 방향은 단위 법선 벡터 $\hat N$과 같고 벡터의 크기는 곡률 $\kappa$ = $1/R$이 된다. 계속해서 단위 벡터 $\hat N$을 $s$에 대해 미분해보면 다음 결과를 얻을 수 있다.

                  (14)

여기서 $a, b$는 적당한 스칼라(scalar)이다. 스칼라 $a$를 구하기 위해 다음 내적 관계식을 미분한다.

                  (15)

식 (6)을 $s$에 대해 미분한 후 식 (14)를 대입하면, $d \hat B/ ds$를 구할 수 있다.

                  (16)

식 (16)의 결과를 분석해 스칼라 $b$를 비틀림률(torsion) $\tau$로 정한다.

                  (17)

[그림 4]를 보면, 단위 벡터 $\hat B$는 $\hat T$와 $\hat N$이 만드는 평면을 표현한다. 그러면 식 (17)에 의해 비틀림률은 곡선 $\bar r(t)$가 $\hat T$와 $\hat N$이 만드는 평면을 빠져나가는 정도를 뜻한다. 예를 들어 곡선이 특정 평면에만 머물러 있으면,  $\hat T$와 $\hat N$이 만드는 표면의 법선 벡터인 $\hat B$는 모든 평면에서 일정해서 비틀어지지 않고 변화도 없다. 그래서 $s$에 대한 $\hat B$의 변화량인 $d\hat B/ds$는 $0$이 되어서 비틀림률 $\tau$ = $0$이 된다. 반대로 [그림 7]처럼 표면이 휘어져 있으면, 곡선을 따라감에 따라 $\hat B$는 방향이 계속 바뀐다. 즉 $d\hat B/ds$ 혹은 $\tau$에 비례한 만큼 표면이 비틀어져 있다.

[그림 7] 비틀림 막대(출처: wikipedia.org)

이상의 결과를 모두 합치면 유명한 프레네–세레 공식(Frenet–Serret formula)을 다음처럼 공식화할 수 있다.

                  (18)

식 (18)은 행렬(matrix) 형태로 표기할 수도 있다.

                  (19)

단위 벡터 TNB에 대한 행렬 관계식인 식 (19)를 쓰면 $s$에 대한 고차 미분을 다음처럼 표현할 수도 있다.

                  (20)

여기서 $\kappa, \tau$는 $s$에 대해 상수라고 가정한다. 식 (17)에 있는 비틀림률의 정의는 너무 복잡하므로 식 (12)와 같은 간략화된 표현식으로 바꾸어본다. 이 경우는 프레네–세레 공식을 쓰면 확실히 편하다. 먼저 식 (17)에 의해 비틀림률은 $t$에 대한 3차 미분과 분명히 관계있다. 그래서 다음과 같은 1차에서 3차까지 위치 벡터의 미분을 계산한다.

                  (21)

식 (21)의 첫째식과 둘째식에 외적을 적용하면 단위 벡터 $\hat B$를 만들 수 있다.

                  (22)

여기서 $s' = |\bar r'(t)|$이다. 식 (22)에 절대값을 취하면, 어렵게 얻은 식 (12)를 더욱 쉽게 유도할 수 있다. 또한 식 (22)와 식 (21)의 셋째식에 대해 내적을 연산하면, 비틀림률 $\tau$와 곡률 $\kappa$의 관계식을 유도할 수 있다.

                  (23)

식 (12)를 식 (23)에 대입해 정리해서 위치 벡터의 미분으로 표현한 비틀림률 $\tau$에 대한 공식을 만든다.

                  (24)

비틀림률 정의를 위해 식 (17)을 사용했지만, 분명히 식 (24)에서 얻은 $\tau$의 정의가 더 간결하고 분명하다. 식 (24)에 의해 $\tau$는 양과 음이 모두 가능하다. 반면에 $\kappa$는 항상 0이거나 양이다.
벡터 미분학에서 매우 중요한 벡터의 미분 연산자는 전역 좌표계 $(x, y, z)$에서 다음처럼 정의할 수 있다.

                  (25: 구배 연산자)

                  (26: 발산 연산자)

                  (27: 회전 연산자)

서로 직교하는 TBN 기저를 사용하면, 전역 좌표계 $(x, y, z)$에 대비되는 국소 좌표계 $(t, n, b)$를 다음처럼 표현할 수 있다.

                  (28)

여기서 식 (7) 기준으로 보면 $t = s$이다. 단위 접선 벡터 $\hat T$는 곡선 $\bar r(t)$의 접선이므로, $t$는 $\bar r(t)$를 따라가는 매개변수이다. 매개변수 $n$은 단위 법선 벡터 $\hat N$ 방향으로 변한다. 곡선 $\bar r(t)$가 정해지면 $t$에 대해 $\hat N$은 딱 하나만 생긴다. 그러면 고정된 혹은 일정한 $t$에 대해, $\hat N$ 방향으로 $n$을 독립적으로 정할 수 있어서 $n$은 새로운 좌표 성분이 된다. 마찬가지로 정해진 $t$와 $n$에 대해 $\hat B$ 방향으로 새로운 좌표 성분 $b$를 설정할 수 있다. 식 (28)을 이용하면 매개변수 $t$를 따라가는 선 적분(line integral)을 다음처럼 정의할 수 있다.

                  (29)

식 (29)를 참고해서 국소 좌표계의 구배 연산자를 만들어본다. 먼저 완전 미분(exact differential)을 도입해서 함수 $f$의 변화율을 쓰면 다음과 같다.

                  (30)

따라서 국소 좌표계의 구배 연산자를 다음처럼 표현한다.

                  (31)

선 적분인 식 (29)를 2차원으로 확장하면 면적 적분(surface integral)을 구성할 수 있다. 단위 벡터 TNB로 만들 수 있는 평면의 기저 쌍은 $(t, n)$, $(n, b)$, $(b, t)$이다. 그래서 각 기저의 쌍이 만드는 미소 평면을 이용해서 임의의 면적 적분을 다음처럼 기술할 수 있다.

                  (32)

면적 적분을 더 확장하면 국소 좌표계에 대한 체적 적분(volume integral)도 쉽게 정의할 수 있다.

                  (33)

식 (33)처럼 체적 적분은 면적 적분에서 사용한 면적 미분소에 독립적인 방향으로 미분소를 곱한 적분이다. 발산 및 회전 연산자를 국소 좌표계로 확장하기 위해 국소 좌표계의 성분 $(t, n, b)$를 다음처럼 고쳐쓴다.

                  (34)

여기서 $h_i$는 척도 인자(尺度因子, scale factor), $u_i$는 오직 $i$에만 종속인 매개변수, $i$ = $t, n, b$이다. 국소 좌표계 $(t, n, b)$는 직교 좌표계(orthogonal coordinate system)이기 때문에, 주변 좌표 성분에 관계없이[혹은 독립적으로] 특정 방향 성분을 조정할 수 있다. 이를 다른 말로 표현하면, 각 좌표계의 기저는 [그림 4]처럼 서로 직교한다는 말이다. 예를 들면 $(\bar A + \Delta t \hat T)\cdot \hat T$ = $t + \Delta t$, $(\bar A + \Delta t \hat T)\cdot \hat N$ = $n$, $(\bar A + \Delta t \hat T)\cdot \hat B$ = $b$가 성립한다. 여기서 $\bar A$ = $(t, n, b)$이다. 이러한 직교성으로 인해 좌표 성분은 식 (34)와 같이 서로 독립인 매개변수 $u_i$와 척도 인자 $h_i$의 곱으로 구성된다. 척도 인자는 서로 독립인 매개변수와 좌표 성분을 연결해주는 비례 상수이다. 척도 인자가 커질수록, 일정한 $u_i$의 변화에 대해 좌표 성분이 더 크게 변한다. 식 (34)의 관계를 이미 알고 있는 발산 정리(divergence theorem)의 우변에 대입한다.

                  (35)

그러면 다음처럼 닫힌 면적 적분을 체적 적분으로 바꿀 수 있다.

                  (36)

식 (35)의 좌변과 식 (36)을 비교해서 국소 좌표계의 발산 연산자를 다음처럼 공식화할 수 있다.

                  (37)

비슷한 방법으로 국소 좌표계에 대한 회전 연산자를 유도할 수 있다. 증명의 간략화를 위해 벡터가 존재하는 영역을 [그림 8]처럼 단위 벡터 TN이 만드는 평면으로 한정한다.

[그림 8] 단위 벡터 TN이 만드는 평면에 있는 벡터의 회전

스토크스의 정리(Stokes' theorem)에 의해 면적 적분과 닫힌 선 적분은 다음 관계를 가진다.

                  (38)

식 (38)의 좌변에 식 (34)를 대입해서 벡터 성분별로 정리한다.

                  (39)

식 (39)에 나타나는 미분의 부호 특성을 쉽게 이해하려면 [그림 8]이 필수적이다. [그림 8]에서 $t$가 커질 때, 벡터 $\hat n$의 방향으로 인해 $A_n dn$은 (-)에서 (+)로 바뀐다. 그래서 $A_n dn$의 변화량은 미분의 정의[미분은 함수값이 커지는 방향, 즉 기울기가 (+)인 경우가 기준이다.]를 잘 따라가므로 (+)가 된다. 벡터 $\hat t$의 방향으로 인해 [그림 8]에서 변수 $n$이 커지는 경우는 $A_t dt$의 부호가 (+)에서 (-)로 변한다. 따라서 $A_t dt$의 변화량은 음의 기울기라서 식 (39)처럼 (-)를 붙여야 한다. 식 (38)에서 얻은 결과를 TNB가 만들 수 있는 모든 평면으로 확장하면, 국소 좌표계의 회전 연산자는 다음과 같다.

                  (40)

직교하는 국소 좌표계 $(t, n, b)$에 대한 미분 연산자는 텐서(tensor) 이론으로 얻은 직교 좌표계의 미분 연산자와 완전히 동일하다.
국소 좌표계를 위한 기저로 [그림 4]에 있는 TNB 벡터를 사용했지만, 더 간단한 방법으로 국소 좌표계의 기저를 구성할 수도 있다. 평행하지 않은 두 벡터 $\bar u, \bar v$를 고려한다. 그러면 식 (6)과 유사한 벡터의 외적을 이용해 간단한 TNB 기저를 다음처럼 구성할 수 있다.

                  (41)

여기서 $\hat T$는 곡선 $\bar r(t)$가 아닌 주어진 벡터인 $\bar u$로 구성하고, $\hat N$은 $\bar v$ 중에서 $\hat T$에 평행한 성분을 제거해서 만들며, $\hat B$는 $\bar u, \bar v$가 만드는 평면의 법선 벡터(normal vector)로 정의한다. 식 (41)을 이용해서 $\bar v$를 다시 쓰면, $\bar v$ = $(\bar v \cdot \hat T) \hat T + (\bar v \cdot \hat N) \hat N$이다.
특정한 곡선 $\bar r(t)$를 입력 변수로 하는 함수 $f[\bar r(t)]$의 미분은 완전 미분(exact differential)에 의해 다음처럼 공식화된다.

                  (42)

여기서 $\bar r(t)$ = $(x(t), y(t), z(t))$이다. 식 (42)에 의해 함수의 미분 $df/dt$는 구배(gradient)와 곡선의 미분 $d \bar r (t)/dt$ 사이의 내적과 같다. 혹은 $df/dt$는 우리가 계산하고 싶은 방향인 접선 벡터 $\hat T$으로 정사영한 구배에 비례한다. 식 (42)에서 미분을 정의한 곡선을 단순화해서 $\bar r(t)$ = $\bar r_0 + \bar v t$라 한다. 여기서 $\bar r_0$과 $\bar v$는 상수 벡터, 매개변수 $t$에 대해서만 변하는 곡선의 궤적은 3차원 공간의 직선이다. 그러면 식 (42)는 방향 $\bar v$를 따라가는 미분인 방향 도함수(directional derivative)가 된다.

                  (43)

여기서 $\bar v$ = $d \bar r(t) / dt$이다. 방향 도함수 정의에 벡터가 쓰이기는 하지만, 내적 연산에 의해 방향 도함수 자체는 스칼라이며 통상적인 미분과 거의 같다. 다만 방향 도함수는 기존 미분과 다르게 미분 방향을 $\bar v$로 자유롭게 설정할 수 있다. 

[다음 읽을거리]

댓글 13개 :

  1. 식 4번에서 내적의 미분연산자 분배법칙의 증명법에 대해 알려주실 수 있을까요??

    답글삭제
    답글
    1. 식 (4)는 단순한 곱에 대한 미분입니다. 앞에 벡터 한 번 미분하고 뒤에 벡터를 다시 미분해서 더하면 됩니다.

      삭제
    2. 답변 해주셔서 감사합니다 지금 댓글을 봤네요 어떻게 보면 식(4)는 같은 함수T헷 에 대해 곱미분 해준거라고 볼 수 있겠네요??

      삭제
  2. 식(11) 에서 우변 미분이 약간 이해가 안되는데요, dt/ds = 1/sqrt(x'^2+y'^2) 에서 s에 대해 한번 더 미분해서 d2t/ds2 를 구하는건데
    매개변수인 t에 대한 미분이 아니고 s에 대해 미분을 하면 마지막에 미분하고 dt/ds 도 붙어줘야 하는거 아닌가요? s가 t의 함수이므로 t도 s의 함수라고 보고? 헷갈리네요.. 그리고 우변항 분모가 제곱으로 되어 있는데 지수가 3/2 제곱 아닌가요?

    답글삭제
    답글
    1. 매개변수 $s$에 대한 2차 미분을 해보세요. 문제는 없어요.
      그리고 $s$는 다시 $t$에 대해 종속하므로, 연쇄 법칙으로 미분할 수 있어요.
      이해를 위해 식 (11) 밑에 수식을 조금 더 추가했어요.

      삭제
  3. 식(17) 밑에 비틀림률은 T와 N이 만드는 평면을 빠져나가느 정도라구 했는데 식이 N 과 d/ds(T x N)의 내적인데 뒤에 식은 단위벡터 N의 차원이라서 T와 N이 만드는 평면 위에 있는데 어떻게 평면을 빠져나가는 정도가 되는지 머릿속에 잘 안그려집니다

    답글삭제
  4. 식(39)에서 두번째로 갈 때 완전미분하는데 왜 두번째 항이 마이너스가 되나요

    답글삭제
    답글
    1. 잘 이해하려면 스토크스의 정리를 봐야 합니다. 아래 링크 참고하세요.

      https://ghebook.blogspot.com/2010/07/curl.html

      삭제
  5. 14번에서 이해안됩니다
    dN/ds를 어떻게 T와 B의 선형관계식으로 나타낼 수 있는건가요? 마치 비유하자면 x,y,z좌표계에서 x를 c에 대해 미분하면 그건 y와 z의 함수로 나타낼 수 있다는건데 그걸 어케아는것이죠?

    답글삭제
    답글
    1. 3차원 공간이라서 기저(basis)는 3개만 가능합니다. 미분 $d \hat N/ds$는 $\hat N$에 수직이기 때문에, 나머지 2개 기저의 선형 결합으로 $d \hat N/ds$를 표현합니다.

      삭제
  6. 안녕하세요 언제나 감사드립니다
    식 (21)의 3차 미분에서 마지막 항에서 (k`s`+3ks``)s`N이 되어야하는게 아닌가요?
    그리고 질문이 더 있습니다!
    식 (29)에서 왜 A와 dl을 내적하는 지와 dl이 dtT로 쓸 수 있는지 궁금합니다.

    답글삭제
    답글
    1. 1. Unknown님, 지적 정말 감사합니다. 수식에 오타가 있었네요. 다시 고쳤습니다.

      2. 식 (29)의 적분은 벡터 적분학에 흔히 나오는 선 적분입니다. 관련 내용은 식 (29) 위에 있는 단락을 참고해주세요.

      삭제

욕설이나 스팸글은 삭제될 수 있습니다. [전파거북이]는 선플운동의 아름다운 인터넷을 지지합니다.