formy.m

可以计算电力系统概率潮流

function Y=FormY(bus, line)
% form admittance matrix for AC load flow study
glo_var

global nline;

nbus=length(bus(:,1));
nline=length(line(:,1));

Y=zeros(nbus, nbus);
for k=1:nline
    frombus=line(k,1);
    tobus=line(k,2);
    frompos=FindPos(bus, frombus);
    topos=FindPos(bus, tobus);
    % nl is number of lines
    nl=line(k, 3);
    % line reactance 
    r=line(k, 5)/nl;
    x=line(k, 6)/nl;
    g=line(k, 7)*nl;
    
    Y(frompos, frompos)=Y(frompos, frompos)+1.0/(r+x*i)+g/2;
    Y(frompos, topos)=-1.0/(r+x*i);
    Y(topos, frompos)=-1.0/(r+x*i);
    Y(topos, topos)=Y(topos, topos)+1.0/(r+x*i)+g/2;
end

return