信息.txt

用MATLAB编的配电网自动化的一些程序

d=[0 0 0 0 0 0 1 0 0             %ͼ3-1(b)p25由D求C
    0 0 1 0 1 0 0 0 0
    0 1 0 0 0 0 1 0 0
    0 0 0 0 0 0 1 0 0
    0 1 0 0 0 0 0 1 1
    0 0 0 0 0 0 0 0 1
    1 0 1 1 0 0 0 0 0
    0 0 0 0 1 0 0 0 0
    0 0 0 0 1 1 0 0 0];
m=[1 0 0 0 0 1 0 0 0];
t=[1 1 0 0 1 1 1 0 1];


m=[0 1 1 0 0 1 0 1 0 0 1 0 0]
d=[0 0 0 0 1 0 0 0 0 1 0 0 0;    %p12ͼ2-2连通系的分解
   0 0 0 0 0 0 0 0 1 0 0 0 0;
   0 0 0 0 0 0 0 0 0 0 0 1 0;
   0 0 0 0 0 0 0 0 0 0 0 1 0;
   1 0 0 0 0 0 0 0 1 0 1 0 0;
   0 0 0 0 0 0 0 0 0 1 0 0 0;
   0 0 0 0 0 0 0 0 0 0 0 1 0;
   0 0 0 0 0 0 0 0 0 0 0 0 1;
   0 1 0 0 1 0 0 0 0 0 0 0 0;
   1 0 0 0 0 1 0 0 0 0 0 0 1;
   0 0 0 0 1 0 0 0 0 0 0 0 0;
   0 0 1 1 0 0 1 0 0 0 0 0 0;
   0 0 0 0 0 0 0 1 0 1 0 0 0]
b=[0 0 0 0 1 0 0 0 0 1 0 1 0];
t=[0 1 1 0 1 1 0 1 0 1 1 1 0]


t=[1 1 1 0 1 1 1 0 0 0]
b=[0 0 1 0 1 0 1 0 0 0]
c=[0 0 0 0 0 0 1 0 0 0;    %ͼ2-4p19最小配电区域分解
   0 0 0 0 1 0 0 0 0 0;
   0 1 0 0 0 0 0 0 0 1;
   0 0 0 0 0 0 0 0 0 0;
   0 0 0 0 0 0 0 1 1 0;
   0 0 0 0 0 0 0 0 1 0;
   0 0 1 1 0 0 0 0 0 0;
   0 0 0 0 0 0 0 0 0 0;
   0 0 0 0 0 0 0 0 0 0;
   0 0 0 0 0 0 0 0 0 0]

z=[1,7;2,3;2,5;3,2;3,7;3,10;4,7;5,8;5,2;5,9;6,9;7,1;7,3;7,4;8,5;9,6;9,5;10,3]  %P19图2-4 一个典型配电网络及其区域
z=[11 6;6 13;13 12;12 14;14 7;7 9;7 16;16 3;3 8;3 17;12 5;5 4;4 15;4 1;5 10;10 2]  %P38图3-4配电网的分层
z=[6 10;10 13;13 8;10 1;1 5;5 11;5 9;9 2;3 12;12 7;12 4]   %p12图2-2一个具有两个连通系的信息矩阵

vp=[vp1;vp2;vp3];
max_1=max(abs(vp1))
for i=1:max_1
    [i2,j2]=find(abs(vp)==i);
    
j2=j2';fprintf('p[%d]',i)
disp(j2);



for i=1:c2
         if(c(a,i)==1&b(i)==0)
             h(i)=mn
             if(length(q)<=1)
             mn=mn+1
             end
             q=[i q]
         elseif(c(a,i)==1&b(i)==1)
             q=[i q]
             
         end
       
         
     end


while sign
  while sign
     [a1,a2]=find(c(a,:)==1)
     [a3,a4]=size(a2)
     for i=1:a4
         if(b(a2(i))==0)
             h(a2(i))=mn
             
         elseif(b(a2(i))==1)
             a=a2(i)
         end
     end
     break
 end
mn=mn+1


                 [bs1,bs2]=find(bs==d(ls(1),:))                     %回退找到T接点的处理办法
                if(length(bs2)~=0)
                    ls=[bs2 ls]
                else
                    ls(1)=[]
                end
                [s1,s2]=find(ls==ls(1))
                if(length(s2)==2)
                ls2{mm}=ls                                            %找到第二类环
                ls(1)=[]
          



  [t1,t2]=find(d(ls(1),:)==1)
     [t3,t4]=size(t2)
     w=0
     for i=1:t4
          [bs1,bs2]=find(bs==t2(i))                     %回退找到T接点的处理办法
        if(length(bs2)==0)
             w=w+1
         end
     end



                              elseif((i2<3)|i2==3&d(ls(1),i)==0)
                                  ls=ls
                                  continue