【正文】 tlabPowerWorld等面積法（s）極限切除角（度）124從表14可知，采用等面積法求出的極限切除角是最小的，其次是Matlab仿真法，最大的是PowerWorld仿真法。如果把等面積法的結(jié)果視為實際值，那Matlab仿真法的相對誤差為PowerWorld仿真法的相對誤差為兩者的相對誤差皆小于5%，所以采用Matlab與PowerWorld仿真計算的結(jié)果精確度較高。具體參數(shù)設(shè)置分別如圖339所示圖3發(fā)電機GENTRA模型參數(shù)設(shè)置圖3勵磁系統(tǒng)BPA_EA模型參數(shù)設(shè)置（1）在5s節(jié)點3時發(fā)生三相對稱短路，所得搖擺曲線如下，即此時系統(tǒng)暫態(tài)穩(wěn)定。圖41（3）在5s時節(jié)點3發(fā)生三相對稱短路，所得搖擺曲線如下，即此時系統(tǒng)暫態(tài)穩(wěn)定。圖43綜上，取tc。 故障方案二：線路L23中點處發(fā)生三相對稱短路發(fā)電機與勵磁系統(tǒng)選取的模型與參數(shù)設(shè)置不變，即圖21。所得搖擺曲線如下，此時系統(tǒng)暫態(tài)穩(wěn)定。圖44（2）5s時線路L23中點處發(fā)生三相對稱短路。圖45（3）5s時線路L23中點處發(fā)生三相對稱短路。圖46（4）5s時線路L23中點處發(fā)生三相對稱短路。圖47綜上，取tc。8 發(fā)電機模型及勵磁調(diào)節(jié)參數(shù)對穩(wěn)定計算結(jié)果的影響 發(fā)電機模型對穩(wěn)定計算結(jié)果的影響 不同發(fā)電機模型對極限切除時間的影響發(fā)電機分別采用二階模型和三階模型時，極限切除時間的計算結(jié)果如表15所示：表1不同發(fā)電機模型的暫定分析模型種類故障方案一故障方案二(s)二階模型三階模型從表15可知，采用三階模型的發(fā)電機的極限切除時間較小。 不同發(fā)電機模型對暫定過渡時間的影響暫定過渡時間是指從發(fā)生故障時刻算起到功角穩(wěn)定于一個定值后的這一段時間，暫態(tài)過渡時間越短意味功角穩(wěn)定地越快，即系統(tǒng)的暫態(tài)穩(wěn)定性越好。對不同的模型，三階模型暫態(tài)過渡時間比二階模型小很多，采用三階模型的發(fā)電機的暫態(tài)過程更短。圖4采用二階模型時暫定過渡時間較大（故障方案一、tc=）圖4采用三階模型時暫定過渡時間較?。ü收戏桨敢弧c=） TJ對穩(wěn)定計算結(jié)果的影響慣性時間常數(shù)TJ是反映發(fā)電機轉(zhuǎn)子機械慣性的重要參數(shù)，由TJ的定義可知，它是轉(zhuǎn)子在額定轉(zhuǎn)速下的動能的兩倍除以基準(zhǔn)功率。下面用PowerWorld軟件來分析TJ對穩(wěn)定計算結(jié)果的影響。計算結(jié)果如表16所示：表1TJ取不同值時系統(tǒng)的極限切除角TJ20406080100(s)(s)隨TJ變化特性如圖30所示，可以發(fā)現(xiàn)，當(dāng)TJ增大時，極限切除角也變大。換言之，當(dāng)發(fā)電機的TJ較小時，通過增大TJ的方式來提高系統(tǒng)的暫態(tài)穩(wěn)定性可以獲得較好的效果。PowerWorld的定性分析如下表表1勵磁參數(shù)對穩(wěn)定計算結(jié)果的影響勵磁系統(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)定。在PowerWorld軟件的勵磁模型對話框中，先后輸入TR=0、120這五個數(shù)據(jù)，然后分別采用試湊法求出極限切除時間，其中故障選取故障方案1：節(jié)點3發(fā)生三相對稱短路。此外，當(dāng)Tr較小時，極限切除角下降較為明顯；當(dāng)Tr，極限切除角基本不變。2 3 +i* 1i* 1 0。3 6 +i* 1i* 1 0。5 6 +i** 1 0]。0 +i* 1 0 0 2。0 +*i 1 0 0 2。 0 1 0 3]。2 0。4 0。6 0]。請輸入節(jié)點數(shù)：n=39。nl=input(39。)。請輸入平衡母線節(jié)點號：isb=39。 pr=input(39。)。請輸入由支路參數(shù)形成的矩陣：B1=39。 %輸入B1 B2=input(39。)。請輸入由節(jié)點號及其對地阻抗形成的矩陣：X=39。%輸入Xna=input(39。)。YI=zeros(n)。f=zeros(1,n)。O=zeros(1,n)。 p=X(i,1)。endendfori=1:nlif B1(i,6)==0 p=B1(i,1)。else p=B1(i,2)。end Y(p,q)=Y(p,q)1./(B1(i,3)*B1(i,5))。Y(q,p)=Y(p,q)。 Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2。 Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2。end %求導(dǎo)納矩陣G=real(Y)。BI=imag(Y)。BI(i,i)=BI(i,i)+B2(i,5)。Q=imag(S)。f(i)=imag(B2(i,3))。endfori=1:nif B2(i,6)==2 V(i)=sqrt(e(i)^2+f(i)^2)。 endendfori=2:n ifi==n B(i,i)=1./B(i,i)。 for j1=IC1:n B(i,j1)=B(i,j1)./B(i,i)。 for k=i+1:n for j1=i+1:n B(k,j1)=B(k,j1)B(k,i)*B(i,j1)。q=0。k=0。 A(p,k)=BI(i,j1)。 else k=i+1。endA(i,i)=1./A(i,i)。endendendendICT2=1。kp=1。K=1。ICT3=1。ICT3=0。for k=1:n C(i)=C(i)+V(k)*(G(i,k)*cos(O(i)O(k))+BI(i,k)*sin(O(i)O(k)))。DP(i)=DP1(i)./V(i)。if DET=pr ICT2=ICT2+1。if ICT2~=0fori=2:nDP(i)=B(i,i)*DP(i)。for k=IC1:nDP(k)=DP(k)B(k,i)*DP(i)。 IC4=L1。DP(I)=DP(I)B(I,L)*DP(L)。endkq=1。fori=1:nif B2(i,6)==2C(i)=0。for k=1:n C(i)=C(i)+V(k)*(G(i,k)*sin(O(i)O(k))BI(i,k)*cos(O(i)O(k)))。DQ(L)=DQ1(i)./V(i)。if DET =pr ICT3=ICT3+1。 ifkq~=0。fori=1:n if B2(i,6)==2C(i)=0。for k=1:n C(i)=C(i)+V(k)*(G(i,k)*sin(O(i)O(k))BI(i,k)*cos(O(i)O(k)))。DQ(L)=DQ1(i)./V(i)。endendendendNq(K)=ICT3。fori=1:naDQ(i)=A(i,i)*DQ(i)。 IC4=L1。DQ(I)=DQ(I)A(I,L)*DQ(L)。for k=IC1:naDQ(k)=DQ(k)A(k,i)*DQ(i)。fori=1:nif B2(i,6)==2 L=L+1。endendkp=1。elsekq=0。endendfori=1:nDy(K1,i)=V(i)。迭代次數(shù)39。disp(K)。每次沒有達(dá)到精度要求的有功功率個數(shù)為39。disp(Np)。每次沒有達(dá)到精度要求的無功功率個數(shù)為39。disp(Nq)。O(k)=O(k)*180./pi。各節(jié)點的電壓標(biāo)幺值E為（節(jié)點號從小到大排）：39。disp(E)。各節(jié)點的電壓V大?。ü?jié)點號從小到大排）為：39。disp(V)。各節(jié)點的電壓相角O（節(jié)點號從小到大排）為：39。disp(O)。for q=1:n C(p)=C(p)+conj(Y(p,q))*conj(E(q))。enddisp(39。)。disp(39。)。q=B1(i,2)。q=B1(i,1)。disp(Si(p,q))。各條支路的末端功率Sj（順序同您輸入B1時一樣）為：39。 fori=1:nlif B1(i,6)==0 p=B1(i,1)。else p=B1(i,2)。 end Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))conj(E(p)))*conj(1./(B1(i,3)*B1(i,5))))。 enddisp(39。)。q=B1(i,2)。q=B1(i,1)。 disp(DS(i))。 for j=1:n Dy(K,j)=Dy(K1,j)。以下是每次迭代后各節(jié)點的電壓值（如圖所示）39。 plot(Cs,Dy),xlabel(39。),ylabel(39。),title(39