電力系統(tǒng)穩(wěn)定器pss設(shè)計與仿真本科畢業(yè)論文

電力系統(tǒng)穩(wěn)定器PSS的設(shè)計與仿真摘要由于電力系統(tǒng)在正常運(yùn)行時會發(fā)生頻率的振蕩，對我們的生產(chǎn)生活帶來了很多的危害，給我國的國民經(jīng)濟(jì)造成了巨大的損失。如果在電力系統(tǒng)中加入PSS后會對系統(tǒng)的穩(wěn)定性提高給予了很大的幫助。電力系統(tǒng)穩(wěn)定器就是為抑制低頻振蕩而研究的一種附加勵磁控制技術(shù)。它在勵磁電壓調(diào)節(jié)器中，引入領(lǐng)先于軸速度的附加信號，產(chǎn)生一個正阻尼轉(zhuǎn)矩，去克服原勵磁電壓調(diào)節(jié)器中產(chǎn)生的負(fù)阻尼轉(zhuǎn)矩作用,用于提高電力系統(tǒng)阻尼，解決低頻振蕩問題，是提高電力系統(tǒng)動態(tài)穩(wěn)定性的重要措施之一。本文首先分析了電力系統(tǒng)穩(wěn)定器的研究背景和國內(nèi)外研究狀況，然后又對無限大系統(tǒng)的數(shù)學(xué)模型進(jìn)行了具體的分析，然后又分析了系統(tǒng)帶PSS和不帶PSS時系統(tǒng)的運(yùn)行情況，最后對兩種情況進(jìn)行了比較，分析了PSS的優(yōu)點(diǎn)和缺點(diǎn)，同時對電力系統(tǒng)的未來提出了更高的挑戰(zhàn)。關(guān)鍵詞：勵磁控制技術(shù);低頻振蕩；電力系統(tǒng)穩(wěn)定器DesignandSimulationofPowerSystemStabilizer(PSS)ABSTRACTThispaperdescribesthepowersystemstabilizer'sdesignprinciplesandtheimpactonthepowersystem,theuseofthebenefitsofpowersystemstabilizer.ThefrequencyofoscillationmaybehappentopowersystemwhenthesystemInthenormaloperation.ItgivesagreathelpifthepowersystemaftertheadditionofPSS.Powersystemstabilizer(PSS)istosuppressalowfrequencyoscillationofadditionalexcitationcontrol.Itistheexcitationvoltageregulator,theintroductionofaxialvelocityaheadofadditionalsignalstoproduceapositivedampingtorquetoovercometheprimaryexcitationvoltageregulatorproducednegativedampingtorqueeffect.Improvingpowersystemdamping,lowingfrequencyoscillationproblemsolvingistoimprovepowersystemdynamicstabilityoftheimportantmeasures.Thispaperanalyzesthepowersystemstabilizerresearchbackgroundandresearchstatusathomeandabroad,andthentotheinfinitesystemofmathematicalmodelsofspecificanalysis,andthenanalyzedthesystemwithPSSandwithoutPSSoperationofthesystem,andfinallyOfthetwocaseswerecomparedandanalyzedtheadvantagesanddisadvantagesofPSS,whilethefutureofthepowersystemofahigherchallenge.Keyword:Excitationcontrol；Low-frequencyoscillation；PowerSystemStabilizer畢業(yè)設(shè)計（論文）原創(chuàng)性聲明和使用授權(quán)說明原創(chuàng)性聲明本人鄭重承諾：所呈交的畢業(yè)設(shè)計（論文），是我個人在指導(dǎo)教師的指導(dǎo)下進(jìn)行的研究工作及取得的成果。盡我所知，除文中特別加以標(biāo)注和致謝的地方外，不包含其他人或組織已經(jīng)發(fā)表或公布過的研究成果，也不包含我為獲得及其它教育機(jī)構(gòu)的學(xué)位或?qū)W歷而使用過的材料。對本研究提供過幫助和做出過貢獻(xiàn)的個人或集體，均已在文中作了明確的說明并表示了謝意。作者簽名：日期：指導(dǎo)教師簽名：日期：使用授權(quán)說明本人完全了解大學(xué)關(guān)于收集、保存、使用畢業(yè)設(shè)計（論文）的規(guī)定，即：按照學(xué)校要求提交畢業(yè)設(shè)計（論文）的印刷本和電子版本；學(xué)校有權(quán)保存畢業(yè)設(shè)計（論文）的印刷本和電子版，并提供目錄檢索與閱覽服務(wù)；學(xué)校可以采用影印、縮印、數(shù)字化或其它復(fù)制手段保存論文；在不以贏利為目的前提下，學(xué)?？梢怨颊撐牡牟糠只蛉績?nèi)容。作者簽名：日期：

學(xué)位論文原創(chuàng)性聲明本人鄭重聲明：所呈交的論文是本人在導(dǎo)師的指導(dǎo)下獨(dú)立進(jìn)行研究所取得的研究成果。除了文中特別加以標(biāo)注引用的內(nèi)容外，本論文不包含任何其他個人或集體已經(jīng)發(fā)表或撰寫的成果作品。對本文的研究做出重要貢獻(xiàn)的個人和集體，均已在文中以明確方式標(biāo)明。本人完全意識到本聲明的法律后果由本人承擔(dān)。作者簽名： 日期：年月日學(xué)位論文版權(quán)使用授權(quán)書本學(xué)位論文作者完全了解學(xué)校有關(guān)保留、使用學(xué)位論文的規(guī)定，同意學(xué)校保留并向國家有關(guān)部門或機(jī)構(gòu)送交論文的復(fù)印件和電子版，允許論文被查閱和借閱。本人授權(quán)大學(xué)可以將本學(xué)位論文的全部或部分內(nèi)容編入有關(guān)數(shù)據(jù)庫進(jìn)行檢索，可以采用影印、縮印或掃描等復(fù)制手段保存和匯編本學(xué)位論文。涉密論文按學(xué)校規(guī)定處理。作者簽名： 日期：年月日導(dǎo)師簽名：日期：年月日

注意事項1.設(shè)計（論文）的內(nèi)容包括：1）封面（按教務(wù)處制定的標(biāo)準(zhǔn)封面格式制作）2）原創(chuàng)性聲明3）中文摘要（300字左右）、關(guān)鍵詞4）外文摘要、關(guān)鍵詞5）目次頁（附件不統(tǒng)一編入）6）論文主體部分：引言（或緒論）、正文、結(jié)論7）參考文獻(xiàn)8）致謝9）附錄（對論文支持必要時）2.論文字?jǐn)?shù)要求：理工類設(shè)計（論文）正文字?jǐn)?shù)不少于1萬字（不包括圖紙、程序清單等），文科類論文正文字?jǐn)?shù)不少于1.2萬字。3.附件包括：任務(wù)書、開題報告、外文譯文、譯文原文（復(fù)印件）。4.文字、圖表要求：1）文字通順，語言流暢，書寫字跡工整，打印字體及大小符合要求，無錯別字，不準(zhǔn)請他人代寫2）工程設(shè)計類題目的圖紙，要求部分用尺規(guī)繪制，部分用計算機(jī)繪制，所有圖紙應(yīng)符合國家技術(shù)標(biāo)準(zhǔn)規(guī)范。圖表整潔，布局合理，文字注釋必須使用工程字書寫，不準(zhǔn)用徒手畫3）畢業(yè)論文須用A4單面打印，論文50頁以上的雙面打印4）圖表應(yīng)繪制于無格子的頁面上5）軟件工程類課題應(yīng)有程序清單，并提供電子文檔5.裝訂順序1）設(shè)計（論文）2）附件：按照任務(wù)書、開題報告、外文譯文、譯文原文（復(fù)印件）次序裝訂3）其它

目錄摘要 )根據(jù)以上的系統(tǒng)模型公式可以運(yùn)用MATLAB編程，得到程序附錄3，所以得到以下的功角曲線：圖5-3功角特性曲線圖5-4變化曲線5.4主要結(jié)論本章主要分析了電力系統(tǒng)的特點(diǎn)，說明了現(xiàn)代的電力系統(tǒng)存在的種種問題，提出提高電力系統(tǒng)穩(wěn)定性的重要性。通過分析了國內(nèi)外多種提高電力系統(tǒng)穩(wěn)定性的方法，總結(jié)出PSS的應(yīng)用范圍比較廣泛，而且取得了不錯的效果。所以我們通過對無限大系統(tǒng)母線在附加PSS和不加PSS兩種狀態(tài)下系統(tǒng)運(yùn)行情況的比較，得出了PSS在抑制系統(tǒng)低頻振蕩方面具有明顯的效果，通過兩幅圖的比較也可以看出在附加PSS的系統(tǒng)振蕩會很快的趨于平緩，而未加PSS的系統(tǒng)則一直振蕩。通過以上的分析我們又總結(jié)出了PSS的一些不足點(diǎn)：首先，PSS的參數(shù)的整定是在某一個低頻振蕩的條件下計算的，如果頻率振蕩更加嚴(yán)重，在低頻振蕩的環(huán)境下計算出來的數(shù)據(jù)有可能在其他的情況下不是很適用；還有PSS的安裝地點(diǎn)應(yīng)該也是一個比較重要的問題；另外PSS的種類繁多，對一個機(jī)組裝PSS時有必要做測試，這樣也浪費(fèi)了大量的時間，從經(jīng)濟(jì)方面考慮又存在這問題。隨著世界經(jīng)濟(jì)的發(fā)展，人們對供電的可靠性和穩(wěn)定性日益提高。同時電力網(wǎng)絡(luò)的開發(fā)和商業(yè)化的運(yùn)營使得電力系統(tǒng)的運(yùn)行越來越趨于極限，電力系統(tǒng)不穩(wěn)定的因素越來越多，特別是電壓的不穩(wěn)和低頻振蕩問題更加嚴(yán)重，所以提出更有效的措施來解決問題可以說是當(dāng)務(wù)之急。致謝經(jīng)過半年的忙碌和工作，本次畢業(yè)設(shè)計已經(jīng)接近尾聲，作為一個本科生的畢業(yè)設(shè)計，由于經(jīng)驗(yàn)的匱乏，難免有許多考慮不周全的地方，如果沒有導(dǎo)師的督促指導(dǎo)，以及一起工作的同學(xué)們的支持，想要完成這個設(shè)計是難以想象的。在這里首先要感謝我的老師向老師。向老師平日里工作繁多，但在我做畢業(yè)設(shè)計的每個階段，從查閱資料到設(shè)計草案的確定和修改，中期檢查，后期詳細(xì)設(shè)計等整個過程中都給予了我悉心的指導(dǎo)。我的設(shè)計較為復(fù)雜煩瑣，但是向老師仍然細(xì)心地糾正模型中的錯誤。除了敬佩向老師的專業(yè)水平外，她的治學(xué)嚴(yán)謹(jǐn)和科學(xué)研究的精神也是我永遠(yuǎn)學(xué)習(xí)的榜樣，并將積極影響我今后的學(xué)習(xí)和工作。其次要感謝我的同組同學(xué)對我無私的幫助，特別是在MATLAB軟件的使用方面，正因?yàn)槿绱宋也拍茼樌耐瓿稍O(shè)計，我要感謝我的母?！骼砉ご髮W(xué)應(yīng)用科學(xué)學(xué)院，是母校給我們提供了優(yōu)良的學(xué)習(xí)環(huán)境；另外，我還要感謝那些曾給我授過課的每一位老師，是你們教會我專業(yè)知識。在此，我再說一次謝謝！謝謝大家?。?！。

參考文獻(xiàn)[1]韓禎祥．電力系統(tǒng)穩(wěn)定[M]．北京．中國電力出版社．1995，56-62．[2]倪以信．動態(tài)電力系統(tǒng)理論與分析[M]．北京．清華大學(xué)出版社．2002，3-6，15-16．[3]孫伯明．高等電力網(wǎng)絡(luò)分析[M]．北京．清華大學(xué)出版社．1996．[4]陳珩．電力系統(tǒng)穩(wěn)態(tài)分析[M]．北京．中國電力出版社．1995．[5]李光琦．電力系統(tǒng)暫態(tài)分析[M]．北京．中國電力出版社．1995．[6]王錫凡．現(xiàn)代電力系統(tǒng)分析[M]．北京．現(xiàn)代出版社．2003，333-342．[7]劉曉鵬．小干擾穩(wěn)定性部分特征值分析的多重變換法[J]．電力系統(tǒng)自動化．1998，22（9）：38-42．[8]蘇金明．MATLAB7.0實(shí)用指南[M]．北京．電子工業(yè)出版社．2004．[9]彭淦榕．MATLAB在電力系統(tǒng)分析中的應(yīng)用[J]．山東電力技術(shù)．2005．2（2）：49-53.[10]S.SUGIMOTO,ThyistorFanltCurrentLimitingSystem.IEEETransVol.11.No.2April1996[11]MAJinl&WANGHaoJing1&ZHANGPu2.RenewedinvestigationonPowerSystemStabilizerdesign.SchoolofElectricalandElectronicEngineering,NorthReceivedApril21,2011;acceptedMay23,2011;publishedonlineAugust7,2011.附錄A外文文獻(xiàn)原文RenewedinvestigationonPowerSystemStabilizer"ScienceandTechnology"inOctober2011Powersystemstabilizer1960swasintroducedtosolvetheproblemoflow-frequencyoscillation,widelyusedinhigh-gainrapidresponseexcitation.Theanglecompensationthebasicdesignofthepowersystemstabilizeristoincreasethedampingtorque,andsincethen,whichhasbecomeanimportantprincipleinthedesignofvariouspowersystemdamping,suchasthesuperiorvenacava,TCSC,UPFC.Althoughmanypaperdesigniscommittedtotheapplicationofthisprinciple,itisinterestingtonotethat,hadtobecarefultosettheparametersofthepowersystemstabilizerinrealindustrialapplications,althoughitsdesigntheoryhasbeenverymature.Therefore,thetheclassicalelectricpowersystemstabilizerdesigntheoryreallymeettheirdesignrequirements?Combinationoffrequencydomainandtimedomainanalysis,thispaperre-investigatethebasicideabehindtheclassicalpowersystemstabilizerdesigntheory.Thispaperclarifiestheconceptofthesynchronizationtorqueanddampingtorque,thatclassicpowersystemstabilizerdesignprinciplesbasedontheconceptofjusttheory.Thisarticlediscussesthebestlinearregulatordesignmethodandanalyzeitsrelationshipwiththeroutineuseoftheclient'sprocessdesign.Todoso,papersrevealthetrueprinciplesofthepowersystemstabilizer,andrecommendedtheuseofmoresystematic,moreadvancedcontroltoolstoimprovecontrollerperformance.Powersystemstabilizer,low-frequencyoscillations,powersystemdynamics.1Introduction.Inthemoderndevelopmentofpowerindustry,inordertomeettheenlargedscalepowersystemdevelopment,largegeneratorsandnetworkoperation.Therefore,thetransientstabilityhasbecomeaveryseriousproblem.Inordertosignificantlyimprovetransientstability,fastresponseandhighgainexcitationhavebeenwidelyusedinpowersystems,whichprovokesanotherseriousproblem,thelow-frequencyoscillation.Powersystemstabilizer(PSS)isanauxiliaryexcitationcontroltechnologyinordertoeliminatelow-frequencyoscillationsinthepowersystem.Itsdesigntheoryhasbeengreatlydevelopedinthe1960s.Usethisasabasisforthetheory,theoccurrenceoflow-frequencyoscillationisweakdampingpowersystemoperationornegativedamping.Sinceusuallysolvethisproblemisanincreaseinthepowersystemdamping.Doingso,andwhereinthesolenoidisbrokendownintoasynchronizingtorque,dampingtorque;fromspeedorelectromagnetictorqueorbothfeedbackisusedtoincreasethedampingtorque.Thisisacriticalstepliesinthethecompensationlaganglebetweentheturn,andtheelectromagneticturn.Doingsoanadditionalthingdampingtorqueisappendedtotheoriginaltorque.Thereisnodoubtthatthelowfrequencyoscillationisusedtoimprovetheweakdampingornegativedampingsystem.However,thequestionofwhatdampingtorque?Itismerelytheoreticaldefinition?Itcanreallyhelptoimprovethesystemtodampingtheso-calleddampingtorquetoimprove?Wefoundtherightwaytoimprovethispaperre-examinetheclassic"powersystemstabilityifwereallywanttodothis?designtheory,toclarifythesynchronizingtorque,dampingtorqueconceptthepurposeofthispaperistorevealthetruenatureofthepowersystemstabilizer,thisisnotbasedonthedampingtorquecompensation.articlefurtheranalysisofthelinearoptimizationdesigntechniquedisclosedintheproportionalintegralthecloserelationshipbetweenthedifferentialcontrolandlinearoptimalcontrol.articleisdividedintoasfollows:Section2describesthedetailedideahowtousetheclassicpowersystemstabilizerdesigntheory;part3,pointedoutthattheywillclarifythesynchronoustorque,thedampingtorquedefects;part4analysistooptimizethedesigntechnology;conclusionsandfutureresearchworkwasintroducedinsection5.2classicpowersystemstabilizerdesignprinciplesClassicpowersystemstabilizerdesigntheorytothestand-alone-theinfinitecapacityofthepowersystemforthestudyoftherelationshipequation1.Therotorofthegeneratorequationcanbedescribedas（1）andsynchronizationgeneratorshaftrotorangle,speeddeviation;synchronousangularvelocityisrad/s;2hgeneratorinertia;saidmechanicaltorqueandtorque,thedampingcoefficient,onbehalfofthefrictionanddamping.IfthegeneratorisexpressedintoinstantaneousvoltageeandtheclassicalmodeloftheinstantaneousresistanceXD.Electromagnetictorquecanbeexpressedbythefollowingformula:（2）（3）（4）Equation5andEquation3isnotdirectlyequaltothedecompositionoftheelectromagnetictorquetothesynchronizingtorqueanddampingtorque.ApplicationstabilitycriterioninEquation3andmustbepositive,and,sinbereplaced.Separatedandafter,wecametothefollowingformula:附錄B外文文獻(xiàn)譯文電力系統(tǒng)穩(wěn)定器的最新研究《科學(xué)技術(shù)》2011年十月電力系統(tǒng)穩(wěn)定器是20世紀(jì)60年代被推出解決低頻振蕩問題，廣泛用于高增益快速反應(yīng)勵磁。電力系統(tǒng)穩(wěn)定器的基本設(shè)計在于角度賠償增加阻尼力矩，從那時起，這已成為一個重要原則在設(shè)計各個電力系統(tǒng)阻尼器，例如腔靜脈，可控串補(bǔ)，統(tǒng)一潮流控制器。盡管許多紙面設(shè)計致力于了應(yīng)用這原則，很有趣注意，在真正的工業(yè)應(yīng)用是電力系統(tǒng)穩(wěn)定器的參數(shù)不得不小心設(shè)置，盡管它的設(shè)計理論已十分成熟。所以古典的電力系統(tǒng)穩(wěn)定器的設(shè)計理論真滿足其設(shè)計需求?結(jié)合頻域和時域分析，這篇論文重新調(diào)查基本的想法落后于古典的電力系統(tǒng)穩(wěn)定器設(shè)計理論。這篇論文澄清了同步轉(zhuǎn)矩以及阻尼力矩的概念，證明經(jīng)典的電力系統(tǒng)穩(wěn)定器設(shè)計原則以這些為依據(jù)概念僅僅由于理論。并且此文討論了線性最佳的調(diào)節(jié)器設(shè)計方法和分析它同常規(guī)使用客戶的進(jìn)程號設(shè)計的關(guān)系。這樣做的話，次論文揭示了電力系統(tǒng)穩(wěn)定器的真正原理，并建議使用更系統(tǒng)的，更先進(jìn)的控制工具提高控制器表現(xiàn)。電力系統(tǒng)穩(wěn)定器，低頻振蕩，電力系統(tǒng)動力學(xué)。1.介紹。在現(xiàn)代電力工業(yè)的發(fā)展中，為了滿足放大的規(guī)模電力系統(tǒng)的發(fā)展，大型發(fā)電機(jī)已經(jīng)并網(wǎng)運(yùn)行。所以暫態(tài)穩(wěn)定已成為非常嚴(yán)重的問題。為了大幅度提高暫態(tài)穩(wěn)定，快速反應(yīng)高增益勵磁開始廣泛用于電力系統(tǒng)中，但這又引起了另一嚴(yán)重問題，低頻振蕩。電力系統(tǒng)穩(wěn)定器(PSS)是一個輔助勵磁控制技術(shù)為了消除發(fā)生在電力系統(tǒng)中的低頻振蕩。它的設(shè)計理論在20世紀(jì)60年代中有了很大的發(fā)展。以這個為依據(jù)理論，低頻振蕩的發(fā)生是由于電力系統(tǒng)運(yùn)行中的弱阻尼或負(fù)阻尼。自通常解決這個問題的辦法是在電力系統(tǒng)中增加阻尼。這樣做，其中電磁被分解成同步轉(zhuǎn)矩，阻尼力矩；并且來自于轉(zhuǎn)速或電磁式轉(zhuǎn)矩或兩者的反饋用來增加阻尼力矩。這關(guān)鍵的一步在于補(bǔ)償在轉(zhuǎn)度和電磁式轉(zhuǎn)之間的滯后角。這樣做的話一個額外的東西阻尼力矩就被附加在原來的力矩上。毫無疑問的是低頻振蕩是用來提高弱阻尼或負(fù)阻尼系統(tǒng)。然而，問題是什么阻尼力矩？它僅僅是理論上的定義?它真的能幫助改善系統(tǒng)阻尼這個所謂阻尼力矩提高?或有我們發(fā)現(xiàn)合適的方式提高它如果我們真想這樣做?這篇論文系統(tǒng)重新研究經(jīng)典的“電力系統(tǒng)穩(wěn)定器的設(shè)計理論，澄清同步轉(zhuǎn)矩，阻尼力矩的概念。此文的目的是揭示電力系統(tǒng)穩(wěn)定器真正的性質(zhì)，這不基于阻尼力矩補(bǔ)償。此文進(jìn)一步分析了線性優(yōu)化設(shè)計技術(shù)，揭示在比例積分微分控制和線性最優(yōu)控制之間的密切關(guān)系。此文分為如下：部分2會說明詳細(xì)思想如何運(yùn)用的經(jīng)典的電力系統(tǒng)穩(wěn)定器設(shè)計理論;部分3會澄清同步轉(zhuǎn)矩，阻尼力矩，指出他們的缺陷;部分4分析優(yōu)化設(shè)計技術(shù);結(jié)論和未來研究工作被介紹在部分5。2經(jīng)典的電力系統(tǒng)穩(wěn)定器設(shè)計原則經(jīng)典的電力系統(tǒng)穩(wěn)定器設(shè)計理論以單機(jī)-無限大容量電力系統(tǒng)為研究對象時的關(guān)系用方程1來說明。發(fā)電機(jī)的轉(zhuǎn)子方程可以被描述為（1）和是發(fā)電機(jī)有關(guān)的同步轉(zhuǎn)軸轉(zhuǎn)子角，轉(zhuǎn)速偏差;是同步角度速度單位是rad/s;2h是發(fā)電機(jī)慣性;其中和分別表示機(jī)械轉(zhuǎn)矩和電磁轉(zhuǎn)矩，是阻尼系數(shù)，代表摩擦力和阻尼。如果發(fā)電機(jī)表示成瞬時電壓e和瞬時電阻xd的經(jīng)典模型。那么電磁轉(zhuǎn)矩可以被如下公式表示：（2）（3）（4）公式5和公式3并不直接相等在分解電磁轉(zhuǎn)矩至同步轉(zhuǎn)矩和阻尼轉(zhuǎn)矩后。應(yīng)用穩(wěn)定性準(zhǔn)則公式3中的和應(yīng)該為正的，并且，sin被換成。分開和后，我們會得出下列公式：附錄C附加圖表附錄D主要源程序RL的程序：clearallcloseallt=0;ff=50;nx=13;fn=18000;h=1/ff/80;x=[-0.5856;1.0136;1;1.2145;0;0;0;0;0.852646;0.5684;0;0;2.9049];z=zeros(nx,fn);z(:,1)=x;ak=zeros(nx,1);ad=[h/2;h/2;h;h];forfi=2:fnak=RLPSSfun(x,nx,t);b=x;forfj=1:3x=z(:,fi-1)+ak.*ad(fj);b=b+ak.*ad(fj+1)/3;ak=RLPSSfun(x,nx,t);endx=b+h*ak/6;z(:,fi)=x;t=t+h;endAngles=zeros(1,fn);Angles(1,:)=z(4,:);plot(Angles')RLfunfunctionak=RLPSSfun(x,nx,t)ak=zeros(nx,1);Te=0.05;Ke=100;Xd=2.2639;Xq=2.2639;Xd1=0.2689;Xq1=0.2689;r=0.01;Td=8.375;Tq=0.8375;Tj=5.2722;D=0;Re=0.06366;Xe=0.48785;Re1=0.122;Xe1=0.77958;Kg=20;Ts=0.3;Tch=0.2;Trh=8;a1=1/3;M0=0.852646;K0=3;T1=0.2415;T2=0.070343;T5=0.291;T6=0.0583;a=0.26027;b=-2.52995;c=-0.087357;d=0.84332;%a=0.392417;b=-5.060425;c=0.161447;d=0.051273;V=0.78731;Vr=1;Wb=314.1592;Wr=1;Ug=Wr-x(3);Ur=Wr-x(3);C4=5.9;C2=1.05;C1=1.05;ift>0&&t<=0.15Id=((1+a*r-b*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*r-b*Xd1)*(1+a*r-b*Xd1)+(a*Xd1+b*r)*(a*Xd1+b*r));Iq=(a*x(1)+b*x(2)+d*V*cos(x(4))-c*V*sin(x(4))-(1+a*r-b*Xd1)*Id)/(a*Xd1+b*r);elseift>0.15&&t<=0.65Id=((r+Re1)*(x(1)+V*sin(x(4)))-(Xd1+Xe1)*(x(2)-V*cos(x(4))))/((r+Re1)*(r+Re1)+(Xd1+Xe1)*(Xd1+Xe1));Iq=(x(1)+V*sin(x(4))-(r+Re1)*Id)/(Xd1+Xe1);elseId=((r+Re)*(x(1)+V*sin(x(4)))-(Xd1+Xe)*(x(2)-V*cos(x(4))))/((r+Re)*(r+Re)+(Xd1+Xe)*(Xd1+Xe));Iq=(x(1)+V*sin(x(4))-(r+Re)*Id)/(Xd1+Xe);endVd=x(1)-r*Id-Xd1*Iq;Vq=x(2)-r*Iq+Xd1*Id;v2=Vd*Vd+Vq*Vq;Vt=sqrt(v2);ifx(13)>C4x(13)=C4;elseifx(13)<-C4x(13)=-C4;endifx(8)+M0>C2x(8)=C2-M0;elseifx(8)+M0<-C1x(8)=-(C1-M0);endak(1)=(-x(1)-(Xq-Xq1)*Iq)/Tq;ak(2)=(x(8)-x(2)+(Xd-Xd1)*Id)/Td;ak(3)=(x(7)+a1*x(6)-D*x(3)-(x(1)*Id+x(2)*Iq-(Xq1-Xd1)*Id*Iq))/Tj;ak(4)=(x(3)-1)*Wb;ak(5)=(Kg*(T5/T6*(T5/T6*(T5/T6*Ug)))-x(5))/Ts;ak(6)=(x(5)+M0-x(6))/Tch;ak(7)=((1-a1)*x(6)-x(7))/Trh;ak(8)=(Ke*(Vr-Vt+T1/T2*(-K0*T1/T2*Ur))-x(8)+2.9049)/Te;說明：RL和RLfun程序應(yīng)該在同一個文件夾中，因?yàn)镽Lfun程序是RL的子程序，放在同一個文件夾中，以備RL調(diào)用。RLPSSclearallcloseallt=0;ff=50;nx=13;fn=18000;h=1/ff/80;x=[-0.5856;1.0136;1;1.2145;0;0;0;0;0.852646;0.5684;0;0;2.9049];z=zeros(nx,fn);z(:,1)=x;ak=zeros(nx,1);ad=[h/2;h/2;h;h];forfi=2:fnak=RLPSSfun(x,nx,t);b=x;forfj=1:3x=z(:,fi-1)+ak.*ad(fj);b=b+ak.*ad(fj+1)/3;ak=RLPSSfun(x,nx,t);endx=b+h*ak/6;z(:,fi)=x;t=t+h;endAngles=zeros(1,fn);Angles(1,:)=z(4,:);plot(Angles')RLPSSfunfunctionak=RLPSSfun(x,nx,t)ak=zeros(nx,1);Te=0.05;Ke=100;Xd=2.2639;Xq=2.2639;Xd1=0.2689;Xq1=0.2689;r=0.01;Td=8.375;Tq=0.8375;Tj=5.2722;D=0;Re=0.06366;Xe=0.48785;Re1=0.122;Xe1=0.77958;Kg=20;Ts=0.3;Tch=0.2;Trh=8;a1=1/3;M0=0.852646;K0=3;T1=0.2415;T2=0.070343;T5=0.291;T6=0.0583;a=0.26027;b=-2.52995;c=-0.087357;d=0.84332;V=0.78731;Vr=1;Wb=314.1592;Wr=1;Ug=Wr-x(3);Ur=Wr-x(3);C4=5.9;C2=1.05;C1=1.05;ift>0&&t<=0.15Id=((1+a*r-b*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*r-b*Xd1)*(1+a*r-b*Xd1)+(a*Xd1+b*r)*(a*Xd1+b*r));Iq=(a*x(1)+b*x(2)+d*V*cos(x(4))-c*V*sin(x(4))-(1+a*r-b*Xd1)*Id)/(a*Xd1+b*r);elseift>0.15&&t<=0.65Id=((r+Re1)*(x(1)+V*sin(x(4)))-(Xd1+Xe1)*(x(2)-V*cos(x(4))))/((r+Re1)*(r+Re1)+(Xd1+Xe1)*(Xd1+Xe1));Iq=(x(1)+V*sin(x(4))-(r+Re1)*Id)/(Xd1+Xe1);elseId=((r+Re)*(x(1)+V*sin(x(4)))-(Xd1+Xe)*(x(2)-V*cos(x(4))))/((r+Re)*(r+Re)+(Xd1+Xe)*(Xd1+Xe));Iq=(x(1)+V*sin(x(4))-(r+Re)*Id)/(Xd1+Xe);endVd=x(1)-r*Id-Xd1*Iq;Vq=x(2)-r*Iq+Xd1*Id;v2=Vd*Vd+Vq*Vq;Vt=sqrt(v2);ifx(13)>C4x(13)=C4;elseifx(13)<-C4x(13)=-C4;endifx(8)+M0>C2x(8)=C2-M0;elseifx(8)+M0<-C1x(8)=-(C1-M0);endak(1)=(-x(1)-(Xq-Xq1)*Iq)/Tq;ak(2)=(x(13)-x(2)+(Xd-Xd1)*Id)/Td;ak(3)=(x(10)+a1*x(9)-D*x(3)-(x(1)*Id+x(2)*Iq-(Xq1-Xd1)*Id*Iq))/Tj;ak(4)=(x(3)-1)*Wb;ak(5)=((1-T5/T6)*Ug-x(5))/T6;ak(6)=((1-T5/T6)*(x(5)+T5/T6*Ug)-x(6))/T6;ak(7)=((1-T5/T6)*(x(6)+T5/T6*(x(5)+T5/T6*Ug))-x(7))/T6;ak(8)=(Kg*(x(7)+T5/T6*(x(6)+T5/T6*(x(5)+T5/T6*Ug)))-x(8))/Ts;ak(9)=(x(8)+M0-x(9))/Tch;ak(10)=((1-a1)*x(9)-x(10))/Trh;ak(11)=(-K0*(1-T1/T2)*Ur-x(11))/T2;ak(12)=((1-T1/T2)*(x(11)-K0*T1/T2*Ur)-x(12))/T2;ak(13)=(Ke*(Vr-Vt+x(12)+T1/T2*(x(11)-K0*T1/T2*Ur))-x(13)+2.9049)/Te;說明：同樣RLPSS和RLPSSfun在同一個文件夾中。5．myfunfunctionfit=myfun(pop)t=0;ff=50;nx=13;fn=18000;h=1/ff/80;x=[-0.5856;1.0136;1;1.2145;0;0;0;0;0.852646;0.5684;0;0;2.9049];z=zeros(nx,fn);z(:,1)=x;ak=zeros(nx,1);ad=[h/2;h/2;h;h];forfi=2:fnak=myfunfun(x,nx,t,pop);b=x;forfj=1:3x=z(:,fi-1)+ak.*ad(fj);b=b+ak.*ad(fj+1)/3;ak=myfunfun(x,nx,t,pop);endx=b+h*ak/6;z(:,fi)=x;t=t+h;endfit=[z(4,:)-1.2241761]*[(z(4,:)-1