ModelPredictiveControl.jl

An open source model predictive control package for Julia.

The package depends on ControlSystemsBase.jl for the linear systems and JuMP.jl for the solving.

Installation

To install the ModelPredictiveControl package, run this command in the Julia REPL:

using Pkg; Pkg.add("ModelPredictiveControl")

Getting Started

To construct model predictive controllers (MPCs), we must first specify a plant model that is typically extracted from input-output data using system identification. The model here is linear with one input, two outputs and a large time delay in the first channel:

\mathbf{G}(s) = \frac{\mathbf{y}(s)}{\mathbf{u}(s)} = 
\begin{bmatrix}
    \frac{2e^{-20s}}{10s + 1} \\[3pt]
    \frac{10}{4s +1}
\end{bmatrix}

We first construct the plant model with a sample time $T_s = 1$ s:

using ModelPredictiveControl, ControlSystemsBase
G = [ tf( 2 , [10, 1])*delay(20)
      tf( 10, [4,  1]) ]
Ts = 1.0
model = LinModel(G, Ts)

Our goal is controlling the first output $y_1$, but the second one $y_2$ should never exceed 35:

mpc = LinMPC(model, Mwt=[1, 0], Nwt=[0.1])
mpc = setconstraint!(mpc, ymax=[Inf, 35])

The keyword arguments Mwt and Nwt are the output setpoint tracking and move suppression weights, respectively. A setpoint step change of five tests mpc controller in closed-loop. The result is displayed with Plots.jl:

using Plots
ry = [5, 0]
res = sim!(mpc, 40, ry)
plot(res, plotry=true, plotymax=true)

See the manual for more detailed examples.

Features

Legend

• [x] implemented feature
• [ ] planned feature

Model Predictive Control Features

• [x] linear and nonlinear plant models exploiting multiple dispatch
• [x] model linearization based on automatic differentiation (exact Jacobians)
• [x] supported objective function terms:
  • [x] output setpoint tracking
  • [x] move suppression
  • [x] input setpoint tracking
  • [x] terminal costs
  • [x] economic costs (economic model predictive control)
• [x] adaptive linear model predictive controller
  • [x] manual model modification
  • [x] automatic successive linearization of a nonlinear model
  • [x] objective function weights and covariance matrices modification
• [x] explicit predictive controller for problems without constraint
• [x] online-tunable soft and hard constraints on:
  • [x] output predictions
  • [x] manipulated inputs
  • [x] manipulated inputs increments
  • [x] terminal states to ensure nominal stability
• [ ] custom manipulated input constraints that are a function of the predictions
• [x] supported feedback strategy:
  • [x] state estimator (see State Estimation features)
  • [x] internal model structure with a custom stochastic model
• [x] automatic model augmentation with integrating states for offset-free tracking
• [x] support for unmeasured model outputs
• [x] feedforward action with measured disturbances that supports direct transmission
• [x] custom predictions for:
  • [x] output setpoints
  • [x] measured disturbances
• [x] easy integration with Plots.jl
• [x] optimization based on JuMP.jl:
  • [x] quickly compare multiple optimizers
  • [x] nonlinear solvers relying on automatic differentiation (exact derivative)
• [x] additional information about the optimum to ease troubleshooting
• [x] implementation that carefully limits allocations for real-time applications

State Estimation Features

• [x] supported state estimators/observers:
  • [x] steady-state Kalman filter
  • [x] Kalman filter
  • [x] Luenberger observer
  • [x] internal model structure
  • [x] extended Kalman filter
  • [x] unscented Kalman filter
  • [x] moving horizon estimator
• [x] easily estimate unmeasured disturbances by adding one or more integrators at the:
  • [x] manipulated inputs
  • [x] measured outputs
• [x] bumpless manual to automatic transfer for control with a proper initial estimate
• [ ] estimators in two possible forms:
  • [x] predictor (or delayed) form to reduce computational load
  • [ ] filter (or current) form to improve accuracy and robustness
• [x] moving horizon estimator in two formulations:
  • [x] linear plant models (quadratic optimization)
  • [x] nonlinear plant models (nonlinear optimization)
• [x] moving horizon estimator online-tunable soft and hard constraints on:
  • [x] state estimates
  • [x] process noise estimates
  • [x] sensor noise estimates