clear;clc;
tic
T=24;
J=7; %cchp7种方案冷热电联供
K=10; %gb 10种方案燃气锅炉
M=10; %ac 10种方案制冷机
P1=[12.50 24.49 36 47.06 57.69 67.92 77.78];%CCHP最大多能输入
P2=[10.87 21.74 32.61 43.48 54.35 65.22 76.09 86.96 97.83 108.70];%gb最大多能输入
P3=[2.5 5 7.5 10 12.5 15 17.5 20 22.5 25];%AC最大多能输入电
Price=xlsread('D:\参考文献\算例\1.1\data.xlsx','sheet3','B4:Y4'); %分时电价
gas=2.5; % 气价
c1=xlsread('D:\参考文献\算例\1.1\data.xlsx','sheet2','B3:H3'); %cchp建设成本
c2=xlsread('D:\参考文献\算例\1.1\data.xlsx','sheet2','B11:K11');%gb建设成本
c3=xlsread('D:\参考文献\算例\1.1\data.xlsx','sheet2','B19:K19');%ac建设成本
P_load=xlsread('D:\参考文献\算例\1.1\data.xlsx','sheet1','A1:A24');%电负荷
C_load=xlsread('D:\参考文献\算例\1.1\data.xlsx','sheet1','B1:B24');%冷负荷
H_load=xlsread('D:\参考文献\算例\1.1\data.xlsx','sheet1','C1:C24');%热负荷
xcchp=binvar(1,J,'full');
xgb=binvar(1,K,'full');
xac=binvar(1,M,'full');
Pcchp=sdpvar(1,T,'full');
Pgb=sdpvar(1,T,'full');
Pac=sdpvar(1,T,'full');
Psub=sdpvar(1,T,'full');
P_load=P_load';
C_load=C_load';
H_load=H_load';
%约束条件
st=[];
st=[0<=Pcchp<=xcchp*P1'
sum(xcchp)<=1
0<=Pgb<=xgb*P2'
sum(xgb)<=1
0<=Pac<=xac*P3'
sum(xac)<=1
0<=Psub<=30
];
%多能负荷平衡约束
for t=1:T
st=[st (P_load(t)+Pac(t))<=(0.4.*Pcchp(t)+Psub(t))
C_load(t)<=4.*Pac(t)+0.432.*Pcchp(t)
H_load(t)<= 0.92.*Pgb(t)+0.324.*Pcchp(t)
];
end不知道你这个问题是否已经解决, 如果还没有解决的话: