32.5 자코비안 행렬의 정의와 기하학적 의미

자코비안 행렬(Jacobian matrix)은 속도 기구학의 중심 수학적 객체이다. 이는 매니퓰레이터의 관절 공간 속도와 작업 공간 속도 사이의 선형 사상(linear mapping)을 기술하며, 동시에 순기구학 사상의 1계 테일러 계수로서 국소 선형 근사의 계수 행렬 역할을 한다. 자코비안은 대수적 객체로서의 행렬 구조와 기하학적 객체로서의 미분 사상 구조를 동시에 가지며, 로봇 공학의 해석과 제어 전반에 관통하는 핵심 도구이다. 본 절에서는 자코비안 행렬의 수학적 정의, 속도 사상으로서의 역할, 열 벡터 해석, 기하학적·해석적 자코비안의 구분, 접공간 관점에서의 의미, 그리고 부분 공간 구조, 계수, 특이값 분해와 조건수 등 자코비안의 기하학적·대수적 의미를 체계적으로 다룬다.

1. 자코비안 행렬의 수학적 정의

관절 변수 벡터 $\vec{q} \in \mathbb{R}^n$ 를 매개변수로 하는 순기구학 사상 $\vec{x} = \vec{f}(\vec{q}) \in \mathbb{R}^m$ 에 대하여, 자코비안 행렬은 $\vec{f}$ 의 $\vec{q}$ 에 대한 1계 편미분 행렬로 정의된다.

$\mathbf{J}(\vec{q}) = \frac{\partial \vec{f}}{\partial \vec{q}} = \begin{bmatrix} \dfrac{\partial f_1}{\partial q_1} & \cdots & \dfrac{\partial f_1}{\partial q_n} \\ \vdots & \ddots & \vdots \\ \dfrac{\partial f_m}{\partial q_1} & \cdots & \dfrac{\partial f_m}{\partial q_n} \end{bmatrix} \in \mathbb{R}^{m \times n}$

행 수 $m$ 은 작업 공간 속도의 차원이며, 열 수 $n$ 은 관절 공간의 차원이다. 공간 운동의 일반적 경우 $m = 6$ 이고, 평면 운동의 경우 $m = 3$ 이며, 위치만 고려하는 응용에서는 $m = 3$ 또는 $m = 2$ 이다.

자코비안은 일반적으로 관절 구성 $\vec{q}$ 의 비선형 함수이다. 즉 매니퓰레이터의 자세가 변하면 자코비안의 성분 값들이 변하며, 이러한 비선형 의존성이 자코비안 기반 기구학 해석의 많은 성질을 결정한다.

32.5.2 속도 사상으로서의 자코비안

자코비안은 속도 수준에서 관절 공간과 작업 공간 사이의 선형 사상을 정의한다. 관절 속도 벡터 $\dot{\vec{q}}$ 와 작업 공간 속도 벡터 $\dot{\vec{x}}$ 사이의 관계는 다음과 같다.

$\dot{\vec{x}} = \mathbf{J}(\vec{q}) \, \dot{\vec{q}}$

이 사상은 고정된 $\vec{q}$ 에서 선형이며, 관절 속도에 대해 중첩 원리가 성립한다. 그러나 $\mathbf{J}(\vec{q})$ 자체는 $\vec{q}$ 에 비선형적으로 의존하므로 매니퓰레이터가 움직이면 사상의 계수 행렬이 변화한다. 이러한 이중적 성질은 자코비안을 “국소적 선형화”(local linearization)의 관점에서 이해하게 만든다.

역방향 문제, 즉 주어진 작업 공간 속도 $\dot{\vec{x}}$ 를 실현하는 관절 속도 $\dot{\vec{q}}$ 를 구하는 문제는 자코비안의 역행렬(정방 가역의 경우) 또는 의사 역행렬(비정방·특이의 경우)을 매개로 수행되며, 이는 역속도 기구학의 학술적 기초이다.

2. 자코비안의 열 벡터 해석

자코비안의 $i$ 번째 열 벡터 $\vec{J}_i(\vec{q}) \in \mathbb{R}^m$ 은 관절 $i$ 의 단위 속도가 엔드 이펙터에 만드는 속도 기여를 나타낸다. 구체적으로 다른 모든 관절이 고정된 상태에서 $\dot{q}_i = 1$ , $\dot{q}_j = 0$ ( $j \neq i$ )일 때의 엔드 이펙터 속도가 $\vec{J}_i$ 이다.

$\vec{J}_i(\vec{q}) = \frac{\partial \vec{f}}{\partial q_i}$

따라서 작업 공간 속도는 자코비안 열 벡터들의 선형 결합으로 표현된다.

$\dot{\vec{x}} = \sum_{i=1}^{n} \vec{J}_i(\vec{q}) \, \dot{q}_i$

이러한 열별 해석은 기하학적 자코비안의 관절별 구성에서 특히 유용하다. 회전 관절의 열 벡터는 관절 축 방향 단위 벡터와 엔드 이펙터까지의 변위 벡터의 외적으로 구성되는 선속도 성분과 관절 축 방향 각속도 성분으로, 직동 관절의 열 벡터는 관절 축 방향 선속도 성분과 영의 각속도 성분으로 결정된다.

3. 기하학적 자코비안과 해석적 자코비안

자코비안은 작업 공간 속도의 표현 방식에 따라 두 가지 형태로 구분된다.

기하학적 자코비안(geometric Jacobian) $\mathbf{J}_g(\vec{q}) \in \mathbb{R}^{6 \times n}$ 은 엔드 이펙터의 물리적 선속도 $\vec{v} \in \mathbb{R}^3$ 와 각속도 $\vec{\omega} \in \mathbb{R}^3$ 의 결합 벡터를 관절 속도와 연결한다.

$\begin{bmatrix} \vec{v} \\ \vec{\omega} \end{bmatrix} = \mathbf{J}_g(\vec{q}) \, \dot{\vec{q}}$

기하학적 자코비안은 선속도 부분 $\mathbf{J}_v \in \mathbb{R}^{3 \times n}$ 과 각속도 부분 $\mathbf{J}_\omega \in \mathbb{R}^{3 \times n}$ 의 수직 결합으로 분해된다. 이 표현은 회전의 매개변수 선택에 독립적이라는 장점을 가진다.

해석적 자코비안(analytical Jacobian) $\mathbf{J}_a(\vec{q}) \in \mathbb{R}^{m \times n}$ 은 위치 $\vec{p}$ 와 방향 매개변수 $\vec{\phi}$ (예: 오일러 각, RPY 각)의 시간 미분을 관절 속도와 연결한다.

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

해석적 자코비안은 순기구학 사상 $\vec{f}$ 의 편미분을 직접 구한 것이며, 방향 매개변수의 선택에 의존한다. 두 자코비안은 방향 매개변수의 시간 미분 $\dot{\vec{\phi}}$ 와 각속도 $\vec{\omega}$ 사이의 변환 행렬 $\mathbf{T}(\vec{\phi})$ 에 의해 연결된다.

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

따라서 기하학적 자코비안과 해석적 자코비안의 관계는 다음과 같이 표현된다.

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

변환 행렬 $\mathbf{T}(\vec{\phi})$ 는 오일러 각 또는 RPY 각 표현의 특이 자세에서 특이해지며, 이는 표현의 특이점(representation singularity)으로 나타난다.

4. 자코비안의 기하학적 의미: 접공간 사상

미분 기하학적 관점에서 자코비안은 관절 공간 다양체 $\mathcal{Q}$ 로부터 작업 공간 다양체 $\mathcal{X}$ 로의 매끄러운 사상 $\vec{f} : \mathcal{Q} \to \mathcal{X}$ 의 접사상(tangent map, differential)이다.

$\mathbf{J}(\vec{q}) = d\vec{f}_{\vec{q}} : T_{\vec{q}} \mathcal{Q} \to T_{\vec{f}(\vec{q})} \mathcal{X}$

즉 자코비안은 각 점 $\vec{q}$ 에서 정의된 관절 공간 접공간 $T_{\vec{q}} \mathcal{Q}$ 를 작업 공간 접공간 $T_{\vec{x}} \mathcal{X}$ 로 보내는 선형 사상이며, 관절 속도 벡터를 작업 공간 속도 벡터로 변환한다.

자코비안은 또한 순기구학 사상의 국소 선형 근사로 해석된다. 임의의 구성 $\vec{q}_0$ 근방에서 테일러 전개를 취하면

$\vec{f}(\vec{q}_0 + \delta \vec{q}) = \vec{f}(\vec{q}_0) + \mathbf{J}(\vec{q}_0) \, \delta \vec{q} + \mathcal{O}(\|\delta \vec{q}\|^2)$

가 성립하여, 자코비안이 $\vec{q}_0$ 에서 비선형 사상 $\vec{f}$ 를 가장 잘 근사하는 선형 사상의 계수 행렬임이 나타난다. 이러한 기하학적·국소 선형 해석은 미분 기구학(differential kinematics)의 학술적 토대를 이룬다.

5. 자코비안의 부분 공간 구조

자코비안은 선형 대수의 기본 정리에 따라 네 개의 부분 공간을 정의한다. 이 중 로봇 공학적으로 특히 중요한 두 부분 공간은 다음과 같다.

상공간(range space) 또는 열 공간은 자코비안이 만드는 작업 공간 속도의 집합이다.

$\mathcal{R}(\mathbf{J}) = \{\mathbf{J}(\vec{q}) \, \dot{\vec{q}} \mid \dot{\vec{q}} \in \mathbb{R}^n\} \subseteq \mathbb{R}^m$

상공간은 현재 구성에서 실현 가능한 엔드 이펙터 속도 방향의 집합을 나타낸다. 상공간의 차원은 자코비안의 계수와 같으며, 상공간이 $\mathbb{R}^m$ 전체를 채우지 못하면 엔드 이펙터가 특정 방향으로 운동할 수 없다.

영 공간(null space) 또는 커널은 엔드 이펙터 속도에 기여하지 않는 관절 속도의 집합이다.

$\mathcal{N}(\mathbf{J}) = \{\dot{\vec{q}} \in \mathbb{R}^n \mid \mathbf{J}(\vec{q}) \, \dot{\vec{q}} = \vec{0}\} \subseteq \mathbb{R}^n$

여유 자유도 매니퓰레이터의 경우 $n > m$ 이므로 영 공간의 차원이 $n - m$ 이상이 되어 비자명하며, 이 공간의 관절 속도는 엔드 이펙터를 정지시킨 채 내부 자기 운동(self-motion)을 수행한다. 이는 부차적 과제(secondary task)의 수행을 가능하게 한다.

상공간과 영 공간은 계수-영차원 정리(rank-nullity theorem)에 의해 다음 관계를 만족한다.

$\dim \mathcal{R}(\mathbf{J}) + \dim \mathcal{N}(\mathbf{J}) = n$

32.5.7 자코비안의 계수와 특이점

자코비안의 계수(rank) $r = \text{rank}(\mathbf{J})$ 는 선형 독립인 열(또는 행)의 최대 개수이다. 계수는 항상 $r \leq \min(m, n)$ 을 만족한다.

완전 계수(full rank)의 경우 $r = \min(m, n)$ 이며, 자코비안은 주어진 차원 내에서 최대의 속도 매핑 능력을 가진다.
계수 감소(rank deficiency)의 경우 $r < \min(m, n)$ 이며, 이러한 구성을 기구학적 특이점(kinematic singularity)이라 한다.

특이점에서는 상공간의 차원이 감소하여 엔드 이펙터가 특정 방향의 순간 속도를 생성할 수 없으며, 동시에 영 공간의 차원이 증가하여 외견상 ‘여유로운’ 관절 속도가 엔드 이펙터 운동에 기여하지 않게 된다. 특이점 근방에서는 자코비안의 역행렬이 수치적으로 불안정해지므로, 역속도 기구학의 수치적 해법에서 감쇠 최소 제곱법 등의 정칙화 기법이 필요해진다.

32.5.8 자코비안의 특이값 분해

자코비안의 기하학적 성질은 특이값 분해(singular value decomposition, SVD)를 통해 정량화된다. $\mathbf{J} \in \mathbb{R}^{m \times n}$ 에 대하여 SVD는 다음과 같이 주어진다.

$\mathbf{J}(\vec{q}) = \mathbf{U}(\vec{q}) \, \mathbf{\Sigma}(\vec{q}) \, \mathbf{V}^\top(\vec{q})$

여기서 $\mathbf{U} \in \mathbb{R}^{m \times m}$ 과 $\mathbf{V} \in \mathbb{R}^{n \times n}$ 은 직교 행렬이며, $\mathbf{\Sigma} \in \mathbb{R}^{m \times n}$ 은 대각 성분에 특이값 $\sigma_1 \geq \sigma_2 \geq \cdots \geq \sigma_r > 0$ 을 갖는 대각 행렬이다. 특이값 $\sigma_i$ 는 다음의 기하학적 의미를 가진다.

좌특이 벡터 $\vec{u}_i$ 는 작업 공간에서 단위 속도를 실현하는 방향이다.
우특이 벡터 $\vec{v}_i$ 는 관절 공간에서 대응하는 단위 속도 방향이다.
특이값 $\sigma_i$ 는 관절 공간 단위 속도가 작업 공간에서 만드는 속도의 크기의 이득(gain)이다.

특이값 중 하나라도 $\sigma_i \to 0$ 이면 자코비안의 계수가 감소하며, 해당 좌특이 벡터 방향이 엔드 이펙터가 순간적으로 도달할 수 없는 방향이 된다. 이는 특이점의 정량적 정의이자 근방 해석의 기초이다.

6. 자코비안의 조건수와 등방성

자코비안의 수치적·기구학적 성능 척도로 조건수(condition number)가 사용된다.

$\kappa(\mathbf{J}) = \frac{\sigma_{\max}}{\sigma_{\min}}$

조건수는 관절 공간 단위 속도가 만들 수 있는 작업 공간 속도의 최대 크기와 최소 크기의 비이다. 조건수의 학술적 의미는 다음과 같다.

$\kappa(\mathbf{J}) = 1$ 인 경우 모든 방향에서 동등한 속도 이득이 얻어지며, 이를 등방 구성(isotropic configuration)이라 한다.
$\kappa(\mathbf{J}) \to \infty$ 인 경우 특이점에 접근하며, 수치적으로도 역사상이 불안정해진다.
조건수의 역수 $1/\kappa(\mathbf{J})$ 는 특이점으로부터의 근접 정도를 측정하는 지표이다.

조건수 외에도 행렬식 기반의 매니퓰러빌리티 측도 $w(\vec{q}) = \sqrt{\det(\mathbf{J} \mathbf{J}^\top)}$ , 최소 특이값 $\sigma_{\min}$ 등 다양한 척도가 매니퓰레이터의 기구학적 성능 평가에 활용된다.

32.5.10 힘-속도 쌍대성

자코비안은 속도 사상뿐만 아니라 힘-토크 사상을 통해 정역학과도 연결된다. 가상 일(virtual work)의 원리에 따라, 엔드 이펙터에 작용하는 외력 렌치 $\vec{F} \in \mathbb{R}^m$ 와 이를 평형시키는 관절 토크 $\vec{\tau} \in \mathbb{R}^n$ 사이에는 다음 쌍대 관계가 성립한다.

$\vec{\tau} = \mathbf{J}^\top(\vec{q}) \, \vec{F}$

이는 속도 사상 $\mathbf{J}$ 와 힘 사상 $\mathbf{J}^\top$ 가 서로 수반(adjoint) 관계에 있음을 나타낸다. 자코비안의 영 공간은 속도 측면에서 자기 운동을 허용하며, 자코비안 전치의 영 공간은 힘 측면에서 관절 토크에 나타나지 않는 외력의 존재 가능성을 설명한다. 이러한 쌍대 구조는 힘 제어, 임피던스 제어, 정역학 해석의 기초를 이룬다.

7. 학술적 의의

본 절에서 다룬 자코비안 행렬의 정의와 기하학적 의미는 로봇 공학의 여러 학술적 주제를 연결하는 근간이다. 첫째, 자코비안은 속도 기구학의 선형 사상 계수로서 관절 공간과 작업 공간의 미분적 관계를 명시한다. 둘째, 열 벡터 해석은 각 관절의 속도 기여에 대한 직관적·기하학적 이해를 제공한다. 셋째, 기하학적 자코비안과 해석적 자코비안의 구분은 방향 표현 선택에 대한 명료한 체계를 제공한다. 넷째, 접공간 사상 관점은 미분 기구학의 현대적 기하학적 언어를 제공한다. 다섯째, 상공간·영 공간·계수·특이값 분해·조건수 등 선형 대수적 성질은 특이점 분석, 매니퓰러빌리티 평가, 여유 자유도 활용의 정량적 기초를 제공한다. 여섯째, 자코비안 전치를 통한 힘-속도 쌍대성은 정역학과 속도 기구학의 통합적 해석을 가능하게 한다.

8. 출처

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.
Yoshikawa, T., “Manipulability of robotic mechanisms”, International Journal of Robotics Research, Vol. 4, No. 2, pp. 3–9, 1985.
Yoshikawa, T., Foundations of Robotics: Analysis and Control, MIT Press, 1990.
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.
Siciliano, B., Sciavicco, L., Villani, L., and Oriolo, G., Robotics: Modelling, Planning and Control, Springer, 2009.
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.
Salisbury, J. K. and Craig, J. J., “Articulated hands: Force control and kinematic issues”, International Journal of Robotics Research, Vol. 1, No. 1, pp. 4–17, 1982.

9. 버전

문서 버전: 1.1
작성일: 2026-04-19