NorthChinaElectricPowerUniversity電力工程系DepartmentofElectricalEngineering電力系統穩定技術研究組

常鮮戎

Email:changxr@電力系統穩定分析和計算內容提要1.電力系統穩定研究的意義2.電力系統穩定分析3.電力系統穩定計算4.電力系統穩定研究的方向電力系統穩定研究的意義1.電力系統穩定破壞的損失

保持電力系統的穩定性是電力系統安全穩定運行的重要任務。由于電力系統在國民經濟中的特殊地位,電網穩定事故所造成的經濟損失十分巨大，帶來的社會后果也十分嚴重。

以1977年7月13日發生在美國紐約的大停電事故為例，美國能源部估計的經濟損失約為35億5千萬美元。另外，在停電持續的25小時內，紐約市發生火災1037起，暴力事件1809起，全部交通中斷。影響了900萬人的工作與生活。前年發生的北美大停電的損失更大，初步估計經濟損失達到300億美元。

電力系統穩定研究的意義2.電力系統穩定研究的目的探討各種穩定問題發生的機理:電壓穩定問題產生的原因低頻振蕩產生的原因次同步諧振的原因

電力系統穩定研究的意義

尋找避免或減輕穩定破壞的方法:尋找穩定域,使系統運行在穩定域內確定提高穩定性的措施,擴大系統的穩定域在系統穩定遭到破壞時,找到使系統回到穩定域的措施電力系統穩定研究的意義電力系統穩定分析1.電力系統小干擾分析法小干擾分析法可用來分析電力系統在小擾動條件下的穩定性,如靜態穩定性，低頻振蕩等。電力系統穩定分析1-1.系統狀態方程為了研究電力系統小干擾穩定性，首先要建立電力系統的狀態方程。電力系統穩定分析如果一系統的所有狀態變量x的變化率都不是時間t的顯函數，則稱該系統為自治系統。此時方程(2-1)可簡化為：

電力系統穩定分析集合｛x1,x2,…,xn｝是系統(2-2)的一個狀態。系統的狀態是描述該系統行為的一組最少信息。當已知系統在任意時刻t0的狀態x0后，就可根據系統t≥t0時的輸入描述該系統t≥t0后的行為，而不需要知道系統t＜t0時的輸入。任意一組n個線性獨立的系統變量都可以用來表示系統的狀態，這些變量稱為狀態變量。系統的任何其它變量都可以通過狀態變量來表示。電力系統穩定定分析系統的狀態變量可以是是該系統的物物理變量，也可以是是描述該系統的純純粹數學變量。盡管管在任意時時刻系統的狀態態是唯一的的，但系統狀態變量的選擇擇不是唯一一的，即描描述系統狀態的的信息不是唯唯一的。描述系統狀態的的n維歐氏空間稱為該該系統的狀態空間。當系統的狀態隨時間間變化時，，在狀態空空間代表系系統狀態的的點將構成成一軌跡，，稱為狀態軌跡。電力系統穩穩定分析1-2系統狀態方方程的線性化設x0，u0分別是非線性系統(2-2)在所關注平衡點的狀態向量和輸入入向量。若此時系統統受到一小小干擾，使得得：電力系統穩穩定分析這個新狀態態也滿足式式(2-2)，因此：將非線性函數在平衡點作Taylor展開。忽略略二次及以以上高次項項后有：電力系統穩穩定分析因此，非線性系統(2-2)的線性化狀狀態方程為：把分別重新記記為有：電力系統穩穩定分析1-3系統狀態方方程的本征征特性特征根與特特征向量設A是線性動態態系統狀態態方程的系系數矩陣,若滿足這個方方程的標量量為矩陣A的特征根。一個動態系系統不同的的狀態方程程有相同的特征根。。電力系統穩穩定分析滿足方程的非零向量量為矩陣A的對于特征根的右特征向量量。滿足方程的非零向量量為矩陣A的對于特征根的左特征向量量。電力系統穩穩定分析對應于不同同特征根的的右特征向向量和左特特征向量是是正交的，即而對應于同同一特征根根的右特征征向量和左左特征向量量有關系電力系統穩穩定分析特征根與系系統穩定性性⑴實數特征根根對應于非非振蕩模式式。負實數數特征根對對應于衰減減模式，特特征根的幅幅值越大，，衰減越快快。正實數數特征根對對應于非周周期失穩。。⑵復數特征根根以共軛形形式出現，，每一對對對應于一個個振蕩模式式。電力系統穩穩定分析模式分布形形態右特征向量量給出了系系統動態模模式的分布布形態右特征向量量的元素的幅值表示示第式在第個狀態變量量中的幅度。個動態模的模大，反映了對的可觀性強。電力系統穩穩定分析左特征向量量的元素的幅值表示示第量個動態模式式中的比重重。個狀態變的模大，反映了的變化可使使有較大變化化,可控性強。。在第電力系統穩穩定分析1-4電力系統低低頻振蕩分分析與電力力系統穩定定器電力系統低低頻振蕩分分析先分析簡單單電力系統統的低頻振振蕩問題。。由于勵磁磁調節系統統在電力系系統低頻振振蕩分析方方面起著很很重要的作作用，因此此在分析電電力系統低低頻振蕩時時，發電機機組的模型型要包括勵勵磁系統的的模型。電力系統穩穩定分析對單機——無窮大系統統，分析低低頻振蕩問問題的模型型為：電力系統穩穩定分析計算表明：：在系統運運行方式變變化時，及都是正數，而而在重負荷即即較大時變為為負數。這一現象在在低頻振蕩蕩分析時是是很重要的的。下面面分分析析發發電電機機轉轉子子繞繞組組及及勵勵磁磁對對低低頻頻振振蕩蕩的的影影響響。。電力力系系統統穩穩定定分分析析⑴設勵勵磁磁系系統統輸輸出出為常常數數。。此時時，，狀態態方程程（（2-5）為為三電力力系系統統穩穩定定分分析析從（（2-6）式可可以推推得：：由（（2-7）式可可得得以以下下結結論論：：①主主要要影影響響振振蕩蕩頻頻率率;②主主要要影影響響振振蕩蕩阻尼尼。。電力力系系統統穩穩定定分分析析此時時所以以發發電電機機勵勵磁磁繞繞組組的的動動態態作作用用有有助于于抑抑制制低低頻頻振振蕩蕩。。⑵勵磁磁系統統對對低頻頻振振蕩蕩的的影響響由（2-5）第第四四式式有有：：電力系系統穩穩定分分析將其代代入（（2-5）第三三式，，整理得:式中為為勵磁磁系統的的放大大倍數數，高高放大大倍數數時，，。與相乘，，將加加速系系統出出現負負阻尼的進程程。電力系系統穩穩定分分析PSS的工作作原理理電力系系統出出現低低頻振振蕩是是由于于勵磁磁調節節系統統產生生了負負阻尼尼，如如果能能在勵勵磁調調節系系統引引入附附加控控制功功能，，使其其產生生正阻阻尼，，抵消消由于于變變負產產生的的負阻阻尼，，就能能抑制制電力力系統統的低低頻振振蕩。。這就就是電電力系系統穩穩定器器（PowerSystemStabilizer簡稱PSS）的設計計思想想。電力系系統穩穩定分分析下面我我們分分析取取為為輸入入信號號的PSS裝置。。將PSS信號引引入勵勵磁調調節通通道，，則發發電機機勵磁磁電勢勢為：：電力系系統穩穩定分分析由(2-8)可推得得:電力系系統穩穩定分分析若要產產生正阻尼尼，則則（2-9）式對對應于于的系數數應為式中：：為為正正實數數。電力系系統穩穩定分分析所以：：應應該為為:由于，隨隨系系統運運行狀狀況變變化，，而且且在現現實情情況下下很難難構造造純超超期環環節，，所以以實際際上取取電力系系統穩穩定分分析將（2-10）式代代入（（2-9）式，分析析得:當，，時時，對應于于的實數數部分分產生生正阻尼尼。電力系系統穩穩定分分析多機系系統低低頻振振蕩分分析若發電電機采采用三三階實實用模模型，，勵磁磁系統統用一一階模模型，，忽略略調速速器動動態，，負荷荷只計計及靜靜態效效應，，則多多機系系統的的線性性化狀狀態方方程如如（2-5）式所所示。只不不過式中相應的的變量量都為為向量量，相相應的的系數數都為為矩陣陣。電力系系統穩穩定分分析多機系系統低低頻振振蕩分分析的的主要要內容容有：：①計算系系統的的特征征根及及左，，右特特征向向量，。一般般用QR法計算算。②從這些些特征征根中中挑出出振蕩蕩頻率率為的特征征根，計算算其與與各狀狀態變變量的的相相關因子子和和機電電回路路相關關比，鑒別別出感感興趣的的機電電振蕩蕩模式式。電力系系統穩穩定分分析③分析機機電振振蕩模模式的的振振蕩頻頻率和和阻尼尼特性性，并并根據據其特特征向向量分分析該該振蕩蕩模式式在各各機觀觀察察時的的相對對振幅幅和相相位，，從而而求出出該模模式發發生在在哪兩兩臺機機組（（或機機群））之間間。④根據相相關因因子判判斷斷機電電模式式同同哪哪些機機組強強相關關，確確定安安裝PSS的地點點。電力系系統穩穩定分分析⑤通過靈靈敏度度分分析，，得到到和和的的相互互關系系，取取為為PSS放大倍倍數時時，可可提供供PSS參數設設置所所需信信息。。電力系系統穩穩定分分析2.電力系系統時時域仿仿真分分析法法時域仿仿真分分析方方法是是利用用計算算機求求解描描述電電力系系統動動態過過程的的微分分方程程組,得到電電力系系統各各變量量在動動態過過程的的數值值解。。電力系系統穩穩定分分析時域仿仿真分分析原原理及及基本本要求求時域仿仿真分分析方方法的的原理理是基基于微微分方方程的的數值值解法法。用用時域域仿真真分析析的方方法分分析電電力系系統的的動態態特性性時,首先要確定定描述電力力系統動態態特性的微微分方程,即電力系統統的數學模型,然后再選擇擇合適的數值計算方方法進行求解,最后,還要對計算算結果進行行輸出格式式的處理,通過良好的的人-機接口界面面輸出結果果,以方便閱讀讀和分析。。電力系統穩穩定分析對時域仿真真分析方法法的要求:①要有足夠的的計算精度度;②要保證數值值計算過程程的數值穩穩定性;③要有盡量快快的計算速速度。電力系統穩穩定分析時域仿真方方法在電力力系統穩定定分析的應應用時域仿真分分析方法在在電力系統統各種穩定定分析計算算中都在應應用,如電力系統統暫態穩定定分析計算算、電力系系統中、長長期動態過過程分析計計算、電力力系統電壓壓穩定分析析計算、電電力系統次次同步諧振振分析計算算等等。電力系統穩穩定分析時域仿真分分析方法需需在以下幾幾方面繼續續完善:①研究合適的的數學模型型;②研究有較高高數值穩定定性及計算算精度的算算法;③研究適合電電力系統時時域仿真分分析的計算算平臺。電力系統穩穩定分析電力系統穩穩定分析中中時域仿真真算法的新新進展電力系統穩穩定的時域域仿真分析析使用的算算法一般為為顯式四階階龍格-庫塔法和隱隱式梯形法法。顯式四四階龍格-庫塔法的計計算精度高高,但數值穩定定性較差。。仿真時間間較長時容容易出現數數值不穩定定現象。隱隱式梯形法法的數值穩穩定性好,但計算精度度較低,只有二階精精度。電力系統穩穩定分析隱式式泰泰勒勒級級數數展展開開法法的的計計算算精精度度高高,而且且精精度度可可人人為為設設置置,使用用時時根根據據需需要要設設置置合合適適的的精精度度。。數數值值穩穩定定性性好好,經適適當當設設置置,隱式式泰泰勒勒級級數數展展開開法法具具有有數數值值計計算算算算法法的的A類穩穩定定性性。。另另外外隱隱式式泰泰勒勒級級數數展展開開法法有有較較快快的的計計算算速速度度,在計計算算精精度度相相同同的的條條件件下下,隱式式泰泰勒勒級級數數展展開開法法的的計計算算速速度度是是顯顯式式四四階階龍龍格格-庫塔塔法法的的三三倍倍。。電力力系系統統穩穩定定分分析析使用用隱隱式式泰泰勒勒級級數數展展開開法法的的電電力力系系統統暫暫態態穩穩定定分分析析軟軟件件可可在在一一般般微微型型計電力力系系統統穩穩定定分分析析3.電力力系系統統穩穩定定性性的的直直接接法法分分析析常規規的的計計算算系系統統穩穩定定性性的的數數值值分分析析計計算算的的方方法法目目前前還還不不能能做做到到實實時時計計算算,只能能采采取取事事先先離離線線計計算算的的方方法法,對系統各種運運行方式下可可能發生的故故障都進行計計算,然后按最嚴重重的情況整定定各種緊急措措施的設備。。由于這些設設備是按離線線整定的方式式運行,一旦電網遇到到的擾動與離離線計算的擾擾動情況不同同,這些緊急措施施的效果就可可能不太理想想。電力系統穩定定分析直接法分析原原理1892年Lyapunov在<<關于運動穩定定性的一般問問題>>的著作中提出出了一種不用用求解微分方方程而直接判判斷運動穩定定性的方法。。電力系統穩定定分析對于任一運動動系統:如果存在標量量函數,其對時間的導導數,則系統在原點點是漸近穩定定的。電力系統穩定定分析直接法可以不不用求解系統統微分方程而而僅根據標量量函數及其導導函數的數值值判斷系統的的穩定性,因而大大提高高了判斷系統統的穩定性的的速度,受到人們廣泛泛的重視。電力系統穩定定分析直接法在電力力系統的應用用將直接法用于電電力系統暫態態穩定分析已已有幾十年歷史。同其它領域域的應用相似似,直接法在電力力系統的應用用研究也一直直致力于構造造一個合理的的Lyapunov函數。電力系統穩定定分析近年來,Lyapunov函數的定義一一般都基于發發電機轉子運運動方程式對對轉子角的一一階積分基礎礎上。在單機機-無窮大系統,構造的能量函函數:電力系統穩定定分析研究的重點則則側重于如何何快速、正確確、可靠地確確定臨界能量這一問題上。。用直接法判斷斷電力系統暫暫態穩定性時時,只需求出擾動動結束時的,,據此計算系統統總能量,并確定臨界能能量。通過比較和和來判斷穩定性性。電力系統穩定定分析為了描述系統統的穩定程度度,可用來度量系統穩穩定度。一般般使用規格化化的穩定度來度量系統穩穩定度。電力系統穩定定分析電力系統使用用的幾種直接接法1.相關不穩定平平衡點法RUEPUEP法的基本假定定是:當系統出現不不穩定時,系統運行點肯肯定要穿過不不穩定平衡點點從穩定區域域進入不穩定定區域。系統統在不穩定平平衡點的勢能能就就是系統統的臨界能量量。只要擾動結束束后系統的總總能量,系統就是穩定定的。電力系統穩定定分析對于多機系統統,有個個失穩穩模式,相應的有個個不穩定平平衡點。應該該用哪個不穩穩定平衡點的的作作為系統臨臨界能量呢呢?一種辦法是計計算所有個個不穩定平平衡點的,然后取最小的的為為系統的臨界界能量。電力系統穩定定分析系統只能以這這個失穩模式中的某一個模模式失去穩定定,與該失穩模式式對應的系統統軌跡只能在在某一個不穩穩定平衡點穿穿過。取該不不穩定平衡點點的作作為系統臨臨界能量才才能反映系系統的實際狀狀況。由于系系統的失穩模模式與系統的運行方方式、故障地地點及故障類類型等運行條條件相隔,故該不穩定平平衡點稱為相相關不穩定平平衡點RUEP。電力系統穩定定分析使用RUEP法判斷系統穩穩定性時要先先判別主導失失穩模式,解出相關不穩穩定平衡點,然后再計算。。該方法法稱為RUEP法。使用RUEP法后,的計算精度和和速度比UEP法有了很大提提高。電力系統穩定定分析設有一n機系統,以各機組的角角度為坐標繪繪出系統勢能能與各機組角角度的關系曲曲面,在勢能曲面上上把勢能相等等的點連起來來,構成一個個閉閉合的等勢能線。在圖上找到到各個不穩定定平衡點的位位置,然后過這些不不穩定平衡點點畫一條曲線線,要求該曲線與與各處的等勢勢能線正交。。這條連接各各不穩定平衡衡點并與各處處的等勢能線線正交的曲線線稱為勢能界面。2.勢能界面法PEBS電力系統穩定定分析PEBS法的基本假定定是:在UEP附近系統勢能能的變化比較較平衡,因此可用持續續故障的軌跡跡與PEBS的交點所所對應的系統統勢能近近似代替UEP點的系系統勢能。電力系統穩定定分析使用PEBS計算極限切除除時間的步驟驟為:①進行故障前的的系統潮流計計算,求出穩定計算算所需的初值值(SEP);②計算持續故障障下的系統軌軌跡;③在每一時段計計算勢能;④重復計算上述述步驟,一直到達達到最大,此時的即即為臨界能量量;電力系統穩定定分析3.擴展等面積法法EEAC等面積法則在在單機-無窮大系統是是判斷系統暫暫態穩定的充充要條件,如能把等面積積法則應用到到多機系統,則可以快速、、準確地判斷斷出系統的穩穩定性。這就就是擴展等面面積法EEAC。電力系統穩定定分析EEAC將觀察系統穩穩定性的觀察察空間和計算算系統軌跡的的積分空間分分開處理,在n維實數域內對對系統的的模型進行積積分運算,而將數值解中中的各發電機機轉角軌跡用用參量形式映映射到一系列列低維觀察空空間中,再將等面積法法則應用到低低維觀察空間間的轉角影像像軌跡來判斷斷中系統統的穩定性。。電力系統穩定定分析多機系統的互互補分割:一個n機系統有且僅有種種不同的的方式將n臺機分割為兩兩個非空的互互補機群和和電力系統穩定定分析對于某個特定定的劃分,把屬于同一群群的機組的轉轉子運動方程程相加,可得到兩個多多元微分方程程:(2-11)電力系統穩定定分析定義1PCOI(n,2)映射是從n維實數空間到到各個個等值兩機空空間的一種線線性變換,記作PCOI(n,2):對某個互補劃劃分來說說,其變換函數為為:其中:電力系統穩定定分析利用此變換,把(2-11)式改寫為兩個個參量化的二二元微分方程程(2-12)電力系統穩定定分析定義2:OMIB變換將各兩機機觀察子空間間映射到相應應的單機平面面,記作OMIB:。其變換函數為為:(2-13)OMIB變換將(2-12)式進一步變換換為電力系統穩定定分析定義3PCOI(n,1)映射是把中中的數值值積分解依次次經過PCOI(n,2)和OMIB兩次變換,而得到映映像。記作作PCOI(n,1):或簡稱PCOI映射。利用PCOI映射把觀察空空間和積分空空間分開,即把多機系統統的暫態穩定定分析問題變變換到一系列列映像兩機系系統的子問題題。而且,判斷的結果在在軌跡跡精度的含義義上是精確的的。電力系系統穩穩定分分析PCOI映射保保留了了的的穩穩定狀狀態和和臨界界條件件,但個個映像像子系系統卻卻是一一個可可觀的的數字字。必必須有有效地地篩選選出極極少數數的候候選,只對它它們進進行觀觀察。。這是是關系系到PCOI映射能能否實實用的的問題題。電力系系統穩穩定分分析按臨界界不穩穩定軌軌跡通通過NRP處的各各機轉轉角來來確定定發電電機臨臨界度度隊列列CMR的方法法稱CMR法。CMR法把需需要觀觀察的的平平面數數從個個減減少到到個個。。電力系系統穩穩定分分析進一步步將個個候選選映像像按對對應的的角度度間隙隙重新新再排排一次次隊,并把注注意力力放在在隊首首的少少量候候選映映像。。“角角間隙隙概念念”進進一步步加快快了CCI過程。。電力系系統穩穩定計計算1.電力系系統數數學模模型⑴同步發發電機機模型型根據同同步發發電機機模擬擬的詳詳盡程程度,同步發發電機機有以以下幾幾種模模型。。模型Ⅰ：按派派克模模型導導出的的精細細模型型，六六繞組組（定定子兩兩繞組組，轉轉子四四繞組組）或或五繞繞組，，考慮慮轉子子勵磁磁繞組組及阻阻尼繞繞組的的暫態態過程程。六六或五五階。。電力系系統穩穩定計計算定子繞繞組方方程式為：電力系系統穩穩定計計算采用定定子繞繞組電電動勢勢來描描述的的轉子子繞組組電磁磁暫態態過程程方程程式為：電力系系統穩穩定計計算轉子運運動方方程式式為：：電力系系統穩穩定計計算電磁力力矩方方程式式為：：電力系系統穩穩定計計算模型Ⅱ:三階，，除轉轉子運運動方方程外外，還還考慮慮轉子子勵磁磁繞組組的暫暫態過過程，，但不不考慮慮轉子子阻尼尼繞組組的暫暫態過過程。。該模型型考慮慮勵磁磁電壓壓的變變化，，但不不考慮慮轉子子阻尼尼繞組組的作作用，，即認認為和和繞組組均開路。。（3-2）式各各參數數可取取為:電力系系統穩穩定計計算定子繞繞組方程式式可寫寫為：：電力力系系統統穩穩定定計計算算轉子子繞繞組組方方程程式簡化化為為：：電力力系系統統穩穩定定計計算算模型型Ⅲ：二二階階，，不不考考慮慮轉轉子子勵勵磁磁繞繞組組及及阻阻尼尼繞繞組組的的暫暫態態過過程程，，即即認認為為恒定定。。該模模型型不不考考慮慮勵勵磁磁電電壓壓變變化化和和轉轉子子阻阻尼尼繞繞組組的的作作用用。。（（3-2）式式各各參參數數取取值值同同模模型型Ⅱ。電力力系系統統穩穩定定計計算算模型型Ⅳ：二二階階，，不不考考慮慮轉轉子子勵勵磁磁繞繞組組及及阻阻尼尼繞繞組組的的暫暫態態過過程程，，后的的恒恒定定。。該模模型型（（3-2）式式各各參參數數可可取取為為：：均取取的值值。。電力力系系統統穩穩定定計計算算該模模型型的的定子子回路路方方程程為為：：電力力系系統統穩穩定定計計算算⑵勵磁磁系系統統及及調調節節器器模模型型⑶原動動機機及及調調速速器器模模型型⑷負荷荷模型型⑸電力力網網絡絡模模型型電力力系系統統穩穩定定計計算算xy-dq兩種種不不同同坐坐標標系系的的變變換換關關系系式式如如下下:電力力系系統統穩穩定定計計算算⑹電力力系系統統模模型型的的一一般般形形式式整個個電電力力系系統統模模型型是是由由一一組組一一階階微微分分方方程程式式，，即即和一一組組代代數數方方程程式式，，即即電力力系系統統穩穩定定計計算算（3-10）式式也也可可表表示示為為一一種種準準線線性性方方程程，，即即電力力系系統統穩穩定定計計算算2.電力力系系統統小干干擾擾穩穩定定性性計計算算電力力系系統統小小干干擾擾穩穩定定性性分分析析是是研研究究電電力力系系統統在在某某個個運運行行工工況況下下受受到到小小干干擾擾后后電電力力系系統統能能否否保保持持同同步步運運行行的的問問題題。。根根據據研研究究的的重重點點和和深深度度不不同同，，所所涉涉及及的的電電力力系系統統各各部部件件的的方方程程也也不不同同。。一一般般有有以以下下方方程程：：電力力系系統統穩穩定定計計算算⑴同步步機機組組轉轉子子運運動動方方程程研究究電電力力系系統統小小干干擾擾穩穩定定性性的的系系統統狀態態方程程必必須須有有能能反反映映同同步步機機組組轉轉速速和和角角度度的的各各同同步步機機組組的的轉轉子子運運動動方方程程：：電力力系系統統穩穩定定計計算算將上式在在運運行行點點線線性性化化，，整整理理得得：：電力力系系統統穩穩定定計計算算⑵原動動機機功功率率方程程分析電力力系統小小干擾穩穩定性時時，通常常有以下下簡化條條件：⑴原動機功功率（轉轉矩）恒恒定，即即⑵用恒定阻阻抗代替替負荷；；⑶不計電力力網絡內內的電磁磁暫態過過程。電力系統統穩定計計算發電機采采用，，模模型時時多機系系統小干干擾穩定定性的計計算步驟驟：⑴確定待分分析的電電力系統統某一運運行方式式并作潮潮流計算算，算出出系統各各節點的的電壓相相量和和各發電機輸輸出功率⑵根據給定定的節點點負荷功功率和和對對應的節節點電壓壓，求出代代替負荷荷功率的的導納。即用恒定導納（阻抗）代代替負荷荷。電力系統統穩定計計算⑶修正網絡絡方程。。且把發電機節點排在在前面。。⑷消去聯絡節點點。整理得由發電