Dynamic Control Introduction

March 11, 2025, at 04:48 PM by 10.35.117.248 -

Changed lines 59-60 from:

~~import matplotlib.pyplot as plt~~

to:

Deleted line 62:

~~#constants~~

Deleted lines 63-64:

~~#Parameters~~

Deleted line 66:

~~#Manipulated variable~~

Changed lines 68-69 from:

~~#Controlled Variable~~

to:

p.STATUS = 1 # allow optimizer to change
p.DCOST = 0.1 # smooth out gas pedal movement
p.DMAX = 20 # slow down change of gas pedal

Changed lines 73-79 from:

to:

v.STATUS = 1 # add the SP to the objective
m.options.CV_TYPE = 2 # sq error
v.SP = 40 # set point
v.TR_INIT = 1 # setpoint trajectory
v.TAU = 5 # time constant of setpoint trajectory

m.options.IMODE = 6 #control

Deleted lines 80-97:

m.options.IMODE = 6 #control

#MV tuning
p.STATUS = 1 #allow optimizer to change
p.DCOST = 0.1 #smooth out gas pedal movement
p.DMAX = 20 #slow down change of gas pedal

#CV tuning
#setpoint
v.STATUS = 1 #add the SP to the objective
m.options.CV_TYPE = 2 #L2 norm
v.SP = 40 #set point
v.TR_INIT = 1 #setpoint trajectory
v.TAU = 5 #time constant of setpoint trajectory

#%% Solve

Added line 84:

import matplotlib.pyplot as plt

Changed lines 91-92 from:

plt.ylabel('velocity')
plt.xlabel('time')

to:

plt.ylabel('velocity'); plt.xlabel('time')

October 18, 2024, at 05:21 PM by 10.35.117.248 -

Changed lines 117-120 from:

%width=400px%Attach:control_simulation.gif

(:toggle hide gekko button show="Show GEKKO (Python) Code":)
(:div id=gekko:)

to:

Model Predictive Control is solved repeatedly at a specified cycle time as new measurements arrive or a new setpoint is specified. The model is updated with initial conditions from the next time step and the horizon is extended by an additional time unit to maintain a constant time horizon for predictive control.

(:toggle hide gekko2 button show="Show GEKKO (Python) Code":)
(:div id=gekko2:)

Added lines 196-197:

%width=400px%Attach:control_simulation.gif

October 18, 2024, at 05:17 PM by 10.35.117.248 -

Deleted line 56:

~~#%%Import packages~~

Deleted line 57:

~~from random import random~~

Deleted lines 60-62:

#%% Build model

#initialize GEKKO model

Deleted lines 61-62:

~~#time~~

Added line 70:

Deleted line 76:

~~#Equations~~

Deleted lines 78-80:

#%% Tuning

#global

Added lines 116-195:

%width=400px%Attach:control_simulation.gif

(:toggle hide gekko button show="Show GEKKO (Python) Code":)
(:div id=gekko:)
(:source lang=python:)
import numpy as np
from gekko import GEKKO
import matplotlib.pyplot as plt
from PIL import Image

# Initialize the GEKKO model
m = GEKKO(remote=False)
m.time = np.linspace(0, 20, 21) # Time grid
mass = 500
b = m.Param(value=50)
K = m.Param(value=0.8)
p = m.MV(value=0, lb=0, ub=100) # Manipulated Variable (gas pedal)
v = m.CV(value=0) # Controlled Variable (velocity)

# Define the dynamic equation
m.Equation(mass * v.dt() == -v * b + K * b * p)

# Set optimization mode for control
m.options.IMODE = 6 # Control mode

# Tuning parameters for MV (gas pedal)
p.STATUS = 1 # Allow optimizer to change
p.FSTATUS = 0
p.DCOST = 0.1 # Smooth out gas pedal movement
p.DMAX = 20 # Limit rate of change

# Tuning parameters for CV (velocity)
v.STATUS = 1 # Include set point in the objective
m.options.CV_TYPE = 2 # L2 norm
v.SP = 40 # Set point for velocity
v.TR_INIT = 1 # Enable set point trajectory
v.TAU = 5 # Time constant of trajectory

# List to store frames for GIF
frames = []

# Run the optimization in a loop and save plots
for i in range(25):
if i==15:
v.SP = 20

print(f'MPC cycle: {i} of 25')
m.solve(disp=False) # Solve the optimization problem

# Create a plot for each iteration
fig, (ax1, ax2) = plt.subplots(2, 1, figsize=(6, 8))

ax1.plot(m.time+i, p.value, 'b-', lw=2)
ax1.set_ylabel('Gas Pedal (%)')
ax1.set_title(f'MPC Cycle {i+1}')
ax1.set_ylim([0,100])

ax2.plot(m.time+i, v.value, 'r--', lw=2)
ax2.plot(m.time+i,np.ones_like(m.time)*v.SP,'k:',lw=2)
ax2.set_ylim([0,50])
ax2.set_ylabel('Velocity (m/s)')
ax2.set_xlabel('Time (s)')
fig.tight_layout()

# Save the plot as an image frame
fig.canvas.draw()
frame = np.frombuffer(fig.canvas.tostring_rgb(), dtype='uint8')
frame = frame.reshape(fig.canvas.get_width_height()[::-1] + (3,))
frames.append(Image.fromarray(frame))

plt.close(fig) # Close the figure to save memory

# Save the frames as a GIF
frames[0].save('control_simulation.gif', format='GIF',
append_images=frames[1:], save_all=True, duration=500, loop=0)

print("GIF saved as 'control_simulation.gif'")
(:sourceend:)
(:divend:)

December 06, 2022, at 10:52 PM by 136.36.4.38 -

Added lines 120-125:

* The first step to solve a model predictive control (MPC) in Python Gekko is to define the MPC problem. This may involve setting the time horizon and decision variables, constraints, objective function, and other parameters.
* Once the problem is defined, the next step is to create the MPC object. This is done by setting ''m.options.IMODE=6''. The parameters of the MPC object can be set in ''m.options''.
* The third step is to solve the MPC problem. This is done by calling the ''m.solve()'' command. This will solve the MPC problem and return the optimal solution.
* The fourth step is to process the solution. This can involve extracting the optimal decisions, evaluating the objective function, and other processing steps.
* Finally, the fifth step is to apply the optimal decisions. This can involve setting the plant state or inputs based on the optimal decisions.

November 17, 2021, at 12:48 AM by 10.35.117.248 -

Changed line 111 from:

plt.plot(m.time,p.value,'b-',~~LineWidth~~=2)

to:

plt.plot(m.time,p.value,'b-',lw=2)

Changed line 114 from:

plt.plot(m.time,v.value,'r--',~~LineWidth~~=2)

to:

plt.plot(m.time,v.value,'r--',lw=2)

January 31, 2018, at 01:52 AM by 10.37.66.71 -

Added lines 53-119:

(:toggle hide gekko button show="Show GEKKO (Python) Code":)
(:div id=gekko:)
(:source lang=python:)
#%%Import packages
import numpy as np
from random import random
from gekko import GEKKO
import matplotlib.pyplot as plt

#%% Build model

#initialize GEKKO model
m = GEKKO()

#time
m.time = np.linspace(0,20,41)

#constants
mass = 500

#Parameters
b = m.Param(value=50)
K = m.Param(value=0.8)
#Manipulated variable
p = m.MV(value=0, lb=0, ub=100)

#Controlled Variable
v = m.CV(value=0)

#Equations
m.Equation(mass*v.dt() == -v*b + K*b*p)

#%% Tuning

#global
m.options.IMODE = 6 #control

#MV tuning
p.STATUS = 1 #allow optimizer to change
p.DCOST = 0.1 #smooth out gas pedal movement
p.DMAX = 20 #slow down change of gas pedal

#CV tuning
#setpoint
v.STATUS = 1 #add the SP to the objective
m.options.CV_TYPE = 2 #L2 norm
v.SP = 40 #set point
v.TR_INIT = 1 #setpoint trajectory
v.TAU = 5 #time constant of setpoint trajectory

#%% Solve

m.solve()

#%% Plot solution
plt.figure()
plt.subplot(2,1,1)
plt.plot(m.time,p.value,'b-',LineWidth=2)
plt.ylabel('gas')
plt.subplot(2,1,2)
plt.plot(m.time,v.value,'r--',LineWidth=2)
plt.ylabel('velocity')
plt.xlabel('time')
plt.show()
(:sourceend:)
(:divend:)

October 31, 2016, at 03:18 PM by 10.5.113.135 -

Changed line 17 from:

'''Objective:''' Implement a model predictive controller that automatically regulates vehicle velocity. Implement the controller in Excel, MATLAB, Python, or Simulink and tune the controller for acceptable performance. Discuss factors that may be ~~import~~ for evaluating controller performance. ''Estimated time: 1 hour''

to:

'''Objective:''' Implement a model predictive controller that automatically regulates vehicle velocity. Implement the controller in Excel, MATLAB, Python, or Simulink and tune the controller for acceptable performance. Discuss factors that may be important for evaluating controller performance. ''Estimated time: 1 hour''

May 18, 2015, at 04:01 PM by 45.56.3.184 -

Changed line 17 from:

'''Objective:''' Implement a model predictive controller that automatically regulates vehicle velocity. Implement the controller in ~~MATLAB~~, ~~Python~~, ~~and~~ Simulink and tune the controller for acceptable performance. Discuss factors that may be import for evaluating controller performance. ''Estimated time: ~~2 hours~~''

to:

'''Objective:''' Implement a model predictive controller that automatically regulates vehicle velocity. Implement the controller in Excel, MATLAB, Python, or Simulink and tune the controller for acceptable performance. Discuss factors that may be import for evaluating controller performance. ''Estimated time: 1 hour''

May 18, 2015, at 03:57 PM by 45.56.3.184 -

Deleted lines 11-12:

~~!!!!Model Predictive Control Example~~

May 18, 2015, at 03:55 PM by 45.56.3.184 -

Changed lines 21-22 from:

# The dynamic relationship between a vehicle gas pedal position (MV) and velocity (CV) is given by the following set of conditions and a single dynamic equation.

to:

The dynamic relationship between a vehicle gas pedal position (MV) and velocity (CV) is given by the following set of conditions and a single dynamic equation.

Changed line 36 from:

# Discuss the controller performance and how it could be tuned to meet multiple objectives including:

to:

Discuss the controller performance and how it could be tuned to meet multiple objectives including:

May 18, 2015, at 03:55 PM by 45.56.3.184 -

Changed lines 15-16 from:

A method to solve dynamic control problems is by numerically integrating the dynamic model at discrete time intervals, much like measuring a physical system at particular time points. The numerical solution is compared to a desired trajectory and the difference is minimized by adjustable parameters in the model that may change at every time step. The first control action is taken and then the entire process is repeated at the next time instance. The process is repeated because objective targets may change or updated measurements may have adjusted parameter or state estimates. Excel, MATLAB, Python, and Simulink are used in the following example to both solve the differential equations that describe the velocity of a vehicle as well as minimize the control objective function.

to:

A method to solve dynamic control problems is by numerically integrating the dynamic model at discrete time intervals, much like measuring a physical system at particular time points. The numerical solution is compared to a desired trajectory and the difference is minimized by adjustable parameters in the model that may change at every time step. The first control action is taken and then the entire process is repeated at the next time instance. The process is repeated because objective targets may change or updated measurements may have adjusted parameter or state estimates.

!!!! Exercise

'''Objective:''' Implement a model predictive controller that automatically regulates vehicle velocity. Implement the controller in MATLAB, Python, and Simulink and tune the controller for acceptable performance. Discuss factors that may be import for evaluating controller performance. ''Estimated time: 2 hours''

# The dynamic relationship between a vehicle gas pedal position (MV) and velocity (CV) is given by the following set of conditions and a single dynamic equation.

Constants
m = 500 ! Mass (kg)
Parameters
b = 50 ! Resistive coefficient (N-s/m)
K = 0.8 ! Gain (m/s-%pedal)
p = 0 >= 0 <= 100 ! Gas pedal position (%)
Variables
v = 0 ! initial condition
Equations
m * $v = -v * b + K * b * p

Implement a model predictive controller that adjusts gas pedal position to regulate velocity. Start at an initial vehicle velocity of 0 m/s and accelerate to a velocity of 40 m/s.

# Discuss the controller performance and how it could be tuned to meet multiple objectives including:

* minimize travel time
* remain within speed limits
* improve vehicle fuel efficiency
* discourage excessive gas pedal adjustments
* do not accelerate excessively

There is no need to implement these advanced objectives in simulation for this second part of the exercise, only discuss the possible competing objectives.

!!!! Solution

Added lines 50-51:

Excel, MATLAB, Python, and Simulink are used in this example to both solve the differential equations that describe the velocity of a vehicle as well as minimize the control objective function.

Deleted line 54:

May 18, 2015, at 03:28 PM by 45.56.3.184 -

Deleted lines 22-45:

!!!! Exercise

'''Objective:''' Set up and solve several [[Attach:Dynamic_Optimization_Benchmarks.pdf|dynamic optimization benchmark problems]]'^1^'. Create a program to optimize and display the results. ''Estimated Time (each): 30 minutes''

* Example 1a - Nonlinear, unconstrained, minimize final state
->Attach:dynopt_1a.png
* Example 1b - Nonlinear, unconstrained, minimize final state with terminal constraint
->Attach:dynopt_1b.png
* Example 2 - Nonlinear, constrained, minimize final state
->Attach:dynopt_2.png
* Example 3 - Tubular reactor with parallel reaction
->Attach:dynopt_3.png
* Example 4 - Batch reactor with consecutive reactions A->B->C
->Attach:dynopt_4.png
Example 5 - Catalytic plug flow reactor with A->B->C
->Attach:dynopt_5.png

!!!! Solution

Attach:download.png [[Attach:dynamic_optimization_benchmarks.zip|Dynamic Optimization Benchmarks in MATLAB and Python]]

!!!! References

# M. Čižniar, M. Fikar, M.A. Latifi: A MATLAB Package for Dynamic Optimisation of Processes, 7th International Scientific – Technical Conference – Process Control 2006, June 13 – 16, 2006, Kouty nad Desnou, Czech Republic. [[Attach:DynOpt_Benchmarks.pdf|Article]]

April 28, 2015, at 05:52 PM by 45.56.3.184 -

Changed lines 25-26 from:

'''Objective:''' Set up and solve several dynamic optimization benchmark problems'^1^'. Create a program to optimize and display the results. ''Estimated Time (each): 30 minutes''

to:

'''Objective:''' Set up and solve several [[Attach:Dynamic_Optimization_Benchmarks.pdf|dynamic optimization benchmark problems]]'^1^'. Create a program to optimize and display the results. ''Estimated Time (each): 30 minutes''

Changed line 42 from:

Attach:download.png [[Attach:~~dynopt~~_benchmarks.zip|Dynamic Optimization Benchmarks in MATLAB and Python]]

to:

Attach:download.png [[Attach:dynamic_optimization_benchmarks.zip|Dynamic Optimization Benchmarks in MATLAB and Python]]

April 28, 2015, at 05:50 PM by 45.56.3.184 -

Changed line 46 from:

# M. Čižniar, M. Fikar, M.A. Latifi: A MATLAB Package for Dynamic Optimisation of Processes, 7th International Scientific – Technical Conference – Process Control 2006, June 13 – 16, 2006, Kouty nad Desnou, Czech Republic.

to:

# M. Čižniar, M. Fikar, M.A. Latifi: A MATLAB Package for Dynamic Optimisation of Processes, 7th International Scientific – Technical Conference – Process Control 2006, June 13 – 16, 2006, Kouty nad Desnou, Czech Republic. [[Attach:DynOpt_Benchmarks.pdf|Article]]

April 27, 2015, at 11:11 PM by 10.5.113.179 -

Added lines 23-46:

!!!! Exercise

'''Objective:''' Set up and solve several dynamic optimization benchmark problems'^1^'. Create a program to optimize and display the results. ''Estimated Time (each): 30 minutes''

* Example 1a - Nonlinear, unconstrained, minimize final state
->Attach:dynopt_1a.png
* Example 1b - Nonlinear, unconstrained, minimize final state with terminal constraint
->Attach:dynopt_1b.png
* Example 2 - Nonlinear, constrained, minimize final state
->Attach:dynopt_2.png
* Example 3 - Tubular reactor with parallel reaction
->Attach:dynopt_3.png
* Example 4 - Batch reactor with consecutive reactions A->B->C
->Attach:dynopt_4.png
Example 5 - Catalytic plug flow reactor with A->B->C
->Attach:dynopt_5.png

!!!! Solution

Attach:download.png [[Attach:dynopt_benchmarks.zip|Dynamic Optimization Benchmarks in MATLAB and Python]]

!!!! References

# M. Čižniar, M. Fikar, M.A. Latifi: A MATLAB Package for Dynamic Optimisation of Processes, 7th International Scientific – Technical Conference – Process Control 2006, June 13 – 16, 2006, Kouty nad Desnou, Czech Republic.

April 04, 2015, at 02:03 PM by 45.56.12.124 -

Changed line 17 from:

* [[Attach:model_predictive_control.zip|MPC in Excel, MATLAB, Python, and Simulink]]

to:

Attach:download.png [[Attach:model_predictive_control.zip|MPC in Excel, MATLAB, Python, and Simulink]]

April 03, 2015, at 01:40 AM by 10.24.17.95 -

Added lines 6-7:

* [[Attach:Intro_Dynamic_Control.pdf|Introduction to Dynamic Control/Optimization (pdf)]]

April 02, 2015, at 11:12 PM by 10.10.145.116 -

Changed line 18 from:

to:

April 02, 2015, at 04:48 PM by 45.56.12.124 -

Added lines 6-9:

(:html:)
<iframe width="560" height="315" src="https://www.youtube.com/embed/DFqOf5wbQtc?rel=0" frameborder="0" allowfullscreen></iframe>
(:htmlend:)

April 02, 2015, at 04:19 PM by 45.56.12.124 -

Changed lines 5-16 from:

Dynamic control is a method to use model predictions to plan an optimized future trajectory for time-varying systems. It is often referred to as Model Predictive Control (MPC) or Dynamic Optimization.

to:

Dynamic control is a method to use model predictions to plan an optimized future trajectory for time-varying systems. It is often referred to as Model Predictive Control (MPC) or Dynamic Optimization.

!!!!Model Predictive Control Example

A method to solve dynamic control problems is by numerically integrating the dynamic model at discrete time intervals, much like measuring a physical system at particular time points. The numerical solution is compared to a desired trajectory and the difference is minimized by adjustable parameters in the model that may change at every time step. The first control action is taken and then the entire process is repeated at the next time instance. The process is repeated because objective targets may change or updated measurements may have adjusted parameter or state estimates. Excel, MATLAB, Python, and Simulink are used in the following example to both solve the differential equations that describe the velocity of a vehicle as well as minimize the control objective function.

* [[Attach:model_predictive_control.zip|MPC in Excel, MATLAB, Python, and Simulink]]

(:html:)
<iframe width="560" height="315" src="https://www.youtube.com/embed/umxAfu44kWo?rel=0" frameborder="0" allowfullscreen></iframe>
(:htmlend:)

March 26, 2015, at 04:08 PM by 127.0.0.1 -

Added lines 1-5:

(:title Dynamic Control Introduction:)
(:keywords control, dynamics, dynamic optimization, simulation, modeling language, differential, algebraic, tutorial:)
(:description Dynamic control in MATLAB and Python for use in real-time or off-line applications:)

Dynamic control is a method to use model predictions to plan an optimized future trajectory for time-varying systems. It is often referred to as Model Predictive Control (MPC) or Dynamic Optimization.

Dynamic Optimization

Dynamic Control Introduction

Main.DynamicControl History

Search

Options:

Dynamic Optimization