32.10 기하학적 자코비안과 해석적 자코비안의 관계
기하학적 자코비안(geometric Jacobian)과 해석적 자코비안(analytical Jacobian)은 동일한 매니퓰레이터의 속도 기구학을 서로 다른 수학적 언어로 표현하는 두 자코비안이다. 두 자코비안은 엔드 이펙터 방향을 물리적 각속도로 기술하느냐, 또는 최소 매개변수화된 방향 표현의 시간 미분으로 기술하느냐의 차이에서 비롯되며, 그 관계는 방향 매개변수의 종류에 의존하는 표현 행렬(representation matrix)을 매개로 명시적으로 도출된다. 본 절에서는 두 자코비안의 정의 재확인, 블록 단위 관계식, 표현 행렬의 역할, 특이점의 비교, 수치적 안정성, 응용에 따른 선택, 상호 변환 알고리즘을 학술적으로 다룬다.
1. 두 자코비안의 정의 재확인
기하학적 자코비안 \mathbf{J}_g(\vec{q}) \in \mathbb{R}^{6 \times n}은 관절 속도 \dot{\vec{q}}로부터 엔드 이펙터의 물리적 선속도 \vec{v}와 물리적 각속도 \vec{\omega}를 결합한 6차원 속도 벡터로의 선형 사상을 정의한다.
\begin{bmatrix} \vec{v} \\ \vec{\omega} \end{bmatrix} = \mathbf{J}_g(\vec{q}) \, \dot{\vec{q}}
해석적 자코비안 \mathbf{J}_a(\vec{q}) \in \mathbb{R}^{(3+k) \times n}은 관절 속도로부터 엔드 이펙터 위치의 시간 미분 \dot{\vec{p}}와 방향 매개변수의 시간 미분 \dot{\vec{\phi}}를 결합한 벡터로의 선형 사상을 정의한다.
\begin{bmatrix} \dot{\vec{p}} \\ \dot{\vec{\phi}} \end{bmatrix} = \mathbf{J}_a(\vec{q}) \, \dot{\vec{q}}
양자의 위치 관련 출력은 동일한 \dot{\vec{p}} = \vec{v}이지만, 방향 관련 출력이 다르다는 점이 두 자코비안의 구조적 차이를 낳는다.
2. 각속도와 방향 매개변수 미분 사이의 변환 관계
각속도 \vec{\omega}와 방향 매개변수의 시간 미분 \dot{\vec{\phi}} 사이에는 방향 매개변수의 선택에 의존하는 표현 행렬 \mathbf{E}(\vec{\phi}) \in \mathbb{R}^{3 \times k}를 통해 다음 관계가 성립한다.
\vec{\omega} = \mathbf{E}(\vec{\phi}) \, \dot{\vec{\phi}}
이 관계의 학술적 근거는 방향 매개변수가 특정 규약의 기본 회전(elementary rotation)들의 순차적 합성으로 엔드 이펙터의 방향을 기술하며, 각 기본 회전이 각속도에 기여하는 순시 성분들의 누적합이 물리적 각속도로 나타난다는 점이다. 세 축 회전의 각속도 기여를 누적하여 얻어지는 \mathbf{E}(\vec{\phi})는 방향 매개변수 규약마다 상이하다.
k = 3인 표준적 3매개변수 표현에서 \mathbf{E}(\vec{\phi})가 가역인 구성에서는 역관계가 성립한다.
\dot{\vec{\phi}} = \mathbf{E}^{-1}(\vec{\phi}) \, \vec{\omega}
그러나 \det \mathbf{E}(\vec{\phi}) = 0이 되는 특정 방향 구성에서는 이 역사상이 정의되지 않으며, 이는 방향 표현의 특이점(representation singularity)으로 불린다. 대표적 예로 ZYZ 오일러 각에서는 \theta = 0, \pi, ZYX 각(RPY)에서는 \beta = \pm \pi/2에서 \mathbf{E}가 특이해진다.
3. 블록 단위 관계식
두 자코비안의 관계는 위치 블록과 방향 블록으로 분리하여 기술할 때 가장 명료하다.
3.1 위치 블록의 일치
엔드 이펙터 위치 \vec{p}는 3차원 유클리드 공간의 벡터로서 매개변수화 선택과 무관하게 정의되므로, 위치 관련 속도 성분에 대응하는 두 자코비안의 블록은 일치한다.
\mathbf{J}_{a,p}(\vec{q}) = \mathbf{J}_{g,v}(\vec{q}) = \frac{\partial \vec{p}}{\partial \vec{q}}
이는 \dot{\vec{p}} = \vec{v}라는 등식이 항상 성립한다는 물리적 사실의 직접적 귀결이다.
32.10.3.2 방향 블록의 표현 행렬을 통한 연결
방향 관련 블록은 표현 행렬 \mathbf{E}(\vec{\phi})의 역행렬을 통해 연결된다. 관계식 \dot{\vec{\phi}} = \mathbf{E}^{-1}(\vec{\phi}) \, \vec{\omega} = \mathbf{E}^{-1}(\vec{\phi}) \, \mathbf{J}_{g,\omega} \dot{\vec{q}}로부터 다음이 얻어진다.
\mathbf{J}_{a,\phi}(\vec{q}) = \mathbf{E}^{-1}(\vec{\phi}(\vec{q})) \, \mathbf{J}_{g,\omega}(\vec{q})
3.2 전체 블록 관계식
두 블록을 결합하면 해석적 자코비안과 기하학적 자코비안의 관계가 블록 대각 변환의 형태로 표현된다.
\mathbf{J}_a(\vec{q}) = \mathbf{T}_a(\vec{\phi}) \, \mathbf{J}_g(\vec{q}), \qquad \mathbf{T}_a(\vec{\phi}) = \begin{bmatrix} \mathbf{I}_3 & \mathbf{0} \\ \mathbf{0} & \mathbf{E}^{-1}(\vec{\phi}) \end{bmatrix}
이 변환 행렬 \mathbf{T}_a(\vec{\phi}) \in \mathbb{R}^{(3+k) \times 6}은 방향 매개변수화의 규약에 의해 결정되며, 위치 블록은 항등 사상, 방향 블록은 \mathbf{E}^{-1}(\vec{\phi})에 의한 사상으로 구성된다. 역으로 해석적 자코비안으로부터 기하학적 자코비안의 방향 블록을 얻으려면 \mathbf{J}_{g,\omega} = \mathbf{E}(\vec{\phi}) \, \mathbf{J}_{a,\phi}로 변환하면 된다.
32.10.4 표현 행렬의 구체적 형태
표현 행렬 \mathbf{E}(\vec{\phi})는 방향 매개변수 규약에 따라 다음과 같이 구체화된다.
ZYZ 오일러 각 (\phi, \theta, \psi)의 경우:
\mathbf{E}_{\text{ZYZ}}(\phi, \theta, \psi) = \begin{bmatrix} 0 & -\sin\phi & \cos\phi \sin\theta \\ 0 & \cos\phi & \sin\phi \sin\theta \\ 1 & 0 & \cos\theta \end{bmatrix}, \qquad \det \mathbf{E}_{\text{ZYZ}} = \sin\theta
ZYX 각(RPY) (\alpha, \beta, \gamma)의 경우(공간 고정축 규약):
\mathbf{E}_{\text{ZYX}}(\alpha, \beta, \gamma) = \begin{bmatrix} \cos\alpha \cos\beta & -\sin\alpha & 0 \\ \sin\alpha \cos\beta & \cos\alpha & 0 \\ -\sin\beta & 0 & 1 \end{bmatrix}, \qquad \det \mathbf{E}_{\text{ZYX}} = \cos\beta
다른 오일러 각 규약(XYZ, ZXZ, XYX 등)에서도 유사한 형태의 표현 행렬이 유도된다. 규약을 명시하지 않은 채 해석적 자코비안을 사용하면 수치적 일관성이 유지되지 않으므로, 학술적·실무적 응용에서 규약의 명시적 선언은 필수적이다.
32.10.5 두 자코비안의 행렬식 관계와 특이점 구분
두 자코비안의 관계식으로부터 행렬식 사이의 직접적 관계가 도출된다. \mathbf{T}_a의 블록 대각 구조에 따라
\det \mathbf{T}_a(\vec{\phi}) = \det \mathbf{I}_3 \cdot \det \mathbf{E}^{-1}(\vec{\phi}) = \frac{1}{\det \mathbf{E}(\vec{\phi})}
이므로, n = 6 정방 자코비안의 경우 다음 관계가 성립한다.
\det \mathbf{J}_a(\vec{q}) = \frac{\det \mathbf{J}_g(\vec{q})}{\det \mathbf{E}(\vec{\phi}(\vec{q}))}
이 관계로부터 두 자코비안의 특이점이 다음과 같이 분해된다.
기구학적 특이점(kinematic singularity)은 \det \mathbf{J}_g = 0이 되는 구성이며, 이는 매니퓰레이터의 본질적 기구학 구조에 기인하고 방향 매개변수의 선택과 무관하다. 이러한 구성에서는 두 자코비안이 모두 특이하다.
표현 특이점(representation singularity)은 \det \mathbf{E}(\vec{\phi}) = 0이 되는 방향 구성이며, 이는 매개변수화 규약에 의해서만 발생한다. 이러한 구성에서 해석적 자코비안은 특이하지만, 기하학적 자코비안은 정상일 수 있다.
따라서 기하학적 자코비안은 오직 기구학적 특이점에서만 특이해지는 반면, 해석적 자코비안은 기구학적 특이점과 표현 특이점의 합집합에서 특이해진다. 수치적 해석과 제어 설계에서는 이 두 유형의 특이점을 분리하여 다루는 것이 중요하다.
32.10.6 수치적 안정성
기하학적 자코비안은 방향 매개변수화에 의존하지 않으므로 표현 특이점의 영향을 받지 않으며, 수치적으로 더 안정적이다. 해석적 자코비안은 짐벌 잠김 근방에서 \det \mathbf{E} \to 0에 따라 \mathbf{E}^{-1}의 성분이 발산하므로 방향 부분의 성분이 수치적으로 불안정해진다.
실무적 구현에서는 다음과 같은 대응이 취해진다. 첫째, 방향 추종이 표현 특이점 근방에서 이루어져야 하는 경우 해석적 자코비안 대신 기하학적 자코비안과 각속도 기반 오차 표현(축-각 또는 사원수 오차)을 사용한다. 둘째, 해석적 자코비안을 사용할 때는 표현 특이점을 회피하는 작업 공간 계획을 수행한다. 셋째, 여러 방향 매개변수 규약을 상황에 따라 전환하여 각 규약의 특이점을 회피한다. 넷째, 사원수 또는 축-각 표현처럼 특이점이 없거나 제한된 표현을 활용하여 특이점의 영향을 최소화한다.
32.10.7 응용에 따른 선택 지침
두 자코비안 중 어느 것을 사용할지는 응용의 성격에 따라 결정된다.
속도 수준 작업 공간 제어: 엔드 이펙터의 선속도와 각속도를 직접 제어 변수로 사용하는 작업 공간 제어, 해상도 모션 제어, 임피던스 제어 등의 응용에서는 기하학적 자코비안이 자연스럽다. 물리적 속도 벡터가 제어 신호에 직접 대응하므로 해석과 설계가 명료하다.
방향 매개변수 기반 궤적 추종: 방향 경로가 오일러 각 또는 RPY 각의 시간 함수로 주어지는 궤적 추종 제어에서는 해석적 자코비안이 직접적 수학적 도구를 제공한다. 방향 매개변수의 오차를 직접 자코비안으로 변환하여 관절 속도 명령을 얻을 수 있다.
힘 제어와 임피던스 제어: 힘-속도 쌍대성이 요구되는 제어 설계에서는 기하학적 자코비안이 선호된다. 자코비안 전치를 통한 힘 사상 \vec{\tau} = \mathbf{J}_g^\top \vec{F}가 물리적 의미가 명료하기 때문이다.
스크류 이론과 리 군 기반 해석: 현대적 기구학의 스크류 이론, 트위스트 표현, 리 대수 기반 해석에서는 기하학적 자코비안이 표준적 수학 언어이다.
교육적·해석적 용도: 학부 수준 교육과 초기 분석에서는 직관적 이해가 용이한 해석적 자코비안이 교육적 가치를 지니며, 전문적 해석에는 기하학적 자코비안이 활용된다.
32.10.8 상호 변환 알고리즘
기하학적 자코비안과 해석적 자코비안 사이의 변환은 다음 절차로 수행된다.
기하학적에서 해석적으로의 변환:
- 순기구학의 회전 행렬 {}^0\mathbf{R}_n(\vec{q})로부터 선택된 규약에 따라 방향 매개변수 \vec{\phi}(\vec{q})를 추출한다.
- 해당 규약의 표현 행렬 \mathbf{E}(\vec{\phi})를 구성한다.
- 위치 블록은 \mathbf{J}_{a,p} = \mathbf{J}_{g,v}로 그대로 사용한다.
- 방향 블록은 선형 연립 방정식 \mathbf{E}(\vec{\phi}) \mathbf{J}_{a,\phi} = \mathbf{J}_{g,\omega}의 해로 구한다. 명시적 \mathbf{E}^{-1} 계산 대신 연립 방정식 풀이가 수치적으로 안정적이다.
해석적에서 기하학적으로의 변환:
- 위치 블록은 \mathbf{J}_{g,v} = \mathbf{J}_{a,p}로 그대로 사용한다.
- 방향 블록은 \mathbf{J}_{g,\omega} = \mathbf{E}(\vec{\phi}) \, \mathbf{J}_{a,\phi}로 직접 계산한다. 이 변환에는 \mathbf{E}의 가역성이 요구되지 않으므로 표현 특이점에서도 수치적으로 안정적이다.
변환의 계산 비용은 방향 매개변수 추출과 표현 행렬 구성의 \mathcal{O}(1) 고정 비용과 3 \times 3 행렬 연산의 \mathcal{O}(n) 비용의 합으로, 실시간 제어에 충분히 효율적이다.
32.10.9 학술적·문헌적 선호의 경향
현대 로봇 공학 교과서와 연구 문헌은 기하학적 자코비안을 중심적 학술 객체로 다루는 경향이 뚜렷하다. 그 이유는 다음과 같다. 첫째, 방향 매개변수화에 독립적이므로 수학적 일반성이 높다. 둘째, 수치적 안정성이 우수하다. 셋째, 스크류 이론, 리 군-리 대수, 수반 사상 등 현대적 수학 구조와 자연스럽게 결합한다. 넷째, 힘-속도 쌍대성의 표현이 명료하다.
반면 해석적 자코비안은 초기 로봇 공학 문헌과 특정 응용 맥락에서 널리 사용되어 왔으며, 현재에도 방향 매개변수 기반 궤적 설계와 교육적 문맥에서 유용하다. 학술적 이해를 위해 두 자코비안의 정의, 관계, 장단점을 모두 파악하고 응용에 적절한 형태를 선택하는 능력이 요구된다.
32.10.10 학술적 의의
본 절에서 정립한 기하학적 자코비안과 해석적 자코비안의 관계는 로봇 속도 기구학에서 다음의 학술적 의의를 가진다. 첫째, 표현 행렬 \mathbf{E}(\vec{\phi})를 매개로 한 명시적 블록 변환 공식 \mathbf{J}_a = \text{diag}(\mathbf{I}_3, \mathbf{E}^{-1}) \mathbf{J}_g는 두 자코비안의 체계적 연결을 제공한다. 둘째, 행렬식 분해 \det \mathbf{J}_a = \det \mathbf{J}_g / \det \mathbf{E}는 기구학적 특이점과 표현 특이점의 독립적 분석을 가능하게 한다. 셋째, 두 자코비안의 수치적 안정성과 응용 적합성의 차이는 제어 설계와 궤적 계획에서 자코비안 선택의 공학적 근거를 제공한다. 넷째, 두 자코비안 사이의 상호 변환 알고리즘은 응용에 따라 표현을 전환할 수 있는 실무적 도구를 제공한다. 다섯째, 본 절의 내용은 후속 절의 오일러 각 기반·RPY 각 기반 해석적 자코비안 변환 논의의 기초가 된다.
출처
- Siciliano, B., Sciavicco, L., Villani, L., and Oriolo, G., Robotics: Modelling, Planning and Control, Springer, 2009.
- Sciavicco, L. and Siciliano, B., Modelling and Control of Robot Manipulators, 2nd edition, Springer, 2000.
- Spong, M. W., Hutchinson, S., and Vidyasagar, M., Robot Modeling and Control, 2nd edition, Wiley, 2020.
- Craig, J. J., Introduction to Robotics: Mechanics and Control, 4th edition, Pearson, 2018.
- Lynch, K. M. and Park, F. C., Modern Robotics: Mechanics, Planning, and Control, Cambridge University Press, 2017.
- Murray, R. M., Li, Z., and Sastry, S. S., A Mathematical Introduction to Robotic Manipulation, CRC Press, 1994.
- Yoshikawa, T., Foundations of Robotics: Analysis and Control, MIT Press, 1990.
- Shuster, M. D., “A survey of attitude representations”, Journal of the Astronautical Sciences, Vol. 41, No. 4, pp. 439–517, 1993.
버전
- 문서 버전: 1.1
- 작성일: 2026-04-19