1 主要内容
程序以IEEE33节点为例，分析含sop的配电网优化，包括sop有功约束、无功约束和容量约束，非线性部分通过转换为旋转锥约束进行编程，并且包括33节点配电网潮流及对应电压、电流约束、OLTC约束等，程序是对《Coordinated Control Method of Voltage and Reactive Power for Active Distribution Networks Based on Soft Open Point》方法实现，里面还有更深一些的知识，包括OLTC和CB等部分内容，具体sop部分和文献《考虑分布式电源运行特性的有源配电网智能软开关sop规划方法》一致。
2 部分代码
%% 定义约束条件 Constraints=[]; %% SOP有功功率限制 (1) Constraints=[Constraints,x_p_sop1(1,:)+x_p_sop1(2,:)+x_p_sop1_loss(1,:)+... x_p_sop1_loss(2,:)==0]; Constraints=[Constraints,x_p_sop2(1,:)+x_p_sop2(2,:)+... x_p_sop2_loss(1,:)+x_p_sop2_loss(2,:)==0]; %% CBs操作限制(27、29) 33 Constraints=[Constraints,Q_CB==N_CB*q_CB]; Constraints=[Constraints,N_CB<=5]; Constraints=[Constraints,N_CB>=0]; %% 有功功率损耗（33） f_loss1=sum(r_ij'*x_Iij_square)*delta_T; Constraints=[Constraints,x_ui_square(1,:)==12.66^2]; %平衡节点每小时电压平方 f_loss2=sum(x_p_sop1_loss)*delta_T; f_loss3=sum(x_p_sop2_loss)*delta_T; Constraints=[Constraints,x_Iij_square>=0]; %% sop运行约束%（40,41） Constraints=[Constraints,0.02*sqrt(x_p_sop1(1,:).^2+x_q_sop1(1,:).^2)<=... x_p_sop1_loss(1,:)]; %Constraints=[Constraints,(x_p_sop1(2,opt_num)^2+x_q_sop1(2,opt_num)^2)<=... % 2*(x_p_sop1_loss(2,opt_num)/(sqrt(2)*0.02))*(x_p_sop1_loss(2,opt_num)/(sqrt(2)* 0.02))]; Constraints=[Constraints,0.02*sqrt(x_p_sop1(2,:).^2+x_q_sop1(2,:).^2)<=... x_p_sop1_loss(2,:)]; % Constraints=[Constraints,x_p_sop1(1,opt_num)+x_p_sop1(2,opt_num)==0]; Constraints=[Constraints,0.02*sqrt(x_p_sop2(1,:).^2+x_q_sop2(1,:).^2)<=... x_p_sop2_loss(1,:)]; Constraints=[Constraints,0.02*sqrt(x_p_sop2(2,:).^2+x_q_sop2(2,:).^2)<=... x_p_sop2_loss(2,:)]; % Constraints=[Constraints,x_p_sop2(1,opt_num)+x_p_sop2(2,opt_num)==0];
程序代码标注了和文献公式的对应性，方便学习！
3 程序结果
