【正文】 由于電力系統(tǒng)規(guī)模的不斷擴(kuò)大，對計(jì)算速度的要求不斷提高，計(jì)算機(jī)的并行計(jì)算技術(shù)也將在潮流計(jì)算中得到廣泛的應(yīng)用，成為重要的研究領(lǐng)域。此外，隨著人工智能理論的發(fā)展，遺傳算法、人工神經(jīng)網(wǎng)絡(luò)、模糊算法也逐漸被引入潮流計(jì)算。另外，為了解決病態(tài)潮流計(jì)算，出現(xiàn)了將潮流計(jì)算表示為一個(gè)無約束非線性規(guī)劃問題的模型，即非線性規(guī)劃潮流算法。牛頓法的特點(diǎn)是將非線性方程線性化。在牛頓法的基礎(chǔ)上，根據(jù)電力系統(tǒng)的特點(diǎn)，抓住主要矛盾，對純數(shù)學(xué)的牛頓法進(jìn)行了改造，得到了PQ分解法。解決電力系統(tǒng)潮流計(jì)算問題是以導(dǎo)納矩陣為基礎(chǔ)的，因此，只要在迭代過程中盡可能保持方程式系數(shù)矩陣的稀疏性，就可以大大提高牛頓潮流程序的計(jì)算效率?？朔杩狗ㄈ秉c(diǎn)的另一途徑是采用牛頓拉夫遜法（以下簡稱牛頓法）。為了克服阻抗法在內(nèi)存和速度方面的缺點(diǎn)，后來發(fā)展了以阻抗矩陣為基礎(chǔ)的分塊阻抗法。但是，阻抗法的主要缺點(diǎn)就是占用計(jì)算機(jī)的內(nèi)存很大，每次迭代的計(jì)算量很大。而且阻抗法每迭代一次都要求順次取阻抗矩陣中的每一個(gè)元素進(jìn)行計(jì)算，因此，每次迭代的計(jì)算量很大。阻抗矩陣是滿矩陣，阻抗法要求計(jì)算機(jī)儲(chǔ)存表征系統(tǒng)接線和參數(shù)的阻抗矩陣。這個(gè)方法的原理比較簡單，要求的數(shù)字計(jì)算機(jī)的內(nèi)存量也比較小，適應(yīng)當(dāng)時(shí)的電子數(shù)字計(jì)算機(jī)制作水平和電力系統(tǒng)理論水平，于是電力系統(tǒng)計(jì)算人員轉(zhuǎn)向以阻抗矩陣為主的逐次代入法（以下簡稱阻抗法）。這種情況促使電力系統(tǒng)的研究人員不斷尋求新的更可靠的計(jì)算方法。非線性代數(shù)方程組的解法離不開迭代，因此，潮流計(jì)算方法首先要求它是能可靠的收斂，并給出正確答案。00 0======線路傳輸功率========== 2to103to14to204to35to206to206to407to507to608to606to906to10011to9010to904to12013to12 14to1230to29 15to1228to8 28to6 請按任意鍵繼續(xù)... 20to19 20to10 22to21 26to25 27to25 28to27 高等電力系統(tǒng)分析 IEEE30節(jié)點(diǎn)潮流程序班級：電研114班姓名：王大偉學(xué)號：2201100151第五篇：電力系統(tǒng)潮流計(jì)算發(fā)展史電力系統(tǒng)潮流計(jì)算發(fā)展史對潮流計(jì)算的要求可以歸納為下面幾點(diǎn)：（1）算法的可靠性或收斂性（2）計(jì)算速度和內(nèi)存占用量（3）計(jì)算的方便性和靈活性電力系統(tǒng)潮流計(jì)算屬于穩(wěn)態(tài)分析范疇，不涉及系統(tǒng)元件的動(dòng)態(tài)特性和過渡過程。。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 。00 0 節(jié)點(diǎn)7。0 0 節(jié)點(diǎn)6。0 00 節(jié)點(diǎn)5。00 0 節(jié)點(diǎn)4。節(jié)點(diǎn)3。====節(jié)點(diǎn)電壓===============發(fā)電機(jī)發(fā)出功率======節(jié)點(diǎn)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)。//顯示一個(gè)復(fù)數(shù)static void PrintfmultiComplex(Complex C,int N)。//求一個(gè)復(fù)數(shù)實(shí)部static double getCompleximage(Complex c1)。//求一個(gè)復(fù)數(shù)共軛static Complex getinverse(Complex c1)。//求兩個(gè)復(fù)數(shù)商static Complex ComDivRea(Complex c1,double r2)。//求兩個(gè)復(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é)點(diǎn)電壓// Bus::PrintfNodeV(bus,N)。gauss::gauss_slove(JacAug,x,NN)。//whether converbence看迭代是否結(jié)束if(icon==1){coutBus::JisuanJacAug(JacAug,X,bus,N)。//計(jì)算節(jié)點(diǎn)功率差值 Bus::PrintfNodeScal(X,bus,N)。//計(jì)算節(jié)點(diǎn)功率Bus::JisuanNodeScal(X,bus,N)。LOOP:Bus::JisuanNodeI(X,bus,N)。x=new double[NN]。for(i=0。double **JacAug。//計(jì)算節(jié)點(diǎn)導(dǎo)納矩陣Bus::PrintfNodeDnz(X,N)。for(i=0。//輸出結(jié)點(diǎn)參數(shù)Complex **X。}Bus::ScanfBusData(bus)。ibus[i].=0。//輸出支路參數(shù)Bus *bus=new Bus[N]。//動(dòng)態(tài)分配支路結(jié)構(gòu)體Line::ScanfLineData(line)。N)==0){return 0。//輸入支路個(gè)數(shù)if((Lamp。//i作為整個(gè)程序的循環(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(Complex::getconj(bus[endtemp1].V), Complex::getconj(bus[statemp1].V))。line[i].stoe=Complex::productCo