【正文】 75．[14]IEEE Power Engmeering Socie哆．IEEE std 421．5．2005 IEEE玎ccOmmended practice for excitation system models for power system stabiI時studies【s】．。End J(:,no)=[]。J(p,q)=X2。end m=p+1。J(m,q)=X4。q=q+1。j=jph m=p+1。q=2*j1。X4=0。X2=C(i)+G(i,j)*e(i)+B(i,j)*f(i)。J(m,q)=X4。q=q+1。for j=1:n m=p+1。switch B2(i4)X4=0。De(k)=DX(i+1)。End k=0。if ph=n DX(no)=0。If ph=n Df(noJ=[]。Df(p)=P(i)P1。IfB2(i4)=2 p=2*i1。D(i)=D(i)+G(i,j)*f(j)+B(i,j)*e(j)。printcut(result1,S,b,c,result2)。for i=1:ma(i)=conj((v(B1(i1))/B1(i5)v(B1(i2))/B1(i3))。B2(ph,1)=S(ph)。[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)。t=0。Q=imag(S)。p=B1(i1)。p=A(i1)；Y(p p)=1./A(i2)。Y=zeros（n）；Y(p,q)=Y(p,q)1./(B1(i3)*B1(i5)。第五篇：電力系統(tǒng)潮流計算程序設(shè)計電力系統(tǒng)潮流計算程序設(shè)計姓名：韋應(yīng)順學(xué)號：2011021052 電力工程學(xué)院牛頓—拉夫遜潮流計算方法具有能夠?qū)⒎蔷€性方程線性化的特點，而使用MATLAB語言是由于MATLAB語言的數(shù)學(xué)邏輯強，易編譯。case2_539。3 5 0 0 0 0 0 0 1360 360。]。5 3 0 0 0 0 1 0 100 1 。%% Bus data% bus_i type PdQdGsBs area VmVabaseKV zone VmaxVmin =[ 1 1 160 80 0 0 1 1 0 100 1 。讀者利用MATPOWER進(jìn)行潮流計算時。六 結(jié)論通過這個任務(wù)，自己在matlab編程，潮流計算，word文檔的編輯功能等方面均有提高，但也暴漏出一些問題：理論知識儲備不足，對matlab的性能和特點還不能有一個全面的把握，對word軟件也不是很熟練，相信通過以后的學(xué)習(xí)能彌補這些不足，達(dá)到一個新的層次。S% 求節(jié)點注入電流I=Y*U39。endk,U,a% 求節(jié)點注入功率i=5。oU(i)=x1(i+4)*U(i)。L(i,i)=L(i,i)+2*(U(i))^2*B(i,i)。endendend% 當(dāng)i=j時，求H,N,M,Lfor i=1:4for j=1:5if i~=j H(i,i)=H(i,i)+U(i)*U(j)*(G(i,j)*sin(a(i)a(j))B(i,j)*cos(a(i)a(j)))。% x2為不平衡量列向量% 求雅克比矩陣% 當(dāng)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)))。oQ(i)=oQ(i)U(i)*U(j)*(G(i,j)*sin(a(i)a(j))B(i,j)*cos(a(i)a(j)))。L(i,j)=0。k=0。U(5)=。S(5)=0。% 輸入原始節(jié)點的給定注入功率S(1)=+。% 求節(jié)點導(dǎo)納矩陣中互導(dǎo)納for i=1:5for j=1:5if i~=jY(i,j)=y(i,j)。y(2,5)=1/(+)。y(1,3)=1/(+)。由課本總結(jié)基本步驟如下：1)形成節(jié)點導(dǎo)納矩陣Y；2)設(shè)各節(jié)點電壓的初值，如果是直角坐標(biāo)的話設(shè)電壓的實部e和虛部f；如果是極坐標(biāo)的話則設(shè)電壓的幅值U和相角a；3)將各個節(jié)點電壓的初值代入公式求修正方程中的不平衡量以及修正方程的系數(shù)矩陣的雅克比矩陣；4)解修正方程式，求各節(jié)點電壓的變化量，即修正量； 5)計算各個節(jié)點電壓的新值，即修正后的值；6)利用新值從第（3）步開始進(jìn)入下一次迭代，直至達(dá)到精度退出循環(huán)； 7)計算平衡節(jié)點的功率和線路功率，輸出最后計算結(jié)果； ① 公式推導(dǎo)② 流程圖三matlab編程代碼clear。牛頓迭代法其最大優(yōu)點是在方程f(x)=0的單根附近時誤差將呈平方減少，而且該法還可以用來求方程的重根、復(fù)根。一 研究內(nèi)容通過一道例題來認(rèn)真分析牛頓拉夫遜法的原理和方法（采用極坐標(biāo)形式的牛拉法），同時掌握潮流計算計算機算法的相關(guān)知識，能看懂并初步使用MATLAB軟件進(jìn)行編程，培養(yǎng)自己電力系統(tǒng)潮流計算機算法編程能力。這次課程設(shè)計對自己意義很大，自己從中獲得很多東西。原始數(shù)據(jù)輸入格式的設(shè)計，主要應(yīng)從使用的角度出發(fā)，原則是簡單明了，便于修改。通過M語言，可以用類似數(shù)學(xué)公式的方式來編寫算法，大大降低了程序所需的難度并節(jié)省了時間，從而可把主要的精力集中在算法的構(gòu)思而不是編程上。MATLAB是一種交互式、面向?qū)ο蟮某绦蛟O(shè)計語言，廣泛應(yīng)用于工業(yè)界與學(xué)術(shù)界，主要用于矩陣運算，同時在數(shù)值分析、自動控制模擬、數(shù)字信號處理、動態(tài)分析、繪圖等方面也具有強大的功能。MATLAB軟件的應(yīng)用MATLAB Compiler是一種編譯工具，它能夠?qū)編寫的函數(shù)文件生成函數(shù)庫或者可執(zhí)行文件COM組件等，以提供給其他高級語言如C++、C等進(jìn)行調(diào)用由此擴(kuò)展MATLAB的應(yīng)用范圍，將MATLAB的開發(fā)效率與其他高級語言的運行結(jié)合起來，取長補短，豐富程序開發(fā)的手段?？偨Y(jié)為在電力系統(tǒng)運行方式和規(guī)劃方案的研究中，都需要進(jìn)行潮流計算以比較運行方式或規(guī)劃供電方案的可行性、可靠性和經(jīng)濟(jì)性。具體表現(xiàn)在以下方面：①在電網(wǎng)規(guī)劃階段,通過潮流計算,合理規(guī)劃電源容量及接入點,合理規(guī)劃網(wǎng)架,選擇無功補償方案,滿足規(guī)劃水平的大、小方式下潮流交換控制、調(diào)峰、調(diào)相、調(diào)壓的要求。另外，論文還討論發(fā)現(xiàn)，當(dāng)變壓器支路阻抗過小時，利用Π型模型會產(chǎn)生數(shù)值巨大的對地導(dǎo)納，由此會導(dǎo)致潮流不收斂。待求的運行狀態(tài)參量包括電網(wǎng)各母線節(jié)點的電壓幅值和相角，以及各支路的功率分布、網(wǎng)絡(luò)的功率損耗等。四、結(jié)語對給定系統(tǒng)通過Matpower進(jìn)行潮流計算后，其運算結(jié)果從節(jié)點、支路等方面與Matlab編程與Simulink仿真的結(jié)果進(jìn)行了對比，得出系統(tǒng)運行穩(wěn)定，且在符合系統(tǒng)要求的情況下，三種方法的潮流計算結(jié)果基本一致。線路的有功損耗逐漸增加，四個變電所低壓側(cè)電壓均在允許范圍內(nèi)，符合課題要求，具體支路損耗見表5。（2）根據(jù)圖6可知，PV節(jié)點10的有功功率和電壓幅值在潮流計算過程中保持不變，是因為它的作用是平衡系統(tǒng)功率。（五）Matpower的M文件的編寫。并設(shè)平衡節(jié)點為節(jié)點1，PV節(jié)點為節(jié)點10，剩余節(jié)點為PQ節(jié)點。②變電所的負(fù)荷如表1所示：③每個變電所的功率因數(shù)均為cosΦ=；④變電所2和變電所4分別配有兩臺容量為75MVA的變壓器，短路損耗414KW，短路電壓（Uk%）=；變電所1和變電所3分別配有兩臺容量為63MVA的變壓器，短路損耗為245KW，短路電壓（Uk%）=；（3）輸電線路資料：發(fā)電廠和變電所之間的輸電線路的電壓等級及長度標(biāo)于圖中，*1