32.40 가중 의사 역행렬(Weighted Pseudoinverse)

32.40 가중 의사 역행렬(Weighted Pseudoinverse)

가중 의사 역행렬은 Moore-Penrose 의사 역행렬을 일반화한 연산자로, 관절 공간과 작업 공간에 각각 대칭 양정 가중 행렬을 도입하여 해의 노름 척도와 잔여 척도를 작업의 물리적 의미에 맞게 변형한다. 가중 의사 역행렬은 관절 가속 에너지 최소화, 동적 일관성, 관절 한계 회피, 단위 이질성의 정규화, 액추에이터 토크 가중 등 다양한 제어 목적에서 활용되며, 자코비안 기반 속도 사상의 단일 형태 위에 작업의 의도를 직접 반영하는 구조적 자유도를 제공한다. 본 절에서는 가중 의사 역행렬의 정의와 그에 대응하는 최적화 문제, 대표적 가중 선택과 응용, 영공간 투영 성질, 그리고 수치적 구현과 한계를 차분하게 기술한다.

1. 정의와 가중 노름 기반 최적화 해석

행렬 \mathbf{J} \in \mathbb{R}^{m \times n}과 두 대칭 양정 가중 행렬 \mathbf{W}_q \in \mathbb{R}^{n \times n}, \mathbf{W}_x \in \mathbb{R}^{m \times m}를 도입한다. 자코비안의 행 계수가 완전한 경우, 자코비안의 가중 의사 역행렬 \mathbf{J}^\#_{\mathbf{W}_q, \mathbf{W}_x}는 다음과 같이 정의된다. \mathbf{J}^\#_{\mathbf{W}_q, \mathbf{W}_x} = \mathbf{W}_q^{-1}\,\mathbf{J}^\top \bigl( \mathbf{J}\,\mathbf{W}_q^{-1}\,\mathbf{J}^\top \bigr)^{-1}로 표현되며, 이는 작업 공간 가중치 \mathbf{W}_x가 단위 행렬인 경우의 표준 형태이다. 일반적인 작업 공간 가중치를 포함하는 형태는 잔여 항에 가중을 부여한 최소 제곱 문제의 해로서 정의되며, 자유도 부족, 비여유, 여유 자유도의 차원 관계에 따라 그 의미가 달라진다.

가중 의사 역행렬의 최적화 해석은 다음의 가중 노름 최적화 문제로 명료하게 정리된다. 자유도 부족, 즉 m > n인 경우 \dot{\vec{q}}^\ast = \arg\min_{\dot{\vec{q}}}\, (\dot{\vec{x}} - \mathbf{J}\,\dot{\vec{q}})^\top\, \mathbf{W}_x\, (\dot{\vec{x}} - \mathbf{J}\,\dot{\vec{q}})의 가중 잔여 최소화 해이며, 가중 정규 방정식 \mathbf{J}^\top\,\mathbf{W}_x\,\mathbf{J}\,\dot{\vec{q}} = \mathbf{J}^\top\,\mathbf{W}_x\,\dot{\vec{x}}의 해로 표현된다. 여유 자유도, 즉 n > m인 경우 \dot{\vec{q}}^\ast = \arg\min_{\dot{\vec{q}}}\, \dot{\vec{q}}^\top\,\mathbf{W}_q\,\dot{\vec{q}}\quad \text{s.t.}\quad \mathbf{J}\,\dot{\vec{q}} = \dot{\vec{x}}의 가중 노름 최소화 해이며, 라그랑주 승수법으로부터 \dot{\vec{q}}^\ast = \mathbf{W}_q^{-1}\,\mathbf{J}^\top\,(\mathbf{J}\,\mathbf{W}_q^{-1}\,\mathbf{J}^\top)^{-1}\,\dot{\vec{x}}가 유도된다. 일반적 형태의 가중 의사 역행렬은 두 최적화의 결합인 가중 최소 노름 최소 제곱 문제의 해로 단일하게 정리된다.

2. 대표적 가중 선택과 물리적 의미

관절 공간 가중 행렬 \mathbf{W}_q의 가장 단순한 선택은 양의 대각 행렬 \mathbf{W}_q = \operatorname{diag}(w_1, w_2, \dots, w_n)이며, 각 w_i는 해당 관절 변수의 가중치를 의미한다. 관절 한계 회피를 위해 Chan과 Dubey의 1995년 논문 “A weighted least-norm solution based scheme for avoiding joint limits for redundant joint manipulators“가 제시한 가중 함수 w_i(q_i) = 1 + |\partial H/\partial q_i|는 관절이 한계에 가까워질수록 가중치가 증가하여 해당 관절의 사용을 자동으로 억제한다. 여기서 H(\vec{q})는 관절 한계까지의 거리를 표현하는 비용 함수이며, 그 부분 도함수의 절댓값이 가중치에 더해진다. 이 접근은 영공간 투영 기반 부차 목표 처리와 달리 비여유 매니퓰레이터에도 직접 적용된다는 이점을 가진다.

관성 행렬 \mathbf{M}(\vec{q})를 가중 행렬로 사용하는 동적 일관 의사 역행렬 \mathbf{J}^\sharp_{\mathbf{M}} = \mathbf{M}^{-1}\,\mathbf{J}^\top (\mathbf{J}\,\mathbf{M}^{-1}\,\mathbf{J}^\top)^{-1}은 Khatib의 1987년 논문 “A unified approach for motion and force control of robot manipulators: The operational space formulation“에서 제시한 작업 공간 동역학 정식화의 핵심 연산자이다. 이 의사 역행렬은 관절 가속 에너지 \dot{\vec{q}}^\top \mathbf{M} \dot{\vec{q}}를 최소화하는 해를 산출하며, 영공간 운동이 작업 공간 가속에 동적으로 결합되지 않는 동적 일관성을 보장한다. 작업 공간 가중 행렬 \mathbf{W}_x의 활용은 선속도와 각속도의 단위 이질성을 정규화하거나, 특정 작업 방향에 더 큰 정밀도 가중을 부여하는 등의 목적에 사용된다. 특성 길이 기반 정규화에서 도입한 가중치는 자코비안의 등방성 평가와 직접 연결된다.

3. 영공간 투영 성질과 일반 해의 구조

가중 의사 역행렬의 영공간 투영자는 가중 행렬에 따라 변형된 형태로 정의된다. 표준 영공간 투영자 \mathbf{N} = \mathbf{I} - \mathbf{J}^+\,\mathbf{J}의 일반화로서, 가중 의사 역행렬에 대응하는 투영자는 \mathbf{N}_{\mathbf{W}_q} = \mathbf{I} - \mathbf{J}^\#_{\mathbf{W}_q}\,\mathbf{J}이며, 이 연산자는 작업 공간 명령에 영향을 미치지 않는 자기 운동 부분 공간을 가중 노름의 직교 보충으로 사영한다. 따라서 일반 해는 \dot{\vec{q}} = \mathbf{J}^\#_{\mathbf{W}_q}\,\dot{\vec{x}} + \mathbf{N}_{\mathbf{W}_q}\,\dot{\vec{q}}_0의 형태로 표현되며, 부차 목표 \dot{\vec{q}}_0는 가중 노름의 의미에서 영공간 부분에 사영된다. 표준 의사 역행렬 해와 달리 가중 의사 역행렬 해는 일반적으로 직교 사영이 아니라 가중 노름 의미의 사영이라는 점에서 구별되며, 이 차이는 작업의 물리적 의미를 반영하는 핵심 자유도가 된다.

가중 의사 역행렬 해의 해석은 여유 자유도 매니퓰레이터의 우선순위 기반 다중 작업 처리에서 특히 중요하다. 주작업의 자코비안에 동적 일관 의사 역행렬을 적용하면 작업 공간 가속 명령이 관절 가속 에너지를 최소화하는 형태로 분해되며, 부작업은 가중 영공간 투영자를 통해 결합된다. 이 다층 구조는 인간형 로봇과 같이 다수의 작업을 동시에 만족시켜야 하는 시스템에서 표준적으로 사용되며, Sentis와 Khatib의 2005년 논문 “Synthesis of whole-body behaviors through hierarchical control of behavioral primitives“는 이러한 가중 우선순위 처리를 일반화된 행동 기반 제어 틀로 확장하였다. 가중 의사 역행렬은 따라서 단일한 자코비안 사상의 형식 안에서 다양한 작업 의도와 동역학적 일관성을 동시에 표현하는 통합적 도구로 자리매김한다.

4. 수치적 구현, 특이 근방 거동, 그리고 한계

가중 의사 역행렬의 수치적 구현은 일반적으로 변환된 자코비안의 표준 의사 역행렬 계산으로 환원된다. \mathbf{W}_q가 양정이면 콜레스키 분해 \mathbf{W}_q = \mathbf{L}_q\,\mathbf{L}_q^\top를 도입하여 변환된 자코비안 \tilde{\mathbf{J}} = \mathbf{J}\,\mathbf{L}_q^{-\top}을 정의하면, 가중 의사 역행렬은 \mathbf{J}^\#_{\mathbf{W}_q} = \mathbf{L}_q^{-\top}\,\tilde{\mathbf{J}}^+로 표현된다. 이 변환을 통해 가중 의사 역행렬 계산은 표준 의사 역행렬 계산 알고리즘에 대한 사전·사후 변환 단계의 결합으로 환원되며, 특이값 분해, QR 분해, LU 분해 등 표준 도구를 그대로 활용할 수 있다. 작업 공간 가중 행렬 \mathbf{W}_x의 도입도 동일한 방식으로 잔여 항을 변환하여 처리된다.

가중 의사 역행렬의 특이 근방 거동은 표준 의사 역행렬과 본질적으로 동일한 한계를 가진다. 즉 자코비안의 최소 특이값이 영에 근접하면 가중 의사 역행렬의 노름도 발산하며, 가중 인자가 이 발산을 제한적으로 완화할 수는 있어도 본질적으로 제거하지는 못한다. 따라서 특이 근방의 수치 안정성을 위해서는 가중 의사 역행렬에 정칙화 항을 결합한 형태, 즉 가중 감쇠 최소 제곱 의사 역행렬 \mathbf{J}^\sharp_{\mathbf{W}_q,\lambda} = \mathbf{W}_q^{-1}\,\mathbf{J}^\top\,(\mathbf{J}\,\mathbf{W}_q^{-1}\,\mathbf{J}^\top + \lambda^2 \mathbf{I})^{-1}이 사용된다. 이 결합 형태는 다음 절에서 다루어질 감쇠 최소 제곱 자코비안의 가중 일반화에 해당하며, 가중 의사 역행렬과 정칙화의 결합이 매니퓰레이터 제어에서 표준적인 자코비안 처리 도구임을 보여 준다. 가중 행렬의 시간 변동, 가중 함수의 미분 가능성, 부차 목표와의 상호 작용 등은 실시간 제어 안정성과 직접 연관되므로, 설계 시 주의 깊은 평가가 요구된다.

5. 본 절의 정리

가중 의사 역행렬은 표준 Moore-Penrose 의사 역행렬을 가중 노름 기반 최적화의 해로 일반화한 연산자이며, 관절 공간 가중 \mathbf{W}_q와 작업 공간 가중 \mathbf{W}_x의 도입을 통해 작업의 물리적 의미를 자코비안 사상에 직접 반영한다. 관절 한계 회피, 동적 일관성, 단위 이질성 정규화, 액추에이터 토크 가중과 같은 다양한 응용은 모두 가중 행렬의 적절한 선택으로 단일 형식 안에서 통합된다. 가중 영공간 투영자를 통한 일반 해의 구조는 여유 자유도 매니퓰레이터의 다중 작업 처리에서 핵심적이며, 우선순위 기반 처리의 표준 도구로 활용된다. 그러나 자코비안의 특이 근방에서의 노름 발산이라는 본질적 한계는 가중 의사 역행렬에서도 그대로 유지되므로, 후속 절에서 다루어질 감쇠 최소 제곱법과의 결합이 표준적 실무 형식으로 자리잡고 있다. 본 절은 자코비안 기반 사상 연산자의 일반화와 그 응용 다양성을 정리하면서, 다음 단계의 정칙화 기반 변형으로 이어지는 가교를 형성한다.

6. 출처

  • Penrose, R., “A generalized inverse for matrices”, Mathematical Proceedings of the Cambridge Philosophical Society, Vol. 51, No. 3, pp. 406–413, 1955.
  • 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.
  • Khatib, O., “A unified approach for motion and force control of robot manipulators: The operational space formulation”, IEEE Journal of Robotics and Automation, Vol. 3, No. 1, pp. 43–53, 1987.
  • Chan, T. F. and Dubey, R. V., “A weighted least-norm solution based scheme for avoiding joint limits for redundant joint manipulators”, IEEE Transactions on Robotics and Automation, Vol. 11, No. 2, pp. 286–292, 1995.
  • Sentis, L. and Khatib, O., “Synthesis of whole-body behaviors through hierarchical control of behavioral primitives”, International Journal of Humanoid Robotics, Vol. 2, No. 4, pp. 505–518, 2005.
  • Sciavicco, L. and Siciliano, B., Modelling and Control of Robot Manipulators, 2nd edition, Springer, 2000.
  • Ben-Israel, A. and Greville, T. N. E., Generalized Inverses: Theory and Applications, 2nd edition, Springer, 2003.

7. 버전

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