[논문 리뷰] Optimal Control of a Parabolic Distributed Parameter System Using a Barycentric Shifted Gegenbauer Pseudospectral Method

Kareem T. Elgindy|arXiv (Cornell University)|2016. 03. 04.

Advanced Numerical Methods in Computational Mathematics참고 문헌 2인용 수 3

한 줄 요약

이 논문은 포물형 분포 매개변수 시스템에 의해 지배되는 최적 제어 문제를 해결하기 위해 새로운 바리센트릭 이동 Gegenbauer 양자법(BCGPM)을 제안한다. 문제를 적분 형태로 재구성하고, 이동 Gegenbauer-Gauss 구분점과 안정적인 바리센트릭 보간을 사용함으로써, 최소한의 구분점 수로도 지수 수렴을 달성하며 높은 정확도를 확보한다. 이는 기존 방법들보다 효율성과 정확도 면에서 뛰어나다.

ABSTRACT

In this paper, we introduce a novel pseudospectral method for the numerical solution of optimal control problems governed by a parabolic distributed parameter system. The infinite-dimensional optimal control problem is reduced into a finite-dimensional nonlinear programming problem through shifted Gegenbauer quadratures constructed using a stable barycentric representation of Lagrange interpolating polynomials and explicit barycentric weights for the shifted Gegenbauer-Gauss (SGG) points. A rigorous error analysis of the method is presented, and a numerical test example is given to show the accuracy and efficiency of the proposed pseudospectral method.

연구 동기 및 목표

• 포물형 편미분방정식에 의해 지배되는 최적 제어 문제를 해결하기 위한 고차수이고 안정적인 수치적 방법을 개발한다.
• 적분 재구성과 스펙트럼 이산화를 통해 무한차원 최적 제어 문제를 잘 조절된 유한차원 비선형 프로그래밍 문제로 축소한다.
• 안정적인 바리센트릭 보간과 이동 Gegenbauer 구분점들을 활용하여 최소한의 구분점 수로 지수 수렴률을 달성한다.
• 기존의 라디얼 기저 함수 방법보다 뛰어난 정확도와 계산 효율성을 제공함으로써, 이에 비해 개선된 성능을 제공한다.
• 일반적인 PDE 제약 최적 제어 문제에 적용 가능한 강력한 프레임워크를 제공한다.

제안 방법

• 상태 및 제어 변수의 연속적 적분을 이용하여 원래의 PDE 제약 최적 제어 문제를 적분 형태로 재구성한다.
• 이동 Gegenbauer-Gauss (SGG) 점을 구분점으로 사용하며, 명시적인 바리센트릭 가중치를 가진 라그랑주 보간 다항식의 안정적인 바리센트릭 표현을 활용한다.
• 이동 Gegenbauer 구분점을 사용하여 적분 행렬을 구성함으로써 적분 연산자를 정확하고 효율적으로 근사한다.
• 무한차원 문제를 선형 제약 조건을 가진 유한차원 비선형 프로그래밍 문제로 이산화하며, 표준 최적화 솔버로 해결할 수 있다.
• Gegenbauer 다항식의 스펙트럼 성질과 바리센트릭 보간의 잘 조절된 성질을 활용하여 고차수 정확도를 확보한다.
• 엄밀한 오차 분석을 수행하여, 구분점 수가 증가함에 따라 오차가 스펙트럼적으로 감소함을 증명한다.

실험 결과

연구 질문

• RQ1 바리센트릭 이동 Gegenbauer 양자법이 포물형 최적 제어 문제에서 지수 수렴을 달성할 수 있는가?
• RQ2 이동 Gegenbauer 다항식의 형태 매개변수 α의 선택이 해의 정확도와 안정성에 어떤 영향을 미치는가?
• RQ3 동일한 문제에 대해 제안된 방법이 라디얼 기저 함수 방법보다 정확도와 계산 비용 면에서 뛰어나게 성능을 높일 수 있는가?
• RQ4 구분점 수가 증가함에 따라 이 방법의 수렴 행동은 어떠한가?
• RQ5 이산화 형식에서 초기 조건 및 경계 조건과 같은 제약 조건을 얼마나 잘 유지하는가?

주요 결과

• BCGPM는 지수 수렴률을 달성하며, 구분점 수가 증가함에 따라 오차가 스펙트럼적으로 감소한다.
• 공간 및 시간 방향 각각 N=5개의 구분점만으로도, 각 방향에 120개의 노드를 사용하는 라디얼 기저 함수 방법과 비교해 유사한 정확도를 달성한다.
• 작은 N 값에서, 음수 또는 0 이하의 α 값(예: α = -0.2)을 가진 SGG 점에서의 이산화가 초기 조건(ψ₁)에서 최대 오차가 가장 작게 나타나, 소규모 문제에 대해 뛰어난 정확도를 보인다.
• 모든 테스트 구성에서 해의 타당성과 최대 경계 조건 위반(ψ₂)은 기계 에프로크스미메이션 수준에 가까이 유지되어 강력한 제약 조건 충족을 나타낸다.
• 모든 테스트된 N 및 α 값에서 근사 최적 비용 기능 Jₙ,ₙ*는 약 15로 안정화되며, 수치적 일관성을 확인한다.
• 그림 2 및 3에서 시각화된 바와 같이, N=12 및 α=-0.2 조건에서도 상태 및 제어 프로파일이 매우 정확하게 생성되며, 최소한의 격자 점으로도 고해상도 근사가 가능함을 보여준다.