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

초록

디스크립터 시스템 표현을 이용해 유리 행렬을 안정적으로 다루는 이론과, 최소 차원 실현, Weierstrass·Kronecker 정규형, 그리고 이를 구현한 최신 수치 소프트웨어 도구들을 종합적으로 소개한다.

상세 분석

본 논문은 선형 시불변(LTI) 시스템을 일반화한 디스크립터 시스템(정의식 E·ẋ = A x + B u, y = C x + D u) 형태를 핵심 모델링 프레임워크로 채택한다. E가 정칙이 아니더라도 A − λE가 정규(pencil)라면 시스템은 완전하게 정의되며, 전통적인 상태‑공간 모델(I = E)보다 더 넓은 클래스의 연속·이산 시간 시스템을 포괄한다.

논문은 특히 유리 행렬 G(λ)=C(A−λE)^{-1}B+D 를 직접 다루는 대신, 해당 행렬을 구현하는 디스크립터 실현 (A−λE, B, C, D) 로 변환하는 ‘실현 문제’를 강조한다. 실현은 유일하지 않으며, U·V 가역 변환에 의해 동등한 실현들이 존재한다(유사 변환). 이러한 자유도는 최소 차원 실현을 찾는 데 활용되며, 최소 실현은 다음 다섯 가지 조건을 만족한다: (i) 전·무한 제어 가능성, (ii) 전·무한 관측 가능성, (iii) E의 영공간이 E의 범위에 포함되는 비동적 모드 부재 등.

핵심 수치 기법으로는 두 종류의 정규형이 제시된다. 첫째, 정규 펜슬 A−λE에 대해 Weierstrass 정규형을 이용해 유한·무한 고유값을 명시적으로 분리한다. 둘째, 일반(특이) 펜슬 M−λN에 대해 Kronecker 정규형을 적용해 좌·우 특이 구조(K_r, K_reg, K_l)와 고유값 구조를 동시에 파악한다. 두 정규형 모두 비직교 변환을 사용하면 수치적으로 불안정하므로, 논문은 완전 직교(또는 정규) 변환만을 이용한 ‘Kronecker‑like’ 형태를 제안한다. 이러한 형태는 SVD·QR 기반 알고리즘으로 구현 가능하며, 고유값·특이값 계산에 내재된 병렬화와 안정성을 확보한다.

또한 정상 차수(normal rank), 극점(poles), 영점(zeros) 등 시스템 특성을 디스크립터 실현을 통해 직접 계산하는 절차를 제시한다. 정상 차수는 시스템 매트릭스 펜슬 S(λ)=