《自动控制原理》课程设计.docx

1、-贮 一. f 一心一卢 - 勹.-_,-自动控制原理课程设计一2-,.-,lll-. _- 一-=-一 作者： 日期：个人收集整理 勿做商业用途名称： 自动控制原理课程设计题目:基于自动控制原理的性能分析设计与校正院系：建筑环境与能源工程系班级： 学生姓名：指导教师：目录一、课程设计的目的与要求-3二、设计内容2。1 控制系统的数学建模-42。2 控制系统的时域分析-6 2.3 控制系统的根轨迹分析-3个人收集整理 勿做商业用途-82。4 控制系统的频域分析-10 2.5 控制系统的校正-12 三、课程设计总结-17 四、参考文献-18一、 课程设计的目的与要求4个人收集整理 勿做商业用途本

2、课程为自动控制原理的课程设计，是课堂的深化。设置自动控制原理课程设计的目的是使 MATLAB 成为学生的基本技能， 熟悉 MATLAB 这一解决具体工程问题的标准软件，能熟练地应用 MATLAB 软件解决控制理论中的复杂和工程实际问题,并给以后的模糊控制理论、最优控制理论和多变量控制理论等奠定基础.使相关专业的本科学生学会应用这一强大的工具，并掌握利用 MATLAB 对控制理论内容进行分析和研究的技能，以达到加深对课堂上所讲内容理解的目的。通过使用这一软件工具把学生从繁琐枯燥的计算负担中解脱出来，而把更多的精力用到思考本质问题和研究解决实际生产问题上去。通过此次计算机辅助设计，学生应达到以下的

3、基本要求： 1.能用 MATLAB 软件分析复杂和实际的控制系统。2。能用 MATLAB 软件设计控制系统以满足具体的性能指标要求。3。能灵活应用 MATLAB 的 CONTROL SYSTEM 工具箱和 SIMULINK 仿真软件，分析系统的性能。5个人收集整理 勿做商业用途二、设计内容2。1 控制系统的数学建模控制系统的分析是以控制系统的数学模型为基础的。数学模型的定义:数学模型是描述系统动态特性及其变量之间关系的数学表达式或其他形式的表达。描述系统变量的各阶导数之间关系的微分方程称为系统的动态模型。在静态条件（描述系统变量的各阶导数为零）下，描述变量之间关系的代数方程称为静态模型。数学模

4、型的特点:相似化和抽象化，尽管组成系统模型参数的物理含义各不相同，但它们数学模型的形式很可能是相同的，从数学观点来看,只要数学模型是相同的,那么它们就应该有相同的运动规律，而不论它们的具体参数含义是什么，具有相同数学模型的不同的具体系数称为相似系统。简化性和精确性， 在建模的时候，要再简化和精确之间作折衷选择，其原则是简化后的数学方程的解的结果必须满足工程实际的要求并留有一定的余地.数学模型的种类:数学模型有多种形式，究竟选用哪一种模型，一般要视采用的分析方法和系统的类型而定，比如：连续系统的单输入/单输出系统的时域分析法可采用微分方程，连续多输入多输出系统的时域分析法可以采用状态方程，离散系

5、统可以采用差分方程等。常用的数学模型有微分方程、传递函数、差分方程、状态方程、结构图、频率特性等。在 MATLAB 中，常用的控制系统数学模型主要包括 TF 模型(多项式模型）、ZPK 模型（零极点模型）和 SS 模型（状态空间模型)。s 2 + 2s + 3例 104 已知系统的传递函数为： G(s) = s3 + 4s 2 + 6s + 9 ,在MATLAB 环境下 获 得 其 连 续 传 递 函 数 形 式 模 型 。 已 知 系 统 的 脉 冲 传 递 函 数z 2 + 2z + 3为: G(s) = z 3 + 4z 2 + 6z + 9 ，在MATLAB 环境下获得其采样时间为 4

6、 秒的传递函数形式模型. 解： num=1 2 3；描述分子多项式系数，按照降幂排列,且最后一位是 0 次幂的系数,中间用空格隔开.den=1 4 6 9;描述分母多项式系数，按降幂，最后一位是 0 次幂。G1=tf（num，den）%定义 G，tf（）表示建立多项式传递函数。6个人收集整理 勿做商业用途运行结果：Transfer function:s2 + 2 s + 3-s3 + 4 s2 + 6 s + 9再输入命令：zpk(G）%得到系统的零极点模型（zpk 模型）:Zero/pole/gain：（s2+ 2s + 3)-（s+3) （s2+ s + 3）7个人收集整理 勿做商业用途2

7、。2 控制系统的时域分析时域分析法就是根据输入、输出微分方程或传递函数数学模型，在时间域中分析控制系统的稳定性、稳态性能、动态性能。时域分析法是一种直接准确的分析方法，易为人么所接受，它可以接受系统时域内的全部信息。时域分析法包括稳定性分析、稳态性能分析（稳态误差）、动态性能分析三大方面。在 MATLAB 软件中稳定性能的分析可以直接求出特征根或用古尔维茨判据判定稳定性，而稳态误差的求取可根据静态误差系数,利用求极限的方法求取（与手算类似不再考虑），也可以输出中直接看出。第三方面动态性能主要是根据系统的各种响应来分析的，所以要学习一下在 MATLAB 软件中如何获取各种响应的命令函数。常用的响

8、应命令函数如下图所示:1.古尔维茨判据设线性定常系统特征方程：01n-1n0D(s) = a sn + a sn-1 + .+ as1 + a = 0（ a 0 ）将特征方程的系数按下列规则则 系统稳定的必要充分条件是：特征方程的各项系数均大于 0，且如下 hurwitz 行列式全部大于 0：aaa135.a2n-1D1aaa0a2a4. a2n-2=1= a3Da1123。.。 Dna= 0aa.a2n-302.000.an古尔维茨判据为：系统稳定的充分必要条件是古尔维茨行列式的各阶主子行列式均大于零，即 Di 0(i = 1,2,.,n) 。【例 1010】 系统闭环特征方程分别如下，试确

9、定特征根在 s 平面的位置，并判断系统闭环稳定性。（1）s4+2s3+3s2+4s+5=0 (2）s3+20s2+9s+100=0试用古尔维茨判据判别系统的稳定性。（1）:解： 运行程序：d1=2;一阶行列式的值a=2 4；1 3;%建立二阶矩阵 ad2=det（a）；%求矩阵 a 所对应的行列式的值b=2 4 0；1 3 5;0 2 4;%建立三阶矩阵 bd3=det（b）；%求矩阵 b 所对应的行列式的值c=2 4 0 0；1 3 5 0；0 2 4 0；0 1 3 5;%建立四阶矩阵 cd4=det(c)；%求矩阵 c 所对应的行列式的8个人收集整理 勿做商业用途值if(（d10）（d2

10、0）(d30）&（d40）)给出条件：如果 d1,d2，d3,d4同时大于 0，则出现如下结果：WARNDLG(The system is stable，Stability Analysis);else如果 d1，d2，d3，d4 不全大于 0，则出现如下对话框:WARNDLG（The system is unstable，Stability Analysis)； end运行结果如下：（2):解：运行程序如下：d=1 20 9 100;描述系统闭环特征方程的系数矩阵，系数中间用空格隔开。r=roots(d）%求解特征方程的根.运行结果： r =19.8005-0。0997 + 2。2451i-

11、0。0997 - 2.2451i由此可见:三个根都在 s 平面的左半部，则系统稳定，且其中一个位于实轴.9个人收集整理 勿做商业用途2。3 控制系统的根轨迹分析在控制系统分析中，为了避开直接求解高阶多项式的根时遇到的困难，在实践中提出了一种图解求根法，即根轨迹法。所谓根轨迹是指当系统的某一个（或几个）参数从到时，闭环特征方程的根在复平面上描绘的一些曲线。应用这些曲线，可以根据某个参数确定相应的特征根。在根轨迹法中，一般取系统的开环放大倍数 K 作为可变参数,利用它来反映出开环系统零极点与闭环系统极点(特征根)之间的关系。根轨迹可以分析系统参数和结构已定的系统的时域响应特性,以及参数变化对时域响

12、应特性的影响，而且还可以根据对时域响应特性的要求确定可变参数及调整开环系统零极点的位置，并改变它们的个数，也就是说根轨迹法可用于解决线性系统的分析与综合问题.10-26 (1)： 负反馈系统开环传递函数如下G (s) = k (s + 4)( s + 8)s 2 (s + 12) 2试绘制k 由 0+变化时其闭环系统的根轨迹并求出单位阶跃响应为衰减,等幅震荡,增幅震荡，单调增幅时的 K 值。解：运行程序:num=conv（1 4,1 8);1 4，1 8分别是分子两个相乘多项式的系数.多项式相乘，conv（多项式 1,多项式 2),最多两个多项式.den1=conv（1 12，1 12）;1

13、12,1 12分别是分母相乘多项式的系数。den2=conv(1 0,1 0）;% 同上。 den=conv（den1，den2)；%分母两个多项式相乘。rlocus(num,den)%rlocus：求系统根轨迹。运行结果：10个人收集整理 勿做商业用途10.80.60.40.20-0.2-0.4-0.6-0.8零极点图-1-12-10-8-6Real Axis-4-20Root Locus86420-2-4-6-8-12-10-8-6Real Axis-4-20因为根轨迹与虚轴交点是开环极点，对应的 K=0。所以当 K0 时响应是衰减，当 K=0 时响应是等幅震荡，没有增幅振荡和单调增幅。1

14、1个人收集整理 勿做商业用途2。4 控制系统的频域分析频域分析法是应用频率特性研究控制系统的一种经典方法。采用这种方法可直观地表达出系统的频率特性,分析方法比较简单，物理概念比较明确，对于诸如防止结构谐振，抑制噪声、改善系统稳定性和暂态性能等问题，都可以从系统的频率特性上明确的看出其物理实质和解决途径。频率分析法主要包括三种方法：Bode 图（幅频/相频特性曲线)、Nyquist 曲线、Nichols 图.1028 已知一振荡环节的传递函数为G (s) =1，求当 T=10， x=0.1、T 2 s 2 + 2xTs + 10.2、0。3、1.2 时的幅相频率特性曲线和对数幅频相频特性曲线.

15、解：运行程序： T=10;a=TT;num=1;for ks=0.1：0.1:1.2%ks（)变化范围是一个等差数列, 公差是 0。1，首末项分别为 0。1、1。2den=a 2*T*ks 1； figure（1）; nyquist(num，den）；hold on%输出图表 1045 奈奎斯特图（幅频相频特性曲线)。figure（2）；bode(num,den);hold on%输出图表 1046 伯德图（对数幅频相频特性曲线）end运行结果如下：得到其曲线分别为图 10-45 和图 1046 所示。12个人收集整理 勿做商业用途Nyquist Diagram_,l/二三/呵夕X/64)勹-

16、.一b-/ 子t/!夕2_l1/x-j/惑二(，, -、f- 一一今,/I, J - l-/夕0，_一-2./J夕夕-i.f-4-6-3-2-10Real Axis12310-45 系统的幅相频率特性曲线Bode Diagram200-20-40-60-800-45-90-135-180-3-2-1011010101010Frequency (rad/sec)131046 系统的对数幅频相频特性曲线个人收集整理 勿做商业用途2.5 控制系统的校正一个完整的自动控制系统的设计包括静态设计和动态设计两个部分，亦称系统的综合.静态设计包括选择执行元件、测量元件、比较元件和放大元件等，即把系统不可变部

17、分确定下来.而由不可变部分组成的控制系统往往不能满足性能指标的要求，甚至不能正常工作。动态设计则是根据性能指标的要求选择校正装置的形式和参数，使校正后系统的性能指标完全满足给定的性能指标的要求,即控制系统的校正.例 10-33已知单位负反馈系统的开环传递函数为：G(s) =10，试设计串联超前校正装置，使系统指标满足单位斜(0.1s +1)(0.001s +1)坡输入信号时稳态误差 ess0。1%，相位裕度g45，穿越频率c150 rad/s。解：G(s) =1000s(0.1s +1)(0.001s +1)根据稳态误差的要求做静态校正，则系统传递函数为：绘制 bode 图，求性能指标运行程序

18、：第一步输入程序：num=1000；描述分子多项式的系数。den=conv（0。1 1 0,0。001 1）; 描述分母多项式的系数。 margin（num,den）给定系统开环系统模型对象。运行结果:14个人收集整理 勿做商业用途100500-50-100-150-90-135-180-225Bode DiagramGm = 0.0864 dB (at 100 rad/sec) , Pm = 0.0584 deg (at 99.5 rad/sec)-270-101010121010345101010Frequency (rad/sec)10-52(a）动态校正前的系统 Bode 图从图 10

19、-52(a)可以看出相位裕度为 0.0584，穿越频率 99.5，都不满足要求，于是设计串联超前校正。第二步:串联超前校正的补偿角为:=45-0+7=52；校正参数 a 由超前网络最大超前角计算公式得： a=(1+sin(*pi/180）/（1-sin(pi/180）=8.43输入程序：syms t；%定义变量 tt=450+7；a=（1+sin(tpi/180)）/（1sin（tpi/180）) 运行结果如下：a =8。4344第三步:超前校正装置应该在系统的中频段，并给一定的裕度，所以取校正后的穿越频率即超前校正网络最大超前角频率为 160。令m =160;15个人收集整理 勿做商业用途可

20、以计算出超前校正装置时间常数 T输入程序： wm=160;T=1/（sqrt（a)wm) 运行结果如下:T =0。0022第四步输入程序：nc=a*T 1；dc= T 1；n=conv （num，nc)；d=conv(den，dc); Margin(n，d)margin：求幅值裕度和相角裕度及对应的转折频率运行结果：Bode DiagramGm = 16.8 dB (at 631 rad/sec) , Pm = 45.2 deg (at 175 rad/sec)100500-50-100-150-90-135-180-225-270-101010121010345101010Frequency

21、 (rad/sec)10-52（b）从图 10-52（b）可以看出相位裕度为 45。2,穿越频率为 175,都满足要求。第五步,验证系统的性能如下： 运行程序：t1=0:0。1:120;G1=tf(num，den); G11=feedback（G1,1）；step(G11,t1）;%step（）；求取系统单位阶跃16个人收集整理 勿做商业用途响应。G2=tf(n,d）;G22=feedback（G2，1）； figure(2); hold on;t2=0:0。1:8;step(G22，t2)运行结果：得到动态校正前后系统的阶跃响应曲线如下图21.81.61.41.210.80.60.4Step

22、 Response0.200204060Time (sec)80100120图 1053 动态校正前后系统的阶跃响应曲线17个人收集整理 勿做商业用途1.41.210.80.60.40.2Step Response001234Time (sec)567818个人收集整理 勿做商业用途三、课程设计总结学习的时间总是过的这么匆忙，两周的自动控制原理课程设计就这样不知不觉的过去了，好快，.在这两周，我从迷茫到想放弃，再到静下心来慢慢琢磨最后总算是还算圆满的完成的任务书上的内容。不想说累，因为充实的感觉来之不易。每次课程设计都是我们既期待又害怕的，期待的是大家可以比较自由的支配时间,而害怕的是，每个课

23、程设计我们总要好好的熬几个夜。还好最后的收获一般都是令人满意的,那种充实的感觉， 会让我们完全忘了经过的辛苦。生活就是这样，汗水预示着结果也见证着收获, 劳动是人类生存生活永恒不变的话题.我们不怕苦累，不畏麻烦与困难，才有机会品偿到丰收的喜悦。经过这两周的设计，我也深刻的体会到了团队的力量与合作的快乐，每个组员都为那个目标而努力琢磨，我们讨论,也争吵,但最终都是为了找到正确的答案。我不会忘了大家静静思索的皱眉，也不会忘了大家恍然明白时的幸福,你会感受到一群人为了同一个目标前进的力量是多么强大.团队需要个人，个人也离不开团队，必须发扬团结协作的精神。某个人的离群都可能导致导致整项工作的失败。设计

24、中只有一个人知道原理是远远不够的，必须让每个人都知道,否则一个人的错误，就有可能导致整个工作失败。团结协作是我们设计成功的一项非常重要的保证.我们需要彼此帮助，也需要彼此支撑,一个人的完成任务并不能算是这个小组的成功,每个人都完美的做出设计，才是我们每个人应该骄傲的事情. 这是最令我欣慰的地方,我们都不自私，我们都很团结，我们都在为其他人做最多的考虑，这是团队精神，也是同学情谊吧。我很希望我们不仅是在设计中，也在以后的生活学习中，都能时刻记着这样的情谊。学习的过程是辛苦的也是快乐的,我们乐在其中便也不会觉得累。谢谢我的组员给我的帮助以及包容，谢谢我的老师给我的指导以及关爱，每个人的成果都不能说

25、是他自己的收获，我们得到了太多人的帮助，从各方各面。在这里我只想表达对我的组员以及老师的感激，谢谢你们，在学习的道路上,有你们真好！19个人收集整理 勿做商业用途四、参考文献1.教材：自动控制原理,于希宁编著,中国电力出版社，2008 年. 2.基于 MATLAB 的系统分析与设计控制系统，楼顺天、于卫编著，西安电子科技大学出版社,1999 年.3。Matlab 工具箱应用指南-控制工程篇，徐昕、李涛、伯晓晨等编著，电子工业出版社，2000 年。4。反馈控制系统设计与分析MATLAB 语言应用,薛定宇 著,清华大学出版社，2000 年。5.自动控制原理,于希宁,孙建平著，中国电力出版社，2009 年。6.MATLAB在控制理论中的应用,陈晓平，李长杰，毛彦新著，中国科学技术大学出版社出版。7。自动控制原理习题精解，刘坤，刘翠响，李妍著,国防工业出版社,2005。2 （2）。8.MATLAB语言与控制系统仿真实训教程,刘振全，杨世凤著。20