控制系统数字仿真大作业(共19页).doc

《控制系统数字仿真大作业(共19页).doc》由会员分享，可在线阅读，更多相关《控制系统数字仿真大作业(共19页).doc（19页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上控制系统数字仿真课程 大作业 姓 名 ： 学 号 ： 班 级 ： 日 期 ： 同组人员 ： 专心-专注-专业目录一、引言MATLAB语言是当今国际控制界最为流行的控制系统计算机辅助设计语言，它的出现为控制系统的计算机辅助分析和设计带来了全新的手段。其中图形交互式的模型输入计算机仿真环境SIMULINK，为MATLAB应用的进一步推广起到了积极的推动作用。现在，MATLAB语言已经风靡全世界，成为控制系统CAD领域最普及、也是最受欢迎的软件环境。随着计算机技术的发展和应用，自动控制理论和技术在宇航、机器人控制、导弹制导及核动力等高新技术领域中的应用也愈来愈深入广泛。不仅

2、如此，自动控制技术的应用范围现在已发展到生物、医学、环境、经济管理和其它许多社会领域中，成为现代社会生活中不可或缺的一部分。随着时代进步和人们生活水平的提高，在人类探知未来，认识和改造自然，建设高度文明和发达社会的活动中，控制理论和技术必将进一步发挥更加重要的作用。作为一个自动化专业的学生，了解和掌握自动控制的有关知识是十分必要的。利用MATLAB软件及其SIMULINK仿真工具来实现对自动控制系统建模、分析与设计、仿真，能够直观、快速地分析系统的动态性能和稳态性能，并且能够灵活的改变系统的结构和参数，通过快速、直观的仿真达到系统的优化设计，以满足特定的设计指标。二、设计方法1、系统数学模型

3、美国卡耐尔基-梅隆大学机器人研究所开发研制了一套用于星际探索的系统，其目标机器人是一个六足步行机器人，如图(a)所示。该机器人单足控制系统结构图如图(b)所示。 要求： （1）建立系统数学模型； （2）绘制系统阶跃响应曲线、根轨迹图、频率特性； （3）分析系统的稳定性，及性能指标； （4）设计控制器Gc(s)，使系统指标满足：ts10s,ess=0,，超调量小于5%。 在不加入任何控制器的情况下，由控制系统结构图可得其开环传递函数为Simulink仿真得在单位阶跃信号作用下系统响应曲线为2、系统性能指标2.1 绘制系统阶跃响应曲线、根轨迹图、频率特性编写M文件如下：%set up system

4、 math modlenum=1;den=1 2 10 0;G0=tf(num,den)G=feedback(G0,1)figure(1);t=0:100;step(G,t);grid;title(单位阶跃响应曲线)y t=step(G,t);Y k=max(y);MATLAB命令窗口中运行后得到结果为：Transfer function: 1-s3 + 2 s2 + 10 s Transfer function: 1-s3 + 2 s2 + 10 s + 1ans = RiseTime: 21.5386 SettlingTime: 38.5480 SettlingMin: 0.9001 Se

5、ttlingMax: 0.9997 Overshoot: 0 Undershoot: 0 PeakTime: 79.1536tp=t(k)ess=1-y;figure(2);plot(t,ess);grid;title(单位阶跃响应误差曲线)figure(3);pzmap(G0);title(零极点分布图)figure(4);margin(G0);title(系统开环对数特性曲线)grid;figure(5);nyquist(G0);title(nyquist曲线)figure(6);rlocus(G0)title(系统根轨迹特性曲线)sgridstepinfo(G)以上各图依次为系统单位阶跃

6、响应曲线、单位阶跃响应误差曲线、零极点分布图、系统开环对数特性曲线、系统根轨迹特性曲线、奈奎斯特曲线。2.2 稳定性分析由闭环零极点分布图可知，系统闭环传递函数的极点全部位于s左半平面，因此系统闭环稳定。2.3 性能指标分析（1） 系统稳态误差=0。由以上分析可知在加入控制器之前，=0，符合设计要求。（2） 系统调节时间10s。在加入控制器之前，调节时间=38s，不满足设计要求，因此需要加入控制器来缩短调节时间，以提高系统响应速度。（3） 系统阶跃响应的超调量%0.95*C)&(y(i)0.98*C)&(y(i)1.02*C); i=i-1;endts=t(i)grid;MATLAB命令窗口中

7、运行后得到结果为：percentovershoot = 0.3270ts = 4.6300综上所述，选择Ki=0.001，Kp=9，Kd=0.2时系统各方面性能都能令人满意，所以可以作为PID控制参数。此时系统的开环传递函数为：系统性能验证： 加入PID之前 加入PID之后MATLAB命令窗口运行结果ans = RiseTime: 21.5386 SettlingTime: 38.5480 SettlingMin: 0.9001 SettlingMax: 0.9997 Overshoot: 0 Undershoot: 0 PeakTime: 79.1536MATLAB命令窗口运行结果ans =

8、 RiseTime: 2.3903 SettlingTime: 4.6403 SettlingMin: 0.9190 SettlingMax: 1.0032 Overshoot: 0.3214 Undershoot: 0 PeakTime: 5.1556四、设计总结通过以上分析，可以看出，从P控制一直到PID控制，系统的性能越来越好。同时，可以发现PID控制所起的作用，不是P、I、D三种作用的简单叠加，而是三种作用的相互促进。增大比例系数Kp一般将加快系统的响应，在有静差的情况下有利于减小静差，但是过大的比例系数会使系统有比较大的超调，并产生振荡，使稳定性变坏。所以调试时将比例参数由小变大，并

9、观察相应的系统响应，直至得到反应快、超调小的响应曲线。如果系统没有静差或静差已经小到允许范围内，并且对响应曲线已经满意，则只需要比例调节器即可。如果在比例调节的基础上系统的静差不能满足设计要求，则必须加入积分环节。增大积分时间I有利于减小超调，减小振荡，使系统的稳定性增加，但是系统静差消除时间变长。如果系统的动态过程反复调整还不能得到满意的结果，则可以加入微分环节。增大微分时间D有利于加快系统的响应速度，使系统超调量减小，稳定性增加，但系统对扰动的抑制能力减弱。在PID参数进行整定时如果能够有理论的方法确定PID参数当然是最理想的方法，但是在实际的应用中，更多的是通过凑试法来确定PID的参数。

10、将原系统以及P、PD、PI、PID四种控制方式下的阶跃响应曲线画在同一坐标系下如右图：五、心得体会这次实验，认识了自动控制领域最常用的PID控制，基本掌握了PID控制的基本规律，同时也认识到自动控制系统的复杂性。在利用MATLAB软件时经常会碰到一些新问题，而我们手头的资料有限，时间和精力有限，并不能解决所有问题。比如在PID控制时，一旦选定了Ki和Kd后，超调量随Kp的变化并不明显，这是我无法理解的，当Kp增加时，系统仅仅提高了响应的快速性，而超调量并没有显著的变化。又如，在PD控制时，当Kd和Kp取值足够大时，便可以使响应曲线完全理想化，即响应时间趋于0，超调量趋于0，在本系统中也满足足够的稳态精度，我就会这样怀疑，并不是所有系统采用PID控制效果一定比其他控制效果要好，等等。所有这些问题将在今后的学习和实验中寻求答案。六、参考文献（1）西安交通大学出版社 反馈控制问题使用MATLAB及其控制系统工具箱 （2）科学出版社自动控制原理第五版 胡寿松主编（3）重庆大学出版社控制系统计算机辅助设计 蔡启仲等编著 （4）西北工业大学出版社基于MATLAB7.x/Simulink/Stateflow系统仿真、分析及设计 贾秋玲、袁冬莉等编著（5）西安电子科技大学出版社MATLAB在自动控制中的应用 吴晓燕主编（6）机械工业出版社MATLAB/Simulink建模与仿真实例精讲