14,085 research outputs found
Order reduction methods for solving large-scale differential matrix Riccati equations
We consider the numerical solution of large-scale symmetric differential
matrix Riccati equations. Under certain hypotheses on the data, reduced order
methods have recently arisen as a promising class of solution strategies, by
forming low-rank approximations to the sought after solution at selected
timesteps. We show that great computational and memory savings are obtained by
a reduction process onto rational Krylov subspaces, as opposed to current
approaches. By specifically addressing the solution of the reduced differential
equation and reliable stopping criteria, we are able to obtain accurate final
approximations at low memory and computational requirements. This is obtained
by employing a two-phase strategy that separately enhances the accuracy of the
algebraic approximation and the time integration. The new method allows us to
numerically solve much larger problems than in the current literature.
Numerical experiments on benchmark problems illustrate the effectiveness of the
procedure with respect to existing solvers
Galerkin approximations for the optimal control of nonlinear delay differential equations
Optimal control problems of nonlinear delay differential equations (DDEs) are
considered for which we propose a general Galerkin approximation scheme built
from Koornwinder polynomials. Error estimates for the resulting
Galerkin-Koornwinder approximations to the optimal control and the value
function, are derived for a broad class of cost functionals and nonlinear DDEs.
The approach is illustrated on a delayed logistic equation set not far away
from its Hopf bifurcation point in the parameter space. In this case, we show
that low-dimensional controls for a standard quadratic cost functional can be
efficiently computed from Galerkin-Koornwinder approximations to reduce at a
nearly optimal cost the oscillation amplitude displayed by the DDE's solution.
Optimal controls computed from the Pontryagin's maximum principle (PMP) and the
Hamilton-Jacobi-Bellman equation (HJB) associated with the corresponding ODE
systems, are shown to provide numerical solutions in good agreement. It is
finally argued that the value function computed from the corresponding reduced
HJB equation provides a good approximation of that obtained from the full HJB
equation.Comment: 29 pages. This is a sequel of the arXiv preprint arXiv:1704.0042
Recommended from our members
Asymptotic and numerical analysis of the optimal investment strategy for an insurer
Combination of direct methods and homotopy in numerical optimal control: application to the optimization of chemotherapy in cancer
We consider a state-constrained optimal control problem of a system of two
non-local partial-differential equations, which is an extension of the one
introduced in a previous work in mathematical oncology. The aim is to minimize
the tumor size through chemotherapy while avoiding the emergence of resistance
to the drugs. The numerical approach to solve the problem was the combination
of direct methods and continuation on discretization parameters, which happen
to be insufficient for the more complicated model, where diffusion is added to
account for mutations. In the present paper, we propose an approach relying on
changing the problem so that it can theoretically be solved thanks to a
Pontryagin Maximum Principle in infinite dimension. This provides an excellent
starting point for a much more reliable and efficient algorithm combining
direct methods and continuations. The global idea is new and can be thought of
as an alternative to other numerical optimal control techniques
Numerical Analysis
Acknowledgements: This article will appear in the forthcoming Princeton Companion to Mathematics, edited by Timothy Gowers with June Barrow-Green, to be published by Princeton University Press.\ud
\ud
In preparing this essay I have benefitted from the advice of many colleagues who corrected a number of errors of fact and emphasis. I have not always followed their advice, however, preferring as one friend put it, to "put my head above the parapet". So I must take full responsibility for errors and omissions here.\ud
\ud
With thanks to: Aurelio Arranz, Alexander Barnett, Carl de Boor, David Bindel, Jean-Marc Blanc, Mike Bochev, Folkmar Bornemann, Richard Brent, Martin Campbell-Kelly, Sam Clark, Tim Davis, Iain Duff, Stan Eisenstat, Don Estep, Janice Giudice, Gene Golub, Nick Gould, Tim Gowers, Anne Greenbaum, Leslie Greengard, Martin Gutknecht, Raphael Hauser, Des Higham, Nick Higham, Ilse Ipsen, Arieh Iserles, David Kincaid, Louis Komzsik, David Knezevic, Dirk Laurie, Randy LeVeque, Bill Morton, John C Nash, Michael Overton, Yoshio Oyanagi, Beresford Parlett, Linda Petzold, Bill Phillips, Mike Powell, Alex Prideaux, Siegfried Rump, Thomas Schmelzer, Thomas Sonar, Hans Stetter, Gil Strang, Endre Süli, Defeng Sun, Mike Sussman, Daniel Szyld, Garry Tee, Dmitry Vasilyev, Andy Wathen, Margaret Wright and Steve Wright
- …