matlab代码翻译成pythonG_Ebalance=[]; for j=1:OperationPeriod G_Ebalance=G_Ebalance+[sum(EG(:,j))+sum(Data.ESDisch(:,j)) + sum(PV(:,j)) == sum(Data.PL(:,j))+sum(Data.ESCh(:,j))]; end %% 约束2：电源出力上下限 %% G_G=[]; for i = 1 : NumG for j = 1 : OperationPeriod G_G=G_G+[Data.G(i,1)<=EG(i,j)]; end end %% 约束3：光伏出力约束 %% G_PV=[]; for i = 1 : NumPV for j = 1 : OperationPeriod if Data.PVForecast(i,j)<=Data.PL(i,j) G_PV = G_PV+[PV(i,j)==Data.PVForecast(i,j)]; else G_PV = G_PV+[PV(i,j)==Data.PL(i,j)]; end end end %% ------------------------目标函数------------------------------- %% %% 成本变量 %% CostG = 0; CostWT = 0; CostESdisch = 0; CostESch = 0; %% 电源成本 CostG = CostG+sum(EG.*Data.PriceE); %% 清洁能源成本 CostWT=CostWT+sum(sum(0.5*PV)); %% 弃光量 PVWAST = sum(sum(Data.PVForecast-PV)); %% 光电总量 PVF=0; PVF=PVF+sum(sum(Data.PVForecast)); %% 储能设备充电成本 %% CostESch = CostESch+sum(sum(ESCh).*Data.PriceE); %% 储能设备售电成本 %% CostESdisch = CostESdisch+sum(sum(ESDisch).*Data.PriceE); %% 分时成本 C2=sum(ESCh).*Data.PriceE-sum(ESDisch).*Data.PriceE; %% 总成本 %% + CostHeatBoiler + CostElectricBoiler + CO2 Cost = (CostESch-CostESdisch); G=[G_Ebalance+G_G+G_PV]; options =sdpsettings('verbose',2,'solver','CPLEX'); optimize(G,Cost,options); WYH.C2=value(C2); WYH.Cost=Cost; end

Translated code in Python:

G_Ebalance = []
for j in range(OperationPeriod):
    G_Ebalance = G_Ebalance + [sum(EG[:, j]) + sum(Data.ESDisch[:, j]) + sum(PV[:, j]) == sum(Data.PL[:, j]) + sum(Data.ESCh[:, j])]
    
G_G = []
for i in range(NumG):
    for j in range(OperationPeriod):
        G_G = G_G + [Data.G[i, 1] <= EG[i, j]]

G_PV = []
for i in range(NumPV):
    for j in range(OperationPeriod):
        if Data.PVForecast[i, j] <= Data.PL[i, j]:
            G_PV = G_PV + [PV[i, j] == Data.PVForecast[i, j]]
        else:
            G_PV = G_PV + [PV[i, j] == Data.PL[i, j]]

CostG = 0
CostWT = 0
CostESdisch = 0
CostESch = 0

CostG = CostG + sum(EG * Data.PriceE)

CostWT = CostWT + sum(sum(0.5 * PV))

PVWAST = sum(sum(Data.PVForecast - PV))

PVF = 0
PVF = PVF + sum(sum(Data.PVForecast))

CostESch = CostESch + sum(sum(ESCh) * Data.PriceE)

CostESdisch = CostESdisch + sum(sum(ESDisch) * Data.PriceE)

C2 = sum(ESCh) * Data.PriceE - sum(ESDisch) * Data.PriceE

Cost = (CostESch - CostESdisch)

G = G_Ebalance + G_G + G_PV

options = sdpsettings(verbose=2, solver='CPLEX')

optimize(G, Cost, options)

WYH.C2 = value(C2)
WYH.Cost = Cost

Note: This code translation is based on my understanding of the original code. Please make sure to confirm the correctness of the translated code before using it.