] u = ["u1"] u_bounds = {"u1": (162.183495794 * 0.5, 162.183495794 * 1.6)} ref_state = {("c_capture", ((), )): 0.52} c = NmpcGen(d_mod=bfb_dae, u=u, states=states, ref_state=ref_state, u_bounds=u_bounds) c.ss.dref = snap c.load_iguess_ss() c.solve_ss() c.load_d_s(c.d1) c.solve_d(c.d1) c.update_state_real() # update the current state c.find_target_ss() c.create_nmpc() c.update_targets_nmpc() c.compute_QR_nmpc(n=-1) c.new_weights_olnmpc(1000, 1e+06) # q_cov = {} # for i in range(1, nfet): # for j in range(1, ntrays + 1): # q_cov[("x", (j,)), ("x", (j,)), i] = 1e-05 # q_cov[("M", (j,)), ("M", (j,)), i] = 0.5 # # m_cov = {}
states = ["Ngb", "Hgb", "Ngc", "Hgc", "Nsc", "Hsc", "Nge", "Hge", "Nse", "Hse", "mom"] u = ["u1"] u_bounds = {"u1":(162.183495794 * 0.0005, 162.183495794 * 10000)} # ref_state = {("c_capture", ((),)): 0.50} ref_state = {("c_capture", ((),)): 0.40} c = NmpcGen(d_mod=bfb_dae, u=u, states=states, ref_state=ref_state, u_bounds=u_bounds) c.ss.dref = snap c.load_iguess_ss() # sys.exit() # c.ss.lydot.display() # sys.exit() c.solve_ss() c.solve_d(c.ss) c.ss.write_nl() # sys.exit() c.load_d_s(c.d1) c.d1.create_bounds() c.solve_d(c.d1, halt_on_ampl_error=True) c.update_state_real() # update the current state c.find_target_ss() with open("mult_boundsss2.txt", "w") as f: c.ss2.ipopt_zL_out.pprint(ostream=f) f.close() c.create_nmpc() c.update_targets_nmpc() c.compute_QR_nmpc(n=-1) c.new_weights_olnmpc(10, 1e+06)
u_bounds = { "u1": (0.0001, 9286.03346463 * 100), "u2": (0.0001, 9286.03346463 * 100) } ref_state = {("c_capture", ((), )): 0.50} # weights = {("c_capture", ((),)): 1E+05} c = NmpcGen(d_mod=bfb_dae, u=u, states=states, ref_state=ref_state, u_bounds=u_bounds) c.load_iguess_ss() c.ss.create_bounds() c.solve_ss() c.load_d_s(c.d1) c.solve_d(c.d1, iter_max=100) c.find_target_ss() c.ss2.write_nl() sys.exit(-1) # c.plant_input_gen(c.ss2, 1) c.create_nmpc() c.update_targets_nmpc() # sys.exit() c.compute_QR_nmpc(n=-1) c.initialize_olnmpc(c.d1) c.olnmpc.display(filename="somefile0.txt") c.create_suffixes() # c.solve_k_aug_nmpc() # c.olnmpc.write_nl()
"x": [(i, ) for i in range(1, ntrays + 1)], "M": [(i, ) for i in range(1, ntrays + 1)] } nfet = 10 tfe = [i for i in range(1, nfet + 1)] c = NmpcGen(d_mod=DistDiehlNegrete, u=u, states=states, ref_state=ref_state, u_bounds=u_bounds) c.load_iguess_ss() c.solve_ss() c.load_d_s(c.d1) c.solve_d(c.d1) c.update_state_real() # update the current state c.find_target_ss() c.create_nmpc() c.update_targets_nmpc() c.compute_QR_nmpc(n=-1) c.new_weights_olnmpc(10000, 1e+06) c.d1.create_bounds() c.create_predictor() c.predictor_step(c.d1, "real") q_cov = {} for j in range(1, ntrays + 1):
from numpy.random import normal as npm from pyutilib.common._exceptions import ApplicationError import sys states = ["Ngb", "Hgb", "Ngc", "Hgc", "Nsc", "Hsc", "Nge", "Hge", "Nse", "Hse", "mom"] u = ["u1"] u_bounds = {"u1":(162.183495794 * 0.5, 162.183495794 * 1.6)} ref_state = {("c_capture", ((),)): 0.52} c = NmpcGen(d_mod=bfb_dae, u=u, states=states, ref_state=ref_state, u_bounds=u_bounds) c.ss.dref = snap c.load_iguess_ss() c.solve_ss() c.load_d_s(c.d1) c.solve_d(c.d1) c.update_state_real() # update the current state c.find_target_ss() c.create_nmpc() c.update_targets_nmpc() c.compute_QR_nmpc(n=-1) c.new_weights_olnmpc(1000, 1e+06) # q_cov = {} # for i in range(1, nfet): # for j in range(1, ntrays + 1): # q_cov[("x", (j,)), ("x", (j,)), i] = 1e-05 # q_cov[("M", (j,)), ("M", (j,)), i] = 0.5 #