32.31 특이점(Singularity)의 수학적 정의

특이점은 매니퓰레이터의 자코비안이 정상적인 사상 성질을 잃어버리는 관절 구성을 가리킨다. 이러한 구성에서는 관절 속도와 작업 공간 속도 사이의 일대일 대응이 깨지며, 일정한 작업 공간 방향에 대한 운동 생성 능력이 순간적으로 상실되거나, 거꾸로 일정한 관절 속도가 어떠한 작업 공간 운동도 만들어 내지 못하는 상태가 발생한다. 본 절에서는 특이점의 정의를 자코비안 행렬의 계수 손실로 형식화하고, 이를 특이값·행렬식·미분 사상의 관점에서 동등하게 기술하며, 미분기하학적 함의와 검출 절차, 그리고 운동학·정역학 양면에서의 의미를 차분하게 정리한다.

1. 자코비안 계수 손실에 의한 특이점 정의

매니퓰레이터의 순기구학 사상을 $\vec{x} = \vec{f}(\vec{q})$ 로 두고, 그 야코비안을 $\mathbf{J}(\vec{q}) = \partial \vec{f}/\partial \vec{q} \in \mathbb{R}^{m \times n}$ 이라 한다. 정상 구성에서는 $\operatorname{rank}\,\mathbf{J}(\vec{q}) = \min(m, n)$ 이 성립하지만, 어떤 관절 구성 $\vec{q}^\ast$ 에서 $\operatorname{rank}\,\mathbf{J}(\vec{q}^\ast) < \min(m, n)$ 이 성립하면 그 구성을 기구학적 특이점(kinematic singularity)이라 정의한다. 이러한 정의는 작업 공간 변수와 관절 변수의 차원에 무관하게 적용되며, 직렬 매니퓰레이터, 병렬 기구, 여유 자유도 매니퓰레이터, 이동 로봇의 운동학 사상 모두에 대해 일관된 형식적 기반을 제공한다. 특이점 집합은 일반적으로 관절 공간 $\mathcal{Q}$ 의 측도 영(zero-measure) 부분집합을 이루며, 정상 구성의 여집합으로서 닫힌 부분 다양체 또는 그 합집합의 형태로 나타난다.

특이점에서는 자코비안의 계수 결손 정도 $r = \min(m, n) - \operatorname{rank}\,\mathbf{J}(\vec{q}^\ast)$ 를 특이점의 차원 결손(corank)이라 부르며, 이 값이 클수록 특이성이 더 깊다고 표현한다. $r = 1$ 인 특이점은 단순 특이점, $r \ge 2$ 인 특이점은 다중 특이점으로 분류된다. 자코비안의 계수가 정수 단계로 떨어지므로 특이 집합의 위상은 일반적으로 매끄러운 부분 다양체와 그 경계의 합집합으로 분해되며, 각 성분은 결손 정도에 따라 위계를 이룬다. 직렬 매니퓰레이터에서 $m = n$ 인 경우 정사각 자코비안의 행렬식 $\det \mathbf{J}(\vec{q})$ 의 영점 집합이 특이점 집합과 일치하므로, 특이성은 결국 다항식 또는 삼각함수 다항식의 영점 집합 분석으로 환원된다.

2. 동치 조건과 대수적 표현

자코비안의 특이값 분해 $\mathbf{J}(\vec{q}) = \mathbf{U}(\vec{q})\,\boldsymbol{\Sigma}(\vec{q})\,\mathbf{V}^\top(\vec{q})$ 를 도입하면, $\vec{q}^\ast$ 가 특이점일 조건은 적어도 하나의 특이값이 영이 되는 조건과 동치이다. 즉 $\sigma_{\min}(\mathbf{J}(\vec{q}^\ast)) = 0$ 이 성립하며, 일반적으로 $\sigma_k(\vec{q}^\ast) = \sigma_{k+1}(\vec{q}^\ast) = \cdots = 0$ 이 동시에 성립할 때 차원 결손이 그만큼 증가한다. 조건수의 관점에서는 $\kappa(\mathbf{J}(\vec{q}^\ast)) = \sigma_{\max}/\sigma_{\min} \to \infty$ 로 발산하므로, 특이점은 조건수가 무한대로 가는 극한 구성이며 등방 구성이 위치한 하한 $1$ 과 정반대 극단에 자리한다.

$m = n$ 인 정사각 자코비안의 경우 특이성 조건은 $\det \mathbf{J}(\vec{q}^\ast) = 0$ 의 단일 대수 방정식으로 표현되며, 매니퓰러빌리티 측도 $w(\vec{q}) = \sqrt{\det(\mathbf{J} \mathbf{J}^\top)}$ 도 동시에 영이 된다. $m \neq n$ 인 직사각 자코비안의 경우 정의는 다음과 같이 정리된다. $\operatorname{rank}\,\mathbf{J}(\vec{q}^\ast) < \min(m, n) \iff \det(\mathbf{J} \mathbf{J}^\top)(\vec{q}^\ast) = 0 \iff \det(\mathbf{J}^\top \mathbf{J})(\vec{q}^\ast) = 0$ . 또한 $\mathbf{J}$ 의 모든 $\min(m,n) \times \min(m,n)$ 부분 행렬의 행렬식이 동시에 영이 되는 조건과도 동치이며, 이는 그라스만 다양체상의 영점 조건으로 일반화된다.

3. 미분기하학적 해석과 임계점 구조

순기구학 사상 $\vec{f}: \mathcal{Q} \to \mathcal{X}$ 의 미분 $d\vec{f}_{\vec{q}}: T_{\vec{q}} \mathcal{Q} \to T_{\vec{f}(\vec{q})} \mathcal{X}$ 를 고려하면, 자코비안 $\mathbf{J}(\vec{q})$ 는 이 미분의 좌표 표현이다. 특이점은 $d\vec{f}_{\vec{q}^\ast}$ 가 그 정의역에서 가능한 최대 계수에 도달하지 못하는 점, 즉 미분 위상수학적 의미의 임계점(critical point)이다. 임계점의 상 $\vec{f}(\vec{q}^\ast)$ 는 임계값(critical value)으로 불리며, Sard 정리에 의해 $\mathcal{X}$ 안에서 측도 영의 부분집합을 이룬다. 따라서 일반적인 작업 공간 위치 $\vec{x}$ 에 대응하는 역기구학 해는 정상 구성으로 구성되지만, 임계값에 가까워질수록 해의 분기 구조가 변화한다.

특이점에서는 미분 사상의 영공간과 상공간이 비자명한 구조를 가진다. 영공간 $\ker \mathbf{J}(\vec{q}^\ast)$ 에 속하는 관절 속도는 엔드 이펙터에 어떠한 운동도 유발하지 않는 내부 운동을 형성하며, 직교 여집합 $\mathrm{range}\,\mathbf{J}^\top(\vec{q}^\ast)$ 바깥의 작업 공간 방향은 어떠한 관절 속도로도 즉시 생성될 수 없는 차단 방향(blocked direction)이 된다. 이 구조는 특이성의 본질이 단순한 수치적 악조건이 아니라 운동 사상의 위상적 변화임을 보여 준다. Whitney와 후속 연구자들이 정립한 일반 위치(generic position)의 임계점 분류 이론은 매끄러운 사상의 임계점이 폴드(fold), 큐스프(cusp) 등의 표준 형태로 분류됨을 보였으며, 이러한 표준형은 매니퓰레이터의 특이점 근방 구조를 분석하는 출발점으로 활용된다.

4. 특이점 검출과 분류의 대수적 절차

직렬 매니퓰레이터의 특이점 검출은 일반적으로 자코비안의 행렬식 또는 매니퓰러빌리티 측도가 관절 변수의 다항식 또는 삼각함수 다항식이 됨을 이용하여 그 영점 집합을 식별하는 절차로 수행된다. $\det \mathbf{J}(\vec{q}) = 0$ 또는 $\det(\mathbf{J} \mathbf{J}^\top)(\vec{q}) = 0$ 의 인수분해를 통해 특이 집합은 여러 성분으로 분해되며, 각 성분은 손목, 팔꿈치, 어깨와 같이 기구학적 의미를 가지는 부분 다양체로 해석된다. 그뢰브너 기저 기법, 종결식(resultant) 기법, 호모토피 연속법 등은 다항계의 영점 집합을 체계적으로 다루는 도구로 사용된다. 병렬 기구의 경우에는 입력 자코비안과 출력 자코비안을 분리하여 특이성을 정의하는 Gosselin과 Angeles의 1990년 논문 “Singularity analysis of closed-loop kinematic chains“의 분류가 표준적 틀로 자리잡았으며, 이는 직접 특이성, 역특이성, 결합 특이성의 세 가지 유형으로 특이 집합을 분해한다.

수치적 검출에서는 자코비안의 최소 특이값 $\sigma_{\min}(\mathbf{J}(\vec{q}))$ , 매니퓰러빌리티 측도 $w(\vec{q})$ , 조건수 $\kappa(\mathbf{J}(\vec{q}))$ 등의 스칼라 지표가 특이성의 근접도를 정량화하는 척도로 사용된다. 부동 소수점 연산의 한계로 인해 엄밀한 영점은 거의 검출되지 않으므로, 사용자는 임계 임곗값 $\epsilon > 0$ 을 도입하여 $\sigma_{\min}(\mathbf{J}(\vec{q})) < \epsilon$ 또는 $w(\vec{q}) < \epsilon$ 을 만족하는 영역을 특이 근방으로 정의한다. 이러한 임곗값 기반 정의는 엄밀한 특이점 정의를 일반화한 실용적 형태이며, 제어 알고리즘이 수치 불안정성을 회피하기 위해 동작 모드를 전환하는 기준으로 활용된다.

5. 운동학적·역학적 함의와 본 절의 정리

특이점에서의 운동학적 결과는 두 측면에서 동시에 발현된다. 한편으로는 작업 공간 운동의 일정 방향이 순간적으로 생성 불가능해지므로, 그 방향에 대한 속도 명령은 어떠한 관절 속도 조합으로도 추종될 수 없다. 다른 한편으로는 영공간이 비자명해지므로 관절 변수가 변화하면서도 엔드 이펙터가 정지 상태에 머무르는 내부 운동이 가능해진다. 이 두 결과는 자코비안의 사상 구조가 가역적이지 않다는 단일한 사실의 두 측면이며, 가상 일의 원리에 따라 정역학적 측면에서도 동시에 나타난다. 즉 작업 공간 힘의 일정 방향이 관절 토크에 의해 더 이상 평형화되지 못하고, 어떤 비자명한 관절 토크 조합은 엔드 이펙터에 어떠한 힘 또는 모멘트도 발생시키지 않는 자기 평형(self-equilibrated) 상태에 놓인다.

본 절에서 다룬 특이점의 수학적 정의는 자코비안의 계수 결손이라는 단일한 핵심 조건으로 요약되며, 이 조건은 특이값 영점, 행렬식 영점, 부분 행렬식의 동시 영점, 매니퓰러빌리티 측도 영점 등 다양한 동치 형태로 변환된다. 미분기하학적 해석은 특이점이 순기구학 사상의 임계점이라는 위상적 의미를 부여하며, 이는 운동학과 정역학의 양면에서 사상 구조의 가역성 상실로 이어진다. 본 절은 이후에 다룰 다양한 특이점의 분류, 근방 거동, 회피 전략, 그리고 자코비안 기반 수치적 처리 기법의 출발점으로 기능한다.

6. 출처

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.
Yoshikawa, T., “Manipulability of robotic mechanisms”, International Journal of Robotics Research, Vol. 4, No. 2, pp. 3–9, 1985.
Gosselin, C. and Angeles, J., “Singularity analysis of closed-loop kinematic chains”, IEEE Transactions on Robotics and Automation, Vol. 6, No. 3, pp. 281–290, 1990.
Park, F. C. and Brockett, R. W., “Kinematic dexterity of robotic mechanisms”, International Journal of Robotics Research, Vol. 13, No. 1, pp. 1–15, 1994.
Murray, R. M., Li, Z., and Sastry, S. S., A Mathematical Introduction to Robotic Manipulation, CRC Press, 1994.
Sciavicco, L. and Siciliano, B., Modelling and Control of Robot Manipulators, 2nd edition, Springer, 2000.
Merlet, J.-P., Parallel Robots, 2nd edition, Springer, 2006.

7. 버전

문서 버전: 2.0
작성일: 2026-04-26