目录
1主要内容
1.1 目标函数
1.2 约束条件
1.3 不确定变量
2部分代码
3程序效果
4 程序链接
1主要内容
程序主要复现的是《考虑寿命损耗的微网电池储能容量优化配置》，解决微网中电源/储能容量优化配置的问题，即风电、光伏、储能以及燃气轮机的容量规划问题，程序考虑了不确定性，采用两阶段鲁棒优化方法进行优化，采用的是KKTCCG算法，一阶段主要决策储能、风电、光伏的容量，二阶段主要规划风光燃储的实际出力变量，最后结果不仅给出了微网电源容量配置的结果，还给出了各个机组的出力。
1.1 目标函数
1.2 约束条件
1.3 不确定变量
不确定变量考虑风光出力和负荷的不确定性，采用盒式不确定的方式。
tips:
1.该程序虽然收敛性良好，但是两阶段鲁棒优化流程还不完善，在一阶段未实现不同迭代次数下变量叠加，只是实现了变量的更新！有兴趣的可以看一下比较规范的两阶段鲁棒优化程序：《6节点电网两阶段鲁棒优化调度matlab》《【重磅】IEEE33配电网两阶段鲁棒优化调度CCG》
2.模型对偶采用kkt方法，还是值得学习的，可结合《Solving two-stage robust optimization problems using a column-and-constraint generation》进行学习。
2代码
%% 1.设参%投资成本参数rp = 0.08;%折现率rbat = 10;rPV = 20;rWT = 15;rG = 15;%折现年数cbat = 1107;cPV = 100;cWT = 300;cG = 2000;%单位容量投资成本pm_max = 500;%联络线功率上限eta = 0.95;%储能充放电效率c_wt_om = 0.0296;c_pv_om = 0.0096;c_g_om = 0.059;c_bat_om = 0.009;%运维成本系数c_fuel = 0.6;%燃料成本系数%% 2.设决策变量p_ch = sdpvar(24,4);%储能充电p_dis = sdpvar(24,4);%储能放电uu_bat = binvar(24,4);%充放电标识uu_m = binvar(24,4);p_buy = sdpvar(24,4);%配网购电p_sell = sdpvar(24,4);%配网售电p_wt = sdpvar(24,4);p_pv = sdpvar(24,4);p_load = sdpvar(24,4);p_g = sdpvar(24,4);%微型燃气轮机%% 3.设变量ee_bat_int = sdpvar(1);%储能容量上限p_pv_int = sdpvar(1);p_wt_int = sdpvar(1);p_g_int = sdpvar(1);yita = sdpvar(1);p_bat_int = ee_bat_int*0.21;%假设储能的功率上限和容量上限有比值关系ee0 = 0.55*ee_bat_int;%储能初始电量%风光出力和电价（以春季典型日为例）p_l = xlsread('四个典型日数据.xlsx','0%','B3:E26')*900;%？max_p_wt = xlsread('四个典型日数据.xlsx','0%','H3:K26')*p_wt_int; max_p_pv = xlsread('四个典型日数据.xlsx','0%','N3:Q26')*p_pv_int; %price=xlsread('四个典型日数据.xlsx','电价','A2:A25');price = [0.48;0.48;0.48;0.48;0.48;0.48;0.48;0.9;1.35;1.35;1.35;0.9;0.9;0.9;0.9;0.9;0.9;0.9;1.35;1.35;1.35;1.35;1.35;0.48];%% 4.设约束C = [];load = p_l';wwt = 0.05;wpv = 0.1;wl = 0.15;%不确定度，缩放比例C = [C, (1 - wwt)*max_p_wt <= p_wt,p_wt <= (1 + wwt)*max_p_wt];%不确定性风C = [C, (1 - wpv)*max_p_pv <= p_pv,p_pv <= (1 + wpv)*max_p_pv];%不确定性光C = [C, (1 - wl)*load' <= p_load,p_load <= (1 + wl)*load'];%不确定性负荷%储能线性化
3程序效果
4 程序链接
这里是《电力程序》集中地，单击直达！
