《完整版(2022年)MATLAB进行控制系统的超前设计—课程设计说明书.docx》由会员分享,可在线阅读,更多相关《完整版(2022年)MATLAB进行控制系统的超前设计—课程设计说明书.docx(27页珍藏版)》请在优知文库上搜索。
1、完整版(2022年)用MATLAStt行控制系统的超前校正设计初始条件:Y已知一单位反馈系统的开环传承函数右二G(三)s(10。5S)(Io.5s)要求系统跟随2rmin的斜坡输入产生的最大稳态误差为2。,45o要求完成的主要任务:(包括课程设计工作量及其技术要求,以及说明书撰写等具体要求)1、用MATLA酢出满足初始条件的最小K值的系统伯德图,计算系统的幅值裕量和相位裕量。2、在系统前向通路中插入一相位超前校正,确定校正网络的传递函数P37用EIKnZAfflE出未校正和已校正系统的根轨迹。4、课程设计说明书中要求写清楚计雏分析的过程,列出MATSB程序和MATLA的出。说明书的格式按蕊处标
2、准书写。时间安排:任务时间(天)审题、查阅相关资料1分析、计算2编写程序2撰写报告2论文答辩s指导教师签名:武汉理工大学自动控制原理课程设计说明书系主任(或责任教师)签名:年月日控制系统的超前校正设计1校正简介以及超前校正的原理方法1.1 校正简介所谓校正就是在系统不可变部分的基础上,加入适当的校正元部件,使系统满足所给定的性能指标。校正环节的形式及其在系统中的位置称为校正方案。一般有:串联校正,并联校正,反馈校正,以及前馈校正,本次课设就是用的串联校正中的超前校正。1.2 超前校正及其特性超前校正就是在前向通道中串联传递函数为:GcsSIaTSlRsaTs1式中:aRlRz1R2RIR2CI
3、ORlR武汉理工大学自动控制原理课程设计说明书通常a为分度系数,T叫时间常数,由式(2-1)可知,采用无方超前网络进行用联校正时,整个系统的开环增益要下降a倍,因此需要提高放大器增益交易补偿.如果对无源超前网络传递函数的衰减由放大器增益所补偿,则aTs1aGcsTs1上式称为超前校正装置的传递函数。无源超前校正网络的对数频率特性如图1。()=图1无源超前校正网络的对数频率特性超前校正对频率在1/aT和1/T之间的输入信号有微分作用,在该频率范围内,输出信号相角比输入信号相角超前,超前网络的名称由此而得。因此超前校正的基本原理就是利用超前相角补偿系统的滞后相角,改善系统的动态性能,如增加相位裕度
4、,提高系统的稳定性等。下面先求取超前校正的最大超前相角fn及取得最大超前相角的频率一则像频特性:C=arctanaT-arctaTdcaTT22d1aT1T当10,则有:marctanaT从而有:1既当m=-r-.三nYaT时,超前相角最大为m=arcsinJ,可以看出a+19只与a有关这一点对于超前校正是相当重要的超前校正RC网络图如图2oRC图2无源超前校正RC网络图1.3参数的一般设计步骤=Y-Y+(1)根据给定的系统性能指标,确定开环增益Ko+利盛务前1中求的的K绘制耒校正系统的伯德图。在伯德图上量取未校正系统的相位裕度和幅值裕度,并计算为使相位裕度达到。其中为给定的相位裕度指标,。为
5、未校正系给定的指标所需补偿的超前相角统的相位裕度,4附加的角度。(4)取1sin1-sin求出a。m(5)取未校正系统的幅值为-10lga(dB)时的频率作为校正后系统的截止频率武汉理工大学自动控制原理课程设计说明书(6)由m1计算出参数T,并写出超前校正传递函数。maT(7)检验指标:绘制系统校正后的伯德图,检验是否满足给定的性能指标。当系统仍不能满足要求时增大值,从步骤3开始重新计算。2超前校正的设计2.1 校正前系统分析由于要求要求系统跟随2rmin的斜坡输入,所以:R=12M60=一R2QsK取K最小值K=6。则待校正的系统的宾标传递函数为$)S10.05s10.5s该最小相位系统的伯
6、德图matlab程序如下:num=6;den=conv(0.051,0.51)0;BodeDiagramFrequency(radsec图3伯德图为了求出校正前的相角裕度和幅值裕度,在matlab中输入程序:G=tf(6,0.0250.5510);margin(G).BodeDiagramo3dB(at6.32rad/sec),Pm-23.3deg(at3.17rad/sec)再输入:3580251-12-27010Frequency(rad/sec)图4相频特性曲线得到:武汉理工大学自动控制原理课程设计说明书wg,wc=margin(G)G=tf(6,0.025相角裕度r=23.29200.
7、55 1O)kg,r,穿越频率Wg=6.3246截止频率wc=3.1654显然45。,需进行超前校正。再画出其跟轨迹,程序如下:n=6;d=0,025,0.55,1,0;rlocus(n,d)得到跟轨迹图如图5o图5根轨迹图2.2 校正系统设计由于给定的相位裕度指标为45 ,未校正系统的相位裕度为23.3 ,不妨设附加角度为15.3。则:37,从而求出a1sn求出a:1-Slm1sinm.a41sinrn设校正后的截所以可得超前网络函数为:C(.1+0.448sGcs加入校正环节之后的传递函数为:()()=GsGs10.112s610.448ss10.05s10.112s10.5s在计算之后还
8、可用其它的方法来进行检验,看所加装置参数的选择是否真的符合题意,满足要求下面用MATLAB来进行验证用MATLAB求校正后的相角裕度和幅值裕度程序为:num1=6*0.4481;den1=conv(0.11210,conv(0.051,0.51);bode(num1,den1)grid6月-50-100-135-180-225-270BodeDiagramGm三14.3dB(at13.1rad/sec)lPm三46.6Cfeg(at4.71rad/sec)-101231010101010Frequency(rad/sec)图6校正后相角裕度和幅值裕度校正后相角裕度46.645,截止频率为:4.
9、71rads满足课设要求。在用matlab画出其伯德图程序为:num1=6*0.4481;den1=0.00280.0860.6621O;bode(num1,de1)grid得到系统的伯德图见(图7卜用Matlabden1gridmp)db三(6。P)*eud10*60BodeDaoramRootLocus0RealAxis图8校正后根轨迹图3校正前后系统比较用Matlab软件作系统校正前的奈奎斯特曲线的程序为:num=6den=0.0250.5510;nyquist(num,den)grid从而得到系统校正前的奈奎斯特曲线(图9)如下。图9校正前的奈奎斯特曲线校正后的奈奎斯特曲线程序如下:n
10、um=2.6886den=0.00280.0860.66210;nyquist(num,den)grid得出校正后的奈奎斯特曲线如图10。图10校正后的奈奎斯特曲线由上图可以看出来,系统开环传递函数无右极点,其奈奎斯特曲线都不包括(,Oj)点,所以闭环系统是稳定的。校正后使开环系统截止频率增大,从而闭环系统带宽也增大,使响应速度加快。运用matlab软件作系统校正前后的单位冲击响应曲线比较,程序为:num1=6;den1=0.025,0.55,1,0;num3=2,526,6;de3=0.0028l0.086,0.662,1,0;t=0:0.02:5;numc1,denc1=cloop(num
11、1,den1);y1=step(numd,dend,t);numc3,denc3=cloop(num3,de3);y3=step(numc3,denc3,t);plot(t,y1,y3);grid;title(,校正前后阶跃响应对比图,);xlabel(t(sec),);ylabel(,c(t),);gtext(,校正前,);gtext(,校正后);得到校正前后阶跃响应对比图(图11)o校正前后阶跃响应对比图图11校正前后阶跃响应对比图从图中可以看出系统在校正后,系统超调量,调节时间都相对下降,系统很快就稳定了。4心得体会这次的自动控制原理课程设计虽然时间不是很充足,但恰好由于这种原因,无形之中就有一种压力,迫使你提高效率。然后将前后知识连贯起来,将问题解决,虽说没有手到擒来的感觉,但也在这种效率之中,找到了自信,找到了乐趣。通过本次课程设计,书本中有关超前校正,开环增益,相位裕度,幅值裕度等基本概念有了更深刻的理解,同时通SMATLAB仿真将奈奎斯特曲线,伯德图,像频特性等曲线画出来,也看到了MATLAB功能的强大。参考文献1胡寿松.自动控制原理(第五版).北京:科学出版社,2007.20032刘金琨.先进PID控制及其MATLA昉真M.北京:电子工业出版社,3孙亮.自动控制原理(第三版).北京:高等教育出版社,2011