电力系统潮流计算高斯-赛德尔迭代法

clear
FilePath=input（‘Please input the dirctory of network raw data:‘‘s‘）;
Network=xlsread（FilePath1）;
NodeNum=max（max（Network（:2））max（Network（:3）））;
data=xlsread（FilePath4）;


VS=zeros（NodeNum3）;

for N=1:NodeNum
    VS（N1）=N;
    VS（N2）=data（N2）*cos（data（N3））+i*data（N2）*sin（data（N3））;
    VS（N3）=（data（N4）-data（N6）+i*（data（N5）-data（N7）））;
end

Y=zeros（NodeNum）;
linkNum=length（Network（:1））;
for L=1:linkNum
    y1=1/（Network（L4）+i*Network（L5））;
    j=Network（L2）;
    k=Network（L3）;
    Y（jj）=Y（jj）+y1+i*Network（L6）;
    Y（jk）=Y（jk）-y1;
    Y（kj）=Y（kj）-y1;
    Y（kk）=Y（kk）+y1+i*Network（16）;
end

Transformer=xlsread（FilePath2）;
TranNum=length（Transformer（:1））;
for T=1:TranNum
    j=Transformer（T2）;
    k=Transformer（T3）;
    for J=1:linkNum
        if （j==Network（J2））
          if （k==Network（J3））
              yl=1/（Network（L4）+i*Network