毕业设计（论文）-matlab在自动控制方面的应用研究.doc

matlab在自动控制中的应用研究学 生 毕 业 论 文（ 2007届）题目（中文） matlab在自动控制中的应用研究 （英文） application of matlab in automatic 系别： 物理与电子工程系 专业： 电气工程及其自动化 班级： 20071261 姓名： 学号: 2007126109 指导教师： 陈燕娟（高级实验师） 诚 信 声 明我声明，所呈交的毕业论文是本人在老师指导下进行的研究工作及取得的研究成果。据我查证，除了文中特别加以标注和致谢的地方外，论文中不包含其他人已经发表或撰写过的研究成果，我承诺，论文中的所有内容均真实、可信。毕业论文作者签名： 签名日期：2010年5月23日摘要：本文介绍了自动控制中matlab的一般运用，以及在自动控制主要系统中的具体应用。利用matlab强大的数值计算和图形绘制能力，作为程序主要设计方向。运用这种功能设计了“自动控制原理理论分析”的计算辅助程序。程序包含了稳定性分析、时域响应分析和根轨迹分析三个主要模块，这些模块改变了传统的工程繁杂运算，充分利用计算机集成软件的表现能力，将抽象问题具体化，而且程序易懂，操作简便，对自动控制理论的应用研究具有一定的实用性关键词：自动控制原理；matlab；稳定性分析；时域响应分析；根轨迹分析abstract: this article mainly introduces the general applications of matlab in the main system of automatic control. besides, with the help of the powerful functions like calculating and charting, the writer design a computer aided process called “theoretical analysis of automatic control”, which has three core modules including stability analysis, time domain analysis and root locus analysis. this process aims to change the traditional engineering complex operations, and to take advantage of computer integrated software to materialize some abstract problems, easing the operation so as to make it much easier for understanding and become more practical for analyzing the theory of automatic control.keywords: theory of automatic control; matlab; stability analysis; time domain analysis; root locus analysis目录1 概论11.1 设计目的及意义11.2 设计基本要求11.3 设计的基本框架和功能22 matlab软件的编程研究32.1 matlab软件的使用32.2 matlab软件在自动控制原理中的应用32.3 matlab在绘图方面的应用43 matlab程序的设计和调试43.1 利用matlab进行系统稳定性判定的设计4311自动控制原理中稳定性的判断4312 利用matlab程序设计进行系统稳定性的判断53.2 利用matlab进行系统时域分析的设计6321自动控制原理中的时域分析7322利用matlab进行系统的时域分析13323利用matlab设计绘制时域分析系统响应曲线的程序143.3 利用matlab进行根轨迹绘制的设计18331自动控制原理中的根轨迹分析19332利用matlab分析系统的根轨迹20333利用matlab设计绘制系统根轨迹图的程序21结束语24参考文献25附录26致谢31matlab在自动控制中的应用研究1 概论1.1 设计目的及意义自动控制原理是自动化学科重要的专业基础，分为经典控制理论和现代控制理论。经典控制理论是以传递函数为基础，利用时域分析法、频率分析法和根轨迹法对线性定常系统进行分析和设计。这是主要研究自动控制系统的一般规律，涉及控制系统的模型建立、系统分析，以及系统设计的基本理论和相关技术。其特点是具有一定的概括性和抽象性，数学知识和运算的应用多，计算复杂。因此，通过计算机的计算表现能力将抽象问题具体化。随着计算机技术的发展和应用，自动控制理论和技术在宇航、机器人控制等高新技术领域中的应用也愈来愈深入广泛。不仅如此，自动控制技术的应用范围现在已扩展到生物、医学、经济管理和其它许多社会生活领域中，成为现代社会生活中不可缺少的一部分。随着时代进步和人们生活水平的提高，建设高度文明和发达社会的活动中，自动控制理论和技术必将进一步发挥更加重要的作用1。而在自动控制工程实际中，可控制系统的结构往往很复杂，如果不借助专用的系统建模软件，则很难准确地把一个控制系统的复杂模型输入计算机，对其进行进一步的分析与仿真。所以需要一种高级科学分析与计算软件，具有强大的数值计算能力，而且速度快，使用方法简便，matlab是为数不多的完全满足这些要求和条件的软件工具，并且凭借它在科学计算方面的天然优势，建立了从设计构思实现到最终设计要求的可视化桥梁，大大弥补了传统设计与开发工具的不足， 是工程控制设计中不可多得的工具，并且在科研、生产和数学等领域具有广泛的应用前景和推广价值。因此，matlab在自动控制方面的应用具有相当广泛的研究意义以及达到更优化设计的目的。1.2 设计基本要求根据所学的matlab的知识来进行编程实现基于matlab的自动控制原理理论进行分析与设计（时域），本次设计预期完成以下功能：（1）运用matlab命令分析控制系统的性能及其稳定性；（2）利用时域分析法对控制系统（一阶系统和典型二阶系统）进行分析，并给出控制系统的瞬态性能指标；（3）利用根轨迹法对控制系统进行分析；（4）提供直观、形象的波形演示；1.3 设计的基本框架和功能matlab在自动控制原理中的应用研究稳定性分析时域响应分析根轨迹分析系统的根稳定性判断一阶系统二阶系统任意阶系统参数计算根轨迹分析图1 程序设计框架图描述各个模块的功能：（1）稳定性分析模块：包含系统的稳定性的判断和显示闭环系统所有的特征根，很直观地显示系统的稳定性判断是否正确，快捷地表示已知系统的稳定性。（2）时域分析响应模块：包括了一阶系统、典型二阶系统、任意阶系统三个主要界面，可以通过此界面绘制时域的响应曲线并且计算相关参数。（3）根轨迹分析模块：包含了绘制根轨迹的函数计算，可以通过此界面绘制系统的根轨迹。2 matlab软件的编程研究2.1 matlab软件的使用matlab环境(中文名是矩阵实验室)是matlab是math works公司推出的种面向工程和科学运算的交互式计算软件，经过近二十年的发展与竞争、完善，现已成为国际公认的最优秀的科技应用软件。matlab有三大特点:一是功能强大，它包括了数值计算和符号计算、计算结果和编程可视化、数学和文字统一处理、离线和在线计算等功能；二是界面友好、语言自然，matlab以复数矩阵为计算单元,指令表达与标准教科书的数学表达式相近；三是开放性强，matlab有很好的可扩充性，可以把它当作一种高级的语言去使用，用它容易地编写各种通用或专用应用程序1。simulink是matlab的一个重要的伴随工具，它通过对真实世界中的各种物理系统建立模型，进而用计算机实现仿真的软件工具。将simulink应用于自动控制系统，可以很容易构建出符合要求的模型，灵活的修改参数，方便的改变系统结构或进行转换模型，同时可以得到大量的有关系统设计的充分的、直观的曲线，这使得它成为国际控制界应用最广的首选计算机工具软件。正是因为matlab具有这些特点，因而被广泛使用，不仅成为世界上最受欢迎的科学与工程计算软件之一，而且成为国际上最流行的控制系统计算机辅助设计的工具。现在的matlab已经成为一种具有广阔应用前景的全新的计算机高级编程语言。2.2 matlab软件在自动控制原理中的应用自动控制是控制理论中理论性较强的技术基础课，在工科的专业培养方向都占有重要地位。在控制理论中涉及许多复杂的计算问题，并且还有时域分析中涉及许多根轨迹作图和伯德制图，而matlab以矩阵和向量为基本数据单位，具有十分突出的矩阵计算能力。同时matlab含有各种可选的工具箱，如神经网络、小波分析、模糊控制、信号处理、鲁棒控制等功能强大的计算辅助工具。在自主探索式学习过程中掌握用matlab来计算与模拟自动控制理论中一些难以理解的问题，使原本比较抽象的问题更加形象化。在自动控制原理理论的研究应用上，我们运用matlab语言解决稳态分析和作图问题的尝试，使计算和作图问题变得简单化，高效化。随着计算机技术的发展和应用，matlab对自动控制理论和技术在宇航、机器人控制等高新技术领域中的应用也愈来愈深入广泛。不仅如此，自动控制技术在生物、医学、经济管理和其它许多社会生活领域的应用越来越依靠matlab强大的功能处理来实现自动控制要求，这些应用成为现代社会生活中不可缺少的一部分。随着时代进步和人们生活水平的提高，建设高度文明和发达社会的活动中，matlab对自动控制理论和技术的应用和研究必将进一步发挥更加重要的作用。2.3 matlab在绘图方面的应用虽然数值计算、符号计算给人们的科学研究带来了很大帮助，很多时候人们难以直接从大量的数据中感受到它们的具体含义和内在规律。人们更喜欢通过图形的方式直观形象地感受科学结果的全局意义和内在本质。而matlab的出现解决了这个问题。matlab有很强的图形功能，可以方便地实现数据的视觉化，这使得matlab在工程方面的提升不仅仅在数据上而且具体化到图形绘制。强大的计算功能与图形功能的结合，为matlab在科学技术和自动控制图形绘制方面的应用提供了更加广阔的天地。3 matlab程序的设计和调试3.1 利用matlab进行系统稳定性判定的设计311自动控制原理中稳定性的判断稳定性是指控制系统在受到扰动信号作用，原有平衡状态被破坏后，经过自动调节能够重新达到平衡状态的性能。当系统在扰动信号作用（如电网电压波动，电动机负载转矩变化等）下偏离了原来的平衡状态时，若系统能通过自身的调节作用使得偏差逐渐见笑，重新回到平衡状态，则系统是稳定的；若偏差不断增加，即使扰动消失，系统也不能回到平衡状态，则这种系统是不稳定的，这表明稳定性是表征系统在扰动消失后的一种恢复能力，它是系统的一种固有特性。系统的稳定性又分为两种：一种是大范围的稳定，即初始偏差可以很大，但系统仍然稳定；另一种是小范围稳定，即初始偏差必须在一定限度内系统才稳定，超出了这个限定值则不稳定。对于线性系统，如果小范围内是稳定的，则它的大范围也是稳定的。而非线性系统不存在类似结论。通常而言，线性定常系统的稳定性表现为其时域响应的收敛性。当把控制系统的响应分为过渡状态和稳定状态来考虑时，若随着时间的推移，其过渡过程逐渐衰减，系统的响应最终收敛到稳定状态，则称该系统是稳定的；如果过渡过程是发散的，则该系统就是不稳定的2。任何一个自动控制系统正常运行的首要条件是，它必须是稳定的。因此，判别系统的稳定性和使系统处于稳定的工作状态，是自动控制的基本问题之一。稳定性是系统去掉扰动以后，系统自身的一种恢复能力，是系统本身所固有的特性。它仅仅取决系统的结构参数，而与初始条件及输入信号无关。根据分析我们知道，如果系统所有的闭环特征根（闭环极点）都分布在s平面左半部，则系统的暂态分量随时间增加逐渐消失为零，这种系统是稳定的。如果有一个或一个以上的闭环特征根是位于s平面右半部或虚轴上，则系统是不稳定的。综上所述，线性系统稳定的充分必要条件是：系统的闭环特征方程式所有的根，均位于s平面的左半部。虽然对于高阶系统，求其闭环特征方程的根是件很麻烦的工作，但是利用劳斯判据来判断需要设计很复杂的界面，又因为matlab具有很强大的计算能力，并且提供了直接求解代数方程根的函数，所以利用matlab来求闭环特征方程的根不是一件很难的事，同时本次设计对软件的运行速度的要求不是很高，只要起到演示的效果即可。经过求闭环特征方程的根和劳师判据来判断系统稳定性的优缺点的对比，我们在本次设计中选择求闭环特征方程的根来判断系统的稳定性。312 利用matlab程序设计进行系统稳定性的判断系统稳定是系统实现预期目标的一个重要因素。因此对系统进行稳定性判定是很有必要的。利用matlab软件对系统稳定性进行分析判定，获得系统传递函数的特征值分布，以此为依据进行判断，是判断过程更加方便、简单、快捷，避免繁琐的分析计算过程。同时，有数据的支持能使得结论的科学性和真实性更加清晰。具体利用matlab实现系统稳定性判定方法：在matlab中编制一个m文件，用于判断系统的闭环特征根是否都位于s平面的左半部，并显示闭环特 征的根，保存的名字为wending.m，并将此文件保存在matlab的work文件夹下。论文已给出wending.m文件，详细程序见附录。a=roots (sys.den1)sum=0;a,b=size(a);c=a*bfor i=1:c if real(a(c)0 sum=sum+1; endend由自动控制原理理论可知系统的传递函数标准型为：输入参数运行程序，运行结果如图2：图2 稳定性判断程序图3.2 利用matlab进行系统时域分析的设计对控制系统而言，其数学模型由微分方程和差分方程给出，因此可以从给定的初始值开始，通过某种算法逐步求出系统某一时刻的响应，从而丝线对控制系统的分析。此外，通过对系统的时域分析，可以求得系统响应的性能指标。在经典控制理论中，时域分析法是一种十分重要的分析和设计控制系统的方法，它包括系统稳定性分析、动态性能和稳态性能指标的计算等内容。时域分析法是通过传递函数、拉氏变换及其反拉氏变换求出系统在典型输入下的输出表达式，从而分析系统的时间响应的全部信息。与其他分析法相比较，时域分析法是一种直接分析法，具有直观和准确的优点，尤其适用于一、二阶系统的分析计算。而对于二阶以上的高阶系统则须采用频率分析法或根轨迹法，因此在本次设计中对于高阶系统的性能指标不进行计算。 与控制系统的其他matlab仿真一样，时域响应matlab的仿真方法也可以在matlab函数的指令方式下进行时域仿真，对于线性系统，matlab控制系统工具箱提供了若干函数完成线性系统的仿真。321自动控制原理中的时域分析（1）典型输入信号控制系统中常用的典型输入信号有：单位阶跃函数、单位斜坡函数、单位加速函数、单位脉冲函数及正弦函数。在典型输入信号作用下，任何一个控制系统的时间响应都由动态过程和稳态过程这两部分组成。相应地，控制系统在典型输入信号作用下的性能指标，通常也由动态性能指标和稳态性能指标这两部分组成。（2）时域性能指标动态过程又称瞬态过程或暂态过程，是指系统在典型输入信号作用下，其输出量从初始状态到最终状态的响应过程。系统在动态过程中所提供的系统响应速度和阻尼情况等用动态性能指标描述。时域中评价系统的暂态性能，通常以系统对单位阶跃输入信号的暂态响应为依据。为了评价系统的暂态性能，则规定其性能指标如图3所示。 图3 单位阶跃响应曲线图通常，在单位阶跃函数作用下，稳定系统的动态过程随时间t变化的指标成为动态性能指标。对于图3所示的单位阶跃响应，通常定义动态性能指标为以下几种2。上升时间（rise time）对于无振荡的系统，定义系统响应从终值的10%上升到90%所需的时间为上升时间；对于有振荡的系统，定义响应从零第一次上升到终值所需的时间为上升时间。缺省情况下，matlab按照第一种定义方式计算上升时间，但可以通过设置得到第二种方式定义的上升时间。峰值时间（peak time）响应超过其终值到达第一个峰值所需的时间定义为峰值时间。超调量（overshoot）响应的最大偏差量与终值的差与终值之比的百分数，定义为超调量，即超调量也称为最大超调量或百分比超调量。调节时间（setting time）响应到达并保持在终值2%或5%内所需的最短时间定义为调节时间。缺省情况，matlab计算动态性能时，取误差范围为2%，可以通过设置得到误差范围为5%时的调节时间。在上述的指标中，上升时间、峰值时间都是表征系统响应初始阶段速度的快慢；调节时间表征过渡时间（暂态时间）的持续时间，从整体上反应了系统的快速性；而超调量表征暂态过程的稳定性；稳态误差反映系统复现输入信号的最终精度。（3）一阶系统时域分析由一阶微分方程描述的系统，称为一阶系统，它的传递函数为：式中的t成为时间常数，它是表征系统惯性的重要参数。所以一阶系统是一个非周期的惯性环节。下面以单位阶跃信号为例来进行一阶系统的时域分析。当输入信号u=1时，u=1/s，系统输出量的拉氏变换为：对上式取拉氏反变换，得单位阶跃响应为： 图4 一阶系统响应图由此可见，一阶系统的阶跃响应是一条初始值为0，按指数规律上升到稳态值1的曲线，见图4。由系统的输出响应可得到如下的性能：单位阶跃响应曲线的初始斜率为这表明一阶系统的单位阶跃响应如果以初始速度上升到稳态值1，所需的时间恰好等t。根据暂态性能指标的定义可以求得调节时间为： （5%的误差带） （2%的误差带）延迟时间为： 上升时间为：其中峰值时间和超调量都为0。（4）二阶系统时域分析凡是可用二阶微分方程描写的系统称为二阶系统。在工程实践中，二阶系统不乏其例。特别是：不少高阶系统在一定条件下可用二阶系统的特性来近似表征。因此，研究典型二阶系统的分析和计算方法，具有较大的实际意义。典型的二阶系统图3.4为典型的二阶系统动态结构图，系统的开环传递函数为：系统的闭环函数为： 上式称为典型二阶系统传递函数，其中为典型二阶系统阻尼比，（或相对阻尼比）为无阻尼振荡频率或称自然角频率。系统闭环传递函数的分母等于零所得方程式称为系统的特征方程式。典型二阶系统的特征方程式为两个特征根为：当01，称为过阻尼状态。特征根为两个不相等的负实数。当=0，称为无阻尼状态。特征根为一对纯虚数。由此可得，和是二阶系统两个非常重要的参数，系统响应特性完全由这两个参数来决定的。二阶系统的阶跃响应3由于本次设计中只以单位阶跃响应为例来分析系统的参数，所以在这里不再像一阶系统那样来详细的分析各种输入信号的输出响应，在这里只以单位阶跃信号为例来分析典型二阶系统的输出响应及其性能指标。在单位阶跃函数作用下，二阶系统输出的拉氏变换为：求y(s)的拉氏变换，可得典型二阶系统单位阶跃响应。由于特征根与系统阻尼比有关。当阻尼比为不同值时，单位阶跃响应有不同的形式，下面只对处于欠阻尼情况下的二阶系统进行分析（因为本次设计的软件主要实现欠阻尼的系统分析）。当系统处于欠阻尼情况下，则00&n0 sum=sum+1; endendif sum=0disp(该系统是稳定的)disp(系统的闭环特征根为：);disp(a)%显示系统闭环特征根elsedisp(该系统是不稳定的)disp(系统的闭环特征根为：);disp(a)%显示系统闭环特征根disp(系统不稳定根的个数：);disp(sum)%显示不稳定根的个数end附录2、一阶系统的时域分析m文件，yijie.m程序内容：k=input(请输入传递函数k值：);t=input(请输入传递函数时间常数t值：);num=k;den=t 1;sys=tf(num,den);impulse(sys);t0=2.2*t;t1=4*t;title(系统的冲激响应);disp(该系统的上升时间为：);disp(t0);disp(该系统的调节时间为：);disp(t1);附录3、二阶系统时域分析m文件，erjie.m程序内容：a=input(请输入传递函数分子值：);b=input(请输入传递函数一阶系数：);c=input(请输入传递函数常数项：);num=a;den=1 b c;sys=tf(num,den);m=sqrt(c);n=b/(2*m);if n0&nq b=zeros(1,n-q),b; gs=n; else nq a=zeros(1,q-n),a;% 两个向量同维数 gs=q; end zm=a-b; gs=q; rt=roots(zm); h,l=size(rt); for i=1:h;ns=poly2sym(num,s); ds=poly2sym(den,s); k(i)=subs(abs(ds/ns),s,rt(i,1);%利用幅值条件求根轨迹增益 m(i)=k(