【正文】 自動化學院 電氣工程基于牛頓拉夫遜法的潮流計算例題編程報學院(系): 專業(yè): 題目: 任課教師 碩士導師 告楊偉 XX2015年6月10號基于牛頓拉夫遜法的潮流計算例題編程報告摘要：電力系統(tǒng)潮流計算的目的在于：確定電力系統(tǒng)的運行方式、檢查系統(tǒng)中各元件是否過壓或者過載、為電力系統(tǒng)繼電保護的整定提供依據(jù)、為電力系統(tǒng)的穩(wěn)定計算提供初值、為電力系統(tǒng)規(guī)劃和經(jīng)濟運行提供分析的基礎(chǔ)。例題如下：用牛頓拉夫遜法計算下圖所示系統(tǒng)的潮流分布，其中系統(tǒng)中5為平衡節(jié)點，節(jié)點5電壓保持U=，其他四個節(jié)點分別為PQ節(jié)點，給定的注入功率如圖所示。電力系統(tǒng)潮流計算，一般來說，各個母線所供負荷的功率是已知的，各個節(jié)點的電壓是未知的（平衡節(jié)點外）可以根據(jù)網(wǎng)絡(luò)結(jié)構(gòu)形成節(jié)點導納矩陣，然后由節(jié)點導納矩陣列寫功率方程，由于功率方程里功率是已知的，電壓的幅值和相角是未知的，這樣潮流計算的問題就轉(zhuǎn)化為求解非線性方程組的問題了。% 如圖所示1，2，3，4為PQ節(jié)點,5為平衡節(jié)點y=0。y(1,4)=1/(+)。y(3,4)=1/(+)。endendend% 求節(jié)點導納矩陣中自導納for i=1:5Y(i,i)=sum(y(i,:))。S(2)=。P=real(S)。a=zeros(1,5)。while max(x2)1e6for i=1:4for j=1:4H(i,j)=0。oP(i)=0。endoP(i)=oP(i)+P(i)。N(i,j)=U(i)*U(j)*(G(i,j)*cos(a(i)a(j))+B(i,j)*sin(a(i)a(j)))。N(i,i)=N(i,i)U(i)*U(j)*(G(i,j)*cos(a(i)a(j))+B(i,j)*sin(a(i)a(j)))。endJ=[H,N。endfor i=1:4a(i)=a(i)+oa(i)。for j=1:5P(i)=U(i)*U(j)*(G(i,j)*cos(a(i)a(j))+B(i,j)*sin(a(i)a(j)))+P(i)。四運行結(jié)果節(jié)點導納矩陣經(jīng)過五次迭代后的雅克比矩陣迭代次數(shù)以及節(jié)點電壓的幅值和相角（弧度數(shù)）節(jié)點注入功率和電流五 結(jié)果分析在這次學習和實際操作過程里：首先，對電力系統(tǒng)分析中潮流計算的部分特別是潮流計算的計算機算法中的牛頓拉夫遜法進行深入的研讀，弄明白了其原理、計算過程、公式推導以及設(shè)計流程。第三篇：電力系統(tǒng)潮流計算程序設(shè)計電力系統(tǒng)潮流計算程序設(shè)計姓名：韋應順學號：2011021052 電力工程學院牛頓—拉夫遜潮流計算方法具有能夠?qū)⒎蔷€性方程線性化的特點，而使用MATLAB語言是由于MATLAB語言的數(shù)學邏輯強，易編譯。Y=zeros（n）；Y(p,q)=Y(p,q)1./(B1(i3)*B1(i5)。p=A(i1)；Y(p p)=1./A(i2)。p=B1(i1)。Q=imag(S)。t=0。[C,D,DF]=xxf(G,B,e,f,P,Q,n,B2,ph,V,no)；J=jacci(Y,G,B,P,Q,e,f,V,C,D,B2,n,ph,no)；[De,Df]=hxf(J,Df,ph,n,no)。B2(ph,1)=S(ph)。for i=1:ma(i)=conj((v(B1(i1))/B1(i5)v(B1(i2))/B1(i3))。printcut(result1,S,b,c,result2)。D(i)=D(i)+G(i,j)*f(j)+B(i,j)*e(j)。IfB2(i4)=2 p=2*i1。Df(p)=P(i)P1。If ph=n Df(noJ=[]。if ph=n DX(no)=0。End k=0。De(k)=DX(i+1)。switch B2(i4)X4=0。for j=1:n m=p+1。q=q+1。J(m,q)=X4。X2=C(i)+G(i,j)*e(i)+B(i,j)*f(i)。X4=0。q=2*j1。j=jph m=p+1。q=q+1。J(m,q)=X4。end m=p+1。J(p,q)=X2。End J(:,no)=[]。75．[14]IEEE Power Engmeering Socie哆．IEEE std 421．5．2005 IEEE玎ccOmmended practice for excitation system models for power system stabiI時studies【s】．第四篇：電力系統(tǒng)潮流計算程序電力系統(tǒng)潮流計算c語言程序，兩行，大家可以看看，仔細研究，然后在這個基礎(chǔ)上修改。//輸入節(jié)點個數(shù)int L=Line::ScanflineNo()。} //如果找不到兩個文件中的任意一個，退出Line *line=new Line[L]。//動態(tài)分配結(jié)點結(jié)構(gòu)體for(int i=0。//輸入節(jié)點參數(shù)Bus::PrintfBusData(bus,N)。iBus::JisuanNodeDnz(X,line,bus,L,N)。JacAug=new double *[NN]。int count=1。//計算節(jié)點功率Bus::JisuanNodeSdelta(bus,N)。//計算雅可比增廣矩陣 // Bus::PrintfJacAug(JacAug,N)。count++。Complex aa,bb,cc,dd,B。endtemp=line[i].end。dd=Complex::CaddC(aa,cc)。cc=Complex::productComplex(bb , Complex::getconj(line[i].Y))。//計算節(jié)點功率for(i=0。}coutfor(i=0。cout}coutfor(i=0。coutComplex::PrintfComplex(Complex::ComDivRea(line[i].stoe,))。int _tmain(int argc, _TCHAR* argv[]){int i。amp。//輸入支路參數(shù)Line::PrintfLineData(line,L)。bus[i].=0。X=new Complex *[N]。//輸出節(jié)點導納矩陣int NN=(N1)*2。idouble *x。//計算節(jié)點注入電流Bus::JisuanNodeScal(X,bus,N)。//輸出節(jié)點功率差值int icon=wehcon1(bus,N)。//解方程組求出電壓差值Bus::ReviseNodeV(bus,x,N)。}else{for(i=0。=line[i].B。bb=Complex::subComplex(Complex::getconj(bus[statemp1].V), Complex::getconj(bus[endtemp1].V))。aa=Complex::productComplex(Complex::getconj(bus[endtemp1].V), B)。line[i].etos=Complex::productComplex(bus[endtemp1].V,dd)。//發(fā)電機功率=注入功率+負荷功率bus[i].= bus[i].+ bus[i].。coutt(bus[i].)。statemp=line[i].start。cout} }return 0。intRecPolar。//求兩個復數(shù)積static Complex divideComplex(Complex c1,Complex c2)。//取倒數(shù)static double getComplexReal(Complex c1)。//顯示多個復數(shù)static void zeroComplex(Complex c1)。Complex Complex::Rec2Polar(Complex c1)//極坐標表示 { Complex Node。return Node。return Node。return Node。return Node。return Node。return Node。return Node。Node=(Complex::divideComplex(Node,c1))。}voidComplex::PrintfComplex(Complex c1)//按直角坐標輸出 { if(==0){ (6)。(6)。} else {coutComplex::zeroComplex(Complex c1)//清零 { =0。}。pivrow=new int[n]。pivrow[k]=k。i{for(j=k。//列}}}if(pivot{coutgetchar()。a[pivrow[k]][j]=a[k][j]。a[i][pivcol[k][1]]=a[i][k]。for(j=k。i=0。 }x[i]=(a[i][n]sum)/a[i][i]。x[pivcol[k][1]]=x[pivcol[k][0]]。節(jié)點2。0節(jié)點26 。0節(jié)點28 。0 20節(jié) 點30 。A、等值電路比 型等值電路更精確；B、采用 等值電路可以減少電力系統(tǒng)等值電路的節(jié)點總數(shù)；C、采用 等值電路可以增加電力系統(tǒng)的節(jié)點總數(shù)；D、電力系統(tǒng)運行方式改變時，采用 等值電路更方便節(jié)點導納矩陣的修改。A、幾何平均距離越大，電納越大；B、幾何平均距離越大，電納越??；C、輸電線路的電納與幾何平均距離無關(guān)；D、改變導線之間的幾何平均距離可以明顯改變線路的電納。A、減小線路電抗； B、增大線路電納；C、減小線路電阻； D、改善輸電線路的電暈條件。對于自耦變壓器，等值電路中各繞組的電阻，下述說法中正確的是（）。1電力系統(tǒng)等值電路中，所有參數(shù)應為歸算到同一電壓等級（基本級）的參數(shù)，關(guān)于基本級的選擇，下述說法中正確的是（）。1關(guān)于電力系統(tǒng)等值電路參數(shù)計算時，變壓器變比的選擇，下述說法中正確的是（）。A、變壓器 的功率大于變壓器 的功率；B、變壓器 的功率小于變壓器 的功率；C、變壓器 和變壓器 的功率相等；D、不能確定，還需其他條件。1在下圖所示的簡單網(wǎng)絡(luò)中，變壓器T中（）。1環(huán)形網(wǎng)絡(luò)中自然功率的分布規(guī)律是（） 在不計網(wǎng)絡(luò)功率損耗的情況下，下圖所示網(wǎng)絡(luò)各段電路中（）。A、網(wǎng)絡(luò)的潮流分布可以調(diào)控； B、網(wǎng)絡(luò)的潮流分布不可以調(diào)控；C、網(wǎng)絡(luò)的潮流分布由線路長度決定； D、網(wǎng)絡(luò)的潮流分布由線路阻抗確定。A、節(jié)點電壓大小和節(jié)點電壓相角；B、節(jié)點電壓大小和發(fā)電機無功功率；C、發(fā)電機有功功率和無功功率；D、節(jié)點電壓相角和發(fā)電機無功功率。（）同步發(fā)電機進相運行時，其運行極限由發(fā)電機并列運行的穩(wěn)定性和端部發(fā)熱條件確定。（）分裂導線多采用2～4分裂，最多不超過6分裂。（）10對于容量比不等于100/100/100的三繞組自耦變壓器，計算變壓器參數(shù)時不僅需要對銘牌給出的短路損耗進行歸算，還需要對銘牌給出的短路電壓進行歸算。（）。（）1均一電網(wǎng)功率的經(jīng)濟分布與其功率的自然分布相同。1如果兩個節(jié)點之間無直接聯(lián)系，則兩個節(jié)點互導納為零，兩個節(jié)點的互阻抗也為零。（）2電力系統(tǒng)潮流計算中，必須設(shè)置，并且只設(shè)置一個平衡節(jié)點。