[논문 리뷰] Optimal H2 order-one reduction by solving eigenproblems for polynomial equations

Bernard Hanzon, J.M. Maciejowski외 1명|ArXiv.org|2007. 06. 13.

Numerical methods for differential equations참고 문헌 32인용 수 4

한 줄 요약

이 논문은 SISO 안정 선형 시스템의 H₂ 순서일치 모델 축소를 위한 구조적 대수 기반 방법을 제시하며, 일阶 최적성 조건에서 유도된 고유값 문제를 해결하여 국소 최솟값을 피한다. 그뢰브너 기저 구성과 선형 대수를 통해 전역 최적을 보장하며, 이전의 대수적 방법에 비해 계산 비용을 크게 감소시켜 고차수 시스템의 효율적 해법을 가능하게 한다.

ABSTRACT

A method is given for solving an optimal H2 approximation problem for SISO linear time-invariant stable systems. The method, based on constructive algebra, guarantees that the global optimum is found; it does not involve any gradient-based search, and hence avoids the usual problems of local minima. We examine mostly the case when the model order is reduced by one, and when the original system has distinct poles. This case exhibits special structure which allows us to provide a complete solution. The problem is converted into linear algebra by exhibiting a finite-dimensional basis for a certain space, and can then be solved by eigenvalue calculations, following the methods developed by Stetter and Moeller. The use of Buchberger's algorithm is avoided by writing the first-order optimality conditions in a special form, from which a Groebner basis is immediately available. Compared with our previous work the method presented here has much smaller time and memory requirements, and can therefore be applied to systems of significantly higher McMillan degree. In addition, some hypotheses which were required in the previous work have been removed. Some examples are included.

연구 동기 및 목표

• 기울기 기반 접근 방식에서 유래하는 국소 최솟값을 피하는 전역 최적 H₂ 모델 축소 방법을 개발하는 것.
• 특히 원래 시스템의 분리된 극점이 있는 경우, SISO 안정 시스템의 차수를 한 단계 감소시킨 H₂ 근사 문제를 해결하는 것.
• 이전 연구에서 제기된 제약 조건(비판점의 유한성 등)을 제거하면서도 전역 최적성을 유지하는 것.
• 기존 대수적 방법보다 계산 효율성이 높은 대안을 제공하여 고차수 맥밀란 차수를 가진 시스템에 적용 가능하게 하는 것.
• 특히 2^N - 1개의 비판점 수에 대한 상한을 포함한, 차수 감소 문제의 구조적 이론적 통찰을 제공하는 것.

제안 방법

• H₂ 근사 문제를 전이 함수 형태의 일阶 최적성 조건에서 유도된 다항방정식계로 변환한다.
• 특정 구조—'대각선 2차' 방정식—을 활용하여, 정의 다항식이 그뢰브너 기저를 이룬다는 점을 이용해 부버거 알고리즘 없이 직접 해를 구한다.
• 해공간이 유한차원임을 보이며, 선형 대수 기법을 통해 고유값 문제로 환원할 수 있음을 보여준다.
• 시스템 데이터로부터 행렬을 구성하고, 이 행렬의 일반화된 고유값 문제를 풀어 최적 근사치를 도출한다.
• 수치적 기울기 최적화를 피하고, 기호 계산과 수치적 고유값 해법에 의존함으로써 효율성과 정확성을 확보한다.
• 예제를 통한 검증을 통해 이론적 방법의 타당성을 입증하였으며, 비판점 행렬 기반 알고리즘과의 비교 분석도 수행하였다.

실험 결과

연구 질문

• RQ1 반복 기울기 기반 최적화에 의존하지 않고도 전역 최적 H₂ 순서일치 축소를 계산할 수 있는가?
• RQ2 H₂ 차수 감소 문제의 비판점 수의 최대값은 얼마이며, 이를 상한으로 둘 수 있는가?
• RQ3 H₂ 모델 축소를 위한 대수적 방법의 계산 복잡도를 전역 최적성을 유지하면서 크게 감소시킬 수 있는가?
• RQ4 알고리즘이 실패하는 조건은 무엇이며, 특히 고유벡터의 조건수와 실수 허리츠 안정성과의 관련성은 무엇인가?
• RQ5 이 방법을 중복 극점을 가진 시스템이나 고차수 축소로 확장할 수 있으며, 어떤 이론적 통찰을 제공하는가?

주요 결과

• 이 방법은 분리된 극점을 가진 SISO 안정 시스템에 대해 H₂ 순서일치 축소에서 전역 최적성을 보장하며, 기울기 기반 방법의 국소 최솟값 문제를 피한다.
• 차수 감소 문제의 비판점 수는 최대 2^N - 1개이며, 이는 유한하고 유계된 해공간을 제공한다.
• N=5일 때, α ∈ [0.38, 0.79] 범위에서 최적 근사치를 성공적으로 계산하며, (0, 0.38) 및 (0.79, 1) 범위에서는 고유벡터의 악조건화 또는 비실수 허리츠 해로 인해 실패한다.
• α=0.78일 때, 상대 오차 1.96%를 가진 단일 타당한 해를 반환하며, 이 근사치는 리라크세이션 시스템으로 확인된다.
• 기존 대수적 접근에 비해 계산 효율성이 뛰어나 고차수 맥밀란 차수를 가진 시스템의 해를 구하는 데 유리하다.
• 고유값 기반 알고리즘과 비판점 행렬 기반 알고리즘 모두 일관된 결과를 도출하지만, 수치 정밀도 측면에서 α=1 근처에서는 고유값 기반 방법이 더 우수한 안정성을 보인다.