自动控制原理课程设计报告.doc

大学自动控制原理课程设计报告课程名称：自动控制原理 设计题目：自动控制原理MATLAB仿真 院 系： 班 级： 设 计 者： 学 号： 指导教师： 设计时间：2015年12月 目 录一.实验目的和意义：3二简述MATLAB语言的特点及其主要功能：3三实训内容41.用matlab语言编制程序，实现以下系统：42.两环节G1，G2串联，求等效的整体传递函数G(s)：53.两环节G1，,G2并联，求等效的整体传递函数G(s)：54已知系统结构如图，求闭环传递函数:55. 求闭环系统的单位阶跃响应曲线，单位脉冲响应曲线：66.典型二阶系统的为自然频率，为阻尼比，试绘制出当=0.5，分别取-2 ,0 ,2 ,4 ,6 ,8 ,10时该系统的单位阶跃响应曲线；分析阻尼比分别为-0.5，-1时系统的稳定性。77试绘制高阶系统开环传递函数的零极点图和闭环根轨迹图：98试绘制单位反馈系统的Bode图和Nyquist曲线：119试绘制开环传递函数系统的开环频率特性曲线，并求出系统的幅值和相位裕量：1310. 在SIMULINK中建立系统，绘制其单位阶跃响应曲线，分析其峰值时间，延迟时间，上升时间，调节时间及超调量：14*11. 给定系统如下图所示，试设计一个串联校正装置，使幅值裕度大于h10分贝、相位裕度45。15四、心得体会19五.参考文献19一.实验目的和意义：1了解matlab软件的基本特点和功能，熟悉其界面、菜单和工具条；掌握线性系统模型的计算机表示方法、变换以及模型间的相互转换。了解控制系统工具箱的组成、特点及应用；掌握求线性定常连续系统输出响应的方法，运用连续系统时域响应函数（impulse，step，lsim），得到系统的时域响应曲线。2掌握使用MATLAB软件作出系统根轨迹；利用根轨迹图对控制系统进行分析；掌握使用MATLAB软件作出开环系统的波特图，奈奎斯特图；观察控制系统的观察开环频率特性，对控制系统的开环频率特性进行分析；3掌握MATLAB软件中simulink工具箱的使用；熟悉simulink中的功能模块，学会使用simulink对系统进行建模；掌握simulink的仿真方法。二简述MATLAB语言的特点及其主要功能：MATLAB是矩阵实验室（Matrix Laboratory）的简称，是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。 MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连 matlab开发工作界面接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。实训内容1.用matlab语言编制程序，实现以下系统：1）G(s)= 程序:num=5 24 0 18; den=1 4 6 2 2;h=tf(num,den)Transfer function: 5 s3 + 24 s2 + 18-s4 + 4 s3 + 6 s2 + 2 s + 22) G(s)= 输入以下程序： n1=4*1 2;n2=1 6 6;n3=1 6 6;num=conv(n1,conv(n2,n3);d1=1 1;d2=1 1;d3=1 1;d4=1 3 2 5;den1=conv(d1,d2);den2=conv(d3,d4);den=den1 den2 0;h=tf(num,den)Transfer function: 4 s5 + 56 s4 + 288 s3 + 672 s2 + 720 s + 288-s8 + 2 s7 + s6 + s5 + 4 s4 + 5 s3 + 7 s2 + 5 s2.两环节G1，G2串联，求等效的整体传递函数G(s)： G1(s)= G2(s)=程序： n1=2;d1=1 3;sys1=tf(n1,d1); n2=7;d2=1 2 1;sys2=tf(n2,d2); sys12=sys1*sys2Transfer function: 14-s3 + 5 s2 + 7 s + 33.两环节G1，,G2并联，求等效的整体传递函数G(s)： G1(s)= G2(s)=输入以下指令： num1=2;den1=1 3;sys1=tf(num1,den1); num2=7;den2=1 2 1;sys2=tf(num2,den2); sys12=sys1+sys2Transfer function: 2 s2 + 11 s + 23-s3 + 5 s2 + 7 s + 34已知系统结构如图，求闭环传递函数:已知系统结构如图，求闭环传递函数。其中的两环节G1,G2分别为 G1(s)= G2(s)=输入以下指令： n1=3 100;d1=1 2 81;n2=2;d2=2 5;s1=tf(n1,d1);s2=tf(n2,d2);sys=feedback(s1,s2)Transfer function: 6 s2 + 215 s + 500-2 s3 + 9 s2 + 178 s + 6055. 求闭环系统的单位阶跃响应曲线，单位脉冲响应曲线： 已知某闭环系统的传递函数为G(s)=,求其单位阶跃响应曲线，单位脉冲响应曲线。1） 单位阶跃响应：输入如下指令：G=tf(10,25,0.16,1.96,10,25)step(G); %计算并绘制系统的单位阶跃响应title;（单位阶跃响应曲线图） 图05.01单位阶跃响应曲线 （2）、单位脉冲响应： 【程序】 G=tf(10,25,0.16,1.96,10,25)impulse(G); %计算并绘制系统的单位脉冲响应title(单位脉冲响应)图05.02单位脉冲响应6.典型二阶系统的为自然频率，为阻尼比，试绘制出当=0.5，分别取-2 ,0 ,2 ,4 ,6 ,8 ,10时该系统的单位阶跃响应曲线；分析阻尼比分别为-0.5，-1时系统的稳定性。典型二阶系统的传递函数为G(s)=,为自然频率，为阻尼比，试绘制出当=0.5，分别取-2 ,0 ,2 ,4 ,6 ,8 ,10时该系统的单位阶跃响应曲线；分析阻尼比分别为-0.5，-1时系统的稳定性。输入如下指令：（1）.当=0.5，分别取-2、0、2、4、6、8、10时：G0=tf(0,4,1,-2,4); % =0.5 Wn=-2G1=tf(0,4,1,2,4); % =0.5 Wn=2G2=tf(0,16,1,4,16); % =0.5 Wn=4G3=tf(0,36,1,6,36); % =0.5 Wn=6G4=tf(0,64,1,8,64); % =0.5 Wn=8G5=tf(0,100,1,10,100); % =0.5 Wn=10step(G0,G1,G2,G3,G4,G5); title(图06.01单位阶跃响应曲线G6=tf(0,100,1,-10,100);%=-0.5,Wn=10时的单位阶跃响应G7=tf(0,100,1,-20,100);%=-1,Wn=10时的单位阶跃响应step(G6,G7);title(=-0.5，=-1时的阶跃响应)图06.02阶跃响应曲线7试绘制高阶系统开环传递函数的零极点图和闭环根轨迹图：设有一高阶系统开环传递函数为G(s)=,试绘制该系统的零极点图和闭环根轨迹图。 1）系统的零极点 输入如下指令： G=tf(0.016,0.218,1.436,9.359,0.06,0.268,0.635,6.271) pzmap(G);title()图07.01系统的零极点2）系统的闭环根轨迹 输入如下指令： num=0.016 0.218 1.436 9.359; den=0.06 0.268 0.635 6.271; rlocus(num,den)图07.02系统的闭环根轨迹8试绘制单位反馈系统的Bode图和Nyquist曲线：单位反馈系统前向通道的传递函数为：G(s)=,试绘制该系统的Bode图和Nyquist曲线，说明软件绘制曲线与手动绘制曲线的异同。 1）绘制该系统的Bode图： 输入如下指令： num=0 0 2 8 12 8 2; den=1 5 10 10 5 1 0; margin(num,den)图08.01 系统的Bode图 2)系统的Nyquist图： 输入如下指令： num=0 0 2 8 12 8 2; den=1 5 10 10 5 1 0; nyquist(num,den)图08.02 系统Nyquist曲线9试绘制开环传递函数系统的开环频率特性曲线，并求出系统的幅值和相位裕量：已知某控制系统的开环传递函数G(s)=,=1.5,试绘制系统的开环频率特性曲线，并求出系统的幅值和相位裕量。 输入如下指令: d1=1 1; d2=1 2 den1=conv(d1,d2); den=den1 0; num=10; bode(num,den) 程序运行结果如图所示：图09.01系统的开环频率特性曲线10. 在SIMULINK中建立系统，绘制其单位阶跃响应曲线，分析其峰值时间，延迟时间，上升时间，调节时间及超调量：在SIMULINK中建立系统，该系统阶跃输入时的连接示意图如下。K为学生学号后三位。绘制其单位阶跃响应曲线，分析其峰值时间，延迟时间，上升时间，调节时间及超调量。图10.01单位阶跃响应曲线输入如下程序： num=238; den=1 9 238; step(num,den); grid;图10.02示波器显示单位阶跃响应曲线延迟时间：=0.2757上升时间:tr=0.094峰值时间:tp=0.392调节时间:ts=1.192超调量:36.24%*11. 给定系统如下图所示，试设计一个串联校正装置，使幅值裕度大于h10分贝、相位裕度45。（1）、原系统的幅值裕度和相位裕度图11.01原系统传函图校正后的图图11.02校正系统传函图G=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)holdbode(G2,r)gridfigureG1c=feedback(G1,1);G2c=feedback(G2,1);step(G1c)holdstep(G2c,r)grid结果如下：Gw = InfPw = 28.0243Wcg = InfWcp = 46.9701Transfer function: 2.5 s + 100-0.0004 s3 + 0.05 s2 + s图11.03校正前系统频率特性曲线图11.04校正后系统频率特性曲线四、心得体会为期一周的自动控制原理课程设计，不仅是对我们在课本上学到的知识的一种检验，而且是对自己动手能力得一种提升。设计过程中我们遇到了很大的困难比如原先学过的知识有些遗忘以及自己对课本知