【正文】 =0) t=t+v[j]*(yg[i][j]*cos(w[i]w[j])+yb[i][j]*sin(w[i]w[j]))。 } pp[i]=pq[i][2]/v[i]t。 } 25 /*計(jì)算 m 個(gè)節(jié)點(diǎn)不平衡無(wú)功并除其電壓 */ for(i=1。i=m。i++) { t=0。 for(j=1。j=n。j++) { if (yb[i][j]!=0) t=t+v[j]*(yg[i][j]*sin(w[i]w[j])yb[i][j]*cos(w[i]w[j]))。 } qq[i]=pq[i][3]/v[i]t。 } for(i=1。i=n1。i++) { if (fabs(pp[i])) { goto a4 。 } } for(i=1。i=m。i++) { if (fabs(qq[i])) { goto a2。 } } goto a3。 /*利用 B1 因子表對(duì)右端項(xiàng)變換 */ a4:for(i=1。i=n2。i++) { t=0。 for(j=1。j=i。j++) { t=t+b1[i+1][j]*pp[j]。 } pp[i+1]=pp[i+1]+t。 } nl=n1。 /* 回代求解 */ pp[nl]=pp[nl]/b1[nl][nl]。 for(k=1。k=nl1。k++) { t=0。 for(j=nlk+1。j=nl。j++) { 26 t=t+b1[nlk][j]*pp[j]。 } pp[nlk]=(pp[nlk]t)/b1[nlk][nlk]。 } /*修正電壓角 */ for(i=1。i=nl。i++) { w[i]=w[i]pp[i]。 } /* 利用 b2 因子表對(duì)右端項(xiàng)變換 */ a2:for(i=1。i=m1。i++) { t=0。 for(j=1。j=i。j++) { t=t+b2[i+1][j]*qq[j]。 } qq[i+1]=qq[i+1]+t。 } qq[m]=qq[m]/b2[m][m]。 /*回代求解 */ for(k=1。k=m1。k++) { t=0。 for(j=nk+1。j=m。j++) { t=t+qq[j]*b2[mk][j]。 } qq[mk]=(qq[mk]t)/b2[mk][mk]。 } /*修正電壓模 */ for(i=1。i=m。i++) { v[i]=v[i]qq[i]。 } kk=kk++。 goto a1。 a3:for(i=1。i=n。i++) { e[i]=v[i]*cos(w[i])。 f[i]=v[i]*sin(w[i])。 ww[i]=w[i]*180/。 } /* 計(jì)算支路功率 */ for(k=1。k=ls。k++) 27 { i=z[k][1]。j=z[k][2]。kb=z[k][5]。 t1=yg[i][j]*(e[i]e[j])+yb[i][j]*(f[i]f[j])。 t2=yb[i][j]*(e[i]e[j])yg[i][j]*(f[i]f[j])。 p1=e[i]*t1+f[i]*t2。/*IJ 節(jié)點(diǎn)間有功的部分 Vi(UiUj)Yij*/ q1=f[i]*t1e[i]*t2。 p2=e[j]*(t1)+f[j]*(t2)。 q2=f[j]*(t1)e[j]*(t2)。 if (kb0) { gg[k][1]=p1+(kb1)*v[i]*v[i]*(yg[i][j])。 tt[k][1]=q1(kb1)*v[i]*v[i]*(yb[i][j])。 gg[k][2]=p2+(1kb)*v[j]*v[j]*(yg[i][j])/kb。 tt[k][2]=q2(1kb)*v[j]*v[j]*(yg[i][j])/kb。 } else { gg[k][1]=p1。gg[k][2]=p2。/*有功 */ tt[k][1]=q1+kb*v[i]*v[i]。/*無(wú)功 */ tt[k][2]=q2+kb*v[j]*v[j]。 } } /*計(jì)算支路功率損耗 */ for(k=1。k=ls。k++) { pw[k]=gg[k][1]+gg[k][2]。 qw[k]=tt[k][1]+tt[k][2]。 } tp=0。tq=0。 for(k=1。k=ls。k++) { tp=tp+pw[k]。 tq=tq+qw[k]。 } printf(迭代 %d 次 電壓實(shí)部 e 電壓虛部 f 電壓模 V 電壓角 (度 )\n,kk)。 for(i=1。i=n。i++) { printf(節(jié)點(diǎn) %6d%%%%\n,i,e[i],f[i],v[i],ww[i])。 } printf(\n)。 printf(支路功率 : 有功 無(wú)功 \n)。 for(k=1。k=ls。k++) { i=z[k][1]。j=z[k][2]。 printf (支路 %d%d%%\n,i,j,gg[k][1],tt[k][1])。 28 printf (支路 %d%d%%\n,j,i,gg[k][2],tt[k][2])。 } printf(支路損耗 : 有功 無(wú)功 \n)。 for(k=1。k=ls。k++) { i=z[k][1]。j=z[k][2]。 printf (支路 %d%d%%\n,i,j,pw[k],qw[k])。 } printf(總網(wǎng)耗 : 有功 無(wú)功 \n)。 printf(%%,tp,tq)。goto begin。 } 29 潮流結(jié)果 N=11 LS=11 M=9 支路數(shù)據(jù) 節(jié)點(diǎn) I 節(jié)點(diǎn) J 電阻 電抗 變比或電納 支路 1 1 3 支路 2 3 10 支路 3 10 5 支路 4 5 6 支路 5 1 2 支路 6 2 4 支路 7 4 9 支路 8 9 8 支路 9 8 7 支路 10 7 6 支路 11 6 11 PQ節(jié)點(diǎn) 注入有功 注入無(wú)功 1 2 3 4 5 6 7 8 9 PV節(jié)點(diǎn) 注入有功 電壓模 10 平衡節(jié)點(diǎn) 電壓模 電壓角 11 迭代 70次 電壓實(shí)部 e 電壓虛部 f 電壓模 V 電壓角 (度 ) 節(jié)點(diǎn) 1 節(jié)點(diǎn) 2 節(jié)點(diǎn) 3 節(jié)點(diǎn) 4 節(jié)點(diǎn) 5 節(jié)點(diǎn) 6 節(jié)點(diǎn) 7 節(jié)點(diǎn) 8 節(jié)點(diǎn) 9 節(jié)點(diǎn) 10 節(jié)點(diǎn) 11 30 支路功率 : 有功 無(wú)功 支路 13 支路 31 支路 310 支路 103 支路 105 支路 510 支路 56 支路 65 支路 12 支路 21 支路 24 支路 42 支路 49 支路 94 支路 98 支路 89 支路 87 支路 78 支路 76 支路 67 支路 611 支路 116 按數(shù)字鍵繼續(xù) :6 支路損耗 : 有功 無(wú)功 支路 13 支路 310 支路 105 支路 56 支路 12 支路 24 支路 49 支路 98 支路 87 支路 76 支路 611 總網(wǎng)耗