電力系統(tǒng)分析課程設(shè)計(jì)

1、電力系統(tǒng)分析課程設(shè)計(jì)學(xué) 院 電氣工程學(xué)院 班 級 學(xué) 號 姓 名 指導(dǎo)教師 時 間 前言 電力系統(tǒng)潮流計(jì)算是研究電力系統(tǒng)穩(wěn)態(tài)運(yùn)行情況的一種計(jì)算，它根據(jù)給定的運(yùn)行條件及系統(tǒng)接線情況確定整個電力系統(tǒng)各部分的運(yùn)行狀態(tài)：各母線的電壓，各元件中流過的功率，系統(tǒng)的功率損耗等等。在電力系統(tǒng)規(guī)劃的設(shè)計(jì)和現(xiàn)有電力系統(tǒng)運(yùn)行方式的研究中，都需要利用潮流計(jì)算來定量地分析比較供電方案或運(yùn)行方式的合理性?？煽啃院徒?jīng)濟(jì)性。此外，電力系統(tǒng)潮流計(jì)算也是計(jì)算系統(tǒng)動態(tài)穩(wěn)定和靜態(tài)穩(wěn)定的基礎(chǔ)。所以潮流計(jì)算是研究電力系統(tǒng)的一種很重要和很基礎(chǔ)的計(jì)算。 隨著科學(xué)技術(shù)的發(fā)展，電力系統(tǒng)變得越來越復(fù)雜，電氣工程師掌握一種好的能對電力系統(tǒng)進(jìn)行仿真的

2、軟件是學(xué)習(xí)和研究的需要。與眾多專門的電力系統(tǒng)仿真軟件相比，MATLAB軟件具有易學(xué)、功能強(qiáng)大和開放性好，是電力系統(tǒng)仿真研究的有力工具。目錄1 設(shè)計(jì)題目31.1系統(tǒng)圖的確定31.2各節(jié)點(diǎn)的初值及阻抗參數(shù)42 潮流計(jì)算52.1潮流計(jì)算概述與發(fā)展52.2復(fù)雜電力系統(tǒng)潮流計(jì)算52.3 MATLAB概述62.4牛頓-拉夫遜法原理72.5牛頓-拉夫遜法解決潮流計(jì)算問題82.6計(jì)算機(jī)潮流計(jì)算的步驟92.7計(jì)算機(jī)潮流計(jì)算流程圖113 手算潮流計(jì)算123.1確定節(jié)點(diǎn)類型123.2求節(jié)點(diǎn)導(dǎo)納矩陣Yb123.3計(jì)算各節(jié)點(diǎn)功率的修正方程的初始值（不平衡量）134 計(jì)算機(jī)算法潮流計(jì)算154.1計(jì)算機(jī)法潮流計(jì)算過程154

3、.2計(jì)算機(jī)法潮流計(jì)算結(jié)果234.3系統(tǒng)功率分布圖25個人心得29附錄：源程序30參考文獻(xiàn)401 設(shè)計(jì)題目1.1系統(tǒng)圖的確定選擇六節(jié)點(diǎn)、環(huán)網(wǎng)、兩電源和多引出的電力系統(tǒng)，簡化電力系統(tǒng)圖如圖1-1所示，等值導(dǎo)納圖如圖1-2所示。運(yùn)用以直角坐標(biāo)表示的牛頓-拉夫遜計(jì)算如圖1-1所示系統(tǒng)中的潮流分布。計(jì)算精度要求各節(jié)點(diǎn)電壓的誤差或修正量不大于。 圖1-1 電力系統(tǒng)圖圖1-2 電力系統(tǒng)等值導(dǎo)納圖1.2各節(jié)點(diǎn)的初值及阻抗參數(shù)該系統(tǒng)中，節(jié)點(diǎn)為平衡節(jié)點(diǎn)，保持=1.05+j0為定值，節(jié)點(diǎn)為PV節(jié)點(diǎn)，其他四個節(jié)點(diǎn)都是PQ節(jié)點(diǎn)。給定的注入電壓標(biāo)幺值如表1-1、線路阻抗標(biāo)幺值如表1-2、輸出功率標(biāo)幺值如表1-3。表1-1

4、 各節(jié)點(diǎn)電壓標(biāo)幺值參數(shù)UUUUUU1.051.001.001.001.001.05表1-2 線路、變壓器阻抗標(biāo)幺值 線路L2L3L4L5T1T2Y/2阻抗0.06+j0.250.04+j0.250.08+j0.300.1+j0.35j0.03j0.015j0.25表1-3 節(jié)點(diǎn)輸出功率節(jié)點(diǎn)功率2+j11.8+j0.41.6+j0.83.7+j1.35注：各PQ節(jié)點(diǎn)的電壓取1是為了方便計(jì)算和最后驗(yàn)證程序的正確性。2 潮流計(jì)算2.1潮流計(jì)算概述與發(fā)展 電力系統(tǒng)潮流計(jì)算也分為離線計(jì)算和在線計(jì)算兩種，前者主要用于系統(tǒng)規(guī)劃設(shè)計(jì)和安排系統(tǒng)的運(yùn)行方式，后者則用于正在運(yùn)行系統(tǒng)的經(jīng)常監(jiān)視及實(shí)時控制。 利用電子數(shù)

5、字計(jì)算機(jī)進(jìn)行電力系統(tǒng)潮流計(jì)算從50年代中期就已經(jīng)開始。在這20年內(nèi)，潮流計(jì)算曾采用了各種不同的方法，這些方法的發(fā)展主要圍繞著對潮流計(jì)算的一些基本要求進(jìn)行的。 牛頓-拉夫遜法作為一種實(shí)用的，有競爭力的電力系統(tǒng)潮流計(jì)算方法，是在應(yīng)用了稀疏矩陣技巧和高斯消去法求修正方程后。牛頓-拉夫遜法是求解非線性代數(shù)方程有效的迭代計(jì)算。2.2復(fù)雜電力系統(tǒng)潮流計(jì)算 電力系統(tǒng)潮流計(jì)算是對復(fù)雜電力系統(tǒng)正常和故障條件下穩(wěn)態(tài)運(yùn)行狀態(tài)的計(jì)算。潮流計(jì)算的目標(biāo)是求取電力系統(tǒng)在給定運(yùn)行方式下的節(jié)點(diǎn)電壓和功率分布，用以檢查系統(tǒng)各元件是否過負(fù)荷、各點(diǎn)電壓是否滿足要求、功率的分布和分配是否合理以及功率損耗等。對現(xiàn)有電力系統(tǒng)的運(yùn)行和擴(kuò)建，

6、對新的電力系統(tǒng)進(jìn)行規(guī)劃設(shè)計(jì)都是以潮流計(jì)算為基礎(chǔ)。 潮流計(jì)算結(jié)果的用途，例如用于電力系統(tǒng)穩(wěn)定研究、安全估計(jì)或最優(yōu)潮流等也對潮流計(jì)算的模型和方法有直接影響。 節(jié)點(diǎn)類型：（1）PV節(jié)點(diǎn)：柱入有功功率P為給定值，電壓也保持在給定數(shù)值。（2）PQ節(jié)點(diǎn)：諸如有功功率和無功功率是給定的。（3）平衡節(jié)點(diǎn)：用來平衡全電網(wǎng)的功率。選一容量足夠大的發(fā)電機(jī)擔(dān)任平衡全電網(wǎng)功率的職責(zé)。平衡節(jié)點(diǎn)的電壓大小與相位是給定的，通常以它的相角為參考量，即取其電壓相角為0。一個獨(dú)立的電力網(wǎng)中只設(shè)一個平衡點(diǎn)。 基本步驟：（1）形成節(jié)點(diǎn)導(dǎo)納矩陣；（2）將各節(jié)點(diǎn)電壓設(shè)初值U；（3）將節(jié)點(diǎn)初值代入相關(guān)求式，求出修正方程式的常數(shù)項(xiàng)向量；（4）

7、將節(jié)點(diǎn)電壓初值代入求式，求出雅可比矩陣元素；（5）求解修正方程，求修正向量；（6）求取節(jié)點(diǎn)電壓的新值；（7）檢查是否收斂，如不收斂，則以各節(jié)點(diǎn)電壓的新值作為初值自第3步重新開始進(jìn)行狹義次迭代，否則轉(zhuǎn)入下一步；（8）計(jì)算支路功率分布，PV節(jié)點(diǎn)無功功率和平衡節(jié)點(diǎn)柱入功率。2.3 MATLAB概述 目前電子計(jì)算機(jī)已廣泛應(yīng)用于電力系統(tǒng)的分析計(jì)算，潮流計(jì)算是其基本應(yīng)用之一?，F(xiàn)有很多潮流計(jì)算方法。對潮流計(jì)算方法有五方面的要求：（1）計(jì)算速度快；（2）內(nèi)存需要少；（3）計(jì)算結(jié)果有良好的可靠性和可信性；（4）適應(yīng)性好，亦即能處理變壓器變比調(diào)整、系統(tǒng)元件的不同描述和與其它程序配合的能力強(qiáng)；（5）簡單。 MATL

8、AB是一種交互式、面向?qū)ο蟮某绦蛟O(shè)計(jì)語言，廣泛應(yīng)用于工業(yè)界與學(xué)術(shù)界，主要用于矩陣運(yùn)算，同時在數(shù)值分析、自動控制模擬、數(shù)字信號處理、動態(tài)分析、繪圖等方面也具有強(qiáng)大的功能。 MATLAB程序設(shè)計(jì)語言結(jié)構(gòu)完整，且具有優(yōu)良的移植性，它的基本數(shù)據(jù)元素是不需要定義的數(shù)組。它可以高效率地解決工業(yè)計(jì)算問題，特別是關(guān)于矩陣和矢量的計(jì)算。MATLAB與C語言和FORTRAN語言相比更容易被掌握。通過M語言，可以用類似數(shù)學(xué)公式的方式來編寫算法，大大降低了程序所需的難度并節(jié)省了時間，從而可把主要的精力集中在算法的構(gòu)思而不是編程上。2.4牛頓-拉夫遜法原理假設(shè)有n個聯(lián)立的非線性代數(shù)方程： 假設(shè)以給出各變量的初值，令其分

9、別為個變量的修正量，使?jié)M足以上方程，所以： 將上式中的n個多元函數(shù)在初始值附近分別展開成泰勒級數(shù)，并略去含有，的二次及以上階次的各項(xiàng)，便得：方程可寫成： 以上方程是對于修正量，的線性方程組，稱為牛頓法的修正方程，可解出，。對初始近似解進(jìn)行修正： （i=1,2,，n）反復(fù)迭代，在進(jìn)行k+1次迭代時，從求解修正方程式：得到修正量,對各量進(jìn)行修正 （i=1,2,，n）迭代過程一直進(jìn)行到滿足收斂判據(jù)2.5牛頓-拉夫遜法解決潮流計(jì)算問題節(jié)點(diǎn)總數(shù)為n；PQ節(jié)點(diǎn)有m，；PV節(jié)點(diǎn)有n-m-1，平衡節(jié)點(diǎn)有1個，節(jié)點(diǎn)編號按照先PQ節(jié)點(diǎn)，再PV節(jié)點(diǎn)，最后平衡節(jié)點(diǎn)的順序進(jìn)行編號，即:1，2，m為PQ節(jié)點(diǎn)；m+1，m+

10、2，n-1為PV節(jié)點(diǎn)；n為平衡節(jié)點(diǎn)?？尚纬山Y(jié)點(diǎn)導(dǎo)納矩陣。導(dǎo)納矩陣元素可表示為，本文中節(jié)點(diǎn)電壓以直角坐標(biāo)形式表示，即。由此下列公式可求出Pi，Qi假設(shè)系統(tǒng)中的第1，2，m號節(jié)點(diǎn)為PQ節(jié)點(diǎn)，第i個節(jié)點(diǎn)的給定功率為和，對該節(jié)點(diǎn)可列方程： 假設(shè)系統(tǒng)中的第m+1，m+2，n-1號節(jié)點(diǎn)為PV節(jié)點(diǎn)，則對其中每一個節(jié)點(diǎn)可列方程: 第n號節(jié)點(diǎn)為平衡節(jié)點(diǎn)，其電壓為是給定的，故不參加迭代。 修正方程可寫成分塊矩陣的形式： 通過反復(fù)求解修正方程，解出各節(jié)點(diǎn)的未知量，再通過收斂判據(jù)判定是否已為真值。從而求得PQ節(jié)點(diǎn)的電壓V及相角的真值，PV節(jié)點(diǎn)的Q、真值，平衡節(jié)點(diǎn)的P、Q真值，以上即為牛頓-拉夫遜迭代法的潮流計(jì)算過程，

11、其優(yōu)點(diǎn)為計(jì)算精確，運(yùn)行速度快。其中的各個環(huán)節(jié)都可通過MATLAB程序來實(shí)現(xiàn)。2.6計(jì)算機(jī)潮流計(jì)算的步驟（1）對電力網(wǎng)絡(luò)的所有參數(shù)設(shè)初值,包括電壓、相角、有功、無功等。（2）處理非標(biāo)準(zhǔn)變比支路,使其變成標(biāo)準(zhǔn)變比為1的變壓器支路。（3）形成節(jié)點(diǎn)導(dǎo)納矩陣Y。（4）計(jì)算有功功率的不平衡量Pi,從而求出。（5）根據(jù)節(jié)點(diǎn)的類型形成J。（6）解修正方程式,求各節(jié)點(diǎn)的電壓的變化量ei(i=1,2,3.n,is)（7）求各節(jié)點(diǎn)相角的新值ei=ei+ei (i=1,2,3. n,is)（8）計(jì)算無功功率的不平衡量Qi,從而求出 (i=1,2,3.n,is)（9）解修正方程式,求各節(jié)點(diǎn)的電壓大小的變化量 (i=1,

12、2,3.,n,is)。（10）求各節(jié)點(diǎn)的電壓大小的新值 (i=1,2,3.,n,is)。（11）運(yùn)用個節(jié)點(diǎn)的電壓的新值自第四步開始下一次迭代。 計(jì)算平衡節(jié)點(diǎn)的功率和線路功率。其中平衡節(jié)點(diǎn)的功率的計(jì)算公式為線路上的功率為: 從而線路上的損耗的功率為:2.7計(jì)算機(jī)潮流計(jì)算流程圖選用牛頓-拉夫遜法進(jìn)行計(jì)算機(jī)潮流計(jì)算。牛頓-拉夫遜法潮流計(jì)算流程圖如圖2-1 。圖2-1 牛頓-拉夫遜法潮流計(jì)算流程圖3 手算潮流計(jì)算3.1確定節(jié)點(diǎn)類型節(jié)點(diǎn)為平衡節(jié)點(diǎn)節(jié)點(diǎn)為PV節(jié)點(diǎn)，其他四個節(jié)點(diǎn)都是PQ節(jié)點(diǎn)。3.2求節(jié)點(diǎn)導(dǎo)納矩陣Yb 由圖1可知，該系統(tǒng)以串聯(lián)支路的導(dǎo)納標(biāo)幺值和對地并聯(lián)導(dǎo)納標(biāo)幺值得等值電路如圖2所示?？傻脠D2的

13、等值節(jié)點(diǎn)導(dǎo)納矩陣。對角線上的元素為：非對角線上的元素為：所以節(jié)點(diǎn)導(dǎo)納矩陣為：Y=3.3計(jì)算各節(jié)點(diǎn)功率的修正方程的初始值（不平衡量）設(shè)各節(jié)點(diǎn)電壓初始值為：根據(jù)上述N-R法的求解過程，將各節(jié)點(diǎn)電壓的初始值代入（1）式和（2）式，進(jìn)行潮流計(jì)算，得節(jié)點(diǎn)功率和節(jié)點(diǎn)電壓：雅可比矩陣：J=修正各節(jié)點(diǎn)電壓：4 計(jì)算機(jī)算法潮流計(jì)算4.1計(jì)算機(jī)法潮流計(jì)算過程導(dǎo)納矩陣 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 1

14、5.0311-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.4686i功率方程第(1)次差值： Columns 1 through 12 0 0 4.2619 -2.1000 0.1000 -1.8000 -0.5500 -1.6000 8.4738 -3.7000 0 5.0000

15、形成的第(1)次Jacobi矩陣： Columns 1 through 11 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -37.3887 14.8252 5.9172 -14.2012 0 0 3.9002 -0.6240 0 0 0 -14.8252 -47.9126 14.2012 5.9172 0 0 0.6240 3.9002 0 0 0 5.9172 -14.2012 -8.0292 15.0311 3.1120 -0.8299 0 0 0 0 0 14.2012 5.9172 -15.0311 -9.0292 0.8299 3.

16、1120 0 0 0 0 0 0 0 3.1120 -0.8299 -5.2535 1.5846 2.6415 -0.7547 0 0 0 0 0 0.8299 3.1120 -1.5846 -5.7535 0.7547 2.6415 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 0 0 0 0 0 0 0 0 0 0 -2.4000 0 0 0 0 0 0 0 0 0 76.1905 0 Columns

17、 12 through 13 0 0 0 0 0 4.2619 0 -2.1000 0 0.1000 0 -1.8000 0 -0.5500 0 -1.6000 0 8.4738 63.4921 -3.7000 0 0 -63.4921 5.0000Jacobi矩陣第(1)次回代運(yùn)算 Columns 1 through 12 0 0 -0.1490 0.1103 -0.0024 0.1657 0.1015 0.3342 -0.1435 0.0297 0 -0.0431各個節(jié)點(diǎn)電壓模 0 1.1543 1.0160 0.9586 1.1439 1.2008功率方程第(2)次差值： Columns

18、 1 through 12 0 0 -0.9196 -0.3298 -0.2847 -0.0127 -0.5652 0.0384 -1.5185 0.1960 -0.0019 -0.3928形成的第(2)次Jacobi矩陣： Columns 1 through 11 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -47.1567 23.2568 5.2333 -16.9700 0 0 4.4126 -1.1470 0 0 0 -20.2170 -47.5883 16.9700 5.2333 0 0 1.1470 4.4126 0 0 0 3.5

19、786 -15.2150 -5.6605 18.1967 2.9819 -1.3474 0 0 0 0 0 15.2150 3.5786 -14.7623 -6.4582 1.3474 2.9819 0 0 0 0 0 0 0 2.5188 -1.7857 -3.5899 4.7795 2.1211 -1.5609 0 0 0 0 0 1.7857 2.5188 -1.7466 -5.2407 1.5609 2.1211 0 0 0 4.4413 -0.8294 0 0 2.9982 -0.9415 -83.4895 7.1529 72.6032 0 0 0.8294 4.4413 0 0 0

20、.9415 2.9982 -0.3333 -83.2844 1.8854 0 0 0 0 0 0 0 0 0 0 -2.4000 0 0 0 0 0 0 0 0 -2.7375 76.1905 -1.8854 Columns 12 through 13 0 0 0 0 0 -0.9196 0 -0.3298 0 -0.2847 0 -0.0127 0 -0.5652 0 0.0384 -1.8854 -1.5185 72.6032 0.1960 -0.0862 -0.0019 -72.6032 -0.3928Jacobi矩陣第(2)次回代運(yùn)算 Columns 1 through 12 0 0

21、0.0441 0.0155 0.0700 -0.0210 0.1790 -0.0275 0.0288 0.0422 -0.001 0.0487各個節(jié)點(diǎn)電壓模 0 1.1120 0.9435 0.7821 1.1170 1.2010功率方程第(3)次差值： Columns 1 through 12 0 0 -0.0417 0.0017 -0.0185 -0.0126 -0.1340 -0.0217 -0.0367 0.0646 -0.0024 -0.0916形成的第(3)次Jacobi矩陣： Columns 1 through 11 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

22、 0 0 0 0 0 0 0 0 -44.1917 23.5230 4.7526 -16.4350 0 0 4.2309 -1.1798 0 0 0 -19.9621 -46.3315 16.4350 4.7526 0 0 1.1798 4.2309 0 0 0 3.4621 -14.0963 -5.0872 17.0579 2.7814 -1.2240 0 0 0 0 0 14.0963 3.4621 -13.4379 -6.4674 1.2240 2.7814 0 0 0 0 0 0 0 1.9845 -1.5515 -1.8991 4.3504 1.6691 -1.3531 0 0 0

23、0 0 1.5515 1.9845 -1.3054 -5.0484 1.3531 1.6691 0 0 0 4.3026 -0.9762 0 0 2.8902 -1.0313 -79.8824 10.0760 70.7751 0 0 0.9762 4.3026 0 0 1.0313 2.8902 -3.4952 -82.5737 4.5677 0 0 0 0 0 0 0 0 0 0 -2.4019 0 0 0 0 0 0 0 0 0.3535 76.2524 -4.5677 Columns 12 through 13 0 0 0 0 0 -0.0417 0 0.0017 0 -0.0185 0

24、 -0.0126 0 -0.1340 0 -0.0217 -4.5677 -0.0367 70.7751 0.0646 0.0111 -0.0024 -70.7751 -0.0916Jacobi矩陣第(3)次回代運(yùn)算 Columns 1 through 12 0 0 0.0078 0.0043 0.0191 -0.0063 0.0768 -0.0055 0.0056 0.0156 0.0011 0.0181各個節(jié)點(diǎn)電壓模 0 1.1048 0.9237 0.7098 1.1126 1.2001功率方程第(4)次差值： Columns 1 through 12 0 0 -0.0004 0.000

25、2 -0.0002 -0.0012 -0.0276 -0.0058 0.0006 0.0016 -0.0003 -0.0053形成的第(4)次Jacobi矩陣： Columns 1 through 11 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -43.7449 23.5927 4.6455 -16.3504 0 0 4.1980 -1.1918 0 0 0 -20.0302 -45.9892 16.3504 4.6455 0 0 1.1918 4.1980 0 0 0 3.4389 -13.7881 -4.9899 16.7680 2.7

26、273 -1.1885 0 0 0 0 0 13.7881 3.4389 -13.0467 -6.4292 1.1885 2.7273 0 0 0 0 0 0 0 1.7502 -1.4706 -1.1217 4.2480 1.4705 -1.2806 0 0 0 0 0 1.4706 1.7502 -1.1039 -4.9984 1.2806 1.4705 0 0 0 4.2712 -1.0337 0 0 2.8637 -1.0684 -79.3737 11.1438 70.4222 0 0 1.0337 4.2712 0 0 1.0684 2.8637 -4.6945 -82.2282 5

27、.5608 0 0 0 0 0 0 0 0 0 0 -2.3998 0 0 0 0 0 0 0 0 1.5031 76.1844 -5.5608 Columns 12 through 13 0 0 0 0 0 -0.0004 0 0.0002 0 -0.0002 0 -0.0012 0 -0.0276 0 -0.0058 -5.5608 0.0006 70.4222 0.0016 0.0473 -0.0003 -70.4222 -0.0053Jacobi矩陣第(4)次回代運(yùn)算 Columns 1 through 12 0 0 0.0021 0.0014 0.0056 -0.0020 0.025

28、9 -0.0018 0.0016 0.0048 0.0002 0.0053各個節(jié)點(diǎn)電壓模 0 1.1029 0.9178 0.6857 1.1114 1.2000功率方程第(5)次差值： Columns 1 through 12 0 0 -0.0000 0.0000 0.0000 -0.0001 -0.0032 -0.0007 -0.0000 0.0001 -0.0000 -0.0005形成的第(5)次Jacobi矩陣： Columns 1 through 11 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -43.6300 23.6210 4

29、.6137 -16.3286 0 0 4.1889 -1.1958 0 0 0 -20.0561 -45.8841 16.3286 4.6137 0 0 1.1958 4.1889 0 0 0 3.4337 -13.6963 -4.9688 16.6803 2.7114 -1.1776 0 0 0 0 0 13.6963 3.4337 -12.9312 -6.4137 1.1776 2.7114 0 0 0 0 0 0 0 1.6711 -1.4437 -0.8564 4.2162 1.4034 -1.2564 0 0 0 0 0 1.4437 1.6711 -1.0344 -4.9843 1

30、.2564 1.4034 0 0 0 4.2621 -1.0517 0 0 2.8559 -1.0801 -79.2374 11.4911 70.3226 0 0 1.0517 4.2621 0 0 1.0801 2.8559 -5.0505 -82.1223 5.8687 0 0 0 0 0 0 0 0 0 0 -2.3993 0 0 0 0 0 0 0 0 1.8419 76.1690 -5.8687Columns 12 through 13 0 0 0 0 0 -0.0000 0 0.0000 0 0.0000 0 -0.0001 0 -0.0032 0 -0.0007 -5.8687

31、-0.0000 70.3226 0.0001 0.0580 -0.0000 -70.3226 -0.0005Jacobi矩陣第(5)次回代運(yùn)算 Columns 1 through 12 0 0 0.0003 0.0002 0.0008 -0.0003 0.0039 -0.0003 0.0002 0.0007 0.0000 0.0008各個節(jié)點(diǎn)電壓模 0 1.1026 0.9169 0.6820 1.1113 1.2000功率方程第(6)次差值： 1.0e-004 * Columns 1 through 12 0 0 -0.0026 0.0007 0.0032 -0.0182 -0.7419 -

32、0.1647 -0.0016 0.0184 -0.0063 -0.1036形成的第(6)次Jacobi矩陣： Columns 1 through 11 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -43.6128 23.6253 4.6089 -16.3254 0 0 4.1876 -1.1964 0 0 0 -20.0599 -45.8683 16.3254 4.6089 0 0 1.1964 4.1876 0 0 0 3.4330 -13.6825 -4.9658 16.6670 2.7090 -1.1760 0 0 0 0 0 13.6

33、825 3.4330 -12.9139 -6.4113 1.1760 2.7090 0 0 0 0 0 0 0 1.6591 -1.4396 -0.8161 4.2114 1.3932 -1.2528 0 0 0 0 0 1.4396 1.6591 -1.0238 -4.9822 1.2528 1.3932 0 0 0 4.2607 -1.0543 0 0 2.8548 -1.0818 -79.2170 11.5431 70.3078 0 0 1.0543 4.2607 0 0 1.0818 2.8548 -5.1036 -82.1067 5.9148 0 0 0 0 0 0 0 0 0 0

34、-2.3993 0 0 0 0 0 0 0 0 1.8924 76.1670 -5.9148 Columns 12 through 13 0 0 0 0 0 -0.0000 0 0.0000 0 0.0000 0 -0.0000 0 -0.0001 0 -0.0000 -5.9148 -0.0000 70.3078 0.0000 0.0596 -0.0000 -70.3078 -0.0000Jacobi矩陣第(6)次回代運(yùn)算 1.0e-004 * Columns 1 through 12 0 0 0.0764 0.0496 0.2048 -0.0728 0.9503 -0.0641 0.0564 0.1753 0.0074 0.1923各