/*最大矩阵阶数*//*迭代次数*//*循环控制变量*//*中间变量*//*节点数*//*支路数*//*PQ 节点数*//*PV 节点数*//*迭代精度*/tt;/*中间变量*//*复数运算函数的返回值*//*节点导纳矩阵中/*雅克比矩阵、不平/*利用牛顿—拉夫逊迭代法(极坐标形式),计算复杂电力系统潮流,具有收敛性好,收敛速度快等优点.所有参数应归算至标幺值下./*可计算最大节点数为 100,可计算 PQ,PV/*可计算非标准变比和平行支路*/#includevstdio。h〉#include〈math。h〉#include〈stdlib.h〉#defineM100#defineNl100inti,j,k,a,b,c;intt,l;doubleP,Q,H,J;intn,m,pq,pv;doubleeps;doubleaa[M],bb[M],cc[M],dd[M],max,rrdoublemo,c1,d1,c2,d2;doubleG[M][M],B[M][M],Y[M][M];的实部、虚部及其模方值*/doubleykb[M][M],D[M],d[M],dU[M];衡量矩阵*/inth,number;structjd{intnum,ty;型*/doublep,q,S,U,zkj,dp,dq,du,dj;模值,电压模值,阻抗角平衡量、电压修正量*/}jd[M];structzl{intnumb;intp1,p2;doublekx;doubler,x;}zl[M];FILE*fp1,*fp2;voiddata(){/*节点结构体*//*num 为节点号,ty 为节点类/*节点有功、无功功率,功率牛顿——拉夫逊中功率不/*支路结构体*//*numb 为支路号*//*支路的两个节点*//*非标准变比*//*支路的电阻与电抗*//*读取数据*/fpl=fopen("input。txt","r");fscanf(fp1,”%d,%d,%d,%d,%lf\n",&n,&m,&pq,&pv,&eps);/*输入节点数,支路数,PQ 节点数,PV 节点数和迭代精度*/for(i=l;i<=n;i++)/*输入节点编号、类型、输入功率和电压初值*/{fscanf(fpl,”%d,%d”,&number,&h);if(h==1)/*类型 h=1是 PQ 节点*/{fscanf(fp1,"%lf,%lf,%lf,%lf\n",&jd[i].p,&jd[i].q,&jd[i]。U,&jd[i].zkj);jd[i]。num=number;jd[i]。ty=h;}if(h==2)/*类型 h=2是 pv 节点*/{fscanf(fp1,”,%lf,%lf,%lf\n”,&jd[i].p,&jd[i]。U,&jd[i].zkj);jd[i].num=number;jd[i].ty=h;jd[i].q=-1。567;}if(h==3)/*类型 h=3是平衡节点*/{fscanf(fp1,",%lf,%lf\n”,&jd[i]。U,&jd[i]。zkj);jd[i].num=number;jd[i]。ty=h;}}for(i=1;i〈=m;i++)/*输入支路阻抗*/fscanf(fp1,”%d,%lf,%d,%d,%lf,%lf\n”,&zl[i]。numb,&zl[i]。kx,&zl[i].p1,&zl[i].p2,&zl[i].r,&zl[i]。x);fclose(fp1);if((fp2=fopen(”output.txt",”w"))==NULL){print...