文档介绍：#include<>#include<>#defineN5//节点数#definen_PQ3//PQ节点数#definen_PV1//PV节点数#[N][N]={{,-,-,0,0},{-,,-,0,0},{-,-,,0,0},{0,0,0,0,0},{0,0,0,0,0}};//导纳矩阵的实部floatB[N,N]={{-,1,,0,0},{1,-,,,0},{,,-,0,},{0,-,0,,0},{0,0,,0,-}};//导纳矩阵的虚部floatU_amp[N]={,,,,,};//电压初值floatU_ang[N]=[0,0,0,0,0]//电压相角初值floatPs[N-1]={,2,,5};//有功初值floatQs[N-2]={,1,};//无功初值floatdU_amp[n_PQ],dU_ang[N-1];//电压幅值和相角修正值floatP1[N],Q1[N],dP[N],dQ[n_PQ],H[N-1,N-1],L[n_PQ-1,n_PQ-1];floatJ[N+n_PQ-1,N+n_PQ-1],M[N+n_PQ-1,1],N[N+n_PQ-1,N+n_PQ-1];//J[N+n_PQ-1-1,N+n_PQ-1]雅可比矩阵,N[N+n_PQ-1-1,N+n_PQ-1]雅可比矩阵逆矩阵floatdS[N+n_PQ-1,1],T,U[N],SS;floatprecision=;//precision计算精度intkk=0;//迭代次数inti,j,k;//主程序===================================================================while(precision>){for(i=1;i<N;i++){U[i]=U_amp[i]*cos(U_ang[i])+U_amp[i]*sin(U_amp[i]*);}for(i=0;i<N;i++){P1[i]=0;Q1[i]=0;for(j=0;j<N;j++){P1[i]=P1[i]+U_amp[i]*(G[i][j]*U_amp[j]*cos(U_ang[j])+B[i][j]*U_amp[j]*sin(U_ang[j]));Q1[i]=Q1[i]+U_amp[i]*(G[i][j]*U_amp[j]*sin(U_ang[j])-B[i][j]*U_amp[j]*cos(U_ang[j]));}dP[i]=Ps[i]-P1[i];if(i<n_PQ)dQ[i]=Qs[i]-Q1[i];}for(i=1;i<N+n_PQ;i++){M(i,1)=dP[i];If(i>N-1)M[i,1]=dQ[i];}for(i=1;i<N,i++){for(j<1;j<N;j++)H[i,j]=0;If(i<N-1&&j<N-1){L[i,j]=0;}}for(i=0;