[논문 리뷰] Convergence rate for a Radau collocation method applied to unconstrained optimal control

William W. Hager, Hongyan Hou외 1명|arXiv (Cornell University)|2015. 08. 16.

Advanced Numerical Methods in Computational Mathematics인용 수 11

한 줄 요약

이 논문은 비제약 최적제어 문제에 Radau 경사법을 적용할 때 지수 수렴 속도를 확립한다. 한 끝점인 t=+1을 경사점으로 포함하는 Radau 구적점들을 사용함으로써, 매끄럽고 강력한 볼록성 조건 하에서 스펙트럼 수렴—즉, 경사점 수에 대해 지수적으로 감소하는 수렴—을 달성한다. 이는 구조적 차이가 있음에도 불구하고 가우스 기반 경사법과 동일한 수렴 속도를 갖는다.

ABSTRACT

A local convergence rate is established for an orthogonal collocation method based on Radau quadrature applied to an unconstrained optimal control problem. If the continuous problem has a sufficiently smooth solution and the Hamiltonian satisfies a strong convexity condition, then the discrete problem possesses a local minimizer in a neighborhood of the continuous solution, and as the number of collocation points increases, the discrete solution convergences exponentially fast in the sup-norm to the continuous solution. An earlier paper analyzes an orthogonal collocation method based on Gauss quadrature, where neither end point of the problem domain is a collocation point. For the Radau quadrature scheme, one end point is a collocation point.

연구 동기 및 목표

• 비제약 최적제어 문제에서 Radau 경사법의 국소 수렴 속도를 확립하기 위해.
• 경사점 수가 증가함에 따라 이산 해가 연속 해로 수렴하는 행동을 분석하기 위해.
• Radau 방법이 구조적 차이가 있음에도 불구하고 가우스 기반 경사법과 동일한 지수 수렴 속도를 달성함을 보여주기 위해.
• 최종 제약 조건의 구현이 유리하기 때문에 hp 이산화 프레임워크에서 Radau 경사법의 사용에 대한 이론적 근거를 제공하기 위해.

제안 방법

• 이 방법은 Radau 구적점 기반의 직교 경사법을 사용하며, 한 끝점(t=+1)을 경사점으로 포함한다.
• 이산 최적제어 문제는 N개의 점에서 경사진 운동 방정식과 초기 조건을 충족시키는 조건 하에서 비용 기능을 최소화하는 방식으로 설정된다.
• 상태와 제어는 최대 차수 N의 다항식으로 근사되며, 운동 방정식은 Radau 점들 τ₁ < … < τ_N = +1 에서 강제된다.
• 이산 수반 시스템이 유도되며, 미분 행렬 D^‡의 구조가 분석되어 역행성과 안정성을 보장한다.
• 이론적 분석은 라그랑주 보간과 직교 경사법의 성질에 기반하며, 특히 미분 행렬의 역행성에 대한 경계를 포함한다.
• 수렴 증명은 해밀토니안의 강력한 볼록성과 해의 매끄러움을 이용하여 균등 노름에서의 지수 감소를 확립한다.

실험 결과

연구 질문

• RQ1 표준적인 매끄러움과 볼록성 가정 하에서 Radau 경사법이 비제약 최적제어 문제에 대해 지수 수렴을 달성하는가?
• RQ2 경사점 분포의 구조적 차이가 있음에도 불구하고 Radau 방법의 수렴 속도가 가우스 경사법과 어떻게 비교되는가?
• RQ3 Radau 방법에서 최종 경사점의 역할은 무엇이며, 이는 최종 제약 조건의 조건화와 구현에 어떻게 영향을 미치는가?
• RQ4 경사점 수가 증가함에 따라 이산 해가 연속 해로 균등 노름에서 국소 수렴하는가?
• RQ5 미분 행렬 D^‡의 역행성에 대한 경계는 무엇이며, 이는 방법의 안정성과 수렴에 어떻게 영향을 미치는가?

주요 결과

• 경사점 수 N이 증가함에 따라 이산 해가 연속 해로 균등 노름에서 지수적으로 수렴한다.
• 연속 해가 η+1개의 연속 도함수를 가지며 η ≥ 3일 경우, 균등 노름에서 오차는 O(N^{2−η})로 경계된다.
• Radau 경사법의 수렴 속도는 Radau 방법이 비대칭 구조를 가지며 종단 경계에 경사점이 있다는 점을 감안할 때, 가우스 기반 경사법과 동일하다.
• 미분 행렬 D^‡의 역행성은 명시적으로 유도되었으며, 라그랑주 기저 함수와 구적 가중치를 통해 표현된다.
• 수치 실험은 관측된 수렴 속도가 이론적 예측과 일치하는 지수 수렴임을 확인한다.
• D^‡⁻¹과 D^‡⁻¹W⁻¹/²의 행 범위는 유계이며, N이 증가함에 따라 각각 2와 √2에 수렴함을 보여, 수치적 안정성을 시사한다.