电力系统稳定器pss设计与仿真本科毕业论文

电力系统稳定器PSS的设计与仿真摘要由于电力系统在正常运行时会发生频率的振荡，对我们的生产生活带来了很多的危害，给我国的国民经济造成了巨大的损失。如果在电力系统中加入PSS后会对系统的稳定性提高给予了很大的帮助。电力系统稳定器就是为抑制低频振荡而研究的一种附加励磁控制技术。它在励磁电压调节器中，引入领先于轴速度的附加信号，产生一个正阻尼转矩，去克服原励磁电压调节器中产生的负阻尼转矩作用,用于提高电力系统阻尼，解决低频振荡问题，是提高电力系统动态稳定性的重要措施之一。本文首先分析了电力系统稳定器的研究背景和国内外研究状况，然后又对无限大系统的数学模型进行了具体的分析，然后又分析了系统带PSS和不带PSS时系统的运行情况，最后对两种情况进行了比较，分析了PSS的优点和缺点，同时对电力系统的未来提出了更高的挑战。关键词：励磁控制技术;低频振荡；电力系统稳定器DesignandSimulationofPowerSystemStabilizer(PSS)ABSTRACTThispaperdescribesthepowersystemstabilizer'sdesignprinciplesandtheimpactonthepowersystem,theuseofthebenefitsofpowersystemstabilizer.ThefrequencyofoscillationmaybehappentopowersystemwhenthesystemInthenormaloperation.ItgivesagreathelpifthepowersystemaftertheadditionofPSS.Powersystemstabilizer(PSS)istosuppressalowfrequencyoscillationofadditionalexcitationcontrol.Itistheexcitationvoltageregulator,theintroductionofaxialvelocityaheadofadditionalsignalstoproduceapositivedampingtorquetoovercometheprimaryexcitationvoltageregulatorproducednegativedampingtorqueeffect.Improvingpowersystemdamping,lowingfrequencyoscillationproblemsolvingistoimprovepowersystemdynamicstabilityoftheimportantmeasures.Thispaperanalyzesthepowersystemstabilizerresearchbackgroundandresearchstatusathomeandabroad,andthentotheinfinitesystemofmathematicalmodelsofspecificanalysis,andthenanalyzedthesystemwithPSSandwithoutPSSoperationofthesystem,andfinallyOfthetwocaseswerecomparedandanalyzedtheadvantagesanddisadvantagesofPSS,whilethefutureofthepowersystemofahigherchallenge.Keyword:Excitationcontrol；Low-frequencyoscillation；PowerSystemStabilizer毕业设计（论文）原创性声明和使用授权说明原创性声明本人郑重承诺：所呈交的毕业设计（论文），是我个人在指导教师的指导下进行的研究工作及取得的成果。尽我所知，除文中特别加以标注和致谢的地方外，不包含其他人或组织已经发表或公布过的研究成果，也不包含我为获得及其它教育机构的学位或学历而使用过的材料。对本研究提供过帮助和做出过贡献的个人或集体，均已在文中作了明确的说明并表示了谢意。作者签名：日期：指导教师签名：日期：使用授权说明本人完全了解大学关于收集、保存、使用毕业设计（论文）的规定，即：按照学校要求提交毕业设计（论文）的印刷本和电子版本；学校有权保存毕业设计（论文）的印刷本和电子版，并提供目录检索与阅览服务；学校可以采用影印、缩印、数字化或其它复制手段保存论文；在不以赢利为目的前提下，学校可以公布论文的部分或全部内容。作者签名：日期：

学位论文原创性声明本人郑重声明：所呈交的论文是本人在导师的指导下独立进行研究所取得的研究成果。除了文中特别加以标注引用的内容外，本论文不包含任何其他个人或集体已经发表或撰写的成果作品。对本文的研究做出重要贡献的个人和集体，均已在文中以明确方式标明。本人完全意识到本声明的法律后果由本人承担。作者签名： 日期：年月日学位论文版权使用授权书本学位论文作者完全了解学校有关保留、使用学位论文的规定，同意学校保留并向国家有关部门或机构送交论文的复印件和电子版，允许论文被查阅和借阅。本人授权大学可以将本学位论文的全部或部分内容编入有关数据库进行检索，可以采用影印、缩印或扫描等复制手段保存和汇编本学位论文。涉密论文按学校规定处理。作者签名： 日期：年月日导师签名：日期：年月日

注意事项1.设计（论文）的内容包括：1）封面（按教务处制定的标准封面格式制作）2）原创性声明3）中文摘要（300字左右）、关键词4）外文摘要、关键词5）目次页（附件不统一编入）6）论文主体部分：引言（或绪论）、正文、结论7）参考文献8）致谢9）附录（对论文支持必要时）2.论文字数要求：理工类设计（论文）正文字数不少于1万字（不包括图纸、程序清单等），文科类论文正文字数不少于1.2万字。3.附件包括：任务书、开题报告、外文译文、译文原文（复印件）。4.文字、图表要求：1）文字通顺，语言流畅，书写字迹工整，打印字体及大小符合要求，无错别字，不准请他人代写2）工程设计类题目的图纸，要求部分用尺规绘制，部分用计算机绘制，所有图纸应符合国家技术标准规范。图表整洁，布局合理，文字注释必须使用工程字书写，不准用徒手画3）毕业论文须用A4单面打印，论文50页以上的双面打印4）图表应绘制于无格子的页面上5）软件工程类课题应有程序清单，并提供电子文档5.装订顺序1）设计（论文）2）附件：按照任务书、开题报告、外文译文、译文原文（复印件）次序装订3）其它

目录摘要 )根据以上的系统模型公式可以运用MATLAB编程，得到程序附录3，所以得到以下的功角曲线：图5-3功角特性曲线图5-4变化曲线5.4主要结论本章主要分析了电力系统的特点，说明了现代的电力系统存在的种种问题，提出提高电力系统稳定性的重要性。通过分析了国内外多种提高电力系统稳定性的方法，总结出PSS的应用范围比较广泛，而且取得了不错的效果。所以我们通过对无限大系统母线在附加PSS和不加PSS两种状态下系统运行情况的比较，得出了PSS在抑制系统低频振荡方面具有明显的效果，通过两幅图的比较也可以看出在附加PSS的系统振荡会很快的趋于平缓，而未加PSS的系统则一直振荡。通过以上的分析我们又总结出了PSS的一些不足点：首先，PSS的参数的整定是在某一个低频振荡的条件下计算的，如果频率振荡更加严重，在低频振荡的环境下计算出来的数据有可能在其他的情况下不是很适用；还有PSS的安装地点应该也是一个比较重要的问题；另外PSS的种类繁多，对一个机组装PSS时有必要做测试，这样也浪费了大量的时间，从经济方面考虑又存在这问题。随着世界经济的发展，人们对供电的可靠性和稳定性日益提高。同时电力网络的开发和商业化的运营使得电力系统的运行越来越趋于极限，电力系统不稳定的因素越来越多，特别是电压的不稳和低频振荡问题更加严重，所以提出更有效的措施来解决问题可以说是当务之急。致谢经过半年的忙碌和工作，本次毕业设计已经接近尾声，作为一个本科生的毕业设计，由于经验的匮乏，难免有许多考虑不周全的地方，如果没有导师的督促指导，以及一起工作的同学们的支持，想要完成这个设计是难以想象的。在这里首先要感谢我的老师向老师。向老师平日里工作繁多，但在我做毕业设计的每个阶段，从查阅资料到设计草案的确定和修改，中期检查，后期详细设计等整个过程中都给予了我悉心的指导。我的设计较为复杂烦琐，但是向老师仍然细心地纠正模型中的错误。除了敬佩向老师的专业水平外，她的治学严谨和科学研究的精神也是我永远学习的榜样，并将积极影响我今后的学习和工作。其次要感谢我的同组同学对我无私的帮助，特别是在MATLAB软件的使用方面，正因为如此我才能顺利的完成设计，我要感谢我的母校——江西理工大学应用科学学院，是母校给我们提供了优良的学习环境；另外，我还要感谢那些曾给我授过课的每一位老师，是你们教会我专业知识。在此，我再说一次谢谢！谢谢大家！！！。

参考文献[1]韩祯祥．电力系统稳定[M]．北京．中国电力出版社．1995，56-62．[2]倪以信．动态电力系统理论与分析[M]．北京．清华大学出版社．2002，3-6，15-16．[3]孙伯明．高等电力网络分析[M]．北京．清华大学出版社．1996．[4]陈珩．电力系统稳态分析[M]．北京．中国电力出版社．1995．[5]李光琦．电力系统暂态分析[M]．北京．中国电力出版社．1995．[6]王锡凡．现代电力系统分析[M]．北京．现代出版社．2003，333-342．[7]刘晓鹏．小干扰稳定性部分特征值分析的多重变换法[J]．电力系统自动化．1998，22（9）：38-42．[8]苏金明．MATLAB7.0实用指南[M]．北京．电子工业出版社．2004．[9]彭淦榕．MATLAB在电力系统分析中的应用[J]．山东电力技术．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外文文献原文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外文文献译文电力系统稳定器的最新研究《科学技术》2011年十月电力系统稳定器是20世纪60年代被推出解决低频振荡问题，广泛用于高增益快速反应励磁。电力系统稳定器的基本设计在于角度赔偿增加阻尼力矩，从那时起，这已成为一个重要原则在设计各个电力系统阻尼器，例如腔静脉，可控串补，统一潮流控制器。尽管许多纸面设计致力于了应用这原则，很有趣注意，在真正的工业应用是电力系统稳定器的参数不得不小心设置，尽管它的设计理论已十分成熟。所以古典的电力系统稳定器的设计理论真满足其设计需求?结合频域和时域分析，这篇论文重新调查基本的想法落后于古典的电力系统稳定器设计理论。这篇论文澄清了同步转矩以及阻尼力矩的概念，证明经典的电力系统稳定器设计原则以这些为依据概念仅仅由于理论。并且此文讨论了线性最佳的调节器设计方法和分析它同常规使用客户的进程号设计的关系。这样做的话，次论文揭示了电力系统稳定器的真正原理，并建议使用更系统的，更先进的控制工具提高控制器表现。电力系统稳定器，低频振荡，电力系统动力学。1.介绍。在现代电力工业的发展中，为了满足放大的规模电力系统的发展，大型发电机已经并网运行。所以暂态稳定已成为非常严重的问题。为了大幅度提高暂态稳定，快速反应高增益励磁开始广泛用于电力系统中，但这又引起了另一严重问题，低频振荡。电力系统稳定器(PSS)是一个辅助励磁控制技术为了消除发生在电力系统中的低频振荡。它的设计理论在20世纪60年代中有了很大的发展。以这个为依据理论，低频振荡的发生是由于电力系统运行中的弱阻尼或负阻尼。自通常解决这个问题的办法是在电力系统中增加阻尼。这样做，其中电磁被分解成同步转矩，阻尼力矩；并且来自于转速或电磁式转矩或两者的反馈用来增加阻尼力矩。这关键的一步在于补偿在转度和电磁式转之间的滞后角。这样做的话一个额外的东西阻尼力矩就被附加在原来的力矩上。毫无疑问的是低频振荡是用来提高弱阻尼或负阻尼系统。然而，问题是什么阻尼力矩？它仅仅是理论上的定义?它真的能帮助改善系统阻尼这个所谓阻尼力矩提高?或有我们发现合适的方式提高它如果我们真想这样做?这篇论文系统重新研究经典的“电力系统稳定器的设计理论，澄清同步转矩，阻尼力矩的概念。此文的目的是揭示电力系统稳定器真正的性质，这不基于阻尼力矩补偿。此文进一步分析了线性优化设计技术，揭示在比例积分微分控制和线性最优控制之间的密切关系。此文分为如下：部分2会说明详细思想如何运用的经典的电力系统稳定器设计理论;部分3会澄清同步转矩，阻尼力矩，指出他们的缺陷;部分4分析优化设计技术;结论和未来研究工作被介绍在部分5。2经典的电力系统稳定器设计原则经典的电力系统稳定器设计理论以单机-无限大容量电力系统为研究对象时的关系用方程1来说明。发电机的转子方程可以被描述为（1）和是发电机有关的同步转轴转子角，转速偏差;是同步角度速度单位是rad/s;2h是发电机惯性;其中和分别表示机械转矩和电磁转矩，是阻尼系数，代表摩擦力和阻尼。如果发电机表示成瞬时电压e和瞬时电阻xd的经典模型。那么电磁转矩可以被如下公式表示：（2）（3）（4）公式5和公式3并不直接相等在分解电磁转矩至同步转矩和阻尼转矩后。应用稳定性准则公式3中的和应该为正的，并且，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程序应该在同一个文件夹中，因为RLfun程序是RL的子程序，放在同一个文件夹中，以备RL调用。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