32.9 해석적 자코비안(Analytical Jacobian)의 유도

해석적 자코비안(analytical Jacobian)은 엔드 이펙터의 자세를 위치 벡터와 최소 매개변수화된 방향 표현(오일러 각, RPY 각, 축-각 표현 등)의 결합으로 기술한 뒤, 이 결합 벡터의 관절 변수에 대한 1계 편미분 행렬로 정의되는 자코비안 표현이다. 기하학적 자코비안이 물리적 각속도 벡터와 관절 속도를 직접 연결하는 반면, 해석적 자코비안은 방향 매개변수의 시간 미분과 관절 속도를 연결한다는 점에서 구조적 차이를 가진다. 본 절에서는 해석적 자코비안의 학술적 정의, 위치 부분과 자세 부분의 유도, 방향 매개변수에 따른 변환 행렬의 구성, 표현 특이점, 그리고 계산 절차를 다룬다.

1. 해석적 자코비안의 학술적 정의

엔드 이펙터의 자세를 위치 벡터 $\vec{p} \in \mathbb{R}^3$ 와 방향 매개변수 벡터 $\vec{\phi} \in \mathbb{R}^k$ (일반적으로 $k = 3$ )의 결합으로 표현한다.

$\vec{x}_a = \begin{bmatrix} \vec{p} \\ \vec{\phi} \end{bmatrix} \in \mathbb{R}^{3 + k}$

관절 변수 벡터 $\vec{q} \in \mathbb{R}^n$ 의 함수로 주어지는 순기구학 사상 $\vec{x}_a = \vec{f}_a(\vec{q})$ 에 대하여, 해석적 자코비안은 이 사상의 1계 편미분 행렬로 정의된다.

$\mathbf{J}_a(\vec{q}) = \frac{\partial \vec{f}_a}{\partial \vec{q}} = \frac{\partial \vec{x}_a}{\partial \vec{q}} \in \mathbb{R}^{(3+k) \times n}$

시간 미분을 취하면 관절 속도와 매개변수화된 작업 공간 속도 사이의 선형 사상이 얻어진다.

$\dot{\vec{x}}_a = \begin{bmatrix} \dot{\vec{p}} \\ \dot{\vec{\phi}} \end{bmatrix} = \mathbf{J}_a(\vec{q}) \, \dot{\vec{q}}$

해석적 자코비안의 핵심적 학술적 특징은 그 구성이 선택된 방향 매개변수의 규약에 명시적으로 의존한다는 점이며, 이로 인해 방향 표현의 특이점(representation singularity)이 자코비안에 직접 반영된다.

32.9.2 블록 구조: 위치 부분과 자세 부분

해석적 자코비안은 자연스럽게 위치 부분과 자세 부분의 두 블록으로 분해된다.

$\mathbf{J}_a(\vec{q}) = \begin{bmatrix} \mathbf{J}_p(\vec{q}) \\ \mathbf{J}_\phi(\vec{q}) \end{bmatrix}$

여기서

$\mathbf{J}_p(\vec{q}) = \frac{\partial \vec{p}}{\partial \vec{q}} \in \mathbb{R}^{3 \times n}, \qquad \mathbf{J}_\phi(\vec{q}) = \frac{\partial \vec{\phi}}{\partial \vec{q}} \in \mathbb{R}^{k \times n}$

이다. 위치 부분 $\mathbf{J}_p$ 는 엔드 이펙터 위치의 관절 변수에 대한 편미분이며, 자세 부분 $\mathbf{J}_\phi$ 는 방향 매개변수의 관절 변수에 대한 편미분이다.

32.9.3 위치 부분의 유도와 기하학적 자코비안과의 관계

엔드 이펙터의 위치 벡터 $\vec{p}(\vec{q})$ 는 순기구학의 동차 변환 ${}^0\mathbf{T}_n(\vec{q})$ 의 병진 부분으로 직접 주어진다. 이에 대한 시간 미분은 연쇄 법칙에 따라 다음과 같다.

$\dot{\vec{p}}(t) = \frac{\partial \vec{p}}{\partial \vec{q}} \, \dot{\vec{q}} = \mathbf{J}_p(\vec{q}) \, \dot{\vec{q}}$

엔드 이펙터의 물리적 선속도 $\vec{v}$ 는 정의상 위치 벡터의 시간 미분 $\vec{v} = \dot{\vec{p}}$ 와 같으므로, 해석적 자코비안의 위치 부분은 기하학적 자코비안의 선속도 부분과 일치한다.

$\mathbf{J}_p(\vec{q}) = \mathbf{J}_v(\vec{q})$

이 등식은 위치가 3차원 유클리드 공간에서 자연스럽게 매개변수화 가능한 가환 가산 구조를 가지기 때문에 성립하며, 매개변수 선택에 관계없이 항상 유효하다.

32.9.4 자세 부분의 유도와 표현 행렬

엔드 이펙터의 물리적 각속도 $\vec{\omega}$ 와 방향 매개변수의 시간 미분 $\dot{\vec{\phi}}$ 는 일반적으로 동일하지 않다. 둘 사이의 관계는 방향 매개변수의 선택에 의존하는 표현 행렬(representation matrix) 또는 역 변환 행렬 $\mathbf{E}(\vec{\phi}) \in \mathbb{R}^{3 \times k}$ 에 의해 다음과 같이 주어진다.

$\vec{\omega} = \mathbf{E}(\vec{\phi}) \, \dot{\vec{\phi}}$

$k = 3$ 인 표준적 3매개변수 표현에서 $\mathbf{E}(\vec{\phi})$ 가 가역일 때, 매개변수의 시간 미분은 각속도로부터 역변환을 통해 얻어진다.

$\dot{\vec{\phi}} = \mathbf{E}^{-1}(\vec{\phi}) \, \vec{\omega}$

이 관계를 기하학적 자코비안의 각속도 부분 $\mathbf{J}_\omega$ 로 표현된 각속도에 적용하면 자세 부분 해석적 자코비안이 얻어진다.

$\dot{\vec{\phi}} = \mathbf{E}^{-1}(\vec{\phi}) \, \mathbf{J}_\omega(\vec{q}) \, \dot{\vec{q}} = \mathbf{J}_\phi(\vec{q}) \, \dot{\vec{q}}$

따라서

$\mathbf{J}_\phi(\vec{q}) = \mathbf{E}^{-1}(\vec{\phi}(\vec{q})) \, \mathbf{J}_\omega(\vec{q})$

가 성립한다. 이를 위치 부분과 결합하면 해석적 자코비안의 구성은 다음 블록 대각 변환의 형태로 표현된다.

$\mathbf{J}_a(\vec{q}) = \begin{bmatrix} \mathbf{I}_3 & \mathbf{0} \\ \mathbf{0} & \mathbf{E}^{-1}(\vec{\phi}) \end{bmatrix} \mathbf{J}_g(\vec{q})$

여기서 $\mathbf{J}_g$ 는 기하학적 자코비안이다. 이 관계는 해석적 자코비안과 기하학적 자코비안의 변환 공식을 제공한다.

2. ZYZ 오일러 각 기반 해석적 자코비안

ZYZ 오일러 각 표현을 사용하는 경우 방향 매개변수는 $\vec{\phi} = (\phi, \theta, \psi)$ 이며, 해당 표현 행렬 $\mathbf{E}(\vec{\phi})$ 은 세 축 회전의 각속도 기여를 누적하여 다음과 같이 유도된다. 첫 번째 회전 $\phi$ 는 공간 $z$ 축 주위, 두 번째 회전 $\theta$ 는 중간 좌표계의 $y'$ 축 주위, 세 번째 회전 $\psi$ 는 최종 중간 좌표계의 $z''$ 축 주위이다.

$\mathbf{E}(\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}$

이 행렬의 행렬식은 $\det \mathbf{E} = \sin \theta$ 로 주어지므로, $\theta = 0$ 또는 $\theta = \pi$ 인 구성에서 $\mathbf{E}$ 가 특이해진다. 이 구성이 바로 ZYZ 오일러 각의 표현 특이점(짐벌 잠김)이다. 이러한 구성에서는 해석적 자코비안의 자세 부분 $\mathbf{J}_\phi = \mathbf{E}^{-1} \mathbf{J}_\omega$ 가 정의되지 않거나 수치적으로 극도로 불안정해진다.

32.9.6 ZYX(RPY) 각 기반 해석적 자코비안

항공 분야와 이동 로봇에서 널리 사용되는 RPY(roll–pitch–yaw) 각 또는 ZYX 오일러 각 표현을 사용하는 경우 방향 매개변수는 $\vec{\phi} = (\alpha, \beta, \gamma)$ 이며(각각 요, 피치, 롤 또는 규약에 따라 다른 순서), 표현 행렬은 다음과 같이 유도된다. 공간 $z$ 축 주위 $\alpha$ 회전, 중간 $y'$ 축 주위 $\beta$ 회전, 최종 $x''$ 축 주위 $\gamma$ 회전을 차례로 적용한 경우

$\mathbf{E}(\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}$

이 행렬의 행렬식은 $\det \mathbf{E} = \cos \beta$ 로 주어지므로, $\beta = \pm \pi/2$ 인 구성에서 특이해진다. 이는 RPY 표현의 표현 특이점이며, 해당 구성에서 요와 롤 축이 일치하여 독립적 매개변수화가 불가능해지는 ‘짐벌 잠김’ 현상에 대응한다.

다양한 오일러 각 규약과 RPY 규약에 따라 표현 행렬의 구체적 형태는 달라지며, 특이점의 발생 위치도 이에 따라 변한다. 규약을 명시하지 않으면 수치적 일관성이 유지되지 않으므로 학술적·실무적 응용에서는 규약을 명시적으로 기술하는 것이 필수적이다.

3. 표현 특이점과 기구학적 특이점의 구분

해석적 자코비안은 두 가지 서로 다른 유형의 특이성을 반영한다.

기구학적 특이점(kinematic singularity)은 기하학적 자코비안 $\mathbf{J}_g$ 의 계수 감소로 정의되는 매니퓰레이터의 본질적 기구학적 특성이다. 이는 로봇의 물리적 구성에 의한 것으로 방향 매개변수의 선택과 무관하다.

표현 특이점(representation singularity)은 방향 매개변수의 시간 미분과 각속도 사이의 변환 행렬 $\mathbf{E}(\vec{\phi})$ 의 특이성에 의한 것이다. 이는 매개변수화 규약의 선택에 따라 나타나는 수학적 특이점이며, 실제 로봇의 물리적 운동 능력에는 영향을 미치지 않는다.

해석적 자코비안은 두 유형의 특이점을 모두 반영한다.

$\det \mathbf{J}_a = \det \left( \begin{bmatrix} \mathbf{I}_3 & \mathbf{0} \\ \mathbf{0} & \mathbf{E}^{-1}(\vec{\phi}) \end{bmatrix} \right) \det \mathbf{J}_g = \frac{\det \mathbf{J}_g}{\det \mathbf{E}(\vec{\phi})}$

따라서 기하학적 자코비안이 완전 계수인 구성에서도 $\det \mathbf{E} = 0$ 이면 해석적 자코비안이 특이해진다. 반대로 기하학적 자코비안이 특이하면 해석적 자코비안도 특이해지므로, 해석적 자코비안의 특이점 해석에는 두 요인을 분리하여 고려할 필요가 있다.

32.9.8 해석적 자코비안과 기하학적 자코비안의 선택

두 자코비안의 학술적·실무적 장단점은 응용 맥락에 따라 상이하다.

해석적 자코비안의 장점은 방향 매개변수의 시간 미분을 직접 제공하므로, 방향 매개변수에 기반한 궤적 생성과 추종 제어에 자연스럽게 대응한다는 점이다. 특히 목표 자세가 오일러 각 또는 RPY 각으로 주어지는 궤적 계획에서 계산이 직접적이다.

기하학적 자코비안의 장점은 물리적 각속도를 직접 다루므로 방향 매개변수의 선택에 독립적이며, 표현 특이점에서 자유롭다는 점이다. 또한 스크류 이론과 리 군-리 대수 기반 현대적 기구학의 자연스러운 수학적 언어이다.

실무적으로 작업 공간 제어(operational-space control)와 임피던스 제어, 그리고 방향 오차를 각속도 또는 축-각 표현으로 처리하는 시스템에서는 기하학적 자코비안이 선호된다. 반면 궤적 추종 시 방향을 오일러 각으로 보간하는 단순한 구현에서는 해석적 자코비안이 직접적으로 사용된다.

32.9.9 해석적 자코비안의 계산 절차

해석적 자코비안의 구성은 다음 체계적 단계로 수행된다.

단계 1: 순기구학 계산. 관절 변수 $\vec{q}$ 에 대하여 엔드 이펙터의 동차 변환 ${}^0\mathbf{T}_n(\vec{q})$ 를 계산하고, 병진 부분에서 $\vec{p}(\vec{q})$ 를, 회전 부분 ${}^0\mathbf{R}_n(\vec{q})$ 로부터 방향 매개변수 $\vec{\phi}(\vec{q})$ 를 추출한다.

단계 2: 위치 부분 구성. $\mathbf{J}_p = \partial \vec{p} / \partial \vec{q}$ 를 기호적 미분 또는 기하학적 자코비안의 선속도 블록 $\mathbf{J}_v$ 로부터 얻는다.

단계 3: 표현 행렬 계산. 선택된 방향 매개변수 규약에 따라 표현 행렬 $\mathbf{E}(\vec{\phi})$ 를 구성한다.

단계 4: 자세 부분 구성. 기하학적 자코비안의 각속도 블록 $\mathbf{J}_\omega$ 로부터 $\mathbf{J}_\phi = \mathbf{E}^{-1}(\vec{\phi}) \, \mathbf{J}_\omega$ 를 계산한다. 또는 자세 매개변수 $\vec{\phi}(\vec{q})$ 에 대한 기호적 편미분을 직접 수행한다.

단계 5: 행렬 조립. 위치 부분과 자세 부분을 수직으로 결합하여 $\mathbf{J}_a = [\mathbf{J}_p^\top, \mathbf{J}_\phi^\top]^\top$ 을 구성한다.

기호 연산 소프트웨어를 활용하면 위 절차가 자동화될 수 있으나, 기호 표현의 크기가 증가하는 단점이 있다. 수치적 구현에서는 단계 4에서 $\mathbf{E}^{-1} \mathbf{J}_\omega$ 를 명시적 역행렬 대신 선형 연립 방정식 $\mathbf{E} \mathbf{J}_\phi = \mathbf{J}_\omega$ 의 해로 구하는 것이 수치적으로 안정적이다.

32.9.10 학술적 의의

본 절에서 정립한 해석적 자코비안의 유도 방법은 로봇 공학에서 다음의 학술적 의의를 가진다. 첫째, 방향의 최소 매개변수 표현과 관절 변수 사이의 명시적 1계 미분 관계를 제공한다. 둘째, 해석적 자코비안과 기하학적 자코비안의 변환 공식 $\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.
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.
Whitney, D. E., “The mathematics of coordinated control of prosthetic arms and manipulators”, Journal of Dynamic Systems, Measurement, and Control, Vol. 94, No. 4, pp. 303–309, 1972.
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