9.31 짐벌 락(Gimbal Lock) 현상의 원인과 분석

1. 짐벌 락의 정의

짐벌 락(gimbal lock)은 3축 짐벌 메커니즘에서 두 회전 축이 서로 평행하게 정렬되어 한 자유도를 상실하는 현상이다. 수학적으로는 오일러 각으로 매개화된 $SO(3)$ 이 특정 자세에서 매개화의 특이점(singularity)을 가지는 현상이며, 세 독립 각도 중 두 개가 동일한 축 주위의 회전으로 퇴화하여 남은 자유도가 2로 감소한다.

“짐벌 락“이라는 용어는 기계적 짐벌의 물리적 잠김 현상에서 유래하지만, 수학적 의미는 기구가 실제로 잠기는 것이 아니라 특정 자세에서 오일러 각 매개화가 더 이상 회전 공간의 3자유도를 올바르게 표현하지 못하는 것이다.

2. 기계적 짐벌에서의 물리적 짐벌 락

3축 짐벌은 세 개의 링이 각각 수직으로 배치된 회전 축을 제공하는 기계 구조이다. 정상적인 배치에서 세 축은 서로 직교하여 강체가 독립적인 세 방향으로 회전할 수 있다. 그러나 두 번째 링의 회전 각도가 특정 값(보통 $\pm \pi/2$ )에 도달하면 세 번째 링의 축이 첫 번째 링의 축과 평행해진다. 이 상태에서는 첫 번째 링의 회전과 세 번째 링의 회전이 동일한 방향의 회전을 만들어내므로, 3자유도가 2자유도로 퇴화한다.

역사적으로 아폴로 11호의 관성 유도 장치(IMU)에서 짐벌 락이 문제가 될 뻔한 사례가 널리 알려져 있으며, 이 사건은 짐벌 락의 실용적 중요성을 상징적으로 보여준다.

3. 수학적 원인

3.1 $ZYX$ 규약에서의 짐벌 락

$ZYX$ 이동 축 규약에서 피치 각 $\theta = \pi/2$ 일 때, 회전 행렬은

$\mathbf{R}(\psi, \pi/2, \phi) = \mathbf{R}_z(\psi)\mathbf{R}_y(\pi/2)\mathbf{R}_x(\phi)$

이 된다. $\mathbf{R}_y(\pi/2)$ 는 $x$ 축을 $-z$ 축 방향으로, $z$ 축을 $x$ 축 방향으로 사상한다. 따라서 $\mathbf{R}_y(\pi/2)\mathbf{R}_x(\phi)$ 는 원래의 $x$ 축이 이제 $-z$ 방향으로 정렬된 상태에서 그 축 주위의 회전을 의미하는데, 이는 결과적으로 초기 $z$ 축 주위의 회전과 동일하다. 즉, 첫 번째 회전 $\mathbf{R}_z(\psi)$ 와 세 번째 회전 $\mathbf{R}_x(\phi)$ 가 같은 축(월드 $z$ 축) 주위의 회전으로 퇴화한다.

구체적으로 전개하면

$\mathbf{R}(\psi, \pi/2, \phi) = \begin{bmatrix} 0 & \sin(\phi - \psi) & \cos(\phi - \psi) \\ 0 & \cos(\phi - \psi) & -\sin(\phi - \psi) \\ -1 & 0 & 0 \end{bmatrix}$

이다. 이 행렬은 $\phi$ 와 $\psi$ 의 선형 결합 $\phi - \psi$ 에만 의존하며, 두 각도의 개별 값에는 의존하지 않는다.

3.2 $ZYZ$ 규약에서의 짐벌 락

$ZYZ$ 고유 오일러 각 규약에서 중간 각 $\beta = 0$ 일 때, 회전 행렬은

$\mathbf{R}(\alpha, 0, \gamma) = \mathbf{R}_z(\alpha)\mathbf{R}_y(0)\mathbf{R}_z(\gamma) = \mathbf{R}_z(\alpha)\mathbf{R}_z(\gamma) = \mathbf{R}_z(\alpha + \gamma)$

이 된다. 첫 번째와 세 번째 회전이 모두 $z$ 축 주위이므로 두 회전은 하나의 회전 $\mathbf{R}_z(\alpha + \gamma)$ 로 병합된다. 이 경우 $\alpha$ 와 $\gamma$ 의 개별 값은 결정되지 않고, 오직 합 $\alpha + \gamma$ 만이 의미를 가진다.

4. 자코비안의 특이성

짐벌 락의 수학적 본질은 오일러 각 매개화의 자코비안(Jacobian)의 특이성이다. 오일러 각 벡터 $\boldsymbol{\eta} = (\alpha, \beta, \gamma)^T$ 의 시간 미분과 각속도 벡터 $\boldsymbol{\omega}$ 사이의 관계는 다음과 같이 표현된다.

$\boldsymbol{\omega} = \mathbf{E}(\boldsymbol{\eta})\dot{\boldsymbol{\eta}}$

여기서 $\mathbf{E}(\boldsymbol{\eta})$ 는 오일러 각-각속도 자코비안이다. $\mathbf{E}$ 는 대부분의 각도에서 가역이지만, 짐벌 락 특이점에서는 행렬식이 0이 되어 가역이 아니다.

4.1 $ZYX$ 규약의 자코비안 예

$ZYX$ 규약에서

$\mathbf{E}_{ZYX}(\psi, \theta, \phi) = \begin{bmatrix} 1 & 0 & -\sin\theta \\ 0 & \cos\phi & \sin\phi\cos\theta \\ 0 & -\sin\phi & \cos\phi\cos\theta \end{bmatrix}$

이고 행렬식은

$\det(\mathbf{E}_{ZYX}) = \cos\theta$

이다. $\theta = \pm \pi/2$ 에서 $\cos\theta = 0$ 이 되어 자코비안이 특이해진다. 이 특이성이 짐벌 락의 수학적 표현이다.

4.2 역변환의 발산

자코비안의 역 $\mathbf{E}^{-1}$ 이 짐벌 락 근처에서 발산하므로, 각속도를 오일러 각 미분으로 변환하는 공식

$\dot{\boldsymbol{\eta}} = \mathbf{E}^{-1}(\boldsymbol{\eta})\boldsymbol{\omega}$

에서 $\dot{\boldsymbol{\eta}}$ 의 성분이 무한대로 발산할 수 있다. 이는 유한한 각속도가 오일러 각의 무한한 시간 변화율을 요구하는 병적 상황이다.

5. 기하학적 해석

5.1 회전 축의 정렬

짐벌 락의 기하학적 본질은 첫 번째와 세 번째 회전 축의 정렬이다. 중간 회전이 $\pm \pi/2$ (테이트-브라이언) 또는 $0/\pi$ (고유 오일러 각)에 도달하면 두 축이 평행하게 되어 동일한 방향의 회전만을 생성할 수 있다.

5.2 회전 공간의 국소 구조

$SO(3)$ 은 3차원 매니폴드이므로 국소적으로 $\mathbb{R}^3$ 과 동형이다. 그러나 유클리드 매개화는 매니폴드의 전역 구조를 올바르게 포착하지 못한다. 이는 “털공 정리(hairy ball theorem)“와 관련된 위상적 결과이며, $S^2$ (구면)을 덮는 연속 단위 벡터 장이 존재하지 않는다는 사실과 유사하다. $SO(3)$ 역시 세 매개변수로 전역적으로 부드럽게 덮을 수 없다.

6. 짐벌 락의 위치

12가지 오일러 각 규약 모두 짐벌 락 특이점을 가지며, 특이점의 위치는 규약에 따라 다르다.

6.1 테이트-브라이언 각 (예: $ZYX$ , $XYZ$ )

중간 각이 $\pm \pi/2$ 일 때 짐벌 락이 발생한다. 이는 기체의 자세가 수직 상승 또는 수직 하강에 해당한다.

6.2 고유 오일러 각 (예: $ZYZ$ , $ZXZ$ )

중간 각이 $0$ 또는 $\pi$ 일 때 짐벌 락이 발생한다. 이는 첫 번째와 세 번째 회전 축이 동일한 방향으로 정렬된 자세이다.

7. 짐벌 락의 영향

7.1 수치 계산의 불안정성

짐벌 락 근처에서 오일러 각 추출이 수치적으로 불안정해진다. 작은 섭동이 오일러 각에 큰 변화를 일으키므로, 센서 잡음이나 부동 소수점 오차가 증폭된다.

7.2 궤적 생성의 단절

오일러 각으로 매개화된 궤적이 짐벌 락을 통과하면 각도가 불연속적으로 점프할 수 있다. 이로 인해 보간된 자세가 의도와 다른 경로를 따를 수 있다.

7.3 제어 시스템의 오류

자세 제어기가 오일러 각 기반으로 설계된 경우, 짐벌 락 근처에서 제어 입력이 비현실적으로 커지거나 발산할 수 있다. 이는 실제 구동기의 포화를 초래하여 제어 성능을 훼손한다.

8. 표면적 현상 vs. 본질적 한계

짐벌 락은 물리적 회전 자체의 한계가 아니라 특정 매개화의 한계임을 명확히 이해해야 한다. 회전 행렬이나 쿼터니언 수준에서는 짐벌 락에 해당하는 자세가 특별하지 않으며, 정상적으로 작동한다. 짐벌 락은 오일러 각 매개화가 $SO(3)$ 을 전역적으로 덮을 수 없다는 위상적 사실의 국소적 표현이다.

9. 로봇 공학적 고려

9.1 작업 영역의 설계

짐벌 락 특이점이 빈번히 발생하지 않도록 로봇의 작업 영역을 설계한다. 예를 들어 지상 이동 로봇이 수직 자세에 도달하지 않으면 $ZYX$ 규약의 짐벌 락을 회피할 수 있다.

9.2 대안 매개화의 선택

내부 계산에서 쿼터니언, 축-각도, 리 대수 표현을 사용하면 짐벌 락이 완전히 제거된다. 오일러 각은 사용자 인터페이스와 해석 출력에만 사용한다.

9.3 예외 처리

부득이하게 오일러 각을 사용하는 경우, 짐벌 락 근처에서 특별한 예외 처리를 구현한다. 자코비안이 특이해지면 의사 역행렬(pseudoinverse) 또는 감쇠 최소 제곱(damped least squares)으로 수치적 안정성을 확보한다.

10. 참고 문헌

Shuster, M. D. (1993). “A Survey of Attitude Representations.” Journal of the Astronautical Sciences, 41(4), 439–517.
Diebel, J. (2006). “Representing Attitude: Euler Angles, Unit Quaternions, and Rotation Vectors.” Stanford University Technical Report.
Stuelpnagel, J. (1964). “On the Parametrization of the Three-Dimensional Rotation Group.” SIAM Review, 6(4), 422–430.
Siciliano, B., Sciavicco, L., Villani, L., & Oriolo, G. (2009). Robotics: Modelling, Planning and Control. Springer.
Craig, J. J. (2018). Introduction to Robotics: Mechanics and Control (4th ed.). Pearson.

version: 1.0