add funs

PNNLCodes/CoSim/func_EDC_CoSim_test.m

@@ -9,9 +9,8 @@
 num_zone = 16; % number of zones
 
 %% load the coefficients a_0 ... a_5
-top_floor = load(strcat('paras/Top_floor/', sample_time, '/top_floor.mat'));
-mid_floor = load(strcat('paras/Mid_floor/', sample_time, '/mid_floor.mat'));
-bot_floor = load(strcat('paras/Bot_floor/', sample_time, '/bot_floor.mat'));
+load(strcat('paras/zone/', sample_time, '/zones.mat'));
+
 
 a_0 = zeros(num_zone, 1);
 a_1 = zeros(num_zone, 1);
@@ -21,38 +20,14 @@
 a_5 = zeros(num_zone, 1);
 
 for i = 1 : num_zone
-    if  1 <= i && i <= 5 % top_floor
-        a_0(i) = eval(strcat('top_floor.zone4_', num2str(i), '_summer_result.a0'));
-        a_1(i) = eval(strcat('top_floor.zone4_', num2str(i), '_summer_result.a1'));
-        a_2(i) = eval(strcat('top_floor.zone4_', num2str(i), '_summer_result.a2'));
-        a_3(i) = eval(strcat('top_floor.zone4_', num2str(i), '_summer_result.a3'));
-        a_4(i) = eval(strcat('top_floor.zone4_', num2str(i), '_summer_result.a4'));
-        a_5(i) = eval(strcat('top_floor.zone4_', num2str(i), '_summer_result.a5'));
-    end
-    if  6 <= i && i <= 10 % mid_floor
-        a_0(i) = eval(strcat('mid_floor.zone4_', num2str(i-5), '_summer_result.a0'));
-        a_1(i) = eval(strcat('mid_floor.zone4_', num2str(i-5), '_summer_result.a1'));
-        a_2(i) = eval(strcat('mid_floor.zone4_', num2str(i-5), '_summer_result.a2'));
-        a_3(i) = eval(strcat('mid_floor.zone4_', num2str(i-5), '_summer_result.a3'));
-        a_4(i) = eval(strcat('mid_floor.zone4_', num2str(i-5), '_summer_result.a4'));
-        a_5(i) = eval(strcat('mid_floor.zone4_', num2str(i-5), '_summer_result.a5'));
-    end
-    if  11 <= i && i <= 15 % bot_floor
-        a_0(i) = eval(strcat('bot_floor.zone4_', num2str(i-10), '_summer_result.a0'));
-        a_1(i) = eval(strcat('bot_floor.zone4_', num2str(i-10), '_summer_result.a1'));
-        a_2(i) = eval(strcat('bot_floor.zone4_', num2str(i-10), '_summer_result.a2'));
-        a_3(i) = eval(strcat('bot_floor.zone4_', num2str(i-10), '_summer_result.a3'));
-        a_4(i) = eval(strcat('bot_floor.zone4_', num2str(i-10), '_summer_result.a4'));
-        a_5(i) = eval(strcat('bot_floor.zone4_', num2str(i-10), '_summer_result.a5'));
-    end
-    if  i == 16 % basement
-        a_0(i) = eval(strcat('bot_floor.zone', num2str(i-10), '_summer_result.a0'));
-        a_1(i) = eval(strcat('bot_floor.zone', num2str(i-10), '_summer_result.a1'));
-        a_2(i) = eval(strcat('bot_floor.zone', num2str(i-10), '_summer_result.a2'));
-        a_3(i) = eval(strcat('bot_floor.zone', num2str(i-10), '_summer_result.a3'));
-        a_4(i) = eval(strcat('bot_floor.zone', num2str(i-10), '_summer_result.a4'));
-        a_5(i) = eval(strcat('bot_floor.zone', num2str(i-10), '_summer_result.a5'));
-    end
+
+        a_0(i) = eval(strcat('z', num2str(i-1), '.a0'));
+        a_1(i) = eval(strcat('z', num2str(i-1), '.a1'));
+        a_2(i) = eval(strcat('z', num2str(i-1), '.a2'));
+        a_3(i) = eval(strcat('z', num2str(i-1), '.a3'));
+        a_4(i) = eval(strcat('z', num2str(i-1), '.a4'));
+        a_5(i) = eval(strcat('z', num2str(i-1), '.a5'));
+
 end
 
 
@@ -268,4 +243,4 @@ expression J(N_schd)
 % SolverStatus.Status = cvx_status;
 % SolverStatus.OptVal = cvx_optval;
 
-Optimal_Temp_Ctrl = T_z(:,1);
+Optimal_Temp_Ctrl = T_z(:,2);

PNNLCodes/CoSim/mpc_run.py

@@ -59,14 +59,14 @@ def write_port_file(port,host):
 
 price_n=np.loadtxt('paras/price_g.csv')
 
-control_horizon_length=3600
+control_horizon_length=int(sys.argv[1])
 
-samp_time=3600
+samp_time=int(sys.argv[2])
 
 
-ele=[float(x) for item in price_e for x in repeat(item, 3600/samp_time)]
+ele=[float(x) for item in price_e for x in repeat(item, 3600/control_horizon_length)]
 
-gas=[float(x) for item in price_n for x in repeat(item, 3600/samp_time)]
+gas=[float(x) for item in price_n for x in repeat(item, 3600/control_horizon_length)]
 
 
 print len(ele)
@@ -76,7 +76,7 @@ def write_port_file(port,host):
 prediction=pd.read_csv('paras/input.csv' )
 
 
-prediction=prediction.groupby(np.arange(len(prediction))//int(samp_time/60)).mean()
+prediction=prediction.groupby(np.arange(len(prediction))//int(control_horizon_length/60)).mean()
 
 tout=prediction['tout'].values.tolist()
 
@@ -99,6 +99,8 @@ def write_port_file(port,host):
 
 
 record=0
+tini=[]
+tini_arr=[]
 while 1:
 
 
@@ -126,7 +128,22 @@ def write_port_file(port,host):
               mssg = '%r %r %r 0 0 %r' % (vers, flag, 38, time) 
               if record<=1439:
                    tset=[[22]]*16
-
+
+
+              if record>0 and record%(int(samp_time/60))==0:
+                     tini=[]                    
+                     for i in range(16):
+                              sum=0
+                              for j in range(len(tini_arr)):
+                                     sum=sum+tini_arr[j][i]
+
+                              tini.append(sum/float(len(tini_arr)))
+                     tini_arr=[]
+              temp=[]
+              for i in range(7,7+16):
+                                  temp.append(float(arry[i]))
+
+              tini_arr.append(temp)  
 #              print mssg
               if record>1439 and index%(int(control_horizon_length/60))==0:
 
@@ -145,9 +162,7 @@ def write_port_file(port,host):
 
                        gas_m=matlab.double(gas[start_t:end_t])
 
-                       tini=[]
-                       for i in range(7,7+16):
-                                  tini.append(float(arry[i]))                                       
+
                        tini_m=matlab.double([tini])         
 
 
@@ -172,7 +187,7 @@ def write_port_file(port,host):
                        set=[]
                        m=[]
                        rh=[]
-                       tset,m,rh= eng.func_EDC_CoSim_test(samp_time/60.,samp_time/60.,ele_m,gas_m,tout_m,qin_m,tini_m, nargout=3)
+                       tset,m,rh= eng.func_EDC_CoSim_test(control_horizon_length/60.,samp_time/60.,ele_m,gas_m,tout_m,qin_m,tini_m, nargout=3)
                        for tt in range(len(m)-1):
                              f1.writelines(str(m[tt][0])+',')
                        f1.writelines(str(m[-1][0])+'\n')