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];
