32.21 자코비안의 영공간(Null Space)의 기하학적 의미

자코비안의 영 공간(null space)은 엔드 이펙터의 순간 속도에 기여하지 않는 관절 속도 벡터의 집합이며, 여유 자유도 매니퓰레이터의 자기 운동(self-motion)을 수학적으로 기술하는 핵심 부분 공간이다. 영 공간의 기하학적 해석은 매니퓰레이터가 엔드 이펙터의 자세를 유지하면서 내부 관절 구성을 변화시키는 운동의 가능성을 정량화하며, 부차적 과제 수행, 장애물 회피, 특이점 회피, 매니퓰러빌리티 최적화 등의 다양한 응용에 직접 활용된다. 본 절에서는 영 공간의 학술적 정의, 기하학적 의미, 계수-영차원 정리를 통한 차원 관계, 특이점에서의 영 공간 확장, 영 공간 투영 연산자, 그리고 자기 운동의 응용을 체계적으로 다룬다.

1. 영 공간의 학술적 정의

자코비안 \mathbf{J}(\vec{q}) \in \mathbb{R}^{m \times n}의 영 공간은 자코비안이 영 벡터로 사상하는 관절 속도 벡터의 집합으로 정의된다.

\mathcal{N}(\mathbf{J}(\vec{q})) = \{\dot{\vec{q}} \in \mathbb{R}^n : \mathbf{J}(\vec{q}) \, \dot{\vec{q}} = \vec{0}\}

영 공간은 \mathbb{R}^n의 선형 부분 공간이며, 그 차원은 계수-영차원 정리(rank-nullity theorem)에 따라 다음과 같이 결정된다.

\dim \mathcal{N}(\mathbf{J}(\vec{q})) = n - \text{rank}(\mathbf{J}(\vec{q}))

영 공간의 원소인 관절 속도 벡터를 영 공간 속도(null-space velocity) 또는 자기 운동 속도(self-motion velocity)라 한다.

2. 영 공간의 기하학적 해석

영 공간의 관절 속도 벡터는 다음의 기하학적 의미를 가진다. 매니퓰레이터가 어떤 구성 \vec{q}_0에서 출발하여 영 공간 속도 \dot{\vec{q}} \in \mathcal{N}(\mathbf{J}(\vec{q}_0))로 순간 운동을 시작하면, 엔드 이펙터의 순간 속도는 영 벡터이므로 엔드 이펙터는 순간적으로 정지 상태를 유지한다. 반면 매니퓰레이터의 관절들은 비영 속도로 움직이므로 내부 구성은 변화한다.

이러한 운동을 자기 운동이라 하며, 엔드 이펙터의 자세(위치와 방향)를 보존하면서 매니퓰레이터의 내부 기하학적 형태를 변경하는 운동이다. 여유 자유도 매니퓰레이터의 경우 동일 엔드 이펙터 자세에 대응하는 관절 구성이 무한히 존재하며, 이러한 구성들이 형성하는 부분 다양체(submanifold)의 접공간이 영 공간이 된다.

3. 계수와 영 공간 차원의 관계

자코비안의 계수와 영 공간의 차원은 다음의 직접적 관계를 가진다.

특이점에서의 영 공간 차원 증가는 기구학적 특이점의 본질적 특징이다. 정방 매니퓰레이터의 경우 완전 계수에서 영 공간 차원이 0이지만, 특이점에 진입하면 계수가 n-1로 감소하여 영 공간 차원이 1로 증가한다. 여유 자유도 매니퓰레이터의 경우 정규 구성에서 영 공간 차원이 n-m이었다면 특이점에서 n-m+1로 증가한다.

4. 영 공간의 기저 구성

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

4.1 SVD 기반 구성

자코비안의 특이값 분해 \mathbf{J} = \mathbf{U} \mathbf{\Sigma} \mathbf{V}^\top에서 영 특이값에 대응하는 우특이 벡터들이 영 공간의 정규 직교 기저를 제공한다. \mathbf{J}의 계수가 r이면 \mathbf{V}의 뒤 n - r개 열 \vec{v}_{r+1}, \ldots, \vec{v}_n이 \mathcal{N}(\mathbf{J})의 기저가 된다.

\mathcal{N}(\mathbf{J}) = \text{span}\{\vec{v}_{r+1}, \ldots, \vec{v}_n\}

이 기저는 정규 직교 기저이며 수치적으로 안정적으로 계산된다.

32.21.4.2 QR 분해 기반 구성

자코비안 전치의 QR 분해 \mathbf{J}^\top = \mathbf{Q} \mathbf{R}에서 \mathbf{R}의 영 행에 대응하는 \mathbf{Q}의 열이 영 공간의 기저를 제공한다. QR 기반 방법은 SVD보다 계산량이 적으나 특이점 근방에서의 수치 정확도는 낮을 수 있다.

32.21.4.3 영 공간 투영 연산자

명시적 기저 없이 영 공간으로의 직교 투영 연산자를 구성할 수 있다. 여유 자유도 매니퓰레이터의 경우 의사 역행렬을 이용한 투영 행렬은 다음과 같다.

\mathbf{P}_N(\vec{q}) = \mathbf{I}_n - \mathbf{J}^+(\vec{q}) \, \mathbf{J}(\vec{q})

여기서 \mathbf{J}^+ = \mathbf{J}^\top (\mathbf{J} \mathbf{J}^\top)^{-1}은 무어-펜로즈 오른 의사 역행렬이다. 투영 행렬 \mathbf{P}_N은 \mathbf{P}_N^2 = \mathbf{P}_N(멱등성)과 \mathbf{P}_N^\top = \mathbf{P}_N(대칭성)을 만족하는 직교 투영 연산자이다.

임의의 관절 속도 벡터 \dot{\vec{q}}_0에 대하여 \mathbf{P}_N \dot{\vec{q}}_0 \in \mathcal{N}(\mathbf{J})이 성립하며, 이는 \dot{\vec{q}}_0의 영 공간 성분을 추출한다.

5. 역속도 기구학에서의 영 공간

여유 자유도 매니퓰레이터의 역속도 기구학 일반 해는 최소 노름 해와 영 공간 성분의 합으로 표현된다.

\dot{\vec{q}} = \mathbf{J}^+(\vec{q}) \, \dot{\vec{x}}_d + \mathbf{P}_N(\vec{q}) \, \dot{\vec{q}}_0

여기서 \dot{\vec{x}}_d는 목표 엔드 이펙터 속도이고, \dot{\vec{q}}_0는 임의의 관절 속도 벡터이다. 첫 번째 항 \mathbf{J}^+ \dot{\vec{x}}_d는 목표 속도를 최소 관절 속도 노름으로 실현하는 해이며, 두 번째 항 \mathbf{P}_N \dot{\vec{q}}_0는 엔드 이펙터 속도에 기여하지 않는 자기 운동 성분이다.

이러한 분해는 Liégeois(1977)의 선구적 연구에서 체계화되었으며, 여유 자유도 매니퓰레이터의 부차적 과제 수행의 수학적 토대이다.

32.21.6 자기 운동 부분 다양체

여유 자유도 매니퓰레이터의 경우, 고정된 엔드 이펙터 자세 \vec{x}_0에 대응하는 관절 구성의 집합은 관절 공간 \mathbb{R}^n 내의 부분 다양체(self-motion manifold)를 형성한다.

\mathcal{M}_{\vec{x}_0} = \{\vec{q} \in \mathbb{R}^n : \vec{f}(\vec{q}) = \vec{x}_0\}

이 부분 다양체의 차원은 일반적으로 n - m이며, 매니퓰레이터의 여유 자유도 수와 같다. 부분 다양체의 한 점 \vec{q}_0에서의 접공간이 바로 \mathcal{N}(\mathbf{J}(\vec{q}_0))이다.

자기 운동 부분 다양체의 위상적 성질(연결성, 경계, 폐쇄성 등)은 매니퓰레이터의 기구학적 특성에 의해 결정된다. 연결 부분 다양체 내에서는 자기 운동을 통해 한 관절 구성에서 다른 구성으로 이동할 수 있으나, 부분 다양체가 여러 연결 성분으로 분리된 경우에는 다른 성분으로의 이동이 자기 운동만으로는 불가능하다. 이러한 성질은 7자유도 인간형 로봇 팔의 기구학 해석에서 중요하게 다루어진다.

6. 특이점 근방의 영 공간 거동

특이점 근방에서 영 공간은 다음의 특징적 거동을 보인다.

영 공간의 확장: 특이점에 접근함에 따라 최소 특이값이 영으로 수렴하고, 대응하는 우특이 벡터가 영 공간의 추가 기저 벡터가 된다. 즉 영 공간 차원이 불연속적으로 증가한다.

영 공간 기저의 회전: 특이점에 접근할 때 기존 영 공간 기저 벡터들도 회전하여 새로운 방향으로 재배치된다. 이러한 기저의 급격한 변화는 영 공간 투영 기반 제어의 불연속성을 초래할 수 있다.

가속도 수준의 폭발: 특이점에 가까워지면 의사 역행렬의 노름이 발산하며, 영 공간 투영 연산자의 시간 미분이 무한대로 발산할 수 있다. 이는 자기 운동 기반 2차 과제 제어의 수치적 안정성을 저해한다.

이러한 특이점 근방 거동의 완화를 위해 감쇠 최소 제곱법(damped least squares), 가중 의사 역행렬, 영 공간 최적화 기반 정칙화 등의 기법이 활용된다.

7. 영 공간을 활용한 부차적 과제 수행

영 공간의 자기 운동은 엔드 이펙터 주된 과제에 영향을 주지 않고 부차적 과제를 수행하는 데 활용된다. 부차적 과제의 전형적 예는 다음과 같다.

7.1 관절 한계 회피

관절 한계로부터의 거리를 최대화하는 부차적 과제. 관절 변수 q_i의 중앙값 q_{i,c}로부터의 편차를 최소화하는 비용 함수 H(\vec{q}) = \sum_i ((q_i - q_{i,c}) / \Delta q_i)^2를 정의하고, 그래디언트 \nabla H의 부호 반대 방향으로 자기 운동을 실행한다.

\dot{\vec{q}} = \mathbf{J}^+ \dot{\vec{x}}_d - \alpha \, \mathbf{P}_N \, \nabla H(\vec{q})

32.21.8.2 매니퓰러빌리티 최적화

매니퓰러빌리티 지수 w(\vec{q}) = \sqrt{\det(\mathbf{J} \mathbf{J}^\top)}를 최대화하는 자기 운동. 이는 특이점으로부터 멀어지는 방향으로의 자기 운동을 생성하며, 제어 성능의 유지에 기여한다.

32.21.8.3 장애물 회피

환경 내 장애물과의 거리를 최대화하는 부차적 과제. 매니퓰레이터의 특정 링크 또는 엔드 이펙터와 장애물 사이의 최소 거리를 목적 함수로 설정하고, 영 공간 성분을 이용하여 엔드 이펙터 경로에 영향을 주지 않으며 회피 운동을 수행한다.

32.21.8.4 에너지 최소화

관절 토크 노름이나 동역학적 에너지를 최소화하는 자기 운동. 이는 에너지 효율적 작업 수행에 유용하며, 특히 배터리 구동 이동 매니퓰레이터에서 중요하다.

32.21.9 여러 우선순위 과제의 영 공간 투영

복수의 부차적 과제가 우선순위에 따라 수행되는 계층적 제어에서 영 공간 투영이 반복적으로 적용된다. 주된 과제 1, 부차적 과제 2, 3차 과제 3의 우선순위 구조에서

\dot{\vec{q}} = \mathbf{J}_1^+ \dot{\vec{x}}_1 + \mathbf{P}_{N,1} \, (\mathbf{J}_2 \mathbf{P}_{N,1})^+ (\dot{\vec{x}}_2 - \mathbf{J}_2 \mathbf{J}_1^+ \dot{\vec{x}}_1) + \cdots

의 형태로 표현된다. 여기서 \mathbf{P}_{N,1}은 과제 1의 영 공간 투영 연산자이고, 연속적 영 공간 투영을 통해 상위 과제를 보존하면서 하위 과제를 수행한다. 이러한 우선순위 기반 영 공간 제어는 Nakamura, Hanafusa, Yoshikawa(1987)의 선구적 연구와 Siciliano, Slotine(1991)의 확장 연구에서 정립되었다.

8. 물리적 해석과 기구학적 함의

영 공간의 물리적 해석은 다음 두 가지 관점에서 통찰을 제공한다.

운동학적 관점: 영 공간은 매니퓰레이터가 엔드 이펙터를 고정한 채 ’재배치’할 수 있는 내부 운동의 공간이다. 인간 팔의 엔드 이펙터(손)를 고정한 채 팔꿈치를 움직이는 운동이 이에 해당한다.

에너지 관점: 영 공간 운동은 엔드 이펙터에 일을 하지 않는 운동이다. 따라서 외부 부하에 대한 일 관점에서는 불필요한 운동으로 보일 수 있으나, 내부 최적화(한계 회피, 장애물 회피 등)의 관점에서는 필수적이다.

리 군 기하학적 관점: 자기 운동 부분 다양체는 엔드 이펙터의 자세를 보존하는 관절 구성의 집합이며, 이는 SE(3)의 단위원 \mathbf{T}_0 \in SE(3)의 역상(inverse image) \vec{f}^{-1}(\mathbf{T}_0)으로 이해된다. 순기구학 사상 \vec{f}의 접사상의 커널로서 영 공간이 자기 운동 부분 다양체의 접공간을 형성한다.

9. 학술적 의의

본 절에서 정립한 자코비안의 영 공간의 기하학적 의미는 로봇 공학에서 다음의 학술적 의의를 가진다. 첫째, 영 공간을 자기 운동의 수학적 기술로 정립함으로써 여유 자유도 매니퓰레이터의 핵심 개념적 토대를 제공한다. 둘째, 계수-영차원 정리를 통한 차원 관계는 정규 구성과 특이점에서의 영 공간 거동을 통일적으로 해석하는 틀을 제공한다. 셋째, SVD 기반·QR 기반·투영 연산자 기반 영 공간 구성 방법은 응용에 따른 수치적 선택 지침을 제공한다. 넷째, 자기 운동 부분 다양체의 개념은 영 공간을 관절 공간의 기하학적 구조로 확장 해석한다. 다섯째, 영 공간 투영을 활용한 부차적 과제 수행의 수학적 틀은 여유 자유도 활용의 실무적 설계 도구이다. 여섯째, 우선순위 기반 영 공간 제어는 복수 과제 동시 수행의 체계적 방법이다. 일곱째, 본 절의 내용은 상공간 해석과 영 공간 투영의 응용, 그리고 매니퓰러빌리티와 특이점 분석의 후속 주제에 대한 공통 기초를 제공한다.

10. 출처

• Liégeois, A., “Automatic supervisory control of the configuration and behavior of multibody mechanisms”, IEEE Transactions on Systems, Man, and Cybernetics, Vol. 7, No. 12, pp. 868–871, 1977.
• Hollerbach, J. M. and Suh, K. C., “Redundancy resolution of manipulators through torque optimization”, IEEE Journal on Robotics and Automation, Vol. 3, No. 4, pp. 308–316, 1987.
• Nakamura, Y., Hanafusa, H., and Yoshikawa, T., “Task-priority based redundancy control of robot manipulators”, International Journal of Robotics Research, Vol. 6, No. 2, pp. 3–15, 1987.
• Siciliano, B. and Slotine, J.-J. E., “A general framework for managing multiple tasks in highly redundant robotic systems”, Proceedings of the 5th International Conference on Advanced Robotics, pp. 1211–1216, 1991.
• Nakamura, Y., Advanced Robotics: Redundancy and Optimization, Addison-Wesley, 1991.
• 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.
• Chiaverini, S., “Singularity-robust task-priority redundancy resolution for real-time kinematic control of robot manipulators”, IEEE Transactions on Robotics and Automation, Vol. 13, No. 3, pp. 398–410, 1997.
• 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.

11. 버전

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