H = A / B,其中H,A和B都是决策变量。H 是要求解的变量,A 和 B 是由连续决策变量构造的。A=HW*SW*sum(k(i,m))+HG*SG,B=SW*sum(k(i,m))+SG)。其中,SW、SG、HW、HG都是常数项,k(i,m)是由0-1决策变量构建的。
% 稳态频率约束
for m=2:NT
for i=1:NC
st=st+[
(SW*sum(k(i,m))+PL*sum(n(i,m))+SG)*K(1,m)==KW*SW*sum(k(i,m))+KL*PL*sum(n(i,m))+KG*SG;
P1(1,m)<=fmax*K(1,m)*(SW*sum(k(i,m))+SG);
P2(1,m)<=fmax*K(1,m)*SW*sum(k(i,m))];
end
end
%对每一时步的01决策变量进行计数
for m=1:NT
for i=1:NC
st=st+[n(i,m)==v(i,m)*NE(i,1);
k(i,m)==u(i,m)*NE(i,1)];
end
end
NE=[5;5;5;6;6];
我通过将H=A/B转换为HB=A进行代码编写,得到的结果H为0,导致我的约束出现矛盾。