【正文】 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)技術(shù)，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ù)字仿真技術(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)：71J(:,no)=[]。J(no:)=[]。end J(p,q)=X1。X4= C(i)+G(i,j)*f(i)B(i,j)*e(i)。X2=C(i)+G(i,j)*e(i)+B(i,j)*f(i)。Else if j=amp。J(p,q)=X2。J(p,q)=X1。p=2*i1。X3=X2。X1=G(i)*f(i)B(i,j)*e(i)。case 1 J(p,q)=X1。%該子程序是用來求取jacci矩陣for i=1:n X3=0。Df(k)=DX(i)。End ElseDX=[DX1,0,0]。x1=length(DX1)。end end end Df=Df’。else p=2*i1。V1=e(i)^2+f(i)^2。For j=i:nC(i)=C(i)+G(i,j)*e(j)B(i,j)*f(j)。end result2=[B1(:,6),B1(:,1),B1(:,2),real(b’),imag(b’),real(c’),imag(c’), real(b’+c’),imag(b’+c’)]。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))]。end S(ph)=sum(hh)*v(ph)。e=e+De；f=f+Df。[De,Di]=hxf(J,D,F,ph,n,no)。end Y(p,p)=Y(p,p)+1./(B1(i3)*B1(i5)^2+B1(i4)./2P=real(S)。For i=1:m S(i)=B2(i1)B2(i2)。if A（i2）=0B=imag(Y)。%它以矩陣形式存貯支路的情況，每行存貯一條支路 %第一列存貯支路的一個端點 %第二列存貯支路的另一個端點 %第三列存貯支路阻抗%第四列存貯支路的對地導納%第五列存貯變壓器的變比，注意支路為1 %第六列存貯支路的序號B2=input(‘n請輸入節(jié)點信息：B2=’)； %第一列為電源側(cè)的功率 %第二列為負荷側(cè)的功率 %第三列為該點的電壓值%第四列為該點的類型：1為PQ，2為PV節(jié)點，3為平衡節(jié)點 A=input（‘n請輸入節(jié)點號及對地阻抗：A=’）； ip=input（‘n請輸入修正值：ip=’）； %ip為修正值)。六 結(jié)論通過這個任務(wù)，自己在matlab編程，潮流計算，word文檔的編輯功能等方面均有提高，但也暴漏出一些問題：理論知識儲備不足，對matlab的性能和特點還不能有一個全面的把握，對word軟件也不是很熟練，相信通過以后的學習能彌補這些不足，達到一個新的層次。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% 當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為不平衡量列向量% 求雅克比矩陣% 當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é)點導納矩陣中互導納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é)點導納矩陣Y；2)設(shè)各節(jié)點電壓的初值，如果是直角坐標的話設(shè)電壓的實部e和虛部f；如果是極坐標的話則設(shè)電壓的幅值U和相角a；3)將各個節(jié)點電壓的初值代入公式求修正方程中的不平衡量以及修正方程的系數(shù)矩陣的雅克比矩陣；4)解修正方程式，求各節(jié)點電壓的變化量，即修正量； 5)計算各個節(jié)點電壓的新值，即修正后的值；6)利用新值從第（3）步開始進入下一次迭代，直至達到精度退出循環(huán)； 7)計算平衡節(jié)點的功率和線路功率，輸出最后計算結(jié)果； ① 公式推導② 流程圖三matlab編程代碼clear。牛頓迭代法其最大優(yōu)點是在方程f(x)=0的單根附近時誤差將呈平方減少，而且該法還可以用來求方程的重根、復根。一 研究內(nèi)容通過一道例題來認真分析牛頓拉夫遜法的原理和方法（采用極坐標形式的牛拉法），同時掌握潮流計算計算機算法的相關(guān)知識，能看懂并初步使用MATLAB軟件進行編程，培養(yǎng)自己電力系統(tǒng)潮流計算機算法編程能力。本文通過對電力系統(tǒng)模型采用PSASP軟件進行潮流計算，可以得出PSASP能夠較好地分析的潮流分布和負荷情況，能夠滿足電力系統(tǒng)運行分析與規(guī)劃的數(shù)字仿真需要，從而保證電網(wǎng)的安全穩(wěn)定運行。點擊“編輯”按鈕填寫有關(guān)數(shù)據(jù)，點擊“刷新”按鈕保存數(shù)據(jù)。圖1電力系統(tǒng)模型 Figure 1 Power System Model單線圖完成后需要進行方案定義，通常定義成常規(guī)型?？赏ㄟ^UD功能實現(xiàn)自動發(fā)電控制、可控串補及嵌套FACTS元件功能； 潮流計算有以下功能：內(nèi)存要求?？朔杩狗ㄈ秉c是另一個途徑是采用牛頓拉夫遜法。這就迫使電力系統(tǒng)的計算人員轉(zhuǎn)向以阻抗矩陣為基礎(chǔ)的逐次代入法。潮流計算結(jié)果可用如電力系統(tǒng)穩(wěn)態(tài)研究，安全估計或最優(yōu)潮流等對潮流計算的模型和方法有直接影響。 潮流計算的意義及其發(fā)展電力系統(tǒng)潮流計算是電力系統(tǒng)分析中的一種最基本的計算，是對復雜電力系統(tǒng)正常和故障條件下穩(wěn)態(tài)運行狀態(tài)的計算。各元件中流過的功率，系統(tǒng)的功率損耗等等。運營方式管理中，潮流是確定電網(wǎng)運行方式的基本出發(fā)點：在規(guī)劃領(lǐng)域，需要進行潮流分析驗證規(guī)劃方案的合理性；在實時運行環(huán)境，調(diào)度員潮流提供了電網(wǎng)在預想操作預想下的電網(wǎng)的潮流分布以及校驗運行的可靠性。電力系統(tǒng)潮流計算概述電力工業(yè)發(fā)展初期，電能是直接在用戶附近的發(fā)電站中生產(chǎn)的，各發(fā)電站孤立運行。手工的潮流計算工作大量且易于出錯。除它自身的重要作用之外，潮流計算還是網(wǎng)損計算、靜態(tài)安全分析、暫態(tài)穩(wěn)定計算、小干擾靜態(tài)穩(wěn)定計算、短路計算、靜態(tài)和動態(tài)等值計算的基礎(chǔ)。本文通過psasp電力系統(tǒng)分析綜合程序軟件進行簡單電力系統(tǒng)潮流計算仿真，利用牛頓法進行計算潮流。潮流計算是電力系統(tǒng)非常重要的分析計算，用以研究系統(tǒng)規(guī)劃和運行中提出的各種問題。psasp0 引言潮流計算是電力系統(tǒng)分析最基本的計算。本文計算采用直角坐標形式下的牛頓拉夫遜法，牛頓拉夫遜法有很好的收斂性，但要求有合適的初值。主要用于電力系統(tǒng)規(guī)劃設(shè)計人員確定經(jīng)濟合理、技術(shù)可行的規(guī)劃設(shè)計方案；運行調(diào)度人員確定系統(tǒng)運行方式、分析系統(tǒng)事故、尋求反事故措施；科研人員研究新設(shè)備、新元件投入系統(tǒng)等新問題以及高等院校用于教學和研究。這種由發(fā)電機、升壓和降壓變電所，送電線路以及用電設(shè)備有機連接起來的整體，即稱為電力系統(tǒng)。電力系統(tǒng)潮流計算是研究電力系統(tǒng)穩(wěn)態(tài)運行情況的一種計算，它根據(jù)給定的運行條件及系統(tǒng)接線情況確定整個電力系統(tǒng)各部分的運行狀態(tài)：各母線的電壓。此外，電力系統(tǒng)的潮流計算也是計算機系統(tǒng)動態(tài)穩(wěn)定和靜態(tài)穩(wěn)定的基礎(chǔ)，所以潮流計算是研究電力系統(tǒng)的一種和重要和基礎(chǔ)的計算。對現(xiàn)有的電力系統(tǒng)的運行和擴建，對新的電力系統(tǒng)進行規(guī)劃設(shè)計以及對電力系統(tǒng)進行靜態(tài)和穩(wěn)態(tài)分析都是以潮流計算為基礎(chǔ)。這個方法的原理比較簡單，要求的數(shù)字計算機內(nèi)存量比較差下，適應50年代電子計算機制造水平和當時電力系統(tǒng)理論水平，但它的收斂性較差，當系統(tǒng)規(guī)模變大時，迭代次數(shù)急劇上升，在計算中往往出現(xiàn)迭代不收斂的情況。這個方法把一個大系統(tǒng)分割為幾個小的地區(qū)系統(tǒng)，在計算機內(nèi)只需要存儲各個地區(qū)系統(tǒng)的阻抗矩陣及它們之間聯(lián)絡(luò)的阻抗，這樣不僅大幅度的節(jié)省了內(nèi)存容量，同時也提高了計算速度。自從60年代中期，牛頓法中利用了最佳順序消去法以后，牛頓法在收斂性。包括:穩(wěn)態(tài)分析的潮流計算、網(wǎng)損分析、最優(yōu)潮流和無功優(yōu)化、靜態(tài)安全分析、諧波分析、靜態(tài)等值等?？煽紤]自動調(diào)節(jié)變壓器和移相器；此系統(tǒng)中共有1個平衡節(jié)點、1個PV節(jié)點及3個PQ節(jié)點。采用牛頓算法，比較容易收斂，迭代次數(shù)設(shè)置為100次。圖4潮流計算摘要信息報表Figure 4 flow calculation summary information report圖5物理母線報表 Figure 5 physical bus report圖6負荷結(jié)果報表 Figure 6 load resulting report圖7結(jié)果綜述報表Figure 7 results summary reports圖8交流線結(jié)果報表 Figure 8 AC line results report 結(jié)束語潮流計算是電力系統(tǒng)非常重要的分析計算，用以研究系統(tǒng)規(guī)劃和運行中提出的各種問題。同時基于MATLAB的計算機算法以雙精度類型進行數(shù)據(jù)的存儲和運算, 數(shù)據(jù)精確度高，能進行潮流計算中的各種矩陣運算，使得傳統(tǒng)潮流計算方法更加優(yōu)化。二 牛頓拉夫遜法潮流計算 1 基本原理牛頓法是取近似解x(k)之后，在這個基礎(chǔ)上，找到比x(k)更接近的方程的根，一步步地迭代，找到盡可能接近方程根的近似根。2 基本步驟和設(shè)計流程圖形成了雅克比矩陣并建立了修正方程式，運用牛頓拉夫遜法計算潮流的核心問題已經(jīng)解決，已有可能列出基本計算步驟并編制流程圖。y(4,5)=0。y(2,3)=1/(+)。endendY=0。B=imag(Y)。S(4)=。% 賦初值，U為節(jié)點電壓的幅值，a為節(jié)點電壓的相位角U=ones(1,5)。x2=ones(8,1)。M(i,j)=0。endend% 求有功、無功功率不平衡量for i=1:4f