自动控制理论课程设计-自动控制系统建模、分析及校正.doc

课程设计（大作业）报告课程名称： 自动控制理论 设计题目： 自动控制系统建模、分析及校正 院 系：自动控制与机械工程学院 班 级：电气工程及其自动化2014级 2班 设 计 者： 学 号： 指导教师： 设计时间：2016.12.122016.12.16 昆明学院课程设计（大作业）任务书姓 名：院（系）：自动控制与机械工程学院专 业：电气工程及其自动化 学 号：201404170146任务起止日期：第十五周（2016年12月12日2016年12月16日）课程设计题目：自动控制系统建模、分析及校正课程设计要求：1了解matlab软件的基本特点和功能,熟悉其界面、菜单和工具条；掌握线性系统模型的计算机表示方法、变换以及模型间的相互转换。了解控制系统工具箱的组成、特点及应用；掌握求线性定常连续系统输出响应的方法，运用连续系统时域响应函数（impulse,step,lsim），得到系统的时域响应曲线。2掌握使用MATLAB软件绘制开环系统的幅相特性曲线、对数频率特性曲线；观察控制系统的观察开环频率特性，对控制系统的开环频率特性进行分析；3掌握MATLAB软件中simulink工具箱的使用；熟悉simulink中的功能模块，学会使用simulink对系统进行建模；掌握simulink的仿真方法。工作计划及安排：在第15周内（2016.12.122016.12.16）内完成规定的题目。1：2016.12.12明确注意事项，完成题目1、2。2：2016.12.13完成题目3、4、5。3：2016.12.14完成题目6、7、8。4：2016.12.15完成题目9、10。5：2016.12.16完成题目11，打印并提交课程设计报告。指导教师签字 2016年 12 月 8 日 课程设计（大作业）成绩学号：201404170146 姓名：张雨 指导教师：杨祖元课程设计题目：自动控制系统建模、分析及校正总结：MATLAB 是一种直观、高效的计算机语言，同时也是一个科学计算平台。它的伴随工具Simulink 是用来对真实世界的动力学系统建模、模拟仿真和分析的软件。我们可将综合性和设计性实验项目通过MATLAB 在计算机上仿真，使系统的观察实验的动态过程。目前，MATLAB 已经成为我们当代大学生必须掌握的基本技能，在设计研究单位和工业部门，MATLAB 已经成为研究和解决各种具体工程问题的一种标准软件。在完成了验证性、综合性和设计性实验后，课程设计必不可少。课程设计是工科实践教学的一个重要的环节，目的是培养我们综合运用理论知识分析和解决实际问题的方法和能力,实现由知识向技能的初步化。所以课程设计是培养我们思维创造能力最有效的途径。指导教师评语：成绩：填表时间：指导教师签名：课程设计（大作业）报告1.用matlab语言编制程序，实现以下系统：num =4 56 288 72 720 288den =1 6 14 21 24 17 5 0Transfer function:Transfer function1） 2）2.两环节G1、G2串联，求等效的整体传递函数G(s) Transfer function: 3.两环节G1、G2并联，求等效的整体传递函数G(s) Transfer function:4.已知系统结构如图，求闭环传递函数。其中的两环节G1、G2分别为（2）负反馈： Transfer function: （1）正反馈： Transfer function: 5.已知某闭环系统的传递函数为，求其单位阶跃响应曲线，单位脉冲响应曲线。（1）单位阶跃响应曲线：（2）单位脉冲响应曲线：6.典型二阶系统的传递函数为 ， 为自然频率， 为阻尼比，试绘出当=0.5，分别取-2、0、2、4、6、8、10时该系统的单位阶跃响应曲线；分析阻尼比分别为0.5、1时系统的稳定性。（1）当=0.5，分别取2、4、6、8、10时w=0:2:10;kosai=0.5;figure(1)hold onfor Wn=wnum=Wn2; den=1,2*kosai*Wn,Wn2; step(num,den) endgrid on;w=-2kosai=0.5;figure(1)Wn=wnum=Wn2;den=1,2*kosai*Wn,Wn2;step(num,den)grid on;title()xlabel()ylabel()w =-2Wn = -2（2）分析阻尼比分别为0.5、1时系统的稳定性1、当=-0.5时：w=-2:2:10; kosai=-0.5;hold onfor Wn=wden=1,2*kosai*Wn,Wn2;v=roots(den)endv = -1.0000 + 1.7321i -1.0000 - 1.7321iv = 0 0v = 1.0000 + 1.7321i 1.0000 - 1.7321iv = 2.0000 + 3.4641i 2.0000 - 3.4641iv = 3.0000 + 5.1962i 3.0000 - 5.1962iv = 4.0000 + 6.9282i 4.0000 - 6.9282iv = 5.0000 + 8.6603i5.0000 - 8.6603i由特征根可得，当w=-2时系统稳定，当w=0时系统临界稳定，当w=2、4、6、8、10时，系统不稳定。2、当=-1时：clearw=-2:2:10;kosai=-1;hold onfor Wn=wden=1,2*kosai*Wn,Wn2;v=roots(den)endv = -2 -2v = 0 0v = 2 2v = 4 4v =6.0000 + 0.0000i 6.0000 - 0.0000iv = 8 8v = 10 10由特征根可得，当w=-2时系统稳定，当w=0时系统临界稳定，当w=2、4、6、8、10时，系统不稳定。7. 设有一高阶系统开环传递函数为，试绘制该系统的零极点图和闭环根轨迹图。（1）系统的零极点：num=0.016 0.218 1.463 9.359;den=0.06 0.268 0.635 6.271;z,p,k=tf2zp(num,den)pzmap(num,den)运行结果：z = -10.2679 -1.6785 + 7.3587i -1.6785 - 7.3587ip = -5.7710 0.6522 + 4.2054i 0.6522 - 4.2054ik = 0.2667（2）系统的闭环根轨迹：num=0.016 0.218 1.463 9.359;den=0.06 0.268 0.635 6.271;rlocus(num,den)8.单位反馈系统前向通道的传递函数为： ，试绘制该系统的Bode图和Nyquist曲线，说明软件绘制曲线与手动绘制曲线的异同。（1）bode图num=0,0,2,8,12,8,2;den=1,5,10,10,5,1,0;bode(num,den)（2）Nyquist曲线num=0,0,2,8,12,8,2;den=1,5,10,10,5,1,0;nyquist(num,den)9已知某控制系统的开环传递函数，试绘制系统的开环频率特性曲线，并求出系统的幅值与相位裕量。G=tf(1.5,conv(conv(1,0,1,1),1,2);bode(G) gridGm,Pm,Wcg,Wcp=margin(G)Gm =4.0000Pm =41.5340Wcg =1.4142Wcp =0.611810.在SIMULINK中建立系统，该系统阶跃输入时的连接示意图如下。k为学生学号后三位。绘制其单位阶跃响应曲线，分析其峰值时间tp、延迟时间td、上升时间tr、调节时间ts及超调量。num1=1;den1=1,0;sys1=tf(num1,den1);num2=146; den2=1,9; sys2=tf(num2,den2); sys12=sys1*sys2; G=feedback(sys12,1);step(G)grid分析其峰值时间tp、延迟时间td、上升时间tr、调节时间ts及超调量,闭环系统的函数：（1） 峰值时间：G1=tf(146,1,9,146);y,t=step(G1);Y,K=max(y);tp=t(K)tp =0.2798（2） 延迟时间y,t=step(G1);c=dcgain(G); i=1;while y(i)0.5*ci=i+1;endtd=t(i)td =0.1049（3） 上升时间G1=tf(146,1,9,146);y,t=step(G1);c=dcgain(G); n=1; while y(n)0.98*c)&(y(i)1.02*c) i=i-1; end ts=t(i) ts = 0.8743（5） 超调量G1=tf(146,1,9,146);y,t=step(G1);c=dcgain(G1);Y,K=max(y);a=(Y-c)/ca = 0.2835*11. 给定系统如下图所示，试设计一个串联校正装置，使截止频率c=40rad/s、相位裕度45G=tf(100,0.04,1,0);Gw,Pw,Wcg,Wcp=margin(G)G1=tf(100,0.04,1,0);G2=tf(100*0.025,1,conv(0.04,1,0,0.01,1)bode(G1)holdbo