复现文章：

微电网两阶段鲁棒优化经济调度方法——刘一欣（中国电机工程学报）

主要内容：

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

部分程序：

%燃气轮机参数设置
pg_max=800; %燃气轮机最大功率限制
pg_min=80; %燃气轮机最小功率限制
a=0.67; %燃气轮机成本系数a,b设置
b=0;

%蓄电池参数设置
ps_max=500; %储能允许最大充放电给功率
Es_max=1800; %蓄电池调度过程中允许的最大剩余容量
Es_min=400; %蓄电池调度过程中允许的最小剩余容量
Es_0=1000; %调度过程中初始容量
Ks=0.38; %折算后充放电成本
yita=0.95; %充放电效率

%需求响应负荷参数设置
K_DR=0.32; %需求响应负荷单位调度成本
D_DR=2940; %需求响应总用电需求
D_DR_min=50; %需求响应用电需求最大值
D_DR_max=200; %需求响应用电需求最小值

%配电网交互功率参数设置
pm_max=1500; %微电网与配电网交互功率最大值

%配电网日前交易电价，为24*1向量
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];

%光伏日前预测，为24*1向量
p_pv_forecast_0=1500*[ 0 0 0 0 0 0.0465 0.1466 0.3135 0.4756 0.5213 0.6563 1.0000 0.7422 0.6817 0.4972 0.4629 0.2808 0.0948 0.0109 0 0 0 0 0]';
%p_pv_forecast=[0; 0; 0; 0; 0; 0; 40; 200; 425; 731; 884; 1180; 900; 830; 600; 510; 340; 50; 0; 0; 0; 0; 0; 0]; %初始最坏数据
%负荷日前预测，为24*1向量
p_l_forecast_0=1500*[ 0.4658 0.4601 0.5574 0.5325 0.5744 0.6061 0.6106 0.6636 0.7410 0.7080 0.7598 0.8766 0.7646 0.7511 0.6721 0.5869 0.6159 0.6378 0.6142 0.6752 0.6397 0.5974 0.5432 0.4803]';
%p_l_forecast=[400; 350; 320; 300; 300; 310; 451; 561; 605; 748; 720; 810; 891; 836; 770; 726; 775.5; 730; 790; 810; 850; 800; 505; 410]; %初始最坏数据
%
C=[];
c=[a*ones(1,24) Ks*yita*ones(1,24) (Ks/yita)*ones(1,24) zeros(1,24) K_DR*ones(1,48) price' -price' zeros(1,48)];
%% 2.变量设置
%对偶变量设置（决策形函数）
gamma=sdpvar(192,1);
lamda=sdpvar(50,1);
miu=sdpvar(192,1);
pai=sdpvar(48,1);

%二元变量B设置
%B=binvar(48,1); %B为子问题初始二元变量，取到1即为最坏情况