Move PID matplotlib demo harnesses out of production modules
temp_controller.py, temp_controller_smith.py, and kalman.py imported matplotlib at module level just to support eyeballed-plot __main__ blocks, coupling the live server's import graph to a GUI plotting lib it never uses at runtime. Relocate those demos (and kalman_eval.py) to scripts/demos/pid/ and strip the now-unused imports/__main__ blocks from the production files. Co-Authored-By: Claude Sonnet 4.6 <noreply@anthropic.com>
This commit is contained in:
Claude Sonnet 4.6
@@ -0,0 +1,56 @@
|
||||
import numpy as np
|
||||
from matplotlib.pyplot import plot, figure, subplot, grid, show
|
||||
from components.pid.kalman import Kalman
|
||||
|
||||
|
||||
if __name__ == '__main__':
|
||||
dt = 1.0
|
||||
|
||||
params = {
|
||||
'var_P' : 1,
|
||||
'var_Q' : 0.0001,
|
||||
'var_R' : 0.5
|
||||
}
|
||||
k = Kalman(dt, params)
|
||||
|
||||
_x1 = np.empty(0)
|
||||
_y1 = np.empty(0)
|
||||
_x2 = np.empty(0)
|
||||
_y2 = np.empty(0)
|
||||
|
||||
N = int(1000/dt)
|
||||
seqn = range(0, N)
|
||||
|
||||
_seqn = range(0, 2*N)
|
||||
X = np.array([dt, 0])
|
||||
for n in seqn:
|
||||
Z = k.process_measurement(X, 0.5)
|
||||
# Kalman.print("Z:", Z)
|
||||
Xp = k.process(Z)
|
||||
|
||||
_x1 = np.append(_x1, Z[0])
|
||||
_x2 = np.append(_x2, Z[1])
|
||||
_y1 = np.append(_y1, Xp[0])
|
||||
_y2 = np.append(_y2, Xp[1])
|
||||
|
||||
for n in seqn:
|
||||
X[0] = X[0] + dt
|
||||
Z = k.process_measurement(X, 0.5)
|
||||
# Kalman.print("Z:", Z)
|
||||
Xp = k.process(Z)
|
||||
|
||||
_x1 = np.append(_x1, Z[0])
|
||||
_x2 = np.append(_x2, Z[1])
|
||||
_y1 = np.append(_y1, Xp[0])
|
||||
_y2 = np.append(_y2, Xp[1])
|
||||
|
||||
figure(1)
|
||||
subplot(2, 1, 1)
|
||||
plot(_seqn, _x1, 'bx', _seqn, _y1, '-r', linewidth=1)
|
||||
grid(True)
|
||||
subplot(2, 1, 2)
|
||||
plot(_seqn, _x2, 'bx', _seqn, _y2, '-r', linewidth=1)
|
||||
grid(True)
|
||||
show()
|
||||
|
||||
print("End of program")
|
||||
@@ -0,0 +1,48 @@
|
||||
import numpy as np
|
||||
from matplotlib.pyplot import plot, figure, subplot, legend, grid, show
|
||||
|
||||
|
||||
T_true = 72
|
||||
T_est = 68
|
||||
E_est = 4
|
||||
E_mea = 2
|
||||
|
||||
_t = np.empty(0)
|
||||
_T_true = np.empty(0)
|
||||
_KG = np.empty(0)
|
||||
_T_est = np.empty(0)
|
||||
_E_est = np.empty(0)
|
||||
|
||||
for t in range(0, 10000):
|
||||
_t = np.append(_t, t)
|
||||
_T_true = np.append(_T_true, T_true)
|
||||
|
||||
T_meas = T_true + E_mea*np.random.randn()
|
||||
|
||||
# 1: Calculate Kalman gain KG
|
||||
KG = E_est / (E_est + E_mea)
|
||||
|
||||
# 2: Calculate estimation T_est
|
||||
T_est = T_est + KG*(T_meas - T_est)
|
||||
|
||||
# 3: Calculate error in estimate E_est
|
||||
E_est = (1-KG)*E_est
|
||||
|
||||
_KG = np.append(_KG, KG)
|
||||
_T_est = np.append(_T_est, T_est)
|
||||
_E_est = np.append(_E_est, E_est)
|
||||
|
||||
figure(1)
|
||||
subplot(3, 1, 1)
|
||||
plot(_t, _T_est, 'b-', _t, _T_true, '-r', linewidth=1)
|
||||
legend(["T_est", "T_true"])
|
||||
grid(True)
|
||||
subplot(3, 1, 2)
|
||||
plot(_t, _E_est, 'b-', linewidth=1)
|
||||
legend(["E_est"])
|
||||
grid(True)
|
||||
subplot(3, 1, 3)
|
||||
plot(_t, _KG, '-r', linewidth=1)
|
||||
legend(["KG"])
|
||||
grid(True)
|
||||
show()
|
||||
@@ -0,0 +1,77 @@
|
||||
import numpy as np
|
||||
from matplotlib.pyplot import plot, figure, subplot, grid, show, legend
|
||||
from components.plant.pot import Pot
|
||||
from components.pid.temp_controller import TempController
|
||||
from components.pid.tc_constants import Test
|
||||
|
||||
|
||||
if __name__ == '__main__':
|
||||
dt = 1.0
|
||||
|
||||
temp_ist = 0
|
||||
temp_soll = 20
|
||||
ctrl = TempController(dt, Test.tc_ctrl_params)
|
||||
plant = Pot(dt, Test.tc_pot_params)
|
||||
_temp_ist = np.empty(0)
|
||||
_temp_soll = np.empty(0)
|
||||
_y = np.empty(0)
|
||||
_fb = np.empty(0)
|
||||
_t = np.empty(0)
|
||||
_heatrate_ist_kalman = np.empty(0)
|
||||
_temp_ist_kalman = np.empty(0)
|
||||
a = 0.5
|
||||
fb = 0
|
||||
rho = 0.02
|
||||
temps = [{'Temp': 20, 'Duration': 1000}, {'Temp': 40, 'Duration': 1000}, {'Temp': 50, 'Duration': 1000}, {'Temp': 60, 'Duration': 1000}, {'Temp': 70, 'Duration': 1000}, {'Temp': 80, 'Duration': 1000}, {'Temp': 78, 'Duration': 1000}]
|
||||
|
||||
t = 0
|
||||
for temp in temps:
|
||||
temp_soll = temp['Temp']
|
||||
hold_counter = temp['Duration']
|
||||
hold = False
|
||||
ctrl.set_theta_soll(temp_soll)
|
||||
ctrl.set_heatrate_soll(1.0)
|
||||
|
||||
while True:
|
||||
if hold:
|
||||
if hold_counter == 0:
|
||||
break
|
||||
hold_counter -= 1
|
||||
plant.process()
|
||||
temp_ist = plant.get_temperature() + 0.0 * np.random.randn()
|
||||
ctrl.set_theta_ist(temp_ist)
|
||||
ctrl.process()
|
||||
|
||||
y = 3500*ctrl.get_power()
|
||||
power = max(0, y)
|
||||
plant.set_power(power)
|
||||
fb = plant.get_power()
|
||||
if abs(temp_ist - temp_soll) < 0.1:
|
||||
hold = True
|
||||
|
||||
temp_ist -= rho
|
||||
_temp_ist = np.append(_temp_ist, temp_ist)
|
||||
_temp_soll = np.append(_temp_soll, temp_soll)
|
||||
_y = np.append(_y, y)
|
||||
_fb = np.append(_fb, fb)
|
||||
_t = np.append(_t, t)
|
||||
_heatrate_ist_kalman = np.append(_heatrate_ist_kalman, max(-1, min(3, ctrl.heatrate_ist)))
|
||||
_temp_ist_kalman = np.append(_temp_ist, ctrl.theta_ist)
|
||||
t += 1
|
||||
|
||||
figure(1)
|
||||
subplot(3, 1, 1)
|
||||
plot(_t, _temp_ist, _t, _temp_soll, 'r-', linewidth=1)
|
||||
legend(["ist", "soll"])
|
||||
grid(True)
|
||||
subplot(3, 1, 2)
|
||||
plot(_t, _y, '-b', _t, _fb, '-r', linewidth=1)
|
||||
legend(["y", "pot"])
|
||||
grid(True)
|
||||
subplot(3, 1, 3)
|
||||
plot(_t, _heatrate_ist_kalman, '-b', linewidth=1)
|
||||
legend(["heatrate"])
|
||||
grid(True)
|
||||
show()
|
||||
|
||||
print("End of program")
|
||||
@@ -0,0 +1,90 @@
|
||||
import numpy as np
|
||||
from matplotlib.pyplot import plot, figure, subplot, grid, show, legend
|
||||
from components.plant.pot import Pot
|
||||
from components.pid.temp_controller_smith import TempController
|
||||
from components.pid.tc_constants import Test
|
||||
|
||||
|
||||
if __name__ == '__main__':
|
||||
dt = 1.0
|
||||
|
||||
temp_soll = 20
|
||||
ctrl = TempController(dt, Test.tc_ctrl_params, Test.tc_model_params)
|
||||
plant = Pot(dt, Test.tc_pot_params)
|
||||
_y = np.empty(0)
|
||||
_fb = np.empty(0)
|
||||
_t = np.empty(0)
|
||||
_temp_soll = np.empty(0)
|
||||
_temp_ist_kalman = np.empty(0)
|
||||
_heatrate_ist_kalman = np.empty(0)
|
||||
_temp_ist_kalman_plant = np.empty(0)
|
||||
_heatrate_ist_kalman_plant = np.empty(0)
|
||||
_temp_ist_kalman_model = np.empty(0)
|
||||
_heatrate_ist_kalman_model = np.empty(0)
|
||||
a = 0.5
|
||||
fb = 0
|
||||
rho = 0.02
|
||||
temps = [{'Temp': 20, 'Duration': 1000}, {'Temp': 40, 'Duration': 1000}, {'Temp': 50, 'Duration': 1000}, {'Temp': 60, 'Duration': 1000}, {'Temp': 70, 'Duration': 1000}, {'Temp': 80, 'Duration': 1000}, {'Temp': 78, 'Duration': 1000}]
|
||||
|
||||
t = 0
|
||||
|
||||
for temp in temps:
|
||||
temp_soll = temp['Temp']
|
||||
hold_counter = temp['Duration']
|
||||
hold = False
|
||||
ctrl.set_theta_soll(temp_soll)
|
||||
ctrl.set_heatrate_soll(1.0)
|
||||
|
||||
while True:
|
||||
if hold:
|
||||
if hold_counter == 0:
|
||||
break
|
||||
hold_counter -= 1
|
||||
plant.process()
|
||||
temp_ist = plant.get_temperature()
|
||||
ctrl.set_theta_ist(temp_ist)
|
||||
ctrl.process()
|
||||
|
||||
y = 3500*ctrl.get_power()
|
||||
power = max(0, y)
|
||||
plant.set_power(power)
|
||||
fb = plant.get_power()
|
||||
if abs(temp_ist - temp_soll) < 0.1:
|
||||
hold = True
|
||||
|
||||
_y = np.append(_y, y)
|
||||
_fb = np.append(_fb, fb)
|
||||
_t = np.append(_t, t)
|
||||
_temp_soll = np.append(_temp_soll, temp_soll)
|
||||
_temp_ist_kalman_plant = np.append(_temp_ist_kalman_plant, ctrl.theta_ist_plant)
|
||||
_heatrate_ist_kalman_plant = np.append(_heatrate_ist_kalman_plant, max(-1, min(3, ctrl.dtheta_ist_plant)))
|
||||
_temp_ist_kalman_model = np.append(_temp_ist_kalman_model, ctrl.theta_ist_model_delay)
|
||||
_heatrate_ist_kalman_model = np.append(_heatrate_ist_kalman_model, max(-1, min(3, ctrl.dtheta_ist_model_delay)))
|
||||
t += 1
|
||||
|
||||
figure(1)
|
||||
subplot(3, 1, 1)
|
||||
plot(_t, _temp_ist_kalman_plant, _t, _temp_soll, 'r-', linewidth=1)
|
||||
legend(["ist", "soll"])
|
||||
grid(True)
|
||||
subplot(3, 1, 2)
|
||||
plot(_t, _y, '-b', _t, _fb, '-r', linewidth=1)
|
||||
legend(["y", "pot"])
|
||||
grid(True)
|
||||
subplot(3, 1, 3)
|
||||
plot(_t, _heatrate_ist_kalman_plant, '-b', linewidth=1)
|
||||
legend(["heatrate"])
|
||||
grid(True)
|
||||
|
||||
figure(2)
|
||||
subplot(2, 1, 1)
|
||||
plot(_t, _temp_ist_kalman_plant, '-b', _t, _temp_ist_kalman_model, 'r-', linewidth=1)
|
||||
legend(["plant", "model"])
|
||||
grid(True)
|
||||
subplot(2, 1, 2)
|
||||
plot(_t, _heatrate_ist_kalman_plant, '-b', _t, _heatrate_ist_kalman_model, '-r', linewidth=1)
|
||||
legend(["plant", "model"])
|
||||
grid(True)
|
||||
show()
|
||||
|
||||
print("End of program")
|
||||
Reference in New Issue
Block a user