典型3阶系统的二阶参考.docx

编 号：_审定成绩：_ 课程设计报告课程设计题目:_ 典型3阶系统的二阶参考_ _ _ 模型设计及仿真研究_ 单 位（系别）：_ 自动化_ 学 生 姓 名：_ 李春斌_ 专 业：_电气工程及其自动化_ _ 班 级：_（1）班_ 学 号：_0511110110_ 指 导 教 师：_汪纪锋_ 填表时间： 2015 年5 月重庆邮电大学移通学院教务处制重庆邮电大学移通学院毕业设计(论文)任务书设计(论文) 题目28 典型3阶系统的二阶参考模型设计及仿真研究 学生姓名 李春斌 系别 自动化系 专业 电气工程及其自动化 班级 05111101指导教师 汪纪锋 职称 教授 联系电话 42871150 教师单位 自动化系 主 要 研 究 内 容 、 方 法 和 要 求研究典型型3阶系统动静态性能特性并完成综合设计，以达到使该系统满足工程实际性能指标的要求。基本要求：1.运用经典控制理论中的频域理论方法，分析给定的典型系统基本特性；2.根据提出的系统实际工艺指标，运用二阶参考模型法综合设计一套满足系统性能指标闭环系统方案；3.应用MATLAB/SIMULINK(或物理模拟)对设计系统进行仿真验证，给出仿真结果图。 进 度 计 划第4周第6周：完成系统建模，性能分析并提出校正设计方案；第4周第9周：按性能指标进行系统的综合设计；第10周第14周：完成原系统及校正后系统仿真研究，并作比较研究；第15周第17周：撰写论文、修改论文，完成答辩。 主 要 参 考 文 献1.汪纪锋、党晓圆，现代控制理论，人民邮电出版社，2013.11；2.郑大钟，线性系统理论（第2版），清华大学出版社，2008.04；3.Chi-Tsong Chen,Linear System Theory and Design,HOLT.RINEHAET AND WINSTON,2001.03；4.熊晓君自动控制原理实验教程（硬件模拟与MATLAB仿真），机械工业出版社，2009.01；5.方水良现代控制理论及其MATLAB实践，浙江大学出版社，2009.6；6.相关学术期刊等。指导教师签字： 汪纪锋 年 月 日教研室主任签字： 年 月 日备注：此任务书由指导教师填写，并于毕业设计(论文)开始前下达给学生。重庆邮电大学移通学院毕业设计任务书（简明）技术资料一、设计题目：题目28 典型3阶系统的二阶参考模型设计及仿真研究二、系统说明：设三阶系统开环结构如下G0sysus=1ss+12s+1三、系统参量：系统输入信号：u（t）；系统输出信号：y（t）；四、设计指标：1.设定：在输入为r(t) = u(t) = a + bt,(其中：a=4 b=1/sec)2.在保证稳态误差essv0.2的前提下，要求动态期望指标：%5%；ts5s五、设计要求：1.分析原系统性能；2.试用频率特性法按二阶参考模型设计满足系统性能指标的闭环系统；3.绘制系统校正前后的物理模拟仿真图。重庆邮电大学移通学院 自动化系指导教师：汪纪锋2014.12摘 要目录前言第一章 控制系统的简介第一节 自动控制系统的性能指标性能指标，是在分析一个控制系统的时候，评价系统性能好坏的标准。对于一个控制系统首要的要求是系统的绝对稳定性。在实际控制系统中，往往由于具体的对象不同或控制任务的指标不同，而对控制系统性能指标的要求有所不同。话虽如此，但仍可以对各种控制系统的性能指标概括为答题三个方面，即，快（快速性）、稳（稳定性）、准（准确性）。具体而言，对于定值控制系统，要能够迅速克服外界干扰的影响，使被控对象能够准确地恢复至给定值。系统性能的描述，可以分为动态性能和稳态性能。被控量不随时间变化的平衡状态称为“稳态”；而被控量随时间变化的不平衡状态称为“动态”。在系统的动态过程中，其被控量是不断变化的，这一随时间变化的过程称为动态过程，也叫做控制过程。在此次设计中，主要研究三阶定常线性系统。其性能指标包括动态和静态指标两部分。静态指标即误差系数，动态性能指标包括上升时间tr、峰值时间tp、超调量Mp、调节时间ts等。第二节 系统概述对于原系统物理仿真结构首先第一个环节将是一个比较环节，将输入信号与由输出反馈回的信号进行比较，以改善放大器的静态和动态性能；第二个是比例环节，对输出信号进行一定量的放大；第三个是积分环节，可以是系统的跟踪能力增强，积分环节是当输入信号为零时，输出信号才能保持不变，而且能保持在任何位置上。在控制系统中，引用积分环节可以消除被控量的偏差。第四个是惯性环节，由于惯性环节系统的阻力，一开始输出并不与输入同步按比例变化，直到过渡过程结束,输出才能与输出保持比例，从而保证了控制过程作无差控制。第五个是惯性比例环节。第六个是反馈环节，根据输入与输出在广义上是否相等来调节系统使之误差减小。第三节 设计基本要求基于频率特性法，按二阶参考模型法（即=0.707）设计满足系统性能指标的闭环系统。第四节 系统基本指标1.设定：在输入为rt=ut=a+bt,(其中：a=4 b=1/sec)2.在保证稳态误差essv0.2的前提下，要求动态期望指标：%5%；ts5s第二章 系统建模第一节 各环节模型建立1.比例环节系统结构图如下图（1），结构框图如图（2）传递函数:G1=-RfR0=-K1 图（1） 图（2）2.积分环节系统结构图如下图（3），结构框图如图（4）传递函数: G2s=-1R0C0S 图（3） 图（4）3.惯性环节系统结构图如下图（5），结构框图如图（6）传递函数： G3s=-R1R01+R1C1S 图（5） 图（6）4.反相器系统结构图如下图（7），结构框图如图（8）传递函数： G4s=-R2R0=-1 图（7） 图（8）第二节 系统数学模型三阶系统结构如下：G0sysus=1ss+12s+1开环传递函数为：G0s=1ss+12s+1闭环传递函数为：Gs=G0s1+G0s=1ss+12s+1+1特征方程为：Ds= ss+12s+1+1=0第三章 系统分析第一节 二阶系统的数学模型二阶系统的开环传递函数为：G0s=2ss+2n闭环传递函数为：Gcs=n2s2+2ns+n2闭环传递函数的分母多项式等于零的代数方程式称为二阶系统的闭环特征方程，即：s2+2ns+n2=0闭环特征方程的两个根称为二阶系统的特征根，即s1,2=-nn2-1上述二阶系统的数学模型中有两个特征参数 和 n，其中 称为二阶系统的阻尼比， n称为二阶系统的无阻尼振荡频率（rad/s）。二阶系统的系统分析和性能描述，基本上是以这两个特征参数来表示的。上述二阶系统的特征根表达式中，随着阻尼比 的不同取值，特征根si有不同类型的值，或者说特征根si在s平面上位于不同的位置，共有以下五种情况。1时，特征根为一对不相等的负实根，位于s平面的负实轴上，使得系统的响应表现为过阻尼的。=1时，特征根为一对相等的负实根，也是位于s平面的负实轴上，系统的响应表现为临界阻尼的。01时，特征根为一对带有负实部的共轭复数根，位于s平面的左半平面上，使得系统的响应表现为欠阻尼的。=0时，特征根为一对纯虚根，位于s平面的虚轴上，系统的响应表现为无阻尼的。 num=1; den=2 3 1 1; G=tf(num,den) Transfer function: 1-2 s3 + 3 s2 + s + 1 p=eig(G)p = -1.3982 -0.0509 + 0.5958i -0.0509 - 0.5958i p1=pole(G)p1 = -1.3982 -0.0509 + 0.5958i -0.0509 - 0.5958i r=roots(den)r = -1.3982 -0.0509 + 0.5958i -0.0509 - 0.5958i分析：系统的特征根全部位于s平面的左半平面，很显然该系统是稳定的。第四节 动态性能分析系统的阶跃响应一般曲线如图所示。图3.4.1 二阶欠阻尼系统的单位阶跃响应根据图中所展示的响应特性，可以定义如下的性能指标。（1）上升时间tr上升时间tr指阶跃响应ct上升至稳态值所需要的时间。考虑到不敏感区或者允许误差，有时取为从稳态值的10%上升至90%时所需要的时间。（2）峰值时间tp峰值时间tp指阶跃响应ct从运动开始至到达第一峰的时间。（3）超调量Mp超调量Mp指系统响应的第一峰值超出稳态值的部分，将其取百分比，可以表示为Mp=ctp-cc100%=ctp-11100%=-1-2100%（4）调节时间ts调节时间ts指阶跃响应ct达到稳定值的时间。评价系统整个过程的响应速度，是速度和阻尼程度的综合指标。（5）稳态误差lss稳态误差是当时间t趋于无穷时，系统希望的输出与实际的输出之差，表示为lss。从上述系统阶跃响应性能指标可以看出，各个时间指标反映了系统的快速性。其中，上升时间tr、峰值时间tp反映了系统的初始快速性；而调节时间ts反映了系统的总体快速性；超调量Mp描述了系统的平稳性；稳态误差lss描述了系统的准确性。3.4.1动态平稳分析这里指超调量Mp是否满足要求。特征方程为：Ds= ss+12s+1+1=0则有：Ds=2s3+3s2+s+1=Ds为闭环特征方程为阻尼比n为系统频率3.4.2动态快速性分析3.4.3原系统的图像分析在频率分析法中，以控制系统的频率特性作为数学模型，以波德图或者其他图表作为分析工具，来研究。分析控制系统的动态性能与温蔼性能。频率特性曲线有两种种表示方法：极坐标图、对数坐标图。（1）极坐标图极坐标图又称为幅相图、奈奎斯特图（Nyquist）.Nyquist曲线。频率特性Gj是输入信号频率w的复变函数，用实坐标与虚坐标可表示为Gj=ReGj+ImGj，当频率从-变到+时， Gj在由实轴与虚轴构成的复平面上走过的根轨迹就称为Gj的极坐标图。MATLAB在绘制奈奎斯特曲线时，频率是从连续变化的，而在自动控制原理的教材中，一般只绘制频率从部分的曲线。可以分析得出，曲线在范围与内，是以横轴为镜像的。极坐标图常用于频域稳定性分析的作图中。（2）对数坐标图对数坐标图又称为波德图。频率特性的矢量表达式为Gj=GjargGj=A其中A称为频率特性的幅频特性，是的函数。当由0+时，它展示了Gj幅值的变化规律。称为频率特性的相频特性，也是的函数。当由0+时，它展示了Gj相位角的变化规律。对数幅频特性L和对视相频特性两条曲线统称对数频率特性，又称波德图。原系统的波德图如下：系统开环结构如下G0sysus=1ss+12s+1使用MATLAB编写程序观察原系统的频率特性及阶跃响应。程序如下： s=tf(s); G=1/(s*(s+1)*(2*s+1) Transfer function: 1-2 s3 + 3 s2 + s bode(G) nyquist(G)程序运行结果得到系统Bode图和Nyquist曲线，分别如下图 3.4.2和3.4.3所示。图 3.4.2 校正前的系统Bode图图3.4.3 校正前的系统Nyquist曲线由原系统波德图分析可知：原系统的相频特性曲线一直位于-270之上，所以稳定裕度c0，且幅频特性曲线Lg0，显然系统是稳定的。由原系统的Nyquist曲线分析可知：原系统不存在右极点，即P=0，且在（-1，j0）点，正穿越和负穿越个数为0，所以系统是稳定的。3.3.4截止频率的计算c令对数幅频特性L过0dB的频率为c，则有20lg2lg1-lgc=-60解得c=1.23.3.5相位裕度计算c定义相位裕度c，计算公式如下c=180+cc=180+1.2c=40第四章 系统的综合设计第一节 超前校正网络具有超前特性，电路图如图所示。图4.1.1超前校正网络传递函数为GDs=1+s1+s其中T=R1C，=R2R1+R20，即相位角总是超前的；（2）轨迹为上半圆。使用相位超前校正可以改善系统的动态品质，减小超调量，缩短响应时间。超前校正的目的是改善系统的动态性能，实现在系统静态性能不受损的前提下，提高系统的动态性能。通过加入超前校正环节，利用其相位超前特性来增大系统的相位裕度，改变系统的开环频率特性。一般使校正环节的最大相位超前角出现在新的穿越频率点。第二节 滞后校正网络具有相位滞后特性，滞后校正网络电路图如图所示。图4.2.1 滞后校正网络传递函数为G1s=1+Ts1+Ts其中T=R2C，=R1+R2R21则频率特性为G1j=1+jT1+jT极坐标起点为A0=10=0终点为A=1=0可以证明：（1）：00，即相位总是滞后的；（2）轨迹为下半圆。在滞后校正中，我们利用的是滞后网络在高频段的衰减特性，而不是相位滞后特性。滞后校正通过加入滞后校正环节，使系统的开环增益有较大幅度增加，同时又使校正后的系统动态指标保持原系统的良好状态。他利用滞后校正环节的低通滤波特性，在不影响校正后系统低频特性的情况下，使校正后系统中高频段增益降低，从而使其穿越频率前移，达到增加系统相位裕度的目的。第三节 滞后超前校正图4.3.1 超前滞后校正（1） 在00时，校正网络起滞后网络的作用；在0时，它起超前网络的作用；（2） 网络的超前部分增加了相位超前角，增加了系统的相位裕量，改善了系统的动特性；（3） 网络的滞后部分，在剪切频率附近和以上，将产生衰减，因而它允许在低频段增加增益，故改善了系统的稳态特性;（4） 因而系统在动态和稳态均不能满足要求时,可以采用滞后-超前校正。（5） 滞后-超前校正的基本原理是利用校正装置的超前部分增大系统的相位裕度，同时利用其滞后部分来改善系统的稳态性能。滞后超前校正适用于对于校正后系统的动态和静态性能有更多更高要求的场合。施加滞后超前校正环节，主要是利用其超前部分增大系统的相位裕度，以改善系统的动态性能；利用其滞后部分改善系统的静态性能。以上3种不同的校正方法的一般性设计步骤如下：1) 根据静态性能指标，计算开环系统的增益。之后求取校正前后系统的频率特性指标。并与设计要求进行比较。2) 确定校正后期望的穿越频率，具体值的选取与所选择的校正方式相适应。3) 根据待设计的校正环节的形式和转折频率，计算相关参数，进而确定校正环节。4) 得出校正后系统。检验系统是否满足设计要求。如果不满足则从第二步从新开始。第四节 校正方案的确定由于调节时间，超调量都不满足要求，所以可以使用相位超前校正方案。第五节 设计校正量与校正装置4.5.1计算最大相位超前角mm=c-c0+520=4.5.2计算低频衰减=1-sinm1+sinm=4.5.3确定新的开环截止频率cL0=-1220log101=4.5.4确定截止频率1=1T=c=2=1T=1c=4.5.6补偿增益KcKc=1=4.5.7校正后的开环频率特性4.5.8校正后的截止频率和稳定裕度4.5.9校正后的系统波德图和Nyquist曲线图第五章 物理模拟仿真第一节 原系统物理仿真物理仿真，又称物理效应仿真，是指仿真的过程以物理性质和几何形状相似为基础，其他性质不变的仿真。控制系统是由各种元件或部件按照一定的要求组合在一起的。从形式结构上看，有各种各样的原件，但从动态性能或数字模型上看，却可以分为几个基本环节，这样对于我们研究系统和分析系统带来很大的便利，便于理解各个部件对于各个环节或是系统的性能影响。第一章已经对于各个环节进行了讲解，由于各环节之间没有负载效应，故各环节之间可以直接连接，并且个环节之间的放大电路可以合并为一个放大电路，从而可得出原系统物理模拟结构图如下图所示；经过校正后的系统物理模拟结构图如下图所示：第二节 MATLAB仿真5.2.1MATLAB的简介Matlab是起源于美国MathWorks公司发布主要面向数值计算、科学数据可视化以及交互式程序设计的高技术计算语言。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境之中，为科学研究、为适应众多专业科技人员的需要；MathWorks同时提供了数十个应用工具箱为科学和工程领域各类特殊问题及应用定制MATLAB运行环境；并为全面解决复杂数值计算问题以及CAD研究等提供了综合解决方案。MATLAB仿真在科学研究中的地位越来越高，如何利用MATLAB仿真出理想的结果，关键在于如何准确的选择MATLAB的仿真。Matlab的应用范围非常广，包括信号和图像处理、通讯、控制系统设计、调试和测量、财务建模和分析以及计算机生物等众多应用领域。附加的工具箱（单独提供的专用Matlab函数集）扩展了Matlab环境，以解决这些应用领域内特定类型的问题。Matlab可以用来进行以下各种工作：1. 数值分析2. 数值和符号计算3. 工程与科学绘图4. 控制系统的设计与仿真5. 数字图象处理6. 数字信号处理7. 通讯系统设计与仿真8. 财务与金融工程5.2.2MATLAB/SIMULINKSimulink是面向系统信号流图进行组态的仿真软件，通过Simulink模块库建立系统的仿真模型，可直观，方便地对系统进行动态仿真。SIMULINK是MATLAB软件的扩展，它是实现动态系统建模和仿真的一个软件包，它与MATLAB语言的主要区别在于，其与用户交互接口是基于Windows的模型化图形输入，其结果是使得用户可以把更多的精力投入到系统模型的构建，而非语言的编程上。具有以下功能；1） 用绘制方框图代替编写程序，结构和流程清晰；2） 智能化得建立和运行仿真，仿真精细，贴近