牛顿-拉夫逊迭代法极坐标潮流计算C语言程序

/＊最大矩阵阶数＊//*迭代次数*//＊循环控制变量*//＊中间变量*//＊节点数*//＊支路数＊//*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...