Time-Varying MPC

When to Use Time-Varying MPC

适应不断变化的操作条件,adaptive MPC supports updating the prediction model and its associated nominal conditions at each control interval. However, the updated model and conditions remain constant over the prediction horizon. If you can predict how the plant and nominal conditions vary in the future, you can use time-varying MPC to specify a model that changes over the prediction horizon. Such a linear time-varying (LTV) model is useful when controlling periodic systems or nonlinear systems that are linearized around a time-varying nominal trajectory.

To use time-varying MPC, specify arrays for thePlantandNominalinput arguments ofmpcmoveAdaptive. For an example of time-varying MPC, seeTime-Varying MPC Control of a Time-Varying Plant.

Time-Varying Prediction Models

Consider the LTV prediction model

$\begin{matrix} x (k + 1) = A (k) x (k) + B_{u} (k) u (k) + B_{v} (k) v (k) \\ y (k) = C (k) x (k) + D_{v} (k) v (k) \end{matrix}$

whereA,B_u,B_v,C, andD是discrete-time state-space matrices that can vary with time. The other model parameters are:

k— Current control interval time index
x— Plant model states
u— Manipulated variables
v— Measured disturbance inputs
y— Measured and unmeasured plant outputs

Since time-varying MPC extends adaptive MPC, the plant model requirements are the same; that is, for each model in thePlantarray:

Sample time (Ts) is constant and identical to the MPC controller sample time.
Any time delays are absorbed as discrete states.
The input and output signal configuration remains constant.
There is no direct feed-through from the manipulated variables to the plant outputs.

For more information, seePlant Model.

The prediction of future trajectories forpsteps into the future, wherepis the prediction horizon, is the same as for the adaptive MPC case:

$[\begin{matrix} y (1) \\ ⋮ \\ y (p) \end{matrix}] = S_{x} x (0) + S_{u 1} u (- 1) + S_{u} [\begin{matrix} Δ u (0) \\ ⋮ \\ Δ u (p - 1) \end{matrix}] + H_{v} [\begin{matrix} v (0) \\ ⋮ \\ v (p) \end{matrix}]$

However, for an LTV prediction model, the matricesS_x,S_u1,S_u, andH_vare:

$\begin{matrix} S_{x} = [\begin{array}{l} C (1) A (0) \\ C (2) A (1) A (0) \\ ⋮ \\ C (p) \prod_{i = 0}^{p - 1} A (i) \end{array}] \\ S_{u 1} = [\begin{array}{l} C (1) B_{u} (0) \\ C (2) [B_{u} (1) + A (1) B_{u} (0)] \\ ⋮ \\ C (p) \sum_{k = 0}^{p - 1} [(\prod_{i = k + 1}^{p - 1} A (i)) B_{u} (k)] \end{array}] \\ S_{u} = [\begin{array}{l} 0 & 0 & \dots & 0 \\ S_{u 1} & C (2) B_{u} (1) & 0 & \dots & 0 \\ ⋮ \\ C (p) \sum_{k = 1}^{p - 1} [(\prod_{i = k + 1}^{p - 1} A (i)) B_{u} (k)] & \dots & \dots & C (p) B_{u} (p - 1) \end{array}] \\ H_{v} = [\begin{array}{l} C (1) B_{v} (0) & D_{v} (1) & 0 & \dots & 0 \\ C (2) A (1) B_{v} (0) & C (2) B_{v} (1) & D_{v} (2) & \dots & 0 \\ ⋮ & ⋮ & ⋮ & ⋮ \\ C (p) (\prod_{i = 1}^{p - 1} A (i)) B_{v} (0) & \dots & \dots & C (p) B_{v} (p - 1) & D_{v} (p) \end{array}] \end{matrix}$

where $\prod_{i = k_{1}}^{k_{2}} A (i) ≜ A (k_{2}) A (k_{2} - 1) \dots A (k_{1})$ if $k_{2} \geq k_{1}$ , orIotherwise.

For more information on the prediction matrices for implicit MPC and adaptive MPC, seeQP Matrices.

Time-Varying Nominal Conditions

Linear models are often obtained by linearizing nonlinear dynamics around time-varying nominal trajectories. For example, consider the following LTI model, obtained by linearizing a nonlinear system at the time-varying nominal offsetsx_off,u_off,v_off, andy_off:

$\begin{matrix} x (k + 1) - x_{o f f} (k) = A (k) (x (k) - x_{o f f} (k)) + B_{u} (k) (u (k) - u_{o f f} (k)) \\ + B_{v} (k) (v (k) - v_{o f f} (k)) + Δ x_{o f f} (k) \\ y (k) - y_{o f f} (k) = C (k) (x (k) - x_{o f f} (k)) + D_{v} (k) (v (k) - v_{o f f} (k)) \end{matrix}$

If we define

$\begin{matrix} \bar{x_{o f f}} ≜ x (0), \bar{u_{o f f}} ≜ u (0) \\ \bar{v_{o f f}} ≜ v (0), \bar{y_{o f f}} ≜ y (0) \end{matrix}$

as standard nominal values that remain constant over the prediction horizon, we can transform the LTI model into the following LTV model:

$\begin{matrix} x (k + 1) - \bar{x_{o f f}} = A (k) (x (k) - \bar{x_{o f f}}) + B_{u} (k) (u (k) - \bar{u_{o f f}}) + B_{v} (k) (v (k) - \bar{v_{o f f}}) + {\bar{B}}_{v} (k) \\ y (k) - \bar{y_{o f f}} = C (k) (x (k) - \bar{x_{o f f}}) + D_{v} (k) (v (k) - \bar{v_{o f f}}) + {\bar{D}}_{v} (k) \end{matrix}$

where

$\begin{matrix} {\bar{B}}_{v} (k) ≜ Δ x_{o f f} (k) + x_{o f f} (k) - \bar{x_{o f f}} + A (k) (\bar{x_{o f f}} - x_{o f f} (k)) + B_{u} (k) (\bar{u_{o f f}} - u_{o f f} (k)) \\ + B_{v} (k) (\bar{v_{o f f}} - v_{o f f} (k)) \\ {\bar{D}}_{v} (k) ≜ y_{o f f} (k) - \bar{y_{o f f}} + C (k) (\bar{x_{o f f}} - x_{o f f} (k)) + D_{v} (k) (\bar{v_{o f f}} - v_{o f f} (k)) \end{matrix}$

If the original linearized model is already LTV, the same transformation applies.

State Estimation

As with adaptive MPC, time-varying MPC uses a time-varying Kalman filter based onA(0),B(0),C(0), andD(0) from the initial prediction step; that is, the current time at which the state is estimated. For more information, seeState Estimation.

Time-Varying MPC

When to Use Time-Varying MPC

Time-Varying Prediction Models

Time-Varying Nominal Conditions

State Estimation

See Also

Functions

Objects

Related Examples

More About