拓扑网络连通性算法.docVIP

网络连通性算法 网络定义 节点与支路的集合，该集合中的节点与支路的连接关系可通过一节点-节点关联矩阵A充分表达： A=[aij]n×n i,j=1,2,…,n 式中：aij= n—网络节点数 连通性算法 理论算法： 称矩阵A为网络一级连通矩阵，A2为二级连通矩阵，…，An-1为n-1级连通矩阵。 A2=AA=[a2ij]n×n i,j=1,2,…,n 式中：a2ij= k=1,2,…,n，k≠i,j …… An-1= =[an-1ij]n×n i,j=1,2,…,n 式中： an-1ij= 矩阵An-1的每一线性无关的行或列中“1”元素对应的节点均处于同一连通子集中。 实际算法： 若矩阵A第i（i=1,2,…,n）行元素与第j（j=i+1,i+2,…,n）行元素中第k列元素aik和ajk同为“1”，则第j行中的其它“1”元素均填入第i行的相应列中。结果矩阵A第i行中所有“1”元素对应的节点处于同一连通子集中。 数据定义 Nc—元件数 Nd—节点数 NOD（Nc,3）—每个元件的节点编号i、j、k KND（Nc）—每个元件的种类（断路器、隔离开关、母线、线路、变压器……） CNT（Nc）—每个开关元件的分、合状态（逻辑型，例如：合为“真”，分为“假”） NDS0（Nd）—每个节点初始所在连通子集编号 NDS（Nd）—每个节点所在连通子集编号 NCT0（Nc）—每个元件初始所在连通子集编号 NCT（Nc）—每个元件所在连通子集编号 NST（Ns，3）—每个原始连通子集内[子集号，子集内节点数，子集内首位节点号] Ns—最大可能连通子集数 RA（Nd）—节点关联矩阵第i行，逻辑型 RB（Nd）—节点关联矩阵第j行，逻辑型 检验第k0个连通子集的连通性子程序CNTS(k0) 初始化 IND=NST(k0,3) 取第k0个连通子集的首位节点号 N=0 连通子集数置0 LOOP1 l=1，Nd l从1至Nd循环 IF (NDS0（l）=k0)， NDS（l）=0 第k0个连通子集的节点l的子集号临时置0 END LOOP1 连通性检验大循环 10 NSUM=1 节点关联矩阵“真”元素计数置为1 LOOP1 l=1，Nd l从1至Nd循环 RA(l)=FALSE 第k0个连通子集中第IND行第l列元素置为“假” END LOOP1 RA(IND)=TRUE 节点关联矩阵第IND行对角元素置1 M=1 节点关联矩阵第IND行 “真”元素计数置为1 形成节点关联矩阵的第IND行RA LOOP1 l=1,Nc l从1至Nc循环 IF(NCT0(l)=k0)，THEN 如果元件l属于初始连通子集k0，则 IF(KND(l)≠‘开关’or (KND(l)=‘开关’and CNT(l)=‘合’)，THEN I=NOD(l,1) J=NOD(l,2) K=NOD(l,3) 取元件l的各端节点号 IF(I=IND)，THEN 如果节点号I等于节点号IND，则 IF(J≠0 and J≠I)，THEN 如果节点号J不等于0和I，则 RA(J)=TRUE 关联矩阵第IND行第J列元素置为“真” M=M+1 关联矩阵第IND行“真”元素计数+1 END IF IF(K≠0 and K≠I and K≠J)，THEN 如果节点号K不等于0和I和J，则 RA(K)=TRUE 关联矩阵第IND行第K列元素置为“真” M=M+1 关联矩阵第IND行“真”元素计数+1 END IF END IF IF(J=IND)，THEN 如果节点号J等于节点号IND，则 IF(I≠0 and I≠J)，THEN 如果节点号I不等于0和I，则 RA(I)=TRUE 关联矩阵第IND行第I列元素置为“真” M=M+1 关联矩阵第IND行“真”元素计数+1 END IF IF(K≠0 and K≠I and K≠J)，THEN 如果节点号K不等于0和I和J，则 RA(K)=TRUE 关联矩阵第IND行第K列元素置为“真” M=M+1 关联矩阵第IND行“真”元素计数+1 END IF END IF IF(K=IND)，THEN 如果节点号K等于节点号IND，则 IF(I≠0 and I≠K)，THEN 如果节点号I不等于0和K，则 RA(I)=TRUE 关联矩阵第IND行第I列元素置为“真”