電力系統(tǒng)課程設(shè)計(jì)(附分段法程序)-資料下載頁(yè)
2025-07-07 15:29本頁(yè)面
【正文】 ：機(jī)械工業(yè)出版社，2011.[2]陳珩 電力系統(tǒng)穩(wěn)態(tài)分析（第三版）[M]. 北京：中國(guó)電力出版社，2007.[3]李光琦 電力系統(tǒng)暫態(tài)分析（第三版）[M]. 北京：中國(guó)電力出版社，2007.[4]于永源 楊綺雯 電力系統(tǒng)分析（第三版）[M].北京：中國(guó)電力出版社，2007附錄部分源程序：S0=1+。 %輸入初始功率 V0=1。 %輸入無(wú)限大系統(tǒng)母線電壓 B=[ ]。 %輸入系統(tǒng)等值電抗矩陣 Tj=。 %輸入慣性時(shí)間常數(shù) N=300。 %輸入故障次數(shù) Ni=29。 %輸入哪個(gè)時(shí)段發(fā)生故障 Nj=290。 %輸入哪個(gè)時(shí)段故障消除dt=。 %輸入每段間隔的時(shí)間for i=1:Ni %故障前 C(i,1)=1。C(i,2)=B(1)。endfor i=Ni:Nj %故障時(shí) C(i,1)=0。C(i,2)=B(2)。endfor i=Nj:N %故障后 C(i,1)=2。C(i,2)=B(3)。endE0=sqrt((V0+imag(S0)*B(1)./V0)^2+(real(S0)*B(1)./V0)^2)。 %求解初始電壓的過(guò)程 dtj0=atan(real(S0)*B(1)./(V0*(V0+imag(S0)*B(1)./V0)))。 %求解功角 Q(1)=dtj0*180./pi。 %將弧度轉(zhuǎn)變?yōu)榻嵌? K1=360*50*dt^2./Tj。 %求解靜態(tài)穩(wěn)定儲(chǔ)備系數(shù)P1m=E0*V0./B(1)。 %故障前功率特性曲線峰值 P2m=E0*V0./B(2)。 %故障時(shí)功率特性曲線峰值 P3m=E0*V0./B(3)。 %故障后功率特性曲線峰值 dtk=piasin(real(S0)./P3m)。dtjx=acos((real(S0)*(dtkdtj0)+P3m*cos(dtk)P2m*cos(dtj0))./(P3mP2m))。%求解極限切除角dtjx=dtjx*180./pi。 %將弧度轉(zhuǎn)變?yōu)榻嵌? for K=2:N+1。 P(K1)=E0*V0./C(K1,2)。 if C(K1,1)==1 %故障前 Peq(K1)=P(K1)*sin(Q(K1)*pi./180)。 %求解輸出電磁功率 DP(K1)=real(S0)Peq(K1)。 %求解不平衡功率 if K==2 %第一時(shí)段 DQ(K)=K1*DP(K1)./2。 else DQ(K)=DQ(K1)+K1*DP(K1)。 %求解功角增量的過(guò)程 end Q(K)=Q(K1)+DQ(K)。 %求解輸出功角 elseif C(K1,1)==0 %故障時(shí) P1(K1)=P(K2)。 P2(K1)=P(K1)。 Peq(K1)=P(K1)*sin(Q(K1)*pi./180)。 %求解輸出電磁功率 DP(K1)=real(S0)Peq(K1)。 %求解不平衡功率 if K==Ni+1 %第一時(shí)段 DQ(K)=DQ(K1)+K1*DP(K1)./2。 else DQ(K)=DQ(K1)+K1*DP(K1)。 %求解功角增量 end Q(K)=Q(K1)+DQ(K)。 %求解輸出功角 elseif C(K1,1)==2 %故障切除 P1(K1)=P(K2)。 P2(K1)=P(K1)。 Peq(K1)=P1(K1)*sin(Q(K1)*pi./180)。 %求輸出電磁功率 DP(K1)=real(S0)Peq(K1)。 %求解不平衡功率 DP2(K1)=real(S0)P2(K1)*sin(Q(K1))。if K==Nj+1 %第一時(shí)段 DQ(K)=DQ(K1)+K1*DP2(K1)./2。 else DQ(K)=DQ(K1)+K1*DP2(K1)。 %求解功角增量 end Q(K)=Q(K1)+DQ(K)。 %求解輸出功角 endenddisp(39。程序計(jì)算結(jié)果如下：39。)。disp(39。輸出不平等功率:DP=39。)。disp(DP)。disp(39。輸出功角增量:DQ=39。)。disp(DQ)。disp(39。輸出功角:Q=39。)。disp(Q)。disp(39。輸出電磁功率:Peq39。)。 disp(Peq)。disp(39。輸出發(fā)電機(jī)轉(zhuǎn)子搖擺曲線如圖Figure :39。)。t=0::dt*N。plot(t,dtjx)。hold on t=0:dt:dt*N。plot(t,Q,t,dtjx)。xlabel(39。時(shí)間(秒)39。)。ylabel(39。功角(度數(shù))39。)。title(39。發(fā)電機(jī)搖擺曲線39。)
外語(yǔ)相關(guān)推薦
鄂ICP備17016276號(hào)-1