Z-bus 法计算潮流。

根据叠加原理，指定母线的电压可以通过根节点在该母线上产生的电压和该母线上的等值注入电流所产生的电压降叠加求得。

%传统Z-bus法潮流计算程序%
brnum=size（branch1）;
UI=zeros（brnum1）;
Uold=UI;
U=UI;
Sl=U;
Il=Sl;
SLi=Sl;
I=Sl;
busm=bus;
busold=bus;
busu=bus;
Y=zeros（brnum）;
for i=2:brnum  %形成导纳矩阵；
    %branch第一列支路号、第二列首节点号、第五列支路阻抗
    I=branch（i1）-1;
    J=branch（i2）-1;
    Y（II）=Y（II）+1/branch（i5）;
    Y（IJ）=-1/branch（i5）;
    Y（JI）=Y（IJ）;
    Y（JJ）=1/branch（i5）;
end
Y（11）=Y（11）+1/branch（15）;
Uso=bus（14）*（cos（bus（15））+1j*sin（bus（15）））;
%bus矩阵第二列节点有功、第三列无功、第四列电压幅值、第五列电压相角%
for i=1:brnum  %计算当松弛节点独立作用于整个配电网的情况下各母线的电压
    SLi（i）=bus（i+12）+1j*bus（i+13）;
    Uold（i）=bus（i+14）*（cos（bus（i+15））+1j*sin（bus（i+15）））;
end
k=0;
while 1  %计算网络中各条母线的等效注入电流、母线电压
    k=k+1;
    I=-conj（SLi）./conj（Uold）;
    UI=Y\I;
    Unew=UI+Uso;
    dU=abs（Unew-Uold）;
    dUmax=max（dU）;
    if（dUmax>0.00001）
        Uold=Unew;
    else
        break;
    end
    U=[Uso;Unew];
    %计算支路功率损耗
    for i=1:brnum
        Il（i）=（U（branch（i1））-U（branch（i2）））/branch（i5）;
        Sl（i）=（U（branch（i1））-U（branch（i2）））*conj（Il（i））;
    end
    %根据叠加原理得到节点电压
    for i=1:brnum
        busm（i+14）=abs（Unew（i））;
        busm（i+15）=angle（Unew（i））;
    end
    dV=abs（busm（:4）-busold（:4））;
    dtheta=abs（busm（:5）-busold（:5））;
    dVmax=max（dV）;
    dthetamax=max（dtheta）;
    dU=max（dVmaxdthetamax）;
    %收敛条件检验
    if（dU<=0.00001）
        break;
    else
        busold=busm;
        busu（:4:5）=busm（:4:5）;
    end
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1606  2018-05-13 10:43  Z_BUS.m