Battery storage -optimisation problem_charge and discharge

17 次查看（过去 30 天）

NN 2021-8-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/891877-battery-storage-optimisation-problem_charge-and-discharge

评论： Walter Roberson 2021-8-22

I have seperate decision variables for charge and discharge cosntraints for battery optimisation problem, but i am getting proper value as optimised outputs.

Below is the code :

% initial conditions, battery constraints 
batteryMinMax.Pmin = -500e3;    
batteryMinMax.Pmax = 500e3; 
battEnergy = 3.6e6*upsBattCap;
batteryMinMax.Emax = 0.8*battEnergy;
batteryMinMax.Emin = 0.2*battEnergy;
%Decision variables
PbattchV = optimvar('PbattchV',N,'LowerBound',0,'UpperBound',batteryMinMax.Pmax);
EbattchV = optimvar('EbattchV',N,'LowerBound',0,'UpperBound',batteryMinMax.Emax);
PbattDchV = optimvar('PbattDchV',N,'LowerBound',0,'UpperBound',batteryMinMax.Pmin);
EbattDchV = optimvar('EbattDchV',N,'LowerBound',0,'UpperBound',batteryMinMax.Emin);

Please advice if i made any mistake.

13 个评论
显示 11更早的评论隐藏 11更早的评论

NN 2021-8-18

在 MATLAB Online 中打开

Please see the below code :It returns a value of zero ,and i am expecting and output of 500e3 for Pbattups.Please advice.

function [Pd,Pbattups,Ebattups] = battSolarOptimizetrial3(N,dt,Ppv,Pload,Einit,Cost,FinalWeight,upsbatteryMinMax)
prob = optimproblem;
% Decision variables
PdV = optimvar('PdV',N);
PbattupsV = optimvar('PbattupsV',N,'LowerBound',0,'UpperBound'500e3);
EbattupsV = optimvar('EbattupsV',N,'LowerBound',upsbatteryMinMax.Emin,'UpperBound',upsbatteryMinMax.Emax);
% Maximise the profit from the usage of battery and PV
prob.ObjectiveSense = 'minimize';
%prob.Objective = var1 -dt*DayaheadCost'*Ppv-FinalWeight*EbattupsV(N)-FinalWeight*Ebatt1V(N)-FinalWeight*Ebatt2V(N);
prob.Objective =dt*Cost'*PbattupsV- dt*Cost'*Ppv-FinalWeight*EbattupsV(N);
% Power input/output to battery
prob.Constraints.energyBalance = optimconstr(N);
prob.Constraints.energyBalance(1) = EbattupsV(1) == Einit;
prob.Constraints.energyBalance(2:N) = EbattupsV(2:N) == EbattupsV(1:N-1) - PbattupsV(1:N-1)*dt
% Satisfy power load with power from PV, grid and battery
prob.Constraints.loadBalance = Ppv + PdV + PbattupsV == Pload;
%prob.Constraints.battpower = PbattupsV<=500e3;
% Solve the linear program
options = optimoptions(prob.optimoptions,'Display','none');
[values,~,exitflag] = solve(prob,'Options',options);
% Parse optmization results
if exitflag <= 0
    Pd = zeros(N,1);
    Pbattups = zeros(N,1);
    Ebattups = zeros(N,1);
else
    Pd = values.PdV;
    Pbattups = values.PbattupsV; 
    Ebattups = values.EbattupsV;
end