转子绕线机控制系统的滞后校正设计

课程设计课程设计任务书学生姓名：化学院题目:转子绕线机控制系统的滞后校正设计初始条件：已知转子绕线机控制系统的开环传递函数为:要求系统的静态速度误差系数151sKv，相位裕度56。要求完成的主要任务:（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求）1、用MATLAB作出满足初始条件的最小K值的系统的伯德图，计算系统的幅值裕度和相位裕度。2、前向通路中插入一相位滞后校正，确定校正网络的传递函数。3、用MATLAB画出未校正和已校正系统的根轨迹。4、用MATLAB对校正前后的系统进行仿真分析，画出阶跃响应曲线，计算其时域性能指标。5、课程设计说明书中要求写清楚计算分析的过程，列出MATLAB程序和MATLAB输出。说明书的格式按照教务处标准书写。时间安排:任务时间（天）审题、查阅相关资料2分析、计算2编写程序1撰写报告2论文答辩1指导教师签名：年月日系主任（或责任教师）签名：年月日)01)(1()(sssKsG目录摘要-----------------------------------------------11设计目的、要求及原理------------------------------21.1设计目的--------------------------------------21.2设计要求-------------------------------------21.3设计原理--------------------------------------22设计分析与计算------------------------------------32.1最小K值的系统频域分析------------------------32.2滞后校正函数计算------------------------------43用MATLAB画校正前后的轨迹-------------------------83.1校正前的根轨迹--------------------------------83.2校正后的根轨迹--------------------------------94用Matlab对校正前后的系统进行仿真分析------------114.1校正前系统-----------------------------------114.2校正后系统-----------------------------------13心得体会------------------------------------------15参考文献------------------------------------------164摘要自动控制技术已广泛应用于制造业、农业、交通、航空及航天等众多产业部门，极大地提高了社会劳动生产率，改善了人们的劳动条件，丰富和提高了人民的生活水平。在今天的社会生活中，自动化装置无所不在，为人类文明进步做出了重要贡献。在控制技术需求推动下，控制理论本身也取得了显著进步。从线性近似到非线性系统的研究取得了新的成就，借助微分几何的固有非线性框架来研究非线性系统的控制，已成为目前重要研究方向之一。为了实现各种复杂的控制任务首先要将被控制对象和控制装置按照一定的方式连接起来，组成一个有机整体，这就是自动控制系统。关键词:自动控制技术、系统分析、MATLAB、校正51设计目的、要求及原理1.1设计目的滞后校正网络具有低通滤波器的特性，因而当它与系统的不可变部分串联相连时，会使系统开环频率特性的中频和高频段增益降低和截止频率Ｗc减小，从而有可能使系统获得足够大的相位裕度，它不影响频率特性的低频段。由此可见，滞后校正在一定的条件下，也能使系统同时满足动态和静态的要求。本设计通过增加一个滞后校正装置，确定其最适合参数来改变系统性能。1.2设计要求要求系统的静态速度误差系数151sKv，相位裕度56。1.3设计原理通过分析系统逐步确定其矫正系统参数，具体步骤如下：1.根据稳态误差要求求出K值；2.画出未校正系统的波特图，并求；3.波特图上绘制出曲线；4.根据稳态误差要求，求出校正系统的截止频率；5.根据公式和,可求出b和t;6.验证已校正系统的相位裕度和幅值裕度；62设计分析与计算2.1最小K值的系统频域分析已知转子绕线机控制系统的开环传递函数是：)01)(1()(sssKsG，静态速度误差系数151sKv，15110/)(0limsKssGsvk所以最小的K值为：K=150故)01)(1(501)(ssssG1求相角裕度：因为1001501011501)(22sssA在穿越频率处)(A=1，解得Wc≈3.69rad/s穿越频率处的相角为：85.111.090)(11ctgtgcc相角裕度为：γ=180+)(c=180-185.1=-5.1deg2求幅值裕度：先求相角穿越频率:1801.090)(11gggtgtg901.011ggtgtg由三角函数关系得:3.16,11.0ggg解得：1.361001501)(22ggggA7所以，幅值裕度为:)(-2.69)(log20dBALgg使用MATLAB软件可直接得到系统的BODE图和相角,幅值裕度。程序的代码如下：n=150d=[1,20,15,0]g1=tf(n,d)[mag,phase,w]=bode(gf)margin(gf)图2-1校正前系统的BODE图Matlab仿真结果为：Gm=6.02dBPm=18.3deg（与理论计算结果相同）82.2滞后校正函数计算求滞后校正的网络函数可以按设计原理所讲述的方法进行求解，但过程比较麻烦，这里介绍使用matlab进行编程求解的方法。操作简单，可快速得到结果。由于按设计要求幅值裕度56。根据要求令相角裕度γ=65并附加6，即取γ=71。设滞后校正器的传递函数为：校正前的开环传递函数为：)01)(1(501)(ssssG用matlab编写滞后校正的程序代码如下：k0=150;n1=1;d1=conv(conv([10],[11]),[110]);Go=tf(k0*n1,d1);[mag,phase,w]=bode(Go);Mag=20*log10(mag);Pm=65;Pm1=Pm+6;Qm=Pm1*pi/180;b=(1-sin(Qm))/(1+sin(Qm));Lcdb=-20*log10(b);wc=spline(Mag,w,Lcdb);T=10/(wc*b);Tz=b*T;Gc=tf([Tz1],[T1])11)(TsbTssGc9图2-2滞后校正求解图得到结果为：使用matlab检验是否符合要求，程序代码为：K=150;n1=1;d1=conv(conv([10],[11]),[110]);s1=tf(K*n1,d1);n2=[25.581];d2=[913.51];s2=tf(n2,d2);sys=s1*s2;[mag,phase,w]=bode(sys);margin(sys)1913.515.582)(sssGc10图2-3校正后系统BODE图Matlab仿真结果为：Gm=32.8dBPm=77.7deg（符合设计要求）113用MATLAB画校正前后的轨迹3.1校正前的根轨迹由于系统未校正前的开环传递函数为：)01)(1(501)(ssssG使用MATLAB画根轨迹代码如下：num=150den=conv(conv([1,0],[1,1]),[1,10])rlocus(num,den)title('控制系统未校正前根轨迹图')图3-1校正前系统根轨迹图123.2校正后的根轨迹系统校正后的开环传递函数为：)1913.5)(01)(1()15.582(501)(ssssssG使用MATLAB画根轨迹代码如下：num=150*[25.58,1]den=conv(conv([1,10],[913.5,1]),[1,1,0])rlocus(num,den)title('控制系统校正后根轨迹图')图3-2校正后的系统根轨迹图13图3-3右方根轨迹放大图144用Matlab对校正前后的系统进行仿真分析4.1校正前系统系统未校正前的开环传递函数为：)01)(1(501)(ssssG单位负反馈闭环传递函数为：5011011501)()()(23ssssRsCs使用MATLAB求校正前系统单位阶跃响应的性能指标代码如下：num=150den=[1,11,10,150]t=0:0.01:20step(num,den,t)[y,x,t]=step(num,den,t)maxy=max(y)yss=y(length(t))pos=100*(maxy-yss)/yssfori=1:2001ify(i)==maxyn=i;endendtp=(n-1)*0.0115y1=1.05*yssy2=0.95*yssi=2001whilei0i=i-1ify(i)=y1|y(i)=y2;m=i;breakendendts=(m-1)*0.01title('stepresponse')Grid图4-1校正前单位阶跃响应曲线图164.2校正后系统系统校正后的开环传递函数为：)1913.5)(01)(1()15.582(501)(ssssssG单位负反馈闭环传递函数为5013847914610049.5913.5)15.582(501)()()(234ssssssRsCs使用MATLAB求校正后系统单位阶跃响应的性能指标代，代码如下：num=150*[25.58,1]den=[913.5,10049.5,9146,3847,150]s1=tf(K*n1,d1);Lsys=tf(num,den);[y,t,x]=step(Lsys);plot(t,y)然后使用matlab中的LTIViewer工具，在MATLAB提示符后,输入ltiview，即可启动该图形软件。从File的下拉菜单中选中→import选项选择需要仿真的系统。选择窗口中的Lsys系统，并用鼠标点击OK。在画面中点击鼠标右键，选择“Characteristics”选项，再选择“PeakTime”项可得阶跃响应曲线中的峰值时间。在画面中点击鼠标右键，选择“Characteristics”选项，再选择“SettlingTime”、“RiseTime”、“SteadyState”选项可得阶跃响应曲线中的调节时间，上升时间,稳态值，单击各个点，图中显示相应的数据，可通过鼠标移动标注的位置。通过点击“Edit”菜单，在弹出的下拉菜单中选择“ViewerPreferences”项，设定阶跃响应的上升时间范围为最终稳态值的0~90%，调节时间的误差带为2%，即可得到系统校正后的单位阶跃响应曲线。17课程设计答辩或质疑记录：1如何求相角裕度？答：先求c180即可。2如何用MATLAB作图？答：在页面中输入程序，点击保存并运行即可。3对此次课设有什么想法或建议？答：课设是很有意义的事情，在自己的努力下完成一项工作，会有成就感，激发学生兴趣，而且课设可以考验学生查找资料，运用计算机和学习协作的能力，是一次难得的锻炼。成绩评定依据：评定项目最高分限评分成绩1．选题合理、目的明确102．设计方案正确、具有可行性、创新性203．设计结果（例如：系统设计程序、仿真程序）204．态度认真、学习刻苦、遵守纪律155．设计报告的规范化、参考文献充分（不少于5篇）106．答辩25总分100最终评定成绩（以优、良、中、及格、不及格评定）指导教师签字：年月日