【正文】 47。電壓迭代次數(shù)曲線39。電壓39。迭代次數(shù)39。)。endenddisp(39。endfori=1:K Cs(i)=i。 endDS(i)=Si(p,q)+Sj(q,p)。else p=B1(i,2)。 fori=1:nl if B1(i,6)==0 p=B1(i,1)。各條支路的功率損耗DS（順序同您輸入B1時(shí)一樣）為：39。disp(Sj(q,p))。q=B1(i,1)。q=B1(i,2)。)。enddisp(39。 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))))。 else p=B1(i,2)。 fori=1:nl if B1(i,6)==0 p=B1(i,1)。各條支路的首端功率Sj（順序同您輸入B1時(shí)一樣）為：39。disp(S)。各節(jié)點(diǎn)的功率S（節(jié)點(diǎn)號從小到大排）為：39。endS(p)=E(p)*C(p)。for p=1:nC(p)=0。)。disp(39。)。disp(39。)。enddisp(39。for k=1:nE(k)=V(k)*cos(O(k))+V(k)*sin(O(k))*j。)。disp(39。)。 disp(39。)。endenddisp(39。ifkp~=0 K=K+1。 K=K+1。V(i)=V(i)DQ(L)。endendend L=0。endendelse IC1=i+1。for MZ=1:IC4 I=IC4+1MZ。ifi==nafor LZ=2:i L=i+2LZ。if ICT3~=0 L=0。 DET=abs(DQ1(i))。endDQ1(i)=Q(i)V(i)*C(i)。L=L+1。 L=0。endendendelsekp=0。 DET=abs(DQ1(i))。endDQ1(i)=Q(i)V(i)*C(i)。L=L+1。L=0。endendendendfori=2:n O(i)=O(i)DP(i)。for MZ=2:IC4 I=IC4+2MZ。 endelsefor LZ=3:i L=i+3LZ。ifi~=n IC1=i+1。endendendNp(K)=ICT2。 DET=abs(DP1(i))。endDP1(i)=P(i)V(i)*C(i)。 fori=1:nifi~=isbC(i)=0。while ICT2~=0|ICT3~=0 ICT2=0。DET=0。kq=1。ICT1=0。for k=i+1:nafor j1=i+1:naA(k,j1)=A(k,j1)A(k,i)*A(i,j1)。for j1=k:naA(i,j1)=A(i,j1)./A(i,i)。endendendendfori=1:na ifi==naA(i,i)=1./A(i,i)。for j1=1:nif B2(j1,6)==2 k=k+1。fori=1:nif B2(i,6)==2 p=p+1。endendendendp=0。endB(i,i)=1./B(i,i)。else IC1=i+1。 O(i)=atan(f(i)./e(i))。V(i)=B2(i,4)。fori=1:n e(i)=real(B2(i,3))。endP=real(S)。fori=1:n S(i)=B2(i,1)B2(i,2)。B=imag(YI)。YI(p,p)=YI(p,p)+1./B1(i,3)。YI(q,q)=YI(q,q)+1./B1(i,3)。YI(q,p)=YI(p,q)。YI(p,q)=YI(p,q)1./B1(i,3)。q=B1(i,1)。q=B1(i,2)。Y(p,p)=1./X(i,2)。fori=1:n if X(i,2)~=0。V=zeros(1,n)。e=zeros(1,n)。 Y=zeros(n)。請輸入PQ節(jié)點(diǎn)數(shù)na=39。)。 %輸入B2X=input(39。請輸入由支路參數(shù)形成的矩陣：B2=39。)。 B1=input(39。請輸入誤差精度：pr=39。)。 isb=input(39。請輸入支路數(shù)：nl=39。)。主程序n=input(39。5 0。3 0。X=[1 0。0 +*i 1 0 0 2。0 +*i 1 0 0 2。B2=[0 0 1 0 1。4 5 +i** 1 0。2 4 +i* 1i* 1 0。參考文獻(xiàn)[1] 《電力系統(tǒng)分析》（上下冊）華中科技大學(xué)出版[2] 《發(fā)電廠電氣部分》中國電力出版社[3] PowerWorld 16使用手冊[4] 基于Matlab/Simulink的系統(tǒng)仿真技術(shù)與應(yīng)用附錄：Matlab編制PQ分解法潮流計(jì)算程序：本文中的實(shí)例數(shù)據(jù)如下：節(jié)點(diǎn)數(shù)為6；支路數(shù)為6；平衡母線節(jié)點(diǎn)號為1；；PQ節(jié)點(diǎn)數(shù)為4；B1=[1 2 +i** 1 0。計(jì)算結(jié)果如表18：表1定量分析TR對穩(wěn)定計(jì)算結(jié)果的影響Tr012(s)圖5(s)隨勵(lì)磁參數(shù)Tr的變化特性從圖30中可知，當(dāng)勵(lì)磁參數(shù)Tr變大，極限切除角也變大。 定量分析TR對暫態(tài)穩(wěn)定的影響三階模型的勵(lì)磁參數(shù)較多，下面以參數(shù)Tr為例，分析定量勵(lì)磁參數(shù)對暫態(tài)穩(wěn)定的影響。圖50、(s)隨TJ變化特性圖5勵(lì)磁環(huán)節(jié)模型原理框圖 定性分析勵(lì)磁參數(shù)對暫態(tài)穩(wěn)定的影響為了分析勵(lì)磁參數(shù)對電力系統(tǒng)穩(wěn)定性的影響，本文基于發(fā)電機(jī)的三階模型，通過單一變量法（每次僅改變一個(gè)參數(shù)）仿真不同勵(lì)磁參數(shù)的暫態(tài)過程，其中故障選取故障方案1：節(jié)點(diǎn)3發(fā)生三相對稱短路。當(dāng)TJ足夠大后，再增大TJ時(shí)，極限切除角變化幅度不大。在PowerWorld軟件里，發(fā)電機(jī)采用三階模型，勵(lì)磁系統(tǒng)采用BPA_EA模型，TJ分別取50時(shí)，分別求出系統(tǒng)的極限切除角，其中故障選取故障方案1：節(jié)點(diǎn)3發(fā)生三相對稱短路。從理論上說，增大TJ可以減少發(fā)電機(jī)受到擾動(dòng)后轉(zhuǎn)子相對動(dòng)能的變化量，即有利于提高暫態(tài)穩(wěn)定。換言之，采用三階模型的發(fā)電機(jī)的暫態(tài)穩(wěn)定性更好。表16為兩種不同的發(fā)電機(jī)模型在不同故障方案下的暫態(tài)過渡時(shí)間統(tǒng)計(jì)；表1不同的發(fā)電機(jī)模型在不同故障方案下的暫態(tài)過渡時(shí)間發(fā)電機(jī)模型故障方案一故障方案二故障切除時(shí)間（s）二階模型暫態(tài)過渡時(shí)間（s）454748424344三階模型暫態(tài)過渡時(shí)間（s）141516131415從表中可以看出，對于同一模型，故障切除時(shí)間越大，暫態(tài)過渡時(shí)間越大。但總的來說兩種發(fā)電機(jī)模型的極限切除時(shí)間相差并不大。lim=。所得搖擺曲線如下，此時(shí)系統(tǒng)暫態(tài)不穩(wěn)定。所得搖擺曲線如下，此時(shí)系統(tǒng)暫態(tài)穩(wěn)定。所得搖擺曲線如下，此時(shí)系統(tǒng)暫態(tài)不穩(wěn)定。則采用PowerWorld軟件仿真結(jié)果如下：（1）5s時(shí)線路L23中點(diǎn)處發(fā)生三相對稱短路。lim=。圖42（4）在5s時(shí)節(jié)點(diǎn)3發(fā)生三相對稱短路，所得搖擺曲線如下，即此時(shí)系統(tǒng)暫態(tài)不穩(wěn)定。圖40（2）在5s時(shí)節(jié)點(diǎn)3發(fā)生三相對稱短路，所得搖擺曲線如下，即此時(shí)系統(tǒng)暫態(tài)不穩(wěn)定。7 基于PowerWorld，分析三階模型的發(fā)電機(jī)的暫態(tài)過程 故障方案一：節(jié)點(diǎn)3發(fā)生三相對稱短路由于發(fā)電機(jī)采用三階模型，在PowerWorld軟件中給接在節(jié)點(diǎn)6的發(fā)電機(jī)添加GENTRA模型，勵(lì)磁系統(tǒng)設(shè)為BPA_EA。從總體上看，三種方法求出的結(jié)果較為接近。 故障方案二由第5部分知，L23中點(diǎn)處發(fā)生三相對稱短路，當(dāng)故障切除后。從總體上看，三種方法求出的結(jié)果較為接近。圖3故障后功角隨時(shí)間的變化曲線6 采用等面積法求解極限切除角 故障方案一由第5部分知，當(dāng)節(jié)點(diǎn)3發(fā)生三相對稱短路，當(dāng)故障切除后。=%5%（6）通過求解故障后功角隨時(shí)間的變化曲線。圖36綜合上面4點(diǎn)的仿真結(jié)果，=。圖34（3）在事件窗口添加：1s時(shí)線路L23中點(diǎn)處發(fā)生三相對稱短路，得出的搖擺曲線如下，此時(shí)系統(tǒng)暫態(tài)穩(wěn)定。（1）在事件窗口添加：1s時(shí)線路L23中點(diǎn)處發(fā)生三相對稱短路，得出的搖擺曲線如下，此時(shí)系統(tǒng)暫態(tài)穩(wěn)定。圖3tc=同理。圖tc=，修改相應(yīng)的Matlab程序，得到的搖擺曲線如下，即此時(shí)發(fā)電暫態(tài)不穩(wěn)定。b39。r39。x2]。t2]。,[tc,20],x02,h_opt,a,b,c)。%短路后的P3[t2,x2]=ode45(39。b=。%求解第一個(gè)區(qū)間[0,tc]x02= x1 (end,:)39。si