(1)電氣08-1班穩(wěn)態(tài)設(shè)計(jì).doc

內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計(jì)說明書內(nèi)蒙古科技大學(xué)本科生課程設(shè)計(jì)說明書題 目：兩機(jī)五節(jié)點(diǎn)網(wǎng)絡(luò)潮流計(jì)算牛拉法學(xué)生姓名：李永清學(xué) 號(hào)：0867130116專 業(yè)：電氣工程及其自動(dòng)化班 級(jí)：電氣20081班指導(dǎo)教師：楊培宏，劉景霞內(nèi)蒙古科技大學(xué)課程設(shè)計(jì)任務(wù)書課程名稱電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計(jì)設(shè)計(jì)題目?jī)蓹C(jī)五節(jié)點(diǎn)網(wǎng)絡(luò)潮流計(jì)算牛拉法指導(dǎo)教師劉景霞時(shí)間1周一、教學(xué)要求電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計(jì)以設(shè)計(jì)和優(yōu)化電力系統(tǒng)的潮流分析為重點(diǎn)，提高學(xué)生綜合能力為目標(biāo)，盡可能結(jié)合實(shí)際工程進(jìn)行。設(shè)計(jì)內(nèi)容的安排要充分考慮學(xué)校現(xiàn)有的設(shè)備，設(shè)計(jì)時(shí)間及工程實(shí)際需要，并使學(xué)生初步學(xué)會(huì)運(yùn)用所學(xué)知識(shí)解決工程中的實(shí)際問題。二、設(shè)計(jì)資料及參數(shù)（一）設(shè)計(jì)原始資料1、待設(shè)計(jì)電氣設(shè)備系統(tǒng)圖2、電力系統(tǒng)網(wǎng)絡(luò)各元件參數(shù)3、電力系統(tǒng)電氣元件的使用規(guī)范4、電力工程電氣設(shè)計(jì)手冊(cè)（二）設(shè)計(jì)參考資料1、電力系統(tǒng)穩(wěn)態(tài)分析，陳珩，中國(guó)電力出版社，2007，第三版2、電力系統(tǒng)分析，韓禎祥，浙江大學(xué)出版社，2005，第三版3、電力系統(tǒng)分析課程實(shí)際設(shè)計(jì)與綜合實(shí)驗(yàn)，祝書萍，中國(guó)電力出版社，2007，第一版三、設(shè)計(jì)要求及成果1.根據(jù)給定的參數(shù)或工程具體要求，收集和查閱資料；學(xué)習(xí)相關(guān)軟件（軟件自選）。2.在給定的電力網(wǎng)絡(luò)上畫出等值電路圖。3.運(yùn)用計(jì)算機(jī)進(jìn)行潮流計(jì)算。4.編寫設(shè)計(jì)說明書。基本要求：1.編寫潮流計(jì)算程序；2.在計(jì)算機(jī)上調(diào)試通過（？）；3.運(yùn)行程序并計(jì)算出正確結(jié)果（？）；4.寫出課程設(shè)計(jì)報(bào)告（包括以下內(nèi)容）（1份）（1）程序框圖；（2）源程序；（3）符號(hào)說明表；（4）算例及計(jì)算結(jié)果5.編寫計(jì)算說明書（1份）。四、進(jìn)度安排根據(jù)給定的參數(shù)或工程具體要求，收集和查閱資料（半天）學(xué)習(xí)軟件（MATLAB或C語言等）（一天半）編程計(jì)算復(fù)雜系統(tǒng)潮流計(jì)算（三天）編寫計(jì)算設(shè)計(jì)書（一天）五、評(píng)分標(biāo)準(zhǔn)課程設(shè)計(jì)成績(jī)?cè)u(píng)定依據(jù)包括以下幾點(diǎn)：1） 工作態(tài)度（占10%）；2） 基本技能的掌握程度（占20%）；3） 程序編寫是否合理是否有運(yùn)行結(jié)果(40%)；4） 課程設(shè)計(jì)說明書編寫水平(占30%)。5） 分為優(yōu)、良、中、合格、不合格五個(gè)等級(jí)?？己朔绞剑涸O(shè)計(jì)期間教師現(xiàn)場(chǎng)檢查；評(píng)閱設(shè)計(jì)報(bào)告系統(tǒng)接線圖其中節(jié)點(diǎn)1為平衡節(jié)點(diǎn)，保持=1.06+j0為定值：節(jié)點(diǎn)2、3、4、5為PQ節(jié)點(diǎn)，給定的注入功率如圖。36內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計(jì)說明書兩機(jī)五節(jié)點(diǎn)網(wǎng)絡(luò)潮流計(jì)算牛拉法摘 要潮流計(jì)算是電力網(wǎng)絡(luò)設(shè)計(jì)及運(yùn)行中最基本的計(jì)算，對(duì)電力網(wǎng)絡(luò)的各種設(shè)計(jì)方案及各種運(yùn)行方式進(jìn)行潮流計(jì)算，可以得到各種電網(wǎng)各節(jié)點(diǎn)的電壓，并求得網(wǎng)絡(luò)的潮流及網(wǎng)絡(luò)中各元件的電力損耗，進(jìn)而求得電能損耗。傳統(tǒng)的潮流計(jì)算程序缺乏圖形用戶界面，結(jié)果顯示不直觀，難于與其他分析功能集成，網(wǎng)絡(luò)原始數(shù)據(jù)輸入工作量大且易于出錯(cuò)。結(jié)合電力系統(tǒng)的特點(diǎn)，對(duì)于復(fù)雜電力系統(tǒng)，根據(jù)定條件，應(yīng)用牛頓-拉夫遜法進(jìn)行計(jì)算，在手工計(jì)算中，由于涉及大量變量、微分方程、矩陣計(jì)算，求解很繁瑣，計(jì)算不同系統(tǒng)時(shí)需要重新計(jì)算。運(yùn)用 MATLAB軟件進(jìn)行仿真潮流計(jì)算，圖形界面直觀，運(yùn)行穩(wěn)定，計(jì)算準(zhǔn)確，提高了計(jì)算速度，各個(gè)類的有效封裝又使程序具有很好的模塊性.可維護(hù)性和可重用性。關(guān)鍵詞：電力系統(tǒng)潮流分析；牛頓-拉夫遜法潮流分析；Matlab仿真；內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計(jì)說明書兩機(jī)五節(jié)點(diǎn)網(wǎng)絡(luò)潮流計(jì)算牛拉法AbstractFlow calculation is the electric power network design and operation of the most basic calculation. For power network of all sorts of design scheme and various operation flow calculation method. Can get different grid of each node voltage. And get network trends and network of each element of power consumption. And getting electricity consumption. The traditional flow calculation program lack of graphical user interface. Results showed that dont intuitive, difficult to with other analysis function integration, network raw data input big workload and easy to go wrong. Combined with the characteristics of electric power system, for complex electric power systems, according to fixed conditions, application Newton - rafer Johnson method in calculating, manual calculation, because it involves a large variables, ordinary differential equations, calculation, solving very tedious, calculating different system needs to be calculated. Using software MATLAB simulation flow calculation, graphic interface intuitive, stable operation, computational accuracy, improves the calculation speed, each kind of effective encapsulation and makes the program has very good modularity. Maintainability and reusability.Key words: Electric power system flow caculation; Newton - Raphson power flow calculation; Matlab simulation內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計(jì)說明書目 錄摘 要IAbstractII第一章 前言11.1 潮流計(jì)算11.1.1潮流計(jì)算概述11.1.2潮流計(jì)算的要求21.1.3 潮流計(jì)算的優(yōu)勢(shì)31.2 MATLAB簡(jiǎn)介31.2.1 MATLAB概述31.2.2 MATLAB語言4第二章 牛頓-拉夫遜法概述32.1 牛頓-拉夫遜基本原理92.2 直角牛頓-拉夫遜法潮流計(jì)算求解過程12第三章 設(shè)計(jì)過程173.1設(shè)計(jì)資料及參數(shù)173.2牛頓拉夫遜法程序流程圖223.3設(shè)計(jì)程序及運(yùn)行結(jié)果23參考文獻(xiàn)27致謝36內(nèi)蒙古科技大學(xué)電力系統(tǒng)穩(wěn)態(tài)分析課程設(shè)計(jì)說明書第一章 前言1.1 潮流計(jì)算1.1.1 潮流計(jì)算概述潮流計(jì)算是研究電力系統(tǒng)穩(wěn)態(tài)運(yùn)行情況的一種基本電氣計(jì)算，常規(guī)潮流計(jì)算的任務(wù)是根據(jù)給定的運(yùn)行條件和網(wǎng)路結(jié)構(gòu)確定整個(gè)系統(tǒng)的運(yùn)行狀態(tài)，如各母線上的電壓（幅值及相角）、網(wǎng)絡(luò)中的功率分布以及功率損耗等。潮流計(jì)算的結(jié)果是電力系統(tǒng)穩(wěn)定計(jì)算和故障分析的基礎(chǔ)。通過潮流計(jì)算可以判斷電網(wǎng)母線電壓、支路電流和功率是否越限，如果有越限，就應(yīng)采取措施，調(diào)整運(yùn)行方式。對(duì)于正在規(guī)劃的電力系統(tǒng)，通過潮流計(jì)算，可以為選擇電網(wǎng)供電方案和電氣設(shè)備提供依據(jù)。潮流計(jì)算還可以為繼電保護(hù)和自動(dòng)裝置定整計(jì)算、電力系統(tǒng)故障計(jì)算和穩(wěn)定計(jì)算等提供原始數(shù)據(jù)。具體表現(xiàn)在以下方面：(1)在電網(wǎng)規(guī)劃階段,通過潮流計(jì)算,合理規(guī)劃電源容量及接入點(diǎn),合理規(guī)劃網(wǎng)架,選擇無功補(bǔ)償方案,滿足規(guī)劃水平的大、小方式下潮流交換控制、調(diào)峰、調(diào)相、調(diào)壓的要求。(2)在編制年運(yùn)行方式時(shí),在預(yù)計(jì)負(fù)荷增長(zhǎng)及新設(shè)備投運(yùn)基礎(chǔ)上,選擇典型方式進(jìn)行潮流計(jì)算,發(fā)現(xiàn)電網(wǎng)中薄弱環(huán)節(jié),供調(diào)度員日常調(diào)度控制參考,并對(duì)規(guī)劃、基建部門提出改進(jìn)網(wǎng)架結(jié)構(gòu),加快基建進(jìn)度的建議。(3)正常檢修及特殊運(yùn)行方式下的潮流計(jì)算,用于日運(yùn)行方式的編制,指導(dǎo)發(fā)電廠開機(jī)方式,有功、無功調(diào)整方案及負(fù)荷調(diào)整方案,滿足線路、變壓器熱穩(wěn)定要求及電壓質(zhì)量要求。(4)預(yù)想事故、設(shè)備退出運(yùn)行對(duì)靜態(tài)安全的影響分析及作出預(yù)想的運(yùn)行方式調(diào)整方案?？偨Y(jié)為在電力系統(tǒng)運(yùn)行方式和規(guī)劃方案的研究中，都需要進(jìn)行潮流計(jì)算以比較運(yùn)行方式或規(guī)劃供電方案的可行性、可靠性和經(jīng)濟(jì)性。同時(shí)，為了實(shí)時(shí)監(jiān)控電力系統(tǒng)的運(yùn)行狀態(tài)，也需要進(jìn)行大量而快速的潮流計(jì)算。因此，潮流計(jì)算是電力系統(tǒng)中應(yīng)用最廣泛、最基本和最重要的一種電氣運(yùn)算。在系統(tǒng)規(guī)劃設(shè)計(jì)和安排系統(tǒng)的運(yùn)行方式時(shí)，采用離線潮流計(jì)算；在電力系統(tǒng)運(yùn)行狀態(tài)的實(shí)時(shí)監(jiān)控中，則采用在線潮流計(jì)算。此外，電力系統(tǒng)潮流計(jì)算也是計(jì)算系統(tǒng)動(dòng)態(tài)穩(wěn)定和靜態(tài)穩(wěn)定的基礎(chǔ)。所以潮流計(jì)算是研究電力系統(tǒng)的一種很重要和基礎(chǔ)的計(jì)算。1.1.2潮流計(jì)算的要求電力系統(tǒng)運(yùn)行必須滿足一定技術(shù)和經(jīng)濟(jì)上的要求。這些要求夠成了潮流問題中某些變量的約束條件，常用的約束條件如下：1. 節(jié)點(diǎn)電壓應(yīng)滿足 式(1-1) 從保證電能質(zhì)量和供電安全的要求來看，電力系統(tǒng)的所有電氣設(shè)備都必須運(yùn)行在額定電壓附近。PU節(jié)點(diǎn)電壓幅值必須按上述條件給定。因此，這一約束條件對(duì)PQ節(jié)點(diǎn)而言。2. 節(jié)點(diǎn)的有功功率和無功功率應(yīng)滿足 式(1-2) PQ節(jié)點(diǎn)的有功功率和無功功率，以及PU節(jié)點(diǎn)的有功功率，在給定是就必須滿足上述條件，因此，對(duì)平衡節(jié)點(diǎn)的P和Q以及PU節(jié)點(diǎn)的Q應(yīng)按上述條件進(jìn)行檢驗(yàn)。3. 節(jié)點(diǎn)之間電壓的相位差應(yīng)滿足 式(1-3) 為了保證系統(tǒng)運(yùn)行的穩(wěn)定性，要求某些輸電線路兩端的電壓相位不超過一定的數(shù)值。這一約束的主要意義就在于此。 因此，潮流計(jì)算可以歸結(jié)為求解一組非線性方程組，并使其解答滿足一定的約束條件。常用的方法是迭代法和牛頓法，在計(jì)算過程中，或得出結(jié)果之后用約束條件進(jìn)行檢驗(yàn)。如果不能滿足要求，則應(yīng)修改某些變量的給定值，甚至修改系統(tǒng)的運(yùn)行方式，重新進(jìn)行計(jì)算。1.1.3潮流計(jì)算的優(yōu)勢(shì)電力系統(tǒng)潮流計(jì)算是電力系統(tǒng)分析中的一種最基本的計(jì)算，是對(duì)復(fù)雜電力系統(tǒng)正常和故障條件下穩(wěn)態(tài)運(yùn)行狀態(tài)的計(jì)算。潮流計(jì)算的目標(biāo)是求取電力系統(tǒng)在給定運(yùn)行狀態(tài)的計(jì)算。即節(jié)點(diǎn)電壓和功率分布，用以檢查系統(tǒng)各元件是否過負(fù)荷。各點(diǎn)電壓是否滿足要求，功率的分布和分配是否合理以及功率損耗等。對(duì)現(xiàn)有電力系統(tǒng)的運(yùn)行和擴(kuò)建，對(duì)新的電力系統(tǒng)進(jìn)行規(guī)劃設(shè)計(jì)以及對(duì)電力系統(tǒng)進(jìn)行靜態(tài)和暫態(tài)穩(wěn)定分析都是以潮流計(jì)算為基礎(chǔ)。潮流計(jì)算結(jié)果可用于電力系統(tǒng)穩(wěn)態(tài)研究，安全估計(jì)或最優(yōu)潮流等，實(shí)際電力系統(tǒng)的潮流計(jì)算主要采用牛頓-拉夫遜法。借由MATLAB可以輕松實(shí)現(xiàn)計(jì)算復(fù)雜的電力系統(tǒng)潮流分布。1.2 MATLAB簡(jiǎn)介1.2.1 MATLAB概述MATLAB的含義是矩陣實(shí)驗(yàn)室（Matrix Laboratory），是美國(guó)MathWork公司于1982推出的一套高性能的數(shù)值計(jì)算可視化軟件，包括MATLAB主程序、SIMULINK 動(dòng)態(tài)系統(tǒng)仿真包和各種專業(yè)工具箱它集數(shù)值分析，矩陣計(jì)算，信號(hào)處理和圖形顯示于一體，構(gòu)成一個(gè)方便的，界面友好的用戶環(huán)境，具有極強(qiáng)大的計(jì)算功能和極高的編程效率，特別適合于科學(xué)計(jì)算、數(shù)值分析、系統(tǒng)仿真和信號(hào)處理等任務(wù)。MATLAB程學(xué)設(shè)計(jì)語言結(jié)構(gòu)完整，且具有優(yōu)良的移植性，它的基本數(shù)據(jù)元素是不需要定義的數(shù)組。它可以高效率的解決工業(yè)計(jì)算問題，特別是關(guān)于矩陣和矢量的計(jì)算。通過MATLAB語言，可以用類似數(shù)學(xué)公式的方式來編寫算法，大大降低了程序需要的難度別難并節(jié)省了時(shí)間，從而可把主要的經(jīng)歷集中在算法的構(gòu)思而不是編程上。學(xué)習(xí)運(yùn)用MATLAB計(jì)算電力系統(tǒng)潮流分布是本次課程設(shè)計(jì)的重點(diǎn)，可以說，作為工科學(xué)生，會(huì)運(yùn)用MATLAB來解決工程問題已成為必須。到目前為止，MATLAB已發(fā)展成為國(guó)際上最優(yōu)秀的科技應(yīng)用軟件之一。它的功能十分強(qiáng)大，不僅僅可以實(shí)現(xiàn)計(jì)算潮流分布，還可以模擬仿真各式各樣的數(shù)值系統(tǒng)，工程。這里將借助MATLAB來完成用直角牛頓-拉夫遜法計(jì)算電力系統(tǒng)潮流分布。1.2.2 MATLAB語言1、MATLAB語言相關(guān)指令表1-1管理命令和函數(shù)help在線幫助文件doc裝入超文本說明whatM、MAT、MEX文件的目錄列表type列出M文件lookfor通過help條目搜索關(guān)鍵字which定位函數(shù)和文件Demo運(yùn)行演示程序Path控制MATLAB的搜索路徑管理變量和工作空間Who列出當(dāng)前變量Whos列出當(dāng)前變量（長(zhǎng)表）Load從磁盤文件中恢復(fù)變量Save保存工作空間變量Clear從內(nèi)存中清除變量和函數(shù)Pack整理工作空間內(nèi)存Size矩陣的尺寸Length向量的長(zhǎng)度disp顯示矩陣或與文件和操作系統(tǒng)有關(guān)的命令cd改變當(dāng)前工作目錄Dir目錄列表Delete刪除文件Getenv獲取環(huán)境變量值!執(zhí)行DOS操作系統(tǒng)命令Unix執(zhí)行UNIX操作系統(tǒng)命令并返回結(jié)果Diary保存MATLAB任務(wù)控制命令窗口Cedit設(shè)置命令行編輯Clc清命令窗口Home光標(biāo)置左上角Format設(shè)置輸出格式Echo底稿文件內(nèi)使用的回顯命令more在命令窗口中控制分頁(yè)輸出啟動(dòng)和退出MATLABQuit退出MATLABStartup引用MATLAB時(shí)所執(zhí)行的M文件Matlabrc主啟動(dòng)M文件一般信息InfoMATLAB系統(tǒng)信息及Mathworks公司信息Subscribe成為MATLAB的訂購(gòu)用戶hostidMATLAB主服務(wù)程序的識(shí)別代號(hào)Whatsnew在說明書中未包含的新信息Ver版本信息操作符和特殊字符+加減*矩陣乘法.*數(shù)組乘法矩陣冪.數(shù)組冪左除或反斜杠/右除或斜杠./數(shù)組除KronKronecker張量積:冒號(hào)( )圓括號(hào) 方括號(hào).小數(shù)點(diǎn).父目錄繼續(xù),逗號(hào);分號(hào)%注釋!感嘆號(hào)轉(zhuǎn)置或引用=賦值= =相等關(guān)系操作符&邏輯與|邏輯或邏輯非xor邏輯異或邏輯函數(shù)Exist檢查變量或函數(shù)是否存在Any向量的任一元為真，則其值為真All向量的所有元為真，則其值為真Find找出非零元素的索引號(hào)三角函數(shù)Sin正弦Sinh雙曲正弦Asin反正弦Asinh0.00001 for m=1:N if m=isb pp=0;qq=0; for n=1:N pp=pp+ei(m)*(G(m,n)*ei(n)-B(m,n)*fi(n)+fi(m)*(G(m,n)*fi(n)+B(m,n)*ei(n); qq=qq+fi(m)*(G(m,n)*ei(n)-B(m,n)*fi(n)-ei(m)*(G(m,n)*fi(n)+B(m,n)*ei(n); end DP(m)=P(m)-pp; %有功功率的不平衡量 DQ(m)=Q(m)-qq; %無功功率的不平衡量 I(m)= (P(m)-Q(m)/conj(ei(m)+fi(m)*i); H(m,m)=-B(m,m)*ei(m)+G(m,m)*fi(m)+imag(I(m); N(m,m)=G(m,m)*ei(m)+B(m,m)*fi(m)+real(I(m); J(m,m)=-G(m,m)*ei(m)-B(m,m)*fi(m)+real(I(m); L(m,m)=-B(m,m)*ei(m)+G(m,m)*fi(m)-imag(I(m); JJ(2*m-3,2*m-3)=H(m,m);JJ(2*m-3,2*m-2)=N(m,m);JJ(2*m-2,2*m-3)=J(m,m);JJ(2*m-2,2*m-2)=L(m,m); end end for m=1:N if m=isb for n=1:N if n=isb if n=m else H(m,n)=-B(m,n)*ei(m)+G(m,n)*fi(m); N(m,n)=G(m,n)*ei(m)+B(m,n)*fi(m); J(m,n)=-B(m,n)*fi(m)-G(m,n)*ei(m); L(m,n)=G(m,n)*fi(m)-B(m,n)*ei(m); JJ(2*m-3,2*n-3)=H(m,n);JJ(2*m-3,2*n-2)=N(m,n);JJ(2*m-2,2*n-3)=J(m,n);JJ(2*m-2,2*n-2)=L(m,n); end end end end end disp(雅克比矩陣JJ:); disp(JJ); for m=1:N if m=isb DW(2*m-3)=DP(m);DW(2*m-2)=DQ(m); %形成不平衡量的矩陣 end end DY=JJDW; %解修正方程式 pr=max(abs(DY); disp(第M次修正方程的解DY:); disp(DY); for n=1:N if n=isb fi(n)=fi(n) +DY(2*n-3); ei(n)=ei(n) +DY(2*n-2); end %計(jì)算新值 end u=ei+fi*i; disp(節(jié)點(diǎn)電壓的第C(k)次近似值:); disp(u); disp(各點(diǎn)的電壓實(shí)部di(單位:V)為(節(jié)點(diǎn)號(hào)從小到大排列):); disp(ei); disp(各點(diǎn)的電壓虛部fi(單位:V)為(節(jié)點(diǎn)號(hào)從小到大排列):); disp(fi); k=k+1; disp(迭代次數(shù):); disp(k);end for m=1:N I(m)=Y(1,m)*u(m);enddisp(平衡節(jié)點(diǎn)的功率);S1=u(1)*sum(conj(I)%計(jì)算平衡節(jié)點(diǎn)的功率for m=1:N for n=1:N S(m,n)=u(m)*(conj(u(m)-conj(u(n)*conj(-Y(m,n);%計(jì)算各支路功率 endenddisp(各支路功率);disp(S) %結(jié)束運(yùn)行結(jié)果導(dǎo)納矩陣為:Y= 6.2500 -18.7500i -5.0000 +15.0000i -1.2500 + 3.7500i 0 0 -5.0000 +15.0000i 10.8333 -32.5000i -1.6667 + 5.0000i -1.6667 + 5.0000i -2.5000 + 7.5000i -1.2500 + 3.7500i -1.6667 + 5.0000i 12.9167 -38.7500i -10.0000 +30.0000i 0 0 -1.6667 + 5.0000i -10.0000 +30.0000i 12.9167 -38.7500i -1.2500 + 3.7500i 0 -2.5000 + 7.5000i 0 -1.2500 + 3.7500i 3.7500 -11.2500i雅克比矩陣JJ: 32.5000 10.8333 -5.0000 -1.6667 -5.0000 -1.6667 -7.5000 -2.5000 -10.8333 32.5000 1.6667 -5.0000 1.6667 -5.0000 2.5000 -7.5000 -5.0000 -1.6667 38.7500 13.2167 -30.0000 -10.0000 0 0 1.6667 -5.0000 -12.6167 38.7500 10.0000 -30.0000 0 0 -5.0000 -1.6667 -30.0000 -10.0000 38.7500 13.2667 -3.7500 -1.2500 1.6667 -5.0000 10.0000 -30.0000 -12.5667 38.7500 1.2500 -3.7500 -7.5000 -2.5000 0 0 -3.7500 -1.2500 11.2500 4.2500 2.5000 -7.5000 0 0 1.2500 -3.7500 -3.2500 11.2500第M次修正方程的解DY: 0.0473 0.0731 0.0847 0.0965 0.0901 0.0971 0.1043 0.0996節(jié)點(diǎn)電壓的第C(k)次近似值: 1.0600 1.4038 1.6891 1.7277 1.8295各點(diǎn)的電壓實(shí)部di(單位:V)為(節(jié)點(diǎn)號(hào)從小到大排列): 1.0600 1.0731 1.0965 1.0971 1.0996各點(diǎn)的電壓虛部fi(單位:V)為(節(jié)點(diǎn)號(hào)從小到大排列): 0 0.0473 0.0847 0.0901 0.1043迭代次數(shù): 1雅克比矩陣JJ: 35.3867 10.0891 -5.4441 -1.5522 -5.4441 -1.5522 -8.1662 -2.3283 -10.0891 35.3867 1.5522 -5.4441 1.5522 -5.4441 2.3283 -8.1662 -5.6238 -1.4043 43.5841 11.0607 -33.7425 -8.4256 0 0 1.4043 -5.6238 -10.7055 43.5841 8.4256 -33.7425 0 0 -5.6355 -1.3780 -33.8127 -8.2680 43.6748 10.8821 -4.2266 -1.0335 1.3780 -5.6355 8.2680 -33.8127 -10.4769 43.6748 1.0335 -4.2266 -8.5074 -1.9669 0 0 -4.2537 -0.9834 12.7612 3.2236 1.9669 -8.5074 0 0 0.9834 -4.2537 -2.6770 12.7612第M次修正方程的解DY: 0.0003 -0.0001 0.0006 -0.0017 0.0005 -0.0015 0.0005 -0.0013節(jié)點(diǎn)電壓的第C(k)次近似值: 1.0600 1.4055 1.6913 1.7299 1.8320各點(diǎn)的電壓實(shí)部di(單位:V)為(節(jié)點(diǎn)號(hào)從小到大排列): 1.0600 1.0729 1.0949 1.0955 1.0982各點(diǎn)的電壓虛部fi(單位:V)為(節(jié)點(diǎn)號(hào)從小到大排列): 0 0.0475 0.0852 0.0906 0.1048迭代次數(shù): 2雅克比矩陣JJ: 35.3849 10.0795 -5.4438 -1.5507 -5.4438 -1.5507 -8.1658 -2.3260 -10.0795 35.3849 1.5507 -5.4438 1.5507 -5.4438 2.3260 -8.1658 -5.6163 -1.3987 43.5261 11.0173 -33.6976 -8.3922 0 0 1.3987 -5.6163 -10.6626 43.5261 8.3922 -33.6976 0 0 -5.6288 -1.3728 -33.7725 -8.2366 43.6228 10.8412 -4.2216 -1.0296 1.3728 -5.6288 8.2366 -33.7725 -10.4366 43.6228 1.0296 -4.2216 -8.4987 -1.9594 0 0 -4.2494 -0.9797 12.7481 3.2120 1.9594 -8.