admittance_matrix.m

程序实现用牛顿－拉夫逊迭代法进行潮流计算。不仅有潮流计算

%*******************************
%Filename:admittance_matrix.m
%Author:Hweel_Zheng(郑奕辉)     
%First created:2008.08.23      
%Last mended:2008.08.25
%******************************
%节点导纳矩阵的形成
function Y=admittance_matrix(bus,line)
%Purpose: build admittance matrix Y from the line data
%
%Input: bus-bus data
%       line-line data
%Output:Y-admittance matrix
[nb,mb]=size(bus);
[nl,ml]=size(line);

Y=zeros(nb,nb);%对导纳矩阵赋初值0，一个nb阶零矩阵

%下面针对各条线路类型修改节点导纳矩阵
for k=1:nl
    I=line(k,1);    %读入线路参数
    J=line(k,2);
    Zt=line(k,3)+j*line(k,4);%Zt=R+jX
    if I~=0&J~=0
        Yt=1/Zt;
    end
    Ym=line(k,5)+j*line(k,6);
    K=line(k,7);%变比
    
    if(K==0)&(J~=0)     %变通线路：K＝0（即没有变压器）
                        %且J！＝0（即不是对地支路）
        Y(I,I)=Y(I,I)+Yt+Ym;
        Y(J,J)=Y(J,J)+Yt+Ym;
        Y(I,J)=Y(I,J)-Yt;
        Y(J,I)=Y(I,J);
    end
    
    if(K==0)&(J==0)     %对地支路
        Y(I,I)=Y(I,I)+Ym;
    end
    
    if K>0      %变压器线路：Zt和Ym为折算到i侧的值，K在j侧
        Y(I,I)=Y(I,I)+Yt+Ym;
        Y(J,J)=Y(J,J)+Yt/K/K;
        Y(I,J)=Y(I,J)-Yt/K;
        Y(J,I)=Y(I,J);
    end
    
    if K<0      %变压器线路：Zt和Ym为折算到i侧的值，K在i侧
        Y(I,I)=Y(I,I)+Yt+Ym;
        Y(J,J)=Y(J,J)+K*K*Yt;
        Y(I,J)=Y(I,J)+K*Yt;
        Y(J,I)=Y(I,J);
    end
end