1、i.num,jdi.p,jdi.num,jdi.U,jdi.num,jdi.q); fprintf(fp2, 平衡节点: 节点%d e%d=%f f%d=%fn, jdn.num,jdn.num,jdn.U,jdn.num,jdn.zkj); fprintf(fp2, -n); for(i=1;i=m;i+) fprintf(fp2, 支路%d: 相关节点:%d,%d 非标准变比:kx=%f R=%f X=%f n, i,zli.p1,zli.p2,zli.kx,zli.r,zli.x); fprintf(fp2, =n); /*-复数运算函数-*/ double mozhi(double a
2、0,double b0) /*复数求模值函数*/ mo=sqrt(a0*a0+b0*b0); return mo; double ji(double a1,double b1,double a2,double b2) /*复数求积函数 a1为电压模值，a2为阻抗角，a3为导纳实部，a4为导纳虚部*/ a1=a1*cos(b1); b1=a1*tan(b1); c1=a1*a2-b1*b2; d1=a1*b2+a2*b1; return c1; return d1; double shang(double a3,double b3,double a4,double b4) /*复数除法求商函数*
3、/ c2=(a3*a4+b3*b4)/(a4*a4+b4*b4); d2=(a4*b3-a3*b4)/(a4*a4+b4*b4); return c2; return d2; /*-计算节点导纳矩阵-*/ void Form_Y() for(i=1;i=n;i+) /*节点导纳矩阵元素附初值*/ for(j=1;j=n;j+) Gij=Bij=0; for(i=1;i=n;i+) /*节点导纳矩阵的主对角线上的元素,非对地导纳加入相应的主对角线元素中(考虑非标准变比）*/ for(j=1;j=m;j+) if(zlj.p1=i) if(zlj.kx=1) mozhi(zlj.r,zlj.x);
4、 if(mo=0) continue; shang(1,0,zlj.r,zlj.x); Gii+=c2; Bii+=d2; else mozhi(zlj.r,zlj.x); if(mo=0) continue; shang(1,0,zlj.r,zlj.x); Gii+=c2/zlj.kx+c2*(1-zlj.kx)/(zlj.kx*zlj.kx); Bii+=d2/zlj.kx+d2*(1-zlj.kx)/(zlj.kx*zlj.kx); else if(zlj.p2=i) if(zlj.kx=1) mozhi(zlj.r,zlj.x); if(mo=0) continue; shang(1,
5、0,zlj.r,zlj.x); Gii+=c2; Bii+=d2; else mozhi(zlj.r,zlj.x); if(mo=0) continue; shang(1,0,zlj.r,zlj.x); Gii+=c2/zlj.kx+c2*(zlj.kx-1)/zlj.kx; Bii+=d2/zlj.kx+d2*(zlj.kx-1)/zlj.kx; for(k=1;k=m;k+) /*节点导纳矩阵非主对角线上（考虑非标准变比）的元素*/ if(zlk.kx=1) i=zlk.p1; j=zlk.p2; mozhi(zlk.r,zlk.x); if(mo=0) continue; shang(1
6、,0,zlk.r,zlk.x); Gij-=c2; Bij-=d2; Gji=Gij; Bji=Bij; else i=zlk.p1; j=zlk.p2; mozhi(zlk.r,zlk.x); if(mo=0) continue; shang(1,0,zlk.r,zlk.x); Gij-=c2/zlk.kx; Bij-=d2/zlk.kx; Gji=Gij; Bji=Bij; /*-输出节点导纳矩阵-*/ fprintf(fp2,nn * 潮流计算过程 *n); fprintf(fp2,n =n); fprintf(fp2,n 节点导纳矩阵为:); for(i=1;i=n;i+) fprin
7、tf(fp2,n ); for(k=1;k=0) fprintf(fp2,+j); fprintf0挟釕0釕0釕0釕0倀鉣0捑刀鉣0捓吀鉣0捕嘀鉣0捗砀鉑0兹稀鉑0养簀鉑0兽縀鉑0兿怀鉣0捡戀鉣0捣搀鉣0捥昀鉣0捧瀀鉣0捱爀鉣0捳琀鉣0捵瘀鉣0捷退鉣0掑鈀鉣0掓鐀鉣0掕阀鉣0掗鉣0掱鉣0掳鉣0掵鉣0掷鉣0揑鉣0揓鉣0揕鉣0揗鉣0握分鉣0揣堆鉣0揥将鉣0揧弆鉤0搙戆鉤0搛朆鉤0搝欆鉤0搟渆鉤0摂理鉤0摄眆鉤0摆簆鉤0揸舆鉣0揺蔆鉣0揼訆鉣0揾贆鉣0鄆釕0阆釕0餆釕0鰆釕0摘夀鉤0摚嬀鉤0摜崀鉤0摞开鉤0摰焀鉤0摲猀鉤0摴甀鉤0摶眀鉤0撈褀鉤0撊謀鉤0撌贀鉤0撎輀鉤0撘餀鉤0撚鬀鉤0撜鴀鉤0撞鼀鉤0撨鉤0撪鉤
