【文章內容簡介】 (inv(J))*x2)。% x1為所求△x的列向量% 求節(jié)點電壓新值，準備下一次迭代for i=1:4oa(i)=x1(i)。oU(i)=x1(i+4)*U(i)。endfor i=1:4a(i)=a(i)+oa(i)。U(i)=U(i)+oU(i)。endk=k+1。endk,U,a% 求節(jié)點注入功率i=5。for j=1:5P(i)=U(i)*U(j)*(G(i,j)*cos(a(i)a(j))+B(i,j)*sin(a(i)a(j)))+P(i)。Q(i)=U(i)*U(j)*(G(i,j)*sin(a(i)a(j))B(i,j)*cos(a(i)a(j)))+Q(i)。endS(5)=P(5)+Q(5)*sqrt(1)。S% 求節(jié)點注入電流I=Y*U39。四運行結果節(jié)點導納矩陣經(jīng)過五次迭代后的雅克比矩陣迭代次數(shù)以及節(jié)點電壓的幅值和相角（弧度數(shù)）節(jié)點注入功率和電流五 結果分析在這次學習和實際操作過程里：首先，對電力系統(tǒng)分析中潮流計算的部分特別是潮流計算的計算機算法中的牛頓拉夫遜法進行深入的研讀，弄明白了其原理、計算過程、公式推導以及設計流程。牛頓拉夫遜法是求解非線性方程的迭代過程，其計算公式為DF=JDX，式中J為所求函數(shù)的雅可比矩陣；DX為需要求的修正值；DF為不平衡的列向量。利用x(*)=x(k+1)+DX(k+1)進行多次迭代，通過迭代判據(jù)得到所需要的精度值即準確值x(*)。六 結論通過這個任務，自己在matlab編程，潮流計算，word文檔的編輯功能等方面均有提高，但也暴漏出一些問題：理論知識儲備不足，對matlab的性能和特點還不能有一個全面的把握，對word軟件也不是很熟練，相信通過以后的學習能彌補這些不足，達到一個新的層次。第三篇：電力系統(tǒng)潮流計算發(fā)展史電力系統(tǒng)潮流計算發(fā)展史對潮流計算的要求可以歸納為下面幾點：（1）算法的可靠性或收斂性（2）計算速度和內存占用量（3）計算的方便性和靈活性電力系統(tǒng)潮流計算屬于穩(wěn)態(tài)分析范疇，不涉及系統(tǒng)元件的動態(tài)特性和過渡過程。因此其數(shù)學模型不包含微分方程，是一組高階非線性方程。非線性代數(shù)方程組的解法離不開迭代，因此，潮流計算方法首先要求它是能可靠的收斂，并給出正確答案。隨著電力系統(tǒng)規(guī)模的不斷擴大，潮流問題的方程式階數(shù)越來越高，目前已達到幾千階甚至上萬階，對這樣規(guī)模的方程式并不是采用任何數(shù)學方法都能保證給出正確答案的。這種情況促使電力系統(tǒng)的研究人員不斷尋求新的更可靠的計算方法。在用數(shù)字計算機求解電力系統(tǒng)潮流問題的開始階段，人們普遍采用以節(jié)點導納矩陣為基礎的高斯賽德爾迭代法（一下簡稱導納法）。這個方法的原理比較簡單，要求的數(shù)字計算機的內存量也比較小，適應當時的電子數(shù)字計算機制作水平和電力系統(tǒng)理論水平，于是電力系統(tǒng)計算人員轉向以阻抗矩陣為主的逐次代入法（以下簡稱阻抗法）。20世紀60年代初，數(shù)字計算機已經(jīng)發(fā)展到第二代，計算機的內存和計算速度發(fā)生了很大的飛躍，從而為阻抗法的采用創(chuàng)造了條件。阻抗矩陣是滿矩陣，阻抗法要求計算機儲存表征系統(tǒng)接線和參數(shù)的阻抗矩陣。這就需要較大的內存量。而且阻抗法每迭代一次都要求順次取阻抗矩陣中的每一個元素進行計算，因此，每次迭代的計算量很大。阻抗法改善了電力系統(tǒng)潮流計算問題的收斂性，解決了導納法無法解決的一些系統(tǒng)的潮流計算，在當時獲得了廣泛的應用，曾為我國電力系統(tǒng)設計、運行和研究作出了很大的貢獻。但是，阻抗法的主要缺點就是占用計算機的內存很大，每次迭代的計算量很大。當系統(tǒng)不斷擴大時，這些缺點就更加突出。為了克服阻抗法在內存和速度方面的缺點，后來發(fā)展了以阻抗矩陣為基礎的分塊阻抗法。這個方法把一個大系統(tǒng)分割為幾個小的地區(qū)系統(tǒng)，在計算機內只需存儲各個地區(qū)系統(tǒng)的阻抗矩陣及它們之間的聯(lián)絡線的阻抗，這樣不僅大幅度的節(jié)省了內存容量，同時也提高了節(jié)省速度?？朔杩狗ㄈ秉c的另一途徑是采用牛頓拉夫遜法（以下簡稱牛頓法）。牛頓法是數(shù)學中求解非線性方程式的典型方法，有較好的收斂性。解決電力系統(tǒng)潮流計算問題是以導納矩陣為基礎的，因此，只要在迭代過程中盡可能保持方程式系數(shù)矩陣的稀疏性，就可以大大提高牛頓潮流程序的計算效率。自從20世紀60年代中期采用了最佳順序消去法以后，牛頓法在收斂性、內存要求、計算速度方面都超過了阻抗法，成為直到目前仍被廣泛采用的方法。在牛頓法的基礎上，根據(jù)電力系統(tǒng)的特點，抓住主要矛盾，對純數(shù)學的牛頓法進行了改造，得到了PQ分解法。PQ分解法在計算速度方面有顯著的提高，迅速得到了推廣。牛頓法的特點是將非線性方程線性化。20世紀70年代后期，有人提出采用更精確的模型，即將泰勒級數(shù)的高階項也包括進來，希望以此提高算法的性能，這便產(chǎn)生了保留非線性的潮流算法。另外，為了解決病態(tài)潮流計算，出現(xiàn)了將潮流計算表示為一個無約束非線性規(guī)劃問題的模型，即非線性規(guī)劃潮流算法。近20多年來，潮流算法的研究仍然非?；钴S，但是大多數(shù)研究都是圍繞改進牛頓法和PQ分解法進行的。此外，隨著人工智能理論的發(fā)展，遺傳算法、人工神經(jīng)網(wǎng)絡、模糊算法也逐漸被引入潮流計算。但是，到目前為止這些新的模型和算法還不能取代牛頓法和PQ分解法的地位。由于電力系統(tǒng)規(guī)模的不斷擴大，對計算速度的要求不斷提高，計算機的并行計算技術也將在潮流計算中得到廣泛的應用，成為重要的研究領域。第四篇：電力系統(tǒng)潮流計算程序設計電力系統(tǒng)潮流計算程序設計姓名：韋應順學號：2011021052 電力工程學院牛頓—拉夫遜潮流計算方法具有能夠將非線性方程線性化的特點，而使用MATLAB語言是由于MATLAB語言的數(shù)學邏輯強，易編譯?！尽俊尽縁unction tisco %這是一個電力系統(tǒng)潮流計算的程序 n=input（‘n請輸入節(jié)點數(shù)：n=’）； m=input(‘請輸入支路數(shù)：m=’)。ph=input(‘n請輸入平衡母線的節(jié)點號：ph=’)； B1=input（‘n請輸入支路信號：B1=’）。%它以矩陣形式存貯支路的情況，每行存貯一條支路 %第一列存貯支路的一個端點 %第二列存貯支路的另一個端點 %第三列存貯支路阻抗%第四列存貯支路的對地導納%第五列存貯變壓器的變比，注意支路為1 %第六列存貯支路的序號B2=input(‘n請輸入節(jié)點信息：B2=’)； %第一列為電源側的功率 %第二列為負荷側的功率 %第三列為該點的電壓值%第四列為該點的類型：1為PQ，2為PV節(jié)點，3為平衡節(jié)點 A=input（‘n請輸入節(jié)點號及對地阻抗：A=’）； ip=input（‘n請輸入修正值：ip=’）； %ip為修正值)。Y=zeros（n）；Y(p,q)=Y(p,q)1./(B1(i3)*B1(i5)。e=zeros（1，n）；Y(p,q)=Y(p,q)。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)。if A（i2）=0B=imag(Y)。p=A(i1)；Y(p p)=1./A(i2)。for i=1：n End e(i)=real(B2(i3))。End f(i)=imag(B2(i3))。For i=1:m S(i)=B2(i1)B2(i2)。p=B1(i1)。V(i)=B2(i3)。p=B1(i2)。end Y(p,p)=Y(p,p)+1./(B1(i3)*B1(i5)^2+B1(i4)./2P=real(S)。Q=imag(S)。[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,Di]=hxf(J,D,F,ph,n,no)。t=0。whilemax(abs(De))ipamp。max(abs(Dfi)ipt=t+1。e=e+De；f=f+Df。[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)。end v=e+f*j。for i=1:n h