1 主要内容
之前分享了很多配电网重构的程序，每个程序针对场景限定性比较大，程序初学者修改起来难度较大，本次分享一个基础程序，针对含sop的配电网重构模型，含风电和光伏，优化了33节点网络电压合理性，任意修改故障线路不仅得到满足辐射和联络性网络重构结果，也能得到完美的节点电压图。
采用IEEE33节点系统，光伏接入7、13、27节点，风机接于5节点，但是由于考虑单时间段（1时刻），光伏出力为0，风机功率为0.8MW，单线路故障得到如下结果图：
10线路故障：
19线路故障：
多线路故障结果：
5、31线路故障：
2 部分程序
Branch=[Branch;xb];
r_ij=Branch(:,4); x_ij=Branch(:,5);
% p_Solar=zeros(33,1);
p_Wind=zeros(33,1);
p_Wind(5,1)=0.8;%5节点加风电
p_Load=p_load;
% p_Solar(5)=0.05;p_Solar(11)=0.05;p_Solar(15)=0.05;p_Solar(22)=0.05;p_Solar(31)=0.05;%分布式功率
q_Solar=p_Solar.*0.05;%无功按照有功的0.05倍
%% 定义决策变量
lamda=sdpvar(33,N,'full');
x_Iij_square=sdpvar(37,N,'full'); x_ui_square=sdpvar(33,N,'full');
x_pij=sdpvar(37,N,'full'); x_qij=sdpvar(37,N,'full');
x_p_sop1=sdpvar(2,N,'full'); x_q_sop1=sdpvar(2,N,'full');
opt_num=1;
beta=binvar(37,37,'full');
Zij=binvar(nl,1);%网架结构
lx=[5,31];%故障线路编号
% Z0=[ones(nl-nc,1);zeros(nc,1)];%初始拓扑
% assign(Zij,Z0);
% x_sop_L=sdpvar(2,N,'full');
%% 定义约束条件
Constraints=[];
%for opt_num=1:N
Constraints = [Constraints, sum(Zij) == 32];
Constraints = [Constraints, 0<= x_pij(33:37) <= 0.1.*Zij(33:37)];
%% 定义目标函数
f=r_ij'*x_Iij_square(:,opt_num)+100*sum(p_load)+100*sum(-lamda.*p_load);
Constraints=[Constraints,x_ui_square(1,opt_num)==12.66^2];%平衡节点电压
Constraints=[Constraints,0<=lamda<=1];%负荷消纳比例
Constraints=[Constraints,Zij(lx)==0];%设置故障线路
Constraints=[Constraints,Zij(35)==0];
Constraints=[Constraints,0<=x_Iij_square<=100*Zij];