I NFINITE -H ORIZON D IFFERENTIABLE M ODEL P REDICTIVE C ONTROL
Anonymous authors
Paper under double-blind review
A BSTRACT
This paper proposes a differentiable linear quadratic Model Predictive Control (MPC) framework for safe imitation learning.
The inﬁnite-horizon cost is enforced using a terminal cost function obtained from the discrete-time algebraic Riccati equation (DARE), so that the learned controller can be proven to be stabilizing in closed-loop.
A central contribution is the derivation of the analytical derivative of the solution of the DARE, thereby allowing the use of differentiation-based learning methods.
A further contribution is the structure of the MPC optimization problem: an augmented Lagrangian method ensures that the MPC optimization is feasible throughout training whilst enforcing hard constraints on state and input, and a prestabilizing controller ensures that the MPC solution and derivatives are accurate at each iteration.
The learning capabilities of the framework are demonstrated in a set of numerical studies.
I NTRODUCTION
Imitation Learning (IL, Osa et al., 2018) aims at reproducing an existing control policy by means of a function approximator and can be used, for instance, to hot-start reinforcement learning.
Effective learning and generalisation to unseen data are paramount to IL success, especially in safety critical applications.
Model Predictive Control (MPC, Maciejowski, 2000; Camacho & Bordons, 2007;
Rawlings & Mayne, 2009; Kouvaritakis & Cannon, 2015; Gallieri, 2016; Borrelli et al., 2017;
Rakovi ´ c & Levine, 2019) is the most successful advanced control methodology for systems with hard safety constraints .
At each time step, a ﬁnite horizon forecast is made from a predictive model of the system and the optimal actions are computed, generally relying on convex constrained Quadratic
Programming (QP, Boyd & Vandenberghe, 2004; Bemporad et al., 2000).
Stability of the MPC in closed loop with the physical system requires the solution of a simpler unconstrained inﬁnite horizon control problem (Mayne et al., 2000) which results in a value function (terminal cost and constraint) and a candidate terminal controller to be accounted for in the MPC forecasting.
For Linear
Time Invariant (LTI) models and quadratic costs, this means solving (ofﬂine) a Riccati equation (Kalman, 2001) or a linear matrix inequality (Boyd et al., 1994).
Under these conditions, an MPC controller will effectively control a system, up to a certain accuracy, provided that uncertainties in the model dynamics are limited (Limon et al., 2009).
Inaccuracies in the MPC predictions can reduce its effectiveness (and robustness) as the forecast diverges from the physical system trajectory over long horizons.
This is particularly critical in applications with both short and long-term dynamics and it is generally addressed, for instance in robust MPC (Richards, 2004; Rakovi´c et al., 2012), by using a controller to pre-stabilise the predictions.
This paper presents an inﬁnite-horizon differentiable linear quadratic MPC that can be learned using gradient-based methods.
In particular, the learning method uses an MPC controller where the terminal cost and terminal policy are the solution of an unconstrained inﬁnite-horizon Linear
Quadratic Regulator (LQR).
A closed-form solution for the derivative of the Discrete-time Algebraic
Riccati Equation (DARE) associated with the LQR is presented so that the stationary solution of the forward pass is fully differentiable.
This method allows analytical results from control theory to be used to determine the stabilizing properties of the learned controller when implemented in closed-loop.
Once the unconstrained LQR is computed, the predictive model is pre-stabilised using a linear state-feedback controller to improve the conditioning of the QP and the numerical accuracy of the MPC solution and gradients.
The proposed algorithm successfully learns an MPC with both local stability and intrinsic robustness guarantees under small model uncertainties.
Contributions
This paper provides a framework for correctly learning an inﬁnite-horizon, LTI quadratic MPC using recent developments in differentiable QPs (Amos & Kolter, 2017) and principles from optimal control (Blanchini & Miani, 2007).
A primary contribution is that the Discrete-time
Algebraic Riccati Equation (DARE) is used to provide inﬁnite-horizon optimality and stability, and an analytical derivative of the solution of the DARE is derived so that differentiation-based optimization can be used for training.
This connects known results on MPC stability (Limon et al., 2003; 2009) and on inﬁnite-horizon optimality (Scokaert & Rawlings, 1998) to imitation learning (Osa et al., 2018).
A further contribution is the MPC control formulation: a pre-stabilizing linear state-feedback controller is implemented from the solution of the DARE, and then the total control input is obtained as a perturbation of the feedback control law from the solution of a convex QP.
The pre-stabilizing controller ensures that the QP is well conditioned and promotes a highly accurate global solution, which in turn ensures that the gradients calculated in the backwards pass are accurate.
Additionally, an augmented Lagrangian penalty method is used to enforce constraints on state and control input.
This method still enforces hard constraints as for a sufﬁciently large penalty term the constraints are strictly enforced, but ensures that the MPC problem is feasible throughout the training process.
These contributions are in contrast to (Amos et al., 2018) which did not consider state constraints, and implemented a differential dynamic programming (Tassa et al., 2014) method to solve the MPC optimization for which convergence could not be guaranteed.
The framework is implemented on a set of second order mass-spring-damper systems and a vehicle platooning model, where it is demonstrated that the inﬁnite horizon cost can be learned and the hard constraints can be guaranteed using a short ﬁnite prediction horizon.
Notation I n := n × n identity matrix.
O m × n := m × n matrix of zeros.
n := a vector of n zeros.n := a vector of n ones.
All inequalities ≤ and ≥ are considered element-wise in the context of vectors.
ρ ( A ) := largest absolute eigenvalue of given matrix A .
vec : R m × n (cid:55)→ R mn is deﬁned as vec ([ c· · · c n ]) :
= ( c 1 · · · c n ) , i.e. the columns of a matrix stacked into a vector.
For a matrix
A ∈ R m × n , the V m,n ∈ R mn × mn permutation matrix is implicitly deﬁned by V m,n vec A := vec A (cid:62) .
The Kronecker product, ⊗ , is deﬁned as in (Magnus & Neudecker, 1999, pp. 440).
D IFFERENTIABLE MPC
Linear quadratic MPC
This paper considers linear time invariant systems of the form
[Equation]
where x t ∈ R n is the system state, u t ∈ R m is the control input, A ∈ R n × n is the state transition matrix, B ∈ R n × m is the input matrix, t ∈ R is the time, and dt ∈ R is the timestep (assumed constant).
The control problem for this system is to determine the sequence of values of u t that achieve a desired level of performance (e.g. stability, frequency response, etc...), and when the system is subject to hard constraints on control input, u t ∈ U , and state, x t ∈ X , (or a combination of both), a well studied framework for controller synthesis is MPC.
The principle of MPC is that the system’s control input and state are optimized over a ﬁnite prediction horizon, then the ﬁrst element of the obtained control sequence is implemented at the current time step and the process is repeated ad inﬁnitum .
For linear MPC it is common to use a quadratic stage cost and box constraints on state and control ( x ≤ x k ≤ x and u ≤ u k ≤ u where u ≤ 0 ≤ u ), so that at each time index t the vector of optimized control variables ˆ u (cid:63) is obtained from ˆ u (cid:63) 0:
N = argmin
ˆ u
[Equation]
s.t.
[Equation]
and
[Equation]
where ˆ u 0: N is the predicted control trajectory, ˆ x is the predicted state trajectory,
R ∈ R m × m (cid:23)represents the stage cost on control input, Q ∈ R n × n (cid:23) 0 represents the stage cost on state,
Q N ∈ R n × n (cid:23) 0 represents the terminal cost on state, N ∈ N is the prediction horizon, r k ∈ R m are slack variables for the control constraint, s k ∈ R n are slack variables for the state constraint, and k u ∈ R > 0 and k x ∈ R > 0 represent the cost of control and state constraint violations.
The variables s and r enforce the box constraints on state and control using the augmented Lagrangian method (Nocedal & Wright, 2006, §17.2), and it can be shown that for sufﬁciently high k x and k u the constraints x ≤ x k ≤ x and u ≤ u k ≤ u can be exactly guaranteed (Kerrigan & Maciejowski, 2000)
(i.e. s = r = 0 ).
The beneﬁt of this approach is that it ensures that the MPC optimization is feasible at each iteration of the learning process, whilst still ensuring that the constraints are ‘hard’.
To close the MPC control loop, at each timestep, t , the ﬁrst element of the optimized control sequence, ˆ u (cid:63) 0 , is implemented as u t .
Pre-stabilised MPC If the control input is decomposed into u t =
Kx t + δu t , where K ∈ R m × n is a linear state-feedback matrix and δu t is a perturbation to the feedback control, system (1) becomes
[Equation]
and problem (2) becomes
[Equation]
and
[Equation]
so that ˆ u (cid:63)
=
Kx t
+ δ ˆ u (cid:63) 0 is implemented as u t .
Using this decomposition, system (3) controlled with the solution of (4) is precisely equal to system (1) controlled with the solution of (2), but problem (4) is preferable from a computational standpoint if A is open-loop unstable (i.e. ρ ( A ) > 1 ) and N is ‘large’, as this can lead to poor conditioning of the matrices deﬁned in Appendix A.
This is important in the context of differentiable MPC, as if A is being learned then there may be no bounds on its eigenvalues at any given iteration.
MPC derivative.
Problems (2) and (4) can be rearranged into the QP form (see Appendix A for details)
[Equation]
When z (cid:63) is uniquely deﬁned by (5), it can also be considered as the solution of an implicit function deﬁned by the Karush-Kuhn-Tucker (KKT) conditions, and in Amos & Kolter (2017) it was demonstrated that it is possible to differentiate through this function to obtain the derivatives of z (cid:63) with respect to the parameters H , q , l , M , and u .
The MPC controller can then be used as a layer in a neural network, and backpropagation can be used to determine the derivatives of an imitation cost function with respect to the MPC parameters Q , R , A , B , u , u , x , x , k x and k u .
Imitation Learning.
A possible use case of the derivative of a model predictive controller is imitation learning, where a subset of { cost function, system dynamics, constraints } are learned from observations of a system being controlled by an ‘expert’.
Imitation learning can be performed by minimizing the loss
[Equation]
where u t is the measured control input,
ˆ u (cid:63) 0: N ( x t ) is the full MPC solution, and β ≥ 0 is a hyperparameter.
It is assumed that both the learning algorithm and MPC controller have completely precise measurements of both the state and control input.
The ﬁrst term of (6) is the control imitation loss, and the second term penalises the one-step ahead prediction error ˆ w t =
Ax t + Bu t − x t + dt .
In practice, the prediction error loss might not be needed for the MPC to be learned correctly, however its use can be instrumental for stability, as discussed in the next section.
T ERMINAL COST FOR INFINITE HORIZON
Terminal cost.
The inﬁnite-horizon discrete-time Linear Quadratic Regulator (LQR, Kalman, 2001) is given with state feedback gain
[Equation]
) where P is obtained as a solution of the DARE
[Equation]
The principle of the approach presented in this paper is the MPC controller (2,4) is implemented with
Q N = P .
Proposition 1 summarises the relevant properties of the proposed MPC, building on classic
MPC results from Scokaert & Rawlings (1998); Limon et al.
Proposition 1.
Consider the MPC problem (4) with Q N = P , where P and K solve (7-8).
Deﬁne
V (cid:63)N ( x ) as the optimal objective in (4) with x t = x .
Denote the optimal stage cost with x t = x as (cid:96) ( x, ˆ u (cid:63) 0 ( x ))
=
x (cid:62) Qx + ˆ u (cid:63)
( x ) (cid:62) R ˆ u (cid:63) 0 ( x ) .
Then, for the closed-loop system, it follows that:
For any ¯ N ≥ 1 , there exists a closed and bounded set, Ω ¯ N , such that
, if x 0 ∈ Ω ¯ N and ˆ w t = 0 , ∀ t ≥ 0 , then the MPC solution is inﬁnite-horizon optimal for any N ≥ ¯ N . 2.
There exist positive scalars d , α , such that, for any N ≥ 1 , if ˆ w t = 0 , ∀ t ≥ [-
-]0 then the MPC constraints are feasible, ∀ t ≥ 0 , and the origin is asymptotically stable ∀ x 0 ∈ Γ N , with
[Equation]
There exist a scalar, µ ≥ 0 , such that, for any N ≥ 1 the MPC constraints are robustly feasible, ∀ t ≥ 0 , and the system is Input-to-State Stable (ISS) ∀ x 0 ∈ Γ N given an additive model error, ˆ w , such that: [Equation]
In other words:
[Equation]
for some strictly increasing, bounded function, σ ( · ) , with σ (0) = 0 .
The QP matrices, H , M and the vector q , in (5), have ﬁnite norms for any N ≥ 1 .
Proof.
Proof of Proposition 1 is given in Appendix C.
Implications.
Proposition 1 has some important implications.
First, point 1 implies that there exists a state-dependant ﬁnite horizon length, ¯ N , which is sufﬁcient to make the MPC problem inﬁnite-horizon optimal.
This ¯ N can be upper bounded for a closed and bounded set of feasible states, Ω ¯ N .
Scokaert & Rawlings (1998) proposed an iterative search that increases the horizon until optimality is veriﬁed; a similar algorithm is discussed in Appendix E where learning is completed with a large horizon and then iteratively reduced afterwards, although it is not implemented in this paper.
Point 2,3 state that MPC that can provide stability and constraints satisfaction, hence safety , if the model error is small.
This also applies to small errors in the QP solution.
Finally, point 4 states that the QP matrices have ﬁnite norm when the system dynamics are pre-stabilised using the LQR gain 2 , so the MPC problem is well conditioned and can be solved reliably to high accuracy, even over long horizons.
If the open-loop system is unstable then the terms of the matrices in Appendix A for the standard form are unbounded, so the QP solution may be poorly conditioned and the result inaccurate for long horizons.
This can in turn invalidate the results of Amos & Kolter (2017) which assumes that the KKT conditions are exactly satisﬁed in order to compute its gradients.
DARE Derivative.
In order to implement Q N = P in a differentiable imitation learning framework such as that presented in Section 2, the solution of the DARE is differentiated as follows.
Proposition 2.
Let P be the stabilizing solution of (8)
[Equation]
and assume that Z − 1 1 and ( R + B (cid:62) PB ) −exist, then the Jacobians of the implicit function deﬁned by (8) are given by
[Equation]
where Z 1 , . .
, Z 5 are deﬁned by
[Equation]
Z
[Equation]
and M 1 , M 2 , M 3 are deﬁned by
[Equation]
Proof.
The proof of Proposition 2 is given in Appendix D.
The sensitivity of the solution of the DARE has been investigated in the context of robustness to perturbations in the input matrices, e.g. Sun (1998); Konstantinov et al.
(1993), and the analytical derivative of the continuous-time algebraic Riccati equation was derived in Brewer (1977) by differentiating the exponential of the Hamiltonian matrix, but to the best of the authors’ knowledge this is the ﬁrst presentation of an analytic derivative of the DARE using the differential calculus approach of
Magnus & Neudecker (1999).
In: M \ S , N > 0 , β ≥ 0 , N epochs > 0 .
Out:
S for i = 0 ...
N epochs do Forward Pass ( K, P ) ← DARE (7-8) solution Q
T ← P ˆ u (cid:63) 0: N ← MPC QP (3-5) solution
L ← Imitation loss (6) Backward Pass
Differentiate loss (6) Differentiate MPC QP solution, ˆ u (cid:63) 0: N , using Appendix B Differentiate DARE, ( P, K ) , using Proposition 2 Update step
S ← Gradient-based step
Algorithm overview Algorithm 1 presents the overall procedure for learning a subset, S , of the MPC controller parameters, M = { A, B, Q, R, x, x, u, u, k u , k x } , with the key steps of the forwards and backwards pass of a gradientbased optimization method.
In each forward pass the MPC terminal cost matrix, Q N , and the prestabilizing controller, K , are set from the solution of the DARE, then the DARE and MPC QP solutions are differentiated in the backward pass so to obtain the gradients.
Note that the horizon, N , is not differentiable, and that learning the entire set
M simultaneously is challenging in general.
N UMERICAL E XPERIMENTS
In this section the performance of the algorithm was demonstrated through numerical experiments in two test cases: ﬁrstly on a set of second order mass-spring-damper models to provide a performance baseline in an easily interpretable setting, and then on vehicle platooning problem to investigate a higher-dimensional real-world application.
4. M ASS -S PRING -D AMPER
Model & Expert Expert data was generated using a mass-spring-damper model parameterized by a mass, m ∈ R > 0 , damping coefﬁcient, c ∈ R , stiffness, k ∈ R , and timestep dt ∈ R > 0 , where
[Equation]
so that x t ∈ R 2 is the position and velocity of the mass, and the u t ∈ R is a force applied to the mass.
[Table]
Seven models were created with m = 1 , k = 1 , and dt = 0 . 2 , and c was varied as shown in Table 1 to affect the open-loop stability of the models ( c > 0 ⇒ stable, c < 0 ⇒ unstable).
The expert data was then generated by simulating each of the systems the initial condition x 0 = (0 , 3) in closed-loop with an inﬁnite-horizon MPC controller (i.e. the horizon was increased until the open-loop state predictions matched the closed-loop response), using Q = diag ([1 , 1]) , R = 2 , ( u, u ) =
( −∞ , 0 . 5) , x = ( − 1 , −∞ ) , and x ) =
(1 , ∞ ) .
The constraint set was chosen so that the constraints on both state and control input were strongly active at the solution whilst ensuring that the expert MPC optimization was feasbile.
The values k u
= k x = 100 were found to be sufﬁcient to enforce the hard constraints and were used for all experiments.
It is important to note that the approach of (Amos et al., 2018) cannot be used reliably for even this simple example as it does not consider state constraints, and when hard constraints are added to the method it fails in general because the optimization problem has become infeasible in the forwards pass at some time t .
Learning The learner and expert shared all system and controller information apart from the state transition matrix A , which was learned, and the MPC horizon length, which was implemented as each of N ∈ { 2 , 3 , 6 } in three separate experiments.
A was initialized with the correct state transition matrix plus a uniformly distributed pseudo-random perturbation in the interval [ − 0 .
5 , 0 . 5] added to each element.
The learner was supplied with the ﬁrst 50 elements of the closed loop state trajectory and corresponding controls as a batch of inputs, and was trained to minimize the imitation loss (6) with β = 0 , i.e. the state dynamics were learned using predicted control trajectories only , and the state transitions are not made available to the learner (this is the same approach used in Amos et al., 2018).
The experiments were implemented in Pytorch 1.2.0 using the built-in Adam optimizer (Kingma & Ba, 2014) for 1000 steps using default parameters.
The MPC optimization problems were solved for the ‘expert’ and ‘learner’ using OSQP (Stellato et al., 2017) with settings (eps_abs=1E-10, eps_rel=1E-10, eps_rim_inf=1E-10, eps_dual_inf=1E-10).
[Figure]
The model loss was not used as part of the training process, and shown only to indicate whether the model is converging correctly.
Results Figure 4 shows the imitation and model loss at each of the 1000 optimization iterations for each of the tested horizon lengths.
It can be seen that for all of the generated systems the imitation loss converges to a low value, although this is a local minimum in general.
In most cases, the learned model converges to a close approximation of the real model, although as the problem is non-convex this cannot be guaranteed, and it is also shown that there are some cases in which the model does not converge correctly.
This occurred exclusively for N = 2 , where neither system 4 nor system 2 converge to the correct dynamics.
Additionally, it can be seen that both the imitation loss and model loss converge faster as the prediction horizon is increased.
This suggests that a longer learning horizon improves the learning capabilities of the methods, but there is not sufﬁcient data to demonstrate this conclusively.
t −  x
N =  t − N =  t −
N =
To test generalization performance, each of the systems was re-initialized with initial condition x 0 = (0 . 5 , 2) and simulated in closed loop using the learned controller for each horizon length.
The results are compared in Figure 2 against the same systems controlled with an inﬁnite horizon MPC controller.
The primary observation is that as the learned MPC horizon is increased to N = 6 , the closed loop trajectories converge to expert trajectories, indicating that the inﬁnite horizon cost has been learned (when using the inﬁnite horizon cost with no model mismatch or disturbance, the predicted MPC trajectory is exactly the same as the closed loop trajectory), and that the state constraints are guaranteed for N ≥ 4 .
Furthermore, it can be seen that the learned controllers are stabilizing, even for the shortest horizon and the most unstable open-loop systems.
This is also the case for systems 2 and 4 where the incorrect dynamics were learned, although in this case the state constraints are not guaranteed for N = 2 .
4. V EHICLE P LATOONING
[Figure]
Model & Expert Vehicle platoon control is a problem that has been studied using control theory (e.g. Zheng et al.
, but here we demonstrate that a safe, stabilizing controller can be learned from examples of vehicles driving in formation.
Figure 3 shows an illustration of a platoon of n v vehicles for which the objective is to stabilize the relative longitudinal positions of each vehicle to the steady-state conditions
y
i − y i −
→
y ss and ˙
y
i − ˙
y i − 1 → 0 ∀ i , subject to the hard constraint that relative position of the vehicles is never lower than a safe threshold y
i − y
i − 1 ≥ y ∀ i , and that the vehicles’ ability to brake and accelerate is constrained by b ≤ ¨ y i ≤ a ∀ i where b < 0 < a (note that only the relative positions and velocities of the vehicles is considered, as the global position and velocity of the platoon can be controlled separately by adding an equal perturbation to each element of ¨ y ).
In appendix F it is shown that this can be modelled as a discrete time LTI system.
Expert data was generated from the model with n v = 10 vehicles so that x t ∈ R 18 and u t ∈ R 10 .
20 instances were generated using random feasible initial conditions with y ss = 30 m and y = 10 m, and then simulated for 20 s in time intervals of dt = 0 .
7 s with an inﬁnite-horizon MPC controller, using Q = {+I n and R = 
+}
I m .
Learning The learner and expert shared all system and controller information apart from the cost matrices Q and R , which were learned, and the MPC horizon length, which was implemented as each of N ∈ { 5 , 10 , 15 , 20 } in four separate experiments.
The matrices Q and R were initialized as completely random diagonal matrices with each element uniformly distributed in the interval [0 , 3] , and the diagonal structure was maintained through training.
500 training iterations were used; otherwise the learning process (loss function, learning rate, etc...) was the same as in Section 4.1.
Results Figure 4 shows the imitation and cost function losses at each of the 500 optimization iterations for each of the tested horizon lengths and initial conditions.
As with the mass-springdamper experiments, it is suggested that a longer prediction horizon improves training as the imitation loss generally converges to a lower value for the examples with N ∈ { 15 , 20 } , but only convergence to a local minimum is achieved in general.
The cost error also does not converge in general (although better convergence is observed again for the longer horizon lengths), however for this learning problem there is a manifold of matrices Q and R with the same minimizing argument, so divergence of the cost error does not necessarily indicate that the learned cost function is ‘incorrect’. Furthermore, in this case the model is known exactly, so the closed-loop inﬁnite-horizon properties can be obtained even without the correct cost function. [Figure] [Figure]
Iterations
[Figure]
[Figure]
C o s t
F un c t i o n
L o ss
N = N =
N =
N =
Figure 5 shows the model simulated from the same initial condition in closed loop using a learned controller for each of the horizon lengths, together with the error between the MPC state predictions and ensuing closed-loop behaviour.
All of the controllers are observed to successfully satisfy the hard constraints on vehicle separation, and all converge to the correct steady-state vehicle separation.
The differences between the prediction capabilities of the controllers is highlighted by the state prediction errors, and it can be seen that for N = 20 the state predictions match the ensuing behaviour, indicating that the inﬁnite horizon cost is being used and that closed-loop stability is guaranteed, even without the use of a terminal constraint set.
It is also demonstrated for N < 20 that the largest errors occur from predictions made at times when the state constraints are active, suggesting that these controllers deviate from their predictions to satisfy the constraints at later intervals.
4. L IMITATIONS
The above approach is limited in scope to LTI systems, and a more comprehensive solution would cover linear time varying systems (for which the MPC is still obtained from the solution of a QP).
In this case the inﬁnite horizon cost cannot be obtained from the solution of the DARE, and the extension of the methods presented in this paper to time varying or non-linear models is non-trivial (see Appendix G for further discussion).
The derivative of the DARE in Proposition 2 involves multiple Kronecker products and matrix inversions (including an n 2 × n 2 matrix) that do not scale well to large state and control dimensions, although the dynamics of physical systems can usually be reasonably approximated with only a few tens of variables, so this may not become an issue in practice.
The algorithm also requires a stabilizing solution of the DARE to exist; theories for the existence of stabilizing solutions are non-trivial (e.g. Ran & Vreugdenhil, 1988), and it is not immediately obvious how to enforce their existence throughout the training process (stabilizibility can be encouraged using the one-step ahead term in (6)). t (s) − S t a t e P r e d i c t i o n
E rr o r ( m )    y ( m )
N =  t (s) −
N =  t (s) − N =  t (s) − N = A PPENDICES A MPC QUADRATIC PROGRAM
Problem (2) is equivalent to
[Equation]

[Equation]

[Equation]

[Equation]
  ∞ u ∞∞
[Equation]
∞ 
 , where z = (cid:34) ˆ urs (cid:35) ,
R =  
R . . .
R 
 ,
Q = 

Q . . .
Q
Q N   , Φ =  
A ...
A N 
 , Ψ =  
B ...
A N − 1 B · · ·
B 
 , x =   x ... x   , x =   x ... x   , u =   u ... u 
 , u =   u ... u   , are of conformal dimensions.
Using the above, problem (4) is then equivalent to z
[Equation]

[Equation]
K
[Equation]
s.t.

[Equation]

[Equation]
−∞ 0 Nm
[Equation]
−∞
[Equation]
∞
[Equation]
 , where now z = (cid:34)
δ ˆ urs (cid:35) ,
Φ =   ( A + BK ) .
( A + BK ) N 
 and Ψ =  
B ...
( A + BK ) N − 1 B · · ·
B 
 , and ˆ Q = (cid:20) 0 n × n
Q (cid:21) ,
[Equation] [Equation]


[Equation]
are of conformal dimensions.
B OSQP DERIVATIVES
OSQP solves quadratic programs of the form (5), and returns values for z , y , and s that satisfy
[Equation]
(Stellato et al., 2017, §2), where C is the set { s : l ≤ s ≤ u } , and N C is the normal cone of C .
The values of y that are returned by the solver can be used to determine whether the constraints are strongly active at the solution, where y i = 0 indicates that the constraints l i ≤ M i z and M i z ≤ u i are inactive, y i > 0 indicates that M i z ≤ u
i is strongly active, and y i < 0 indicates that l i ≤ M i z is strongly active.
The solution can therefore be completely characterised by the KKT system
[Equation]
where U = { i : y i > 0 } and L
=
{ i : y
i < 0 } , and the notation M S indicates a matrix consisting of the i ∈ S columns of given matrix M , and v S indicates a vector consisting of the i ∈ S elements of given vector v .
Equation (10) can then be differentiated using the techniques detailed in (Amos &
Kolter, 2017, §3)
C P ROOF OF P ROPOSITION
Proof.
( Proposition 1 )
The ﬁrst point follows from (Scokaert & Rawlings, 1998).
The next two points of Proposition 1 stem from the results in (Limon et al., 2003; 2009).
In particular, the closedloop is Lipschitz since the model is linear and the controller is the solution of a strictly convex QP.
Moreover, the LQR provides a contractive terminal set.
The ﬁnal point follows from the fact that ( A + BK ) N has eigenvalues in the unit circle, ∀ N ≥ 1 .
Proof of point 4 is concluded by inspection of the QP matrices (Appendix A) and by application of Theorem 5.6.12, page 298 of Horn & Johnson (2012) which states that, given a bound, ρ , on the spectral radius, then there exists a matrix norm which is also less than ρ .
D P ROOF OF P ROPOSITION
Proof.
( Proposition 2 )
If a stabilizing solution ( ρ ( A + BK ) ≤ 1 ) to (8) exists, it is unique (Ionescu & Weiss, 1992, Proposition 1), and the DARE can therefore be considered an implicit function of A ,
B , Q , and R .
Using the assumption that ( R + B (cid:62) PB ) − 1 exists, it can be concluded that Z 1 , . . .
, Zand M 1 , M 2 , M 3 exist (the Kronecker product and matrix addition, subtraction, and multiplication always exist).
Equation (8) can be given by
[Equation]
which is differentiable, and M 1 , M 2 , M 3 are also differentiable.
Differentials are taken for (11) and each of M 1 , M 2 , M 3 as
[Equation]
then these can be combined using the differential chain rule (Magnus & Neudecker, 1999, Theorem 18.2) to obtain
[Equation]
The Jacobians, as deﬁned in Proposition 2, therefore exist if Z − 1 1 exists.
E V ERIFICATION AND REDUCTION OF THE PREDICTION HORIZON
In: N > 0 , X 0 ⊆ X , M , ( P, K ) from (7-8), (cid:15) > 0 , n s > 0 , η ∈ (0 , 1) .
Out: ¯ N,
X X ← X 0 while X ⊃ ∅ do ¯ N ← N while ¯ N > 0 do
X sample ← n s uniform state samples, s.t.:
x ∈ X δ ˆ u (cid:63) ← Solution of MPC QP (3-5), ∀ x ∈ X sample if (cid:107) δ ˆ u
(cid:63)k ( x ) (cid:107) ≤ (cid:15), ∀ k ≥ ¯ N, ∀ x ∈ X sample then return TRUE
¯
N ← ¯ N −
X ← η X
Procedure failed N ← N +
Go to Algorithm
A method is proposed for the reduction of the MPC prediction horizon after imitation learning.
The idea is to be able to reproduce the inﬁnite-horizon optimal
MPC up to a tolerance (cid:15) with high probability.
Do do so, we check that, for a candidate horizon ¯ N , the MPC action deltas, δ ˆ u
(cid:63)k , satisfy (cid:107) δ ˆ u
(cid:63)k (cid:107) ≤ (cid:15) , for all k ≥ ¯ N .
This means that the optimal action is equal to the LQR up to a tolerance (cid:15) .
In order to provide a high probability guarantee of this condition, we propose the use of a probabilistic veriﬁcation approach, similar to Bobiti (2017).
This is described in Algorithm 2.
In particular, the condition is checked on a high number, n s , of initial states.
These states are sampled uniformly from a set of interest X , which can be either the state constraints X or an estimate of the region of attraction, Γ N .
If veriﬁed, this set is a region of attraction for the system with high probability.
The relationship between the number of samples and the veriﬁcation probability is discussed in (Bobiti, 2017, Chapter 5).
The algorithm also checks whether the inﬁnite horizon condition has been reached for the N used during training.
Finally, a line search for a suitable X is proposed using a scaling factor η ∈ (0 , 1) .
In particular, the initial set is downscaled until either an horizon is found or the set becomes empty.
In latter case the search fails and the procedure returns to the training algorithm with an increased N .
Noticeably, the proposed algorithm does not require to explicitly compute the terminal set in which the LQR is invariant and it could be used also for non-linear MPC if an inﬁnite-horizon (or a stabilising) terminal controller is available.
F P LATOON M ODEL D ERIVATION
The problem described in Section 4.2 can be decomposed into the regulation problem
[Equation]
subject to the constraints
[Equation]
and 1 b ≤ ¨ y ≤ 1 a.
If each vehicle is modelled as a mass then a continuous-time LTI state space model can be formed as  
[Equation]

[Equation]

[Equation]

[Equation]

[Equation]
z n − 1 ˙ z− ˙ z 1 .
˙
z n − ˙ z n −   (cid:124) (cid:123)(cid:122) (cid:125) x +   − ...
− 1  (cid:124) (cid:123)(cid:122) (cid:125)
B c   ¨ y 1 .
¨ y n 
 (cid:124) (cid:123)(cid:122) (cid:125) u , (12) which can then be given as ˙ x = (cid:20)
I  (cid:21) x
+ (cid:20)ˆ B (cid:21) u.
If it is assumed that the control input is constant between sampling intervals t and t + d t , then this can be given in discrete time as
[Equation]
13) where x t ∈ R 2( n v − 1) , and u ∈ R n v and are subject to the constraints
[Equation]
G N ONLINEAR MODELS
As discussed in the main paper, our approach is currently limited to Linear Time Invariant (LTI) systems.
In general, conditions for inﬁnite-horizon optimality of systems that are not LTI are nontrivial.
Some of the results on MPC stability could however be maintained, for example in the case when the LQR value function, x (cid:62) Px , is a local control Lyapunov function (Khalil, 2014; Mayne et al., 2000).
In this case, the stability and intrinsic robustness results are maintained (see Limon et al., 2003; 2009).
For these system, it would be possible to use our method, for instance in combination with Amos et al.
(2018), to provide a stable Non-linear MPC.
This is however a big assumptions for systems that are very non-linear.
Assessing this LQR controllability condition could be done, for instance, by training a local linear model around the target equilibrium (origin) and then checking whether the DARE is solvable.
This should be performed before starting the imitation learning.
We leave the study of more general systems to future work.