适用平台：Matlab+Yalmip+Cplex

参考文献：

《考虑机组禁止运行区间的含风电鲁棒机组组合》-中国电机工程学报

《微电网两阶段鲁棒优化经济调度方法》-中国电机工程学报

程序提出了微电网中电源容量的两阶段鲁棒优化配置模型，第一阶段主要决策储能、风电、光伏的容量（容量配置），第二阶段主要规划风电、光伏、储能和燃气轮机的实际出力变量（功率平衡），最后给出了微网电源容量配置的结果和各个机组的出力。程序算例丰富、注释清晰、干货满满，可扩展性和创新性很高！码超所值！下面对文章和程序做简要介绍！

程序创新点：

1）由于所提模型为非线性，程序先引入整数变量将其转换为MILP模型，随后提出双层循环C&CG算法以求解此两阶段鲁棒问题，其中外层循环为传统的C&CG算法，内层循环为嵌套式列生成算法。

2）针对微电网内可再生能源和负荷的不确定性，建立了min-max-min结构的两阶段鲁棒优化能源配置模型，可得到最恶劣场景下微网不同电源的配置容量方案。

主要工作：

针对微电网内可再生能源和负荷的不确定性，建立了min-max-min结构的两阶段鲁棒优化模型，目的在于寻求不确定工况下微网中不同电源的最优配置容量。模型中考虑了储能、需求侧负荷及可控分布式电源等的运行约束和协调控制，并引入了不确定性调节参数，可灵活调整调度方案的保守性。基于列约束生成算法和强对偶理论，可将原问题分解为具有混合整数线性特征的主问题和子问题进行交替求解，从而得到原问题的最优解。最终通过仿真分析验证了所建模型和求解算法的有效性，同时给出了分时电价机制下微电网对储能进行调度的边界条件。

主要框架：

两阶段鲁棒模型

程序结果：

部分程序：

%% 开始运行
%先运行一次，得到UB-LB
[yita,LB,ee_bat_int, p_wt_int,p_pv_int,p_g_int] = MP;
[p_wt,p_pv,p_load,x,UB] = SP(ee_bat_int,p_wt_int,p_pv_int,p_g_int,LB,yita);
UB1 = UB;p(1)= UB - LB;
pub(1)=0;plb(1)=0;
%开始迭代
for k=1:10[yita,LB,ee_bat_int,p_wt_int,p_pv_int,p_g_int] = MP2(p_wt,p_pv,p_load);%MP迭代[p_wt,p_pv,p_load,x,UB] = SP(ee_bat_int,p_wt_int,p_pv_int,p_g_int,LB,yita);%SP迭代UB = min(UB1,UB);%取UB较小值pub(k+1)=UB;plb(k+1)=LB;p(k+1) = UB-LB;
%%绘图版块：主要绘制了各微网的日运行计划，容量配置结果，迭代过程等等
figure(1)plot(x(1:24),'-*')
xlim([1 24])grid
hold on plot(x(25:48),'-*')
bar(x(49:72))plot(x(73:96),'-d')
plot(x(97:120),'-d')
title('典型日1场景下微网运行计划')

部分内容源自网络，侵权联系删除！

欢迎感兴趣的小伙伴关注并获取完整版代码，小编会不定期更新高质量的学习资料、文章和程序代码，为您的科研加油助力！