33.2 자코비안 행렬의 계수 결손과 특이점의 관계

자코비안 행렬의 계수 결손(rank deficiency)은 특이점 개념의 수학적 핵심이다. 관절 공간에서 작업 공간으로의 미분 사상이 선형 대수적으로 퇴화하는 정도를 정량화함으로써, 특이점의 심각성, 유형, 회피 전략의 공통 기초가 마련된다. 본 절에서는 계수 결손의 엄밀한 정의, 등가 조건, 차수 분류, 선형 대수 관점의 해석, 수치적 판별, 결손의 기구학적 원인, 복구 기법을 체계적으로 해라체로 기술한다.

1. 계수와 계수 결손의 기초

1.1 계수의 정의

자코비안 $\mathbf{J}(\vec{q}) \in \mathbb{R}^{m \times n}$ 의 계수(rank)는 선형 독립인 열 벡터(또는 행 벡터)의 최대 수로 정의되며, 열 공간 $\operatorname{range}(\mathbf{J})$ 와 행 공간 $\operatorname{range}(\mathbf{J}^\top)$ 의 차원과 일치한다. 열 계수와 행 계수는 항상 같다.

1.2 완전 계수 조건

자코비안이 완전 계수(full rank)란 $\operatorname{rank}(\mathbf{J}) = \min(m, n)$ 인 경우를 말한다. 행 완전 계수( $m \le n$ )는 $\mathbf{J} \mathbf{J}^\top$ 이 가역임과 동치이고, 열 완전 계수( $m \ge n$ )는 $\mathbf{J}^\top \mathbf{J}$ 이 가역임과 동치이다.

1.3 계수 결손의 정의

계수 결손은 $\operatorname{rank}(\mathbf{J}(\vec{q})) < \min(m, n)$ 인 상태로, 특이점의 정의와 동일하다.

$\vec{q} \in \mathcal{S} \iff \operatorname{rank}\bigl( \mathbf{J}(\vec{q}) \bigr) < \min(m, n)$

계수 결손은 관절 구성에 의존하는 국소적 성질이며, 관절 공간의 일반점(generic point)에서는 완전 계수가 성립한다.

33.2.1.4 일반 계수

일반적 구성(generic configuration)에서 매니퓰레이터 자코비안은 $\min(m, n)$ 의 일반 계수를 가지며, 이 일반 계수를 달성하지 못하는 구성이 특이 구성이다. Strang의 Introduction to Linear Algebra 제5판(Wellesley-Cambridge Press, 2016)과 같은 선형 대수 표준 교과서가 계수 이론의 기본 틀을 제공한다.

33.2.2 계수 결손의 차수

33.2.2.1 결손 차수의 정의

계수 결손의 차수(corank)는 다음과 같이 정의된다.

$r_d(\vec{q}) = \min(m, n) - \operatorname{rank}\bigl( \mathbf{J}(\vec{q}) \bigr)$

결손 차수는 특이점의 심각성을 정량화하는 기본 지표이다.

1.4 1차 결손

$r_d = 1$ 인 경우를 단순 특이점(simple singularity) 또는 1차 결손이라 부른다. 엔드 이펙터가 작업 공간의 한 방향으로만 운동할 수 없는 상태이며, 대부분의 기구학적 특이점이 이 범주에 속한다.

1.5 고차 결손

$r_d \ge 2$ 인 경우는 고차 특이점(higher-order singularity)이다. 엔드 이펙터가 다수의 작업 공간 방향으로 순간 운동할 수 없는 상태이며, 작업 공간 해석 관점에서 더 심각한 현상을 유발한다.

1.6 최고 차수 결손

이론적으로 가능한 최대 결손 차수는 $\min(m, n)$ 이나, 실제 매니퓰레이터에서는 링크 구조와 관절 배치의 제약으로 인하여 일반적으로 1 또는 2를 초과하지 않는다.

1.7 결손 차수의 여차원

관절 공간 내에서 결손 차수 $r_d$ 를 갖는 구성의 집합은 일반적으로 여차원(codimension) $r_d^2$ 의 부분 다양체를 이룬다. 이는 대수 다양체 이론의 결과이며, 고차 결손이 드물게 발생하는 이유를 설명한다.

2. 등가 수학적 조건

2.1 행렬식 기반 조건

정사각 자코비안( $m = n$ )의 경우, 계수 결손은 다음 조건과 동치이다.

$\det \mathbf{J}(\vec{q}) = 0$

33.2.3.2 그람 행렬식 조건

여유 자유도( $n > m$ ) 매니퓰레이터의 경우, 계수 결손은 그람 행렬식 조건과 동치이다.

$\det\bigl( \mathbf{J}(\vec{q}) \mathbf{J}^\top(\vec{q}) \bigr) = 0$

자유도 부족( $m > n$ )의 경우에는 $\det(\mathbf{J}^\top \mathbf{J}) = 0$ 이 적용된다.

2.2 특이값 조건

계수 결손은 특이값 분해 $\mathbf{J} = \mathbf{U} \boldsymbol{\Sigma} \mathbf{V}^\top$ 에서 영인 특이값의 수로 정량화된다. 영인 특이값의 개수가 결손 차수 $r_d$ 이다.

2.3 소행렬식 조건

계수가 $r$ 이하임은 $\mathbf{J}$ 의 모든 $(r+1) \times (r+1)$ 소행렬식이 영임과 동치이다. 이 조건은 대수 기하학적 관점에서 특이 다양체의 정의 방정식을 제공한다.

2.4 매니퓰러빌리티 조건

$r_d \ge 1$ 은 매니퓰러빌리티 지수 $w(\vec{q}) = \sqrt{\det(\mathbf{J}\mathbf{J}^\top)} = 0$ 과 동치이다.

3. 선형 대수 관점의 해석

3.1 행 공간과 열 공간의 차원

완전 계수에서 행 공간과 열 공간의 차원은 모두 $\min(m, n)$ 이다. 계수 결손 시 두 공간의 차원이 동시에 $r_d$ 만큼 감소한다.

3.2 상공간의 축소

자코비안의 상공간 $\operatorname{range}(\mathbf{J}) \subseteq \mathbb{R}^m$ 의 차원이 감소한다. 이는 엔드 이펙터의 순간 운동이 가능한 작업 공간 방향의 집합이 축소됨을 의미한다.

3.3 영공간의 확장

자코비안의 영공간 $\operatorname{null}(\mathbf{J}) \subseteq \mathbb{R}^n$ 의 차원이 확장된다. 차원 정리에 의하여

$\dim \operatorname{null}(\mathbf{J}) = n - \operatorname{rank}(\mathbf{J})$

이 성립하므로, 계수 저하 시 영공간 차원이 증가한다. 비여유 자유도 매니퓰레이터에서 정상 구성에는 영공간이 자명( $\{\vec{0}\}$ )이나, 특이점에서는 비자명한 영공간이 출현한다.

33.2.4.4 직교 분해

관절 공간 $\mathbb{R}^n$ 은 $\operatorname{null}(\mathbf{J})$ 과 $\operatorname{range}(\mathbf{J}^\top)$ 의 직교합으로 분해된다.

$\mathbb{R}^n = \operatorname{null}(\mathbf{J}) \oplus \operatorname{range}(\mathbf{J}^\top)$

작업 공간 $\mathbb{R}^m$ 은 $\operatorname{range}(\mathbf{J})$ 과 $\operatorname{null}(\mathbf{J}^\top)$ 의 직교합으로 분해된다. 계수 결손 시 각 부분 공간의 차원이 달라진다.

4. 특이값과 결손의 관계

4.1 특이값 순서

특이값은 $\sigma_1 \ge \sigma_2 \ge \cdots \ge \sigma_{\min(m,n)} \ge 0$ 으로 순서화된다. 완전 계수 구성에서 모든 특이값이 양수이며, 결손 차수 $r_d$ 의 구성에서 가장 작은 $r_d$ 개의 특이값이 영이 된다.

4.2 손실 방향

영 특이값에 대응하는 좌특이벡터 $\vec{u}_i$ 는 작업 공간의 손실 운동 방향(엔드 이펙터가 순간적으로 운동할 수 없는 방향)이다.

4.3 영공간 생성 벡터

영 특이값에 대응하는 우특이벡터 $\vec{v}_i$ 는 $\operatorname{null}(\mathbf{J})$ 의 기저 벡터이며, 엔드 이펙터 운동에 기여하지 않는 관절 속도 방향을 제공한다.

4.4 근접 결손 상태

실무에서는 엄밀한 계수 결손과 함께 근접 결손 상태(near-rank-deficient state)가 중요하다. 최소 특이값이 영에 가까운(그러나 영은 아닌) 상태에서 수치적으로 계수 결손과 유사한 거동이 나타난다.

5. 결손의 기구학적 원인

5.1 열 벡터 선형 종속

자코비안 열 벡터 $\vec{J}_i = \mathbf{J} \vec{e}_i$ 의 선형 종속은 여러 관절 속도가 동일한 엔드 이펙터 속도 성분을 생성함을 의미한다. 이는 공선 정렬, 축 중첩, 축 교차와 같은 기구학적 조건에서 발생한다.

5.2 행 벡터 선형 종속

자코비안 행 벡터의 선형 종속은 특정 작업 공간 방향이 관절 공간 운동으로 생성되지 않음을 의미한다. 스크류 이론적으로는 해당 방향이 모든 관절 트위스트와 직교하는 상호 렌치에 해당한다.

5.3 공선 정렬

인접 관절 축 또는 링크가 공선 정렬을 이루면 자코비안의 열들이 동일 방향 성분을 공유하여 계수 저하가 발생한다.

5.4 축 교차

특정 관절 축이 엔드 이펙터 위치(또는 손목 중심)와 교차하면 해당 관절의 위치 기여 벡터가 영이 되어 결손이 발생한다.

5.5 축 중첩

구형 손목의 제1 축과 제3 축이 공선 정렬되면 자세 자코비안 블록의 계수가 저하된다. 이는 손목 특이점의 대표적 기전이다.

6. 수치적 판별

6.1 특이값 임계치

수치적 환경에서는 엄밀한 영 특이값 대신 임계치 $\epsilon$ 이하의 특이값을 영으로 간주한다.

$\operatorname{rank}_\epsilon(\mathbf{J}) = \#\{ \sigma_i : \sigma_i > \epsilon \}$

임계치 $\epsilon$ 의 선택은 기계 엡실론 $\epsilon_{\mathrm{mach}}$ 와 행렬 크기의 곱에 비례하여 설정된다.

33.2.7.2 조건수 기반 판별

조건수 $\kappa(\mathbf{J}) = \sigma_{\max}/\sigma_{\min}$ 이 큰 경우 수치적으로 계수 결손과 유사한 거동을 보인다. 실무 임계치로는 $\kappa > 10^6$ 또는 $\kappa > 10^8$ 이 사용된다.

33.2.7.3 실무적 지표

매니퓰러빌리티 지수 $w$ , 최소 특이값 $\sigma_{\min}$ , 조건수 $\kappa$ 가 근접 결손 상태의 실시간 감지에 사용된다. 각 지표는 계산 비용과 해석적 명확성 면에서 차이를 보인다.

33.2.7.4 Golub과 Van Loan의 수치 해석

Golub와 Van Loan의 Matrix Computations 제4판(Johns Hopkins University Press, 2013)은 수치적 계수 판별, 근접 결손 분석, 조건수 평가의 표준 이론을 제공한다.

33.2.8 결손의 정역학적 쌍대성

33.2.8.1 자코비안 전치의 계수 결손

자코비안 전치의 계수도 동시에 감소한다. $\operatorname{rank}(\mathbf{J}) = \operatorname{rank}(\mathbf{J}^\top)$ 이므로, 속도 사상의 계수 결손과 힘 사상의 계수 결손이 같은 관절 구성에서 발생한다.

33.2.8.2 힘 생성 불가 방향

$\operatorname{null}(\mathbf{J}^\top)$ 의 원소는 관절 토크를 생성하지 않는 렌치 방향이며, 특이점에서 이 방향 외부 힘은 구조적으로 지지된다.

33.2.8.3 속도-힘 쌍대

속도 사상에서 손실된 방향과 힘 사상에서 생성 불가한 방향은 동일 벡터이다. 이는 가상 일 원리에 의한 직접적 귀결이다.

33.2.9 결손 분류와 특이 유형

33.2.9.1 경계 특이점의 결손

경계 특이점은 일반적으로 $r_d = 1$ 의 결손을 가지며, 작업 공간 경계 법선 방향의 손실 운동 방향으로 특징지어진다.

33.2.9.2 내부 특이점의 결손

내부 특이점의 결손 차수는 기구학적 구조에 따라 다양하다. 손목 특이점과 어깨 특이점은 일반적으로 $r_d = 1$ 이며, 특수한 교차 조건에서 $r_d = 2$ 의 고차 결손이 나타날 수 있다.

33.2.9.3 결합 특이점

어깨와 팔꿈치가 동시에 특이 구성에 있거나, 어깨와 손목이 동시에 특이 구성에 있으면 결합 특이점이 발생하며, $r_d = 2$ 이상의 고차 결손을 유발한다.

33.2.9.4 병렬 기구 Type II의 결손

병렬 기구의 Type II 특이점은 순기구학 자코비안의 계수 저하로 정의되며, 플랫폼이 순간적으로 추가 자유도를 획득한다. 이는 결손 차수와 얻는 자유도 수가 같음을 보여준다.

33.2.10 결손 해소 기법

33.2.10.1 관절 구성 변경

관절 구성 $\vec{q}$ 를 특이점으로부터 이동시키면 계수 결손이 해소된다. 이는 경로 계획 또는 자세 재배치를 통해 실현된다.

33.2.10.2 여유 자유도 자체 운동

여유 자유도 매니퓰레이터에서 영공간 자체 운동을 이용하여 엔드 이펙터 자세를 유지한 채 관절 구성을 특이점에서 벗어나게 한다.

33.2.10.3 구조적 설계

링크 길이, 관절 오프셋, 축 배치의 재설계는 특이점 집합의 구조를 근본적으로 변경한다. 오프셋 손목, 오프셋 어깨 등이 대표적 예이다.

33.2.10.4 강건 역기구학

감쇠 최소 제곱, 선택적 감쇠, 가중 의사 역행렬 등 정규화 기법은 근접 결손 상태에서 수치 안정적 해를 제공한다. 이들은 결손 자체를 해소하지는 않으나 결손의 효과를 완화한다.

33.2.11 결손의 해석 절차

33.2.11.1 대수적 절차

자코비안의 해석적 형태를 유도한다.
행렬식 또는 그람 행렬식을 관절 변수 함수로 표현한다.
인수 분해를 통해 특이 조건을 식별한다.
각 인수의 영점 집합을 결정하여 특이 다양체를 기술한다.

33.2.11.2 수치적 절차

현재 관절 구성의 자코비안을 계산한다.
SVD 또는 QR 분해를 통해 특이값을 평가한다.
최소 특이값, 조건수, 매니퓰러빌리티 지수를 계산한다.
임계치 초과 여부로 근접 결손 상태를 판정한다.

33.2.11.3 기하학적 절차

관절 트위스트를 스크류 이론으로 구성한다.
트위스트들이 생성하는 스크류 시스템의 차수를 평가한다.
상호 스크류를 식별하여 손실 방향의 기하학적 의미를 해석한다.

33.2.12 본 절의 학술적 정리

본 절에서 다룬 자코비안 계수 결손과 특이점의 관계는 특이점 이론의 선형 대수적 기초이다. 계수 결손은 행렬식 영점, 특이값 영점, 매니퓰러빌리티 영점, 스크류 시스템 차수 저하 등 다양한 등가 조건으로 기술되며, 결손의 차수는 특이점의 심각성을 정량화한다. 기구학적으로는 공선 정렬, 축 교차, 축 중첩이 결손의 대표적 원인이며, 정역학적으로는 속도 사상과 힘 사상에서 동시에 나타나는 쌍대 현상이다. 수치적 판별에서는 특이값 임계치, 조건수, 매니퓰러빌리티 지수가 활용되며, 실무적으로는 경로 재계획, 자체 운동, 구조 설계, 강건 역기구학을 통해 결손의 영향이 관리된다. 본 절의 분석은 후속 절에서 다룰 특이점 판별, 분류, 회피의 수학적 전제를 제공한다.

출처

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.
Siciliano, B., Sciavicco, L., Villani, L., and Oriolo, G., Robotics: Modelling, Planning and Control, Springer, 2009.
Golub, G. H. and Van Loan, C. F., Matrix Computations, 4th edition, Johns Hopkins University Press, 2013.
Strang, G., Introduction to Linear Algebra, 5th edition, Wellesley-Cambridge Press, 2016.
Lynch, K. M. and Park, F. C., Modern Robotics: Mechanics, Planning, and Control, Cambridge University Press, 2017.
Spong, M. W., Hutchinson, S., and Vidyasagar, M., Robot Modeling and Control, 2nd edition, Wiley, 2020.

버전

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