[논문 리뷰] Energy-Preserving and Passivity-Consistent Numerical Discretization of Port-Hamiltonian Systems

Elena Celledoni, Eirik Hoel Høiseth|arXiv (Cornell University)|2017. 06. 26.

Control and Stability of Dynamical Systems참고 문헌 17인용 수 15

한 줄 요약

이 논문은 포트-해밀토니안 시스템을 이산화할 때 이산적 수동성과 에너지 균형을 유지하는 두 가지 구조 유지 수치 방법—이산 그래디언트 및 스플리팅 방법—을 제안한다. 주요 기여는 수성 안정성 보장이 가능한 이산 포트-해밀토니안 시스템으로, 수성 일致 제어를 통해 전역적 점근적 안정성을 확보한 것이다. 수치 실험을 통해 기존 동일 차수의 적분기들과 비교해 뛰어난 성능을 입증하였다.

ABSTRACT

In this paper we design discrete port-Hamiltonian systems systematically in two different ways, by applying discrete gradient methods and splitting methods respectively. The discrete port-Hamiltonian systems we get satisfy a discrete notion of passivity, which lets us, by choosing the input appropriately, make them globally asymptotically stable with respect to an equilibrium point. We test methods designed using the discrete gradient approach in numerical experiments, and the results are encouraging when compared to relevant existing integrators of identical order.

연구 동기 및 목표

• 연속적 포트-해밀토니안 시스템의 본질적 에너지 균형 및 수성 성질을 유지하는 수치 이산화 기법을 개발하기 위해.
• 이산 수성과 리아푸노프 분석을 활용하여 이산 포트-해밀토니안 시스템의 전역 점근적 안정성을 보장하기 위해.
• 원래 시스템의 기하학적 및 상호연결 성질을 수치적 시간 적분 하에서도 유지하는 구조 유지 적분기 설계를 위해.
• 기존의 고차수 적분기들과 동일한 차수에서의 성능 비교를 통해 제안된 방법을 수치 실험으로 검증하기 위해.

제안 방법

• 이산 그래디언트 방법은 이산 해밀토니안을 구성하고, 라그랑주 보간을 활용한 콜로케이션 기반 적분을 통해 단계 값과 제어 입력을 정의한다.
• 양방향 규칙을 사용한 적분 공식을 통해 이산 에너지 균형 식을 유도하며, 양의 가중치를 갖는다. 이는 $ H(x_{n+1}) - H(x_n) = h \sum_{j=1}^s b_j y_j^T u_j $ 를 만족시켜 이산 수성을 보장한다.
• 스플리팅 방법을 사용하여 시스템을 간단한 구성 요소로 분해함으로써, 기저가 되는 다이어크 구조를 유지하고 상호연결 하에서 수성을 보장한다.
• 이산 시스템이 반경 방향 무한대 저장 함수 $ H $ 를 갖는 수성임을 보여주며, 리아푸노프 함수를 통한 안정성 분석이 가능하다.
• 제어 입력은 $ u_j = -\phi(y_j) $ 로 선택되며, 여기서 $ \phi(0) = 0 $ 이고 $ y^T \phi(y) > 0 $ 이다. $ y \neq 0 $ 인 경우 에너지 변화가 음의 정부호가 되어 점근적 안정성을 보장한다.
• 이론적 분석을 통해 모든 단계에서 출력이 0인 유일한 해가 자명한 해 $ x_n = 0 $ 인 경우에만 원점이 전역적으로 점근적으로 안정하다고 증명한다.

실험 결과

연구 질문

• RQ1 이산 그래디언트 방법을 체계적으로 적용하여 에너지 보존 및 수성 일치 이산 포트-해밀토니안 시스템을 구성할 수 있는가?
• RQ2 스플리팅 방법을 어떻게 활용하여 이산 시간 영역에서 수성과 구조를 유지할 수 있는가?
• RQ3 제안된 이산 시스템은 수성 제어 입력 설계 하에 전역 점근적 안정성을 보장하는가?
• RQ4 콜로케이션 기반 시간 이산화의 맥락에서 이산 에너지 균형과 수성 간의 관계는 무엇인가?
• RQ5 기존의 동일 차수 고차수 적분기들과 비교해 제안된 방법의 성능은 어떠한가?

주요 결과

• 이산 그래디언트 방법을 통해 구성된 이산 포트-해밀토니안 시스템은 이산 에너지 균형 식 $ H(x_{n+1}) - H(x_n) = h \sum_{j=1}^s b_j y_j^T u_j $ 를 만족하며, 이는 수성을 보장한다.
• 제어 입력 $ u_j = -\phi(y_j) $ 가 $ y^T \phi(y) > 0 $ 를 만족하는 경우, 이산 시스템은 $ \Delta H(x_n) \leq 0 $ 를 보장하며, 이로 인해 $ H $ 가 리아푸노프 함수가 된다.
• 모든 단계에서 출력이 0인 해가 자명한 해 $ x_n = 0 $ 인 경우에만 원점의 전역 점근적 안정성이 증명된다.
• 수치 실험을 통해 제안된 방법이 동일 차수의 기존 적분기들보다 에너지 보존 및 안정성 측면에서 뛰어난 성능을 보임을 입증하였다.
• 이 방법은 일반적인 입력-상태-출력 포트-해밀토니안 시스템에 적용 가능하며, 상호연결 하에서도 다이어크 구조를 유지한다.
• 추가 자료에서 이 방법의 고차수 일반화가 유도되었으며, 이는 고차수 정확도를 갖는 방법으로의 적용 가능성을 확장한다.