微❤关注"电气仔推送"获得资料（专享优惠）
参考文档：
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,x*ui*square(:,opt_num)\<=(1.05\*12.66)\^2\]; **运行结果：**    