【正文】 功潮流方程變成式（2）所示線性的直流潮流方程。n階矩陣，不包括平衡節(jié)點(diǎn)，其元素是1236。B(i,j)=239。3）根據(jù)注入功率的情況，形成矩陣Pi（除平衡節(jié)點(diǎn)）。2．了解和掌握PSDBPA電力系統(tǒng)分析程序穩(wěn)態(tài)分析方法（即潮流計(jì)算）。采用什么算法以及迭代的最大步數(shù)可以由用戶指定。打開(kāi)TextEdit應(yīng)用程序時(shí)先讀入該文件。如為可調(diào)的，則要附加填寫(xiě)R卡。分析報(bào)告列表：并根據(jù)用戶的要求（通過(guò)控制語(yǔ)句控制），輸出各種潮流分析報(bào)告。2．實(shí)驗(yàn)程序以及實(shí)驗(yàn)結(jié)果。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)。p=B1(i2)。max(abs(Dfi)ipt=t+1。jd=angle(v)*180/p。D(i)=0。p=p+1。DX1=DX。For i=1:2:x K=k+1。continue q=2*j1。J(p,q)=X2。X3=0。J(p,q)=X1。J(p,q)=X2。if ph=n q=q+1。cess dyn鋤ic Voltage stabil時(shí)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)安全分析的全過(guò)程動(dò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ā)電機(jī)勵(lì)磁系統(tǒng)數(shù)學(xué)模型及參 數(shù)對(duì)電網(wǎng)動(dò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時(shí)【J】．Relay，2007，35(2 1)：2227．[7] 方思立，朱方．快速勵(lì)磁系統(tǒng)對(duì)系統(tǒng)穩(wěn)定的影響[J】．中 國(guó)電機(jī)工程學(xué)報(bào)，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] 劉取．電力系統(tǒng)穩(wěn)定性及發(fā)電機(jī)勵(lì)磁控制[M】．北京： 中國(guó)電力出版社，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] 陳利芳，陳天祿．淺談自并勵(lì)勵(lì)磁系統(tǒng)在大容量機(jī)組 中的應(yī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)=[]。end m=p+1。q=q+1。q=2*j1。X2=C(i)+G(i,j)*e(i)+B(i,j)*f(i)。q=q+1。switch B2(i4)X4=0。End k=0。If ph=n Df(noJ=[]。IfB2(i4)=2 p=2*i1。printcut(result1,S,b,c,result2)。B2(ph,1)=S(ph)。t=0。p=B1(i1)。Y=zeros（n）；Y(p,q)=Y(p,q)1./(B1(i3)*B1(i5)。在PSDBPA單線圖程序上，建立IEEE9節(jié)點(diǎn)系統(tǒng)單線圖。5）迭代過(guò)程。三繞組變壓器先按常規(guī)方法化為三臺(tái)兩繞組變壓器后再用此卡模擬。 定義潮流程序卡片格式文件，用戶可更改及調(diào)整該文件。清晰認(rèn)知電力系統(tǒng)分析的意義。實(shí)驗(yàn)數(shù)據(jù)如圖所示三母線電力系統(tǒng)中，支路電抗和節(jié)點(diǎn)注入的攻入如圖所示，編寫(xiě)程序，求個(gè)各條支路的有功潮流分布。三、直流潮流算法計(jì)算步驟1）選擇平衡節(jié)點(diǎn)。237。i,j185。相當(dāng)于注入的有功功率。在一些應(yīng)用場(chǎng)所，如輸電網(wǎng)絡(luò)中，只要考慮的是電力系統(tǒng)中有功功率的分布，而不需要計(jì)算各個(gè)節(jié)點(diǎn)電壓幅值，且要計(jì)算速度要快，這勢(shì)必要對(duì)潮流計(jì)算進(jìn)行簡(jiǎn)化處理，本節(jié)實(shí)驗(yàn)就是研究直流潮流計(jì)算，編程與調(diào)試，獲得電力系統(tǒng)中各支路的有功分布，為進(jìn)一步進(jìn)行電力系統(tǒng)分析作準(zhǔn)備。利用x(*)=x(k+1)+DX(k+1)進(jìn)行多次迭代，通過(guò)迭代判據(jù)得到所需要的精度值即準(zhǔn)確值x(*)。endk=k+1。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)。endx2=[oP,oQ]39。M(i,j)=0。% 賦初值，U為節(jié)點(diǎn)電壓的幅值，a為節(jié)點(diǎn)電壓的相位角U=ones(1,5)。B=imag(Y)。y(2,3)=1/(+)。2 基本步驟和設(shè)計(jì)流程圖形成了雅克比矩陣并建立了修正方程式，運(yùn)用牛頓拉夫遜法計(jì)算潮流的核心問(wèn)題已經(jīng)解決，已有可能列出基本計(jì)算步驟并編制流程圖。同時(shí)基于MATLAB的計(jì)算機(jī)算法以雙精度類型進(jìn)行數(shù)據(jù)的存儲(chǔ)和運(yùn)算, 數(shù)據(jù)精確度高，能進(jìn)行潮流計(jì)算中的各種矩陣運(yùn)算，使得傳統(tǒng)潮流計(jì)算方法更加優(yōu)化。平衡節(jié)點(diǎn)（松弛節(jié)點(diǎn)、參考節(jié)點(diǎn)、基準(zhǔn)節(jié)點(diǎn)、緩沖節(jié)點(diǎn)）：給定Vi、di為0，求Pi、Qi，一般假設(shè)第n個(gè)節(jié)點(diǎn)為平衡節(jié)點(diǎn)。)。sum1=sum1+c。end V2=A([3],[1,2])。a=op./V1。B=imag(Y)。end end end for m=1:n for j=1:n if m==j Y(m,j)=sum(y(m,:))+sum(yd(m,:))。yd(4,2)=*i。z(1,3)=*i。please enter the short value n:39。Vj(Gijcosqij+Bijsinqij)i=(1,2,3Ln)（15）DQi=QisVi229。MM247。231。247。=231。231。MM247。231。D247。=231。V1231。247。247。231。230。247。231。247。MM231。247。231。232。V2Dq2247。及B1n246。0247。M247。247。0231。0V0Vn248。(10)247。231。230。MM231。1(6)此外，與系統(tǒng)各節(jié)點(diǎn)無(wú)功功率相應(yīng)的導(dǎo)納BLi必定遠(yuǎn)遠(yuǎn)小于該節(jié)點(diǎn)自導(dǎo)納的虛部，即：BLi=Qi=Bii Vi2Gijsinqij=Bij因此，Qi=Vi2Bii(7)考慮到以上關(guān)系后，式(5)中系數(shù)矩陣中的元素表達(dá)式可以化簡(jiǎn)為：Hii=Vi2BiiHij=VVijBijLii=ViBiiLij=VVijBij2(8)這樣，式(5)中系數(shù)矩陣可以表示為：230。DQ二、實(shí)驗(yàn)器材：計(jì)算機(jī)、軟件（已安裝，包括各類編程軟件C語(yǔ)言、C++、VB、VC等、應(yīng)用軟件MATLAB等）、移動(dòng)存儲(chǔ)設(shè)備（學(xué)生自備，軟盤(pán)、U盤(pán)等）三、實(shí)驗(yàn)內(nèi)容：：PQ分解法潮流計(jì)算基本思想是：把節(jié)點(diǎn)功率表示為電壓向量的極坐標(biāo)方程式，抓住主要矛盾，以有功功率誤差作為修正電壓向量角度的依據(jù)，以無(wú)功功率誤差作為修正電壓幅值的依據(jù)，把有功功率和無(wú)功功率迭代分開(kāi)來(lái)進(jìn)行。233。234?；蛘归_(kāi)為：DP=HDd+NDV/VDQ=JDd+LDV/V(4)電力系統(tǒng)中有功功率主要與各節(jié)點(diǎn)電壓向量的角度有關(guān)，無(wú)功功率則主要受各節(jié)點(diǎn)電壓幅值的影響。247。VVB247。231。247。247。232。B11B12L230。2247。231。231。247。247。V1246。231。O231。0Vn247。248。231。Bnn248。247。247。01/Vn247。232。B11B12LP2231。231。Vn248。231。Bnn248。B11B12L231。231。n248。231。Bnn248。(2)根據(jù)(15)計(jì)算各節(jié)點(diǎn)有功功率誤差DPi，并求出DPi/Vi。z=zeros(n,n)。yd(1,2)=*i。y(j,m)=y(m,j)。end end end Y A=[,0。for k1=0:100 for m=1:(n1)sum=0。A([4],[1,2,3])=A([4],[1,2,3])+os。V2=V2+b。end A(1,4)=sum。實(shí)驗(yàn)結(jié)果：五、思考討論題或體會(huì)或?qū)Ω倪M(jìn)實(shí)驗(yàn)的建議