디스크립터 시스템 기반 정밀 수치 해석과 소프트웨어 도구

본 논문은 디스크립터 시스템(Descriptor System)이라는 일반화된 상태‑공간 모델을 중심으로, 유리 행렬(Rational Matrix)의 신뢰성 있는 수치 연산을 위한 이론적 기반과 실용적인 소프트웨어 도구들을 포괄적으로 정리한다. 1. **디스크립터 시스템 정의 및 필요성** - 연속시간 시스템은 E·ẋ(t)=A·x(t)+B·u(t), y(t)=C·x(t)+D·u(t) 형태이며, 이산시간은 E·x(k+1)=A·x(k)+B·u(k) 로 표현된다. 여기서 E는 정방 행렬이며, 정칙이 아니어도 A−λE가 정규(pencil)라면 시스템은 완전하게 정의된다. - E=I인 표준 상태‑공간 모델은 디스크립터 시스템의 특수 경우이며, 복합적인 미분·대수 방정식, 제약 메카니즘, 경제 모델 등에서 자연스럽게 등장한다. 2. **유리 행렬과 실현 문제** - 전이 함수 행렬 G(λ)=C(A−λE)^{-1}B+D 은 복소 변수 λ에 대한 유리 행렬이다. G(λ)를 직접 다루면 다항식 기반 표현의 수치적 민감도가 크게 증가한다. 따라서 G(λ)를 구현하는 디스크립터 실현 (A−λE, B, C, D)을 찾는 것이 핵심이다. - 실현은 유일하지 않으며, U·V 가역 변환에 의해 동등한 실현들이 존재한다. 이러한 자유도는 최소 차원 실현을 찾는 데 활용된다. 3. **최소 실현의 정의와 판정** - 최소 실현은 상태 차원 n이 가능한 최소값인 경우를 말한다. Theorem 1에 따르면, 최소 실현은 (i) 전·무한 제어 가능성, (ii) 전·무한 관측 가능성, (iii) 비동적 모드 부재(A·N(E)⊆R(E))를 만족해야 한다. - 조건 (i)·(ii)는 각각 rank