caculating.m

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

function [X,EXITF,OPT,FEV1,FEV2]=caculating(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,FEV1,FEV2,E_tes,E_dis,E_off,Q_tes,Q_dis,Q_off)



f1=3.6/LHV*p_ng;f2=p_ele;
f=[f2;f1;f2;0;0];
EXITF=[0,0];OPT=[0,0];
X=zeros(length(f),2);
Aeq=[1 a_p 0 -1 0                    %无综合楼电空调耗电 新大楼空调挂在微网母线
    0 0 COP_dis COP_off 1];
A=[0 -COP_ac*a_q 0 0 1];%0 PER_rate-a_p 0 0 -1b_p;
b=[COP_ac*b_q];



beq1 = E_tes+Q_tes/COP_tes+E_dis*0.5-b_p;
beq2 = Q_dis+Q_off;
beq = [beq1;beq2];
lb = [0;ceil((80-b_p)/a_p);0;0;0];
ub = [inf;1800;Q_dis/COP_dis;Q_off/COP_off;1092.5];

[x,feval,exitflag,output]=linprog(f,A,b,Aeq,beq,lb,ub);
X(:,1)=x;
EXITF(1)=exitflag;
OPT(1)=feval;
FEV1(i)=feval;

f3=0.77;
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,A,b,Aeq,beq,lb,ub);
X(:,2)=x;
EXITF(2)=exitflag;
OPT(2)=feval;
FEV2(i)=feval;