【正文】 amp。zx[ip])。swap(amp。zyk[i], amp。zyk[ip])。//**** 節(jié)點優(yōu)化 ****//int ii1, ii2, zls2, nomax。int i, j, l, k1, k。int temp。zls2 = zls + zls。for(i = 1。i old[i] = nnew[i] = 0。//先清零。由此可知：NSNS必須大于2*zls。for(i = 1。i {old[i] = iabs(izl[i])。old[i + zls] = iabs(jzl[i])。} //變壓器節(jié)點號由正變負，old[]前zls個為左節(jié)點號，后zls個為右節(jié)點號。for(i = 1。i {k1 = i + 1。for(j = k1。j if(old[i] old[j])iswap(amp。old[i], amp。old[j])。//交換整數(shù)old[i]、old[j]。小節(jié)點號排在支路左側(cè)。}nomax = old[zls2]。//nomax 即是最大節(jié)點號。 l = 1。for(i = 1。i {ii1 = old[l]。old[i] = ii1。for(j = l。j {ii2 = old[j]。if(ii1!= ii2){l = j。break。}nnew[i] ++。} } for(i = 1。i //即為第 i 行的非零元素個數(shù)。for(k = i1。k if(kq0!= 0){ai =yb[k]。bi = yg[k]。}//yb[]存放導(dǎo)納陣非對角元的虛部(yg[] + j*yb[])else{ai = yg[k]。bi = yb[k]。}//yg[]存放導(dǎo)納陣非對角元的實部(yg[] + j*yb[])j = iy[k]。//iy[]存放的是列足碼。vj = vi * v[j]。// v[]存放節(jié)點電壓幅值。ci = va[i]bi。//非對角元素 j 的功率} }void bbhl(int kq0){//**** 本函數(shù)計算各節(jié)點的功率誤差，求最大功率誤差dsm ****// //**** 和常數(shù)項b[i]。(程序框圖見P423)****//int i。double vi, vj。double pl0, pg0。dsm = 。// dsm 即為最大功率誤差。for(i = 1。i vi = v[i]。// v[]存放節(jié)點電壓幅值。vj = v0[i]。// v0[]存放節(jié)點初始電壓幅值。v0[]存放的是最后一個節(jié)點數(shù)據(jù)。// vi[] // nopt[] =1: pv節(jié)點，對于pv節(jié)點，即為該點應(yīng)維持的電壓值。// nopt[] = 0: pq節(jié)點，對于pq節(jié)點，即為系統(tǒng)平均電壓值。// vj 此時表示的是節(jié)點正常電壓的會定值。if(kq0 == 0){ pl0 = pl[i]。pg0 = pg[i]。}// 負荷p,發(fā)電機pelse{pl0 = ql[i]。pg0 = qg[i]。}// 負荷q,發(fā)電機q if(kk2 == 0)pl0 = pl0 *((vipl0。// pv節(jié)點 if(i == mpj amp。amp。 kq0 == 0)pg[i] = w[i]w[i]。//pv節(jié)點(nobt[] =1)和平衡節(jié)點(mpj)不參與 //求最大功率誤差和常數(shù)項的運算if(((kq0!= 1)||(nobt[i]!=1))amp。amp。(i!= mpj)){ if(fabs(b[i]) fabs(dsm)){dsm = b[i]。dsd = i。}// dsm 即為最大功率誤差,dsd存放其對應(yīng)的節(jié)點號。b[i] = b[i] / vi。// 計算修正方程式的常數(shù)項。} elseb[i]=。// pv節(jié)點(nobt[] =1)和平衡節(jié)點(mpj)不參與// 求最大功率誤差和常數(shù)項的運算。} }node_output()//{ //**** 輸出節(jié)點數(shù)據(jù)和最小電壓幅值、相角(角度)及其節(jié)點號 ****// //****(程序框圖見p426 F116)****//double vi, ci。int i, j, oldnumber。printc(39。+39。, 72)。fprintf(fp2, “n%5s%8s%10s%11s%11s%11s%11sn”, “I”，“V”, “CA”, “PL”, “QL”, “PG”, “QG”)。vmin = v[1]。dsd = 1。for(i = 1。i j = nnew[i]。oldnumber = newsort[i]。//轉(zhuǎn)化為相應(yīng)舊號ci = va[j] * 。//弧度轉(zhuǎn)化為角度。vi = v[j]。if(vi vmin){vmin = vi。dsd = j。}// vmin即為最小電壓,dsd存放其對應(yīng)的新節(jié)點號。fprintf(fp2, “n%5d%%”, oldnumber, vi, ci)。fprintf(fp2, “%%%%”, pl[j], ql[j]，pg[j], qg[j])。} printc(39。39。, 72)。} void branch_output(){//**** 本函數(shù)輸出支路數(shù)據(jù)。(程序框圖見p428 F117)int ii, jj。double r, x, yk, zf, vi, vj, ci, cj。int i, j, l。double de, df, ei, ej, fi, fj, fii, fir, pij, pji, qij, qji。****// dph = 。//統(tǒng)計系統(tǒng)有功網(wǎng)損。dqh = 。//統(tǒng)計系統(tǒng)無功網(wǎng)損。fprintf(fp2, “n%5s%5s%10s%12s%12s%12sn ”, “I”, “J”, “PIJ”，“QIJ”, “PJI”, “QJI”)。for(i =1。i {j = idk[i]。dkk[i] = v[j] * v[j] / dkk[i]。} for(l = 1。l {ii = iabs(izl[l])。// izl[]: 支路左節(jié)點號。jj = iabs(jzl[l])。// jzl[]: 支路右節(jié)點號。i=nnew[ii]。j=nnew[jj]。//轉(zhuǎn)換為新節(jié)點號。ii = newsort[ii]。jj = newsort[jj]。//轉(zhuǎn)化為相應(yīng)舊號 r = zr[l]。x = zx[l]。yk = zyk[l]。vi = v[i]。//v[]: 電壓幅值。ci = va[i]。//va[]: 電壓相角。vj = v[j]。cj = va[j]。//支路左、右節(jié)點電壓值由極坐標轉(zhuǎn)換為直角坐標 ei = vi * cos(ci)。fi = vi * sin(ci)。// ei: 支路左節(jié)點電壓實部，fi: 支路左節(jié)點電壓虛部。ej = vj * cos(cj)。fj = vj * sin(cj)。// ej: 支路右節(jié)點電壓實部，fj: 支路右節(jié)點電壓虛部。if((izl[l] 0)||(jzl[l] 0))// 變壓器支路。{ if(izl[l] 0){ei = ei / yk。fi = fi / yk。}// yk=zyk[l]else{ej = ej / yk。fj = fj / yk。}yk = 。}de = eifj。zf = r * r + x * x。fii =(de * r + df * x)/ zf。fir =(df * rfir * ei。pji =fii * ejvi * vi * yk。qji = qji1。i ++){ip = i。k1 = iabs(izl[i])。k3 = iabs(jzl[i])。for(j = i + 1。j {k2 = iabs(izl[j])。k4 = iabs(jzl[j])。if(k2 k1 ||(k2 == k1 amp。amp。 k4 k3)){ip = j。k1 = k2。k3 = k4。}}if(i!= ip){iswap(amp。izl[i], amp。izl[ip])。iswap(amp。jzl[i], amp。jzl[ip])。swap(amp。zr[i], amp。zr[ip])。swap(amp。zx[i], amp。zx[ip])。swap(amp。zyk[i], amp。zyk[ip])。}} for(i = 1。i k = nob[i]。nob[i] = newsort[k]。printf(“nob[%d] = %dn”, i, nob[i])。} for(i = 1。i for(j = i+1。j {if(nob[i] nob[j]){k = nob[i]。nob[i] = nob[j]。nob[j] = k。}} } printf(“nob[%d] = %dn”, i, nob[i])。}void yy1(){//**** 本函數(shù)形成節(jié)點導(dǎo)納陣(不包括接地支路)****// int j1。double r, x, yk, zf, gij, bij。int i, j, i1, l, ll。for(i = 1。i {gii[i] = 。bii[i] = 。}// 導(dǎo)納陣對角元(與節(jié)點一一對應(yīng))先清零。for(i = 1。i yg[i] = 。yb[i] = 。} //導(dǎo)納陣非零非對角元(與支路一一對應(yīng))先清零。ll = 1。for(l = 1。l i1 = izl[l]。j1 = jzl[l]。i = iabs(i1)。j = iabs(j1)。if(i == j)// 支路左節(jié)點號。// 支路右節(jié)點號。// 變壓器支路有一節(jié)點號為負值。continue。//排除左、右節(jié)點號相等的情況。r = zr[l]。x = zx[l]。yk = zyk[l]。// zr[],zx[],zyk[]:支路三參數(shù)。zf = r * r + x * x。gij = r / zf。//bij =x / zf。bij =1/x。yg[ll] = yg[ll]bij。gii[i] = gii[i] + gij。bii[i] = bii[i] + bij。gii[j] = gii[j] + gij。bii[j] = bii[j] + bij。if((i!= iabs(izl[l + 1]))||(j!= iabs(jzl[l + 1])))ll++。} // 打印導(dǎo)納矩陣。對角元實部為gii，虛部為bii，// 非零非對角元實部為yb[]，虛部為yb[]，列足碼為iy[]。fprintf(fp2, “*******GII(1),BII(1)********n”)。printf2(gii,bii,n)。}void y3(){//**** 本函數(shù)形成節(jié)點導(dǎo)納陣，追加接地支路 ****//int j1。double r, x, yk, zf, gij, bij。int i, j, i1, l, ll, kk = 0。for(i = 1。i {j = idk[i]。bii[j] = bii[j]gij / yk。yb[ll] = yb[ll]bij / yk。}if((i!= iabs(izl[l + 1]))||(j!= iabs(jzl[l + 1]))){ll++。kk = 0。}elsekk = 1。} // 打印導(dǎo)納矩陣。對角元實部為gii，虛部為bii，// 非零非對角元實部為yb[]，虛部為yb[]，列足碼為iy[]。fprintf(fp2, “*******GII,BII********”)。printf2(gii,bii,n)。fprintf(fp2, “n*******YYYYY********”)。for(i = 1。i {if(i % 2 == 1)fprintf(fp2, “n”)。fprintf(fp2, “%%%8d”, yg[i], yb[i], iy[i])。} }第四篇：電力系統(tǒng)潮流計算程序include include include “” include “” include “”void main(){ coutPdQd。coutUd。double data[300][8]。//分別為0前面有功，1前面無功，2壓降，3中間有功，4中間無功,5變比,6后面有功，7后面無功int ii。ii=0。data[ii][0]=。data[ii][1]=。data[ii][2]=Ud。data[ii][3]=Pd。data[ii][4]=Qd。data[ii][5]=。data[ii][6]=。data[ii][7]=。ii++。coutwhile(1){double a。double b。double c。double d。double e。double f。int x。coutcinx。if(x!=0amp。amp。x!=1amp。amp。x!=2amp。amp。x!=3){coutcoutcinx。}if(x==0)break。if(x==1){coutcinab。Node N[100]。int i=0。N[i].Pn=a。N[i].Qn=b。data[ii][0]=。data[ii][1]=。data[ii][2]=。data[ii][3]=N[i].Pn。data[ii][4]=N[i].Qn。data[ii][5]=。data[ii][6]=。data