【正文】 for(i=0。i{bus[i]. = bus[i]. + bus[i].。//發(fā)電機(jī)功率=注入功率+負(fù)荷功率bus[i].= bus[i].+ bus[i].。bus[i].V=Complex::Rec2Polar(bus[i].V)。}coutfor(i=0。i{coutComplex::PrintfComplex(bus[i].V)。coutt(bus[i].)。coutt(bus[i].)。cout}coutfor(i=0。i{int statemp,endtemp。statemp=line[i].start。endtemp=line[i].end。coutComplex::PrintfComplex(Complex::ComDivRea(line[i].stoe,))。Complex::PrintfComplex(Complex::ComDivRea(line[i].etos,))。cout} }return 0。}include using namespace std。class Complex//定義復(fù)數(shù)類 { public: double real。double image。intRecPolar。//0表示直角坐標(biāo)，1表示極坐標(biāo)static Complex CaddC(Complex c1,Complex c2)。//求兩個(gè)復(fù)數(shù)和static Complex subComplex(Complex c1,Complex c2)。//求兩個(gè)復(fù)數(shù)差static Complex productComplex(Complex c1,Complex c2)。//求兩個(gè)復(fù)數(shù)積static Complex divideComplex(Complex c1,Complex c2)。//求兩個(gè)復(fù)數(shù)商static Complex ComDivRea(Complex c1,double r2)。//除數(shù)static Complex getconj(Complex c1)。//求一個(gè)復(fù)數(shù)共軛static Complex getinverse(Complex c1)。//取倒數(shù)static double getComplexReal(Complex c1)。//求一個(gè)復(fù)數(shù)實(shí)部static double getCompleximage(Complex c1)。//求一個(gè)復(fù)數(shù)虛部static void PrintfComplex(Complex c1)。//顯示一個(gè)復(fù)數(shù)static void PrintfmultiComplex(Complex C,int N)。//顯示多個(gè)復(fù)數(shù)static void zeroComplex(Complex c1)。//將復(fù)數(shù)復(fù)零static Complex Rec2Polar(Complex c1)。//取極坐標(biāo)Complex(){RecPolar=0。} }。Complex Complex::Rec2Polar(Complex c1)//極坐標(biāo)表示 { Complex Node。=sqrt(*+*)。=atan2(,)*180/。=1。return Node。}Complex Complex::CaddC(Complex c1,Complex c2)//復(fù)數(shù)加法 {Complex Node。=+。=+。return Node。}Complex Complex::subComplex(Complex c1,Complex c2)//復(fù)數(shù)減法 {Complex Node。=。=。return Node。}Complex Complex::productComplex(Complex c1,Complex c2)//復(fù)數(shù)乘法 {Complex Node。=**。=*+*。return Node。} Complex Complex::divideComplex(Complex c1,Complex c2)//復(fù)數(shù)除法 {Complex Node。=(*+*)/(pow(,2)+pow(,2))。=(**)/(pow(,2)+pow(,2))。return Node。} Complex Complex::ComDivRea(Complex c1,double r1)//復(fù)數(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)//取實(shí)部 {return 。}doubleComplex::getCompleximage(Complex c1)//取虛部 {return 。}voidComplex::PrintfComplex(Complex c1)//按直角坐標(biāo)輸出 { 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é)點(diǎn)功率計(jì)算值==== icon=1進(jìn)行第2次迭代 節(jié)點(diǎn)電壓修正量================= ====節(jié)點(diǎn)功率計(jì)算值==== icon=1進(jìn)行第3次迭代 節(jié)點(diǎn)電壓修正量================= ====節(jié)點(diǎn)功率計(jì)算值==== icon=0,迭代結(jié)束。====節(jié)點(diǎn)電壓===============發(fā)電機(jī)發(fā)出功率======節(jié)點(diǎn)10。節(jié)點(diǎn)2。節(jié)點(diǎn)3。0節(jié)點(diǎn)25 。00 0 節(jié)點(diǎn)4。0節(jié)點(diǎn)26 。0 00 節(jié)點(diǎn)5。節(jié) 點(diǎn)27 。0 0 節(jié)點(diǎn)6。0節(jié)點(diǎn)28 。00 0 節(jié)點(diǎn)7。0節(jié)點(diǎn)29 。00 節(jié)點(diǎn)8 節(jié)點(diǎn)9 0 節(jié)點(diǎn)10 0 節(jié)點(diǎn)11 節(jié)點(diǎn)12 0 節(jié)點(diǎn)13 節(jié)點(diǎn)14 0 節(jié)點(diǎn)15 0 節(jié)點(diǎn)16 0 節(jié)點(diǎn)17 0 節(jié)點(diǎn)18 0 節(jié)點(diǎn)19 0 節(jié)點(diǎn)20 0 節(jié)點(diǎn)21 0 節(jié)點(diǎn)22 0 節(jié)點(diǎn)23 0 節(jié)點(diǎn)24 0 。。。0 20節(jié) 點(diǎn)30 。00 0======線路傳輸功率========== 2to103to14to204to35to206to206to407to507to608to606to906to10011to9010to904to12013to12 14to1230to29 15to1228to8 28to6 請(qǐng)按任意鍵繼續(xù)... 20to19 20to10 22to21 26to25 27to25 28to27 高等電力系統(tǒng)分析 IEEE30節(jié)點(diǎn)潮流程序班級(jí)：電研114班姓名：王大偉學(xué)號(hào)：2201100151