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