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