【正文】 。} } for(i = 1。ii2)){ii1 ++。}k = jzl[i]。ii2)){ii1 ++。i k = izl[i]。ii1 = 0。for(i = 1。}void oldtonew(){ //**** 本函數(shù)將輸入數(shù)據(jù)中的節(jié)點號變成從1開始的連續(xù)節(jié)點號 ****//int i, j, k, ii1, ii2, zls2, k1, k2, k3, k4, ip。return 1。for(i = 1。return a。//輸出節(jié)點和支路數(shù)據(jù)return 1。//初始化pqflow()。void newtoold()。void yy1()。void yzb(int t, int* iu, double* u, double* di, int* nfd)。void printf1(double* aa, int n)。void calc(int* iu, double* u, double* di, int* nfd, double* b)。void swap()。void printc()。void branch_output()。void calc()。void jdgl(int kq0)。void ya0()。void printy()。void zlsort(int* nnew)。void dataio()。void pqflow()。int newsort[NS4]。// 支路電阻zr[] 支路電抗zx[] 輸電線路充電容納zyk[](y0/2)// 接地電抗dkk[] 對角元實部gii[] 對角元虛部 // 非對角元實部yg[] 非對角元虛部yb[] // pg[],qg[],pl[],ql[]:發(fā)電機(jī)，負(fù)荷功率實、虛部 // v[]是電壓幅值,va[]是電壓相角。double zr[ZS], zx[ZS], zyk[ZS], dkk[DKS], gii[NS], bii[NS], yg[ZS2], yb[ZS2]。double v00。dph,dqh:系統(tǒng)有、無功損耗。double eps, dsm, vmin, dph, dqh, af[3]。// nnew[i]中為i對應(yīng)的新號 // nob[]存放的是節(jié)點號。int nnew[NS4], old[NS], nob[NS], nobt[NS]。// yds[]存放各行非零非對角元素的個數(shù)。// 節(jié)點總數(shù)n(包括聯(lián)絡(luò)節(jié)點)支路數(shù)(回路數(shù))zls 節(jié)點數(shù)nb(發(fā)電機(jī)和負(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]。char inname[12], outname[12]。00 0======線路傳輸功率========== 2to103to14to204to35to206to206to407to507to608to606to906to10011to9010to904to12013to12 14to1230to29 15to1228to8 28to6 請按任意鍵繼續(xù)... 20to19 20to10 22to21 26to25 27to25 28to27 高等電力系統(tǒng)分析 IEEE30節(jié)點潮流程序班級：電研114班姓名：王大偉學(xué)號：2201100151第三篇：電力系統(tǒng)潮流計算的C語言實現(xiàn)////////////////////////////////////////////////////////////////////// //PQ分解法潮流// //文件輸入格式：節(jié)點總數(shù)n(包括聯(lián)絡(luò)節(jié)點),支路數(shù)zls//節(jié)點數(shù)(發(fā)電機(jī)和負(fù)荷)nb,接地電抗數(shù)mdk,迭代精度eps//考慮負(fù)荷靜特性標(biāo)志kk2(0考慮),平衡節(jié)點號，優(yōu)化標(biāo)志(0不優(yōu)化)//最大迭代次數(shù)it1,支路左右節(jié)點號izl[],jzl[],支路電阻zr[],電抗zx[] // //支路容納zyk[],節(jié)點號nob[]及標(biāo)志nobt[](0PQ,1PV)//發(fā)電機(jī)和負(fù)荷有功、無功pg[],qg[],pl[],ql[]//電壓v0[](pv節(jié)點輸入實際值，PQ節(jié)點任輸入一值)//電抗節(jié)點號idk[],電抗值dkk[]//////////////////////////////////////////////////////////////////////include “” include “”define NS 2000 //最大節(jié)點數(shù) define NS2 NS * 2 define NS4 define ZS define ZS2 1000 3000 ZS * 2//NSNS必須大于2*zls。。00 節(jié)點8 節(jié)點9 0 節(jié)點10 0 節(jié)點11 節(jié)點12 0 節(jié)點13 節(jié)點14 0 節(jié)點15 0 節(jié)點16 0 節(jié)點17 0 節(jié)點18 0 節(jié)點19 0 節(jié)點20 0 節(jié)點21 0 節(jié)點22 0 節(jié)點23 0 節(jié)點24 0 。00 0 節(jié)點7。0 0 節(jié)點6。0 00 節(jié)點5。00 0 節(jié)點4。節(jié)點3。====節(jié)點電壓===============發(fā)電機(jī)發(fā)出功率======節(jié)點10。x[pivcol[k][0]]=al。k){al=x[pivcol[k][1]]。} for(k=(n2)。j{sum +=a[i][j]*x[j]。i){sum=0。//解方程for(i=(n2)。j{a[i][j]=aik*a[k][j]/a[k][k]。i{aik=a[i][k]。a[i][k]=al。i{al=a[i][pivcol[k][1]]。a[k][j]=al。j{al=a[pivrow[k]][j]。exit(0)。//行pivcol[k][1]=j。j{if(pivot{pivot=fabs(a[i][j])。//列n*2矩陣for(i=k。//行pivcol[k][0]=k。ipivot= fabs(a[k][k])。pivcol=new int *[n]。n=NN。void gauss::gauss_slove(double **a,double *x,int NN){int n,i,j,k,*pivrow,**pivcol。static void gauss_output()。=0。cout”。(8)。cout”。(8)。}doubleComplex::getCompleximage(Complex c1)//取虛部 {return 。return Node。=0。}Complex Complex::getinverse(Complex c1)//取倒數(shù) { Complex Node。=。} Complex Complex::getconj(Complex c1)//取共軛 {Complex Node。=(r1)。} Complex Complex::ComDivRea(Complex c1,double r1)//復(fù)數(shù)除數(shù) { Complex Node。=(**)/(pow(,2)+pow(,2))。} Complex Complex::divideComplex(Complex c1,Complex c2)//復(fù)數(shù)除法 {Complex Node。=*+*。}Complex Complex::productComplex(Complex c1,Complex c2)//復(fù)數(shù)乘法 {Complex Node。=。}Complex Complex::subComplex(Complex c1,Complex c2)//復(fù)數(shù)減法 {Complex Node。=+。}Complex Complex::CaddC(Complex c1,Complex c2)//復(fù)數(shù)加法 {Complex Node。=1。=sqrt(*+*)。} }。//將復(fù)數(shù)復(fù)零static Complex Rec2Polar(Complex c1)。//顯示一個復(fù)數(shù)static void PrintfmultiComplex(Complex C,int N)。//求一個復(fù)數(shù)實部static double getCompleximage(Complex c1)。//求一個復(fù)數(shù)共軛static Complex getinverse(Complex c1)。//求兩個復(fù)數(shù)商static Complex ComDivRea(Complex c1,double r2)。//求兩個復(fù)數(shù)差static Complex productComplex(Complex c1,Complex c2)。//0表示直角坐標(biāo)，1表示極坐標(biāo)static Complex CaddC(Complex c1,Complex c2)。double image。}include using namespace std。Complex::PrintfComplex(Complex::ComDivRea(line[i].etos,))。endtemp=line[i].end。i{int statemp,endtemp。coutt(bus[i].)。i{coutComplex::PrintfComplex(bus[i].V)。bus[i].V=Complex::Rec2Polar(bus[i].V)。i{bus[i]. = bus[i]. + bus[i].。}coutBus::JisuanNodeScal(X,bus,N)。dd=Complex::CaddC(aa,cc)。bb=Complex::subComplex(Complex::getconj(bus[endtemp1].V), Complex::getconj(bus[statemp1].V))。line[i].stoe=Complex::productComplex(bus[statemp1].V,dd)。cc=Complex::productComplex(bb , Complex::getconj(line[i].Y))。aa=Complex::productComplex(Complex::getconj(bus[statemp1].V), B)。statemp=line[i].start。=0。i{int statemp,endtemp。goto LOOP。//修正節(jié)點電壓// Bus::PrintfNodeV(bus,N)。gauss::gauss_slove(JacAug,x,NN)。//whether converbence看迭代是否結(jié)束if(icon==1){coutBus::JisuanJacAug(JacAug,X,bus,N)。//計算節(jié)點功率差值 Bus::PrintfNodeScal(X,bus,N)。//計算節(jié)點功率Bus::JisuanNodeScal(X,bus,N)。LOOP:Bus::JisuanNodeI(X,bus,N)。x=new double[NN]。for(i=0。double **JacAug。//計算節(jié)點導(dǎo)納矩陣Bus::PrintfNodeDnz(X,N)。for(i=0。//輸出結(jié)點參數(shù)Complex **X。}Bus::ScanfBusData(bus)。ibus[i].=0。//輸出支路參數(shù)Bus *bus=new Bus[N]。//動態(tài)分配支路結(jié)構(gòu)體Line::ScanfLineData(line)。N)==0){return 0。//輸入支路個數(shù)if((Lamp。//i作為整個程序的循環(huán)變量int N=Bus::ScanfBusNo()。}include “” include include include include“” include“” include“” using namespace std。Complex::PrintfComplex(Complex::ComDivRea(line[i].etos,))。endtemp=line[i].end。i{int statemp,endtemp。coutt(bus[i].)。i{coutComplex::PrintfComplex(bus[i].V)。bus[i].V=Complex::Rec2Polar(bus[i].V)。i{bus[i]. = bus[i]. + bus[i].。}coutBus::JisuanNodeScal(X,bus,N)。dd=Complex::CaddC(aa,cc)。bb=Complex::subComplex(Compl