%傳統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