【正文】 %P186公式741狀態(tài)方程各個參數(shù) K1=Eqq*U*cos(delta*pi/180)/Xdd_all+U*U*(Xdd_allXd_all)*cos(2*delta*pi/180)/(Xd_all*Xdd_all)。 UGd=U*Xd*sin(delta*pi/180)/Xd_all。 b=2*(Ke^2*Xd*Xe*U*cos(delta*pi/180)/(Xd_all^2)+Eq0+Ke*UG0)。 %系統(tǒng)電抗(包含暫態(tài)電抗)Tj=10。 %線路電抗值Xd=。 U=1。發(fā)電機2339。%由序分量合成相分量 考慮Yd11聯(lián)結組別Iabc_effective=abs(Iabc) %支路25%mn=39。I1_branch=(UU1(1)UU1(4))/z1_branch(4)I2_branch=(UU2(1)UU2(4))/z2_branch(4)I0_branch=(U0(1)0)/(j*)Iabc=T*[I1_branch I2_branch I0_branch].39。strUU1=U2_Node1UU2=U2_Node2m=4U4abc=T*[U2_Node1(4,1)*exp(j*pi/6) U2_Node2(4,1)*exp(j*pi/6) 0].39。U2_Node2=U2fault_node2+U2normal_node2。U2_Node1=zeros(bus_Num2,1)。U2fault_node2=YY2\I2fault_node2。 endend%故障分量U2fault_node1=zeros(bus_Num2,1)。 I2fault_node0(m,1)=1*I_Fault0。I2fault_node2=zeros(bus_Num2,1)。strI1= I_branch1I2= I_branch2I0=I_branch0str=39。 I_branch2(i,1)=(U1fault_node2(m,1)+U1normal_node2(m,1)U1fault_node2(n,1)U1normal_node2(n,1))/z2_branch(i)。I_branch0=zeros(bus_Num1,1)。39。U1_node1=U1fault_node1+U1normal_node1。U1normal_node0=zeros(bus_Num1,1)。U1fault_node2=Y2\I1fault_node2。 endend%故障分量U1fault_node1=zeros(bus_Num1,1)。 I1fault_node0(m,1)=1*I_Fault0。I1fault_node2=zeros(bus_Num1,1)。strI1=I_Fault1I2=I_Fault2I0=I_Fault0str=39。 I_Fault0=I_Fault1*Z2(Fault_Node,Fault_Node)/((Z2(Fault_Node,Fault_Node)+Z0(Fault_Node,Fault_Node)))。 I_Fault2=I_Fault1。 elseif Fault_Type==2 %a相經(jīng)10歐姆電阻接地短路各序電流計算 I_Fault1=1/(Z1(Fault_Node,Fault_Node)+Z2(Fault_Node,Fault_Node)+Z0(Fault_Node,Fault_Node)+3*10*60/(115^2))。 I_Fault0=0。) a=+j*sqrt(3)/2。輸入短路類型 。輸入短路點編號 。Z0=inv(Y0)。 YY0(c,b)=a。 b=branch2_No1(i)。YY0=zeros(bus_Num2)。 YY1(c,c)=YY1(c,c)+a。 c=branch2_No2(i)。YY1(4,4)=1/(j*)。end %第二步：考慮發(fā)電機節(jié)點計算節(jié)點導納矩陣%%節(jié)點導納矩陣，YY1表示計及發(fā)電機節(jié)點的正序網(wǎng)絡節(jié)點導納陣，YY2表示計及發(fā)電機節(jié)點的負序網(wǎng)絡節(jié)點導納陣z1_branch(4)=j*。 Y0(b,b)=Y0(b,b)+a。Y0(2,2)=1/(j*)。 Y1(c,c)=Y1(c,c)+a。 c=branch1_No2(i)。Y1(1,1)=1/(j*)。 z0_branch(1)=j*。 %輸入支路各序阻抗，z1_branch表示支路正序阻抗，z2_branch表示支路負序阻抗，z0_branch表示支路零序阻抗z1_branch(1)=j*。 %包括發(fā)電機支路的支路數(shù)branch1_No1=[1,1,2]。 %輸入已知條件%bus_Num1=3。 減小元件的電抗252。，暫態(tài)穩(wěn)定的概念？答：靜態(tài)穩(wěn)定是指電力系統(tǒng)受到小干擾之后，不發(fā)生非周期性失步或自發(fā)震蕩，自動恢復到初始運行狀態(tài)的能力。 潮流計算時需要迭代計算出結果，而短路計算一次就可以計算出短路電流。而用計算機時可以不用簡化網(wǎng)絡,求出各序導納矩陣后利用I=YU來計算各節(jié)點正負零各序電流后再轉換成三相短路電流。 Ia= Ib= 0 Ic= 0b. 各節(jié)點三序電壓、三相電壓：節(jié)點U1U2U0UaUbUc1 + 2 + 3 + 4 0 + + 5 0 + + c. 各支路三序電流、三相電流：支路I1I2I0IaIbIc12 + + + + 13 23 + + 14 + + + + + + 25 + + + + 發(fā)電機1 + + 0 + 發(fā)電機23 + + 0 + (3) 節(jié)點3發(fā)生b、c兩相短路故障a. 故障點三序電流、三相電流：I1= 0 I2= 0 + I0= 0 。四、參考文獻1. 《電力系統(tǒng)暫態(tài)分析》，李光琦，中國電力出版社，2006年，第三版；2. 《電力系統(tǒng)分析》（上、下），何仰贊，華中科技大學出版社，1998年，第二版。通過此次學習，了解到電力系統(tǒng)的故障時，大部分電磁量將隨時間變化，描述其特性的是微分方程，這給分析計算帶來一定困難。在短路電流的計算過程中，先通過計算系統(tǒng)各個節(jié)點的導納矩陣，再通過判斷系統(tǒng)的故障節(jié)點號和故障類型代入相應的公式進行計算。二、設計正文（詳細內容見附錄,用A4紙，頁數(shù)不限）1． 對稱短路計算過程流程圖和計算結果2． 不對稱短路計算過程流程圖和計算結果3． 靜態(tài)穩(wěn)定性計算過程流程圖和計算結果4． 暫態(tài)穩(wěn)定性計算過程流程圖和計算結果 思考題 實驗程序三、課程設計總結或結論本次電力系統(tǒng)暫態(tài)上機主要討論了計算各種類型短路故障下系統(tǒng)網(wǎng)絡中的電壓電流分布以及電力系統(tǒng)的靜態(tài)穩(wěn)定、暫態(tài)穩(wěn)定問題。在暫態(tài)穩(wěn)定性的計算中，設置了相應的切除時間和分段時間，每一個分段時間進行一次計算，在切除故障前和后分別進行計算，進而得出相應曲線，并通過多次嘗試的除最大切除時間大約為2s。另外，通過此次課程設計，我學習并掌握了一定的Matlab知識，對簡單的選擇語句、循環(huán)語句、數(shù)組的使用等有了一些了解，為今后學習、工作中再次學習使用此軟件打下了基礎，獲益匪淺。 Ia= 0 Ib= 0 Ic= 0 b. 各節(jié)點三序電壓、三相電壓：節(jié)點U1U2U0UaUbUc1 + 2 + 3 +