ybus.m

来自「电力系统计算节点数不多的潮流程序,数度比较快」· M 代码 · 共 36 行

%节点导纳矩阵的形成 
function     Y=ybus(line,nb,nl)
Y=zeros(nb,nb);         % 对导纳矩阵赋初值0 
for k=1:nl,
    I=line(k,1);
    J=line(k,2);
    Zt=line(k,3)+j*line(k,4);
    Yt=1/Zt;
    Ym=line(k,5)+j*line(k,6);
    K=line(k,7);
    
    if (K==0)&(J~=0)                 % 普通线路: K=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)               % 对地支路: K=0,J=0,R=X=0;
        Y(I,I)=Y(I,I)+Ym;
    end
    if K>0                        % 变压器线路: Zt和Ym为折算到1侧的值,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为折算到k侧的值,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