32.19 공간 자코비안(Space Jacobian)과 물체 자코비안(Body Jacobian)
공간 자코비안(space Jacobian)과 물체 자코비안(body Jacobian)은 동일한 매니퓰레이터의 속도 기구학을 서로 다른 기준 좌표계에서 표현한 두 자코비안이다. 공간 자코비안은 기저(공간) 좌표계에서 표현된 트위스트를 관절 속도와 연결하며, 물체 자코비안은 엔드 이펙터(물체) 좌표계에서 표현된 트위스트를 관절 속도와 연결한다. 두 자코비안은 동차 변환의 수반 사상(adjoint map)에 의해 직접 연결되며, 응용 맥락에 따라 선택적으로 활용된다. 본 절에서는 두 자코비안의 학술적 정의, 제품 지수 공식으로부터의 유도, 수반 사상을 통한 상호 변환, 계산 절차, 응용 지침을 체계적으로 다룬다.
1. 공간 트위스트와 공간 자코비안의 정의
엔드 이펙터의 동차 변환 \mathbf{T}(t) \in SE(3)의 시간 미분으로부터 공간 트위스트 \mathcal{V}^s가 정의된다.
[\mathcal{V}^s(t)] = \dot{\mathbf{T}}(t) \, \mathbf{T}^{-1}(t)
여기서 [\mathcal{V}^s] \in \mathfrak{se}(3)은 트위스트의 4 \times 4 반대칭 행렬 표현이다. 이를 6차원 벡터 표현으로 풀어 쓰면
\mathcal{V}^s = \begin{bmatrix} \vec{\omega}^s \\ \vec{v}^s \end{bmatrix} \in \mathbb{R}^6
여기서 \vec{\omega}^s는 공간 좌표계에서 표현된 엔드 이펙터의 각속도이고, \vec{v}^s는 공간 좌표계에서 표현된 엔드 이펙터 원점의 선속도이다. 공간 트위스트의 기하학적 해석은 관성 공간에 고정된 관찰자의 시점에서 본 엔드 이펙터 강체의 순간 운동이다.
공간 자코비안 \mathbf{J}_s(\vec{q}) \in \mathbb{R}^{6 \times n}은 관절 속도 \dot{\vec{q}}와 공간 트위스트를 연결하는 선형 사상의 계수 행렬이다.
\mathcal{V}^s = \mathbf{J}_s(\vec{q}) \, \dot{\vec{q}}
32.19.2 물체 트위스트와 물체 자코비안의 정의
엔드 이펙터의 동차 변환 \mathbf{T}(t)로부터 물체 트위스트 \mathcal{V}^b가 다음과 같이 정의된다.
[\mathcal{V}^b(t)] = \mathbf{T}^{-1}(t) \, \dot{\mathbf{T}}(t)
물체 트위스트의 6차원 벡터 표현은
\mathcal{V}^b = \begin{bmatrix} \vec{\omega}^b \\ \vec{v}^b \end{bmatrix} \in \mathbb{R}^6
여기서 \vec{\omega}^b는 엔드 이펙터(물체) 좌표계에서 표현된 각속도이고, \vec{v}^b는 엔드 이펙터 좌표계에서 표현된 엔드 이펙터 원점의 선속도이다. 물체 트위스트의 기하학적 해석은 엔드 이펙터 자체에 부착된 관찰자의 시점에서 본 순간 운동이다.
물체 자코비안 \mathbf{J}_b(\vec{q}) \in \mathbb{R}^{6 \times n}은 관절 속도와 물체 트위스트를 연결한다.
\mathcal{V}^b = \mathbf{J}_b(\vec{q}) \, \dot{\vec{q}}
2. 제품 지수 공식으로부터의 유도
제품 지수 공식 \mathbf{T}(\vec{q}) = e^{[\mathcal{S}_1] q_1} \cdots e^{[\mathcal{S}_n] q_n} \, \mathbf{M}의 시간 미분과 \dot{\mathbf{T}} \mathbf{T}^{-1}의 계산을 통해 공간 자코비안의 각 열이 다음과 같이 유도된다.
\vec{J}_{s,i}(\vec{q}) = [\text{Ad}_{e^{[\mathcal{S}_1] q_1} \cdots e^{[\mathcal{S}_{i-1}] q_{i-1}}}] \, \mathcal{S}_i
즉 공간 자코비안의 i번째 열은 영 구성에서의 공간 관절 스크류 \mathcal{S}_i를 관절 i 이전의 모든 관절 운동에 의한 수반 변환으로 사상한 결과이다. 특히 \vec{J}_{s,1} = \mathcal{S}_1이다.
같은 방식으로 \mathbf{T}^{-1} \dot{\mathbf{T}}의 계산을 통해 물체 자코비안의 각 열이 유도된다.
\vec{J}_{b,i}(\vec{q}) = [\text{Ad}_{e^{-[\mathcal{B}_n] q_n} \cdots e^{-[\mathcal{B}_{i+1}] q_{i+1}}}] \, \mathcal{B}_i
여기서 \mathcal{B}_i는 엔드 이펙터(물체) 좌표계에서 표현된 관절 스크류이며, 공간 좌표계 관절 스크류와는
\mathcal{B}_i = [\text{Ad}_{\mathbf{M}^{-1}}] \, \mathcal{S}_i
의 관계를 가진다. 물체 자코비안의 i번째 열은 물체 좌표계 관절 스크류 \mathcal{B}_i를 관절 i 이후의 모든 관절 운동의 역변환으로 사상한 결과이다. 특히 \vec{J}_{b,n} = \mathcal{B}_n이다.
32.19.4 수반 사상을 통한 상호 변환
공간 자코비안과 물체 자코비안은 엔드 이펙터의 현재 동차 변환 \mathbf{T}(\vec{q})의 수반 사상으로 직접 연결된다.
\mathbf{J}_s(\vec{q}) = [\text{Ad}_{\mathbf{T}(\vec{q})}] \, \mathbf{J}_b(\vec{q})
역변환은
\mathbf{J}_b(\vec{q}) = [\text{Ad}_{\mathbf{T}^{-1}(\vec{q})}] \, \mathbf{J}_s(\vec{q})
이 관계는 공간 트위스트와 물체 트위스트의 변환 관계 \mathcal{V}^s = [\text{Ad}_{\mathbf{T}}] \mathcal{V}^b로부터 직접 유도된다.
수반 사상의 블록 구조는 [\vec{\omega}^\top, \vec{v}^\top]^\top 규약에서 다음과 같다.
[\text{Ad}_\mathbf{T}] = \begin{bmatrix} \mathbf{R} & \mathbf{0} \\ [\vec{p}]_\times \mathbf{R} & \mathbf{R} \end{bmatrix}, \qquad \mathbf{T} = (\mathbf{R}, \vec{p})
이 블록 삼각 구조는 각속도 블록이 단순 회전 변환으로 사상되는 반면, 선속도 블록은 병진 오프셋에 기인한 각속도 결합 항을 포함함을 보인다. 즉 두 좌표계 사이에 병진 오프셋이 존재하면 각속도가 선속도 블록에 영향을 미친다.
3. 공간 자코비안과 기하학적 자코비안의 관계
기저 좌표계에서 표현된 기하학적 자코비안 \mathbf{J}_g와 공간 자코비안 \mathbf{J}_s는 본질적으로 동일한 수학적 객체의 서로 다른 성분 배치이다. 기하학적 자코비안이 선속도를 상단, 각속도를 하단에 배치하는 반면, 공간 자코비안은 각속도를 상단, 선속도를 하단에 배치한다(본 절의 규약).
또한 공간 자코비안의 선속도 성분은 공간 좌표계 원점 기준이며, 기하학적 자코비안의 선속도 성분은 엔드 이펙터 원점 기준이라는 차이가 있다. 두 표현은 위치 벡터 \vec{p}의 외적 항을 통해 연결된다. 구체적으로 공간 자코비안의 공간 선속도 \vec{v}^s와 엔드 이펙터 원점의 선속도 \vec{v}_n는 다음 관계를 가진다.
\vec{v}^s = \vec{v}_n + \vec{p}_n \times \vec{\omega}^s = \vec{v}_n - \vec{\omega}^s \times \vec{p}_n
따라서 공간 자코비안을 기하학적 자코비안으로 변환하려면 각 열에 대해 이러한 보정을 적용한다.
물체 자코비안은 엔드 이펙터 좌표계에서 표현된 기하학적 자코비안과 성분 순서만 차이를 가진다. 물체 자코비안의 선속도 블록 \vec{v}^b는 엔드 이펙터 원점의 엔드 이펙터 좌표계 표현이므로 기하학적 의미가 직접적이다.
32.19.6 힘-속도 쌍대성과 두 자코비안
공간 자코비안과 물체 자코비안은 각각 공간 좌표계와 물체 좌표계에서 표현된 렌치와 관절 토크의 쌍대 관계를 제공한다.
공간 자코비안 전치 관계:
\vec{\tau} = \mathbf{J}_s^\top(\vec{q}) \, \mathcal{F}^s
여기서 \mathcal{F}^s는 공간 좌표계에서 표현된 렌치이다.
물체 자코비안 전치 관계:
\vec{\tau} = \mathbf{J}_b^\top(\vec{q}) \, \mathcal{F}^b
여기서 \mathcal{F}^b는 물체 좌표계에서 표현된 렌치이다.
두 렌치는 수반 사상의 전치에 의해 연결된다. \mathcal{F}^s = [\text{Ad}_{\mathbf{T}^{-1}}]^\top \mathcal{F}^b. 이 관계는 동일한 물리적 렌치가 좌표계에 따라 서로 다른 성분 표현을 가짐을 반영한다.
32.19.7 계산 절차의 비교
두 자코비안의 수치적 계산은 다음과 같이 수행된다.
32.19.7.1 공간 자코비안의 계산
- 영 구성 관절 스크류 \mathcal{S}_1, \ldots, \mathcal{S}_n을 매니퓰레이터 설계로부터 사전 결정한다.
- 누적 수반 변환을 전방으로 재귀 계산한다: \mathbf{T}_0 = \mathbf{I}, \mathbf{T}_i = \mathbf{T}_{i-1} \cdot e^{[\mathcal{S}_i] q_i}.
- 각 열을 계산한다: \vec{J}_{s,1} = \mathcal{S}_1, \vec{J}_{s,i} = [\text{Ad}_{\mathbf{T}_{i-1}}] \mathcal{S}_i (i \geq 2).
32.19.7.2 물체 자코비안의 계산
- 물체 관절 스크류 \mathcal{B}_i = [\text{Ad}_{\mathbf{M}^{-1}}] \mathcal{S}_i를 사전 계산한다.
- 누적 수반 변환을 후방으로 재귀 계산한다: \mathbf{T}'_{n+1} = \mathbf{I}, \mathbf{T}'_i = e^{-[\mathcal{B}_i] q_i} \cdot \mathbf{T}'_{i+1}.
- 각 열을 계산한다: \vec{J}_{b,n} = \mathcal{B}_n, \vec{J}_{b,i} = [\text{Ad}_{\mathbf{T}'_{i+1}}] \mathcal{B}_i (i < n).
두 방법 모두 \mathcal{O}(n) 시간 복잡도를 가지며, 기하학적 자코비안의 재귀 계산과 동등한 효율을 가진다.
32.19.7.3 한 자코비안으로부터 다른 자코비안의 변환
한 자코비안이 이미 계산된 경우 다른 자코비안은 수반 변환 [\text{Ad}_{\mathbf{T}(\vec{q})}] 또는 그 역의 행렬 곱으로 얻어진다.
\mathbf{J}_s = [\text{Ad}_{\mathbf{T}}] \, \mathbf{J}_b, \qquad \mathbf{J}_b = [\text{Ad}_{\mathbf{T}^{-1}}] \, \mathbf{J}_s
6 \times 6 행렬과 6 \times n 행렬의 곱이므로 계산량은 \mathcal{O}(n)이다.
4. 응용에 따른 선택 지침
공간 자코비안과 물체 자코비안의 선택은 응용 맥락에 의해 결정된다.
공간 자코비안의 주 응용:
- 관성 공간 기준의 작업 공간 제어. 로봇 주위 고정 기준계에서 엔드 이펙터 속도를 지령하는 제어.
- 외부 센서(비전 시스템, 모션 캡처)에서 얻은 공간 좌표계 기준 목표의 추종 제어.
- 시각화와 사용자 인터페이스에서 세계 좌표계 기반의 운동 표현.
물체 자코비안의 주 응용:
- 엔드 이펙터 기준 좌표계 제어. 엔드 이펙터에 부착된 도구의 시점에서의 힘 제어.
- 접촉 작업에서 접촉면 기준의 힘 제어와 순응 제어.
- 힘 센서의 감지 좌표계가 엔드 이펙터와 일치하는 경우의 자연스러운 수식 표현.
- 임피던스 제어에서 엔드 이펙터 관점의 임피던스 설계.
일반적으로 힘 제어와 접촉 작업에서는 물체 자코비안이 유리하고, 자유 공간 운동과 관성 기준 추종에서는 공간 자코비안이 유리하다. 두 자코비안은 수반 사상으로 직접 변환 가능하므로, 응용에 맞게 전환하여 사용하는 것이 실무적이다.
5. 좌표계의 일관성과 수치적 유의점
공간 자코비안과 물체 자코비안을 혼용하는 경우 좌표계의 일관성이 엄격히 유지되어야 한다. 다음의 조합은 올바른 물리적 관계를 만족한다.
- 공간 자코비안과 공간 트위스트: \mathcal{V}^s = \mathbf{J}_s \dot{\vec{q}}.
- 물체 자코비안과 물체 트위스트: \mathcal{V}^b = \mathbf{J}_b \dot{\vec{q}}.
- 공간 자코비안 전치와 공간 렌치: \vec{\tau} = \mathbf{J}_s^\top \mathcal{F}^s.
- 물체 자코비안 전치와 물체 렌치: \vec{\tau} = \mathbf{J}_b^\top \mathcal{F}^b.
공간 자코비안과 물체 트위스트를 혼용하거나 역으로 하면 수반 변환의 누락으로 오류가 발생한다. 수치 구현에서 이러한 일관성 검사는 단위 테스트와 시뮬레이션 검증에서 필수적이다.
수반 사상의 계산은 회전 행렬과 반대칭 행렬의 곱을 포함하므로 수치적으로 안정적이다. 그러나 회전 행렬의 직교성이 부동소수점 오차로 소실되는 경우 수반 사상의 정확성이 저하될 수 있으므로, 정기적 재직교화가 권장된다.
6. 학술적 의의
본 절에서 정립한 공간 자코비안과 물체 자코비안의 체계는 로봇 공학에서 다음의 학술적 의의를 가진다. 첫째, 공간 트위스트와 물체 트위스트의 명확한 구분은 좌표계 선택의 수학적 언어를 제공한다. 둘째, 제품 지수 공식과 결합된 두 자코비안의 닫힌 형식 유도는 데나빗-하르텐버그 매개변수 없는 현대적 기구학 해석의 표준 도구이다. 셋째, 수반 사상 [\text{Ad}_\mathbf{T}]을 매개로 한 두 자코비안의 직접 변환 공식은 좌표계 전환의 체계적 수단을 제공한다. 넷째, 힘-속도 쌍대성의 두 자코비안 형태는 공간 좌표계 기반 제어와 물체 좌표계 기반 제어 각각에 일관된 수학 틀을 제공한다. 다섯째, 두 자코비안의 응용별 선택 지침은 작업 공간 제어, 힘 제어, 임피던스 제어 등의 다양한 제어 설계에서 실무적 판단의 기준을 제시한다. 여섯째, 본 절의 내용은 현대 로봇 공학 소프트웨어(Pinocchio, Drake, OpenRAVE 등)의 표준 구현 기반이며, 리 군 기반 기구학 이론의 핵심적 귀결이다.
7. 출처
- Brockett, R. W., “Robotic manipulators and the product of exponentials formula”, Mathematical Theory of Networks and Systems, Springer, pp. 120–129, 1984.
- Murray, R. M., Li, Z., and Sastry, S. S., A Mathematical Introduction to Robotic Manipulation, CRC Press, 1994.
- Lynch, K. M. and Park, F. C., Modern Robotics: Mechanics, Planning, and Control, Cambridge University Press, 2017.
- Park, F. C., “Computational aspects of the product-of-exponentials formula for robot kinematics”, IEEE Transactions on Automatic Control, Vol. 39, No. 3, pp. 643–647, 1994.
- Park, F. C., Bobrow, J. E., and Ploen, S. R., “A Lie group formulation of robot dynamics”, International Journal of Robotics Research, Vol. 14, No. 6, pp. 609–618, 1995.
- Selig, J. M., Geometric Fundamentals of Robotics, 2nd edition, Springer, 2005.
- Featherstone, R., Rigid Body Dynamics Algorithms, Springer, 2008.
- Spong, M. W., Hutchinson, S., and Vidyasagar, M., Robot Modeling and Control, 2nd edition, Wiley, 2020.
8. 버전
- 문서 버전: 1.1
- 작성일: 2026-04-19