【正文】 ),title(39。 for j=1:n Dy(K,j)=Dy(K1,j)。)。 fori=1:nlif B1(i,6)==0 p=B1(i,1)。q=B1(i,2)。enddisp(39。disp(V)。O(k)=O(k)*180./pi。每次沒有達到精度要求的有功功率個數(shù)為39。elsekq=0。DQ(I)=DQ(I)A(I,L)*DQ(L)。DQ(L)=DQ1(i)./V(i)。if DET =pr ICT3=ICT3+1。endkq=1。if ICT2~=0fori=2:nDP(i)=B(i,i)*DP(i)。ICT3=0。endendendendICT2=1。k=0。 endendfori=2:n ifi==n B(i,i)=1./B(i,i)。BI(i,i)=BI(i,i)+B2(i,5)。 Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2。endendfori=1:nlif B1(i,6)==0 p=B1(i,1)。YI=zeros(n)。)。 pr=input(39。請輸入節(jié)點數(shù)：n=39。 0 1 0 3]。3 6 +i* 1i* 1 0。PowerWorld的定性分析如下表表1勵磁參數(shù)對穩(wěn)定計算結果的影響勵磁系統(tǒng)參數(shù)TRKATATEKESEKFTF參數(shù)初值02010參數(shù)變化趨勢增大增大增大增大增大增大增大增大減小不變減小減小增大減小不變不變從表17知，可以通過改變勵磁系統(tǒng)的參數(shù)來提高系統(tǒng)的暫態(tài)穩(wěn)定性，即減小TR、減小TA、減小TE，增大KE，減小SE有利于改善系統(tǒng)的暫態(tài)穩(wěn)定。圖4采用二階模型時暫定過渡時間較大（故障方案一、tc=）圖4采用三階模型時暫定過渡時間較?。ü收戏桨敢弧c=） TJ對穩(wěn)定計算結果的影響慣性時間常數(shù)TJ是反映發(fā)電機轉(zhuǎn)子機械慣性的重要參數(shù)，由TJ的定義可知，它是轉(zhuǎn)子在額定轉(zhuǎn)速下的動能的兩倍除以基準功率。圖47綜上，取tc。所得搖擺曲線如下，此時系統(tǒng)暫態(tài)穩(wěn)定。具體參數(shù)設置分別如圖339所示圖3發(fā)電機GENTRA模型參數(shù)設置圖3勵磁系統(tǒng)BPA_EA模型參數(shù)設置（1）在5s節(jié)點3時發(fā)生三相對稱短路，所得搖擺曲線如下，即此時系統(tǒng)暫態(tài)穩(wěn)定。臨界角為，滿足式，即（20）求得設極限切除角為，由等面積定則有（21）將參數(shù)代入后解方程得，極限切除角幾種方法求出來的結果統(tǒng)計如下：表1幾種方法求出來的結果統(tǒng)計計算方法MatlabPowerWorld等面積法（s）極限切除角（度）從表13可知，采用等面積法求出的極限切除角是最小的，其次是Matlab仿真法，最大的是PowerWorld仿真法。圖33（2）在事件窗口添加：1s時線路L23中點處發(fā)生三相對稱短路，得出的搖擺曲線如下，此時系統(tǒng)暫態(tài)不穩(wěn)定。,t,x(:,2),39。simb39。%短路后的P2[t1,x1]=ode45(39。%定義轉(zhuǎn)子運動方程組主程序:w0=1。故障后功角隨時間的變化曲線如圖29所示，從圖中可以讀出極限切除角=。圖25（2）在事件窗口添加：1s時節(jié)點3發(fā)生三相對稱短路，得出的搖擺曲線如下，此時系統(tǒng)暫態(tài)不穩(wěn)定。,[0,],x0,h_opt,a,b,c)。h_opt=odeset。b39。,[tc,20],x02,h_opt,a,b,c)。simb39。tc=。當節(jié)點3發(fā)生三相對稱短路時，通過網(wǎng)絡變換法同理可得，轉(zhuǎn)移阻抗Z17= + =，輸入阻抗Z77= + =。網(wǎng)絡變換時，先將這些并聯(lián)電容消去。所以，兩種計算潮流得出的支路潮流具有很好的一致性。最后，用輸電線路將各節(jié)點連接起來。此外，系統(tǒng)模型可以通過使用仿真軟件的圖形編輯工具很容易地進行修改，用戶只需輕輕點擊幾下鼠標就可以在檢修期間切換線路、增加新的線路或發(fā)電機、確定新的交易容量。 節(jié)點分類標號：1為平衡節(jié)點（應為1號節(jié)點）；2為PQ節(jié)點；3為PV節(jié)點；根據(jù)前文對待求解電網(wǎng)的線路參數(shù)的歸算和標幺化，所構造的支路信息矩陣B1和構造節(jié)點信息矩陣B2如下所示（MATLAB界面截圖表示），圖MATLAB界面截圖2）接下來，根據(jù)支路信息矩陣B1和節(jié)點信息矩陣B2，程序先形成待求解電網(wǎng)的節(jié)點導納矩陣：先對空矩陣加入對角線元素，即先加入每個節(jié)點的對地導納y0；根據(jù)支路信息矩陣B1，對導納矩陣再加入互導納ypq；再將互導納ypq加入到導納矩陣的對角線元素，形成自導納。1 潮流的參數(shù)計算和標幺化 線路參數(shù)計算與標幺化圖電網(wǎng)的等效線路圖圖1是待求解電網(wǎng)的等效線路圖，由已知取基準功率100MW，基準電壓220kV：則有 線路正序參數(shù)： 線路阻抗標幺值計算公式： 線路導納標幺值計算公式： 由于電網(wǎng)線路為雙回路線路，即阻抗值為單回路的一半，導納值為單回路的兩倍。對上面23種故障方案，采用等面積定則計算故障的極限切除角，并分析比較與時域仿真法計算結果的差別；用Matlab編制穩(wěn)定計算程序，發(fā)電機采用三階模型，并對第6步的23種故障方案進行穩(wěn)定計算，給出搖擺曲線，并計算故障的極限切除時間；與PowerWorld軟件的分析結果進行比較校核。發(fā)電廠按PV方式運行。F/km；零序參數(shù)： r0 = , x0 = , C0 = 181。F/km；40186?？紤]該電廠開機三臺，所有發(fā)電機保留10%的功率裕度。比較兩種模型的仿真結果，分析發(fā)電機模型選擇對于穩(wěn)定計算結果的影響，并且分析勵磁調(diào)節(jié)系統(tǒng)參數(shù)變化對于穩(wěn)定計算結果的影響；編制課程設計報告。L12的線路長度為27km，則Z12=*（+）*27/484Ω=(+) ΩB12=2*2**50**106*27*484S=L23的線路長度為17km，則Z23=*（+）*17/484Ω=(+) ΩB23=2*2**50**106*17*484S=L24的線路長度為6km，則Z24=*（+）*6/484Ω=(+) ΩB24=2*2**50**106*6*484S=L45的線路長度為5km，則Z45=*（+）*5/484Ω=(+) ΩB45=2*2**50**106*5*484S=L36的線路長度為18km，則Z36=*（+）*18/484Ω=(+) ΩB36=2*2**50**106*18*484S=L56的線路長度為24km，則Z56=*（+）*24/484Ω=(+) ΩB56=2*2**50**106*24*484S=綜上，電網(wǎng)線路參數(shù)的標幺值如表1所示：表線路參數(shù)的標幺值統(tǒng)計線路L12L23L24L45L36L56阻抗Z/Ω++++++導納B/S 節(jié)點信息統(tǒng)計表節(jié)點給定參數(shù)統(tǒng)計編號123456節(jié)點類型平衡節(jié)點PQ節(jié)點PQ節(jié)點PQ節(jié)點PQ節(jié)點PV節(jié)點P(MW)無250380280420Q(MVar)無40908595無V(kV)230無無無無由圖1可知，網(wǎng)絡中共有6個節(jié)點，其中節(jié)點1為平衡節(jié)點，視為無窮大母線，電壓穩(wěn)定在230kV。此部分程序可在MATLAB平臺中獨立寫成m文件，提高效率；3）初始化各節(jié)點電壓（幅值和相角）的存儲矩陣和各節(jié)點注入功率（實部和虛部）的存儲矩陣，以存放迭代過程中的中間值和輸出最終值；4）構造各節(jié)點的潮流方程，進入潮流計算程序；5）基于MATLAB平臺，得出計算結果，采用誤差精度為10e5。仿真器廣泛地使用了圖形和動畫功能，大大地增強了用戶對系統(tǒng)特性、問題和約束的理解，以便于用戶對系統(tǒng)進行維護。采用powerworld畫出的單相接線圖如圖5所示。說明兩種計算方法都是準確可行的。對于節(jié)點2，負荷ZL2=+，其中B12= ，B23= ， B24= = 即= + 對于節(jié)點3，負荷ZL3=+，其中B23=，B36== 即= + 對于節(jié)點4，負荷ZL4=+，其中B24=，B45== 即= + 對于節(jié)點5，負荷ZL5=+，其中B45=，B56== 即= + 對于節(jié)點6，不帶負荷，其中B36=，B56==即=j（）綜上，經(jīng)過上述簡化的電網(wǎng)接線圖如圖7所示：圖講過上述簡化后的模型表1圖7中各元件的阻抗標幺值Z12Z23Z24Z45Z36Z56++++++Z2Z3Z4Z5Z6Xd’+XT + + + + j（） 采用網(wǎng)絡變換法的具體步驟 消去節(jié)點4在圖7的基礎上通過星三角變換進行第一步簡化，可將節(jié)點4消去，如圖9示：圖8圖9參數(shù)計算如下：= + = + = + 將并聯(lián)的兩對阻抗合并后得：= + = + 另： 消去節(jié)點5在圖9的基礎上通過星三角變換進行第二步簡化，可將節(jié)點5消去，如圖11所示：圖10圖11參數(shù)計算如下：= + =