17.13 관성 행렬의 대칭성과 양정치성

1. 개요

매니퓰레이터 동역학 방정식에서 관성 행렬 \mathbf{M}(\mathbf{q})는 단순한 수치적 계수가 아니라, 시스템의 운동 에너지를 유도하는 이차 형식의 계수이며, 물리 법칙으로부터 유도되는 두 가지 핵심 대수적 성질을 지닌다. 첫째, 관성 행렬은 대칭(symmetric)이며, 둘째, 양의 정부호(positive definite)이다. 이 두 성질은 단순한 수학적 특성이 아니라, 운동 에너지의 물리적 의미와 직접적으로 연결되며, 수치 해석의 안정성, 제어기 설계의 유효성, 그리고 안정성 증명의 근간이 되는 조건이다.

본 절에서는 관성 행렬의 대칭성과 양정치성을 물리적, 수학적, 그리고 수치적 관점에서 학술적으로 기술한다. 또한 이 성질들이 함의하는 주요 결과와 수치 구현에서 이들을 유지하기 위한 실용적 고려를 다룬다.

2. 대칭성의 정의와 물리적 기원

대칭성이란 관성 행렬의 전치가 그 자신과 같다는 성질, 즉 \mathbf{M}(\mathbf{q}) = \mathbf{M}^T(\mathbf{q})임을 의미한다. 원소별로는 M_{ij}(\mathbf{q}) = M_{ji}(\mathbf{q})가 모든 i, j와 모든 자세 \mathbf{q}에 대하여 성립한다. 이 성질은 다음의 물리적 근거로부터 유도된다.

첫째, 운동 에너지는 관절 속도의 실수 값 이차 형식이며, 임의의 실수 이차 형식은 대칭 행렬로 유일하게 표현될 수 있다. 구체적으로 K = \frac{1}{2} \dot{\mathbf{q}}^T \mathbf{M}(\mathbf{q}) \dot{\mathbf{q}}에서 \mathbf{M}의 비대칭 부분은 이차 형식에 기여하지 않으며, 따라서 물리적 의미를 가지는 부분은 대칭 성분에 한정된다.

둘째, 질량과 관성 텐서의 자연스러운 정의로부터 유도되는 관성 행렬은 본래적으로 대칭 형식을 갖는다. 각 링크의 관성 텐서 \mathbf{I}_{c,i}는 대칭 행렬이며, 자코비안을 통한 변환에서 대칭성이 보존된다. 관성 행렬의 표준 구성 식

\mathbf{M}(\mathbf{q}) = \sum_{i=1}^{n} \left[ m_i \, \mathbf{J}_{v,i}^T \mathbf{J}_{v,i} + \mathbf{J}_{\omega,i}^T \mathbf{I}_{c,i}^{(0)} \mathbf{J}_{\omega,i} \right]

에서 각 항은 \mathbf{A}^T \mathbf{B} \mathbf{A}의 형식이며, \mathbf{B}가 대칭일 때 이 항 자체가 대칭이다. 따라서 각 링크의 기여가 대칭 행렬이 되며, 이들의 합인 \mathbf{M}(\mathbf{q})도 대칭이다.

셋째, 라그랑주 역학의 관점에서, 운동 에너지의 헤시안(Hessian)으로 정의되는 관성 행렬은 정확히 \partial^2 K / \partial \dot{q}_i \partial \dot{q}_j로 표현되며, 편미분의 순서 교환 가능성(Young의 정리)에 의하여 자동으로 대칭이 된다.

이러한 다중 근거는 대칭성이 우연적 특성이 아니라 본질적 물리 성질임을 나타낸다.

양정치성의 정의와 물리적 기원

양의 정부호성이란 관성 행렬이 대칭이면서, 임의의 영이 아닌 벡터 \mathbf{v} \in \mathbb{R}^n에 대하여 다음을 만족함을 의미한다.

\mathbf{v}^T \mathbf{M}(\mathbf{q}) \mathbf{v} > 0

이 성질은 다음의 물리적 근거로부터 유도된다. 시스템의 운동 에너지는 속도가 영이 아닌 어떠한 상태에서도 음이 될 수 없으며, 속도가 영이 아닌 한 양이어야 한다. 이는 에너지 개념 자체의 정의에서 나온 것이며, 질량이 양이고 관성 텐서가 양의 정부호라는 사실의 자연스러운 귀결이다. 질량은 물리적 의미에서 항상 양이며, 관성 텐서 \mathbf{I}_c는 연속적 질량 분포의 이차 모멘트로 정의되어 본래 양의 정부호 성질을 가진다.

수식적으로 운동 에너지 식 K = \frac{1}{2} \dot{\mathbf{q}}^T \mathbf{M}(\mathbf{q}) \dot{\mathbf{q}}에서 \dot{\mathbf{q}} \neq \mathbf{0}일 때 반드시 K > 0임을 요구하면, 이는 \mathbf{M}의 양정치성과 동등하다. 예외적으로 K = 0이 되는 상태는 모든 링크가 동시에 정지한 상태, 즉 \dot{\mathbf{q}} = \mathbf{0}인 경우에 한정된다.

3. 양정치성의 수학적 특성화

대칭 행렬 \mathbf{M}이 양의 정부호일 필요충분조건은 다음 중 어느 하나의 조건과 동등하다. 첫째, 모든 고유값이 양이다. 둘째, 행렬의 모든 주부분 행렬식(principal minor)이 양이다(실베스터 기준). 셋째, 행렬이 유일한 대칭 양의 정부호 제곱근 \mathbf{M}^{1/2}을 가진다. 넷째, 촐레스키 분해가 존재하여 \mathbf{M} = \mathbf{L} \mathbf{L}^T로 표현되며, 여기서 \mathbf{L}은 양의 대각 원소를 갖는 하삼각 행렬이다.

이러한 특성화는 관성 행렬이 양의 정부호임을 수치적으로 검증하거나, 양정치성을 명시적으로 활용하는 알고리즘을 구성할 때 각각 다른 관점에서 유용하다. 특히 촐레스키 분해는 선형 시스템 \mathbf{M} \mathbf{x} = \mathbf{b}의 풀이에서 일반 가우스 소거법보다 효율적이며, 수치적으로 안정적이다.

4. 역행렬과 제곱근의 존재

관성 행렬이 양의 정부호이므로 그 역행렬 \mathbf{M}^{-1}(\mathbf{q})이 모든 자세에서 존재하며, 역시 대칭이고 양의 정부호이다. 이는 순동역학 계산에서 가속도를 유일하게 결정할 수 있게 한다.

\ddot{\mathbf{q}} = \mathbf{M}^{-1}(\mathbf{q}) \left[ \boldsymbol{\tau} + \mathbf{J}^T \mathbf{F}_e - \mathbf{C}(\mathbf{q}, \dot{\mathbf{q}}) \dot{\mathbf{q}} - \mathbf{g}(\mathbf{q}) - \mathbf{f}(\dot{\mathbf{q}}) \right]

역행렬의 존재가 보장되지 않는다면 순동역학 해가 정의되지 않거나 특이 상황이 발생할 수 있으며, 관성 행렬의 양정치성은 이러한 문제가 결코 발생하지 않음을 보증한다.

또한 양의 정부호 대칭 행렬은 유일한 대칭 제곱근을 가지므로, \mathbf{M}^{1/2}을 이용한 변환을 통해 동역학 방정식을 단위 관성 공간에서 분석하거나, 가중 노름을 정의하는 데 사용된다. 예를 들어 가속도의 가중 최소화 문제에서 \mathbf{M}-가중 노름은 자연스러운 물리적 해석(운동 에너지)을 가진다.

자세 의존적 대칭성과 양정치성

관성 행렬은 자세 \mathbf{q}의 함수이지만, 대칭성과 양정치성은 모든 자세에서 균일하게 성립한다. 즉 관성 행렬의 이 두 성질은 작업 공간 전역에 걸쳐 유지되는 구조적 불변성이다. 이는 물리 법칙이 자세에 무관하게 성립한다는 사실의 반영이다.

구체적으로 양정치성이 모든 자세에서 성립함을 보이려면, 이차 형식 \mathbf{v}^T \mathbf{M}(\mathbf{q}) \mathbf{v}가 모든 영이 아닌 \mathbf{v}와 모든 \mathbf{q}에 대하여 양임을 확인하면 된다. 운동 에너지 표현을 이용하면 이는 자명하게 따라오며, 추가적 조건이 필요하지 않다.

고유값과 조건수의 자세 의존성

관성 행렬의 대칭성과 양정치성이 모든 자세에서 성립하지만, 구체적 고유값은 자세에 따라 크게 변한다. 가장 큰 고유값 \lambda_{\max}(\mathbf{M})과 가장 작은 고유값 \lambda_{\min}(\mathbf{M})의 비율인 조건수

\kappa(\mathbf{M}) = \frac{\lambda_{\max}(\mathbf{M}(\mathbf{q}))}{\lambda_{\min}(\mathbf{M}(\mathbf{q}))}

는 자세에 의존하는 값이며, 특정 자세에서 매우 커질 수 있다. 조건수가 큰 자세에서는 역행렬 계산의 수치 오차가 증폭되며, 관성적 결합이 강한 방향과 약한 방향 사이의 스케일 차이가 크다는 것을 의미한다.

조건수의 최솟값은 항상 1 이상이며, 이는 관성 행렬이 직교 행렬일 때(즉 모든 고유값이 같을 때) 달성된다. 실제 매니퓰레이터에서는 링크 질량 분포의 비대칭과 다양한 결합으로 인하여 1보다 상당히 큰 조건수를 가진다. 조건수가 크면 수치적 어려움이 발생하며, 이를 완화하기 위한 스케일링 및 정칙화 기법이 사용된다.

5. 에너지 부등식과 고유값 경계

양의 정부호성은 이차 형식에 대한 다음의 부등식으로 표현될 수 있다.

\lambda_{\min}(\mathbf{M}(\mathbf{q})) \|\dot{\mathbf{q}}\|^2 \leq \dot{\mathbf{q}}^T \mathbf{M}(\mathbf{q}) \dot{\mathbf{q}} = 2 K \leq \lambda_{\max}(\mathbf{M}(\mathbf{q})) \|\dot{\mathbf{q}}\|^2

만약 매니퓰레이터의 전 작업 공간에서 고유값이 일정한 경계 0 < \mu_1 \leq \lambda_{\min}(\mathbf{M}) \leq \lambda_{\max}(\mathbf{M}) \leq \mu_2 내부에 머무른다면, 운동 에너지는 관절 속도 크기에 의하여 양쪽에서 균일하게 경계된다. 이러한 균일 경계는 리아푸노프 안정성 증명에서 중요한 조건이며, 많은 제어기의 수학적 정당화에 필수적이다.

실제 매니퓰레이터의 경우, 관절 제한(joint limit) 내부에서 이러한 경계는 일반적으로 성립한다. 그러나 관절이 회전 제한을 갖지 않는 경우에도 질량과 관성의 유한성으로 인하여 경계는 여전히 존재한다.

대칭성과 양정치성이 보장하는 제어 성능

대칭성과 양정치성은 다양한 제어 기법의 안정성 증명에 핵심적으로 활용된다. 대표적 결과는 다음과 같다.

첫째, 중력 보상이 적용된 PD 제어기는 임의의 초기 조건에서 목표 자세로 점근적으로 수렴한다. 이 증명에서 관성 행렬의 양정치성은 운동 에너지 기반 리아푸노프 함수의 양의 정부호성을 보장하며, 시스템의 총 에너지가 시간에 따라 감소함을 보이는 데 사용된다.

둘째, 계산 토크 제어(computed torque control)는 피드백 선형화를 통하여 관성 행렬의 역수를 이용한다. 양정치성은 역행렬이 항상 존재함을 보장하여 피드백 선형화가 전 작업 공간에서 유효함을 담보한다.

셋째, 적응 제어는 관성 행렬의 선형 매개변수 형식과 양정치성을 함께 활용하여 수렴 증명을 수행한다.

넷째, 수동성 기반 제어(passivity-based control)는 운동 에너지를 저장 함수(storage function)로 채택하며, 양정치성이 이 함수의 적법성을 보장한다.

수치 구현에서의 대칭성과 양정치성 유지

이론적으로는 항상 성립하는 대칭성과 양정치성이 수치 계산에서는 반올림 오차에 의하여 미세하게 위배될 수 있다. 이러한 위배를 막기 위하여 다음의 실용적 고려가 적용된다.

첫째, 관성 행렬을 구성할 때 (\mathbf{M} + \mathbf{M}^T) / 2 형식의 명시적 대칭화 단계를 포함시킨다. 이는 비대칭 수치 오차를 제거한다.

둘째, 관성 행렬의 고유값이 예상 범위를 크게 벗어나는 경우 스케일링이나 단위 변환을 재검토한다. 자세에 따라 고유값이 급격히 변할 때에는 자세 의존적 스케일링을 적용하여 수치적 조건을 개선한다.

셋째, 선형 시스템 \mathbf{M} \mathbf{x} = \mathbf{b}의 풀이에는 촐레스키 분해를 우선적으로 사용한다. 촐레스키 분해는 양의 정부호성이 유지되어야만 실패 없이 진행되며, 계산 효율과 수치 안정성에서 가우스 소거법보다 우수하다.

넷째, 수치 시뮬레이션의 장시간 적분에서 적분 오차에 의한 에너지 드리프트가 발생할 수 있다. 이 경우 구조 보존형 적분기(예: 심플렉틱 적분기, 변분 적분기)를 사용하여 에너지 보존 성질을 유지하는 것이 유리하다.

다섯째, 관성 매개변수 식별 과정에서 식별된 매개변수가 물리적으로 일관되지 않으면 관성 행렬이 양의 정부호를 잃을 수 있다. 이를 방지하기 위하여 물리적 일관성 제약(physical consistency constraints)을 식별 최적화 문제에 추가하는 기법이 제안되어 있다.

물리적 일관성과 식별 가능한 매개변수

관성 행렬이 양의 정부호라는 사실은 기초 물리 법칙에서 유도되지만, 매개변수 식별 과정에서 통계적 방법만을 사용하면 식별된 관성 매개변수가 이 조건을 위배하는 비물리적 값을 가질 수 있다. 이러한 경우 식별 결과로 구성된 관성 행렬이 일부 자세에서 양의 정부호를 상실할 수 있다.

이를 방지하기 위하여 식별 과정에 각 링크의 관성 텐서가 양의 정부호이고 삼각 부등식(triangle inequality)을 만족해야 한다는 제약이 추가된다. 구체적으로 각 링크의 관성 매개변수가 물리적으로 실현 가능한 질량 분포에 대응하는지를 검사하는 알고리즘이 개발되어 있으며, 식별 결과의 물리적 일관성을 보장한다.

본 절의 의의

본 절에서 기술한 관성 행렬의 대칭성과 양정치성은 매니퓰레이터 동역학의 가장 기본적인 대수적 성질이며, 이 성질들이 없으면 대부분의 이론적 분석과 수치적 구현이 붕괴된다. 대칭성은 운동 에너지의 이차 형식 구조에서 자연스럽게 나오며, 양정치성은 에너지가 양이라는 기본 물리 법칙의 반영이다.

이 두 성질은 제어기 설계, 안정성 증명, 수치 시뮬레이션, 그리고 매개변수 식별의 전 과정에 걸쳐 반복적으로 활용되며, 본 장의 후속 절들에서 다루어지는 다양한 주제의 수학적 기반을 이룬다.

학습 권장사항

본 절의 내용을 충분히 이해하기 위하여 대칭 행렬과 이차 형식, 고유값 분해, 촐레스키 분해, 그리고 리아푸노프 안정성 이론에 대한 선행 학습이 권장된다. 본 절의 내용을 심화 학습하기 위해서는 수치 선형 대수, 구조 보존형 적분 기법, 그리고 물리적으로 일관된 매개변수 식별 이론에 대한 추가 학습이 권장된다.

참고 문헌

• Spong, M. W., Hutchinson, S., and Vidyasagar, M. (2020). Robot Modeling and Control (2nd ed.). Wiley.
• Siciliano, B., Sciavicco, L., Villani, L., and Oriolo, G. (2009). Robotics: Modelling, Planning and Control. Springer.
• Featherstone, R. (2008). Rigid Body Dynamics Algorithms. Springer.
• Murray, R. M., Li, Z., and Sastry, S. S. (1994). A Mathematical Introduction to Robotic Manipulation. CRC Press.
• Golub, G. H., and Van Loan, C. F. (2013). Matrix Computations (4th ed.). Johns Hopkins University Press.
• Khalil, H. K. (2015). Nonlinear Systems (3rd ed.). Pearson.
• Traversaro, S., Brossette, S., Escande, A., and Nori, F. (2016). Identification of Fully Physical Consistent Inertial Parameters Using Optimization on Manifolds. In Proceedings of the IEEE/RSJ International Conference on Intelligent Robots and Systems.

version: 1.0