电网参数:
1 1 2 0.10000 0.41000 0.01528 1.00000
2 1 3 0.00000 0.30000 0.00000 1.10000
3 1 4 0.12000 0.50000 0.01920 1.00000
4 2 4 0.08000 0.40000 0.01413 1.00000
节点导纳矩阵形成过程
for(l=0; l<nL; l++){
i=sL[l].NumI;
j=sL[l].NumJ;
r=sL[l].R;
x=sL[l].X;
d1=r*r+x*x;
g=r/d1;
b=-x/d1;
if(fabs(sL[l].K-1.0)<0.000001)
{/*Normal lines or transformers*/
YG[i][i]=YG[i][i]+g;
YG[j][j]=YG[j][j]+g;
YB[i][i]=YB[i][i]+b-sL[l].B;
YB[j][j]=YB[j][j]+b-sL[l].B;
YG[i][j]=YG[i][j]-g;
YG[j][i]=YG[j][i]-g;
YB[i][j]=YB[i][j]-b;
YB[j][i]=YB[j][i]-b;
}
else
{/*abnormal transformer ratio*/
YG[i][i]=YG[i][i]+g/(sL[l].K)/(sL[l].K);
YG[j][j]=YG[j][j]