1、 PAGE47 / NUMPAGES47科技大學電力系統穩態分析課程設計題 目：基于MATLAB的電力系統復雜潮流分析學生：建峰學 號：1167130207專 業：電氣工程與其自動化班 級：電氣20112班指導教師：景霞摘要電力系統潮流計算是電網分析的基礎應用，是對復雜電力系統正常和故障條件下穩態運行狀態的計算。給定電力系統的網絡結構、參數和決定電力系統運行狀況的邊界條件，確定電力系統運行的方法之一是朝流計算。MATLAB是一種交互式、面向對象的程序設計語言，廣泛應用于工業界與學術界，主要用于矩陣運算采用迭代法，通過建立矩陣的修正方程來依次迭代，逐步逼近真值來計算出電力網的電壓，功率分布。PQ

2、分解法是極坐標形式牛頓-拉潮流計算的一種簡化計算方法，。PQ分解法通過對電力系統具體特點的分析，對牛頓法修正方程式的雅可比矩陣進行了有效的簡化和改進。由于這些簡化只涉與修正方程式的系數矩陣，并未改變節點功率平衡方程和收斂判據，因不會降低計算結果的精度。用手算和計算機算法對其進行設計。使用MATLAB軟件進行編程，在很大程度上節省了存，減少了計算量。通過對本題計算我們了解了一些工程計算和解決工程問題的方法。關鍵詞：潮流計算，PQ分解法，MATLAB Electrical power system complex tidal current analysis based on MATLABPowe

3、r Flow Analysis Grid computing is the basis of applications, the complex power system under normal and fault conditions for the calculation of steady state operation. Given the power system network structure, parameters and decisions operation of the power system boundary conditions, to determine th

4、e method of operation of the power system is one of North Korea flow calculation.MATLAB is an interactive, object-oriented programming language, widely used in industry and academia, mainly for matrix calculation. Using iteration, the amendment through the establishment of matrix iterative equation

5、to turn, gradually moving towards a true value to calculate the voltage electricity grid, power distribution.PQ decomposition method is the form of polar coordinates Newton - the widening trend of a simplified calculation method. P-Q decomposition method adopted by the specific characteristics of th

6、e power system analysis, Newtons Law of the Jacobian matrix formula has effectively simplified and improved. As a result of these simplified formula that involves only the coefficient matrix, the balance of power has not changed node equations and the convergence criterion, because the results will

7、not reduce the accuracy.Use MATLAB software programming, saving memory to a large extent, reduce the amount of computation. By this calculation we understand that a number of engineering calculation and solve engineering problemsKeywords : The trend, the PQ decomposition method, MATLAB目錄 TOC o 1-3 h

8、 z u HYPERLINK l _Toc376010692科技大學課程設計任務書 PAGEREF _Toc376010692 h 5HYPERLINK l _Toc376010693第一章引言 PAGEREF _Toc376010693 h 7HYPERLINK l _Toc3760106941.1 研究背景與意義 PAGEREF _Toc376010694 h 7HYPERLINK l _Toc3760106951.2潮流計算的意義 PAGEREF _Toc376010695 h 8HYPERLINK l _Toc3760106961.3電力系統穩態分析潮流計算總結 PAGEREF _To

9、c376010696 h 9HYPERLINK l _Toc3760106971.4 MATLAB的概述 PAGEREF _Toc376010697 h 9HYPERLINK l _Toc376010698第二章理論計算 PAGEREF _Toc376010698 h 11HYPERLINK l _Toc3760106992.1 P-Q法潮流計算的基本步驟 PAGEREF _Toc376010699 h 11HYPERLINK l _Toc3760107002.2 PQ分解法潮流計算流程圖 PAGEREF _Toc376010700 h 12HYPERLINK l _Toc3760107012

10、.3 兩機五節點網絡潮流計算 PAGEREF _Toc376010701 h 13HYPERLINK l _Toc376010702第三章程序設計 PAGEREF _Toc376010702 h 18HYPERLINK l _Toc3760107033.1 設計流程 PAGEREF _Toc376010703 h 18HYPERLINK l _Toc37601070432 程序設計 PAGEREF _Toc376010704 h 18HYPERLINK l _Toc3760107053.3 程序運行結果 PAGEREF _Toc376010705 h 22HYPERLINK l _Toc376

11、010706第四章設計感想 PAGEREF _Toc376010706 h 44HYPERLINK l _Toc376010707參考文獻 PAGEREF _Toc376010707 h 45科技大學課程設計任務書課程名稱電力系統穩態分析課程設計設計題目兩機五節點網絡潮流計算 P-Q法指導教師景霞時間1周一、教學要求電力系統穩態分析課程設計以設計和優化電力系統的潮流分析為重點，提高學生綜合能力為目標，盡可能結合實際工程進行。設計容的安排要充分考慮學校現有的設備，設計時間與工程實際需要，并使學生初步學會運用所學知識解決工程中的實際問題。二、設計資料與參數（一）設計原始資料1、待設計電氣設備系統圖

12、2、電力系統網絡各元件參數3、電力系統電氣元件的使用規4、電力工程電氣設計手冊（二）設計參考資料1、電力系統穩態分析，珩，中國電力，2007，第三版2、電力系統分析，禎祥，大學，2005，第三版3、電力系統分析課程實際設計與綜合實驗，祝書萍，中國電力，2007，第一版三、設計要求與成果1.根據給定的參數或工程具體要求，收集和查閱資料；學習相關軟件（軟件自選）。2.在給定的電力網絡上畫出等值電路圖。3.運用計算機進行潮流計算。4.編寫設計說明書。基本要求：1.編寫潮流計算程序；2.在計算機上調試通過（？）；3.運行程序并計算出正確結果（？）；4.寫出課程設計報告（包括以下容）（1份）（1）程序框

13、圖；（2）源程序；（3）符號說明表；（4）算例與計算結果5.編寫計算說明書（1份）。四、進度安排根據給定的參數或工程具體要求，收集和查閱資料（半天）學習軟件（MATLAB或C語言等）（一天半）編程計算復雜系統潮流計算（三天）編寫計算設計書（一天）五、評分標準課程設計成績評定依據包括以下幾點：工作態度（占10%）；基本技能的掌握程度（占20%）；程序編寫是否合理是否有運行結果(40%)；課程設計說明書編寫水平(占30%)。分為優、良、中、合格、不合格五個等級。考核方式：設計期間教師現場檢查；評閱設計報告第一章 引言1.1 研究背景與意義電力系統在運行時，在電源電勢激勵作用下，電流或功率從電源通過

14、系統各元件流入負荷，分布于電力網各處，稱為潮流分布。潮流計算是電力系統分析中的一種最基本的計算，即根據給定電力系統的網絡結構和參數，在滿足電力系統運行狀態參數的邊界條件情況下，確定電力系統穩態運行狀態的最基本的方法。它的任務是在給定的接線方式和運行條件下，確定系統的運行狀態，如各母線上的電壓（幅值和相角）、網絡中的功率分布與功率損耗等，是電力系統的穩態計算。潮流計算是對電力系統正常運行狀況的分析和計算，即電力系統中的電壓、電流、功率的計算，即潮流計算；潮流計算方法很多：高斯塞德爾法、牛頓拉夫遜法、P-Q分解法、直流潮流法，以與由高斯塞德爾法、牛頓拉夫遜法演變的各種潮流計算方法。潮流計算可以用傳

15、統的手工方式進行，也可以計算機為工具通過軟件完成。兩種方法各有優缺點。前者物理概念清晰，可用來計算一些接線較簡單的電力網，但若將其用于接線復雜的電力網則計算量過大，難于保證計算準確性。后者從數學上看可歸結為用數值方法解非線性代數方程，數學邏輯簡單完整，借助計算機可快速精確地完成計算，但其缺點是物理概念不明顯，物理規律被埋沒在循環往復的數值求解過程中。潮流計算是電力系統非常重要的分析計算，用以研究系統規劃和運行中提出的各種問題。對規劃中的電力系統，通過潮流計算可以檢驗所提出的電力系統規劃方案能否滿足各種運行方式的要求：對運行中的電力系統，通過潮流計算可以預知各種負荷變化和網絡結構的改變會不會危與

16、系統的安全，系統中所有母線的電壓是否在允許的圍以，系統中各元件（線路、變壓器等）是否會出現過負荷，以與可能出現過負荷時應事先采取哪些預防措施等。因此潮流計算的目的是： 為電力系統規劃設計提供接線、電氣設備選擇和導線截面選擇的依據。 提供電力線運行方式和制定檢修計劃的依據。 提供繼電保護、自動裝置設計和整定計算的依據。 為調壓計算、經濟運行計算、短路和穩定計算提供必要的數據。潮流計算的研究是從20世紀50年代隨著電網的產生而開始的，涌現出各種算法都是圍繞著算法的可靠性或收斂性；對計算機存量的要求和計算速度；計算的方便性和靈活性等。1.2潮流計算的意義潮流計算是研究電力系統穩態運行情況的一種基本電

17、氣計算，常規潮流計算的任務是根據給定的運行條件和網路結構確定整個系統的運行狀態，如各母線上的電壓（幅值與相角）、網絡中的功率分布以與功率損耗等。潮流計算的結果是電力系統穩定計算和故障分析的基礎。通過潮流計算可以判斷電網母線電壓、支路電流和功率是否越限，如果有越限，就應采取措施，調整運行方式。對于正在規劃的電力系統，通過潮流計算，可以為選擇電網供電方案和電氣設備提供依據。潮流計算還可以為繼電保護和自動裝置定整計算、電力系統故障計算和穩定計算等提供原始數據。具體表現在以下方面：(1)在電網規劃階段,通過潮流計算,合理規劃電源容量與接入點,合理規劃網架,選擇無功補償方案,滿足規劃水平的大、小方式下潮

18、流交換控制、調峰、調相、調壓的要求。(2)在編制年運行方式時,在預計負荷增長與新設備投運基礎上,選擇典型方式進行潮流計算,發現電網中薄弱環節,供調度員日常調度控制參考,并對規劃、基建部門提出改進網架結構,加快基建進度的建議。(3)正常檢修與特殊運行方式下的潮流計算,用于日運行方式的編制,指導發電廠開機方式,有功、無功調整方案與負荷調整方案,滿足線路、變壓器熱穩定要求與電壓質量要求。(4)預想事故、設備退出運行對靜態安全的影響分析與作出預想的運行方式調整方案。總結為在電力系統運行方式和規劃方案的研究中，都需要進行潮流計算以比較運行方式或規劃供電方案的可行性、可靠性和經濟性。同時，為了實時監控電力

19、系統的運行狀態，也需要進行大量而快速的潮流計算。因此，潮流計算是電力系統中應用最廣泛、最基本和最重要的一種電氣運算。在系統規劃設計和安排系統的運行方式時，采用離線潮流計算；在電力系統運行狀態的實時監控中，則采用在線潮流計算。1.3電力系統穩態分析潮流計算總結潮流計算針對電力革統各正常運行方式，而靜態安全分析則要研究各種運行方式下個別系統元件退出運行后系統的狀況。其目的是校驗系統是否能安全運行，即是否有過負荷的元件或電壓過低的母線等。原則上講，靜態安全分析也可U用潮流計算來代替。但是一般靜態安全分析需要校驗的狀態數非常多，用嚴格的潮流計算來分析這些狀態往往計算量過大，因此不得不尋求一些特殊的算法

20、以滿足要求。利用電子數字計算機進行電力系統潮流計算從20 世紀50 年代中期就己開始，此后，潮流計算曾采用了各種不同的方法，這些方法的發展主要是圍繞著對潮流計算的一些基本要求進行的，對潮流計算的要求可以歸納為下面幾點：(1)計算方法的可靠性或收斂性(2)對計算速度和存量的要求(3)計算的方便性和靈活性一般情況下，采用P-Q分解法計算時要求的迭代次數較采用牛頓-拉夫遜法時多，但每次迭代所需的時間則較采用牛頓-拉夫遜法時少，以致總的計算速度仍是P-Q分解法快。1.4 MATLAB的概述目前電子計算機已廣泛應用于電力系統的分析計算，潮流計算是其基本應用軟件之一。現有很多潮流計算方法。對潮流計算方法有

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

22、時間，從而可以把主要的精力集中在算法的構思而不是編程上。另外，MATLAB提供了一種特殊的工具：工具箱（TOOLBOXES）.這些工具箱主要包括：信號處理（SIGNAL PROCESSING）、控制系統（CONTROL SYSTEMS）、神經網絡（NEURAL NETWORKS）、模糊邏輯(FUZZY LOGIC)、小波(WAVELETS)和模擬（SIMULATION）等等。不同領域、不同層次的用戶通過相應工具的學習和應用，可以方便地進行計算、分析與設計工作。 MATLAB設計中，原始數據的填寫格式是很關鍵的一個環節，它與程序使用的方便性和靈活性有著直接的關系。原始數據輸入格式的設計，主要應從

23、使用的角度出發，原則是簡單明了，便于修改。第二章 理論計算2.1 P-Q法潮流計算的基本步驟運用P-Q分解法計算潮流分布時的基本步驟是：形成系數矩陣B、B，并求其逆陣。設各節點電壓的初值（i=1,2，n，is）和（i=1,2，m，is）。按式=-計算有功功率的不平衡量，從而求出（i=1，2，ns）。解修正方程式 ()=BU，求各節點電壓相位角的變量（i=1,2，n，is）。求各節點電壓相位角的新值=+（i=1,2，n，is）。按式=- 計算無功功率的不平衡量，從而求出（i=1,2，m，is）。解修正方程式()=BU（i=1,2，m，is）。求各節點電壓大小的新值=+（i=1,2，m，is）。運

24、用各節點電壓的新值自第三步開始進入下一次迭代。計算平衡節點功率和線路功率。2.2 PQ分解法潮流計算流程圖 啟動輸入原始數據形成節點導納矩陣形成BB并求其逆陣計算 計算電壓新值新值是否滿足要求滿足 停止置迭代次數k=0不滿足k+12.3 兩機五節點網絡潮流計算系統接線圖解： 1,.形成系數矩陣B、B并求他們的逆陣由于節點1為平衡節點，其他的節點均為PQ節點，系數矩陣，階數一樣。又應對該等值網絡，不存在除去與有功功率和電壓相位或無功功率和電壓大小關系較小的因素的可能性，這兩個矩陣B、B完全一樣。他們就由導納矩陣的虛部部分中除第一行第一列外的各個元素所組成，即Y=6.25-18.75j -5+15

25、j -1.25+3.75j 0 0 -5+15j 10.834-32.5j -1.667+5j -1.667+5j -2.5+7.5j -1.25+3.75j -1.667+5j 12.917-38.75j -10+30j0 0 -1.667+5j -10+30j 12.917-38.75j -1.25+3.75j0 -2.5+7.5j 0 -1.25+3.75j 3.75-11.25jB=B= -32.5005.0005.0007.5005.000-38.75030.00005.000 30.000-38.7503.750 7.5000 3.750-11.250 由此可見，網絡的節點導納矩陣

26、雖為奇異矩陣，但它的虛數部分的子陣B或B則是非奇異矩陣，可以求逆，其逆陣為= -0.056190-0.041905-0.44762-0.052381-0.041905-0.099048-0.087619 -0.057143-0.044762 -0.087619-0.105714-0.065079 -0.052381-0。057143 -0.065079 -0.145503 222 2.計算各節點有功功率的不平衡量 取=1.06，=0；=1.0；=0，按下式計算各節點有功功率不平衡量=-=0.201.01.06（-5.000cos0+15.000sin0）1.01.0（10.384cos032.

27、500sin0）1.01.0（-1.667cos0+5.000sin0）1.01.（-1.667cos0+5.000sin0）1.01.0(-2.500cos0+7.500sin0)=0.50000相似地可得=-0.375000；=-0.400000；=-0.600000 3.計算各節點電壓的相位角（弧度） 由方程式 ()=U可得=0.036952；=0.085524；=0.091810；=0.108571；所以=0.036952； =0.085524； =0.091810； =0.108571； 4.計算各節點無功功率不平衡量按下式計算各節點無功功率不平衡量=-=0.201.01.065.0

28、0sin（-0.0369520）15.000cos（-0.0369520）1.01.010.834sin0+32.5cos01.01.0-1.667sin（-0.036952+0.085524）5.000cos（-0.036952+0.085524）1.01.0-1.667sin（0.036952+0,091810）5.000cos（-0,036952+0.091810）1.01.0-2.500sin（0.036952+0.108571）7.500cos（-0.036952+0.108571）=1.211930相似地可得=0.077279； =0.191947； =0.319599； 5.計算

29、各節點電壓的大小 由方程式 ()=U可得=1.039528； =0.008050； =0.006386； =0.000072；所以=1.039528； =0.008050； =0.006386； =0.000072；求得各節點電壓的新值后，就開始第二次迭代。每次迭代所得適于表1表3。由表1可見經6次迭代就可滿足的要求。表1 迭代過程中各節點功率的不平衡量K0-0.5000001.211930-0.375000-0.0772790.400000-0.191947-0.600000-0.3195991-0.411720-0.1364970.0498020.0183000.0762480.02484

30、90.1208990.04021820.0432230.014875-0.005899-0.001936-0.007856-0.002664-0.012588-0.0042653-0.004833-0.001658-0.0006500.0002150.0008720.0002950.0013760.00046640.0005370.000184-0.000072-0,000024-0.000096-0.000033-0.000150-0.0000515-0.000060-0.0000210.0000080.0000030.0000110.0000040.0000160.00000660.000

31、0070.000000-0.0000010.000000-0.0000010.0000000-0.0000020.000000表2 迭代過程中各節點電壓的修正量k0-0.0369520.039528-0.0855240.008050-0.0918100.006386-0.1085710.0000721-0.010063-0.0034060.0018280.0007570.0024610.0009410.0045970.00161820.0010540.000385-0.000230-0.000064-0.000285-0.000082-0.000484-0.0001493-0.000120-0

32、.0000430.0000230.000006-0.0000280.0000080.0000490.00001540.0000130.000005-0.000002-0.000001-0.000003-0.000001-0.000005-0.0000025-0.000002-0.0000010.0000000.0000000.0000000.0000000.0000010.0000006表3 迭代過程中各節點電壓k00.0000001.0000000.0000001.0000000.0000001.0000000.0000001.0000001-0.0369521.039528-0.08552

33、41.008050-0.0918101.006386-0,1085711.0000722-0.0470161.036122-0.0836961.008808-0,0893481.007327-0,1039741.0016903-0.0459621.036507-0.0839261.008744-0.0896331.007245-0.1044581.0015404-0.0460821.036463-0.0839031.008750-0.0896051.007253-0.1044101.0015555-0.0460681.036468-0.0839061.008750-0.0896081.0072

34、52-0.1044151.0015546-0.0460701.036468-0,0839061.008750-0.0896081.007252-0.1044141.001554=1.036468； =1.008750=1.007252； =1.001554 6.計算平衡節點功率和線路功率 迭代收斂后，就可以計算平衡節點功率和線路功率。結果如下=1.298162+j0.244472各線路功率如表4表4 各線路功率j i1234510.889505+j0.1386620.408657+j0.1058102-0.875079-j0.0953850.246884+j0.0814570.279319+j

35、0.0806060.548870+j0.133213-0.395969-j0.067747-0.243108-j0.0701330.189079-j0.0121204-0.274598-j0.066445-0.188726+j0.0131780.063325+j0.0032675-0.536990-j0.097684-0.063008-j0.002316計算完畢。第三章 程序設計3.1 設計流程32 程序設計針對兩機五節點網絡潮流計算基于MATLAB的程序如下：Clc %清除當前屏幕上顯示的容，但并不清楚工作空間的數據Clear %清除工作空間的所有變量disp(節點總數為：); %顯示容N=

36、5disp(平衡節點為：);1Y=input(請輸入導納矩陣：Y=);pr=input(請輸入誤差精度：pr=);G=real(Y);B=imag(Y); %real實部，image虛部B1=B(2:5, 2:5); %去掉平衡節點B2=B1; b1=inv(B1); %取逆b2=inv(B2);ure=input(請輸入平衡節點1電壓初值的實部，Re=);dim=input(請輸入平衡節點1電壓初值的虛部，Im=);for m=2:5dt(m)=0; %電壓虛部初值u(m)=1.0; %電壓實部初值endfor m=2:5p(m)=input(請輸入PQ節點的有功功率：p=);q(m)=in

37、put(請輸入PQ節點的無功功率：q=);endk=0;wucha=1; %迭代次數，while wucha(1)pr;u(1)=ure;dt(1)=dim; %平衡節點電壓實部for m=2:5 %平衡節點電壓虛部for n=1:5pt(n)=u(m)*u(n)*(G(m,n)*cos(dt(m)-dt(n)+B(m,n)*sin(dt(m)-dt(n); %enddisp(各節點有功率不平衡量為：) dp(m)=p(m)-sum(pt) %=-endfor m=2:5u1=diag(u); %對角矩陣u2=inv(u1(2:5,2:5); %取逆u3=flipud(rot90(u(2:5)

38、; %逆時針90上下對稱u4=u3(1:5-1,:);dp1=flipud(rot90(dp(2:5); % ()=Udp2=dp1(1:5-1,:);dtt=(-b1*dp2/u4)*u2;dtt=dtt(:,1:1);dtt=flipud(rot90(dtt)disp(各節點電壓的相位角dt為：) dt(m)=dt(m)+dtt(m-1)endfor m=2:5for n=1:5 qt(n)=u(m)*u(n)*(G(m,n)*sin(dt(m)-dt(n)-B(m,n)*cos(dt(m)-dt(n);end %disp(各節點無功率不平衡量為：) (m)=q(m)-sum(qt) %

39、=-endfor m=2:5 1=flipud(rot90( (2:5); 2= 1(1:5-1,:);ut=-b2* 2/u4;disp( 各節點電壓的修正量：) utdisp( 各節點電壓的大小為：) u(m)=u(m)+ut(m-1) endfor n=1:5 U(n)=u(n)*(cos(dt(n)+j*sin(dt(n); %電壓endfor n=1:5 I(n)=Y(1,n)*U(n);enddisp(各節點功率s為：)S1=U(1)*sum(conj(I) %電流for m=1:5 for n=1:5 S(m,n)=U(m)*(conj(U(m)-conj(U(n)*conj(-

40、Y(m,n); %功率endenddisp(各節點功率sij為：)Swucha=max(abs(ut)k=k+1 %迭代+1end3.3 程序運行結果節點總數為：N = 5平衡節點為：ans = 1請輸入導納矩陣：Y=6.25-18.75i -5+15i -1.25+3.75i 0 0 -5+15i 10.834-32.5i -1.667+5i -1.667+5i -2.5+7.5i -1.25+3.75i -1.667+5i 12.917-38.75i -10+30i 0 0 -1.667+5i -10+30i 12.917-38.75i -1.25+3.75i0 -2.5+7.5i 0 -

41、1.25+3.75i 3.75-11.25i請輸入誤差精度：pr=0.00001請輸入平衡節點1電壓初值的實部，Re=1.06請輸入平衡節點1電壓初值的虛部，Im=0請輸入PQ節點的有功功率：p=0.2請輸入PQ節點的無功功率：q=0.2請輸入PQ節點的有功功率：p=-0.45請輸入PQ節點的無功功率：q=-0.15請輸入PQ節點的有功功率：p=-0.4請輸入PQ節點的無功功率：q=-0.05請輸入PQ節點的有功功率：p=-0.6請輸入PQ節點的無功功率：q=-0.1各節點有功率不平衡量為：dp = 0 0.5000各節點有功率不平衡量為：dp = 0 0.5000 -0.3750各節點有功率

42、不平衡量為：dp = 0 0.5000 -0.3750 -0.4000各節點有功率不平衡量為：dp = 0 0.5000 -0.3750 -0.4000 -0.6000dtt = -0.0370 -0.0855 -0.0918 -0.1086各節點電壓的相位角dt為：dt = 0 -0.0370 0 0 0dtt = -0.0370 -0.0855 -0.0918 -0.1086各節點電壓的相位角dt為：dt = 0 -0.0370 -0.0855 0 0dtt = -0.0370 -0.0855 -0.0918 -0.1086各節點電壓的相位角dt為：dt = 0 -0.0370 -0.08

43、55 -0.0918 0dtt = -0.0370 -0.0855 -0.0918 -0.1086各節點電壓的相位角dt為：dt = 0 -0.0370 -0.0855 -0.0918 -0.1086各節點無功率不平衡量為： = 0 1.2119各節點無功率不平衡量為： = 0 1.2119 -0.0773各節點無功率不平衡量為： = 0 1.2119 -0.0773 -0.1919各節點無功率不平衡量為： = 0 1.2119 -0.0773 -0.1919 -0.3196各節點電壓的修正量：ut = 0.0395 0 0 0 0.0081 0 0 0 0.0064 0 0 0 0.0001

44、 0 0 0各節點電壓的大小為：u = 1.0600 1.0395 1.0000 1.0000 1.0000各節點電壓的修正量：ut = 0.0395 0 0 0 0.0081 0 0 0 0.0064 0 0 0 0.0001 0 0 0各節點電壓的大小為：u = 1.0600 1.0395 1.0081 1.0000 1.0000各節點電壓的修正量：ut = 0.0395 0 0 0 0.0081 0 0 0 0.0064 0 0 0 0.0001 0 0 0各節點電壓的大小為：u = 1.0600 1.0395 1.0081 1.0064 1.0000各節點電壓的修正量：ut = 0.0

45、395 0 0 0 0.0081 0 0 0 0.0064 0 0 0 0.0001 0 0 0各節點電壓的大小為：u = 1.0600 1.0395 1.0081 1.0064 1.0001各節點功率s為：S1 = 1.1389 + 0.2403i各節點功率sij為：S =0 0.7229 + 0.1332i 0.4160 + 0.1071i 0 0 -0.7133 - 0.1044i 0 0.3110 + 0.0850i 0.3469 + 0.0845i 0.6671 + 0.1416i-0.4029 - 0.0676i -0.3052 - 0.0677i 0 0.2083 - 0.012

46、8i 0 0 -0.3398 - 0.0633i -0.2078 + 0.0141i 0 0.0714 + 0.0033i0 -0.6499 - 0.0900i 0 -0.0710 - 0.0021i 0 wucha = 0.0395 0 0 0k = 1各節點有功率不平衡量為：dp = 0 -0.4117 -0.3750 -0.4000 -0.6000各節點有功率不平衡量為：dp = 0 -0.4117 0.0498 -0.4000 -0.6000各節點有功率不平衡量為：dp = 0 -0.4117 0.0498 0.0762 -0.6000各節點有功率不平衡量為：dp = 0 -0.411

47、7 0.0498 0.0762 0.1209dtt = -0.0105 0.0012 0.0017 0.0035各節點電壓的相位角dt為：dt = 0 -0.0474 -0.0855 -0.0918 -0.1086dtt = -0.0105 0.0012 0.0017 0.0035各節點電壓的相位角dt為：dt = 0 -0.0474 -0.0843 -0.0918 -0.1086dtt = -0.0105 0.0012 0.0017 0.0035各節點電壓的相位角dt為：dt = 0 -0.0474 -0.0843 -0.0901 -0.1086dtt = -0.0105 0.0012 0.

48、0017 0.0035各節點電壓的相位角dt為：dt = 0 -0.0474 -0.0843 -0.0901 -0.1050各節點無功率不平衡量為： = 0 -0.1367 -0.0773 -0.1919 -0.3196各節點無功率不平衡量為： = 0 -0.1367 0.0175 -0.1919 -0.3196各節點無功率不平衡量為： = 0 -0.1367 0.0175 0.0237 -0.3196各節點無功率不平衡量為： = 0 -0.1367 0.0175 0.0237 0.0378各節點電壓的修正量：ut = -0.0038 0 0 0 0.0002 0 0 0 0.0004 0 0

49、 0 0.0009 0 0 0各節點電壓的大小為：u = 1.0600 1.0358 1.0081 1.0064 1.0001各節點電壓的修正量：ut = -0.0038 0 0 0 0.0002 0 0 0 0.0004 0 0 0 0.0009 0 0 0各節點電壓的大小為：u = 1.0600 1.0358 1.0083 1.0064 1.0001各節點電壓的修正量：ut = -0.0038 0 0 0 0.0002 0 0 0 0.0004 0 0 0 0.0009 0 0 0各節點電壓的大小為：u = 1.0600 1.0358 1.0083 1.0068 1.0001各節點電壓的修

50、正量：ut = -0.0038 0 0 0 0.0002 0 0 0 0.0004 0 0 0 0.0009 0 0 0各節點電壓的大小為：u = 1.0600 1.0358 1.0083 1.0068 1.0009各節點功率s為：S1 = 1.3260 + 0.2508i各節點功率sij為：S =0 0.9151 + 0.1435i 0.4109 + 0.1072i 0 0 -0.8998 - 0.0977i 0.2415 + 0.0816i 0.2741 + 0.0809i 0.5423 + 0.1344i -0.3981 - 0.0687i -0.2378 - 0.0707i 0 0.1

51、903 - 0.0114i 0 0 -0.2695 - 0.0672i -0.1899 + 0.0125i 0 0.0639 + 0.0036i0 -0.5306 - 0.0995i 0 -0.0636 - 0.0026i 0 wucha = 0.0038 0 0 0k = 2各節點有功率不平衡量為：dp = 0 0.0419 0.0498 0.0762 0.1209各節點有功率不平衡量為：dp = 0 0.0419 -0.0044 0.0762 0.1209各節點有功率不平衡量為：dp = 0 0.0419 -0.0044 -0.0044 0.1209各節點有功率不平衡量為：dp = 0 0

52、.0419 -0.0044 -0.0044 -0.0058dtt = 0.0016 0.0006 0.0006 0.0008各節點電壓的相位角dt為：dt = 0 -0.0459 -0.0843 -0.0901 -0.1050dtt = 0.0016 0.0006 0.0006 0.0008各節點電壓的相位角dt為：dt = 0 -0.0459 -0.0838 -0.0901 -0.1050dtt = 0.0016 0.0006 0.0006 0.0008各節點電壓的相位角dt為：dt = 0 -0.0459 -0.0838 -0.0895 -0.1050dtt = 0.0016 0.0006

53、 0.0006 0.0008各節點電壓的相位角dt為：dt = 0 -0.0459 -0.0838 -0.0895 -0.1043各節點無功率不平衡量為： = 0 0.0153 0.0175 0.0237 0.0378各節點無功率不平衡量為： = 0 0.0153 -0.0005 0.0237 0.0378各節點無功率不平衡量為： = 0 0.0153 -0.0005 -0.0006 0.0378各節點無功率不平衡量為： = 0 0.0153 -0.0005 -0.0006 -0.0001各節點電壓的修正量：ut = 1.0e-003 * 0.7776 0 0 0 0.5118 0 0 0 0

54、.5491 0 0 0 0.6921 0 0 0各節點電壓的大小為：u = 1.0600 1.0366 1.0083 1.0068 1.0009各節點電壓的修正量：ut = 1.0e-003 * 0.7776 0 0 0 0.5118 0 0 0 0.5491 0 0 0 0.6921 0 0 0各節點電壓的大小為：u = 1.0600 1.0366 1.0088 1.0068 1.0009各節點電壓的修正量：ut = 1.0e-003 * 0.7776 0 0 0 0.5118 0 0 0 0.5491 0 0 0 0.6921 0 0 0各節點電壓的大小為：u = 1.0600 1.036

55、6 1.0088 1.0073 1.0009各節點電壓的修正量：ut = 1.0e-003 * 0.7776 0 0 0 0.5118 0 0 0 0.5491 0 0 0 0.6921 0 0 0各節點電壓的大小為：u = 1.0600 1.0366 1.0088 1.0073 1.0016各節點功率s為：S1 = 1.2936 + 0.2441i各節點功率sij為：S =0 0.8855 + 0.1383i 0.4081 + 0.1057i 0 0 -0.8712 - 0.0955i 0 0.2475 + 0.0815i 0.2799 + 0.0806i 0.5495 + 0.1333i-

56、0.3954 - 0.0678i -0.2437 - 0.0701i 0 0.1889 - 0.0122i 0 0 -0.2751 - 0.0664i -0.1885 + 0.0132i 0 0.0632 + 0.0033i0 -0.5376 - 0.0976i 0 -0.0629 - 0.0023i 0 wucha = 1.0e-003 * 0.7776 0 0 0k = 3各節點有功率不平衡量為：dp = 0 -0.0057 -0.0044 -0.0044 -0.0058各節點有功率不平衡量為：dp = 0 -0.0057 0.0003 -0.0044 -0.0058各節點有功率不平衡量為

57、：dp = 0 -0.0057 0.0003 0.0004 -0.0058各節點有功率不平衡量為：dp = 0 -0.0057 0.0003 0.0004 0.0005dtt = 1.0e-003 * -0.2442 -0.1346 -0.1413 -0.1674各節點電壓的相位角dt為：dt = 0 -0.0461 -0.0838 -0.0895 -0.1043dtt = 1.0e-003 * -0.2442 -0.1346 -0.1413 -0.1674各節點電壓的相位角dt為：dt = 0 -0.0461 -0.0839 -0.0895 -0.1043dtt = 1.0e-003 * -

58、0.2442 -0.1346 -0.1413 -0.1674各節點電壓的相位角dt為：dt = 0 -0.0461 -0.0839 -0.0896 -0.1043dtt = 1.0e-003 * -0.2442 -0.1346 -0.1413 -0.1674各節點電壓的相位角dt為：dt = 0 -0.0461 -0.0839 -0.0896 -0.1044各節點無功率不平衡量為： = 0 -0.0020 -0.0005 -0.0006 -0.0001各節點無功率不平衡量為： = 0 -0.0020 0.0001 -0.0006 -0.0001各節點無功率不平衡量為： = 0 -0.0020

59、0.0001 0.0001 -0.0001各節點無功率不平衡量為： = 0 -0.0020 0.0001 0.0001 0.0001各節點電壓的修正量：ut = 1.0e-004 * -0.9362 0 0 0 -0.5617 0 0 0 -0.6002 0 0 0 -0.7119 0 0 0各節點電壓的大小為：u = 1.0600 1.0365 1.0088 1.0073 1.0016各節點電壓的修正量：ut = 1.0e-004 * -0.9362 0 0 0 -0.5617 0 0 0 -0.6002 0 0 0 -0.7119 0 0 0各節點電壓的大小為：u = 1.0600 1.0

60、365 1.0087 1.0073 1.0016各節點電壓的修正量：ut = 1.0e-004 * -0.9362 0 0 0 -0.5617 0 0 0 -0.6002 0 0 0 -0.7119 0 0 0各節點電壓的大小為：u = 1.0600 1.0365 1.0087 1.0072 1.0016各節點電壓的修正量：ut = 1.0e-004 * -0.9362 0 0 0 -0.5617 0 0 0 -0.6002 0 0 0 -0.7119 0 0 0各節點電壓的大小為：u = 1.0600 1.0365 1.0087 1.0072 1.0015各節點功率s為：S1 = 1.298