y.m

来自「利用该MATLAB程序可实现对一般电力系统的节点、支路及其损耗的潮流计算。」· M 代码 · 共 46 行

function Y=y(bus,line)
%目的：根据支路参数建立节点导纳矩阵
%输入：节点参数矩阵--bus；支路参数矩阵--line
%输出：节点导纳矩阵--Y
[nb,mb]=size(bus);
[nl,ml]=size(line);

Y=zeros(nb,nb);
for k=1:nl
    I=line(k,1);
    J=line(k,2);
    Zt=line(k,3)+j*line(k,4);
    if Zt==0
        disp('对地支路');
        Yt=inf;
    else
        Yt=1/Zt;
    end
    Ym=line(k,5)+j*line(k,6);
    K=line(k,7);
    
    if(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)   %对地支路K=J=0，R=X=0
        Y(I,I)=Y(I,I)+Ym;
    end
    if K>0  %K>0时变压器支路
        Y(I,I)=Y(I,I)+Yt+Ym;
        Y(J,J)=Y(J,J)+Yt/K^2;
        Y(I,J)=Y(I,J)-Yt/K;
        Y(J,I)=Y(I,J);
    end
    if K<0   %K<0时变压器支路
        Y(I,I)=Y(I,I)+Yt+Ym;
        Y(J,J)=Y(J,J)+Yt*K^2;
        Y(I,J)=Y(I,J)+Yt*K;
        Y(J,I)=Y(I,J);
    end
end