【正文】 / v0[i]。im jj1 = nfd[im]。di[i] = / b[i]。for(j = i1。x++。}nfd[i+1] = nfd[i] + fd[i]。fprintf(fp2, “%%8i”, u[i], iu[i])。for(i = 1。}void calc(int* iu, double* u, double* di, int* nfd, double* b){//**** 本函數(shù)利用因子表解線形方程組。i {bi = b[i]。j i2。j = i1。k1 = iabs(izl[i])。k4 = iabs(jzl[j])。k1 = k2。izl[ip])。swap(amp。zx[i], amp。zyk[ip])。zls2 = zls + zls。由此可知：NSNS必須大于2*zls。} //變壓器節(jié)點號由正變負(fù)，old[]前zls個為左節(jié)點號，后zls個為右節(jié)點號。j if(old[i] old[j])iswap(amp。小節(jié)點號排在支路左側(cè)。for(i = 1。j {ii2 = old[j]。} } for(i = 1。bi = yg[k]。//iy[]存放的是列足碼。//非對角元素 j 的功率} }void bbhl(int kq0){//**** 本函數(shù)計算各節(jié)點的功率誤差，求最大功率誤差dsm ****// //**** 和常數(shù)項b[i]。dsm = 。// v[]存放節(jié)點電壓幅值。// vi[] // nopt[] =1: pv節(jié)點，對于pv節(jié)點，即為該點應(yīng)維持的電壓值。pg0 = pg[i]。// pv節(jié)點 if(i == mpj amp。amp。b[i] = b[i] / vi。} }node_output()//{ //**** 輸出節(jié)點數(shù)據(jù)和最小電壓幅值、相角(角度)及其節(jié)點號 ****// //****(程序框圖見p426 F116)****//double vi, ci。, 72)。for(i = 1。//弧度轉(zhuǎn)化為角度。}// vmin即為最小電壓,dsd存放其對應(yīng)的新節(jié)點號。39。double r, x, yk, zf, vi, vj, ci, cj。//統(tǒng)計系統(tǒng)有功網(wǎng)損。for(i =1。l {ii = iabs(izl[l])。i=nnew[ii]。jj = newsort[jj]。vi = v[i]。vj = v[j]。// ei: 支路左節(jié)點電壓實部，fi: 支路左節(jié)點電壓虛部。if((izl[l] 0)||(jzl[l] 0))// 變壓器支路。fj = fj / yk。fii =(de * r + df * x)/ zf。i ++){ip = i。j {k2 = iabs(izl[j])。 k4 k3)){ip = j。izl[i], amp。jzl[ip])。swap(amp。zyk[i], amp。nob[i] = newsort[k]。j {if(nob[i] nob[j]){k = nob[i]。}void yy1(){//**** 本函數(shù)形成節(jié)點導(dǎo)納陣(不包括接地支路)****// int j1。i {gii[i] = 。i yg[i] = 。for(l = 1。j = iabs(j1)。continue。yk = zyk[l]。//bij =x / zf。bii[i] = bii[i] + bij。} // 打印導(dǎo)納矩陣。}void y3(){//**** 本函數(shù)形成節(jié)點導(dǎo)納陣，追加接地支路 ****//int j1。i {j = idk[i]。kk = 0。fprintf(fp2, “*******GII,BII********”)。i {if(i % 2 == 1)fprintf(fp2, “n”)。double data[300][8]。data[ii][1]=。data[ii][5]=。coutwhile(1){double a。double e。if(x!=0amp。x!=2amp。if(x==1){coutcinab。N[i].Qn=b。data[ii][3]=N[i].Pn。data[ii。data[ii][5]=。data[ii][1]=。int i=0。x!=3){coutcoutcinx。x!=1amp。int x。double c。data[ii][7]=。data[ii][3]=Pd。ii=0。} }第四篇：電力系統(tǒng)潮流計算程序include include include “” include “” include “”void main(){ coutPdQd。fprintf(fp2, “n*******YYYYY********”)。} // 打印導(dǎo)納矩陣。yb[ll] = yb[ll]bij / yk。int i, j, i1, l, ll, kk = 0。fprintf(fp2, “*******GII(1),BII(1)********n”)。bii[j] = bii[j] + bij。yg[ll] = yg[ll]bij。zf = r * r + x * x。r = zr[l]。// 支路右節(jié)點號。j1 = jzl[l]。} //導(dǎo)納陣非零非對角元(與支路一一對應(yīng))先清零。}// 導(dǎo)納陣對角元(與節(jié)點一一對應(yīng))先清零。int i, j, i1, l, ll。nob[j] = k。} for(i = 1。}} for(i = 1。zx[ip])。zr[i], amp。iswap(amp。k3 = k4。if(k2 k1 ||(k2 == k1 amp。k3 = iabs(jzl[i])。pji =fii * ejvi * vi * yk。}de = eifj。fi = fi / yk。fj = vj * sin(cj)。//支路左、右節(jié)點電壓值由極坐標(biāo)轉(zhuǎn)換為直角坐標(biāo) ei = vi * cos(ci)。ci = va[i]。x = zx[l]。//轉(zhuǎn)換為新節(jié)點號。jj = iabs(jzl[l])。dkk[i] = v[j] * v[j] / dkk[i]。//統(tǒng)計系統(tǒng)無功網(wǎng)損。double de, df, ei, ej, fi, fj, fii, fir, pij, pji, qij, qji。} void branch_output(){//**** 本函數(shù)輸出支路數(shù)據(jù)。fprintf(fp2, “%%%%”, pl[j], ql[j]，pg[j], qg[j])。if(vi vmin){vmin = vi。oldnumber = newsort[i]。vmin = v[1]。printc(39。} elseb[i]=。dsd = i。 kq0 == 0)pg[i] = w[i]w[i]。pg0 = qg[i]。// vj 此時表示的是節(jié)點正常電壓的會定值。// v0[]存放節(jié)點初始電壓幅值。for(i = 1。double vi, vj。// v[]存放節(jié)點電壓幅值。bi = yb[k]。for(k = i1。break。old[i] = ii1。//nomax 即是最大節(jié)點號。old[j])。i {k1 = i + 1。i {old[i] = iabs(izl[i])。i old[i] = nnew[i] = 0。int i, j, l, k1, k。swap(amp。zr[ip])。jzl[i], amp。}}if(i!= ip){iswap(amp。amp。for(j = i + 1。bi = bi1。b[k] = b[k]1。i2 = nfd[i + 1]。int i, j, k, i1, i2。fprintf(fp2, “%”, aa[i])。printf1(di, n)。for(i = 1。}else{fd[i] = 0。iu[x] = j。k = 0。}break。i2 = ydz[i + 1]1。amp。//|b[i] = bii[i]。if(((t!= 1)||(nobt[i]!=1))amp。double ai, b[NS]。,t=1求B39。b[k1] = aa[i]。for(i = 1。*39。//branch data printc(39。}void out(){//**** 本函數(shù)輸出節(jié)點和支路數(shù)據(jù) ****//zlsort(old)。t ++。i if(kq0 == 0)va[i] = va[i]b[i]。if(fabs(dsm) eps){kg[kq0]=1。do { jdgl(kq0)。yzb(1, iu2, u2, di2, nfd2)。double u1[N2], u2[N2], di1[NS], di2[NS]。dkk[i])。v0[i])。pg[i], amp。i fscanf(fp1, “%d %d”, amp。zr[i], amp。i {fscanf(fp1, “%d %d”, amp。mpj，amp。mdk)。fscanf(fp1, “%d %d %d %d”, amp。fp1 = fopen(inname, “r”)。//af[0]和af[1]分別是負(fù)荷有功功率、無功功率靜態(tài)特性系數(shù)。fprintf(fp2, “%8d%8d%”, idk[i], old[idk[i]], dkk[i])。, 78)。fprintf(fp2, “%8d%8d%8d”, nob[i],old[nob[i]], nobt[i])。, 78)。fprintf(fp2, “%8d%8d%8d%8d”, izl[i], jzl[i], old[abs(izl[i])], old[abs(jzl[i])])。, 78)。fprintf(fp2, “n******AF AND V0 ******n”)。printo()。for(i = 1。newval(qg)。yds[k1] = nobt[i]。// exchange the node before and after sort for(i = 1。} for(i = 1。//mpj:平衡節(jié)點 zlsort(nnew)。i {old[i] = i。dataio()。} } //電抗器節(jié)點號轉(zhuǎn)變 for(j = 1。} printf(“nob[%d] = %dn”, i, nob[i])。printf(“izl[%d] = %d, jzl[%d] = %dn”, i, izl[i], i, jzl[i])。k = jzl[i]。i k = izl[i]。j {if(newsort[i] newsort[j]){k = newsort[i]。ii2)){ii1 ++。ii2)){ii1 ++。ii1 = 0。}void oldtonew(){ //**** 本函數(shù)將輸入數(shù)據(jù)中的節(jié)點號變成從1開始的連續(xù)節(jié)點號 ****//int i, j, k, ii1, ii2, zls2, k1, k2, k3, k4, ip。for(i = 1。//輸出節(jié)點和支路數(shù)據(jù)return 1。void newtoold()。void yzb(int t, int* iu, double* u, double* di, int* nfd)。void calc(int* iu, double* u, double* di, int* nfd, double* b)。void printc()。void calc()。void ya0()。void zlsort(int* nnew)。void pqflow()。// 支路電阻zr[] 支路電抗zx[] 輸電線路充電容納zyk[](y0/2)// 接地電抗dkk[] 對角元實部gii[] 對角元虛部 // 非對角元實部yg[] 非對角元虛部yb[] // pg[],qg[],pl[],ql[]:發(fā)電機，負(fù)荷功率實、虛部 // v[]是電壓幅值,va[]是電壓相角。double v00。double eps, dsm, vmin, dph, dqh, af[3]。int nnew[NS4], old[NS], nob[NS], nobt[NS]。// 節(jié)點總數(shù)n(包括聯(lián)絡(luò)節(jié)點)支路數(shù)(回路數(shù))zls 節(jié)點數(shù)nb(發(fā)電機和負(fù)荷)// 接地電抗數(shù)mdk 精度eps平衡節(jié)點號mpj// 節(jié)點優(yōu)化(標(biāo)志)bnsopton(=0節(jié)點不優(yōu)化，!=0節(jié)點優(yōu)化)// 最大迭代次數(shù)it1 最低電壓或最大功率誤差節(jié)點號dsd // 負(fù)荷靜特性標(biāo)志(=0考慮負(fù)荷靜特性)// 支路數(shù)(雙回線算一條支路)int izl[ZS], jzl[ZS], idk[DKS], yds[NS], ydz[NS], iy[ZS2]。00 0======線路傳輸功率========== 2to103to14to204to35to206to206to407to507to608to606to906to10011to9010to904to12013to12 14to1230to29 15to1228to8 28