OS Numerical Optimization: Certified Reduced-Order Model Predictive Control for Linear Switched Evolution Equations
Time
Tuesday, 12. November 2024
15:15 - 16:45
Location
F426
Organizer
S. Volkwein
Speaker:
Michael Kartmann
On 12th November 2024 at 15:15, Michael Kartmann from the University of Konstanz will give a talk.
Abstract: In this talk, we are concerned with the solution of infinite-horizon optimal control problems of the form
\begin{equation}
\begin{aligned}
\min_{u \in L^2((0,\infty),\mathbb{R}^\rho)} \mathcal{J}(u) = \int_{0}^{\infty} \tfrac{1}{2} \|y(t)-y_d(t)\|^2_{\mathbb{R}^{\rho}}+ \tfrac{\lambda}{2}\|u(t)\|_{\mathscr{U}_n}^2 + \mu \|u(t) \|_{L^1(\Omega)}\,\mathrm{d}t,
\end{aligned}
\end{equation}
subject to $(u, y)$ solving the following linear switched input-output system
\begin{equation}
\left\{\quad
\begin{aligned}
\mathcal{M}_{\sigma(t)}\tfrac{\mathrm{d}}{\mathrm{dt}} \theta(t) + \mathcal{A}_{\sigma(t)}\theta(t) &= \mathcal{B}_{\sigma(t)}u(t) &&t\geq 0\\
y(t) &= \mathcal{C}_{\sigma(t)}\theta(t) && t\geq 0\\
\theta(0) &= \theta_{\circ}. && \\
\end{aligned}\right.
\end{equation}
and control constraints $u\in \mathscr{U}_{\mathsf{ad}}$. Here $\sigma:[0,\infty) \to \{1,\ldots,L\}$ is a switching signal, that switches through different system operators $\mathcal{M}_i,\mathcal{A}_i,\mathcal{B}_i,\mathcal{C}_i$ for $i=1,\ldots,L$. \\
To approximate the solution of $(1)$, we apply Model Predictive Control (MPC): the optimal control problem is solved over smaller, receding time intervals $(t_n,t_n+T)$ for some prediction horizon $T>0$ and the solutions are concatenated in the sampling interval $(t_n,t_{{n+1}})$ for $0<t_{{n+1}}<t_n+T$. First, we derive optimality conditions for these small-horizon problems and discuss their suboptimality w.r.t. $(1)$. The difficulty here is that the cost functional $\mathcal{J}$ is not differentiable in the classical sense, due to the presence of the $L^1$-regularization. Second, the repeated solution of small-horizon optimal control problems motivates model reduction: we consider (Petrov-)Galerkin reduced-order models for $(2)$ to speed up the MPC process. To quantify the error, we do a full a posteriori error analysis for the optimal control, optimal state, and optimal value function of the small-horizon problems, which allows us to control the evolving error through the MPC iterations. These estimates are then used to construct two certified ROM-MPC algorithms for the solution of $(1)$, that are up to $10$ times faster than the MPC relying on the full-order model. This is joint work with Stefan Volkwein (U. Konstanz), Mattia Manucci, and Benjamin Unger (U. Stuttgart).