《自动化专业方向设计》报告

1、西南科技大学自动化专业方向设计报告设计名称： 基于Matlab的控制系统分析 姓 名： 学 号: 班 级： 指导教师： 起止日期： 2014.11.5-2014.12.5 西南科技大学信息工程学院制15方 向 设 计 任 务 书学生班级： 学生姓名： 学号： 设计名称： 起止日期： 2014.11.5-2014.12.5 指导教师： 设计要求：电加热的炉温控制系统如图所示，其目的是保持电路内温度恒定，被控量为电路内温度T，控制量为给定电位计的电压炉温控制系统被控对象电加热炉，炉内变化量对控制电压变化量之间的传递函数为，其中为电加热炉时间常数，控制器其它部分还包括放大器控制、电机调控、减速器控制

2、部分（齿轮）、测速发电机转化控制、热电偶反馈控制，取，试求出系统开环和闭环传递函数，并用MATLAB绘制Bode图和Nyquist图，判断系统是否稳定。1、设计内容1） 学习自动控制系统的开环、闭环函数，Bode图及Nyquist图2） 学习Matlab软件，应用Matlab绘制系统Bode图和Nyquist图方 向 设 计 学 生 日 志时间设计内容2014.11.5获得设计题目并与老师联系2014.11.6-11.13完成设计的大致方案2014.11.14-11.21完成传递函数的求解并写成报告形式2014.11.22-11.27Matlab程序编写、图形绘制、系统性能分析以及报告的完成2

3、014.11.28-12.3交予老师检查签字2014.12.5答辩基于Matlab的控制系统分析摘要 稳定是控制系统的重要性能，也是系统能够正常运行的首要条件，否则系统就会无法正常工作，甚至可能导致设备毁坏，造成重大损失。随着科学技术的不断发展，必将对控制理论、技术、系统与应用提出越来越多、越来越高的要求，因此有必要进一步加深加强对稳定性方面的研究，而MATLAB为此提供了可能。实践表明它的确是一个功能强大、形象逼真、便于操作的软件工具。本文就主要是利用MATLAB来对电加热的炉温控制系统进行稳定性分析。关键词：控制系统分析 MATLAB 炉温控制系统Analysis of control s

4、ystem based on MATLABAbstract:Stability is the important performance of control system, also the first condition of a system to normal operation. Otherwise the system will not be able to work properly, and may even lead to equipment damage, resulting in significant lose. With the continuous developm

5、ent of science and technology, the requirements of the control theory, technology, systems and applications is bound to higher and higher. So it is necessary to take further steps in the research of stability, and MATLAB makes

6、0;it possible. The practice shows that it is indeed a powerful, vivid, easily to operate software tool. This paper is mainly about a stability analysis of an electric-controled heating system based on MATLAB. Keywords: analysis of control system MATLAB electric-controled heating system 一、设计目的和意

7、义在科学技术取得突飞猛进的发展的过程中，自动控制始终担负着重要的角色。自动控制作为一种技术手段已广泛地应用与工业、农业、国防乃至日常生活和社会科学的许多领域。自动控制技术的广泛应用，不仅可以改善工作条件，减少劳动强度和提高生产效率，而且在人类征服自然、探知未来、建设高度文明的新社会等方面有着重要的意义。随着经济以及科技、国防事业的发展和人民生活水平的提高，自动控制技术所起的作用越来越重要，自动控制技术本身也将会得到进一步的发展。作为一个学自动化的人员，了解和掌握自动控制方面的知识是十分必要的，如自动控制的基础知识、MATLAB的使用、控制系统的稳定性判断等。2、 主要工作 本设计基于电加热的炉

8、温控制系统，要求利用MATLAB绘制Bode图和Nyquist图来对系统进行稳定性分析，通过本次设计掌握自动控制系统的基本概念和稳定性分析。 已知：电加热的炉温控制系统如图所示，其目的是保持电路内温度恒定，被控量为电路内温度T，控制量为给定电位计的电压图1.炉温控制系统被控对象电加热炉，炉内变化量对控制电压变化量之间的传递函数为，其中为电加热炉时间常数，控制器其它部分还包括放大器控制、电机调控、减速器控制部分（齿轮）、测速发电机转化控制、热电偶反馈控制，取，试求出系统开环和闭环传递函数，并用MATLAB绘制Bode图和Nyquist图，判断系统是否稳定。具体要求如下： 1.学习自动控制系统的开

9、环、闭环函数，Bode图及Nyquist图； 2.学习Matlab软件； 3.应用Matlab绘制系统Bode图和Nyquist图； 4.完成专业方向设计报告。3、 开环闭环传递函数 1.放大器控制部分：将给定温度电压信号与热电偶测量的实际炉温比较后的偏差信号进行放大，其传递函数为。 2.电动机调控部分：偏差电压信号驱动电动机正向或反向旋转产生角速度，其传递函数为。 3.齿轮作为减速器控制部分：其齿数比（转速比）为，其传递函数为。 4.测速发电机转化控制部分：测速发电机的输出电压与其转速成正比，其传递函数为。 5.电炉（被控对象）加热控制部分：作为被控对象，将电热丝散发的热量用以控制炉子的温度

10、，其传递函数为。 6.热电偶作为反馈控制部分：热电偶作为测量元件测量炉内实际温度并与给定温度进行比较产生偏差信号，其传递函数为。 综合以上各控制部分得出系统控制模型，如下图所示：图 2.系统控制模型经计算，系统前向通道传递函数为： 开环传递函数为： 闭环传递函数为： 根据已知可以求得：前向通道传递函数；开环传递函数为；闭环传递函数为。4、 时域分析方法判断稳定性 由系统开环闭环传递函数可知该系统为线性系统，而线性系统稳定的充分必要条件是闭环系统特征方程的所有根都具有负实部或者闭环传递函数的极点均位于左半平面。对此，判别系统稳定性的基本方法有：1) 劳斯-赫尔维茨判据；2) 根轨迹法；3) 奈奎

11、斯特判据；4) 李雅普诺夫第二方法。 利用时域方法判断稳定性主要依据劳斯-赫尔维茨判据。它根据系统特征方程的各项系数，直接利用代数方法判断特征方程的根是否全部位于左半平面，以此判断系统的稳定性。系统的开环传递函数为，可以得出系统特征方程为。列出劳斯表： 可见，劳斯表第一列中所有元素都大于零，故系统稳定。5、 Bode图与Nyquist图 Bode图与Nyquist图均是线性系统的频率特性的几何表示方法。Nyquist曲线又称为幅相频率特性曲线，在复平面上以极坐标的形式表示。Bode图又称为对数频率特性曲线，由对数幅频和对数相频两条曲线组成。对数频率特性曲线的横坐标是频率，并按对数分度，单位是弧

12、度/秒（rad/s）。对数幅频特性曲线的纵坐标表示对数幅频特性的函数值，均匀分度，单位是分贝（dB）。对数相频特性曲线的纵坐标表示相频特性的函数值，均匀分度，单位是度。系统的Bode图与Nyquist图如下：图 3.炉温控制系统Bode图图 4.炉温控制系统Nyquist图6、 利用频域分析方法判断稳定性（1）利用Bode图判断稳定性 由系统开环传递函数可知，系统有3个开环极点，0个开环零点。系统的3个极点分别为-1、-、-，全部位于左半平面，故该系统为最小相位系统。由幅值裕度和相角裕度图（图5）可知： 幅值裕度=2.85dB>1; 相角裕度=11deg>0。图 5.幅值裕度与相角

13、裕度 又，对于最小相位系统，相角裕度大于零，幅值裕度大于1，系统稳定。所以，该系统稳定。（2）利用Nyquist图判断系统稳定性根据Nyquist图判断系统稳定性主要根据奈奎斯特稳定判据。由上可知，系统的开环传递函数在平面的原点及虚轴上无极点，此时，奈奎斯特判据可表述为：反馈控制系统稳定的充要条件是当从变化时，奈氏曲线反时针包围临界点的圈数等于开环传递函数在右半平面的极点数，即；否则闭环系统不稳定。由图6可得出系统Nyquist曲线与实轴的交点为且不包围点，故Nyquist曲线反时针包围临界点的圈数；又系统开环传递函数在在右半平面的极点数。由此可得出，该系统稳定。图 6.Nyquist曲线与左

14、半s平面交点图七、系统的单位阶跃响应性能分析 炉温控制系统的单位阶跃响应曲线如图7。由图可知系统的稳态值。图 7.炉温控制系统单位阶跃响应曲线由MATLAB程序算出系统响应动态性能指标相应参数如下：峰值时间：；上升时间：；调节时间：；超调量：%；综合以上各项动态指标分析可知，调节反应时间足够灵活，峰值很快达到，故系统动态性能良好。八、结论综上所述，此炉温控制系统稳定且动态性能良好。9、 Matlab程序编写、（1） 系统阶跃响应曲线程序： num=12 12;den=6 11 6 8.2;y,x,t=step(num,den);peak,k=max(y) format compactp0=y(

15、length(t)overshoot=(peak-p0)/p0tp=t(k)n=1;while y(n)<1n=n+1;endtr=y(n)m=length(t)while (y(m)>0.98)&(y(m)<1.02)m=m-1;endts=t(m)step(num,den,t)grid ontitle('electric heating control system');xlabel('t(s)');ylabel('T(t)');运行结果为：peak = 2.8911k = 30p0 = 1.4623overshoot = 0.9771tp = 3.1751tr = 1.0166m = 851ts = 93.0644（2） Bode图程序编写： >> num=7.2;>> den=6 11 6 1;>> bode(num,den);>> grid;（3） Nyquist图程序编写：>> num=7.2;>> den=6 11 6 1;>> nyquist(num,den);>> grid;