【正文】 plex c2)。//求兩個復數(shù)差static Complex productComplex(Complex c1,Complex c2)。//求兩個復數(shù)積static Complex divideComplex(Complex c1,Complex c2)。//求兩個復數(shù)商static Complex ComDivRea(Complex c1,double r2)。//除數(shù)static Complex getconj(Complex c1)。//求一個復數(shù)共軛static Complex getinverse(Complex c1)。//取倒數(shù)static double getComplexReal(Complex c1)。//求一個復數(shù)實部static double getCompleximage(Complex c1)。//求一個復數(shù)虛部static void PrintfComplex(Complex c1)。//顯示一個復數(shù)static void PrintfmultiComplex(Complex C,int N)。//顯示多個復數(shù)static void zeroComplex(Complex c1)。//將復數(shù)復零static Complex Rec2Polar(Complex c1)。//取極坐標Complex(){RecPolar=0。} }。Complex Complex::Rec2Polar(Complex c1)//極坐標表示 { Complex Node。=sqrt(*+*)。=atan2(,)*180/。=1。return Node。}Complex Complex::CaddC(Complex c1,Complex c2)//復數(shù)加法 {Complex Node。=+。=+。return Node。}Complex Complex::subComplex(Complex c1,Complex c2)//復數(shù)減法 {Complex Node。=。=。return Node。}Complex Complex::productComplex(Complex c1,Complex c2)//復數(shù)乘法 {Complex Node。=**。=*+*。return Node。} Complex Complex::divideComplex(Complex c1,Complex c2)//復數(shù)除法 {Complex Node。=(*+*)/(pow(,2)+pow(,2))。=(**)/(pow(,2)+pow(,2))。return Node。} Complex Complex::ComDivRea(Complex c1,double r1)//復數(shù)除數(shù) { Complex Node。=(r1)。=(r1)。return Node。} Complex Complex::getconj(Complex c1)//取共軛 {Complex Node。=。=。return Node。}Complex Complex::getinverse(Complex c1)//取倒數(shù) { Complex Node。=1。=0。Node=(Complex::divideComplex(Node,c1))。return Node。}double Complex::getComplexReal(Complex c1)//取實部 {return 。}doubleComplex::getCompleximage(Complex c1)//取虛部 {return 。}voidComplex::PrintfComplex(Complex c1)//按直角坐標輸出 { if(==0){ (6)。(8)。(ios::right)。cout”。(6)。(8)。(ios::left)。cout”。} else {coutComplex::zeroComplex(Complex c1)//清零 { =0。=0。}class gauss { public: static void gauss_slove(double **a,double *x,int NN)。static void gauss_output()。}。void gauss::gauss_slove(double **a,double *x,int NN){int n,i,j,k,*pivrow,**pivcol。double eps,pivot,sum,aik,al。n=NN。pivrow=new int[n]。pivcol=new int *[n]。for(i=0。ipivot= fabs(a[k][k])。pivrow[k]=k。//行pivcol[k][0]=k。pivcol[k][1]=k。//列n*2矩陣for(i=k。i{for(j=k。j{if(pivot{pivot=fabs(a[i][j])。pivrow[k]=i。//行pivcol[k][1]=j。//列}}}if(pivot{coutgetchar()。exit(0)。}if(pivrow[k]!=k)//行變換{for(j=k。j{al=a[pivrow[k]][j]。a[pivrow[k]][j]=a[k][j]。a[k][j]=al。}}if(pivcol[k][1]!=k)//列變換{for(i=0。i{al=a[i][pivcol[k][1]]。a[i][pivcol[k][1]]=a[i][k]。a[i][k]=al。}}if(k!=(n1))//將矩陣化為上三角形式{for(i=(k+1)。i{aik=a[i][k]。for(j=k。j{a[i][j]=aik*a[k][j]/a[k][k]。}}} } x[n1]=a[n1][n]/a[n1][n1]。//解方程for(i=(n2)。i=0。i){sum=0。for(j=(i+1)。j{sum +=a[i][j]*x[j]。 }x[i]=(a[i][n]sum)/a[i][i]。} for(k=(n2)。k=0。k){al=x[pivcol[k][1]]。x[pivcol[k][1]]=x[pivcol[k][0]]。x[pivcol[k][0]]=al。}coutcout}====節(jié)點功率計算值==== icon=1進行第2次迭代 節(jié)點電壓修正量================= ====節(jié)點功率計算值==== icon=1進行第3次迭代 節(jié)點電壓修正量================= ====節(jié)點功率計算值==== icon=0,迭代結束。====節(jié)點電壓===============發(fā)電機發(fā)出功率======節(jié)點10。節(jié)點2。節(jié)點3。0節(jié)點25 。00 0 節(jié)點4。0節(jié)點26 。0 00 節(jié)點5。節(jié) 點27 。0 0 節(jié)點6。0節(jié)點28 。00 0 節(jié)點7。0節(jié)點29 。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 。。。0 20節(jié) 點30 。00 0======線路傳輸功率========== 2to103to14to204to35to206to206to407to507to608to606to906to10011to9010to904to12013to12 14to1230to29 15to1228to8 28to6 請按任意鍵繼續(xù)... 20to19 20to10 22to21 26to25 27to25 28to27 高等電力系統(tǒng)分析 IEEE30節(jié)點潮流程序班級：電研114班姓名：王大偉學號：2201100151第二篇：電力系統(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[ii][7]=。i++。ii++。}if(x==2){coutcinabc。Transmission_line Tm[100]。int j=0。Tm[j].Rl=a。Tm[j].Xl=b。Tm[j].Bl=c。data[ii][0]=。data[ii][1]=Tm[j].GetQf(data[ii1][2])。data[ii][6]=。data[ii][7]=Tm[j].GetQb(data[ii1][2])。data[ii][2]=Tm[j].GetUl(data[ii1][3]+data[ii1][1]+data[ii][7],data[ii1][4]+data[ii1][0]+data[ii][6],data[ii1][2])。//參數(shù)分別為前節(jié)點的有功、無功、電壓data[ii][3]=Tm[j].GetPl(data[ii1][3]+data[ii1][1]+dat