MATLAB自动控制课程设计报告.docx

《MATLAB自动控制课程设计报告.docx》由会员分享，可在线阅读，更多相关《MATLAB自动控制课程设计报告.docx（13页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、MATLAB自动控制课程设计报告 一、实训目的： 了解matlab 软件的基本特点和功能，熟悉其界面，菜单和工具条；掌握线性系统模型的计算机表示方法，变换以及模型间的相互转换。了解控制系统工具箱的组成，特点及应用；掌握求线性定常连续系统输出响应的方法，运用连续系统时域响应函数（impulse,step,lsim）,得到系统的时域响应曲线。 掌握使用MATLAB软件作出系统根轨迹；利用根轨迹图对控制系统进行分析；掌握使用MA TLAB软件作出开环系统的波特图，奈奎斯图；观察控制系统的开环频率特性，对控制系统的开环频率特性进行分析。 掌握MATLAB软件中simulink工具箱的使用；熟悉simu

2、link中的功能模块，学会使用simulink对系统进行建模；掌握simulink的仿真方法。 二、软件介绍： MATLAB是Mathworks公司开发的一种集数值计算、符号计算和图形可视化三大基本功能于一体的功能强大、操作简单的优秀工程计算应用软件。MATLAB不仅可以处理代数问题和数值分析问题，而且还具有强大的图形处理及仿真模拟等功能。从而能够很好的帮助工程师及科学家解决实际的技术问题。 MATLAB是矩阵实验室（Matrix Laboratory）的简称，是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括

3、MATLAB和Simulink两大部分。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。 MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MA TLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连matlab开发工作界面接其他编程语言的程序等，主要

4、应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。 MATLAB最重要的特点是易于扩展。它允许用户自行建立完成指定功能的扩展MATLAB函数（称为M文件），从而构成适合于其它领域的工具箱，大大扩展了MATLAB的应用范围。目前，MATLAB已成为国际控制界最流行的软件，控制界很多学者将自己擅长的CAD方法用MATLAB加以实现，出现了大量的MATLAB配套工具箱，如控制系统工具箱（control systems toolbox），系统识别工具箱（system identification toolbox），鲁棒控制工具箱（robust control to

5、olbox），信号处理工具箱（signal processing toolbox）以及仿真环境SIMULINK 等。 三、实训原理 在自动控制系统中，对于线性定常系统，可以用常系数线性激分方程加以描述。当给定输入的时间函数时，通过解微分方程，可以得出系统的输出响应。根据输出响应的数学表达式可以画出时间响应曲线，直观地反映出系统工作的动态过程。通常采用传递函数这种与微分方程等价的数学模型。 1、数学模型的建立 （1）、由结构图求传递函数，主要的函数有：parallel ，series ，feedback 等 （2）、模型变换 主要函数有：tf2ss ，ss2tf ，ss2zp ，tf2zp ，z

6、p2ss ，zp2tf 可以完成传递函数与状态空间模型间的变换。函数c2d ，c2dt ，c2dm ，d2c ，d2cm 可进行连续系统与离散系统间的变换。函数ctrbf ，obsvf 可对状态空间实现进行可控、可观测性分解。函数minereal 可完成传递函数的零、极相消及状态空间模型的最小实现。 2、时域响应 主要函数有：roots ，step ，impulse ，lsim ，initial 3、根轨迹 主要函数有：pzmap ，rlocus ，rlocfind ，sgrid ，zgrid 。 4、频率法 主要函数有：bode ，nyquist ，margin 当系统复杂连接时其等效的整体

7、传递函数可用simulink 来实现传递函数的求取。 四、实训内容 （1）、用matlab 语言编制程序，实现以下系统： 1）、G(s)= 2 26418 24523423+s s s s s s 输入程序: num=5 24 0 18; den=1 4 6 2 2; G=tf(num,den) 运行结果： Transfer function: 5 s3 + 24 s2 + 18 - s4 + 4 s3 + 6 s2 + 2 s + 2 2)、 G(s)= ) 523()1()66)(242332 2+s s s s s s s s 输入程序： n1=4*1 2; n2=1 6 6; n3=1

8、 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=conv(den1,den2),1 0); h=tf(num,den) 运行结果： Transfer function: 4 s 5 + 5 6 s4 + 288 s3 + 672 s2 + 720 s + 288 - s 7 + 6 s6 + 14 s5 + 21 s4 + 24 s3 + 17 s2 + 5 s （2）、两环节G1，G2串联，求等效的整体传递函数G(s) G1(s

9、)=32 +s G2(s)=1 227+s s 输入程序： num1=2; den1=1 3; sys1=tf(num1,den1); num2=7; den2=1 2 1; sys2=tf(num2,den2); G=sys1*sys2 运行结果： Transfer function: 14 - s3 + 5 s2 + 7 s + 3 （3）、两环节G1，,G2并联，求等效的整体传递函数G(s) G1(s)=32 +s G2(s)=1 227+s s 输入程序： num1=2; den1=1 3; sys1=tf(num1,den1); num2=7; den2=1 2 1; sys2=tf

10、(num2,den2); G=sys1+sys2 运行结果： Transfer function: 2 s2 + 11 s + 23 - s3 + 5 s2 + 7 s + 3 （4）、已知系统结构如图，求闭环传递函数。其中的两环节G1,G2分别为G1(s)=81 2100 s 32+s s G2(s)= 5 22 +s 输入程序： num1=3 100; den1=1 2 81; sys1=tf(num1,den1); num2=2; den2=2 5; sys2=tf(num2,den2); G=feedback(sys1,sys2) 运行结果： Transfer function: 6

11、s2 + 215 s + 500 - 2 s3 + 9 s2 + 178 s + 605 （5）、 已知某闭环系统的传递函数为G(s)=25 10296.1316.025 10+s s s s ,求其单位阶跃响 应曲线，单位脉冲响应曲线。 1）、单位阶跃响应 输入程序： num=10 25; den=0.16 1.96 10 25; G=tf(num,den) y=step(G); plot(y); figure(1) grid ; title(单位阶跃响应); xlabel(时间) ylabel(振幅) 运行结果： 图5.1.1系统的阶跃响应曲线2）、单位脉冲响应 输入程序： num=10

12、25; den=0.16 1.96 10 25; G=tf(num,den) y=impulse(G); figure(2) plot(y); grid ; title(单位脉冲响应) xlabel(时间) ylabel(振幅) 图5.1.2系统的脉冲响应曲线 （6）、典型二阶系统的传递函数为 2 2 2 2 ) ( n n n w s w s w s G + + = ,n w 为自然频率， 为阻尼比， 试绘制出当 =0.5，n w 分别取-2 ,0 ,2 ,4 ,6 ,8 ,10时该系统的单位阶跃响应曲线；分析阻尼比 分别为-0.5，-1时系统的稳定性。 输入程序： i=0; t=0:0.1

13、:10; for wn=-2:10 i=i+1; zeta=0.5; num=wn2;den=1,2*zeta*wn,wn2; sys=tf(num,den); y(:,i)=step(sys,t); end plot(t,y(:,1), t,y(:,2), t,y(:,3), t,y(:,4), t,y(:,5), t,y(:,6), t,y(:,7); grid xlabel(time(s) ylabel(step response) 图6.1.1 系统的单位跃阶响应曲线 输入程序：假设wn=10， num=100;zeta=-0.5,-1;t=0:0.1:10;y=zeros(lengt

14、h(t),2); for i=1:2 den=1,2*zeta(i)*10,100; sys=tf(num,den); y(:,i)=step(sys,t); end plot(t,y(:,1), t,y(:,2); f=1,2*zeta(i)*10,100; p=roots(f) 运行结果为： p = 10 10 所以系统稳定。 （7）、设有一高阶系统开环传递函数为 271.6635.0268.006.0359 .9436.1218.0016.0)(2323+= s s s s s s s G ,试绘制该系统的零极点图和闭环根轨迹图。 1）系统的零极点 输入程序： num=0.016 0.2

15、18 1.436 9.359; den=0.06 0.268 0.635 6.271; sys=tf(num,den); Pzmap(sys) Z P K=Pzmap(sys) 运行结果： Z = -10.4027 -1.6111 + 7.3235i -1.6111 - 7.3235i P = -5.7710 0.6522 + 4.2054i 0.6522 - 4.2054i K = 0.2667 图7.1.1 系统的零极点图 2）系统的闭环根轨迹 输入程序： num=0.016 0.218 1.436 9.359; den=0.06 0.268 0.635 6.271; sys=tf(num

16、,den); rlocus(sys) 运行结果: 图7.1.2 系统的闭环根轨迹图 (8)、单位反馈系统前向通道的传递函数为： s s s s s s s s s s s G +=23456 2345101052 81282)(,试绘制该系统的Bode 图和Nyquist 曲线，说明软件绘制曲线与手动绘制曲线的异同。 1）、绘制该系统的Bode 图： 输入程序： num=2 8 12 8 2; den=1 5 10 10 5 1 0; sys=tf(num,den); margin(sys) 运行结果： 图8.1.1 系统的Bode图 2)、系统的Nyquist图 输入程序： num=2 8 12 8 2; den=1 5 10 10 5 1 0; sys=tf(num,den); nyquist(sys)