32.22 자코비안의 상공간(Range Space) 해석

자코비안의 상공간(range space) 또는 열 공간은 관절 속도의 모든 가능한 값에 대해 자코비안이 생성하는 엔드 이펙터 속도 벡터의 집합이다. 상공간의 기하학적 해석은 매니퓰레이터가 현재 구성에서 순간적으로 실현 가능한 엔드 이펙터 속도의 방향과 크기 범위를 정량화하며, 특이점의 식별, 매니퓰러빌리티 평가, 도달 가능성 해석, 역속도 기구학의 해의 존재 조건 판정에 핵심적 역할을 한다. 본 절에서는 상공간의 학술적 정의, 기하학적 의미, 직교 여공간과의 관계, 특이점에서의 거동, 상공간 투영 연산자, 그리고 매니퓰러빌리티와의 관계를 체계적으로 다룬다.

1. 상공간의 학술적 정의

자코비안 $\mathbf{J}(\vec{q}) \in \mathbb{R}^{m \times n}$ 의 상공간은 관절 속도 벡터 $\dot{\vec{q}} \in \mathbb{R}^n$ 의 모든 값에 대한 엔드 이펙터 속도 $\dot{\vec{x}} = \mathbf{J}(\vec{q}) \dot{\vec{q}}$ 의 집합으로 정의된다.

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

상공간은 $\mathbb{R}^m$ 의 선형 부분 공간이며, 자코비안의 열 벡터들이 생성(span)하는 공간과 일치한다.

$\mathcal{R}(\mathbf{J}) = \text{span}\{\vec{J}_1, \vec{J}_2, \ldots, \vec{J}_n\}$

상공간의 차원은 자코비안의 계수와 같다.

$\dim \mathcal{R}(\mathbf{J}(\vec{q})) = \text{rank}(\mathbf{J}(\vec{q}))$

32.22.2 상공간의 기하학적 의미

상공간의 기하학적 해석은 다음과 같다. 현재 관절 구성 $\vec{q}$ 에서 엔드 이펙터가 순간적으로 실현할 수 있는 속도 벡터는 정확히 상공간 $\mathcal{R}(\mathbf{J}(\vec{q}))$ 에 포함된다. 즉 상공간 외부의 속도 벡터는 어떠한 관절 속도 선택으로도 실현할 수 없다.

공간 운동의 일반적 경우 $m = 6$ 이며 상공간은 $\mathbb{R}^6$ 의 부분 공간이다. 완전 계수 구성에서 상공간은 $\mathbb{R}^6$ 전체를 채우며, 엔드 이펙터는 모든 6차원 속도 방향을 실현할 수 있다. 특이점이나 자유도 부족 구성에서 상공간의 차원이 6보다 작으면 엔드 이펙터가 실현할 수 없는 속도 방향이 존재한다.

상공간의 기저 벡터는 엔드 이펙터 속도 공간에서의 기하학적 축을 제공하며, 이들 축 방향의 선형 결합만이 실현 가능한 속도이다. 이는 매니퓰레이터의 순간 운동학적 능력의 방향별 분포를 정량화하는 기초가 된다.

32.22.3 상공간의 기저 구성

상공간의 기저는 다음의 방법으로 구성된다.

32.22.3.1 SVD 기반 구성

자코비안의 특이값 분해 $\mathbf{J} = \mathbf{U} \mathbf{\Sigma} \mathbf{V}^\top$ 에서 영이 아닌 특이값에 대응하는 좌특이 벡터들이 상공간의 정규 직교 기저를 제공한다. $\mathbf{J}$ 의 계수가 $r$ 이면 $\mathbf{U}$ 의 처음 $r$ 개 열 $\vec{u}_1, \ldots, \vec{u}_r$ 이 $\mathcal{R}(\mathbf{J})$ 의 기저이다.

$\mathcal{R}(\mathbf{J}) = \text{span}\{\vec{u}_1, \ldots, \vec{u}_r\}$

이 기저는 정규 직교 기저이며, 각 $\vec{u}_i$ 는 관절 공간 단위 속도 $\vec{v}_i$ (우특이 벡터)가 만드는 엔드 이펙터 속도의 방향이다. 대응하는 특이값 $\sigma_i$ 가 그 방향의 속도 이득을 나타낸다.

1.1 자코비안 열 벡터의 선형 독립 부분 집합

자코비안의 열 벡터 $\vec{J}_1, \ldots, \vec{J}_n$ 중에서 선형 독립한 최대 부분 집합을 선택하여 상공간의 기저로 사용할 수도 있다. 이 방법은 관절별 물리적 의미와 연결되어 해석에 유리하나, 정규 직교성이 보장되지 않으므로 수치 해석에는 덜 안정적이다.

1.2 QR 분해 기반 구성

자코비안의 QR 분해 $\mathbf{J} = \mathbf{Q} \mathbf{R}$ 에서 비영 피벗에 대응하는 $\mathbf{Q}$ 의 열들이 상공간의 정규 직교 기저를 제공한다. QR 분해는 SVD보다 계산량이 적어 실시간 응용에 유리하다.

2. 직교 여공간: 실현 불가능한 속도 방향

상공간 $\mathcal{R}(\mathbf{J})$ 의 직교 여공간(orthogonal complement)은 상공간에 직교하는 모든 벡터의 집합이다.

$\mathcal{R}(\mathbf{J})^\perp = \{\vec{w} \in \mathbb{R}^m : \vec{w}^\top \vec{u} = 0, \; \forall \vec{u} \in \mathcal{R}(\mathbf{J})\}$

선형 대수의 기본 정리에 따라 $\mathcal{R}(\mathbf{J})^\perp$ 는 자코비안 전치의 영 공간과 일치한다.

$\mathcal{R}(\mathbf{J})^\perp = \mathcal{N}(\mathbf{J}^\top)$

$\mathcal{R}(\mathbf{J})^\perp$ 의 원소는 매니퓰레이터가 현재 구성에서 순간적으로 실현할 수 없는 엔드 이펙터 속도 방향이다. 정규 구성에서 이 공간의 차원은 $m - \text{rank}(\mathbf{J}) = 0$ (정방 완전 계수) 또는 $m - \text{rank}(\mathbf{J}) = m - n$ (자유도 부족의 경우)이다.

특이점에서는 $\mathcal{R}(\mathbf{J})^\perp$ 의 차원이 증가하며, 이는 엔드 이펙터가 실현할 수 없는 속도 방향이 늘어남을 의미한다. 이러한 실현 불가능한 방향은 특이점 해석의 핵심 대상이 된다.

3. 상공간과 영 공간의 직교 분해

자코비안 $\mathbf{J}$ 와 전치 $\mathbf{J}^\top$ 에 대응하는 네 개의 기본 부분 공간 $\mathcal{R}(\mathbf{J}), \mathcal{N}(\mathbf{J}), \mathcal{R}(\mathbf{J}^\top), \mathcal{N}(\mathbf{J}^\top)$ 는 선형 대수의 기본 정리에 의해 다음 직교 분해 관계를 만족한다.

$\mathbb{R}^m = \mathcal{R}(\mathbf{J}) \oplus \mathcal{N}(\mathbf{J}^\top)$

$\mathbb{R}^n = \mathcal{R}(\mathbf{J}^\top) \oplus \mathcal{N}(\mathbf{J})$

첫 번째 분해는 엔드 이펙터 속도 공간 $\mathbb{R}^m$ 이 실현 가능한 속도의 공간 $\mathcal{R}(\mathbf{J})$ 와 실현 불가능한 속도의 공간 $\mathcal{N}(\mathbf{J}^\top)$ 의 직교 합으로 분해됨을 나타낸다.

두 번째 분해는 관절 속도 공간 $\mathbb{R}^n$ 이 엔드 이펙터에 영향을 주는 방향 $\mathcal{R}(\mathbf{J}^\top)$ 과 자기 운동 방향 $\mathcal{N}(\mathbf{J})$ 의 직교 합으로 분해됨을 나타낸다. 이 관점에서 $\mathcal{R}(\mathbf{J}^\top)$ 의 관절 속도만이 엔드 이펙터 속도를 생성하며, $\mathcal{N}(\mathbf{J})$ 의 관절 속도는 자기 운동만을 만든다.

4. 상공간 투영 연산자

상공간으로의 직교 투영 연산자는 다음과 같이 구성된다.

$\mathbf{P}_R(\vec{q}) = \mathbf{J}(\vec{q}) \, \mathbf{J}^+(\vec{q})$

여기서 $\mathbf{J}^+$ 는 무어-펜로즈 의사 역행렬이다. 여유 자유도 매니퓰레이터의 경우 $\mathbf{J}^+ = \mathbf{J}^\top (\mathbf{J} \mathbf{J}^\top)^{-1}$ 이므로 $\mathbf{P}_R = \mathbf{J} \mathbf{J}^\top (\mathbf{J} \mathbf{J}^\top)^{-1} = \mathbf{I}_m$ 이 되며, 이는 상공간이 $\mathbb{R}^m$ 전체임을 반영한다(완전 행 계수 가정).

자유도 부족 매니퓰레이터의 경우 $\mathbf{J}^+ = (\mathbf{J}^\top \mathbf{J})^{-1} \mathbf{J}^\top$ 이고 $\mathbf{P}_R = \mathbf{J} (\mathbf{J}^\top \mathbf{J})^{-1} \mathbf{J}^\top$ 은 $\mathcal{R}(\mathbf{J}) \subsetneq \mathbb{R}^m$ 으로의 직교 투영이 된다. 이 투영 연산자는 $\mathbf{P}_R^2 = \mathbf{P}_R$ (멱등성)과 $\mathbf{P}_R^\top = \mathbf{P}_R$ (대칭성)을 만족한다.

임의의 목표 속도 벡터 $\dot{\vec{x}}_d$ 에 대해 $\mathbf{P}_R \dot{\vec{x}}_d$ 는 상공간으로의 투영이며, 이 투영된 속도가 매니퓰레이터가 실현할 수 있는 가장 가까운 속도이다.

32.22.7 역속도 기구학 해의 존재성

자코비안의 상공간은 역속도 기구학 해의 존재성에 대한 판정 기준을 제공한다.

해의 존재 조건: 목표 엔드 이펙터 속도 $\dot{\vec{x}}_d$ 에 대한 역속도 기구학 해가 존재하려면 $\dot{\vec{x}}_d \in \mathcal{R}(\mathbf{J})$ 이어야 한다.

완전 행 계수 구성( $\mathcal{R}(\mathbf{J}) = \mathbb{R}^m$ )에서는 모든 목표 속도에 대해 해가 존재한다.
자유도 부족 또는 특이 구성에서 $\dot{\vec{x}}_d \notin \mathcal{R}(\mathbf{J})$ 인 경우 엄밀한 해가 존재하지 않으며, 의사 역행렬은 최소 제곱 의미의 근사 해 $\dot{\vec{q}} = \mathbf{J}^+ \dot{\vec{x}}_d$ 를 제공한다. 이때 실제 실현되는 엔드 이펙터 속도는 $\mathbf{J} \mathbf{J}^+ \dot{\vec{x}}_d = \mathbf{P}_R \dot{\vec{x}}_d$ 이며, 이는 $\dot{\vec{x}}_d$ 의 상공간 투영이다.

목표 속도와 실현 속도의 차이 $\dot{\vec{x}}_d - \mathbf{P}_R \dot{\vec{x}}_d = (\mathbf{I} - \mathbf{P}_R) \dot{\vec{x}}_d$ 는 실현 불가능한 성분이며, 이는 $\mathcal{N}(\mathbf{J}^\top)$ 의 원소이다.

32.22.8 특이점에서의 상공간 수축

매니퓰레이터가 특이점에 접근함에 따라 상공간의 구조는 다음의 변화를 겪는다.

차원 감소: 특이점 진입 시 자코비안의 계수가 감소하고 상공간의 차원이 작아진다. 예를 들어 6자유도 정방 매니퓰레이터가 1계 특이점에 도달하면 상공간 차원이 $6$ 에서 $5$ 로 감소한다.

방향의 소실: 감소하는 차원에 대응하는 엔드 이펙터 속도 방향이 실현 불가능한 영역 $\mathcal{N}(\mathbf{J}^\top)$ 으로 이동한다. 이 방향은 특이점의 ‘닫힌 방향’(lost direction)으로 해석된다.

기하학적 해석: 매니퓰러빌리티 타원체(manipulability ellipsoid)가 특이점 근방에서 납작해지고, 특이점에서는 한 축 이상이 영으로 수축하여 저차원 다양체로 변한다.

이러한 상공간의 수축은 매니퓰레이터의 순간 운동학적 능력의 국부적 저하를 반영하며, 특이점 해석의 핵심 현상이다.

32.22.9 상공간과 매니퓰러빌리티 타원체

자코비안의 상공간은 매니퓰러빌리티 타원체의 기하학적 형상과 직접 연결된다. 관절 속도 단위 구 $\{\dot{\vec{q}} : \|\dot{\vec{q}}\| \leq 1\}$ 의 자코비안 상 $\{\mathbf{J} \dot{\vec{q}} : \|\dot{\vec{q}}\| \leq 1\}$ 는 엔드 이펙터 속도 공간의 타원체이다.

이 타원체의 주축은 자코비안의 좌특이 벡터들이며, 주축 길이는 대응 특이값이다. 상공간이 $\mathbb{R}^m$ 전체일 때 타원체는 $m$ 차원 정규 타원체이며, 상공간이 $r < m$ 차원인 경우 타원체는 $r$ 차원 축소 타원체로서 상공간에 포함된다.

매니퓰러빌리티 지수 $w = \sqrt{\det(\mathbf{J} \mathbf{J}^\top)}$ 는 매니퓰러빌리티 타원체의 부피와 비례하며, 특이점에서 $w \to 0$ 이 되어 타원체가 저차원으로 수축한다. 이러한 상공간-매니퓰러빌리티의 기하학적 연결은 매니퓰레이터의 순간 운동학적 성능 평가의 정량적 토대를 제공한다.

32.22.10 힘 영역에서의 상공간 해석

자코비안 전치를 통한 힘 사상 $\vec{\tau} = \mathbf{J}^\top \mathcal{F}$ 에서 상공간 해석은 다음과 같이 적용된다.

자코비안 전치의 상공간 $\mathcal{R}(\mathbf{J}^\top) \subseteq \mathbb{R}^n$ : 엔드 이펙터 힘이 만들 수 있는 관절 토크의 집합.

자코비안 전치의 영 공간 $\mathcal{N}(\mathbf{J}^\top) \subseteq \mathbb{R}^m$ : 관절 토크를 만들지 않는 엔드 이펙터 힘의 집합. 이는 상공간 $\mathcal{R}(\mathbf{J})$ 의 직교 여공간이며, 구조적 구속 힘(constraint force)의 공간으로 해석된다.

힘 제어 설계에서 $\mathcal{N}(\mathbf{J}^\top)$ 의 성분은 관절 토크에 반영되지 않으므로 제어 시스템이 감지할 수 없으며, 이러한 힘은 매니퓰레이터의 구조적 강성에 의해 지탱된다. 이러한 성질은 접촉 작업에서 힘 제어의 한계를 결정한다.

32.22.11 학술적 의의

본 절에서 정립한 자코비안의 상공간 해석은 로봇 공학에서 다음의 학술적 의의를 가진다. 첫째, 상공간을 엔드 이펙터가 실현 가능한 속도의 집합으로 정립함으로써 매니퓰레이터의 순간 운동학적 능력을 정량화하는 기본 개념을 제공한다. 둘째, 영 공간과의 직교 분해를 통해 관절 속도와 엔드 이펙터 속도 공간의 구조를 통일적으로 이해하는 틀을 제공한다. 셋째, SVD·QR 기반 상공간 기저 구성과 상공간 투영 연산자의 수치적 구현은 실시간 해석의 실무적 도구이다. 넷째, 역속도 기구학 해의 존재성 판정 기준은 제어 시스템의 실행 가능성 검토의 근거를 제공한다. 다섯째, 특이점에서의 상공간 수축은 기구학적 성능 저하의 기하학적 해석을 제공한다. 여섯째, 매니퓰러빌리티 타원체와의 직접적 연결은 정량적 성능 평가의 토대이다. 일곱째, 자코비안 전치의 상공간·영 공간 해석은 힘 제어 설계의 이론적 한계를 규정한다. 여덟째, 본 절의 내용은 후속 절의 영 공간 투영을 활용한 여유 자유도 제어, 매니퓰러빌리티 타원체의 상세 분석, 특이점 분류의 이론적 기초가 된다.

출처

Strang, G., Introduction to Linear Algebra, 5th edition, Wellesley-Cambridge Press, 2016.
Golub, G. H. and Van Loan, C. F., Matrix Computations, 4th edition, Johns Hopkins University Press, 2013.
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.
Khatib, O., “A unified approach for motion and force control of robot manipulators: The operational space formulation”, IEEE Journal on Robotics and Automation, Vol. 3, No. 1, pp. 43–53, 1987.
Nakamura, Y., Advanced Robotics: Redundancy and Optimization, Addison-Wesley, 1991.
Spong, M. W., Hutchinson, S., and Vidyasagar, M., Robot Modeling and Control, 2nd edition, Wiley, 2020.
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.

버전

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