【正文】 ak(13)=(Ke*(VrVt+x(12)+T1/T2*(x(11)pop*T1/T2*Ur))x(13)+)/Te。ak(11)=(pop*(1T1/T2)*Urx(11))/T2。ak(9)=(x(8)+M0x(9))/Tch。ak(7)=((1T5/T6)*(x(6)+T5/T6*(x(5)+T5/T6*Ug))x(7))/T6。ak(5)=((1T5/T6)*Ugx(5))/T6。ak(3)=(x(10)+a1*x(9)D*x(3)(x(1)*Id+x(2)*Iq(Xq1Xd1)*Id*Iq))/Tj。endak(1)=(x(1)(XqXq1)*Iq)/Tq。end if x(8)+M0C2 x(8)=C2M0。 if x(13)C4 x(13)=C4。v2=Vd*Vd+Vq*Vq。end Vd=x(1)r*IdXd1*Iq。 else Id=((r+Re)*(x(1)+V*sin(x(4)))(Xd1+Xe)*(x(2)V*cos(x(4))))/((r+Re)*(r+Re)+(Xd1+Xe)*(Xd1+Xe))。t= Id=((r+Re1)*(x(1)+V*sin(x(4)))(Xd1+Xe1)*(x(2)V*cos(x(4))))/((r+Re1)*(r+Re1)+(Xd1+Xe1)*(Xd1+Xe1))。elseif tamp。 t= Id=((1+a*rb*Xd1)*(a*x(1)+b*x(2)+d*V*cos(x(4))c*V*sin(x(4)))(a*Xd1+b*r)*(b*x(1)+a*x(2)+c*V*cos(x(4))+d*V*sin(x(4))))/((1+a*rb*Xd1)*(1+a*rb*Xd1)+(a*Xd1+b*r)*(a*Xd1+b*r))。 if t0 amp。C2=。Ur=Wrx(3)。Wr=1。Vr=1。d=。b=。d=。b=。T6=。T2=。M0=。Trh=8。Ts=。Xe1=。Xe=。D=0。Tq=。r=。Xd1=。Xd=。Te=。)fit=[z(4,:)]*[(z(4,:))]39。,39。,39。)%save(39。,39。,39。 %save(39。 %dA1=Angles(1,:)39。%Angles(1,:)=z(1,:)。 %plot(Angles39。 end %Angles=zeros(1,fn)。 z(:,fi)=x。 ak=myfunfun(x,nx,t,pop)。 for fj=1:3 x=z(:,fi1)+ak.*ad(fj)。 for fi=2:fn ak=myfunfun(x,nx,t,pop)。 h。ad=[h/2。 z(:,1)=x。]。0。0。0。x=[。 fn=18000。 ff=50。說明：同樣RLPSS和RLPSSfun在同一個文件夾中。ak(12)=((1T1/T2)*(x(11)K0*T1/T2*Ur)x(12))/T2。ak(10)=((1a1)*x(9)x(10))/Trh。ak(8)=(Kg*(x(7)+T5/T6*(x(6)+T5/T6*(x(5)+T5/T6*Ug)))x(8))/Ts。ak(6)=((1T5/T6)*(x(5)+T5/T6*Ug)x(6))/T6。ak(4)=(x(3)1)*Wb。ak(2)=(x(13)x(2)+(XdXd1)*Id)/Td。elseif x(8)+M0C1 x(8)=(C1M0)。elseif x(13)C4 x(13)=C4。Vt=sqrt(v2)。Vq=x(2)r*Iq+Xd1*Id。 Iq=(x(1)+V*sin(x(4))(r+Re)*Id)/(Xd1+Xe)。 Iq=(x(1)+V*sin(x(4))(r+Re1)*Id)/(Xd1+Xe1)。amp。 Iq=(a*x(1)+b*x(2)+d*V*cos(x(4))c*V*sin(x(4))(1+a*rb*Xd1)*Id)/(a*Xd1+b*r)。amp。C1=。C4=。Ug=Wrx(3)。Wb=。V=。c=。%a=。c=。a=。T5=。T1=。M0=。Trh=8。Ts=。Xe1=。Xe=。D=0。Tq=。r=。Xd1=。Xd=。Te=。ASCII39。dA239。d:\39。ASCII39。dA139。d:\39。%dA2=Angles(2,:)39。%Angles(2,:)=z(2,:)。) %Angles=zeros(2,fn)。Angles(1,:)=z(4,:)。 t=t+h。 end x=b+h*ak/6。 b=b+ak.*ad(fj+1)/3。 b=x。 h]。 h/2。 ak=zeros(nx,1)。 z=zeros(nx,fn)。0。0。0。1。h=1/ff/80。 nx=13。3. RLPSSclear allclose all%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%t=0。ak(8)=(Ke*(VrVt+T1/T2*(K0*T1/T2*Ur))x(8)+)/Te。ak(6)=(x(5)+M0x(6))/Tch。ak(4)=(x(3)1)*Wb。ak(2)=(x(8)x(2)+(XdXd1)*Id)/Td。elseif x(8)+M0C1 x(8)=(C1M0)。elseif x(13)C4 x(13)=C4。Vt=sqrt(v2)。Vq=x(2)r*Iq+Xd1*Id。 Iq=(x(1)+V*sin(x(4))(r+Re)*Id)/(Xd1+Xe)。 Iq=(x(1)+V*sin(x(4))(r+Re1)*Id)/(Xd1+Xe1)。amp。 Iq=(a*x(1)+b*x(2)+d*V*cos(x(4))c*V*sin(x(4))(1+a*rb*Xd1)*Id)/(a*Xd1+b*r)。amp。C1=。C4=。Ug=Wrx(3)。Wb=。V=。c=。%a=。c=。a=。T5=。T1=。M0=。Trh=8。Ts=。Xe1=。Xe=。D=0。Tq=。r=。Xd1=。Xd=。Te=。ASCII39。dA239。d:\39。ASCII39。dA139。d:\39。%dA2=Angles(2,:)39。%Angles(2,:)=z(2,:)。) %Angles=zeros(2,fn)。Angles(1,:)=z(4,:)。 t=t+h。 end x=b+h*ak/6。 b=b+ak.*ad(fj+1)/3。 b=x。 h]。 h/2。 ak=zeros(nx,1)。 z=zeros(nx,fn)。0。0。0。1。h=1/ff/80。 nx=13。附 錄1. RL的程序：clear allclose all%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%t=0。 展望未來隨著世界經濟的發(fā)展，人們對供電的可靠性和穩(wěn)定性日益提高。所以我們通過對無限大系統(tǒng)母線在附加PSS和不加PSS兩種狀態(tài)下系統(tǒng)運行情況的比較，得出了PSS在抑制系統(tǒng)低頻振蕩方面具有明顯的效果，通過兩幅圖的比較也可以看出在附加PSS的系統(tǒng)振蕩會很快的趨于平緩，而未加PSS的系統(tǒng)則一直振蕩。 求矩陣的特征值將的值帶入下面的矩陣得到上述矩陣的全部特征值： 畫出上述特征值的曲線：圖52 特征值根曲線 在含有PSS的情況下系統(tǒng)的狀態(tài)狀態(tài)方程的詳細表示式（非線性表達式）為：根據(jù)以上的系統(tǒng)模型公式可以運用MATLAB編程，得到程序附錄3，所以得到以下的功角曲線：圖53 功角特性曲線圖54 變化曲線 主要結論和展望 主要結論本文主要分析了電力系統(tǒng)的特點，說明了現(xiàn)代的電力系統(tǒng)存在的種種問題，提出提高電力系統(tǒng)穩(wěn)定性的重要性。 幾個參數(shù)的確定 PSS網(wǎng)絡的計算同步功率系數(shù)固有振蕩頻率當時， 汽輪機及其調節(jié)系統(tǒng)的超前補償網(wǎng)絡的計算通路的相位滯后角度為取，那么這里的確定我們是運用的遺傳算法，我們知道遺傳算法具有取最優(yōu)解的特點，本課題我們把（范圍是1）帶入到myfun（附錄5）程序中，然后通過計算求出，然后比較所求出的，我們的要求是求出的最小值，然后GA工具箱會返回給我們一個最小值，這個時候就是我們要求的最優(yōu)值。發(fā)現(xiàn)在那種情況下，系統(tǒng)的性能很差，不能很穩(wěn)定的運行。( 8) 重復步驟( 3) 至( 7) , 直至參數(shù)收斂達到預定的指標。( 6) 隨機交叉, 兩個個體按一定的交叉概率進行交叉操作, 并產生兩個新的子個體。( 4) 對收斂條件進行判斷, 如果已經找到最佳個體, 則停止,否則繼續(xù)進行遺傳操作。( 2) 隨機產生包含各個個體的初始種群。變異算子是對個體的某一個或幾位按某一較小的概率進行反轉其二進制字符, 模擬自然界的基因突變現(xiàn)象。選擇算子的作用是用來提高群體的平均適應度值。遺傳操作包括三個算子: 選擇、重組和變異。在遺傳算法中, 首先將空間問題中的決策變量通過一定的編碼表示成遺傳空間的一個個體, 它是一個基因型串結構數(shù)據(jù)。除內部函數(shù)外，所有MATLAB主包文件和各種工具包都是可讀可修改的文件，用戶通過對源程序的修改或加入自己編寫程序構造新的專用工具包。學科工具包是專業(yè)性比較強的工具包，控制工具包，信號處理工具包，通信工具包等都屬于此類。工具箱又可以分為功能性工具箱和學科工具箱。（2） 在MATLAB的命令窗口中輸入“quit”命令（3） 直接單擊MATLAB命令窗口的按鈕。 MATLAB的啟動與退出（1）MATLAB的啟動，在windows操作系統(tǒng)中，選擇“開始”“程序”“MATLAB” “R2007b”“MATLAB”，就可以啟動MATLAB的工作界面了。工具箱是MATLAB函數(shù)的子程序庫，每一個工具箱都是為某一類學科專業(yè)和應用而定制的，主要包括信號處理、控制系統(tǒng)、神經網(wǎng)絡、模糊邏輯、小波分析和系統(tǒng)仿真等方面的應用。另外，MATLAB網(wǎng)頁服務程序還容許在Web應用中使用自己的MATLAB數(shù)學和圖形程序。（6）實用的程序接口和發(fā)布平臺新版本的MATLAB可以利用MATLAB編譯器和C/C++數(shù)學庫和圖形庫，將自己的MATLAB程序自動轉換為獨立于MATLAB運行的C和C++代碼。一般來說，它們都是由特定領域的專家開發(fā)的，用戶可以直接使用工具箱學習、應用和評估不同的方法而不需要自己編寫代碼。另外新版本的MATLAB還著重在圖形用戶界面（GUI）的制作上作了很大的改善，對這方面有特殊要求的用戶也可以得到滿足。新版本的MATLAB對整個圖形處理功能作了很大的改進和完善，使它不僅在一般數(shù)據(jù)可視化軟件都具有的功能（例如二維曲線和三維曲面的繪制和處理等）方面更加完善，而且對于一些其他軟件所沒有的功能（例如圖形的光照處理、色度處理以及四維數(shù)據(jù)的表現(xiàn)等），MATLAB同樣表現(xiàn)了出色的處理能力。高層次的作圖包括二維和三維的可視化、圖象處理、動畫和表達式作圖。函數(shù)所能解決的問題其大致包括矩陣運算和線性方程組的求解、微分方程及偏微分方程的組的求解、符號運算、傅立葉變換和數(shù)據(jù)的統(tǒng)計