subcaculating.m

基于长期负荷预测的线性规划燃机出力能量优化控制

function [x,exitflag,feval]=subcaculating(a_p,a_q,b_p,b_q,i,p_ng,p_ele,COP_ac,COP_tes,COP_dis,COP_off,LHV,PER_rate,sta)

f1=3.6/LHV*p_ng;f2=p_ele;f3=0.77;

Aeq=[1 a_p 0 -1 0                     %无综合楼电空调耗电 新大楼空调挂在微网母线
    0 0 COP_dis COP_off 1];
% A=[0 PER_rate-COP_ac*a_q-a_p 0 0];
% %    0 -COP_ac*a_q 0 0 1];
% b=[COP_ac*b_q+b_p];

[E_T,E_D,E_O,Q_T,Q_D,Q_O]=creat(sta);
E_tes=E_T(i); E_dis=E_D(i); E_off=E_O(i);% P_gt=0.4*E_gt-40;   %代入可得     0.4*E_gt+P_ep-E_re_dis-E_re_off=P_loads+E_re_tes+40
Q_tes=Q_T(i); Q_dis=Q_D(i); Q_off=Q_O(i); % Q_gt=0.3*E_gt+40;   %代入可得    E_re_dis*COP_dis+E_re_off*COP_off+Q_ac=Q_dis+Q_off

beq1 = E_tes+Q_tes/COP_tes+E_dis*0.5-b_p;
beq2 = Q_dis+Q_off;
beq = [beq1;beq2];

f=[f3;f1;f2;0;0];
lb = [-inf;ceil((80-b_p)/a_p);0;0;0];
ub = [0;1800;Q_dis/COP_dis;Q_off/COP_off;1092.5];
[x,feval,exitflag,output]=linprog(f,[],[],Aeq,beq,lb,ub);
% x;