32.27 매니퓰러빌리티 지수(Manipulability Index)의 정의
매니퓰러빌리티 지수(manipulability index)는 Yoshikawa가 1985년 논문 “Manipulability of robotic mechanisms“에서 제안한 스칼라 지표로, 매니퓰레이터의 국소 운동학적 능력을 단일 수치로 정량화한다. 이 지수는 매니퓰러빌리티 타원체의 부피와 비례하며, 매니퓰레이터의 성능 평가, 특이점 회피, 구성 최적화, 구조 설계에서 핵심적 역할을 한다. 본 절에서는 매니퓰러빌리티 지수의 학술적 정의, 수학적 성질, 기하학적 의미, 특이값·타원체 부피와의 관계, 정방 자코비안과 비정방 자코비안에서의 특수 형태, 그리고 관련 지표와의 비교를 체계적으로 다룬다.
1. 매니퓰러빌리티 지수의 학술적 정의
자코비안 \mathbf{J}(\vec{q}) \in \mathbb{R}^{m \times n}을 가진 매니퓰레이터의 구성 \vec{q}에서의 Yoshikawa 매니퓰러빌리티 지수는 다음과 같이 정의된다.
w(\vec{q}) = \sqrt{\det(\mathbf{J}(\vec{q}) \, \mathbf{J}^\top(\vec{q}))}
이 정의는 완전 행 계수 구성(\text{rank}(\mathbf{J}) = m)에서 양수 값을 가지며, 특이점에서는 \mathbf{J} \mathbf{J}^\top의 행렬식이 영이 되어 w = 0이 된다.
매니퓰러빌리티 지수는 다음의 수학적 성질을 만족한다.
비음성: w(\vec{q}) \geq 0. \mathbf{J} \mathbf{J}^\top이 양반정치 행렬이므로 행렬식이 비음수이다.
특이점에서의 영: w(\vec{q}) = 0 \iff \text{rank}(\mathbf{J}(\vec{q})) < m. 즉 매니퓰러빌리티 지수가 영인 조건은 기구학적 특이점의 정의와 동치이다.
연속성과 미분 가능성: 순기구학 사상이 매끄러우면 w는 관절 변수의 매끄러운 함수이다. 이러한 미분 가능성은 경사법 기반 최적화에서 중요하다.
관절 변수 불변성: 매니퓰러빌리티 지수는 일반화된 기구학적 불변량이 아니며, 관절 변수의 스케일링 등 기구학적 매개변수화에 의존한다.
32.27.2 특이값과 매니퓰러빌리티 지수의 관계
자코비안의 특이값 분해 \mathbf{J} = \mathbf{U} \mathbf{\Sigma} \mathbf{V}^\top에서 특이값 \sigma_1 \geq \sigma_2 \geq \cdots \geq \sigma_m \geq 0을 사용하면 매니퓰러빌리티 지수는 특이값의 곱으로 표현된다.
w(\vec{q}) = \prod_{i=1}^{m} \sigma_i(\vec{q})
이 결과는 \det(\mathbf{J} \mathbf{J}^\top) = \det(\mathbf{U} \mathbf{\Sigma} \mathbf{\Sigma}^\top \mathbf{U}^\top) = \det(\mathbf{\Sigma} \mathbf{\Sigma}^\top) = \prod_{i=1}^{m} \sigma_i^2로부터 w = \sqrt{\prod \sigma_i^2} = \prod \sigma_i가 성립함을 통해 유도된다.
이 결과는 다음의 기하학적 해석을 제공한다.
모든 특이값이 크면 매니퓰러빌리티가 크다: 각 방향으로의 속도 이득이 균일하게 크면 매니퓰러빌리티 지수가 크다.
한 특이값이라도 작으면 매니퓰러빌리티가 작다: 하나의 특이값이 영에 가까우면 곱의 형태로 매니퓰러빌리티가 급격히 감소한다. 이는 매니퓰러빌리티 지수가 특이점에 민감하게 반응하는 ’병목 지표’의 성격을 가짐을 의미한다.
특이점에서의 영: 하나 이상의 특이값이 정확히 영이 되면 w = 0이며, 이는 특이점의 수학적 정의이다.
2. 매니퓰러빌리티 타원체의 부피
매니퓰러빌리티 지수는 매니퓰러빌리티 타원체의 m차원 부피와 비례한다. 구체적으로 m차원 타원체의 부피 공식에 따라
\text{Vol}(\mathcal{E}_v(\vec{q})) = V_m \cdot \prod_{i=1}^{m} \sigma_i = V_m \cdot w(\vec{q})
여기서 V_m은 단위 구의 m차원 부피이다(V_2 = \pi, V_3 = \frac{4}{3} \pi 등). 따라서 매니퓰러빌리티 지수는 상수 배수로 타원체 부피와 일치한다.
이러한 부피 해석은 매니퓰러빌리티 지수의 직관적 이해를 제공한다. 큰 w는 매니퓰레이터가 “다양한 방향으로 큰 속도를 생성할 수 있는” 구성임을 의미하며, 작은 w는 운동 능력이 저하된 구성을 나타낸다.
32.27.4 정방 자코비안에서의 단순화
관절 수와 작업 공간 차원이 일치하는 정방 매니퓰레이터(n = m)의 경우 자코비안이 정방 행렬이고, 매니퓰러빌리티 지수는 자코비안의 행렬식의 절댓값과 같다.
w(\vec{q}) = |\det \mathbf{J}(\vec{q})|
이 단순화는 \det(\mathbf{J} \mathbf{J}^\top) = (\det \mathbf{J})^2의 관계로부터 직접 유도된다.
정방 자코비안의 경우 매니퓰러빌리티 지수는 계산이 단순하며, 해석적 조사도 용이하다. 예를 들어 평면 2자유도 매니퓰레이터의 매니퓰러빌리티 지수는
w(\theta_1, \theta_2) = |\det \mathbf{J}(\theta_1, \theta_2)| = a_1 \, a_2 \, |\sin \theta_2|
로 구해진다(링크 길이 a_1, a_2, 두 번째 관절 각 \theta_2). 이는 두 번째 관절이 \theta_2 = 0 또는 \theta_2 = \pi(팔이 완전히 펴지거나 접힌 상태)에서 영이 됨을 보이며, 이는 평면 2자유도 매니퓰레이터의 특이점 조건과 일치한다.
32.27.5 비정방 자코비안에서의 의미
여유 자유도 매니퓰레이터(n > m)의 경우 자코비안이 비정방이지만 \mathbf{J} \mathbf{J}^\top \in \mathbb{R}^{m \times m}은 정방이며, 완전 행 계수에서 양정치이다. 매니퓰러빌리티 지수 w = \sqrt{\det(\mathbf{J} \mathbf{J}^\top)}는 동일한 의미를 유지한다.
여유 자유도의 존재는 동일 엔드 이펙터 자세에 대해 여러 관절 구성이 가능함을 의미하므로, 매니퓰러빌리티 지수는 각 구성에 따라 다른 값을 가진다. 이러한 구성 의존성은 자기 운동을 활용한 매니퓰러빌리티 최대화의 이론적 기초가 된다.
자유도 부족 매니퓰레이터(m > n)의 경우 \mathbf{J} \mathbf{J}^\top의 계수가 n이므로 항상 \det(\mathbf{J} \mathbf{J}^\top) = 0이 되어 w = 0이다. 따라서 Yoshikawa 지수는 자유도 부족 매니퓰레이터에는 직접 적용되지 않으며, \sqrt{\det(\mathbf{J}^\top \mathbf{J})}의 변형 형태가 사용된다.
32.27.6 매니퓰러빌리티 지수의 미분 가능성
매니퓰러빌리티 지수는 정규 구성에서 관절 변수에 대해 미분 가능하며, 경사는 다음과 같이 계산된다.
\frac{\partial w}{\partial q_j}(\vec{q}) = w(\vec{q}) \cdot \text{tr}\!\left[ (\mathbf{J} \mathbf{J}^\top)^{-1} \, \frac{\partial (\mathbf{J} \mathbf{J}^\top)}{\partial q_j} \right]
여기서 \text{tr}은 행렬의 대각합이다. \partial (\mathbf{J} \mathbf{J}^\top) / \partial q_j = (\partial \mathbf{J} / \partial q_j) \mathbf{J}^\top + \mathbf{J} (\partial \mathbf{J}^\top / \partial q_j)이므로, 자코비안의 관절 변수에 대한 미분이 경사 계산에 필요하다.
이 경사는 자코비안의 시간 미분 계산 알고리즘과 유사한 방식으로 효율적으로 계산 가능하며, 매니퓰러빌리티 최대화 자기 운동의 방향을 제공한다.
특이점 근방에서는 w 자체가 영에 접근하므로 경사의 크기가 급격히 변화하며, 수치적으로 불안정해질 수 있다. 이러한 경우 로그 매니퓰러빌리티 \log w의 경사를 사용하거나 감쇠 정칙화를 적용한다.
3. 선속도·각속도 매니퓰러빌리티 지수
기하학적 자코비안의 선속도 블록과 각속도 블록으로부터 분리된 매니퓰러빌리티 지수가 정의된다.
선속도 매니퓰러빌리티 지수:
w_v(\vec{q}) = \sqrt{\det(\mathbf{J}_v(\vec{q}) \mathbf{J}_v^\top(\vec{q}))}
각속도 매니퓰러빌리티 지수:
w_\omega(\vec{q}) = \sqrt{\det(\mathbf{J}_\omega(\vec{q}) \mathbf{J}_\omega^\top(\vec{q}))}
분리된 지수는 위치 매니퓰러빌리티와 방향 매니퓰러빌리티를 단위 일관성 있게 평가한다. 선속도 지수의 단위는 m³(3차원 부피), 각속도 지수의 단위는 무차원(rad가 무차원으로 간주됨)이다.
혼합 응용에서 두 지수의 결합 지표는 다음과 같은 형태가 사용된다.
- 곱의 형태: w_{\text{combined}} = w_v^{\alpha_v} \cdot w_\omega^{\alpha_\omega} (가중 기하 평균).
- 가중 평균: w_{\text{combined}} = \alpha_v w_v + \alpha_\omega w_\omega (가중 산술 평균).
가중치 \alpha_v, \alpha_\omega는 응용 맥락에 따라 결정된다.
4. 관련 지표와의 비교
매니퓰러빌리티 지수 외에도 유사한 성능 평가 지표들이 제안되어 왔다.
4.1 조건수
조건수 \kappa(\mathbf{J}) = \sigma_{\max} / \sigma_{\min}는 이방성을 측정하며, 매니퓰러빌리티 지수와 상호 보완적이다. 매니퓰러빌리티는 타원체 부피를 측정하고 조건수는 타원체 형상의 이방성을 측정한다. 두 지표를 함께 사용하면 성능 평가가 보다 포괄적이다.
4.2 최소 특이값
최소 특이값 \sigma_{\min}(\mathbf{J})은 특이점으로부터의 ’거리’를 측정하는 또 다른 지표이다. 매니퓰러빌리티는 모든 특이값의 곱이므로 여러 특이값이 모두 영으로 수렴할 때 급격히 감소하지만, 최소 특이값은 단일 특이값만 추적한다.
4.3 손재주 지표
Klein과 Blaho(1987)의 손재주 지표(dexterity measure)는 매니퓰러빌리티의 변형으로, 작업 공간의 전역적 성능 평가에 초점을 둔 지표이다.
4.4 국소 조건 지표
Angeles와 López-Cajún(1992)의 조건 지표는 조건수의 역수를 전역 작업 공간에 걸쳐 평균한 지표이며, 매니퓰레이터 설계 최적화에 사용된다.
4.5 동역학적 매니퓰러빌리티
Yoshikawa(1985)의 또 다른 연구인 동역학적 매니퓰러빌리티는 기구학적 매니퓰러빌리티를 관성 행렬로 가중하여 토크 생성 능력까지 반영한다. 이는 동역학적 성능을 평가하는 지표이다.
5. 매니퓰러빌리티 지수의 응용
매니퓰러빌리티 지수는 다음의 실무 응용에 활용된다.
부차적 과제 설계: 영 공간 투영 프레임워크에서 \nabla w의 방향으로 자기 운동을 실행하여 매니퓰러빌리티를 최대화.
경로 최적화: 엔드 이펙터 경로 상에서 매니퓰러빌리티의 평균 또는 최소값을 최대화하는 궤적 선택.
매니퓰레이터 구조 설계: 링크 길이와 관절 배치의 설계 변수를 전역 매니퓰러빌리티 기준으로 최적화.
작업 배치 설계: 공작물 또는 작업 대상의 위치를 매니퓰러빌리티가 높은 영역에 배치.
성능 비교: 서로 다른 매니퓰레이터 구조의 성능 비교에서 표준 지표로 활용(링크 길이 정규화 후).
여유 자유도 활용: 7자유도 매니퓰레이터의 팔 각도 선택 기준.
6. 한계와 주의 사항
매니퓰러빌리티 지수의 유용성에도 불구하고 다음의 한계와 주의 사항이 존재한다.
이방성 무시: 매니퓰러빌리티 지수는 타원체의 부피만을 측정하며, 이방성을 반영하지 않는다. 크게 이방적인 타원체라도 부피가 크면 매니퓰러빌리티가 높게 나타난다. 따라서 조건수와 함께 사용해야 한다.
단위 일관성: 기하학적 자코비안의 전체 블록 기반 매니퓰러빌리티는 선속도와 각속도의 단위 혼재 문제가 있다. 분리된 지수 또는 단위 정규화가 요구된다.
설계 변수 스케일 의존성: 매니퓰러빌리티 지수는 링크 길이 등의 설계 변수의 스케일에 민감하므로, 비교 시 정규화가 필요하다.
국소적 지표: 매니퓰러빌리티는 특정 구성에서의 국소 성능을 기술하며, 작업 공간 전체의 성능은 별도로 평가해야 한다.
동역학적 효과 미반영: 기구학 기반 지표이므로 관성, 중력, 토크 한계 등의 동역학적 효과를 반영하지 않는다.
7. 학술적 의의
본 절에서 정립한 매니퓰러빌리티 지수의 정의는 로봇 공학에서 다음의 학술적 의의를 가진다. 첫째, Yoshikawa의 고전적 정의 w = \sqrt{\det(\mathbf{J} \mathbf{J}^\top)}는 매니퓰레이터의 국소 운동학적 능력을 단일 수치로 정량화하는 표준 지표이다. 둘째, 특이값의 곱으로의 분해는 매니퓰러빌리티의 병목 지표적 성격을 명시하며, 특이점 식별의 정량적 기준을 제공한다. 셋째, 매니퓰러빌리티 타원체 부피와의 직접적 연결은 기하학적 해석을 제공한다. 넷째, 미분 가능성과 명시적 경사 공식은 경사 기반 최적화 응용의 기초이다. 다섯째, 선속도·각속도 분리 지수는 단위 일관성을 확보하는 실무적 도구이다. 여섯째, 조건수·최소 특이값 등 관련 지표와의 비교는 성능 평가의 다차원적 접근을 가능하게 한다. 일곱째, 매니퓰러빌리티 지수를 활용한 여유 자유도 제어, 경로 최적화, 구조 설계는 현대 로봇 공학의 핵심 응용이며, 본 절의 정의가 이러한 응용의 수학적 토대를 이룬다. 여덟째, 본 절은 후속 절의 매니퓰러빌리티 지수를 이용한 성능 평가와 조건수·등방성 해석의 공통 기초를 제공한다.
8. 출처
- Yoshikawa, T., “Manipulability of robotic mechanisms”, International Journal of Robotics Research, Vol. 4, No. 2, pp. 3–9, 1985.
- Yoshikawa, T., “Dynamic manipulability of robot manipulators”, Journal of Robotic Systems, Vol. 2, No. 1, pp. 113–124, 1985.
- Yoshikawa, T., Foundations of Robotics: Analysis and Control, MIT Press, 1990.
- Klein, C. A. and Blaho, B. E., “Dexterity measures for the design and control of kinematically redundant manipulators”, International Journal of Robotics Research, Vol. 6, No. 2, pp. 72–83, 1987.
- Angeles, J. and López-Cajún, C. S., “Kinematic isotropy and the conditioning index of serial robotic manipulators”, International Journal of Robotics Research, Vol. 11, No. 6, pp. 560–571, 1992.
- 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.
- Chiu, S. L., “Task compatibility of manipulator postures”, International Journal of Robotics Research, Vol. 7, No. 5, pp. 13–21, 1988.
- 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.
9. 버전
- 문서 버전: 1.1
- 작성일: 2026-04-19