【正文】 data[ii][0]=。double f。//分別為0前面有功，1前面無(wú)功，2壓降，3中間有功，4中間無(wú)功,5變比,6后面有功，7后面無(wú)功int ii。bii[j] = bii[j]gij / yk。bij =1/x。l i1 = izl[l]。nob[i] = nob[j]。swap(amp。k1 = iabs(izl[i])。ej = vj * cos(cj)。j=nnew[jj]。int i, j, l。i j = nnew[i]。(i!= mpj)){ if(fabs(b[i]) fabs(dsm)){dsm = b[i]。vj = v0[i]。}//yb[]存放導(dǎo)納陣非對(duì)角元的虛部(yg[] + j*yb[])else{ai = yg[k]。}nomax = old[zls2]。for(i = 1。}} }void bnsopt(){ iswap(amp。j){k = iu[j]。i {if(i % 5 == 1)fprintf(fp2, “n”)。j {if(fabs(b[j]) ){u[x] = b[j] * ai。(t == 1)amp。}void yzb(int t, int* iu, double* u, double* di, int* nfd){ //**** 本函數(shù)求因子表 ****////參數(shù)1為標(biāo)志(t=0 求B39。// node data branch_output()。elseprintf(“Q: %dt%dt%fn”, t, dsd, dsm)。idk[i], amp。fscanf(fp1, “%lf %lf %lf ”, amp。fp2 = fopen(outname, “w”)。39。39。newval(pg)。//節(jié)點(diǎn)優(yōu)化mpj = nnew[mpj]。break。i for(j = i+1。} return 0。int find(int k, int a[], int* z)。void y2()。// af[0]和af[1]分別是負(fù)荷有功功率、無(wú)功功率靜態(tài)特性系數(shù)。0 20節(jié) 點(diǎn)30 。x[pivcol[k][1]]=x[pivcol[k][0]]。a[i][pivcol[k][1]]=a[i][k]。pivrow[k]=k。(6)。return Node。return Node。//取倒數(shù)static double getComplexReal(Complex c1)。statemp=line[i].start。aa=Complex::productComplex(Complex::getconj(bus[endtemp1].V), B)。//解方程組求出電壓差值Bus::ReviseNodeV(bus,x,N)。//輸出節(jié)點(diǎn)導(dǎo)納矩陣int NN=(N1)*2。amp。}coutfor(i=0。endtemp=line[i].end。//計(jì)算節(jié)點(diǎn)功率Bus::JisuanNodeSdelta(bus,N)。//輸入節(jié)點(diǎn)參數(shù)Bus::PrintfBusData(bus,N)。}第二篇：電力系統(tǒng)潮流計(jì)算程序電力系統(tǒng)潮流計(jì)算c語(yǔ)言程序，兩行，大家可以看看，仔細(xì)研究，然后在這個(gè)基礎(chǔ)上修改。}//全網(wǎng)數(shù)據(jù)// int high=1,low=1。QZL[i][j]=(x1*y2x2*y1)。}}//平衡節(jié)點(diǎn)//for(i=0。//調(diào)用Jacobi函數(shù)和高斯函數(shù)//for(int t=1。K[i][j]=1。ZLs++。B[i].Pd=m3[i]。Generator*G。Ja[2*i+1][2*j+1]=2*f[i]。jif(i!=j){if(type[i]==1){Ja[2*i][2*j]=(YG[i][j]*e[i]+YB[i][j]*f[i])。//deltaP deltaQ deltaV^2// void fun1(double YG[][50],double YB[][50],double e[],double f[],int type[],int N,double W[],double P[],double Q[],double V[]){ double dP=0,dQ=0,dV=0。//節(jié)點(diǎn)號(hào)類型負(fù)荷有功負(fù)荷無(wú)功母線數(shù)據(jù)（類型 1＝PV節(jié)點(diǎn)，2＝PQ節(jié)點(diǎn)，3＝平衡節(jié)點(diǎn)）struct BUS { int busno。float B。} } //Jacobi矩陣// void Jacobi(double YG[][50],double YB[][50],double e[50],double f[50],int type[50],int N ,double Ja[100][101]){int i,j。if(type[i]==1){Ja[2*i+1][2*j]=b[i]+YB[i][i]*e[i]YG[i][i]*f[i]。jvoid main(){ ifstream fin。for(i=0。i++){finm1[i]。iYB[i][j]=0。iC[2*i]=1。for(i=0。y1=f[i]/L[k].kf[j]。pr[j][i]=I2*L[k].R。ki=L[k].busi1。ibus[i].=0。LOOP:Bus::JisuanNodeI(X,bus,N)。=0。i{bus[i]. = bus[i]. + bus[i].。//i作為整個(gè)程序的循環(huán)變量int N=Bus::ScanfBusNo()。for(i=0。//whether converbence看迭代是否結(jié)束if(icon==1){coutBus::JisuanJacAug(JacAug,X,bus,N)。cc=Complex::productComplex(bb , Complex::getconj(line[i].Y))。coutt(bus[i].)。//求兩個(gè)復(fù)數(shù)商static Complex ComDivRea(Complex c1,double r2)。}Complex Complex::CaddC(Complex c1,Complex c2)//復(fù)數(shù)加法 {Complex Node。} Complex Complex::ComDivRea(Complex c1,double r1)//復(fù)數(shù)除數(shù) { Complex Node。(8)。pivcol=new int *[n]。a[k][j]=al。} for(k=(n2)。00 0 節(jié)點(diǎn)7。double eps, dsm, vmin, dph, dqh, af[3]。void zlsort(int* nnew)。void calc(int* iu, double* u, double* di, int* nfd, double* b)。for(i = 1。ii2)){ii1 ++。printf(“izl[%d] = %d, jzl[%d] = %dn”, i, izl[i], i, jzl[i])。i {old[i] = i。yds[k1] = nobt[i]。fprintf(fp2, “n******AF AND V0 ******n”)。fprintf(fp2, “%8d%8d%8d”, nob[i],old[nob[i]], nobt[i])。fp1 = fopen(inname, “r”)。i {fscanf(fp1, “%d %d”, amp。v0[i])。do { jdgl(kq0)。}void out(){//**** 本函數(shù)輸出節(jié)點(diǎn)和支路數(shù)據(jù) ****//zlsort(old)。b[k1] = aa[i]。//|b[i] = bii[i]。k = 0。printf1(di, n)。b[k] = b[k]1。}}if(i!= ip){iswap(amp。int i, j, l, k1, k。old[j])。for(k = i1。for(i = 1。 kq0 == 0)pg[i] = w[i]w[i]。vmin = v[1]。} void branch_output(){//**** 本函數(shù)輸出支路數(shù)據(jù)。jj = iabs(jzl[l])。//支路左、右節(jié)點(diǎn)電壓值由極坐標(biāo)轉(zhuǎn)換為直角坐標(biāo) ei = vi * cos(ci)。pji =fii * ejvi * vi * yk。iswap(amp。} for(i = 1。} //導(dǎo)納陣非零非對(duì)角元(與支路一一對(duì)應(yīng))先清零。zf = r * r + x * x。int i, j, i1, l, ll, kk = 0。} }第四篇：電力系統(tǒng)潮流計(jì)算程序include include include “” include “” include “”void main(){ coutPdQd。double c。int i=0。data[ii][3]=N[i].Pn。if(x!=0amp。data[ii][1]=。kk = 0。bii[i] = bii[i] + bij。j = iabs(j1)。}void yy1(){//**** 本函數(shù)形成節(jié)點(diǎn)導(dǎo)納陣(不包括接地支路)****// int j1。swap(amp。j {k2 = iabs(izl[j])。if((izl[l] 0)||(jzl[l] 0))// 變壓器支路。jj = newsort[jj]。//統(tǒng)計(jì)系統(tǒng)有功網(wǎng)損。//弧度轉(zhuǎn)化為角度。b[i] = b[i] / vi。// vi[] // nopt[] =1: pv節(jié)點(diǎn)，對(duì)于pv節(jié)點(diǎn)，即為該點(diǎn)應(yīng)維持的電壓值。//iy[]存放的是列足碼。for(i = 1。由此可知：NSNS必須大于2*zls。swap(amp。k1 = iabs(izl[i])。}void calc(int* iu, double* u, double* di, int* nfd, double* b){//**** 本函數(shù)利用因子表解線形方程組。x++。b[i] = b[i] + af[1] * ql[i] / v0[i] / v0[i]。)// //參數(shù)2因子表上三角矩陣非零非對(duì)角元素的列足碼 //參數(shù)3因子表上三角矩陣非零非對(duì)角元素的數(shù)值 //參數(shù)4因子表上三角矩陣對(duì)角元素 //參數(shù)5因子表上三角各行非零元素個(gè)數(shù)int i, j, k, i1, i2。, 78)。if(kq0 == 1)calc(iu2, u2, di2, nfd2, b)。}void pqflow(){//**** PQ分解法計(jì)算潮流,程序框圖見(jiàn)P164圖316(從第 7 步起)****//int kq0, iu1[N2], nfd1[NS], iu2[N2], nfd2[NS]。zyk[i])。zls, amp。for(i = 1。for(i = 1。newval(ql)。i k1 = idk[i]。j mpj = j。newsort[j] = k。for(i = 1。void yy1()。void jdgl(int kq0)。double zr[ZS], zx[ZS], zyk[ZS], dkk[DKS], gii[NS], bii[NS], yg[ZS2], yb[ZS2]。char inname[12], outname[12]。====節(jié)點(diǎn)電壓===============發(fā)電機(jī)發(fā)出功率======節(jié)點(diǎn)10。i{aik=a[i][k]。//列n*2矩陣for(i=k。cout”。=。=。//顯示一個(gè)復(fù)數(shù)static void PrintfmultiComplex(Complex C,int N)。Complex::PrintfComplex(Complex::ComDivRea(line[i].etos,))。dd=Complex::CaddC(aa,cc)。goto LOOP。for(i=0。//動(dòng)態(tài)分配支路結(jié)構(gòu)體Line::ScanfLineData(line)。coutt(bus[i].)。cc=Complex::productComplex(bb , Complex::getconj(line[i].Y))。//whether converbence看迭代是否結(jié)束if(icon==1){coutBus::JisuanJacAug(JacAug,X,bus,N)。for(i=0。//i作為整個(gè)程序的循環(huán)變量int N=Bus::ScanfBusNo()。kVh){Vh=sqrt(e[k]*e[k]+f[k]*f[k])。pr[i][j]=I2*L[k].R。ja+=(YG[i][j]*e[j]YB[i][j]*f[j])。f[i]=C[2*i+1]。ifor(j=i。iL[i].busi=m1[i]。i++){finm1[i]。if(!fin){coutgetchar()。for(k=0。Ja[2*i+1][2*j+1]=Ja[2*i][2*j]。idouble A=0,B=0。float Qd。int busj。dQ=Q[i]f[i]*A+e[i]*B。b[i]+=(YG[i][k]*f[k]+YB[i][k]*e[k])。for(k=n2。finm2[i]m3[i]m4[i]。G[i].Pg=m4[i]。()。l// 求A=e+f// double e[50]={0},f[50]={0}。for(i=0。ki=L[k].busi1。QZL[j][i]=(x1*y2x2*y1)。ifout}for(j=0。//輸入支路參