experiment_model_2 = Model(thermostat = therm.HarshCurtailedThermostat(strategy = greedy_strategy), time_step_in_hours = timestep)
experiment_model_3 = Model(thermostat = therm.UpsideCurtailedThermostat(strategy = greedy_strategy), time_step_in_hours = timestep)
total_cost_control = 0
total_cost_experiment = 0
total_cost_experiment_2 = 0
total_cost_experiment_3 = 0
max_temp_control = 0
max_temp_exp = 0
max_temp_exp_2 = 0
max_temp_exp_3 = 0
for i in range(1, tempdata.timespan):
ext_temp = tempdata.external_temp(i)
no_cooling_model.iterate(ext_temp)
control_model.iterate(ext_temp)
d = tempdata.steptime_to_date(i)
greedy_strategy.iterate(d)
experiment_model.iterate(ext_temp)
experiment_model_2.iterate(ext_temp)
experiment_model_3.iterate(ext_temp)
price_mwh = pricedata.price(d)
control_cost = control_model.power_cost(price_mwh)
experiment_cost = experiment_model.power_cost(price_mwh)
experiment_cost_2 = experiment_model_2.power_cost(price_mwh)
experiment_cost_3 = experiment_model_3.power_cost(price_mwh)
if greedy_strategy.current_swing > 0:
total_cost_control += control_cost
total_cost_experiment += experiment_cost
total_cost_experiment_2 += experiment_cost_2
timestep = 1/60.0
tempdata = Tempdata("hourly-temp.csv", timestep)
no_cooling_model = Model(thermostat = NoThermostat(), time_step_in_hours = timestep)
control_model = Model(thermostat = ControlThermostat(deadband_delta_celcius = 0.5), time_step_in_hours = timestep)
control_model_2 = Model(thermostat = ControlThermostat(deadband_delta_celcius = 1.0), time_step_in_hours = timestep)
control_model.temp = util.f_to_c(78)
control_model.therm.state = 1
i = 0
while control_model.therm.state:
#print i, util.c_to_f(control_model.temp)
i+=1
control_model.iterate(util.f_to_c(85))
no_cooling_model = Model(thermostat = NoThermostat(), time_step_in_hours = timestep, internal_temp_celcius = util.f_to_c(68))
no_cooling_model_2 = Model(thermostat = NoThermostat(), time_step_in_hours = timestep, internal_temp_celcius = util.f_to_c(68))
no_cooling_model_3 = Model(thermostat = NoThermostat(), time_step_in_hours = timestep, internal_temp_celcius = util.f_to_c(68))
i = 0
while no_cooling_model.temp < util.f_to_c(85) and i < 420:
print i, util.c_to_f(no_cooling_model.temp)
i+=1
no_cooling_model.iterate(util.f_to_c(85))
#for i in range(0, 200):
# print "%d, %s, %s, %s" % (i, util.c_to_f(no_cooling_model.temp), util.c_to_f(no_cooling_model_2.temp), util.c_to_f(no_cooling_model_3.temp))
# no_cooling_model.iterate(util.f_to_c(85))