Multiscale differential riccati equations for linear quadratic regulator problems
Journal article, 2018

We consider approximations to the solutions of differential Riccati equations in the context of linear quadratic regulator problems, where the state equation is governed by a multiscale operator. Similarly to elliptic and parabolic problems, standard finite element discretizations perform poorly in this setting unless the grid resolves the fine-scale features of the problem. This results in unfeasible amounts of computation and high memory requirements. In this paper, we demonstrate how the localized orthogonal decomposition method may be used to acquire accurate results also for coarse discretizations, at the low cost of solving a series of small, localized elliptic problems. We prove second-order convergence (except for a logarithmic factor) in the L2operator norm and first-order convergence in the corresponding energy norm. These results are both independent of the multiscale variations in the state equation. In addition, we provide a detailed derivation of the fully discrete matrix-valued equations and show how they can be handled in a low-rank setting for large-scale computations. In connection to this, we also show how to efficiently compute the relevant operator-norm errors. Finally, our theoretical results are validated by several numerical experiments.

Differential Riccati equations


Linear quadratic regulator problems

Finite elements

Localized orthogonal decomposition


Axel Målqvist

Chalmers, Mathematical Sciences, Applied Mathematics and Statistics

Anna Persson

Chalmers, Mathematical Sciences, Applied Mathematics and Statistics

Tony Stillfjord

Chalmers, Mathematical Sciences, Applied Mathematics and Statistics

SIAM Journal of Scientific Computing

1064-8275 (ISSN) 1095-7197 (eISSN)

Vol. 40 4 A2406-A2426

Subject Categories

Computational Mathematics

Control Engineering

Mathematical Analysis



More information

Latest update