Physics-informed Gaussian Processes as Linear Model Predictive Controller

Jörn Tebbe,Andreas Besginow,Markus Lange-Hegermann
2024-12-02
Abstract:We introduce a novel algorithm for controlling linear time invariant systems in a tracking problem. The controller is based on a Gaussian Process (GP) whose realizations satisfy a system of linear ordinary differential equations with constant coefficients. Control inputs for tracking are determined by conditioning the prior GP on the setpoints, i.e. control as inference. The resulting Model Predictive Control scheme incorporates pointwise soft constraints by introducing virtual setpoints to the posterior Gaussian process. We show theoretically that our controller satisfies asymptotical stability for the optimal control problem by leveraging general results from Bayesian inference and demonstrate this result in a numerical example.
Optimization and Control,Machine Learning,Systems and Control
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: How to design an accurate and stable tracking controller for a linear time - invariant system in the presence of input and state constraints. Specifically, the author proposes a model predictive control (MPC) method based on Physics - informed Gaussian Processes (LODE - GP) to achieve the following goals: 1. **Satisfy the strict constraints of the linear ordinary differential equation (ODE) system**: Ensure that the trajectory generated by the controller strictly follows the given linear ODE system. 2. **Handle hard and soft constraints**: Handle hard constraints (such as initial and final states) by introducing virtual set - points, and ensure the feasibility of the optimization problem through soft constraints (such as the bounds of states and control inputs). 3. **Guarantee asymptotic stability**: Theoretically prove that the proposed controller can make the system asymptotically stable, and verify this conclusion in numerical experiments. ### Specific description of the problem The control problem discussed in the paper is the tracking control problem for linear time - invariant (LTI) systems. The goal of such problems is to minimize the difference between the system state and the reference trajectory while satisfying the dynamic equations and various constraints of the system. The specific mathematical description is as follows: Consider a controllable linear ODE system, whose state \(x\in\mathbb{R}^{n_x}\), control input \(u\in\mathbb{R}^{n_u}\), output \(y\in\mathbb{R}^{n_y}\), is described by the following equations: \[ \dot{x}=Ax + Bu\quad(1) \] \[ y = Cx\quad(2) \] where \(A\in\mathbb{R}^{n_x\times n_x}\), \(B\in\mathbb{R}^{n_x\times n_u}\), \(C\in\mathbb{R}^{n_y\times n_x}\) are the system, control, and output matrices respectively. The state and control input of the system need to satisfy the following constraints: \[ x_{\min}\leq x(t)\leq x_{\max}\quad(3) \] \[ u_{\min}\leq u(t)\leq u_{\max}\quad(4) \] for the time interval \(t\in[t_0,t_T]\). The general tracking control task is to minimize the difference between the given reference trajectory \(x_{\text{ref}}\in\mathbb{R}^{n_x}\) and the system state \(x\) within a given time range: \[ \min_{u(t)}\int_{t_0}^{t_T}\|x_{\text{ref}}-x(t)\|dt\quad(5) \] The tracking control task in discrete time steps can be expressed as: \[ \min_{u(t)}\sum_{i = 0}^T\|x_{\text{ref}}-x(t_i)\|+\|u\|\quad(6a) \] \[ \text{s.t.}\quad\dot{x}=Ax + Bu,\quad x(t_0)=x_0,\quad(6b),(6c) \] \[ x_{\min}\leq x(t)\leq x_{\max},\quad\forall t\in[t_0,t_T],\quad(6d) \] \[ u_{\min}\leq u(t)\leq u_{\max},\quad\forall t\in[t_0,t_T].\quad(6e) \] ### Proposed method The paper proposes an MPC method based on LODE - GP to solve the above problems through the following steps: 1. **Construct L**