1、ify(i)gt*yss,n=;end;end;end;t=t(i);cs=length(t);j=cs+;n=;whilen==;j=j;ifkey==;ifj==;n=;elseify(j)t;ts=t;end;elseiftgtt;iftltt;ts=t;elsets=t;endend函数dist的源程序如下:fuction[detac,t,tv]=dist(key,y,t)%MATLABFUNCTIONPROGRAMdistm%%Countsgmandt[m,t]=min(y);cs=length(t);yss=y(cs);detac=mt=t(tf)tm=max(t);%Counttsi=cs+;n=;whilen==,i=i;ifkey==,ifi==,n=;elseify(i)gt=,n=;end;elseifkey==,ifi==,n=;elseify(i)gt=,n=;end;endend;t=t(i)ics=length(t);j=cs+;n=;whilen==,j=j;ifkey==ifj==,n=;elseify(j)lt=(),n=;end;elseifkey==,ifj==,n=;elseify(j)ttv=tendelseiftgtt,iftltt,ift~=tv=telsetv=tendelseift~=tv=telsetv=tendendend图形对象句柄命令figure命令figure功能创建一个新的图形对象。图形对象为在屏幕上单独的窗口,在窗口中可以输出图形。用法figure用缺省的属性值创建一个新的图形对象。figure('ProertyName',ProertyValue,)对指定的属性ProertyName用指定的属性值,ProertyValue(属性名与属性值成对出现)创建一个新的图形窗口,对于那些没有指定的属性调量反映系统的相对稳定性。超调量越小,相对稳定性越好,即动态响应比较平稳。。
2、对比两个阶跃扰动响应性能指标值可知,对作用在点的扰动信号的抗扰作用比对作用在点的扰动信号的抗扰作用强,再次表明扰动作用点离被调量越近,转速环对扰动信号的抗扰能力越好。单位冲激信号扰动的响应曲线在图A中两个扰动信号作用点分别施加单位冲激信号,绘制其扰动响应曲线。在程序文件方式下执行以下MATLAB程序Lfm:%MATLABPROGRAMLfm[a,b,c,d]=linmod('mxb');figure();imulse(a,b,c,d);holdon[a,b,c,d]=linmod('mxc');figure();imulse(a,b,c,d)程序执行后,可得转速环在扰动作用点与扰动作用点施加单位冲激信号时的扰动响应曲线如图(a)与(b)所示。(a)(b)图转速环单位冲激信号扰动响应曲线Figrotationloounitimulsesignalturbulenceresonsecurve转速环单位冲激扰动响应过程如图(a)所示呈单调衰减或者如图(b)所示呈一次衰减振荡在s内就结束了,抗扰动过程亦非常短暂,系统转速外环已经感觉到扰动,但经s扰动响应的过程已经结束。说明系统的抗扰性能好转速环频域分析的MATLAB计算及仿真频域分析是运用闭环系统的开环频率特性曲线来分析闭环系统的响应与性能的。频域分析的主要内容是绘制Bode图与计算频域性能指标。转速闭环系统的开环结构图如图a(未画出转速环给定平衡滤波环节),此图即为mxdmdl动态结构图,对mxdmdl进行仿真操作如下:图a转速闭环系统的开环结构图Figarotationclosedloosystemslitloostructurediagram根据题目要求,用函数命令linmod()与margin()给出MATLAB指令如下。在程序文件方式下执行以下MATLAB程序Lgm:%MATLABPROGRAMLgm[a,b,c,d]=linmod('mxd');sys=ss(a,b,c,d)。
3、标数据可见转速环的阶跃响应超调量σ%lt%,峰值时间tlts,调节时间stlts,表明系统的相对稳定性好,动态响应快,这样的系统响应是非常理想的。绘制单位阶跃信号扰动响应曲线并计算其性能指标在图A中两个扰动信号作用点分别施加单位阶跃信号,绘制其扰动响应曲线,并求其最大动态降落与最大动态降落的时间及恢复时间。为进行仿真,在作用点与施加扰动信号时分别构成MATLAB里的动态结构图mxbmdl与mxcmdl。这两个结构图模型中的数据和图A相同,只是扰动信号作用点不同。在程序文件方式下执行以下用linmod()与ste()函数并调用dist()的程序Lem:%MATLABPROGRAMLem[a,b,c,d]=linmod('mxb');s=ss(a,b,c,d);t=[::];[y,t]=ste(s,t);figure();ste(s,t);holdon[detac,t,tv]=dist(,y,t)[detac,t,tv]=dist(,y,t)[a,b,c,d]=linmod('mxc');s=ss(a,b,c,d);t=[::];[y,t]=ste(s,t);figure();ste(s,t);[detac,t,tv]=dist(,y,t)[detac,t,tv]=dist(,y,t)程序执行后,可得转速环的单位阶跃扰动响应曲线如图所示。(a)与(b)图分别对应着扰动信号作用于与两个不同的点。对于图(a),计算其阶跃扰动响应性能指标:最大动态降落:detac=;最大动态降落时间:t=s;基准值%范围的恢复时间:vt=s;基准值%范围的恢复时间:vt=s。对于图(b)计算其阶跃扰动响应性能指标:最大动态降落:detac=;最大动态降落时间:t=s;基准值%范围的恢复时间:vt=s;基准值%范围的恢复时间:vt=s。由阶跃扰动响应性能指标值可知,该系统对扰动信号具有良好的动态抗扰作用,对比两个阶跃扰动响应性能指标值可知,对作用在点的扰动信。
4、时间非常短,系统几乎察觉不出扰动的作用。说明电流环的抗扰性能好。电流环频域分析的MATLAB计算及仿真根据自动控制原理,频域分析的特点是运用闭环系统的开环频率特性曲线来分析闭环系统的响应及其性能。频域分析的主要内容是画Bode图与计算频域性能指标。电流闭环系统的开环结构图如图a所示,它对应着Simulink动态结构图模型mxdmdl图a电流闭环系统的开环结构图Figacurrentclosedloosystemslitloostructurediagram在程序文件方式下执行以下的MATLAB程序Lcm:%MATLABPROGRAMLcmn=;d=[];s=tf(n,d);n=[];d=[];s=tf(n,d);n=;d=[];s=tf(n,d);n=;d=[];s=tf(n,d);n=;d=[];s=tf(n,d);sys=s*s*s*s*s;margin(sys)执行语句后,可得电流环的Bode图如图b所示,在图上就附有经计算出电流环的频域性能指标:模稳定裕度hL=dBπ穿越频率gW=rads相稳定裕度γ=剪切频率cW=rads图b电流环的Bode图FigbcurrentlooofBodediagram工程上,一般要求模稳定裕度hL≥dB,相稳定裕度γ≥。可见电流环有足够的稳定裕量,其频域性能是优良的,反映了电流环具有良好的相对稳定性。转速环的MATLAB计算及仿真转速环校正前后给定阶跃响应的MATLAB计算及仿真与电流环添加低通滤波器措施一样,在转速环反馈通道与给定信号通道都添加了滤波惯性环节,其传递函数为:sTon按照把转速环设计成II型系统的要求,根据自动控制理论,转速调节器也应选择比例积分调节器,其传递函数为:ASRW(s)=nKsTsTnn,nT转速环的校正主要是对转速调节器的参数nK以及对晶闸管整流与移相触发装置的放大倍数sK进行校正。校正前nK=,sK=,转速调节器的传递函数为:ASRW(s)=nKs。
5、的抗扰作用比对作用在点的扰动信号的抗扰作用强,再次表明扰动作用点离被调量越近,转速环对扰动信号的则用缺省值。Figure(h)若h为一个已经存在的图形的句柄,则figure(h)使由h标记的图形成为当前图形,使它可见,且在屏幕上把它显示到所有图形之前。当前图形为图像输出的地方。若h不是已经存在图形的句柄,但是为一整数,则该命令生成一图形窗口,同时把该窗口的句柄赋值为h;若h不是一图形窗口的句柄,也不是一整数,则返回一错误信息。h=figure()返回图形窗口对象的句柄给h。参考文献:[]陈伯时电力拖动自动控制系统(第版)机械工业出版社,[]杨庚辰自动控制原理西安电子科技大学出版社,[]黄忠霖,周向明控制系统MATLAB计算及仿真国防工业出版社,[]黄忠霖,周向明控制系统MATLAB计算及仿真实训国防工业出版社,致谢在此特别感谢我的指导老师在整个毕业设计过程中对我的精心指导和帮助,使我能够顺利按时完成这次毕业设计。在整个毕业设计期间,导师对我的研究工作给予了细致耐心的指导,他谦虚严谨朴实的学风让我看到了自己身上的很多不足,其渊博的知识、敏锐的思维以及丰富的经验给我留下了极其深刻的印象。值此论文完成之际,谨向辛勤培育我的导师致以衷心的感谢和最诚挚的祝福。扰能力越好。(a)(b)图转速环的单位阶跃扰动响应曲线Figrotationloounitsteturbulenceresonsecurve(a)与(b)图分别对应着扰动信号作用于与两个不同的点。对于图(a),计算其阶跃扰动响应性能指标:最大动态降落:detac=;最大动态降落时间:t=s;基准值%范围的恢复时间:vt=s;基准值%范围的恢复时间:vt=s。对于图(b)计算其阶跃扰动响应性能指标:最大动态降落:detac=;最大动态降落时间:t=s;基准值%范围的恢复时间:vt=s;基准值%范围的恢复时间:vt=s。由阶跃扰动响应性能指标值可知,该系统对扰动信号具有良好的动态抗扰作用,。
6、margin(sys)语句指令执行后可得Bode图如图b所示。图b电流环的Bode图FigotationlooofBodediagram从由图上可得转速环频域性能指标:模稳定裕度hL=dBπ穿越频率gW=rads相稳定裕度γ=剪切频率cW=rads数据显示,模稳定裕度hLgtdB,相稳定裕度γgt,说明控制系统有较大的稳定裕度,反映了该系统具有良好的相对稳定性。总结通过对电流环内的触发整流装置的放大系数sK的校正,利用MATLAB及其中的仿真工具Simulink,对所设计的电流环和转速环的阶跃信号进行了仿真计算,使电流环阶跃响应Simulink曲线的峰值时间t缩短,提高系统的加速起动性能,通过对转速环内PI调节器比例部分放大系数nK及触发整流装置的放大系数sK的校正,使转速环阶跃响应性能指标的超调量σ%dB,相稳定裕度γgt,电流环和转速环都有足够的稳定裕度,符合实际工程的设计要求,该系统有良好的相对稳定性。由仿真计算结果表明,利用MATLAB的Simulink对双闭环直流调速系统进行仿真设计,可以迅速直观地分析出系统的跟随性能、抗扰性能及稳定性,使得对系统进行分析、设计及校正变得更简单方便,大大缩短了系统的调试周期,提高了开发系统的效率。对于调速系统的设计,MATLAB的Simulink确实是个经济、简单、快速、高效的工具。附录部分函数的源程序函数erf的源程序如下:function[sigma,t,ts]=erf(key,y,t)%MATLABFUNCTIONPROGRAMerfm%%Countsigmaandt[m,tf]=max(y);cs=length(t);yss=y(cs);sigma=*(myss)ysst=t(tf);%Counttsi=cs+;n=;whilen==;i=i;ifkey==;ifi==;n=;elseify(i)gt*yss,n=;end;elseifkey==;ifi==;n=;else。
7、TsTnn=s,nT构成动态结构图模型mxmdl;校正后nK=,sK=,转速调节器的传递函数为:ASRW(s)=nKsTsTnn=ss,nT构成动态结构图模型mxamdl。其他参数不变,校正前、后的动态结构图模型只是nK与sK的值不一样,所以在此只给出校正后的mxamdl文件的动态结构图的模型,如图A所示。图A带参数双闭环系数的Simulink动态结构图模型mxmdlFigAthebeltarameterdoubleclosedloocoefficientSimulinkdynamicstructurediagrammodelmxmdl用linmod()与ste()函数命令并调用函数erf()编写求其阶跃响应与性能指标的MATLAB程序Ldm如下。在程序文件方式下执行以下的MATLAB程序Ldm:%MATLABPROGRAMLdm[a,b,c,d]=linmod('mx');s=ss(a,b,c,d);figure();ste(s,'k')holdon[y,t]=ste(s);[a,b,c,d]=linmod('mxa');s=ss(a,b,c,d);figure();ste(s,'k')[y,t]=ste(s);[sigma,t,ts]=erf(,y,t);[sigma,t,ts]=erf(,y,t);程序执行后,可得转速环校正前后的单位阶跃响应曲线如图B所示。(a)(b)图B转速环阶阶跃响应Simulink曲线FigBtherotationloosteresonsecurve校正前转速环的单位阶跃响应仿真曲线如图B(a)所示,此时响应产生了衰减振荡,表明系统稳定性差。程序运行后,可得到图B(a)阶跃响应性能指标:超调量:σ%=%;峰值时间:t=s;调节时间:st=s。同时可得到图B(b)阶跃响应性能指标:超调量:σ%=%;峰值时间:t=s;调节时间:st=s。调节时间数值均对应于%的误差带。由计算图B(b)的性能指。