[논문 리뷰] A Pseudospectral Approach to High Index DAE Optimal Control Problems

Harleigh C. Marsh, Mark Karpenko외 1명|arXiv (Cornell University)|2018. 11. 30.

Numerical methods for differential equations참고 문헌 24인용 수 5

한 줄 요약

이 논문은 특수스펙트럼(PS) 방법, 특히 DIDO MATLAB 도구상자 via를 통해 고차수 미분대수방정식(DAE) 최적제어 문제—예를 들어 인덱스-3 펜듈럼 문제—를 인덱스 감소 없이 직접 해결할 수 있음을 보여준다. 이는 오랫동안 이러한 문제들이 본질적으로 어려운 것으로 여겨져 온 믿음을 도전한다. 핵심 기여는 PS 방법이 전통적인 계산적 장벽을 우회함으로써, DAE의 어려움에 대한 새로운 이론이 필요하다는 것을 시사한다는 점이다.

ABSTRACT

Historically, solving optimal control problems with high index differential algebraic equations (DAEs) has been considered extremely hard. Computational experience with Runge-Kutta (RK) methods confirms the difficulties. High index DAE problems occur quite naturally in many practical engineering applications. Over the last two decades, a vast number of real-world problems have been solved routinely using pseudospectral (PS) optimal control techniques. In view of this, we solve a "provably hard," index-three problem using the PS method implemented in DIDO, a state-of-the-art MATLAB optimal control toolbox. In contrast to RK-type solution techniques, no laborious index-reduction process was used to generate the PS solution. The PS solution is independently verified and validated using standard industry practices. It turns out that proper PS methods can indeed be used to "directly" solve high index DAE optimal control problems. In view of this, it is proposed that a new theory of difficulty for DAEs be put forth.

연구 동기 및 목표

• 고차수 DAE 최적제어 문제의 본질적 어려움 여부 또는 특정 계산 방법에 기인한 어려움 여부를 조사한다.
• 기존에 인덱스 감소 없이 룬지-쿠타(RK) 방법으로 해결되지 않았던 기준 인덱스-3 DAE 문제에 대해 특수스펙트럼(PS) 방법의 효용성을 시험한다.
• 실제 공학 문제에서 관찰되는 특이 궤적과 대수적 제약 조건을 포함한 고차수 DAE를 PS 방법이 직접 처리할 수 있음을 검증한다.
• 고차수 DAE 최적제어 문제를 해결하기 위해 인덱스 감소가 반드시 필요하다는 전통적 믿음을 도전한다.
• 문제의 내재적 복잡성보다는 계산 방법에 의존하는 DAE의 어려움에 대한 새로운 이론을 제안한다.

제안 방법

• 이 연구는 특수스펙트럼 최적제어를 코베이터 매핑 원리(CMP)를 통해 필수 조건을 유도하는 데 사용하는 DIDO MATLAB 최적제어 도구상자를 사용한다.
• 인덱스-3 DAE 최적제어 문제는 상태 변수 $x_1, x_2, x_3, x_4, x_5$ 를 사용하여 표준 상태공간 형태로 기술되며, $x_5$ 는 대수적 제약 조건 $x_1^2 + x_3^2 = L^2$ 의 라그랑주 승수를 나타낸다.
• 이차 특수스펙트럼 미분 행렬을 적용하여 상태 및 코스테이트 궤적을 스펙트럼 정확도로 근사함으로써 고정밀도 최적제어 문제 해결이 가능해진다.
• 해결책은 필수 조건—종단 전이 조건 및 제어 조건 (식 9 및 10)—이 수치적으로 만족되는지 확인하기 위해 코스테이트 및 제어 표현의 도표를 통해 검증된다.
• 독립적 검증은 표준 산업 관행을 사용하여 수행되었으며, 이는 인덱스 감소나 문제 변환 없이도 PS 해의 신뢰성을 확인한다.
• PS 방법의 수렴 이론은 스톤-바이어슈트라스 정리와 아르체라-아스콜리 정리에 기반하며, 인덱스나 특이성 가정에 의존하지 않는다.

실험 결과

연구 질문

• RQ1 특수스펙트럼 방법은 인덱스 감소나 문제 변환 없이 고차수 DAE 최적제어 문제를 직접 해결할 수 있는가?
• RQ2 왜 룬지-쿠타 기반 방법은 특수스펙트럼 방법이 성공하는 고차수 DAE 문제에서는 실패하는가?
• RQ3 고차수 DAE 최적제어 문제의 어려움은 문제 자체에 내재된 것이며, 수치적 방법의 선택에 기인한 것인가?
• RQ4 특이 궤적과 고차수 DAE는 특수스펙트럼 방법이 수렴성이나 정확도를 저하시키지 않고 효과적으로 처리할 수 있는가?
• RQ5 복잡한 실제 DAE 문제를 해결하는 데 있어 특수스펙트럼 방법의 성공은 현재 DAE 어려움 이론을 재고할 필요가 있는가?

주요 결과

• DIDO에 구현된 특수스펙트럼 방법은 인덱스 감소나 문제 변환 없이도 기준 인덱스-3 DAE 최적제어 문제를 성공적으로 해결했다.
• 해결책은 종단 전이 조건 및 제어 최적성 조건을 모두 만족하며, 이는 코스테이트 및 제어 표현의 수치적 도표를 통해 확인되었다.
• DIDO에서 생성된 코스테이트는 종단 시간에 0이 되어 종단 전이 조건이 만족됨을 확인했다.
• 필수 조건에서 유도된 제어 법칙 (식 9) 은 수치적으로 확인되었으며, 그림 8에서 양변이 일치함을 보였다.
• 특이 제어 조건 (식 10) 도 수치적으로 만족되었으며, 그림 9에서 이를 확인하여 해가 최적성 조건과 일관됨을 입증했다.
• 이 연구는 고차수 DAE 문제를 적절히 구현된 특수스펙트럼 방법을 사용하면 본질적으로 어려운 것이 아니며, 문제의 내재적 복잡성보다는 계산 방법 의존성에 기반한 DAE 어려움에 대한 새로운 이론이 필요하다고 결론 내린다.