9.37 반대칭 행렬(Skew-Symmetric Matrix)과 외적 연산

9.37 반대칭 행렬(Skew-Symmetric Matrix)과 외적 연산

1. 반대칭 행렬의 정의

실수 정방 행렬 \mathbf{S} \in \mathbb{R}^{n \times n}이 다음 조건을 만족할 때 반대칭 행렬(skew-symmetric matrix) 또는 **반대칭 행렬(anti-symmetric matrix)**이라 한다.

\mathbf{S}^T = -\mathbf{S}

이는 \mathbf{S}의 원소가 다음을 만족함을 의미한다.

s_{ij} = -s_{ji}

특히 대각선 원소는 s_{ii} = -s_{ii}이므로 s_{ii} = 0이다.

2. 차원 반대칭 행렬의 형태

3차원 반대칭 행렬 \mathbf{S} \in \mathbb{R}^{3 \times 3}은 세 개의 독립 성분으로 결정된다.

\mathbf{S} = \begin{bmatrix}0 & -s_3 & s_2 \\ s_3 & 0 & -s_1 \\ -s_2 & s_1 & 0\end{bmatrix}

여기서 (s_1, s_2, s_3)는 세 독립 매개변수이다. 이 세 매개변수는 정확히 3차원 벡터 \mathbf{s} = (s_1, s_2, s_3)^T의 성분과 대응된다. 따라서 3차원 반대칭 행렬 공간과 \mathbb{R}^3는 자연스러운 동형 관계를 갖는다.

3. 외적 연산과의 동일시

3차원 벡터 \mathbf{a} = (a_1, a_2, a_3)^T와 임의 벡터 \mathbf{b} = (b_1, b_2, b_3)^T의 외적(cross product)은 다음과 같이 정의된다.

\mathbf{a} \times \mathbf{b} = \begin{bmatrix}a_2 b_3 - a_3 b_2 \\ a_3 b_1 - a_1 b_3 \\ a_1 b_2 - a_2 b_1\end{bmatrix}

이 외적은 다음의 행렬 곱으로 표현할 수 있다.

\mathbf{a} \times \mathbf{b} = [\mathbf{a}]_\times \mathbf{b}

여기서

[\mathbf{a}]_\times = \begin{bmatrix}0 & -a_3 & a_2 \\ a_3 & 0 & -a_1 \\ -a_2 & a_1 & 0\end{bmatrix}

\mathbf{a}반대칭 행렬 표현 또는 **교차곱 행렬(cross product matrix)**이라 한다.

4. 해트 연산자와 벡터 연산자

두 가지 연산자를 정의한다.

4.1 해트 연산자 (\cdot)^\wedge

벡터에서 반대칭 행렬로의 사상이다.

\mathbf{a}^\wedge = [\mathbf{a}]_\times

4.2 벡터 연산자 (\cdot)^\vee

반대칭 행렬에서 벡터로의 역사상이다.

([\mathbf{a}]_\times)^\vee = \mathbf{a}

두 연산자는 서로의 역함수이다. 이들은 3차원 반대칭 행렬 공간과 \mathbb{R}^3 사이의 선형 동형 사상을 제공한다.

5. 반대칭 행렬의 연산적 성질

5.1 합과 스칼라 곱

반대칭 행렬의 합과 스칼라 배는 다시 반대칭 행렬이다.

[\mathbf{a}]_\times + [\mathbf{b}]_\times = [\mathbf{a} + \mathbf{b}]_\times

c[\mathbf{a}]_\times = [c\mathbf{a}]_\times

이는 해트 연산자가 선형임을 보인다.

5.2 전치

[\mathbf{a}]_\times^T = -[\mathbf{a}]_\times이므로 전치는 부호를 반전시킨다.

5.3 거듭제곱

3차원 반대칭 행렬의 거듭제곱은 다음의 패턴을 따른다.

[\mathbf{a}]_\times^2 = \mathbf{a}\mathbf{a}^T - \lVert \mathbf{a} \rVert^2 \mathbf{I}

[\mathbf{a}]_\times^3 = -\lVert \mathbf{a} \rVert^2 [\mathbf{a}]_\times

\hat{\mathbf{u}}가 단위 벡터인 경우 [\hat{\mathbf{u}}]_\times^3 = -[\hat{\mathbf{u}}]_\times로 간소화되며, 이 항등식은 로드리게스 공식의 핵심이 된다.

5.4 행렬식과 고유값

3차원 반대칭 행렬의 행렬식은 항상 0이다.

\det([\mathbf{a}]_\times) = 0

이는 [\mathbf{a}]_\times\mathbf{a} = \mathbf{a} \times \mathbf{a} = \mathbf{0}으로부터 \mathbf{a}가 영공간에 속함을 통해 확인된다. 고유값은 \{0, i\lVert \mathbf{a} \rVert, -i\lVert \mathbf{a} \rVert\}이며, 영이 아닌 두 고유값은 켤레 복소수이다.

6. 외적의 대수적 성질

반대칭 행렬 표현을 이용하면 외적의 여러 성질을 행렬 연산으로 환원할 수 있다.

6.1 반대칭성

\mathbf{a} \times \mathbf{b} = -\mathbf{b} \times \mathbf{a}

이는 [\mathbf{a}]_\times\mathbf{b} = -[\mathbf{b}]_\times\mathbf{a}에 대응한다.

6.2 분배 법칙

\mathbf{a} \times (\mathbf{b} + \mathbf{c}) = \mathbf{a} \times \mathbf{b} + \mathbf{a} \times \mathbf{c}

행렬 곱의 분배 법칙이다.

6.3 야코비 항등식

\mathbf{a} \times (\mathbf{b} \times \mathbf{c}) + \mathbf{b} \times (\mathbf{c} \times \mathbf{a}) + \mathbf{c} \times (\mathbf{a} \times \mathbf{b}) = \mathbf{0}

이는 리 대수의 야코비 항등식에 해당하며, 반대칭 행렬의 교환자 [\mathbf{A}, \mathbf{B}] = \mathbf{A}\mathbf{B} - \mathbf{B}\mathbf{A}와 외적의 대응 관계로 이해된다.

6.4 라그랑주 항등식(BAC-CAB)

\mathbf{a} \times (\mathbf{b} \times \mathbf{c}) = \mathbf{b}(\mathbf{a}^T\mathbf{c}) - \mathbf{c}(\mathbf{a}^T\mathbf{b})

이 공식은 반대칭 행렬 표현으로 [\mathbf{a}]_\times[\mathbf{b}]_\times의 전개로부터 유도된다.

7. \mathfrak{so}(3) 리 대수와의 관계

3차원 반대칭 행렬의 집합은 리 대수 \mathfrak{so}(3)을 이룬다.

\mathfrak{so}(3) = \{\mathbf{S} \in \mathbb{R}^{3 \times 3} \mid \mathbf{S}^T = -\mathbf{S}\}

이 공간은 3차원 벡터 공간이며, 리 브라켓은 행렬 교환자로 정의된다.

[\mathbf{A}, \mathbf{B}] = \mathbf{A}\mathbf{B} - \mathbf{B}\mathbf{A}

해트 연산자를 통해 \mathfrak{so}(3)\mathbb{R}^3가 리 대수로서 동형임이 다음의 관계로 나타난다.

[[\mathbf{a}]_\times, [\mathbf{b}]_\times] = [\mathbf{a} \times \mathbf{b}]_\times

즉, 반대칭 행렬의 교환자와 벡터의 외적이 정확히 대응한다. 이는 (\mathbb{R}^3, \times)가 리 대수이며 \mathfrak{so}(3)과 동형임을 보여준다.

8. 표준 기저

\mathfrak{so}(3)의 표준 기저는 다음과 같다.

\mathbf{J}_x = [\hat{\mathbf{e}}_x]_\times = \begin{bmatrix}0 & 0 & 0 \\ 0 & 0 & -1 \\ 0 & 1 & 0\end{bmatrix}

\mathbf{J}_y = [\hat{\mathbf{e}}_y]_\times = \begin{bmatrix}0 & 0 & 1 \\ 0 & 0 & 0 \\ -1 & 0 & 0\end{bmatrix}

\mathbf{J}_z = [\hat{\mathbf{e}}_z]_\times = \begin{bmatrix}0 & -1 & 0 \\ 1 & 0 & 0 \\ 0 & 0 & 0\end{bmatrix}

이 세 기저 행렬은 다음의 교환 관계를 만족한다.

[\mathbf{J}_x, \mathbf{J}_y] = \mathbf{J}_z, \quad [\mathbf{J}_y, \mathbf{J}_z] = \mathbf{J}_x, \quad [\mathbf{J}_z, \mathbf{J}_x] = \mathbf{J}_y

이 관계는 외적의 순환 규칙과 정확히 일치한다.

9. 회전 운동학에서의 응용

9.1 각속도의 반대칭 행렬 표현

각속도 벡터 \boldsymbol{\omega}의 회전 행렬 미분에 대한 작용은 반대칭 행렬로 표현된다.

\dot{\mathbf{R}} = [\boldsymbol{\omega}]_\times \mathbf{R}

또는 본체 좌표계의 각속도의 경우

\dot{\mathbf{R}} = \mathbf{R}[\boldsymbol{\omega}_{\text{body}}]_\times

이 관계는 회전 행렬이 단위 구면에서의 운동을 기술하는 \mathfrak{so}(3) 생성기에 의해 생성됨을 나타낸다.

9.2 모멘트와 외력의 교차

강체 동역학에서 외력 \mathbf{f}에 의한 모멘트는 \mathbf{r} \times \mathbf{f} = [\mathbf{r}]_\times\mathbf{f}로 표현된다. 반대칭 행렬은 힘을 모멘트로 변환하는 선형 연산자이다.

9.3 자이로스코프 효과

자이로스코프에 작용하는 회전 운동의 커플링 항은 반대칭 행렬 곱의 형태로 나타난다. 예를 들어 각운동량 \mathbf{L}과 각속도 \boldsymbol{\omega}의 상호 작용이 \boldsymbol{\omega} \times \mathbf{L} = [\boldsymbol{\omega}]_\times\mathbf{L}의 형태로 동역학 방정식에 등장한다.

10. 일반 차원의 반대칭 행렬

n차원 반대칭 행렬 공간의 차원은 n(n-1)/2이다.

  • n = 2: 차원 1 (단일 매개변수의 평면 회전)
  • n = 3: 차원 3 (3차원 벡터 \mathbb{R}^3와 동형)
  • n = 4: 차원 6

3차원이 특별한 이유는 차원이 벡터 공간 \mathbb{R}^3의 차원과 일치하기 때문이며, 이는 외적이 3차원에서만 자연스럽게 정의되는 이유이다. 고차원에서는 반대칭 행렬이 벡터가 아니라 더 일반적인 “바이벡터(bivector)“로 해석된다.

11. 로봇 공학적 활용 정리

반대칭 행렬과 외적 연산은 로봇 공학에서 다음과 같이 광범위하게 사용된다.

  • 로드리게스 공식: 축-각도로부터 회전 행렬 계산
  • 각속도와 회전 행렬 미분: 자세 운동학의 핵심 관계
  • 지수 사상: \mathfrak{so}(3)에서 SO(3)으로의 변환
  • 비선형 필터링: 확장 칼만 필터에서의 회전 오차 전파
  • 힘-모멘트 변환: 강체 동역학과 제어
  • 기구학 자코비안: 관절 각속도와 말단 장치 속도 사이의 변환

12. 참고 문헌

  • Murray, R. M., Li, Z., & Sastry, S. S. (1994). A Mathematical Introduction to Robotic Manipulation. CRC Press.
  • Hall, B. C. (2015). Lie Groups, Lie Algebras, and Representations (2nd ed.). Springer.
  • Siciliano, B., Sciavicco, L., Villani, L., & Oriolo, G. (2009). Robotics: Modelling, Planning and Control. Springer.
  • Barfoot, T. D. (2017). State Estimation for Robotics. Cambridge University Press.
  • Sola, J., Deray, J., & Atchuthan, D. (2018). “A Micro Lie Theory for State Estimation in Robotics.” arXiv:1812.01537.

version: 1.0