【正文】 磁參數(shù)的暫態(tài)過程，其中故障選取故障方案1：節(jié)點3發(fā)生三相對稱短路。計算結(jié)果如表18：表1定量分析TR對穩(wěn)定計算結(jié)果的影響Tr012(s)圖5(s)隨勵磁參數(shù)Tr的變化特性從圖30中可知，當勵磁參數(shù)Tr變大，極限切除角也變大。2 4 +i* 1i* 1 0。B2=[0 0 1 0 1。0 +*i 1 0 0 2。3 0。主程序n=input(39。請輸入支路數(shù)：nl=39。)。 B1=input(39。請輸入由支路參數(shù)形成的矩陣：B2=39。)。 Y=zeros(n)。V=zeros(1,n)。Y(p,p)=1./X(i,2)。q=B1(i,1)。YI(q,p)=YI(p,q)。YI(p,p)=YI(p,p)+1./B1(i,3)。fori=1:n S(i)=B2(i,1)B2(i,2)。fori=1:n e(i)=real(B2(i,3))。 O(i)=atan(f(i)./e(i))。endB(i,i)=1./B(i,i)。fori=1:nif B2(i,6)==2 p=p+1。endendendendfori=1:na ifi==naA(i,i)=1./A(i,i)。for k=i+1:nafor j1=i+1:naA(k,j1)=A(k,j1)A(k,i)*A(i,j1)。kq=1。while ICT2~=0|ICT3~=0 ICT2=0。endDP1(i)=P(i)V(i)*C(i)。endendendNp(K)=ICT2。 endelsefor LZ=3:i L=i+3LZ。endendendendfori=2:n O(i)=O(i)DP(i)。L=L+1。 DET=abs(DQ1(i))。 L=0。endDQ1(i)=Q(i)V(i)*C(i)。if ICT3~=0 L=0。for MZ=1:IC4 I=IC4+1MZ。endendend L=0。 K=K+1。endenddisp(39。 disp(39。disp(39。for k=1:nE(k)=V(k)*cos(O(k))+V(k)*sin(O(k))*j。)。)。)。endS(p)=E(p)*C(p)。disp(S)。 fori=1:nl if B1(i,6)==0 p=B1(i,1)。 end Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))conj(E(q)))*conj(1./(B1(i,3)*B1(i,5))))。)。q=B1(i,1)。各條支路的功率損耗DS（順序同您輸入B1時一樣）為：39。else p=B1(i,2)。endfori=1:K Cs(i)=i。)。電壓39。47。電壓迭代次數(shù)曲線39。迭代次數(shù)39。endenddisp(39。 endDS(i)=Si(p,q)+Sj(q,p)。 fori=1:nl if B1(i,6)==0 p=B1(i,1)。disp(Sj(q,p))。q=B1(i,2)。enddisp(39。 else p=B1(i,2)。各條支路的首端功率Sj（順序同您輸入B1時一樣）為：39。各節(jié)點的功率S（節(jié)點號從小到大排）為：39。for p=1:nC(p)=0。disp(39。disp(39。enddisp(39。)。)。)。ifkp~=0 K=K+1。V(i)=V(i)DQ(L)。endendelse IC1=i+1。ifi==nafor LZ=2:i L=i+2LZ。 DET=abs(DQ1(i))。L=L+1。endendendelsekp=0。endDQ1(i)=Q(i)V(i)*C(i)。L=0。for MZ=2:IC4 I=IC4+2MZ。ifi~=n IC1=i+1。 DET=abs(DP1(i))。 fori=1:nifi~=isbC(i)=0。DET=0。ICT1=0。for j1=k:naA(i,j1)=A(i,j1)./A(i,i)。for j1=1:nif B2(j1,6)==2 k=k+1。endendendendp=0。else IC1=i+1。V(i)=B2(i,4)。endP=real(S)。B=imag(YI)。YI(q,q)=YI(q,q)+1./B1(i,3)。YI(p,q)=YI(p,q)1./B1(i,3)。q=B1(i,2)。fori=1:n if X(i,2)~=0。e=zeros(1,n)。請輸入PQ節(jié)點數(shù)na=39。 %輸入B2X=input(39。)。請輸入誤差精度：pr=39。 isb=input(39。)。5 0。X=[1 0。0 +*i 1 0 0 2。4 5 +i** 1 0。參考文獻[1] 《電力系統(tǒng)分析》（上下冊）華中科技大學(xué)出版[2] 《發(fā)電廠電氣部分》中國電力出版社[3] PowerWorld 16使用手冊[4] 基于Matlab/Simulink的系統(tǒng)仿真技術(shù)與應(yīng)用附錄：Matlab編制PQ分解法潮流計算程序：本文中的實例數(shù)據(jù)如下：節(jié)點數(shù)為6；支路數(shù)為6；平衡母線節(jié)點號為1；；PQ節(jié)點數(shù)為4；B1=[1 2 +i** 1 0。 定量分析TR對暫態(tài)穩(wěn)定的影響三階模型的勵磁參數(shù)較多，下面以參數(shù)Tr為例，分析定量勵磁參數(shù)對暫態(tài)穩(wěn)定的影響。當TJ足夠大后，再增大TJ時，極限切除角變化幅度不大。從理論上說，增大TJ可以減少發(fā)電機受到擾動后轉(zhuǎn)子相對動能的變化量，即有利于提高暫態(tài)穩(wěn)定。表16為兩種不同的發(fā)電機模型在不同故障方案下的暫態(tài)過渡時間統(tǒng)計；表1不同的發(fā)電機模型在不同故障方案下的暫態(tài)過渡時間發(fā)電機模型故障方案一故障方案二故障切除時間（s）二階模型暫態(tài)過渡時間（s）454748424344三階模型暫態(tài)過渡時間（s）141516131415從表中可以看出，對于同一模型，故障切除時間越大，暫態(tài)過渡時間越大。lim=。所得搖擺曲線如下，此時系統(tǒng)暫態(tài)穩(wěn)定。lim=。圖40（2）在5s時節(jié)點3發(fā)生三相對稱短路，所得搖擺曲線如下，即此時系統(tǒng)暫態(tài)不穩(wěn)定。從總體上看，三種方法求出的結(jié)果較為接近。從總體上看，三種方法求出的結(jié)果較為接近。=%5%（6）通過求解故障后功角隨時間的變化曲線。圖34（3）在事件窗口添加：1s時線路L23中點處發(fā)生三相對稱短路，得出的搖擺曲線如下，此時系統(tǒng)暫態(tài)穩(wěn)定。圖3tc=同理。b39。x2]。,[tc,20],x02,h_opt,a,b,c)。b=。simb39。%給定微分方程組的初值a=。tc=。代入式13后得：（18）當節(jié)點3的故障切除后，通過網(wǎng)絡(luò)變換法同理可得，轉(zhuǎn)移阻抗Z17= + =，輸入阻抗Z77= + =。圖2故障后功角隨時間的變化曲線 選定故障2：線路L23中點處發(fā)生三相對稱短路故障說明：因為線路L23是雙回路線路，如果單純是單回線路發(fā)生對稱短路，由于還有另外一回線路正常工作，并聯(lián)后的轉(zhuǎn)移阻抗并不會很大，即故障后的電磁功率不會明顯減小。由于采用Matlab仿真時需要計算各種狀態(tài)下的轉(zhuǎn)移阻抗與輸入阻抗，而且為復(fù)數(shù)運算，所以存在一點的誤差。圖26（3）在事件窗口添加：1s時節(jié)點3發(fā)生三相對稱短路，得出的搖擺曲線如下，此時系統(tǒng)暫態(tài)穩(wěn)定。點擊圖中交點后，輸出x =，y =，即==*180/pi=。%求解區(qū)間為[0,]plot(t,x(:,1),39。c=*pi/180。x0=[*pi/180。其中的Matlab求解程序如下：（子程序與穩(wěn)定計算程序一樣）functionxdot=simb(t,x,flag,a,b,c)xdot=[(x(2)1)*。)%畫出搖擺曲線運行程序后得搖擺曲線如下，功角逐漸穩(wěn)定，可知發(fā)電機暫態(tài)穩(wěn)定。plot(t,x(:,1),39。%求解第二個區(qū)間[tc,20s]t=[t1。c=*pi/180。,[0,tc],x0,h_opt,a,b,c)。b=。%給定角速度初值為1，h_opt=odeset。假設(shè)在5s時刻節(jié)點3發(fā)生了三相對稱短路，則采用Matlab編制的程序如下所示：子函數(shù):functionxdot=simb(t,x,flag,a,b,c)xdot=[(x(2)1)*。代入式12后得：（15）當節(jié)點3的故障切除后，通過網(wǎng)絡(luò)變換法同理可得，轉(zhuǎn)移阻抗Z17= + =，輸入阻抗Z77= + =。