1、編號 課 程 設 計（ 2014級本科）題 目： 系（部）院： 物理與機電工程學院 專 業： 電氣工程及其自動化 作者姓名： 指導教師： 劉永科 職稱： 副教授 完成日期： 2014 年 7 月 1 日1.設計題目 復雜網絡N-R法潮流分析與計算的設計1.1基礎資料1.2系統圖的確定 選擇六節點、環網、兩電源和多引出的電力系統，簡化電力系統圖如1-1所示，等值阻抗圖如1-2所示，運用直角坐標表示的牛頓-拉夫遜計算如圖1-1所示系統中的潮流分布。計算精度要求各節點電壓的誤差或修正量不大于 .1.3各節點的初值及阻抗參數 該系統中，節點為平衡節點，保持 QUOTE 為定值，節點為PV節點，其他四個

2、節點為PQ節點，給定的注入電壓標幺值、線路阻抗標幺值、輸出功率標幺值分別為表1-3、1-4、1-5中的數據。線路對地導納標幺值之半 QUOTE 及線路阻抗標幺值、輸出功率標幺值和變壓器標幺值如圖1-2所示的注釋。表1-1 各節點電壓標幺值參數1.051.001.001.001.001.05表1-2 線路、變壓器阻抗標幺值線路T1L2L3L4L5T6阻抗j0.030.06+j0.250.04+j0.250.08+j0.300.1+j0.35J0.015表1-3 節點輸出功率節點功率2+j11.8+j0.401.6+j0.83.7+j1.35注：各PQ節點的電壓取1是為了方便計算和最后驗證程序的準

3、確性。設計的基本要求2.1設計及計算說明書（1）說明書要求書寫整齊，條理分明，表達正確、語言正確。（2）計算書內容：為各設計內容最終成果、確定提供依據進行的技術分析、論證和定量計算，如。（3）計算書要求：計算無誤，分析論證過程簡單明了，各設計內容列表匯總。2.2圖紙（1）繪制分析所需的必要圖紙（2）圖紙要求：用標準符號繪制，布置均勻，設備符號大小合適，清晰美觀。3.論文(設計)進度安排階段論文（設計）各階段名稱起止日期1熟悉設計任務書、設計題目及設計背景資料2查閱有關資料3閱讀設計要求必讀的參考資料4書寫設計說明書5上交設計成果4.需收集和閱讀的資料及參考文獻（指導教師指定）1: 陳珩.電力系

4、統穩態分析（第三版）M，北京，中國電力出版社，20072：何仰贊，溫增銀.電力系統分析第三版M，武漢,華中科技大學出版社，20023：陳悅.電氣工程畢業設計指南電力系統分冊M，北京，中國水利水電出版社，20084: 5: 教 研 室 意 見負責人簽名： 年 月 日學 院 意 見負責人簽名： 年 月 日摘 要電力系統穩態分析包括潮流計算和靜態安全分析。本文主要運用的事潮流計算，潮流計算是電力網絡設計與運行中最基本的運算，對電力網絡的各種設計方案及各種運行方式進行潮流計算，可以得到各種電網各節點的電壓，并求得網絡的潮流及網絡中的各元件的電力損耗，進而求得電能損耗。本位就是運用潮流計算具體分析，并有

5、MATLAB仿真。關鍵詞： 電力系統 潮流計算 MATLAB仿真AbstractElectric power system steady flow calculation and analysis of the static safety analysis. This paper, by means of the calculation, flow calculation is the trend of the power network design and operation of the most basic operations of electric power network, v

6、arious design scheme and the operation ways to tide computation, can get all kinds of each node of the power grid voltage and seek the trend of the network and the network of the components of the power loss, and getting electric power. The standard is to use the power flow calculation and analysis,

7、 the specific have MATLAB simulation.Key words: Power system; Flow calculation; MATLAB simulation第一章 電力系統的基本概念電力系統：發電機把機械能轉化為電能，電能經變壓器和電力線路輸送并分配到用戶，在那里經電動機、電爐和電燈等設備又將電能轉化為機械能、熱能和光能等。這些生產、變換、輸送、分配、消費電能的發電機、變壓器、變換器、電力線路及各種用電設備等聯系在一起組成的統一整體稱為電力系統。電力網：電力系統中除發電機和用電設備外的部分。動力系統：電力系統和“動力部分”的總和。電力系統潮流計算是研究電力

8、系統穩態運行情況的一種計算，它根據給定的運行條件及系統接線情況確定整個電力系統各部分的運行狀態：各母線的電壓，各元件中流過的功率，系統的功率損耗等等。在電力系統規劃的設計和現有電力系統運行方式的研究中，都需要利用潮流計算來定量地分析比較供電方案或運行方式的合理性。可靠性和經濟性。此外，電力系統潮流計算也是計算系統動態穩定和靜態穩定的基礎。所以潮流計算是研究電力系統的一種很重要和很基礎的計算。隨著科學技術的發展，電力系統變得越來越復雜，電氣工程師掌握一種好的能對電力系統進行仿真的軟件是學習和研究的需要。文章簡要介紹了MATLAB發展歷史、組成和強大的功能，并用簡單例子分別就編程和仿真兩方面分析了

9、MATIAB軟件在電力系統研究中的具體應用。采取等效電路法，能對特殊、復雜地電力系統進行高效仿真研究，因此，掌握編程和仿真是學好MATLAB的基礎。與眾多專門的電力系統仿真軟件相比，MATLAB軟件具有易學、功能強大和開放性好，是電力系統仿真研究的有力工具。第二章 潮流計算2.1潮流計算概述與發展電力系統潮流計算也分為離線計算和在線計算兩種，前者主要用于系統規劃設計和安排系統的運行方式，后者則用于正在運行系統的經常監視及實時控制。利用電子數字計算機進行電力系統潮流計算從50年代中期就已經開始。在這20年內，潮流計算曾采用了各種不同的方法，這些方法的發展主要圍繞著對潮流計算的一些基本要求進行的。

10、對潮流計算的要求可以歸納為下面幾點：1）計算方法的可靠性或收斂性；2）對計算機內存量的要求；3）計算速度；4）計算的方便性和靈活性。電力系統潮流計算問題在數學上是一組多元非線性方程式求解問題，其解法都離不開迭代。因此，對潮流計算方法，首先要求它能可靠地收斂，并給出正確答案。由于電力系統結構及參數的一些特點，并且隨著電力系統不斷擴大，潮流問題的方程式階數越來越高，對這樣的方程式并不是任何數學方法都能保證給出正確答案的。這種情況成為促使電力系統計算人員不斷尋求新的更可靠方法的重要因素。在用數字計算機解電力系統潮流問題的開始階段，普遍采取以節點導納矩陣為基礎的逐次代入法。這個方法的原理比較簡單，要求

11、的數字計算機內存量比較下，適應50年代電子計算機制造水平和當時電力系統理論水平。但它的收斂性較差，當系統規模變大時，迭代次數急劇上升，在計算中往往出現迭代不收斂的情況。這就迫使電力系統計算人員轉向以阻抗矩陣為基礎的逐次代入法。60年代初，數字計算機已發展到第二代，計算機的內存和速度發生了很大的飛躍，從而為阻抗法的采用創造了條件。阻抗法要求數字計算機儲存表征系統接線和參數的阻抗矩陣，這就需要較大的內存量。而且阻抗法每迭代一次都要求順次取阻抗矩陣中的每一個元素進行運算，因此，每次迭代的運算量很大。這兩種情況是過去電子管數字計算機無法適應的。阻抗法改善了系統潮流計算問題的收斂性，解決了導納法無法求解

12、的一些系統的潮流計算，在60年代獲得了廣泛的應用，曾為我國電力系統設計.運行和研究作出了很大的貢獻。目前，我國電力工業中仍有一些單位采用阻抗法計算潮流。阻抗法的主要缺點是占用計算機內存大，每次迭代的計算量大。當系統不斷擴大時，這些缺點就更加突出。一個內存16K的計算機在采用阻抗法時只能計算100以下的系統，32K內存的計算機也只能計算150個節點以下的系統。這樣，我國很多電力系統為了采用阻抗法計算潮流就不得不予先對系統進行相當的簡化工作。為了克服阻抗法在內存和速度方面的缺點，60年代中期發展了以阻抗矩陣為基礎的分塊阻抗法。這個方法把一個大系統分割為幾個小的地區系統，在計算機內只需要存儲各個地區

13、系統的阻抗矩陣及它們之間聯絡線的阻抗，這樣不僅大幅度地節省了內存容量，同時也提高了計算速度。克服阻抗法缺點的另一途徑是采用牛頓-拉夫遜法。這是數學中解決非線性方程式的典型方法，有較好的收斂性。在解決電力系統潮流計算問題時，是以導納矩陣為基礎的，因此，只要我們能在迭代過程中盡可能保持方程式系數矩陣的稀疏性，就可以大大提高牛頓法潮流程序的效率。自從60年代中期，在牛頓法中利用了最佳順序消去法以后，牛頓法在收斂性.內存要求.速度方面都超過了阻抗法，成為60年代末期以后廣泛采用的優秀方法。潮流計算靈活性和方便性的要求，對數字計算機的應用也是一個很關鍵的問題。過去在很長時間內，電力系統潮流計算是借助于交

14、流臺進行的。交流臺模擬了電力系統，因此在交流計算臺上計算潮流時，計算人員可以隨時監視系統各部分運行狀態是否滿足要求，如發現某些部分運行不合理，則可以立即進行調整。這樣，計算的過程就相當于運算人員去系統進行操作.調整的過程，非常直觀，物理概念也很清楚。當利用數字計算機進行潮流計算時，就失去了這種直觀性。為了彌補這個缺點，潮流程序的編制必須盡可能使計算人員在計算機計算的過程中加強對計算機過程的監視和控制，并便于作各種修改和調整。電力系統潮流計算問題并不是單純的計算問題，把它當作一個運行方式的調整問題可能更為確切。為了得到一個合理的運行方式，往往需要不斷根據計算結果，修改原始數據。在這個意義上，我們

15、在編制潮流計算程序時，對使用的方便性和靈活性必須予以足夠的重視。因此，除了要求計算方法盡可能適應各種修改.調整以外，還要注意輸入和輸出的方便性和靈活性，加強人機聯系，以便使計算人員能及時監視計算過程并適當地控制計算的進行。潮流計算的目標是求取電力系統在給定運行狀態的計算。即節點電壓和功率分布，用以檢查系統各元件是否過負荷.各點電壓是否滿足要求，功率的分布和分配是否合理以及功率損耗等。對現有電力系統的運行和擴建，對新的電力系統進行規劃設計以及對電力系統進行靜態和暫態穩定分析都是以潮流計算為基礎。潮流計算結果可用如電力系統穩態研究，安全估計或最優潮流等對潮流計算的模型和方法有直接影響。實際電力系統

16、的潮流技術那主要采用牛頓-拉夫遜法。在運行方式管理中，潮流是確定電網運行方式的基本出發點；在規劃領域，需要進行潮流分析驗證規劃方案的合理性；在實時運行環境，調度員潮流提供了電完個在預想操作情況下電網的潮流分布以校驗運行可靠性。在電力系統調度運行的多個領域都涉及到電網潮流計算。潮流是確定電力網絡運行狀態的基本因素，潮流問題是研究電力系統穩態問題的基礎和前提。牛頓-拉夫遜法作為一種實用的，有競爭力的電力系統潮流計算方法，是在應用了稀疏矩陣技巧和高斯消去法求修正方程后。牛頓-拉夫遜法是求解非線性代數方程有效的迭代計算。2.2復雜電力系統潮流計算電力系統潮流計算是對復雜電力系統正常和故障條件下穩態運行

17、狀態的計算。潮流計算的目標是求取電力系統在給定運行方式下的節點電壓和功率分布，用以檢查系統各元件是否過負荷、各點電壓是否滿足要求、功率的分布和分配是否合理以及功率損耗等。對現有電力系統的運行和擴建，對新的電力系統進行規劃設計都是以潮流計算為基礎。潮流計算結果的用途，例如用于電力系統穩定研究、安全估計或最優潮流等也對潮流計算的模型和方法有直接影響。節點類型：1）PV節點：柱入有功功率P為給定值，電壓也保持在給定數值。2）PQ節點：諸如有功功率和無功功率是給定的。3）平衡節點：用來平衡全電網的功率。選一容量足夠大的發電機擔任平衡全電網功率的職責。平衡節點的電壓大小與相位是給定的，通常以它的相角為參

18、考量，即取其電壓相角為0。一個獨立的電力網中只設一個平衡點。基本步驟：1）形成節點導納矩陣；2）將各節點電壓設初值U；3）將節點初值代入相關求式，求出修正方程式的常數項向量；4）將節點電壓初值代入求式，求出雅可比矩陣元素；5）求解修正方程，求修正向量；6）求取節點電壓的新值；7）檢查是否收斂，如不收斂，則以各節點電壓的新值作為初值自第3步重新開始進行狹義次迭代，否則轉入下一步；8）計算支路功率分布，PV節點無功功率和平衡節點柱入功率。2.3潮流計算的方法及優、缺點潮流計算法有，簡化梯度法、二次規劃法、牛頓-拉夫遜法(Newton-Raphson)等。簡化梯度法是采用梯度法進行搜索，用罰函數處理

19、違約的不等式約束。該方法程序編制簡便，所需存儲量小，對初始點無特殊要求，曾獲得普遍重視，成為第一種有效的優化潮流方法。簡化梯度法的缺點：迭代過程中，尤其是在接近最優點附近會出現鋸齒現象，收斂性較差，收斂速度很慢；每次迭代都要重新計算潮流，計算量很大，耗時較多。二次規劃法是二階的方法，解決最優潮流問題收斂精度較好，能很好地解決耦合的最優潮流問題，但缺點是計算Lagrange函數的二階偏導數，計算量大、計算復雜。2.4 MATLAB概述目前電子計算機已廣泛應用于電力系統的分析計算，潮流計算是其基本應用之一。現有很多潮流計算方法。對潮流計算方法有五方面的要求：（1）計算速度快；（2）內存需要少；（3

20、）計算結果有良好的可靠性和可信性；（4）適應性好，亦即能處理變壓器變比調整、系統元件的不同描述和與其它程序配合的能力強；（5）簡單。MATLAB是一種交互式、面向對象的程序設計語言，廣泛應用于工業界與學術界，主要用于矩陣運算，同時在數值分析、自動控制模擬、數字信號處理、動態分析、繪圖等方面也具有強大的功能。MATLAB程序設計語言結構完整，且具有優良的移植性，它的基本數據元素是不需要定義的數組。它可以高效率地解決工業計算問題，特別是關于矩陣和矢量的計算。MATLAB與C語言和FORTRAN語言相比更容易被掌握。通過M語言，可以用類似數學公式的方式來編寫算法，大大降低了程序所需的難度并節省了時間

21、，從而可把主要的精力集中在算法的構思而不是編程上。2.6牛頓-拉夫遜法原理假設有n個聯立的非線性代數方程：假設以給出各變量的初值，令其分別為個變量的修正量，使滿足以上方程，所以：將上式中的n個多元函數在初始值附近分別展開成泰勒級數，并略去含有，的二次及以上階次的各項，便得：方程可寫成：以上方程是對于修正量，的線性方程組，稱為牛頓法的修正方程，可解出，。對初始近似解進行修正： （i=1,2,，n）反復迭代，在進行k+1次迭代時，從求解修正方程式：得到修正量,對各量進行修正 （i=1,2,，n）迭代過程一直進行到滿足收斂判據圖2.1 牛頓法的幾何解釋2.7牛頓-拉夫遜法解決潮流計算問題節點總數為n

22、；PQ節點有m，；PV節點有n-m-1，平衡節點有1個，節點編號按照先PQ節點，再PV節點，最后平衡節點的順序進行編號，即:1，2，m為PQ節點；m+1，m+2，n-1為PV節點；n為平衡節點。可形成結點導納矩陣。導納矩陣元素可表示為，本文中節點電壓以直角坐標形式表示，即。由此下列公式可求出Pi，Qi假設系統中的第1，2，m號節點為PQ節點，第i個節點的給定功率為和，對該節點可列方程：假設系統中的第m+1，m+2，n-1號節點為PV節點，則對其中每一個節點可列方程: 第n號節點為平衡節點，其電壓為是給定的，故不參加迭代。修正方程可寫成分塊矩陣的形式：通過反復求解修正方程，解出各節點的未知量，再

23、通過收斂判據判定是否已為真值。從而求得PQ節點的電壓V及相角的真值，PV節點的Q、真值，平衡節點的P、Q真值，以上即為牛頓-拉夫遜迭代法的潮流計算過程，其優點為計算精確，運行速度快。其中的各個環節都可通過MATLAB程序來實現。2.8計算機潮流計算的步驟（1）對電力網絡的所有參數設初值,包括電壓、相角、有功、無功等。（2）處理非標準變比支路,使其變成標準變比為1的變壓器支路。（3）形成節點導納矩陣Y。（4）計算有功功率的不平衡量Pi,從而求出。（5）根據節點的類型形成J。（6）解修正方程式,求各節點的電壓的變化量ei(i=1,2,3.n,is)（7）求各節點相角的新值ei=ei+ei (i=1

24、,2,3. n,is)（8）計算無功功率的不平衡量Qi,從而求出 (i=1,2,3.n,is)（9）解修正方程式,求各節點的電壓大小的變化量 (i=1,2,3.,n,is)。（10）求各節點的電壓大小的新值 (i=1,2,3.,n,is)。（11）運用個節點的電壓的新值自第四步開始下一次迭代。計算平衡節點的功率和線路功率。其中平衡節點的功率的計算公式為線路上的功率為:從而線路上的損耗的功率為:2.9計算機程序的實現導納矩陣的形成節點導納矩陣是方陣,其階數等于網絡中出參考節點外的節點數n。節點導納矩陣是稀疏矩陣,其各行非零非對角元數就等于該行相對應節點所連接的不接地支路數。節點導納矩陣的對角元就

25、等于各該節點所連接導納的總和。節點導納矩陣的非對角元等于連接節點i,j支路導納的負值。點導納矩陣一般是對稱矩陣。對于支路中有非標準變比變壓器的支路來說,利用下面的公式來計算它的導納。變壓器的變比,假如已知非標準變比支路i,j上的阻抗(以下沒有特殊說明所有的參數都用標幺值)為,則線路導納為,線路上的對地半導納為。J的形成Y是由最終形成的導納矩陣的虛部組成的,但是pv節點以及平衡節點不參加QV迭代,因此Y中不包含與這些節點有關的元素。迭代條件和約束方程迭代條件就是如果Q時就停止迭代。對節點的約束條件分為三類:即對節點注入功率的約束、對節點電壓大小的約束和對相角的約束。其中對節點注入功率的約束,主要

26、是對電源注入功率的約束條件不能滿足時,將威脅到發電機的安全運行。對電壓大小的約束不能滿足時,將影響電能的質量,嚴重時將影響系統運行的穩定性。對相對相角的約束條件不能滿足時,也將危及系統運行的穩定性。圖2.2 程序流程圖圖3.7 線路潮流大小分布情況3.0 手算潮流計算用圖1和圖2的數據和等值網絡形成節點導納矩陣1. 節點導納矩陣導納矩陣 Y=0 -30.2343i 0 +31.7460i 0 0 0 0 0 +31.7460i 14.8252 -42.6506i -14.2012 + 5.9172i 0 -0.6240 + 3.9002i 0 0 -14.2012+5.9172i 15.031

27、1-8.5292i -0.8299+3.1120i 0 0 0 0 -0.8299 + 3.1120i 1.5846 - 5.5035i -0.7547 + 2.6415i 0 0 -0.6240 + 3.9002i 0 -0.7547 + 2.6415i 1.3787 -72.9583i 0 +63.4921i0 0 0 0 0 +63.4921i 0 -60.4686iB1 =1.0000 2.0000 0 + 0.0300i 0 1.0500 1.0000 1.0000 2.0000 3.0000 0.0600 + 0.0250i 0 + 0.5000i 1.0000 0 0 2.000

28、0 5.0000 0.0400 + 0.2500i 0 + 0.5000i 1.0000 0 0 3.0000 4.0000 0.0800 + 0.3000i 0 + 0.5000i 1.0000 0 0 4.0000 5.0000 0.1000 + 0.3500i 0 1.0000 0 0 6.0000 5.0000 0 + 0.0150i 0 1.0500 1.0000 1.0000 B2 =0 0 1.2000 0 1.0000 0 2.1000 + 1.0000i 1.0000 0 2.0000 0 1.8000 + 0.4000i 1.0000 0 2.0000 0 1.6000 +

29、 0.8000i 1.0000 0 2.0000 0 3.7000 + 1.3000i 1.0000 0 2.0000 0 -5.0000 1.2000 0 3.0000 功率方程第(1)次差值： 0 0 4.2619 -2.1000 0.1000 -1.8000 -0.5500 -1.6000 8.4738 -3.7000 0 5.0000形成的第(1)次Jacobi矩陣： Columns 1 through 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 00 0 -37.3887 14.8252 5.9172 -14.2012 0 0 3

30、.9002 -0.6240 0 00 0 -14.8252 -47.9126 14.2012 5.9172 0 0 0.6240 3.9002 0 00 0 5.9172 -14.2012 -8.0292 15.0311 3.1120 -0.8299 0 0 0 0 0 0 14.2012 5.9172 -15.0311 -9.0292 0.8299 3.1120 0 0 0 00 0 0 0 3.1120 -0.8299 -5.2535 1.5846 2.6415 -0.7547 0 0 0 0 0 0 0.8299 3.1120 -1.5846 -5.7535 0.7547 2.6415

31、0 0 0 0 3.9002 -0.6240 0 0 2.6415 -0.7547 -63.1845 1.3787 63.4921 0 0 0 .6240 3.9002 0 0 0.7547 2.6415 -1.3787 -82.7321 0 63.49210 0 0 0 0 0 0 0 0 0 -2.4000 00 0 0 0 0 0 0 0 0 76.1905 0 -63.4921 Column 13 0 0 4.2619 -2.1000 0.1000 -1.8000 -0.5500 -1.6000 8.4738 -3.7000 0 5.0000Jacobi矩陣第(1)次消去運算 Colu

32、mns 1 through 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1.0000 -0.3965 -0.1583 0.3798 0 0 -0.1043 0.0167 0 0 0 0 0 1.0000 -0.2204 -0.2147 0 0 0.0171 -0.0771 0 0 0 0 0 0 1.0000 -1.0549 -0.3206 0.0855 -0.0845 0.1044 0 0 0 0 0 0 0 1.0000 0.1079 -0.1753 -0.0184 -0.0757 0 0 0 0 0 0 0 0 1.000

33、0 -0.3868 -0.6525 0.1977 0 0 0 0 0 0 0 0 0 1.0000 0.0419 -0.5520 0 0 0 0 0 0 0 0 0 0 1.0000 -0.0135 -1.0440 0 0 0 0 0 0 0 0 0 0 1.0000 0.0106 -0.7856 0 0 0 0 0 0 0 0 0 0 1.0000 0 0 0 0 0 0 0 0 0 0 0 0 1.0000 Column 13 0 0 -0.1140 0.0705 -0.1659 0.1184 0.0717 0.3118 -0.1439 0.0636 0 -0.0431Jacobi矩陣第(

34、1)次回代運算 0 0 -0.1490 0.1103 -0.0024 0.1657 0.1015 0.3342 -0.1435 0.0297 0 -0.0431各個節點電壓模 0 1.1543 1.0160 0.9586 1.1439 1.2008功率方程第(2)次差值： 0 0 -0.9196 -0.3298 -0.2847 -0.0127 -0.5652 0.0384 -1.5185 0.1960 -0.0019 -0.3928形成的第(2)次Jacobi矩陣。 Column 13 0 0 -0.9196 -0.3298 -0.2847 -0.0127 -0.5652 0.0384 -1.

35、5185 0.1960 -0.0019 -0.3928Jacobi矩陣第(2)次消去運算。 Column 13 0 0 0.0195 -0.0011 0.0425 -0.0053 0.2265 -0.0488 0.0271 -0.0013 0.0008 0.0487Jacobi矩陣第(2)次回代運算 0 0 0.0441 0.0155 0.0700 -0.0210 0.1790 -0.0275 0.0288 0.0422 -0.0010 0.0487各個節點電壓模 0 1.1120 0.9435 0.7821 1.1170 1.2010功率方程第(3)次差值： 0 0 -0.0417 0.00

36、17 -0.0185 -0.0126 -0.1340 -0.0217 -0.0367 0.0646 -0.0024 -0.0916形成的第(3)次Jacobi矩陣。 Column 13 0 0 -0.0417 0.0017 -0.0185 -0.0126 -0.1340 -0.0217 -0.0367 0.0646 -0.0024 -0.0916Jacobi矩陣第(3)次消去運算。 Column 13 0 0 0.0009 -0.0004 0.0034 -0.0004 0.1046 -0.0114 0.0036 -0.0002 0.0010 0.0181Jacobi矩陣第(3)次回代運算 0

37、0 0.0078 0.0043 0.0191 -0.0063 0.0768 -0.0055 0.0056 0.0156 0.0011 0.0181各個節點電壓模 0 1.1048 0.9237 0.7098 1.1126 1.2001功率方程第(4)次差值： 0 0 -0.0004 0.0002 -0.0002 -0.0012 -0.0276 -0.0058 0.0006 0.0016 -0.0003 -0.0053形成的第(4)次Jacobi矩陣。 Column 13 0 0 -0.0004 0.0002 -0.0002 -0.0012 -0.0276 -0.0058 0.0006 0.00

38、16 -0.0003 -0.0053Jacobi矩陣第(4)次消去運算。 Column 13 0 0 0.0000 -0.0000 0.0000 0.0000 0.0450 -0.0033 0.0010 0.0002 0.0001 0.0053Jacobi矩陣第(4)次回代運算 0 0 0.0021 0.0014 0.0056 -0.0020 0.0259 -0.0018 0.0016 0.0048 0.0002 0.0053各個節點電壓模 0 1.1029 0.9178 0.6857 1.1114 1.2000功率方程第(5)次差值： 0 0 -0.0000 0.0000 0.0000 -0

39、.0001 -0.0032 -0.0007 -0.0000 0.0001 -0.0000 -0.0005形成的第(5)次Jacobi矩陣。 Column 13 0 0 -0.0000 0.0000 0.0000 -0.0001 -0.0032 -0.0007 -0.0000 0.0001 -0.0000 -0.0005Jacobi矩陣第(5)次消去運算。 Column 13 0 0 0.0000 -0.0000 -0.0000 0.0000 0.0087 -0.0005 0.0002 0.0000 0.0000 0.0008Jacobi矩陣第(5)次回代運算 0 0 0.0003 0.0002

40、 0.0008 -0.0003 0.0039 -0.0003 0.0002 0.0007 0.0000 0.0008各個節點電壓模 0 1.1026 0.9169 0.6820 1.1113 1.2000功率方程第(6)次差值： 1.0e-004 * 0 0 -0.0026 0.0007 0.0032 -0.0182 -0.7419 -0.1647 -0.0016 0.0184 -0.0063 -0.1036形成的第(6)次Jacobi矩陣。 Column 13 0 0 -0.0000 0.0000 0.0000 -0.0000 -0.0001 -0.0000 -0.0000 0.0000 -

41、0.0000 -0.0000Jacobi矩陣第(6)次消去運算。 Column 13 0 0 0.0000 -0.0000 -0.0000 0.0000 0.0002 -0.0000 0.0000 0.0000 0.0000 0.0000Jacobi矩陣第(6)次回代運算 1.0e-004 * 0 0 0.0764 0.0496 0.2048 -0.0728 0.9503 -0.0641 0.0564 0.1753 0.0074 0.1923各個節點電壓模 0 1.1026 0.9169 0.6819 1.1113 1.2000功率方程第(7)次差值： 1.0e-007 * 0 0 -0.0015 0.0004 0.0019 -0.0106 -0.4343