k=0;
while
%解析来是程序1.m
for i=1:bus_num
PI=0.0;
QI=0.0;
for M=IFD(i):(IFD(i+1)-1)
J=JF(M);
VG=V(i)*V(J)*GF(M);
VB=V(i)*V(J)*BF(M);
HIJ=H(i)-H(J);
PI=PI+VG*(cos(HIJ))+VB*(sin(HIJ));
QI=QI+VG*(sin(HIJ))-VB*(cos(HIJ));
end
PI=PI+GD(i)*V(i)*V(i);
QI=QI-BD(i)*V(i)*V(i);
if bus_Type(i)==3
DP(i)=0;
DQ(i)=0;
G=1;
if i==NG(G)
PG(G)=PI;
QG(G)=QI;
else
G=G+1
if G>GM
break
end
end
else G=1;
if i==abs(NG(G))
if NG(G)>0
DP(i)=(PG(G)-PI)/V(i);
DQ(i)=(QG(G)-QI)/V(i);
else DP(i)=(PG(G)-PI)/V(i);
DQ(i)=0;