开源代码分享(19)-配电网孤岛优化划分方法matlab代码

function sx=fig(Zij) % AA=6; A1=[1,3];A2=[2,3];A3=[3,3];A4=[4,3];A5=[5,3];A6=[6,3]; A7=[7,3];A8=[8,3];A9=[9,3];A10=[10,3];A11=[11,3];A12=[12,3]; A13=[13,3];A14=[14,3];A15=[15,3];A16=[16,3];A17=[17,3];A18=[18,3]; A19=[2,1];A20=[3,1];A21=[4,1];A22=[5,1];A23=[3,5];A24=[4,5]; A25=[5,5];A26=[6,5];A27=[7,5];A28=[8,5];A29=[9,5];A30=[10,5]; A31=[11,5];A32=[12,5];A33=[13,5];A34=[9,2];A35=[15,2];A36=[5,6];A37=[9,6];A38=[12,1]; for p = 1:38 % 放进一个矩阵内方便操作 c = num2str(p); s = ['A(p,:) = A' c ';']; eval(s); end % scatter(A(:,1),A(:,2)); for p = 1:33 c = num2str(p); plot(A(p,1),A(p,2),'ko','MarkerFaceColor','k','LineWidth',1); hold on axis([0 20 0 7]); text(A(p,1)+0.05,A(p,2)-0.1,c); end set(gca,'xtick',0:1:20,'ytick',0:1:7); sx = 1:18; for p = 1:length(sx)-1 P1 = A(sx(p),:); P2 = A(sx(p+1),:); line([P1(1) P2(1)],[P1(2) P2(2)],'color','k','LineWidth',1); end sx = [2 19 20 21 22]; for p = 1:length(sx)-1 P1 = A(sx(p),:); P2 = A(sx(p+1),:); line([P1(1) P2(1)],[P1(2) P2(2)],'color','k','LineWidth',1); end sx = [3 23 24 25]; for p = 1:length(sx)-1 P1 = A(sx(p),:); P2 = A(sx(p+1),:); line([P1(1) P2(1)],[P1(2) P2(2)],'color','k','LineWidth',1); end sx = [6 26 27 28 29 30 31 32 33]; for p = 1:length(sx)-1 P1 = A(sx(p),:); P2 = A(sx(p+1),:); line([P1(1) P2(1)],[P1(2) P2(2)],'color','k','LineWidth',1); end % plot(A(AA,1)+0.5,A(AA,2),'p','MarkerSize',10,'MarkerFaceColor','r','MarkerEdgeColor','r','LineWidth',1.5); sx = [21 8];%33 for p = 1:length(sx)-1 P1 = A(sx(p),:); P2 = A(sx(p+1),:); line([P1(1) P2(1)],[P1(2) P2(2)],'color','k','linestyle','-','LineWidth',1); end sx = [9 34 35 15];%34 for p = 1:length(sx)-1 P1 = A(sx(p),:); P2 = A(sx(p+1),:); line([P1(1) P2(1)],[P1(2) P2(2)],'color','k','linestyle','-','LineWidth',1); end sx = [12 38 22];%35 for p = 1:length(sx)-1 P1 = A(sx(p),:); P2 = A(sx(p+1),:); line([P1(1) P2(1)],[P1(2) P2(2)],'color','k','linestyle','-','LineWidth',1); end sx = [18 33];%36 for p = 1:length(sx)-1 P1 = A(sx(p),:); P2 = A(sx(p+1),:); line([P1(1) P2(1)],[P1(2) P2(2)],'color','k','linestyle','-','LineWidth',1); end sx = [25 36 37 29];%37 for p = 1:length(sx)-1 P1 = A(sx(p),:); P2 = A(sx(p+1),:); line([P1(1) P2(1)],[P1(2) P2(2)],'color','k','linestyle','-','LineWidth',1); end grid on for p=1:17 if Zij(p)==0 sx=1:18; P1 = A(sx(p),:); P2 = A(sx(p+1),:); line([P1(1) P2(1)],[P1(2) P2(2)],'color','r','LineWidth',1); end %plot(A(jieru(p),1),A(jieru(p),2),'p','MarkerSize',10,'MarkerFaceColor','r','MarkerEdgeColor','r','LineWidth',1.5); end for p=1:4 sx = [2 19 20 21 22]; if Zij(p+17)==0 P1 = A(sx(p),:); P2 = A(sx(p+1),:); line([P1(1) P2(1)],[P1(2) P2(2)],'color','r','LineWidth',1); end end for p=1:3 sx = [3 23 24 25]; if Zij(p+21)==0 P1 = A(sx(p),:); P2 = A(sx(p+1),:); line([P1(1) P2(1)],[P1(2) P2(2)],'color','r','LineWidth',1); end end for p=1:8 sx = [6 26 27 28 29 30 31 32 33]; if Zij(p+24)==0 P1 = A(sx(p),:); P2 = A(sx(p+1),:); line([P1(1) P2(1)],[P1(2) P2(2)],'color','r','LineWidth',1); end end if Zij(34)==0 sx = [21 8];%33 for p = 1:length(sx)-1 P1 = A(sx(p),:); P2 = A(sx(p+1),:); line([P1(1) P2(1)],[P1(2) P2(2)],'color','r','LineWidth',1); end end if Zij(35)==0 sx = [9 34 35 15];%34 for p = 1:length(sx)-1 P1 = A(sx(p),:); P2 = A(sx(p+1),:); line([P1(1) P2(1)],[P1(2) P2(2)],'color','r','LineWidth',1); end end if Zij(33)==0 sx = [12 38 22];%35 for p = 1:length(sx)-1 P1 = A(sx(p),:); P2 = A(sx(p+1),:); line([P1(1) P2(1)],[P1(2) P2(2)],'color','r','linestyle','-','LineWidth',1); end % text(A(22,1)+2.2,A(22,2)+0.2,'SOP'); end if Zij(37)==0 sx = [18 33];%36 for p = 1:length(sx)-1 P1 = A(sx(p),:); P2 = A(sx(p+1),:); line([P1(1) P2(1)],[P1(2) P2(2)],'color','r','LineWidth',1); end end if Zij(36)==0 sx = [25 36 37 29];%37 for p = 1:length(sx)-1 P1 = A(sx(p),:); P2 = A(sx(p+1),:); line([P1(1) P2(1)],[P1(2) P2(2)],'color','r','LineWidth',1); end end hold off axis off title('故障处理网络重构图');