Solving Singular Control Problems in Mathematical Biology, Using PASA
Abstract
In this paper, we will demonstrate how to use a nonlinear polyhedral constrained optimization solver called the Polyhedral Active Set Algorithm (PASA) for solving a general singular control problem. We present methods of discretizing a general optimal control problem that involves the use of the gradient of the Lagrangian for computing the gradient of the cost functional so that PASA can be applied. When a numerical solution contains artifacts that resemble "chattering'', a phenomenon where the control oscillates wildly along the singular region, we recommend a method of regularizing the singular control problem by adding a term to the cost functional that measures a scalar multiple of the total variation of the control, where the scalar is viewed as a tuning parameter. We then demonstrate PASA's performance on three singular control problems that give rise to different applications of mathematical biology. We also provide some exposition on the heuristics that we use in determining an appropriate size for the tuning parameter.
Keywords
Cite
@article{arxiv.2010.06744,
title = {Solving Singular Control Problems in Mathematical Biology, Using PASA},
author = {Summer Atkins and Mahya Aghaee and Maia Martcheva and William Hager},
journal= {arXiv preprint arXiv:2010.06744},
year = {2020}
}