【正文】 s[statemp1].V), Complex::getconj(bus[endtemp1].V))。endtemp=line[i].end。=line[i].B。Complex aa,bb,cc,dd,B。}else{for(i=0。count++。//解方程組求出電壓差值Bus::ReviseNodeV(bus,x,N)。//計算雅可比增廣矩陣 // Bus::PrintfJacAug(JacAug,N)。//輸出節(jié)點功率差值int icon=wehcon1(bus,N)。//計算節(jié)點功率Bus::JisuanNodeSdelta(bus,N)。//計算節(jié)點注入電流Bus::JisuanNodeScal(X,bus,N)。int count=1。idouble *x。JacAug=new double *[NN]。//輸出節(jié)點導納矩陣int NN=(N1)*2。iBus::JisuanNodeDnz(X,line,bus,L,N)。X=new Complex *[N]。//輸入節(jié)點參數(shù)Bus::PrintfBusData(bus,N)。bus[i].=0。//動態(tài)分配結(jié)點結(jié)構(gòu)體for(int i=0。//輸入支路參數(shù)Line::PrintfLineData(line,L)。} //如果找不到兩個文件中的任意一個，退出Line *line=new Line[L]。amp。//輸入節(jié)點個數(shù)int L=Line::ScanflineNo()。int _tmain(int argc, _TCHAR* argv[]){int i。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ǔ)上修改。cess dyn鋤ic Voltage stabil時of hybrid AC／DC poWer tmnsmission System舶m Yu衄an proVince to G啪gdong province【J】．Electric Power，2008，4l(10)： l4．[4] 朱新立，湯涌，等．大電網(wǎng)安全分析的全過程動態(tài)仿 真技術(shù)[J】．電網(wǎng)技術(shù)，2008，32(22)：23—28． SONG Xin—Ii，TANG Yof唱，et a1． Full dyn鋤ic simulation for the stabilhy a眥lysis of large power system【J】．Power System融IlrIolo影，2008，32(22)： 23．28．[5]Roytelm鋤I，Shallidehpour S M．A prehsivc long teml dynaIIlic simulation for powcr system recoVery【J】． IEEE Transactions 0n Power Systems，1994，9(3)． [6] 石雩梅，汪志宏，等．發(fā)電機勵磁系統(tǒng)數(shù)學模型及參 數(shù)對電網(wǎng)動態(tài)穩(wěn)定性分析結(jié)果影響的研究[J】．繼電 器，2007，35(21)：2227．SHI Xue．mei，WANG Zllihon舀et a1．Iksearch on the innuence of g鋤e翰to璐baScd ∞de詛iled excitation system models柚d parameterS t0 power鏟id dyn鋤ic stabil時【J】．Relay，2007，35(2 1)：2227．[7] 方思立，朱方．快速勵磁系統(tǒng)對系統(tǒng)穩(wěn)定的影響[J】．中 國電機工程學報，1986，6(1)：20．28．FANG Si．1i，ZHU Fang．The effbct of f弧t．respon∞excitation system on the stability of power netwofk【J】． Proceedings ofthe CSEE，1986，6(1)：2028．[8] 劉?。娏ο到y(tǒng)穩(wěn)定性及發(fā)電機勵磁控制[M】．北京： 中國電力出版社，2007．LIU Qu．Power system S詛bility鋤d generator excitation control【M】．BeUing：ChiIla Electric Powef Press，2007． [9] Dallachy J L，Anderson T．EXperience with rcplacing ro詛ting exciters wim static exciters【J】．1k InStitution of Electrical Engineers，1 996．[10] 陳利芳，陳天祿．淺談自并勵勵磁系統(tǒng)在大容量機組 中的應用【J】．繼電器，2007，35(1)：8l培4． CHEN Lif抽島CHEN Tian—lIL Application of 辯l仁exci組tion mode in large capacity髫memtor unit【J】． ReIay39。End J(:,no)=[]。J(m,q)=X4。J(p,q)=X2。if ph=n q=q+1。end m=p+1。end q=2*j1。J(m,q)=X4。J(p,q)=X2。q=q+1。J(m,q)=X3。j=jph m=p+1。J(p,q)=X1。q=2*j1。P=2*i1。X4=0。X3=0。X2=C(i)+G(i,j)*e(i)+B(i,j)*f(i)。X1=D(i)+G(i,j)*f(i)B(i,j)*e(i)。J(m,q)=X4。J(p,q)=X2。q=q+1。J=ph J(m,q)=X3。for j=1:n m=p+1。continue q=2*j1。switch B2(i4)X4=0。X2=G(i,j)*e(i)+B(i,j)*f(i)。De(k)=DX(i+1)。For i=1:2:x K=k+1。End k=0。For i=(no+2):(x1+2)DX(i)=DX1(i2)。if ph=n DX(no)=0。DX1=DX。If ph=n Df(noJ=[]。Df(p)=Q(i)Q1。Df(p)=P(i)P1。p=p+1。IfB2(i4)=2 p=2*i1。Q1=C(i)*f(i)D(i)*e(i)。D(i)=D(i)+G(i,j)*f(j)+B(i,j)*e(j)。D(i)=0。printcut(result1,S,b,c,result2)。c(i)=v(B1(i2))*a(i)j*B1(i4)*v(B1(i2))^2/2。for i=1:ma(i)=conj((v(B1(i1))/B1(i5)v(B1(i2))/B1(i3))。jd=angle(v)*180/p。B2(ph,1)=S(ph)。for i=1:n hh(i)=conj(Y(ph,i)*v(i))。[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)。max(abs(Dfi)ipt=t+1。t=0。J=jacci(Y,G,B,P,Q,e,f,V,C,D,B2,n,ph,no)。Q=imag(S)。p=B1(i2)。p=B1(i1)。End f(i)=imag(B2(i3))。p=A(i1)；Y(p p)=1./A(i2)。f=zeros（1，n）；no=2*ph=1； Y(q,q)=Y(q,q)+1./B1(i3)+B1(i4)/2；End for i=1：nG=real(Y)。Y=zeros（n）；Y(p,q)=Y(p,q)1./(B1(i3)*B1(i5)。ph=input(‘n請輸入平衡母線的節(jié)點號：ph=’)； B1=input（‘n請輸入支路信號：B1=’）。第三篇：電力系統(tǒng)潮流計算程序設(shè)計電力系統(tǒng)潮流計算程序設(shè)計姓名：韋應順學號：2011021052 電力工程學院牛頓—拉夫遜潮流計算方法具有能夠?qū)⒎蔷€性方程線性化的特點，而使用MATLAB語言是由于MATLAB語言的數(shù)學邏輯強，易編譯。利用x(*)=x(k+1)+DX(k+1)進行多次迭代，通過迭代判據(jù)得到所需要的精度值即準確值x(*)。四運行結(jié)果節(jié)點導納矩陣經(jīng)過五次迭代后的雅克比矩陣迭代次數(shù)以及節(jié)點電壓的幅值和相角（弧度數(shù)）節(jié)點注入功率和電流五 結(jié)果分析在這次學習和實際操作過程里：首先，對電力系統(tǒng)分析中潮流計算的部分特別是潮流計算的計算機算法中的牛頓拉夫遜法進行深入的研讀，弄明白了其原理、計算過程、公式推導以及設(shè)計流程。endS(5)=P(5)+Q(5)*sqrt(1)。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)。endk=k+1。endfor i=1:4a(i)=a(i)+oa(i)。% x1為所求△x的列向量% 求節(jié)點電壓新值，準備下一次迭代for i=1:4oa(i)=x1(i)。endJ=[H,N。L(i,i)=L(i,i)U(i)*U(j)*(G(i,j)*sin(a(i)a(j))B(i,j)*cos(a(i)a(j)))endendN(i,i)=N(i,i)2*(U(i))^2*G(i,i)。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)))。M(i,j)=N(i,j)。N(i,j)=U(i)*U(j)*(G(i,j)*cos(a(i)a(j))+B(i,j)*sin(a(i)a(j)))。end