【正文】 統(tǒng)潮流計(jì)算c語(yǔ)言程序，兩行，大家可以看看，仔細(xì)研究，然后在這個(gè)基礎(chǔ)上修改。//輸入節(jié)點(diǎn)個(gè)數(shù)int L=Line::ScanflineNo()。} //如果找不到兩個(gè)文件中的任意一個(gè)，退出Line *line=new Line[L]。//動(dòng)態(tài)分配結(jié)點(diǎn)結(jié)構(gòu)體for(int i=0。//輸入節(jié)點(diǎn)參數(shù)Bus::PrintfBusData(bus,N)。iBus::JisuanNodeDnz(X,line,bus,L,N)。JacAug=new double *[NN]。int count=1。//計(jì)算節(jié)點(diǎn)功率Bus::JisuanNodeSdelta(bus,N)。//計(jì)算雅可比增廣矩陣 // Bus::PrintfJacAug(JacAug,N)。count++。Complex aa,bb,cc,dd,B。endtemp=line[i].end。dd=Complex::CaddC(aa,cc)。cc=Complex::productComplex(bb , Complex::getconj(line[i].Y))。//計(jì)算節(jié)點(diǎn)功率for(i=0。}coutfor(i=0。cout}coutfor(i=0。coutComplex::PrintfComplex(Complex::ComDivRea(line[i].stoe,))。int _tmain(int argc, _TCHAR* argv[]){int i。amp。//輸入支路參數(shù)Line::PrintfLineData(line,L)。bus[i].=0。X=new Complex *[N]。//輸出節(jié)點(diǎn)導(dǎo)納矩陣int NN=(N1)*2。idouble *x。//計(jì)算節(jié)點(diǎn)注入電流Bus::JisuanNodeScal(X,bus,N)。//輸出節(jié)點(diǎn)功率差值int icon=wehcon1(bus,N)。//解方程組求出電壓差值Bus::ReviseNodeV(bus,x,N)。}else{for(i=0。=line[i].B。bb=Complex::subComplex(Complex::getconj(bus[statemp1].V), Complex::getconj(bus[endtemp1].V))。aa=Complex::productComplex(Complex::getconj(bus[endtemp1].V), B)。line[i].etos=Complex::productComplex(bus[endtemp1].V,dd)。//發(fā)電機(jī)功率=注入功率+負(fù)荷功率bus[i].= bus[i].+ bus[i].。coutt(bus[i].)。statemp=line[i].start。cout} }return 0。intRecPolar。//求兩個(gè)復(fù)數(shù)積static Complex divideComplex(Complex c1,Complex c2)。//取倒數(shù)static double getComplexReal(Complex c1)。//顯示多個(gè)復(fù)數(shù)static void zeroComplex(Complex c1)。Complex Complex::Rec2Polar(Complex c1)//極坐標(biāo)表示 { Complex Node。return Node。return Node。return Node。return Node。return Node。return Node。return Node。Node=(Complex::divideComplex(Node,c1))。}voidComplex::PrintfComplex(Complex c1)//按直角坐標(biāo)輸出 { if(==0){ (6)。(6)。} else {coutComplex::zeroComplex(Complex c1)//清零 { =0。}。pivrow=new int[n]。pivrow[k]=k。i{for(j=k。//列}}}if(pivot{coutgetchar()。a[pivrow[k]][j]=a[k][j]。a[i][pivcol[k][1]]=a[i][k]。for(j=k。i=0。 }x[i]=(a[i][n]sum)/a[i][i]。x[pivcol[k][1]]=x[pivcol[k][0]]。節(jié)點(diǎn)2。0節(jié)點(diǎn)26 。0節(jié)點(diǎn)28 。0 20節(jié) 點(diǎn)30 。// fp1輸入數(shù)據(jù)文件指針 fp2輸出文件指針// inname[]輸入數(shù)據(jù)文件名 outname[]輸出數(shù)據(jù)文件名 int n, zls, nb, mdk, mpj, bnsopton, it1, dsd, kk2, nzls。// ydz[i]是第 i 行第一個(gè)非零非對(duì)角元素的首地址，// 即在所有非零非對(duì)角元素中的次序號(hào) // iy[]存放列足碼。nobt[]存放的是節(jié)點(diǎn)類型, 0: pq節(jié)點(diǎn),1: pv節(jié)點(diǎn)。// af[0]和af[1]分別是負(fù)荷有功功率、無(wú)功功率靜態(tài)特性系數(shù)。double pg[NS], qg[NS], pl[NS], ql[NS], v0[NS], v[NS], va[NS]。// newsort[i]存放i對(duì)應(yīng)的老號(hào)void initial()。void bnsopt()。void y2()。void bbhl(int kq0)。void newval(double* aa)。void printf2(double* aa, double* bb, int n)。int find(int k, int a[], int* z)。void y3()。//pq潮流計(jì)算out()。}int find(int k, int a[], int* z){****// //**** 本函數(shù)查找a[]中是否有fabs(k)有則返回0，無(wú)則返回1 ****// //參數(shù)1為待查找量，參數(shù)2待搜索數(shù)組，參數(shù)3返回k在a[]中的次序號(hào)//int i。} return 0。i newsort[i] = 0。if(!find(k, newsort, amp。if(!find(k, newsort, amp。i for(j = i+1。}} } for(i = 1。} } elseprintf(“error!”)。} elseprintf(“error!”)。break。break。} } } void initial(){//**** 本函數(shù)進(jìn)行初始化工作 ****//int i, k1。for(i = 1。//節(jié)點(diǎn)優(yōu)化mpj = nnew[mpj]。idk[i] = nnew[k1]。} // (v00)，電壓相角初值都為0。i {k1 = nnew[nob[i]]。newval(pg)。newval(v0)。//nobt[] = 0: pq節(jié)點(diǎn)，v0[]存放的是最后一個(gè)節(jié)點(diǎn)數(shù)據(jù)，//對(duì)于pq節(jié)點(diǎn)，即為系統(tǒng)平均電壓值。(P407)}void printo(){ //**** 輸出af[]、v00和節(jié)點(diǎn)排序后的支路、節(jié)點(diǎn)和接地電抗數(shù)據(jù) ****//int i。39。i {fprintf(fp2, “n”)。39。i {fprintf(fp2, “n”)。39。i {fprintf(fp2, “n”)。af[1] = 。scanf(“%s”, inname)。fp2 = fopen(outname, “w”)。nb, amp。kk2, amp。//precision, swing node,sort the node,iteration numbers for(i = 1。fscanf(fp1, “%lf %lf %lf ”, amp。} for(i = 1。fscanf(fp1, “%lf %lf %lf %lf %lf”, amp。ql[i], amp。idk[i], amp。int i, t。//form the B matrix of P0 iteration y2()。kg[0] = kg[1] = 1。elseprintf(“Q: %dt%dt%fn”, t, dsd, dsm)。for(i = 1。else{kq0 = 0。fprintf(fp2, “n%s%d”, “times = ”, t)。// node data branch_output()。printc(39。}void newval(double* aa){//**** 本函數(shù)將舊號(hào)換成新號(hào) ****//int i, k1。i {k1 = nnew[nob[i]]。}void yzb(int t, int* iu, double* u, double* di, int* nfd){ //**** 本函數(shù)求因子表 ****////參數(shù)1為標(biāo)志(t=0 求B39。int jj, jj1, jj2, im, x, fd[NS]。i //nobt[] 存放的是節(jié)點(diǎn)類型, 0: pq節(jié)點(diǎn),1: pv節(jié)點(diǎn)。j //|b[j] = 。(t == 1)amp。//af[1] i1 = ydz[i]。jj2 = nfd[im + 1]ai * u[k]。ai = di[i]。j {if(fabs(b[j]) ){u[x] = b[j] * ai。}}fd[i] = k。} fprintf(fp2, “n********U*********”)。} fprintf(fp2, “n********DI********”)。i {if(i % 5 == 1)fprintf(fp2, “n”)。(P417圖F19)double bi。i1 = nfd[i]。j ++){k = iu[j]。j){k = iu[j]。k3 = iabs(jzl[i])。if(k2 k1 ||(k2 == k1 amp。k3 = k4。}} }void bnsopt(){ iswap(amp。zr[i], amp。zx[ip])。//**** 節(jié)點(diǎn)優(yōu)化 ****//int ii1, ii2, zls2, nomax。for(i = 1。for(i = 1。for(i = 1。old[i], amp。}nomax = old[zls2]。i {ii1 = old[l]。if(ii1!= ii2){l = j。i //即為第 i 行的非零元素個(gè)數(shù)。}//yb[]存放導(dǎo)納陣非對(duì)角元的虛部(yg[] + j*yb[])else{ai = yg[k]。vj = vi * v[j]。(程序框圖見P423)****//int i。// dsm 即為最大功率誤差。vj = v0[i]。// nopt[] = 0: pq節(jié)點(diǎn)，對(duì)于pq節(jié)點(diǎn)，即為系統(tǒng)平均電壓值。}// 負(fù)荷p,發(fā)電機(jī)pelse{pl0 = ql[i]。amp。(i!= mpj)){ if(fabs(b[i]) fabs(dsm)){dsm = b[i]。// 計(jì)算修正方程式的常數(shù)項(xiàng)。int i, j, oldnumber。fprintf(fp2, “n%5s%8s%10s%11s%11s%11s%11sn”, “I”，“V”, “CA”, “PL”, “QL”, “PG”, “QG”)。i j = nnew[i]。vi = v[j]。fprintf(fp2, “n%5d%%”, oldnumber, vi, ci)。, 72)。int i, j, l。dqh = 。i {j = idk[i]。// izl[]: 支路左節(jié)點(diǎn)號(hào)。j=nnew[jj]。//轉(zhuǎn)化為相應(yīng)舊號(hào) r = zr[l]。//v[]: 電壓幅值。cj = va[j]。ej = vj * cos(cj)。{ if(izl[l] 0){ei = ei / yk。}yk = 。fir =(df * rfir * ei。k1 = iabs(izl[i])。k4 = iabs(jzl[j])。k1 = k2。izl[ip])。swap(amp。zx[i], amp。zyk[ip])。printf(“nob[%d] = %dn”, i, nob[i])。nob[i] = nob[j]。double r, x, yk, zf, gij, bij。bii[i] = 。yb[i] = 。l i1 = izl[l]。if(i == j)// 支路左節(jié)點(diǎn)號(hào)。//排除左、右節(jié)點(diǎn)號(hào)相等的情況。// zr[],zx[],zyk[]:支路三參數(shù)。bij =1/x。gii[j] = gii[j] + gij。對(duì)角元實(shí)部為gii，虛部為bii，// 非零非對(duì)角元實(shí)部為yb[]，虛部為yb[]，列足碼為iy[]。double r, x, yk, zf, gij, bij。bii[j] = bii[j]gij / yk。}elsekk = 1。printf2(gii,bii,n)。fprintf(fp2, “%%%8d”, yg[i], yb[i], iy[i])。//分別為0前面有功，1前面無(wú)功，2壓降，3中間有功，4中間無(wú)功,5變比,6后面有功，7后面無(wú)功int ii。data[ii][2]=Ud。data[ii][6]=。double b。double f。amp。amp。Node N[100]。data[ii][0]=。data[ii][4]=N[i].Q