【正文】 plot(Cs,Dy),xlabel(39。 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))))。disp(39。disp(E)。迭代次數(shù)39。fori=1:naDQ(i)=A(i,i)*DQ(i)。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)))。DP(i)=DP1(i)./V(i)。 else k=i+1。f(i)=imag(B2(i,3))。end Y(p,q)=Y(p,q)1./(B1(i,3)*B1(i,5))。%輸入Xna=input(39。)。0 +i* 1 0 0 2。計(jì)算結(jié)果如表16所示：表1TJ取不同值時(shí)系統(tǒng)的極限切除角TJ20406080100(s)(s)隨TJ變化特性如圖30所示，可以發(fā)現(xiàn)，當(dāng)TJ增大時(shí)，極限切除角也變大。圖45（3）5s時(shí)線路L23中點(diǎn)處發(fā)生三相對(duì)稱短路。臨界角為，滿足式，即（22）求得設(shè)極限切除角為，由等面積定則有（23）將參數(shù)代入后解方程得，極限切除角幾種方法求出來(lái)的結(jié)果統(tǒng)計(jì)如下：表1幾種方法求出來(lái)的結(jié)果統(tǒng)計(jì)計(jì)算方法MatlabPowerWorld等面積法（s）極限切除角（度）124從表14可知，采用等面積法求出的極限切除角是最小的，其次是Matlab仿真法，最大的是PowerWorld仿真法。=。%給第二個(gè)計(jì)算區(qū)間賦初值a=。 參數(shù)的求取線路L23中點(diǎn)處發(fā)生三相對(duì)稱短路，通過(guò)網(wǎng)絡(luò)變換法同理可得，轉(zhuǎn)移阻抗Z17==，輸入阻抗Z77= + =。在主窗口中輸入：[x,y]=ginput(1)。圖tc=（穩(wěn)定）圖2tc=（不穩(wěn)定）（2）采用Matlab求極限切除角通過(guò)求解故障時(shí)發(fā)電機(jī)轉(zhuǎn)子運(yùn)動(dòng)方程來(lái)確定功角隨時(shí)間變化的特性，進(jìn)而可以求出極限切除角。b=。此外，為了加快功角的收斂速度，添加了阻尼部分，阻尼系數(shù)D設(shè)為10。則歸算到全網(wǎng)基準(zhǔn)后的發(fā)電機(jī)參數(shù)可按式10來(lái)計(jì)算： （10）代入?yún)?shù)，求得xd=，xq=，xd’=.為了簡(jiǎn)化網(wǎng)絡(luò)，將三臺(tái)發(fā)電機(jī)并聯(lián)等值成一臺(tái)機(jī)組。具體參數(shù)設(shè)置如圖3所示，機(jī)組的有功出力是不確定的，可在框內(nèi)暫時(shí)填為100MW。綜上，各節(jié)點(diǎn)的參數(shù)如表2所示：2 采用MATLAB計(jì)算電網(wǎng)潮流通過(guò)對(duì)待求解電網(wǎng)的等效線路參數(shù)的計(jì)算與標(biāo)幺化，得出待求解電網(wǎng)的各PQ節(jié)點(diǎn)的有功功率和無(wú)功功率，各PV節(jié)點(diǎn)的有功功率和節(jié)點(diǎn)電壓。負(fù)荷采用恒阻抗模型。電力系統(tǒng)穩(wěn)定分析和計(jì)算設(shè)計(jì)報(bào)告電力系統(tǒng)穩(wěn)定分析和計(jì)算設(shè)計(jì)報(bào)告專業(yè)：電氣工程及其自動(dòng)化班級(jí)： 09電氣（4）班學(xué)號(hào)：32號(hào)、34號(hào)姓名：禤培正、陳子程華南理工大學(xué)電力學(xué)院2012229目錄課程設(shè)計(jì)題目：電力系統(tǒng)穩(wěn)定分析和計(jì)算 11 潮流的參數(shù)計(jì)算和標(biāo)幺化 4 線路參數(shù)計(jì)算與標(biāo)幺化 4 節(jié)點(diǎn)信息統(tǒng)計(jì) 52 采用MATLAB計(jì)算電網(wǎng)潮流 6 PQ分解法潮流計(jì)算 6 MATLAB平臺(tái)潮流計(jì)算結(jié)果分析 73 采用powerworld仿真軟件計(jì)算潮流 7 PowerWorld仿真軟件簡(jiǎn)介 7 PowerWorld軟件建模與潮流計(jì)算 8 PowerWorld仿真結(jié)果與Matlab計(jì)算結(jié)果對(duì)比分析 114 網(wǎng)絡(luò)變換法求解輸入阻抗與轉(zhuǎn)移阻抗 12 負(fù)荷、發(fā)電機(jī)、變壓器、線路等效電容的處理 12 采用網(wǎng)絡(luò)變換法的具體步驟 15 求解輸入阻抗和轉(zhuǎn)移阻抗 19 求解有功功率傳輸特性 195 故障方案仿真 20 選定故障1：節(jié)點(diǎn)3發(fā)生三相對(duì)稱短路 20 選定故障2：線路L23中點(diǎn)處發(fā)生三相對(duì)稱短路 266 采用等面積法求解極限切除角 31 故障方案一 31 故障方案二 327 基于PowerWorld，分析三階模型的發(fā)電機(jī)的暫態(tài)過(guò)程 33 故障方案一：節(jié)點(diǎn)3發(fā)生三相對(duì)稱短路 33 故障方案二：線路L23中點(diǎn)處發(fā)生三相對(duì)稱短路 358 發(fā)電機(jī)模型及勵(lì)磁調(diào)節(jié)參數(shù)對(duì)穩(wěn)定計(jì)算結(jié)果的影響 36 發(fā)電機(jī)模型對(duì)穩(wěn)定計(jì)算結(jié)果的影響 36 TJ對(duì)穩(wěn)定計(jì)算結(jié)果的影響 38 勵(lì)磁參數(shù)對(duì)穩(wěn)定計(jì)算結(jié)果的影響 39參考文獻(xiàn) 40附錄： 40[鍵入文字]課程設(shè)計(jì)題目：電力系統(tǒng)穩(wěn)定分析和計(jì)算姓名：指導(dǎo)教師：一、 一個(gè)220kV分網(wǎng)結(jié)構(gòu)和參數(shù)如下：500kV站（1）的220kV母線視為無(wú)窮大母線，電壓恒定在230kV。二、 設(shè)計(jì)的主要內(nèi)容：進(jìn)行參數(shù)計(jì)算和標(biāo)幺化，形成潮流計(jì)算參數(shù)；用Matlab編制潮流計(jì)算程序，要求采用PQ分解法潮流計(jì)算方法?；贛ATLAB平臺(tái)編寫(xiě)程序，我們對(duì)待求解電網(wǎng)進(jìn)行潮流計(jì)算，從而得出各PQ節(jié)點(diǎn)的節(jié)點(diǎn)電壓以及各線路的線路損耗。圖接在平衡節(jié)點(diǎn)處的發(fā)電機(jī)參數(shù)設(shè)置節(jié)點(diǎn)6為PV節(jié)點(diǎn)，P=，V=，為了方便后面的穩(wěn)定計(jì)算，按題目中要求，節(jié)點(diǎn)6接3臺(tái)參數(shù)一樣的機(jī)組，其中正常運(yùn)行時(shí)各機(jī)組發(fā)出的有功功率P1=P2=P3=。等效后的機(jī)組出力P=，Xd=，Xq=，Xd’=發(fā)電機(jī)采用二階經(jīng)典模型，即E’恒定模型。假設(shè)在5s時(shí)刻節(jié)點(diǎn)3發(fā)生了三相對(duì)稱短路，則采用Matlab編制的程序如下所示：子函數(shù):functionxdot=simb(t,x,flag,a,b,c)xdot=[(x(2)1)*。c=*pi/180。其中的Matlab求解程序如下：（子程序與穩(wěn)定計(jì)算程序一樣）functionxdot=simb(t,x,flag,a,b,c)xdot=[(x(2)1)*。點(diǎn)擊圖中交點(diǎn)后，輸出x =，y =，即==*180/pi=。代入式13后得：（18）當(dāng)節(jié)點(diǎn)3的故障切除后，通過(guò)網(wǎng)絡(luò)變換法同理可得，轉(zhuǎn)移阻抗Z17= + =，輸入阻抗Z77= + =。b=。圖3tc=同理。從總體上看，三種方法求出的結(jié)果較為接近。所得搖擺曲線如下，此時(shí)系統(tǒng)暫態(tài)穩(wěn)定。當(dāng)TJ足夠大后，再增大TJ時(shí)，極限切除角變化幅度不大。0 +*i 1 0 0 2。 isb=input(39。請(qǐng)輸入PQ節(jié)點(diǎn)數(shù)na=39。YI(p,q)=YI(p,q)1./B1(i,3)。V(i)=B2(i,4)。for j1=k:naA(i,j1)=A(i,j1)./A(i,i)。 DET=abs(DP1(i))。endDQ1(i)=Q(i)V(i)*C(i)。ifi==nafor LZ=2:i L=i+2LZ。)。disp(39。各條支路的首端功率Sj（順序同您輸入B1時(shí)一樣）為：39。disp(Sj(q,p))。迭代次數(shù)39。)。q=B1(i,1)。disp(S)。)。endenddisp(39。if ICT3~=0 L=0。L=L+1。endDP1(i)=P(i)V(i)*C(i)。endendendendfori=1:na ifi==naA(i,i)=1./A(i,i)。fori=1:n e(i)=real(B2(i,3))。q=B1(i,1)。)。請(qǐng)輸入支路數(shù)：nl=39。B2=[0 0 1 0 1。在PowerWorld軟件里，發(fā)電機(jī)采用三階模型，勵(lì)磁系統(tǒng)采用BPA_EA模型，TJ分別取50時(shí)，分別求出系統(tǒng)的極限切除角，其中故障選取故障方案1：節(jié)點(diǎn)3發(fā)生三相對(duì)稱短路。所得搖擺曲線如下，此時(shí)系統(tǒng)暫態(tài)不穩(wěn)定。 故障方案二由第5部分知，L23中點(diǎn)處發(fā)生三相對(duì)稱短路，當(dāng)故障切除后。圖tc=，修改相應(yīng)的Matlab程序，得到的搖擺曲線如下，即此時(shí)發(fā)電暫態(tài)不穩(wěn)定。%求解第一個(gè)區(qū)間[0,tc]x02= x1 (end,:)39。所以，為了檢驗(yàn)系統(tǒng)的暫態(tài)穩(wěn)定特性，文本所提到的線路故障皆為雙回路在同一點(diǎn)發(fā)生短路，即線路L23的兩個(gè)回路同時(shí)在中點(diǎn)處發(fā)生三相對(duì)稱短路。)%畫(huà)出搖擺曲線line([,],[0,3]) %畫(huà)t=圖2故障時(shí)功角隨時(shí)間變化的特性已知極限切除時(shí)間tc=，則可由圖22的曲線上找到對(duì)應(yīng)的切除角，即極限切除角。圖1tc=(1) 采用試湊法不斷改變切除時(shí)間，即給tc賦值，找到穩(wěn)定與不穩(wěn)定的臨界時(shí)間，如圖21所示，取tc=。%給第二個(gè)計(jì)算區(qū)間賦初值a=。由于Pe的表達(dá)式在不同的狀態(tài)下（故障前、故障后、切除后）是不一樣的，所以需要Matlab定義多個(gè)時(shí)間段進(jìn)行仿真。在本例中，正常運(yùn)行下負(fù)荷所在節(jié)點(diǎn)電壓近似為VN，由式9求得各負(fù)荷的等效阻抗如表6所示：表1各負(fù)荷的等效阻抗負(fù)荷所在節(jié)點(diǎn)節(jié)點(diǎn)電壓V（kV）負(fù)荷大小S（MVA）等效阻抗ZL的標(biāo)幺值2220250+j40+3220380+j90+4220280+j85+5220420+j95+ 發(fā)電機(jī)的處理接在節(jié)點(diǎn)6的有三臺(tái)發(fā)電機(jī)，發(fā)電機(jī)的容量SG(N)=300/=，VG(N)=，而基準(zhǔn)容量SB=100MW，為了使變壓器的標(biāo)幺變比k*=1，在10kV電壓等級(jí)下的基準(zhǔn)電壓VB2=*220/242kV=。 PowerWorld軟件建模與潮流計(jì)算 采用PowerWorld仿真軟件建立潮流模型在軟件的編輯模式下