32.2 관절 공간 속도와 작업 공간 속도의 관계

32.2 관절 공간 속도와 작업 공간 속도의 관계

관절 공간 속도와 작업 공간 속도의 관계는 속도 기구학의 중심 주제이다. 두 공간은 로봇 매니퓰레이터의 서로 다른 운동학적 변수 집합을 기술하며, 자코비안 행렬(Jacobian matrix)을 매개로 1계 시간 미분 수준에서 선형 사상(linear mapping)을 형성한다. 본 절에서는 관절 공간과 작업 공간의 학술적 정의, 두 공간 속도의 미분적 대응 관계, 기하학적 해석, 그리고 역매핑의 수학적 구조를 다룬다.

1. 관절 공간과 작업 공간의 학술적 정의

관절 공간(joint space) \mathcal{Q}는 로봇의 모든 관절 변수 집합으로 구성되는 공간이다. n 자유도의 직렬 매니퓰레이터에 대하여 관절 변수 벡터는 다음과 같이 정의된다.

\vec{q} = [q_1, q_2, \ldots, q_n]^\top \in \mathcal{Q}

회전 관절의 경우 q_i = \theta_i이고 그 자연스러운 공간은 원 S^1이며, 따라서 모든 회전 관절로 구성된 매니퓰레이터의 관절 공간은 원환면 \mathbb{T}^n이다. 직동 관절의 경우 q_i = d_i이며 공간은 \mathbb{R}이다. 해석의 편의를 위해 관절 한계 이내의 국소적 영역에서 관절 공간을 \mathbb{R}^n의 부분 집합으로 간주하는 것이 일반적이다.

작업 공간(task space 또는 operational space) \mathcal{X}는 엔드 이펙터의 자세(pose)가 놓이는 공간이다. 공간 운동의 경우 자세는 위치 \vec{p} \in \mathbb{R}^3와 방향 \mathbf{R} \in SO(3)의 쌍으로 구성되며, 자연스러운 공간은 특수 유클리드 군 SE(3)이다. 방향을 오일러 각, RPY 각, 축-각 표현(axis-angle representation) 등으로 매개변수화하면 자세를 \mathbb{R}^6의 벡터로 표현할 수 있다.

\vec{x} = \begin{bmatrix} \vec{p} \\ \vec{\phi} \end{bmatrix} \in \mathbb{R}^6

여기서 \vec{\phi}는 방향의 최소 매개변수 표현이다. 이러한 매개변수화는 국소적으로 유효하나 전역적으로는 표현의 특이점(representation singularity)을 수반하므로 주의가 요구된다.

두 공간의 차원 관계에 따라 매니퓰레이터의 운동학적 특성이 결정된다. n = 6이면 공간 자세 전 영역에 대응하는 표준 매니퓰레이터이며, n > 6이면 운동학적 여유 자유도(kinematically redundant) 매니퓰레이터, n < 6이면 자유도 부족(under-actuated) 매니퓰레이터이다.

2. 순기구학 사상과 시간 미분

관절 공간으로부터 작업 공간으로의 매핑인 순기구학(forward kinematics)은 일반적으로 비선형 함수이다.

\vec{x} = \vec{f}(\vec{q}), \qquad \vec{f} : \mathcal{Q} \to \mathcal{X}

비선형성의 주된 원인은 회전 변환에서 등장하는 삼각함수 항과 동차 변환 행렬(homogeneous transformation matrix)의 누적 곱셈이다.

관절 변수가 시간의 연속 함수 \vec{q}(t)로 주어질 때, 순기구학 사상을 통하여 작업 공간 변수 역시 시간의 함수가 된다. 양변을 시간에 대해 미분하고 연쇄 법칙(chain rule)을 적용하면 다음의 속도 관계가 얻어진다.

\dot{\vec{x}}(t) = \frac{d \vec{x}}{dt} = \frac{\partial \vec{f}}{\partial \vec{q}} \, \frac{d \vec{q}}{dt} = \mathbf{J}(\vec{q}) \, \dot{\vec{q}}

여기서

\mathbf{J}(\vec{q}) = \frac{\partial \vec{f}}{\partial \vec{q}} \in \mathbb{R}^{m \times n}

이 자코비안 행렬이다. 이 식은 관절 공간 속도 \dot{\vec{q}}와 작업 공간 속도 \dot{\vec{x}} 사이의 핵심 관계식으로, 본 절의 중심 결과이다.

32.2.3 선형 사상으로서의 자코비안

자코비안 \mathbf{J}(\vec{q})는 고정된 관절 구성 \vec{q}에서 관절 공간의 접공간(tangent space) T_{\vec{q}}\mathcal{Q}로부터 작업 공간의 접공간 T_{\vec{x}}\mathcal{X}로의 선형 사상이다.

\mathbf{J}(\vec{q}) : T_{\vec{q}}\mathcal{Q} \to T_{\vec{x}}\mathcal{X}

구체적으로 관절 속도 \dot{\vec{q}} \in T_{\vec{q}}\mathcal{Q}가 주어지면 작업 공간 속도 \dot{\vec{x}} = \mathbf{J}(\vec{q}) \dot{\vec{q}} \in T_{\vec{x}}\mathcal{X}가 결정된다. 이 사상은 고정된 \vec{q}에서 선형이지만, 자코비안 자체는 \vec{q}의 비선형 함수이므로 매니퓰레이터의 자세가 변하면 사상의 계수도 변한다.

고정된 \vec{q}_0 근방에서 순기구학 사상의 1계 테일러 전개를 취하면 다음이 성립한다.

\vec{x} \approx \vec{f}(\vec{q}_0) + \mathbf{J}(\vec{q}_0) \, (\vec{q} - \vec{q}_0)

따라서 자코비안은 순기구학 사상의 국소 선형 근사(local linear approximation)의 계수 행렬이다. 이러한 국소 선형성은 역기구학의 뉴턴-랩슨(Newton–Raphson) 반복법, 경로 계획의 국소 최적화, 제어 설계의 선형 해석 등 다양한 응용의 기초가 된다.

32.2.4 자코비안의 차원과 구조

자코비안 \mathbf{J}(\vec{q}) \in \mathbb{R}^{m \times n}의 행 수 m은 작업 공간 속도의 차원이며, 열 수 n은 관절 공간 속도의 차원이다.

  • 작업 공간이 3차원 위치만을 포함하는 경우 m = 3이다.
  • 공간 운동의 일반적 경우 m = 6이며, 선속도 3차원과 각속도 3차원을 포함한다.
  • 평면 운동의 경우 m = 3이며, 평면 내 선속도 2차원과 평면 법선 방향 각속도 1차원을 포함한다.

자코비안의 각 열 \vec{J}_i(\vec{q}) \in \mathbb{R}^m은 i번째 관절 속도 \dot{q}_i가 작업 공간 속도에 기여하는 성분을 나타낸다.

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

이러한 열별 해석은 기하학적 자코비안(geometric Jacobian)의 열 구성에서 핵심적 역할을 수행한다.

3. 역속도 사상

작업 공간 속도 \dot{\vec{x}}_d가 주어질 때 이를 실현하는 관절 속도 \dot{\vec{q}}를 구하는 문제가 역속도 기구학(inverse velocity kinematics)이다. 자코비안의 구조에 따라 해의 성격이 달라진다.

자코비안이 정방이고 가역인 경우 (m = n이고 \det \mathbf{J} \neq 0): 유일한 해가 존재한다.

\dot{\vec{q}} = \mathbf{J}^{-1}(\vec{q}) \, \dot{\vec{x}}_d

여유 자유도의 경우 (n > m이고 \mathbf{J}의 행 계수가 m): 해가 무한히 존재하며, 최소 노름(minimum norm) 해는 무어-펜로즈 의사 역행렬 \mathbf{J}^+ = \mathbf{J}^\top(\mathbf{J}\mathbf{J}^\top)^{-1}를 이용하여 얻어진다. 일반 해는 다음과 같이 영 공간 성분을 포함한다.

\dot{\vec{q}} = \mathbf{J}^+(\vec{q}) \, \dot{\vec{x}}_d + (\mathbf{I}_n - \mathbf{J}^+(\vec{q}) \, \mathbf{J}(\vec{q})) \, \dot{\vec{q}}_0

여기서 \dot{\vec{q}}_0는 임의의 관절 속도 벡터이며, 투영 행렬 (\mathbf{I}_n - \mathbf{J}^+ \mathbf{J})는 자코비안의 영 공간으로의 직교 투영(orthogonal projection)을 수행한다. 이 두 번째 항은 엔드 이펙터 속도에 기여하지 않는 자기 운동(self-motion)이며, 부차적 과제 수행에 활용된다.

자유도 부족의 경우 (n < m이고 \mathbf{J}의 열 계수가 n): 일반적으로 해가 존재하지 않으며, 최소 제곱 의미(least-squares sense)의 근사 해가 \mathbf{J}^+ = (\mathbf{J}^\top \mathbf{J})^{-1} \mathbf{J}^\top에 의해 얻어진다.

특이 구성의 경우: 자코비안의 계수가 최대 계수보다 작아지며, 이 경우 의사 역행렬의 수치적 조건수(condition number)가 발산하므로 감쇠 최소 제곱법(damped least squares) 등의 정칙화(regularization) 기법이 요구된다.

4. 제어 관점에서의 두 공간

관절 공간 속도와 작업 공간 속도의 관계는 로봇 제어 구조에도 직접적 영향을 미친다.

관절 공간 제어(joint-space control)는 관절 변수 \vec{q} 또는 \dot{\vec{q}}를 직접 추종하는 방식이다. 각 관절의 독립적 서보 제어기가 기본 구조를 이룬다.

작업 공간 제어(task-space control 또는 operational-space control)는 엔드 이펙터 자세 \vec{x} 또는 속도 \dot{\vec{x}}를 직접 추종하는 방식이다. 실제 구동은 관절에서 이루어지므로 역속도 사상을 통해 목표 관절 속도 \dot{\vec{q}} = \mathbf{J}^+(\vec{q}) \, \dot{\vec{x}}_d를 계산한 뒤 관절 서보에 전달한다. 이 구조는 Whitney의 해상도 모션 제어(resolved-motion rate control)가 정립한 고전적 틀이다.

경로 계획은 일반적으로 작업 공간에서 수행되며, 실행은 관절 공간에서 이루어지므로 두 공간 사이의 속도 변환이 실시간 제어 루프에서 반복적으로 수행된다.

5. 수치적 계산과 실시간 관점

자코비안과 그 역사상의 계산은 실시간 제어에서 수 kHz 수준의 주기로 반복되어야 하므로 수치적 효율성과 안정성이 중요하다. 자코비안의 구성은 순기구학의 편미분을 해석적으로 유도하거나, 재귀적 기하학 알고리즘(recursive geometric algorithm)을 이용하여 관절별 열 벡터로 직접 구성하는 방법이 존재한다.

역속도 사상의 계산은 일반적으로 자코비안의 명시적 역행렬을 구하기보다는 선형 연립 방정식

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

를 LU 분해, QR 분해, 또는 특이값 분해(singular value decomposition, SVD) 기반으로 풀어 수치적 안정성을 확보한다. 특이 구성 근방에서는 감쇠 최소 제곱법을 적용하여 다음과 같은 정칙화된 해를 사용한다.

\dot{\vec{q}} = \mathbf{J}^\top(\vec{q}) \, \big(\mathbf{J}(\vec{q}) \, \mathbf{J}^\top(\vec{q}) + \lambda^2 \mathbf{I}_m\big)^{-1} \, \dot{\vec{x}}_d

여기서 \lambda > 0은 감쇠 계수이다.

6. 기하학적 해석과 미분 다양체

두 공간의 속도 관계는 미분 다양체(differentiable manifold) 관점에서 보면, 관절 공간 다양체와 작업 공간 다양체 사이의 매끄러운 사상 \vec{f}의 접사상(tangent map, differential) d\vec{f}_{\vec{q}}가 자코비안에 해당한다. 즉 자코비안은 \vec{f}의 미분 구조를 점별로 나타내는 선형 연산자이며, 그 계수와 영 공간은 매니퓰레이터의 국소 운동학적 성질을 결정한다.

이러한 기하학적 관점은 특이점의 정의(자코비안 계수 감소), 매니퓰러빌리티 해석(자코비안 이미지의 타원체), 여유 자유도의 활용(영 공간 투영) 등으로 자연스럽게 확장된다.

7. 학술적 의의

관절 공간 속도와 작업 공간 속도의 관계는 로봇 공학에서 다음과 같은 학술적 의의를 가진다. 첫째, 매니퓰레이터의 비선형 순기구학을 속도 수준의 선형 관계로 환원하여 해석 가능성을 높인다. 둘째, 역기구학의 수치적 반복 해법의 이론적 기초를 제공한다. 셋째, 특이점·여유 자유도·매니퓰러빌리티 등 운동학적 특성 해석의 공통 틀을 마련한다. 넷째, 자코비안 전치를 매개로 하는 힘-토크 쌍대성을 통하여 정적 역학과의 연결을 제공한다. 다섯째, 작업 공간 제어와 관절 공간 제어의 체계적 통합을 가능하게 한다.

본 절에서 정립한 두 공간 사이의 속도 관계는 이후 자코비안의 구체적 유도, 기하학적 해석, 특이점 분석, 역속도 기구학의 다양한 해법, 그리고 매니퓰러빌리티 해석 등의 학술적 전개의 기초가 된다.

8. 출처

  • Whitney, D. E., “Resolved motion rate control of manipulators and human prostheses”, IEEE Transactions on Man-Machine Systems, Vol. 10, No. 2, pp. 47–53, 1969.
  • 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.
  • 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.
  • 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.
  • 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.

9. 버전

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