三相不平衡配电网的潮流计算

隐式高斯法求解三相不平衡配电网的潮流计算。

%%本程序为求解节点导纳矩阵

originalData36
t1=menu（‘选择系统是否含有变压器；是选1，否选2；‘‘1‘‘2‘）
 switch t1
     case 1
         disp（‘系统含有变压器‘）
         linedata=oriLineData（）;
         transformdata=transformdata（）;
         I = linedata（:2）;      % 起节点
         J = linedata（:3）;      % 终节点
         l = linedata（:4）;      %线路长度
         nbus = max（max（I）max（J））;    % 节点数
         nbranch = length（J）;          % 支路数
         Y=zeros（nbus*3）;
         nline=size（linedata1）;
         ntrans=size（transformdata1）;
        if nline>=1
         for m=1:nbranch
            I = linedata（m2）;      % 起节点
            J = linedata（m3）;      % 终节点
            l = linedata（m4）;      %线路长度为英尺feet
            L=l/5280;               %线路长度转换成mile
            Zl=oriZl*L*SB/UB^2/1000;
            Yl=inv（Zl）;
           Y（3*（I-1）+1:3*（I-1）+33*（J-1）+1:3*（J-1）+3）=Y（3*（I-1）+1:3*（I-1）+33*（J-1）+1:3*（J-1）+3）-Yl;
           Y（3*（J-1）+1:3*（J-1）+33*（I-1）+1:3*（I-1）+3）=Y（3*（I-1）+1:3*（I-1）+33*（J-1）+1:3*（J-1）+3）;
           Y（3*（I-1）+1:3*（I-1）+33*（I-1）+1:3*（I-1）+3）=Y（3*（I-1）+1:3*（I-1）+33*（I-1）+1:3*（I-1）+3）+Yl;
           Y（3*（J-1）+1:3*（J-1）+33*（J-1）+1:3*（J-1）+3）=Y（3*（J-1）+1:3*（J-1）+33*（J-1）+1:3*（J-1）+3）+Yl;  
         end
        
    if ntrans>=1
        typeT=transformdata（:4）;
       for m=1:ntrans 
         I = transformdata（m2）;      % 起节点
         J = transformdata（m3）;      % 终节点
         U1=transformdata（m6）;
         U2=transformdata（m7）;
         k=（U1*UB2）/（U2*UB）;  
         SN=transformdata（m5）;
         R_100=transformdata（m8）;
         X_100=transformdata（m9）;
         ZT=U1^2/SN*1000;
         RT=R_100/100*ZT*SB/UB^2/1000;
         XT=X_100/100*ZT*SB/UB^2/1000;  
         yt=1.0/（RT+j*XT）;
          if typeT（m）==1  %类型-1表示变压器为YNyn接线方式
             y11=[yt/k^2 0 0 ;0 yt/k^2 0;0 0 yt/k^2];
             y12=[yt/k 0 0;0 yt/k 0;0 0 yt/k];
             y22=[yt 0 0;0 yt 0;0 0 yt];
          end
          if typeT（m）==2 %类型-2表示变压器为YNd11接线方式
              y11=[yt/k^2 0 0 ;0 yt/k^2 0;0 0 yt/k^2];
              y12=[yt/k/sqrt（3） 0 0;0 yt/k/sqrt（3） 0;0 0 yt/k/sqrt（3）];
              y22=[yt/3 0 0;0 yt/3 0;0 0 yt/3];
          end
              Y（3*（I-1）+1:3*（I-1）+33*（J-1）+1:3*（J-1）+3）=Y（3*（I-1）+1:3*（I-1）+33*（J-1）+1:3*（J-1）+3）-y12;
              Y（3*（J-1）+1:3*（J-1）+33*（I-1）+1:3*（I-1）+3）=Y（3*（I-1）+1:3*（I-1）+33*（J-1）+1:3*（J-1）+3）;
              Y（3*（I-1）+1:3*（I-1）+33*（I-1）+1:3*（I-1）+3）=Y（3*（I-1）+1:3*（I-1）+33*（I-1）+1:3*（I-1）+3）+y11;
              Y（3*（J-1）+1:3*（J-1）+33*（J-1）+1:3*（J-1）+3）=Y（3*（J-1）+1:3*（J-1）+33*（J-1）+1:3*（J-1）+3）+y22;         
       end
    end
 end
 
     case 2
         disp（‘系统没有变压器‘）
         linedata=oriLineData（）;
         I = linedata（:2）;      % 起节点
         J = linedata（:3）;      % 终节点
         l = linedata（:4）;      %线路长度
         nbus = max（max（I）max（J））; % 节点数
         nbranch = length（J）;  % 支路数
         Y=zeros（nbus*3）;
    for m=1:nbranch
           I = linedata（

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       4310  2013-05-08 12:54  yinshigaosi_method\jiediandaona.asv

     文件       3792  2013-05-15 21:39  yinshigaosi_method\jiediandaona.m

     文件       2911  2013-05-08 14:41  yinshigaosi_method\originalData12.asv

     文件       2716  2013-05-08 14:44  yinshigaosi_method\originalData12.m

     文件       6802  2013-05-08 15:19  yinshigaosi_method\originalData36.m

     文件       1972  2013-05-08 14:31  yinshigaosi_method\originalData6.asv

     文件       1955  2013-05-08 14:36  yinshigaosi_method\originalData6.m

     文件       3809  2013-05-15 15:20  yinshigaosi_method\yinshigaosifa.asv

     文件       3809  2013-05-15 21:44  yinshigaosi_method\yinshigaosifa.m

     目录          0  2013-05-16 13:33  yinshigaosi_method

----------- ---------  ---------- -----  ----

                32076                    10