【正文】 ,3)=1/(+)。2 基本步驟和設計流程圖形成了雅克比矩陣并建立了修正方程式，運用牛頓拉夫遜法計算潮流的核心問題已經(jīng)解決，已有可能列出基本計算步驟并編制流程圖。同時基于MATLAB的計算機算法以雙精度類型進行數(shù)據(jù)的存儲和運算, 數(shù)據(jù)精確度高，能進行潮流計算中的各種矩陣運算，使得傳統(tǒng)潮流計算方法更加優(yōu)化。20世紀70年代后期，有人提出采用更精確的模型，即將泰勒級數(shù)的高階項也包括進來，希望以此提高算法的性能，這便產(chǎn)生了保留非線性的潮流算法。這個方法把一個大系統(tǒng)分割為幾個小的地區(qū)系統(tǒng)，在計算機內(nèi)只需存儲各個地區(qū)系統(tǒng)的阻抗矩陣及它們之間的聯(lián)絡線的阻抗，這樣不僅大幅度的節(jié)省了內(nèi)存容量，同時也提高了節(jié)省速度。20世紀60年代初，數(shù)字計算機已經(jīng)發(fā)展到第二代，計算機的內(nèi)存和計算速度發(fā)生了很大的飛躍，從而為阻抗法的采用創(chuàng)造了條件。第一篇：電力系統(tǒng)潮流計算發(fā)展史電力系統(tǒng)潮流計算發(fā)展史對潮流計算的要求可以歸納為下面幾點：（1）算法的可靠性或收斂性（2）計算速度和內(nèi)存占用量（3）計算的方便性和靈活性電力系統(tǒng)潮流計算屬于穩(wěn)態(tài)分析范疇，不涉及系統(tǒng)元件的動態(tài)特性和過渡過程。阻抗矩陣是滿矩陣，阻抗法要求計算機儲存表征系統(tǒng)接線和參數(shù)的阻抗矩陣。克服阻抗法缺點的另一途徑是采用牛頓拉夫遜法（以下簡稱牛頓法）。另外，為了解決病態(tài)潮流計算，出現(xiàn)了將潮流計算表示為一個無約束非線性規(guī)劃問題的模型，即非線性規(guī)劃潮流算法。一 研究內(nèi)容通過一道例題來認真分析牛頓拉夫遜法的原理和方法（采用極坐標形式的牛拉法），同時掌握潮流計算計算機算法的相關知識，能看懂并初步使用MATLAB軟件進行編程，培養(yǎng)自己電力系統(tǒng)潮流計算機算法編程能力。由課本總結基本步驟如下：1)形成節(jié)點導納矩陣Y；2)設各節(jié)點電壓的初值，如果是直角坐標的話設電壓的實部e和虛部f；如果是極坐標的話則設電壓的幅值U和相角a；3)將各個節(jié)點電壓的初值代入公式求修正方程中的不平衡量以及修正方程的系數(shù)矩陣的雅克比矩陣；4)解修正方程式，求各節(jié)點電壓的變化量，即修正量； 5)計算各個節(jié)點電壓的新值，即修正后的值；6)利用新值從第（3）步開始進入下一次迭代，直至達到精度退出循環(huán)； 7)計算平衡節(jié)點的功率和線路功率，輸出最后計算結果； ① 公式推導② 流程圖三matlab編程代碼clear。y(2,5)=1/(+)。% 輸入原始節(jié)點的給定注入功率S(1)=+。U(5)=。L(i,j)=0。% x2為不平衡量列向量% 求雅克比矩陣% 當i~=j時，求H,N,M,Lfor i=1:4for j=1:4if i~=jH(i,j)=U(i)*U(j)*(G(i,j)*sin(a(i)a(j))B(i,j)*cos(a(i)a(j)))。L(i,i)=L(i,i)+2*(U(i))^2*B(i,i)。endk,U,a% 求節(jié)點注入功率i=5。六 結論通過這個任務，自己在matlab編程，潮流計算，word文檔的編輯功能等方面均有提高，但也暴漏出一些問題：理論知識儲備不足，對matlab的性能和特點還不能有一個全面的把握，對word軟件也不是很熟練，相信通過以后的學習能彌補這些不足，達到一個新的層次。if A（i2）=0B=imag(Y)。end Y(p,p)=Y(p,p)+1./(B1(i3)*B1(i5)^2+B1(i4)./2P=real(S)。e=e+De；f=f+Df。resulte1=[A(:,1),real(v),imag(v),V,jd,real(S’)，imag(S’),real(B2(:1)),imag(B2(:1)),real(B2(:2)),imag(B2(:,2))]。For j=i:nC(i)=C(i)+G(i,j)*e(j)B(i,j)*f(j)。else p=2*i1。x1=length(DX1)。Df(k)=DX(i)。case 1 J(p,q)=X1。X3=X2。J(p,q)=X1。Else if j=amp。X4= C(i)+G(i,j)*f(i)B(i,j)*e(i)。J(no:)=[]。2007，35(1)：8184．[11] 方思立，劉增煌，孟慶和．大型汽輪發(fā)電機自并勵勵 磁系統(tǒ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進行長持續(xù)時間過程 的仿真研究【J】．電網(wǎng)技術，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]王卉，陳楷，彭哲，等．數(shù)字仿真技術在電力系統(tǒ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)：71N)==0){return 0。}Bus::ScanfBusData(bus)。double **JacAug。//計算節(jié)點功率Bus::JisuanNodeScal(X,bus,N)。//修正節(jié)點電壓// Bus::PrintfNodeV(bus,N)。statemp=line[i].start。bb=Complex::subComplex(Complex::getconj(bus[endtemp1].V), Complex::getconj(bus[statemp1].V))。bus[i].V=Complex::Rec2Polar(bus[i].V)。endtemp=line[i].end。//輸入支路個數(shù)if((Lamp。ibus[i].=0。//計算節(jié)點導納矩陣Bus::PrintfNodeDnz(X,N)。LOOP:Bus::JisuanNodeI(X,bus,N)。gauss::gauss_slove(JacAug,x,NN)。=0。line[i].stoe=Complex::productComplex(bus[statemp1].V,dd)。i{bus[i]. = bus[i]. + bus[i].。i{int statemp,endtemp。double image。//求一個復數(shù)共軛static Complex getinverse(Complex c1)。} }。=+。=*+*。=(r1)。=0。cout”。static void gauss_output()。ipivot= fabs(a[k][k])。//行pivcol[k][1]=j。i{al=a[i][pivcol[k][1]]。//解方程for(i=(n2)。k){al=x[pivcol[k][1]]。00 0 節(jié)點4。00 節(jié)點8 節(jié)點9 0 節(jié)點10 0 節(jié)點11 節(jié)點12 0 節(jié)點13 節(jié)點14 0 節(jié)點15 0 節(jié)點16 0 節(jié)點17 0 節(jié)點18 0 節(jié)點19 0 節(jié)點20 0 節(jié)點21 0 節(jié)點22 0 節(jié)點23 0 節(jié)點24 0 。架空輸電線路的電納和導線之間幾何平均距離的關系為（）。A、采用 型等值電路比采用T型等值電路精確；B、采用 型等值電路在變壓器變比改變時，便于電力系統(tǒng)節(jié)點導納矩陣的修改； C、采用采用 型等值電路可以減少電力系統(tǒng)的節(jié)點總數(shù)；D、采用采用 型等值電路可以增加電力系統(tǒng)的節(jié)點總數(shù)。A、選擇功率和電壓基準值； B、選擇功率和電流基準值；C、選擇功率和阻抗基準值； D、選擇電壓和阻抗基準值。D、有功功率和無功功率都從末端流向首端。2對于下圖所示的放射性網(wǎng)絡，下述說法中正確的是（）。二、判斷題（下述說法中，對于你認為正確的請選擇“Y”，錯誤的選擇“N”，每小題2分，共50分）同步發(fā)電機降低功率因數(shù)運行時，其運行極限由額定勵磁電流確定。（）對于容量比不等于100/100/100的普通三繞組變壓器，計算變壓器參數(shù)時需要對銘牌給出的短路損耗進行歸算，但銘牌給出的短路電壓不需歸算。（）1任何多電壓等級環(huán)網(wǎng)中都存在循環(huán)功率。（）2未裝設無功補償裝置的變電所母線為PQ節(jié)點。（）2高斯塞德爾潮流計算法，由于收斂速度慢，在電力系統(tǒng)潮流計算中很少單獨使用。（）1在環(huán)形電力網(wǎng)中串聯(lián)縱向串聯(lián)加壓器主要改變電網(wǎng)的有功功率分布。（）1同一電壓等級電力系統(tǒng)中，所有設備的額定電壓都相同。（）電力系統(tǒng)穩(wěn)態(tài)分析時，對于與無限大電力系統(tǒng)并列運行的定出力發(fā)電機，其數(shù)學模型為，約束條件為。2電力系統(tǒng)潮流計算采用的數(shù)學模型是（）。A、有功率通過； B、無功率通過； C、不能確定； D、僅有有功功率通過。A、精確計算時采用實際變比，近似計算時采用平均額定變比；B、近似計算時，采用實際變比；精確計算時采用平均額定變比C、不管是精確計算還是近似計算均應采用額定變比；D、不管是精確計算還是近似計算均應采用平均額定變比。A、等值電路中各繞組的電阻，是各繞組實際電阻按照變壓器變比歸算到同一電壓等級的電阻值；B、等值電路中各繞組的電阻就是各繞組的實際電阻；C、等值電路中各繞組的電阻是各繞組的等效電阻歸算到同一電壓等級的電阻值；D、等值電路中各繞組的電阻一定為正值，因為繞組總有電阻存在。在輸電線路參數(shù)中屬于耗能參數(shù)的是（）。00 0======線路傳輸功率========== 2to103to14to204to35to206to206to407to507to608to606to906to10011to9010to904to12013to12 14to1230to29 15to1228to8 28to6 請按任意鍵繼續(xù)... 20to19 20to10 22to21 26to25 27to25 28to27 高等電力系統(tǒng)分析 IEEE30節(jié)點潮流程序班級：電研114班姓名：王大偉學號：2201100151第五篇：電力系統(tǒng)潮流計算自測題（二）電力系統(tǒng)潮流計算與