【正文】 )。),ylabel(39。以下是每次迭代后各節(jié)點(diǎn)的電壓值（如圖所示）39。 disp(DS(i))。q=B1(i,2)。 enddisp(39。else p=B1(i,2)。各條支路的末端功率Sj（順序同您輸入B1時(shí)一樣）為：39。q=B1(i,1)。)。)。for q=1:n C(p)=C(p)+conj(Y(p,q))*conj(E(q))。各節(jié)點(diǎn)的電壓相角O（節(jié)點(diǎn)號從小到大排）為：39。各節(jié)點(diǎn)的電壓V大?。ü?jié)點(diǎn)號從小到大排）為：39。各節(jié)點(diǎn)的電壓標(biāo)幺值E為（節(jié)點(diǎn)號從小到大排）：39。disp(Nq)。disp(Np)。disp(K)。endendfori=1:nDy(K1,i)=V(i)。endendkp=1。for k=IC1:naDQ(k)=DQ(k)A(k,i)*DQ(i)。 IC4=L1。endendendendNq(K)=ICT3。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)))。 ifkq~=0。DQ(L)=DQ1(i)./V(i)。fori=1:nif B2(i,6)==2C(i)=0。DP(I)=DP(I)B(I,L)*DP(L)。for k=IC1:nDP(k)=DP(k)B(k,i)*DP(i)。if DET=pr ICT2=ICT2+1。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)))。ICT3=1。kp=1。endA(i,i)=1./A(i,i)。 A(p,k)=BI(i,j1)。q=0。 for j1=IC1:n B(i,j1)=B(i,j1)./B(i,i)。endfori=1:nif B2(i,6)==2 V(i)=sqrt(e(i)^2+f(i)^2)。Q=imag(S)。BI=imag(Y)。 Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2。Y(q,p)=Y(p,q)。else p=B1(i,2)。 p=X(i,1)。f=zeros(1,n)。)。請輸入由節(jié)點(diǎn)號及其對地阻抗形成的矩陣：X=39。 %輸入B1 B2=input(39。)。請輸入平衡母線節(jié)點(diǎn)號：isb=39。nl=input(39。6 0]。2 0。0 +*i 1 0 0 2。5 6 +i** 1 0]。2 3 +i* 1i* 1 0。在PowerWorld軟件的勵(lì)磁模型對話框中，先后輸入TR=0、120這五個(gè)數(shù)據(jù)，然后分別采用試湊法求出極限切除時(shí)間，其中故障選取故障方案1：節(jié)點(diǎn)3發(fā)生三相對稱短路。換言之，當(dāng)發(fā)電機(jī)的TJ較小時(shí)，通過增大TJ的方式來提高系統(tǒng)的暫態(tài)穩(wěn)定性可以獲得較好的效果。下面用PowerWorld軟件來分析TJ對穩(wěn)定計(jì)算結(jié)果的影響。對不同的模型，三階模型暫態(tài)過渡時(shí)間比二階模型小很多，采用三階模型的發(fā)電機(jī)的暫態(tài)過程更短。8 發(fā)電機(jī)模型及勵(lì)磁調(diào)節(jié)參數(shù)對穩(wěn)定計(jì)算結(jié)果的影響 發(fā)電機(jī)模型對穩(wěn)定計(jì)算結(jié)果的影響 不同發(fā)電機(jī)模型對極限切除時(shí)間的影響發(fā)電機(jī)分別采用二階模型和三階模型時(shí)，極限切除時(shí)間的計(jì)算結(jié)果如表15所示：表1不同發(fā)電機(jī)模型的暫定分析模型種類故障方案一故障方案二(s)二階模型三階模型從表15可知，采用三階模型的發(fā)電機(jī)的極限切除時(shí)間較小。圖46（4）5s時(shí)線路L23中點(diǎn)處發(fā)生三相對稱短路。圖44（2）5s時(shí)線路L23中點(diǎn)處發(fā)生三相對稱短路。 故障方案二：線路L23中點(diǎn)處發(fā)生三相對稱短路發(fā)電機(jī)與勵(lì)磁系統(tǒng)選取的模型與參數(shù)設(shè)置不變，即圖21。圖41（3）在5s時(shí)節(jié)點(diǎn)3發(fā)生三相對稱短路，所得搖擺曲線如下，即此時(shí)系統(tǒng)暫態(tài)穩(wěn)定。如果把等面積法的結(jié)果視為實(shí)際值，那Matlab仿真法的相對誤差為PowerWorld仿真法的相對誤差為兩者的相對誤差皆小于5%，所以采用Matlab與PowerWorld仿真計(jì)算的結(jié)果精確度較高。如果把等面積法的結(jié)果視為實(shí)際值，那Matlab仿真法的相對誤差為PowerWorld仿真法的相對誤差為兩者的相對誤差皆小于5%，所以采用Matlab與PowerWorld仿真計(jì)算的結(jié)果精確度較高。故障后功角隨時(shí)間的變化曲線如圖29所示，從圖中可以讀出極限切除角=。圖35（4）在事件窗口添加：1s時(shí)線路L23中點(diǎn)處發(fā)生三相對稱短路，得出的搖擺曲線如下，此時(shí)系統(tǒng)暫態(tài)不穩(wěn)定。如圖32所示：圖3故障時(shí)功角隨時(shí)間變化的特性 采用 PowerWorld軟件分析結(jié)果發(fā)電機(jī)仍是采用經(jīng)典二階模型（GENCLS）。)%畫出搖擺曲線運(yùn)行程序后得搖擺曲線如下，功角逐漸穩(wěn)定，可知發(fā)電機(jī)暫態(tài)穩(wěn)定。plot(t,x(:,1),39。%求解第二個(gè)區(qū)間[tc,20s]t=[t1。c=*pi/180。,[0,tc],x0,h_opt,a,b,c)。b=。%給定角速度初值為1，h_opt=odeset。代入式12后得：（19） Matlab編制穩(wěn)定計(jì)算程序假設(shè)在5s時(shí)刻節(jié)點(diǎn)3發(fā)生了三相對稱短路，則采用Matlab編制的程序如下所示：子函數(shù):functionxdot=simb(t,x,flag,a,b,c)xdot=[(x(2)1)*。換言之，單回路發(fā)生故障時(shí)，對系統(tǒng)的沖擊有限。此外兩個(gè)軟件采用的算法不同，對精確的要求也不一樣，所以也會造成一定的計(jì)算差距。圖27（4）在事件窗口添加：1s時(shí)節(jié)點(diǎn)3發(fā)生三相對稱短路，得出的搖擺曲線如下，此時(shí)系統(tǒng)暫態(tài)不穩(wěn)定。 采用 PowerWorld軟件分析結(jié)果在PowerWorld軟件里給發(fā)電機(jī)建立暫態(tài)穩(wěn)定模型，其中平衡節(jié)點(diǎn)用一臺大發(fā)電機(jī)代替，選定GENPWTwoAxis模型，把其中的H值設(shè)得非常大（），其他都用默認(rèn)參數(shù)。r39。%短路后的電磁功率[t,x]=ode45(39。w0]。 1/*(*sin(x(1)+c)10*(x(2)1))]。圖1tc=，得到的搖擺曲線如下，此時(shí)發(fā)電暫態(tài)不穩(wěn)定。r39。t2]。%短路后的P3，如式15[t2,x2]=ode45(39。%求解第一個(gè)區(qū)間[0,tc]x02= x1 (end,:)39。c=*pi/180。x0=[*pi/180。 1/*(*sin(x(1)+c)10*(x(2)1))]。代入式12后得：（16） Matlab編制穩(wěn)定計(jì)算程序轉(zhuǎn)子運(yùn)動方程如下（17）而Pe的表達(dá)式分別為式1115所示，故可以采用Matlab軟件求解這組微分方程。由文中第4部分知，E’=，V=，穩(wěn)定運(yùn)行時(shí)的轉(zhuǎn)移阻抗Z17= +=，輸入阻抗Z77= +j =。則有（12）有功功率輸入特性按式12來計(jì)算，其中，為的阻抗角，為的阻抗角。在網(wǎng)絡(luò)變換中，用電壓為E’、內(nèi)阻為Xd’的電壓源來等效原來的三臺機(jī)組。4 網(wǎng)絡(luò)變換法求解輸入阻抗與轉(zhuǎn)移阻抗 負(fù)荷、發(fā)電機(jī)、變壓器、線路等效電容的處理 負(fù)荷的處理由題目要求，負(fù)荷采用恒阻抗模型，計(jì)算公式為（9）其中，U是負(fù)荷所在節(jié)點(diǎn)的電壓，是負(fù)荷的共軛值。說明兩種計(jì)算方法都是準(zhǔn)確可行的。圖采用powerworld畫出的單相接線圖 采用powerworld仿真軟件的潮流計(jì)算結(jié)果在powerworld軟件的運(yùn)行模型下，進(jìn)行潮流計(jì)算，電網(wǎng)的潮流分布如圖6圖電網(wǎng)的潮流分布示意圖（1）各節(jié)點(diǎn)電壓如表6所示表各節(jié)點(diǎn)電壓節(jié)點(diǎn)編號基準(zhǔn)電壓kV標(biāo)幺電壓實(shí)際電壓(kV)相角 (度)1220022203220422052206220（2）各發(fā)電機(jī)狀態(tài)如表7所示表各發(fā)電機(jī)運(yùn)行狀態(tài)節(jié)點(diǎn)編號發(fā)電機(jī)編號有功出力（MW）無功出力（MW）設(shè)定電壓11616263（3）各支路狀態(tài)如表8所示表各支路狀態(tài)首端節(jié)點(diǎn)編號末端節(jié)點(diǎn)編號首端有功功率（MW）首端無功功率（Mvar）首端視在功率（MVA）有功損耗（MW）無功損耗（Mvar）122324364556150（4）電網(wǎng)網(wǎng)損率的計(jì)算網(wǎng)損率按以下公式計(jì)算： 其中，發(fā)電機(jī)輸出功率包括發(fā)電廠和平衡節(jié)點(diǎn)的功率，發(fā)電機(jī)輸出功率和負(fù)荷均只取有功分量。具體參數(shù)設(shè)置如圖4所示。它基本的工具包括經(jīng)濟(jì)調(diào)度、區(qū)域功率經(jīng)濟(jì)分配分析、功率傳輸分配因子計(jì)算（PTDF）、短路分析以及事故分析等功能的工具。它能夠有效地計(jì)算高達(dá)10,0000個(gè)節(jié)點(diǎn)的電力網(wǎng)絡(luò)，因此當(dāng)它作為一個(gè)獨(dú)立的潮流分析軟件包時(shí)，性非常實(shí)用。由結(jié)果可知，共用7次迭代便能滿足精度要求，； MATLAB平臺潮流計(jì)算結(jié)果分析采用MATLAB平臺計(jì)算電網(wǎng)潮流的輸出結(jié)果如下所示，（1）各節(jié)點(diǎn)電壓如表3所示：表各節(jié)點(diǎn)電壓節(jié)點(diǎn)編號基準(zhǔn)電壓kV標(biāo)幺電壓實(shí)際電壓(kV)相角 (度)12200