【正文】 參考文獻(xiàn)[1]陳珩．電力系統(tǒng)穩(wěn)定分析[M]．北京：中國電力出版社，2002：139—187．[2]鄭阿奇．MATLAB實(shí)用教程[M]．北京：電子工業(yè)出版社，2005：1243．[3] 束洪春，孫士云，等．云電送粵交商流混聯(lián)系統(tǒng)全過 程動(dòng)態(tài)電壓研究[J】．中國電力，2008，4l(10)：l4． SHU Hong—ch吼，SUN Shiyun，et a1．Research on fun prc39。謝謝include “” include include include include“” include“” include“” using namespace std。//輸入支路個(gè)數(shù)if((Lamp。//動(dòng)態(tài)分配支路結(jié)構(gòu)體Line::ScanfLineData(line)。ibus[i].=0。//輸出結(jié)點(diǎn)參數(shù)Complex **X。//計(jì)算節(jié)點(diǎn)導(dǎo)納矩陣Bus::PrintfNodeDnz(X,N)。for(i=0。LOOP:Bus::JisuanNodeI(X,bus,N)。//計(jì)算節(jié)點(diǎn)功率差值Bus::PrintfNodeScal(X,bus,N)。gauss::gauss_slove(JacAug,x,NN)。goto LOOP。=0。aa=Complex::productComplex(Complex::getconj(bus[statemp1].V), B)。line[i].stoe=Complex::productComplex(bus[statemp1].V,dd)。dd=Complex::CaddC(aa,cc)。i{bus[i]. = bus[i]. + bus[i].。i{coutComplex::PrintfComplex(bus[i].V)。i{int statemp,endtemp。Complex::PrintfComplex(Complex::ComDivRea(line[i].etos,))。//i作為整個(gè)程序的循環(huán)變量int N=Bus::ScanfBusNo()。N)==0){return 0。//輸出支路參數(shù)Bus *bus=new Bus[N]。}Bus::ScanfBusData(bus)。for(i=0。double **JacAug。x=new double[NN]。//計(jì)算節(jié)點(diǎn)功率Bus::JisuanNodeScal(X,bus,N)。//whether converbence看迭代是否結(jié)束if(icon==1){coutBus::JisuanJacAug(JacAug,X,bus,N)。//修正節(jié)點(diǎn)電壓// Bus::PrintfNodeV(bus,N)。i{int statemp,endtemp。statemp=line[i].start。cc=Complex::productComplex(bb , Complex::getconj(line[i].Y))。bb=Complex::subComplex(Complex::getconj(bus[endtemp1].V), Complex::getconj(bus[statemp1].V))。}coutBus::JisuanNodeScal(X,bus,N)。bus[i].V=Complex::Rec2Polar(bus[i].V)。coutt(bus[i].)。endtemp=line[i].end。}include using namespace std。//0表示直角坐標(biāo)，1表示極坐標(biāo)static Complex CaddC(Complex c1,Complex c2)。//求兩個(gè)復(fù)數(shù)商static Complex ComDivRea(Complex c1,double r2)。//求一個(gè)復(fù)數(shù)實(shí)部static double getCompleximage(Complex c1)。//將復(fù)數(shù)復(fù)零static Complex Rec2Polar(Complex c1)。=sqrt(*+*)。}Complex Complex::CaddC(Complex c1,Complex c2)//復(fù)數(shù)加法 {Complex Node。}Complex Complex::subComplex(Complex c1,Complex c2)//復(fù)數(shù)減法 {Complex Node。}Complex Complex::productComplex(Complex c1,Complex c2)//復(fù)數(shù)乘法 {Complex Node。} Complex Complex::divideComplex(Complex c1,Complex c2)//復(fù)數(shù)除法 {Complex Node。} Complex Complex::ComDivRea(Complex c1,double r1)//復(fù)數(shù)除數(shù) { Complex Node。} Complex Complex::getconj(Complex c1)//取共軛 {Complex Node。}Complex Complex::getinverse(Complex c1)//取倒數(shù) { Complex Node。return Node。(8)。(8)。=0。void gauss::gauss_slove(double **a,double *x,int NN){int n,i,j,k,*pivrow,**pivcol。pivcol=new int *[n]。//行pivcol[k][0]=k。j{if(pivot{pivot=fabs(a[i][j])。exit(0)。a[k][j]=al。a[i][k]=al。j{a[i][j]=aik*a[k][j]/a[k][k]。i){sum=0。} for(k=(n2)。x[pivcol[k][0]]=al。節(jié)點(diǎn)3。0 00 節(jié)點(diǎn)5。00 0 節(jié)點(diǎn)7。00 0======線路傳輸功率========== 2to103to14to204to35to206to206to407to507to608to606to906to10011to9010to904to12013to12 14to1230to29 15to1228to8 28to6 請(qǐng)按任意鍵繼續(xù)... 20to19 20to10 2。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 0 節(jié)點(diǎn)6。00 0 節(jié)點(diǎn)4。====節(jié)點(diǎn)電壓===============發(fā)電機(jī)發(fā)出功率======節(jié)點(diǎn)10。k){al=x[pivcol[k][1]]。j{sum +=a[i][j]*x[j]。//解方程for(i=(n2)。i{aik=a[i][k]。i{al=a[i][pivcol[k][1]]。j{al=a[pivrow[k]][j]。//行pivcol[k][1]=j。//列n*2矩陣for(i=k。ipivot= fabs(a[k][k])。n=NN。static void gauss_output()。cout”。cout”。}doubleComplex::getCompleximage(Complex c1)//取虛部 {return 。=0。=。=(r1)。=(**)/(pow(,2)+pow(,2))。=*+*。=。=+。=1。} }。//顯示一個(gè)復(fù)數(shù)static void PrintfmultiComplex(Complex C,int N)。//求一個(gè)復(fù)數(shù)共軛static Complex getinverse(Complex c1)。//求兩個(gè)復(fù)數(shù)差static Complex productComplex(Complex c1,Complex c2)。double image。Complex::PrintfComplex(Complex::ComDivRea(line[i].etos,))。i{int statemp,endtemp。i{coutComplex::PrintfComplex(bus[i].V)。i{bus[i]. = bus[i]. + bus[i].。dd=Complex::CaddC(aa,cc)。line[i].stoe=Complex::productComplex(bus[statemp1].V,dd)。aa=Complex::productComplex(Complex::getconj(bus[statemp1].V), B)。=0。goto LOOP。gauss::gauss_slove(JacAug,x,NN)。//計(jì)算節(jié)點(diǎn)功率差值 Bus::PrintfNodeScal(X,bus,N)。LOOP:Bus::JisuanNodeI(X,bus,N)。for(i=0。//計(jì)算節(jié)點(diǎn)導(dǎo)納矩陣Bus::PrintfNodeDnz(X,N)。//輸出結(jié)點(diǎn)參數(shù)Complex **X。ibus[i].=0。//動(dòng)態(tài)分配支路結(jié)構(gòu)體Line::ScanfLineData(line)。//輸入支路個(gè)數(shù)if((Lamp。}include “” include include include include“” include“” include“” using namespace std。endtemp=line[i].end。coutt(bus[i].)。bus[i].V=Complex::Rec2Polar(bus[i].V)。}coutBus::JisuanNodeScal(X,bus,N)。bb=Complex::subComplex(Complex::getconj(bus[endtemp1].V), Complex::getconj(bus[statemp1].V))。cc=Complex::productComplex(bb , Complex::getconj(line[i].Y))。statemp=line[i].start。i{int statemp,endtemp。//修正節(jié)點(diǎn)電壓// Bus::PrintfNodeV(bus,N)。//whether converbence看迭代是否結(jié)束if(icon==1){coutBus::JisuanJacAug(JacAug,X,bus,N)。//計(jì)算節(jié)點(diǎn)功率Bus::JisuanNodeScal(X,bus,N)。x=new double[NN]。double **JacAug。for(i=0。}Bus::ScanfBusData(bus)。//輸出支路參數(shù)Bus *bus=new Bus[N]。N)==0){return 0。//i作為整個(gè)程序的循環(huán)變量int N=Bus::ScanfBusNo()。2007，35(1)：8184．[11] 方思立，劉增煌，孟慶和．大型汽輪發(fā)電機(jī)自并勵(lì)勵(lì) 磁系統(tǒng)的應(yīng)用條件【J]．中國電力，1994，27(12)：61．63． FANG Si．Ii，LIU Zenghu鋤g，MENG Qin爭hc．m application conditions of large turbine generator selfexcitation system【J】．Electric Powef，1994，27(12)： 61．63．[12]梁小冰，黃方能．利用EMTDC進(jìn)行長持續(xù)時(shí)間過程 的仿真研究【J】．電網(wǎng)技術(shù)，2002，26(9)：55．57． LIANG Xiaobing，HUANG Fan爭眥ng．How to cany out simulalion of long dul‘a(chǎn)tion processes by use of EMTDC【J】．Power System 11echnology，2002，26(9)： 5557．[13]王卉，陳楷，彭哲，等．?dāng)?shù)字仿真技術(shù)在電力系統(tǒng)中 的應(yīng)用及常用的幾種數(shù)字仿真工具【J】．繼電器，2004，32(21)：7l一75．wANG Hui，CHEN Kai，PENG zhe，et a1．Application of digital simulation眥hniques棚d severaJ simulation tools in power system[J】．Relay，2004，32(21)：71J(:,no)=[]。J(no:)=[]。end J(p,q)=X1。X4= C(i)+G(i,j)*f(i)B(i,j)*e(i)。X2=C(i)+G(i,j)*e(i)+B(i,j)*f(i)。Else if j=amp。J(p,q)=X2。J(p,q)=X1。p=2*i1。X3=X2。X1=G(i)*f(