c# 节点导纳矩阵形成

电网参数：

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]