【正文】 Ii，LIU Zenghu鋤g，MENG Qin爭(zhēng)hc．m application conditions of large turbine generator selfexcitation system【J】．Electric Powef，1994，27(12)： 61．63．[12]梁小冰，黃方能．利用EMTDC進(jìn)行長(zhǎng)持續(xù)時(shí)間過(guò)程 的仿真研究【J】．電網(wǎng)技術(shù)，2002，26(9)：55．57． LIANG Xiaobing，HUANG Fan爭(zhēng)眥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]王卉，陳楷，彭哲，等．?dāng)?shù)字仿真技術(shù)在電力系統(tǒng)中 的應(yī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)：71J(no:)=[]。P=2*i1。X1=D(i)+G(i,j)*f(i)B(i,j)*e(i)。m=p+1。X4=X1。ifJ=amp。End End case 2 Function for j=1:n J=jacci(Y,G,B,PQ,e,f,V,C,D,B2,n,ph,no)X1=G(i,j)*f(i)B(i,j)*e(i)。DX(no+1)=0。p=p+1。endP1=C(i)*e(i)+D(i)*f(i)。b(i)=v(B1(i1))*a(i)j*B1(i4)*v(B1(i))^2/2。end v=e+f*j。[C,D,DF]=xxf(G,B,e,f,P,Q,n,B2,ph,V,no)。for i=1：n End e(i)=real(B2(i3))?！尽俊尽縁unction tisco %這是一個(gè)電力系統(tǒng)潮流計(jì)算的程序 n=input（‘n請(qǐng)輸入節(jié)點(diǎn)數(shù)：n=’）； m=input(‘請(qǐng)輸入支路數(shù)：m=’)。3 2 0 0 0 0 0 1360 360。4 2 0 0 0 0 1 0 100 1 。僅供學(xué)習(xí)電力系統(tǒng)專(zhuān)業(yè)的讀者參考。endS(5)=P(5)+Q(5)*sqrt(1)。% x1為所求△x的列向量% 求節(jié)點(diǎn)電壓新值，準(zhǔn)備下一次迭代for i=1:4oa(i)=x1(i)。M(i,j)=N(i,j)。endend% 求有功、無(wú)功功率不平衡量for i=1:4for j=1:5oP(i)=oP(i)U(i)*U(j)*(G(i,j)*cos(a(i)a(j))+B(i,j)*sin(a(i)a(j)))。x2=ones(8,1)。S(4)=。endendY=0。y(4,5)=0。二 牛頓拉夫遜法潮流計(jì)算 1 基本原理牛頓法是取近似解x(k)之后，在這個(gè)基礎(chǔ)上，找到比x(k)更接近的方程的根，一步步地迭代，找到盡可能接近方程根的近似根?？吹搅俗约豪碚撀?lián)系實(shí)際的能力還需提高，也知道了自己以后學(xué)習(xí)的方向和目的。MATLAB與C語(yǔ)言和FORTRAN語(yǔ)言相比更容易被掌握。在系統(tǒng)規(guī)劃設(shè)計(jì)和安排系統(tǒng)的運(yùn)行方式時(shí)，采用離線(xiàn)潮流計(jì)算；在電力系統(tǒng)運(yùn)行狀態(tài)的實(shí)時(shí)監(jiān)控中，則采用在線(xiàn)潮流計(jì)算。所以考留到經(jīng)及調(diào)度、電網(wǎng)規(guī)劃、電力系統(tǒng)可靠性分析。通常給定的運(yùn)行條件有系統(tǒng)中各電源和負(fù)荷點(diǎn)的功率、樞紐點(diǎn)電壓、平衡點(diǎn)的電壓和相位角。（5）系統(tǒng)節(jié)點(diǎn)電壓如表4所示，―，變電所低壓母線(xiàn)35KV在35―36KV之間，經(jīng)過(guò)折算發(fā)現(xiàn)節(jié)點(diǎn)9不在指定范圍內(nèi)，需要進(jìn)行電壓調(diào)整，電壓調(diào)整后折算發(fā)現(xiàn)9節(jié)點(diǎn)已滿(mǎn)足要求，且相角隨電壓幅值而變化。對(duì)于變電所低壓母線(xiàn)為35K變壓器，非標(biāo)準(zhǔn)變比的算法與10KV的相同。（2）變電所資料：①變電所4低壓母線(xiàn)的電壓等級(jí)分別為：10KV 35KV 10KV 35KV。Matpower多用于小型電力系統(tǒng)的潮流計(jì)算分析，它運(yùn)行較為穩(wěn)定，計(jì)算速度快，運(yùn)行結(jié)果全面、直觀(guān)易懂，且準(zhǔn)確度高。最后把Matpower潮流計(jì)算的最終結(jié)果與Matlab編程和Simulink仿真兩種方法的結(jié)果做比較。根據(jù)系統(tǒng)圖及發(fā)電廠(chǎng)、變電所、輸電線(xiàn)路等參數(shù)，按照設(shè)計(jì)內(nèi)容對(duì)系統(tǒng)進(jìn)行潮流計(jì)算，并分析計(jì)算結(jié)果。（三）電網(wǎng)的節(jié)點(diǎn)設(shè)置與分類(lèi)。二、潮流計(jì)算結(jié)果分析通過(guò)運(yùn)行M文件，可得系統(tǒng)潮流計(jì)算的部分結(jié)果如圖2所示。（2）隨機(jī)抽取5條電壓調(diào)整后的支路功率進(jìn)行對(duì)比，如表9所示。該算法將網(wǎng)絡(luò)支路按層次進(jìn)行分類(lèi)，并分層并行計(jì)算各層次的支路功率損耗和電壓損耗，因而可大幅度提高配電網(wǎng)潮流的計(jì)算速度。③正常檢修及特殊運(yùn)行方式下的潮流計(jì)算,用于日運(yùn)行方式的編制,指導(dǎo)發(fā)電廠(chǎng)開(kāi)機(jī)方式,有功、無(wú)功調(diào)整方案及負(fù)荷調(diào)整方案,滿(mǎn)足線(xiàn)路、變壓器熱穩(wěn)定要求及電壓質(zhì)量要求?，F(xiàn)有很多潮流計(jì)算方法。不同領(lǐng)域、不同層次的用戶(hù)通過(guò)相應(yīng)工具的學(xué)習(xí)和應(yīng)用，可以方便地進(jìn)行計(jì)算、分析及設(shè)計(jì)工作。潮流計(jì)算的計(jì)算機(jī)算法包含高斯—賽德?tīng)柕ā⑴ｎD拉夫遜法和P—Q分解法等，其中牛拉法計(jì)算原理較簡(jiǎn)單、計(jì)算過(guò)程也不復(fù)雜，而且由于人們引入泰勒級(jí)數(shù)和非線(xiàn)性代數(shù)方程等在算法里從而進(jìn)一步提高了算法的收斂性和計(jì)算速度。為了便于用迭代法解方程組，需要將上述功率方程改寫(xiě)成功率平衡方程，并對(duì)功率平衡方程求偏導(dǎo)，得出對(duì)應(yīng)的雅可比矩陣，給未知節(jié)點(diǎn)賦電壓初值，將初值帶入功率平衡方程，得到功率不平衡量，這樣由功率不平衡量、雅可比矩陣、節(jié)點(diǎn)電壓不平衡量（未知的）構(gòu)成了誤差方程，解誤差方程，得到節(jié)點(diǎn)電壓不平衡量，節(jié)點(diǎn)電壓加上節(jié)點(diǎn)電壓不平衡量構(gòu)成節(jié)點(diǎn)電壓新的初值，將新的初值帶入原來(lái)的功率平衡方程，并重新形成雅可比矩陣，然后計(jì)算新的電壓不平衡量，這樣不斷迭代，不斷修正，一般迭代三到五次就能收斂。y(1,5)=1/(+)。endY% Y為導(dǎo)納矩陣G=real(Y)。Q=imag(S)。N(i,j)=0。oQ(i)=oQ(i)+Q(i)。M(i,i)=M(i,i)U(i)*U(j)*(G(i,j)*cos(a(i)a(j))+B(i,j)*sin(a(i)a(j)))。U(i)=U(i)+oU(i)。牛頓拉夫遜法是求解非線(xiàn)性方程的迭代過(guò)程，其計(jì)算公式為DF=JDX，式中J為所求函數(shù)的雅可比矩陣；DX為需要求的修正值；DF為不平衡的列向量。239。%% generator data% bus PgQgQmaxQmin Vg mbase status PmaxPmin =[ 4 500 0 999909999 100 1 600 0。]。ph=input(‘n請(qǐng)輸入平衡母線(xiàn)的節(jié)點(diǎn)號(hào)：ph=’)； B1=input（‘n請(qǐng)輸入支路信號(hào)：B1=’）。End f(i)=imag(B2(i3))。J=jacci(Y,G,B,P,Q,e,f,V,C,D,B2,n,ph,no)。for i=1:n hh(i)=conj(Y(ph,i)*v(i))。c(i)=v(B1(i2))*a(i)j*B1(i4)*v(B1(i2))^2/2。Q1=C(i)*f(i)D(i)*e(i)。Df(p)=Q(i)Q1。For i=(no+2):(x1+2)DX(i)=DX1(i2)。X2=G(i,j)*e(i)+B(i,j)*f(i)。J=ph J(m,q)=X3。X1=D(i)+G(i,j)*f(i)B(i,j)*e(i)。P=2*i1。J(m,q)=X3。end q=2*j1。J(m,q)