1、课程设计(大作业)报告课程名称: 自动控制理论 设计题目: 自动控制系统建模、分析及校正院系:自动控制与机械工程学院 班级:电气工程及其自动化 2014 级 2 班设 计 者: 学号: 指导教师: 设计时间:2016.12.122016.12.16昆明学院课程设计(大作业)任务书姓名:院(系):自动控制与机械工程学院专业:电气工程及其自动化学号:201404170146任务起止日期:第十五周(2016 年 12 月 12 日2016 年 12 月 16 日)课程设计题目:自动控制系统建模、分析及校正课程设计要求:1. 了解 matlab 软件的基本特点和功能,熟悉其界面、菜单和工具条;掌握线性
2、系统模型的计算机表示方法、变换以及模型间的相互转换。了解控制系统工具箱的组成、特点及应用;掌握求线性定常连续系统输出响应的方法,运用连续系统时域响应函数(impulse,step,lsim),得到系统的时域响应曲线。2. 掌握使用 MATLAB 软件绘制开环系统的幅相特性曲线、对数频率特性曲线;观察控制系统的观察开环频率特性,对控制系统的开环频率特性进行分析;3. 掌握 MATLAB 软件中 simulink 工具箱的使用;熟悉 simulink 中的功能模块,学会使用 simulink 对系统进行建模;掌握 simulink 的仿真方法。工作计划及安排:在第 15 周内(2016.12.12
3、2016.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 是用来对真实世界的动力学系统建模、
4、模拟仿真和分析的软件。我们可将综合性和设计性实验项目通过 MATLAB 在计算机上仿真,使系统的观察实验的动态过程。目前,MATLAB已经成为我们当代大学生必须掌握的基本技能,在设计研究单位和工业部门,MATLAB已经成为研究和解决各种具体工程问题的一种标准软件。在完成了验证性、综合性和设计性实验后,课程设计必不可少。课程设计是工科实践教学的一个重要的环节,目的是培养我们综合运用理论知识分析和解决实际问题的方法和能力,实现由知识向技能的初步化。所以课程设计是培养我们思维创造能力最有效的途径。指导教师评语:成绩:填表时间:指导教师签名:课程设计(大作业)报告1. 用 matlab 语言编制程序,
5、实现以下系统:5s3 + 24s 2 +181) G(s) =s 4 + 4s3 + 6s 2 + 2s + 24(s + 2)(s 2 + 6s + 6)22) G(s) =s(s +1)3 (s3 + 3s 2 + 2s + 5)num = 5, 24, 0,18 den = 1, 4, 6, 2, 2 sys = tf (num, den)Transferfunction5s3 + 24s2 + 18s4 + 4s3 + 6s2 + 2s+2num = 4 * conv (1, 2, conv (1, 6, 6,1, 6, 6)den = conv (1, 0, conv (1,1, c
6、onv (1,1, conv (1,1,1, 3, 2,5)sys = tf (num, den)num =45628872720288den =1614212417504s5 + 56s4 + 288s3 + 672s2 + 720s + 288Transfer function: s7 + 6s6 + 14s5 + 21s4 + 24s3 + 17s2 + 5s2. 两环节 G1、G2 串联,求等效的整体传递函数 G(s)G1 (s) =2s + 3G2 (s) =7s 2 + 2s + 1g1 = tf (2,1, 3)g2 = tf (7,1, 2,1)g = g1* g2Transf
7、er function:14s3 + 5s2 + 7s + 33. 两环节 G1、G2 并联,求等效的整体传递函数 G(s)G1 (s) =2s + 3G2 (s) =7s2 + 2s +1g1 = tf (2,1, 3)g2 = tf (7,1, 2,1)g = g1 + g22s2 + 11s + 23Transfer function: s3 + 5s2 + 7s + 34. 已知系统结构如图,求闭环传递函数。其中的两环节1G1、G2 分别G为(s) =3s +100s2 + 2s + 81G2 (s) =22s + 5(1) 正反馈:g1 = tf (3, 100 , 1, 2, 81
8、)(2) 负反馈:g1 = tf (3, 100 , 1, 2, 81)g2 = tf (2 , 2, 5)gf = feedback(g1, g2, 1)g2 = tf (2 , 2, 5)gf = feedback(g1, g2)Transfer function:Transfer function:6s2+ 215s + 5006s2+ 215s + 5002s3+ 9s2+ 166s + 2052s3+ 9s2+ 178s + 6055. 已知某闭环系统的传递函数为G(s) =曲线,单位脉冲响应曲线。(1) 单位阶跃响应曲线:10s + 250.16s3 + 1.96s 2 + 10s
9、 + 251.4,求其单位阶跃响应Step Responsenum= 10, 25den = 0. 16, 1. 96, 10, 251.2step(num,den) grid1Amplitude0.80.60.40.2(2) 单位脉冲响应曲线:000.511.522.5Time (sec)num= 10, 25den = 0. 16, 1. 96, 10, 25 5Impulse Responseimpulse(num,den)grid43Amplitude210-100.511.522.5Time (sec)w26. 典型二阶系统的传递函数为G(s) = ns2 + 2xws +w2,wn
10、为自然频率,x 为阻尼比,nn试绘出当x=0.5,wn分别取-2、0、2、4、6、8、10 时该系统的单位阶跃响应曲线;分析阻尼比分别为0.5、1 时系统的稳定性。(1)当x=0.5,wnw=0:2:10;kosai=0.5; figure(1) hold on for Wn=w num=Wn2;分别取 2、4、6、8、10 时1.41.21Step Responseden=1,2*kosai*Wn,Wn2; step(num,den)endgrid on;0.8Amplitude0.60.40.2w=-2kosai=0.5; figure(1) Wn=w num=Wn2;den=1,2*ko
11、sai*Wn,Wn2; step(num,den)grid on; title() xlabel() ylabel()w =-2Wn = -200123456Time (sec)x 1012位 位 位 位 位 位864 位20-2-4051015202530 位 (sec)(2)分析阻尼比分别为0.5、1 时系统的稳定性1、当x=-0.5 时:w=-2:2:10;kosai=-0.5; hold onfor Wn=w den=1,2*kosai*Wn,Wn2;v=roots(den) endv =-1.0000 + 1.7321i-1.0000 - 1.7321iv =00v=1.0000+1
12、.7321i1.0000-1.7321iv=2.0000+3.4641i2.0000-3.4641iv=3.0000+5.1962i3.0000-5.1962iv=4.0000+6.9282i4.0000-6.9282iv=5.0000+8.6603i5.0000-8.6603i由特征根可得,当 w=-2 时系统稳定,当 w=0 时系统临界稳定,当 w=2、4、6、8、10时,系统不稳定。2、当x=-1 时:clearw=-2:2:10;kosai=-1; hold on for Wn=wden=1,2*kosai*Wn,Wn2; v=roots(den)endv = -2-2v =00v =22v =44v =6.0000 + 0.0000i 6.0000 - 0.0000iv =88v = 1010由特征根可得,当 w=-2 时系统稳定,当 w=0 时系统临界稳定,当 w=2、4、6、8、10时,系统不稳定。7. 设有一高阶系统开环传递函数为G(s) =0.016s3 + 0.218s2 +1.436s + 9.3590.06s3 +