Lecture-IntroductiontoModelPredictiveControl
Lecture-IntroductiontoModelPredictiveControl
net/publication/259124422
CITATIONS READS
32 673
2 authors:
All content following this page was uploaded by Donal O'Donovan on 15 October 2021.
Presented by:
Donal O’Donovan
Munster Technological University
([email protected])
What Can I Expect From Today’s Lecture?
• Gain an intuitive feel for what Model Predictive Control (MPC) is
• Understand the benefits of MPC
• Where and when MPC may be applied
• Appreciate the implementation requirements
• Why isn’t MPC everywhere?
• Is there an alternative MPC that is accessible for ‘smaller’ plants?
• Predictive Functional Control
Unit 1: Distributed Control System (PID) Unit 2: Distributed Control System (PID)
Basic Dynamic Control
(Every second)
FC TC FC TC
PC LC PC LC
Source: Wikipedia
Source: [8] Donal O'Donovan - Munster Technological University
Chess Game
(𝒌 + 𝑯)
Analogy
Etc. Horizon (H) (𝒌)
(𝒌 + 𝑯 + 𝟏) Horizon (H) (𝒌 + 𝟏)
Opponent’s Move My Move
(Your Move) (Plant) (Controller)
(My Move) (Your Move) Plant
Controller Plant 𝑦(𝑘 + 2)
⋮
⋮
⋮
⋮
𝑘+𝑁 𝑘+𝑖 (My Move)
⋮ Controller
⋮
⋮
⋮ u(𝑘)
⋮
⋮ ⋮ Opponent’s
My Move
⋮ Move
⋮ ⋮ (Plant) (Controller)
⋮ ⋮ ⋮
⋮
⋮
⋮ ⋮ ⋮ (My Move)
⋮ ⋮ Controller
⋮ ⋮
⋮
⋮
⋮
⋮
Convergence 10 s faster
Peak reduced by 50%
Donal O'Donovan - Munster Technological University
Example: Temperature setpoint tracking in Air Conditioner
• The goal of the control equation is to estimate the projection weighting's of the functional basis
set on to the future 𝑢(𝑛 + 𝑖)
• i.e. 𝑀𝑉(0) is the weighted sum of the functional basis set
• For example:
• If the setpoint is a ramp, 𝑢(𝑛) should contain step and ramp basis function components
• The step is needed to deal with constant offsets
Donal O'Donovan - Munster Technological University
Future MV Hypotheses - Basis Functions
• The weighting coefficients of the functional basis are computed as the unknown variables of
the control equation
• The 𝑢(𝑛 + 𝑖) is then, by extension, projected on to a polynomial basis: 𝑃𝑗 𝑛 , where 𝜇𝑗 are
the unknown weighting coefficients of the polynomial 𝑃, which is a function of the future
time 𝑖, i.e.
𝑴𝑽 𝒏 + 𝒊 = 𝝁𝒋 𝑷𝒋 𝒊 with 𝒋 = 𝟏, 𝟐, … , 𝒋𝒎𝒂𝒙
𝒋
• For instance,
𝑴𝑽(𝒏 + 𝒊) = 𝝁𝟏 𝟏 + 𝝁𝟐 𝒊
where 𝑷𝟏 = 𝟏 i.e. step
𝑷𝟐 = 𝒊 i.e. a ramp
𝝁𝒋 are the unknown coefficients to be determined
Donal O'Donovan - Munster Technological University
Future MV Hypotheses - Basis Functions
• If two basis functions are selected, then there should be two equations and, at least,
two coincidence points
• This implies that it is possible to structure the future 𝑢(𝑛 + 𝑖) for 𝑛 = 1, … , ℎ using a
limited number of weighted functions of the functional basis set
• Also, it is not necessary to compute all future 𝑢(𝑛 + 𝑖) over the coincidence horizon
• For instance, in the step and ramp case, it is sufficient to choose only two coincidence
points ℎ1 and ℎ2 over the horizon
• If the functional basis is limited to ramp, step, and parabola, a simple algebraic ‘one-
step’ solver may be used – no iterations required
• In summary, the eigenfunctions, chosen as the functional basis set for the MV, dictate
the ability of track of the controller to track the setpoint
• i.e. It dictates the ‘accuracy’ of the closed-loop control
Donal O'Donovan - Munster Technological University
Implicit Integrator and Time Delay Compensation
• A steady-state analysis of the basic PFC equation shows that there is no
steady-state error present
• It may be easily shown that there is an implicit integrator present in the controller
• If the process possesses a pure time delay 𝜃 the reference trajectory is simply
initialized, using the predicted value of the process output 𝑦ො𝑝 at the future
time 𝑛 + 𝜃 to compensate
[9] H. Bouhenchir, M. Cabassud, M.V. Le Lannc, Computers and Chemical Engineering (2006)
Setpoint Profile of Chemical Batch Reactor
• Test the robustness of the control system when there are
dynamic changeovers, due to the heat release, during
constant setpoint stage
• The reactor was fed with a solution of hydroxide sodium
(NaOH/water = 57.14 g/485.71 g) at temperature of 22 ℃
• The set-point temperature profile is composed of:
• Stage 1: Heating from 22 to 45 ℃ for 1500 s (0.92 ℃ / min)
• Stage 2: Constant temperature at 45 ℃ during 3500 s in
which solution of hydrochloric acid (HCl/water = 44.03
g/74.97 g) is fed during 2000 s
• Stage 3: Cooling from 45 to 30 ◦C during 1500 s (−0.6 ℃/
min) and fourth stage: maintain at 30 ◦C during 500 s
• In order to eliminate discontinuities resulting from change of
stage, the set-point profile was filtered by a procedure called
“docking procedure” which transforms the temperature
profile as shown
Results for PFC-controlled Chemical Batch Reactor
• Set point (Tr,set)
• Adjusted set point (Td,set)
• Reaction mixture (Tr) temperatures (left y-axis)
• The manipulated variable (β) computed by the slave controllers (right y-axis)
• The manipulated variable (β) denotes the fraction (compared to its
maximum) of the electrical power value or the mono-fluid flow rate fraction
dispatched to one of the two PHE
Results with Model Mismatch
The results show the performance of cascade control when there are differences in the dynamics of the different heating/
cooling system elements.
Results with No Model Mismatch