微❤关注“电气仔推送”获得资料（专享优惠）
参考文档：
A New Model for Resilient Distribution Systems by Microgrids Formation；
主动配电网故障恢复的重构与孤岛划分统一模型；
同时考虑孤岛与重构的配电网故障恢复运行策略。
运行环境：
matlab+yalmip+cplex
主要内容：
程序主要以33节点为研究对象，编制配电网故障重构模型，采用虚拟潮流约束保证配电网的连通性和辐射性，以网损和负荷损失作为目标函数，包括潮流约束、电压电流约束、辐射性约束等。
主动配电网故障恢复的重构与孤岛划分统一模型 随着分布式电源和储能装置大量接入配电网，配电网 在发生故障之后可以进行重构和孤岛划分，提高故障恢复的水平。建立包含多类型分布式电源、柔性负荷和储能的多 时间段故障动态恢复模型，并计及故障恢复时间和检修次序。利用二阶锥技术将 所建模型转换成混合整数二阶锥规划模型。结果表明，与传统方法相比，所提故障恢复方法计算速度更快，故障恢复率更高。
部分代码：
%% 检修策略约束
for i=1:nl
if ismember(i,F)
Constraints = [Constraints, Zij(i,1)==0];%故障线路初始状态
for t=2:N
Constraints = [Constraints, Zij(i,t)>=Zij(i,t-1)];%故障线路下一个时刻状态要大于等于上一个时刻
end
end
end
[~,na]=size(F);%故障线路条数
for t=1+Tf:N
Constraints = [Constraints, sum(Zij(F,t))-sum(Zij(F,t-Tf))<=h];%每个时段最多只能检修h条故障线路
end
for i=1:6
for t=2:N
Rep(i,t-1)=Zij(F(i),t)-Zij(F(i),t-1);%线路维修状态和线路状态关系
end
%% 定义目标函数
f=sum(sum(repmat(r_ij,1,N).*x_Iij_square))+100*sum(sum(p_Load1))+100*sum(sum(-lamda.*p_Load1));
Constraints=[Constraints,0<=lamda<=1];
%% 有功与无功平衡
for opt_num=1:N
Constraints=[Constraints,x_ui_square(1,opt_num)==12.66^2];
for k=2:33
Constraints = [Constraints, -0.01.*Zij(33:37,opt_num)<= x_pij(33:37,opt_num) <= 0.01.*Zij(33:37,opt_num)];%线路有功约束
Constraints = [Constraints, -4.*Zij(1:32,opt_num) <= x_pij(1:32,opt_num) <= 4.*Zij(1:32,opt_num)];%线路有功约束
%% 支路欧姆定律
%欧姆定律约束
m = (1.05*1.05 - 0.95*0.95)*12.66^2;
M = (ones(nl,N) - Zij)*m;
%% 运行电压电流约束
Constraints=[Constraints,x_ui_square(:,opt_num)>=(0.95*12.66)^2];
Constraints=[Constraints,xuisquare(:,opt_num)<=(1.05*12.66)^2];
运行结果：