自动控制系统的时域频域分析

1、目录摘要第一章 绪论11.1 自动控制理论发展概述11.2 Matlab简介2第二章 控制系统的时域分析与校正22.1 概述22.2 一阶系统的时间响应及动态性能32.3 二阶系统的时间响应及动态性能42.4高阶系统的阶跃响应、动态性能及近似11第三章 控制系统的频域分析与校正133.1 概述133.2 频率特性的表示方法143.3 频率特性的性能指标153.4 典型环节的频率特性17第四章 结论23课程设计总结24参考文献25附录26 摘要系统利用Matlab进行控制系统时域与频域的分析与设计，对控制系统的给定数学模型，研究系统性能与系统结构、参数之间的关系。其仿真过程是以某种算法从初态出发

2、，逐步计算系统的响应，最后绘制出系统的响应曲线，即可分析系统的性能。自动控制系统的计算机仿真是一门涉及到计算机技术、计算数学与控制理论、系统辨识、控制工程以及系统科学的综合性学科。控制系统仿真就是以控制系统的模型为基础，主要用数学模型代替实际的控制系统，以计算机为工具，对控制系统进行实验和研究的一种方法。控制系统最常用的时域分析法，就是在输入信号的作用下，求出系统的输出响应。系统采用单位阶跃响应为输入信号，求出各典型环节（一阶、二阶及高阶）的输出响应，分析各响应在阻尼比和固有频率变化时对输出响应的影响，从而可以选择最优方案，提高系统的快速性。而频域分析法是应用频率特性研究控制系统的一种经典方法

3、，以此可直观的表达出系统的频率特性，其主要方法有Bode图、Nyquist曲线、Nichols图，由于编写M文件时三种方法只需改变固定的命令，所以系统主要研究Bode图。同样是研究响应的典型环节，及比例、微分、积分、惯性、二阶振荡与高阶环节，分析其对数幅频特性与对数相频特性。经过对两种分析方法的对比与分析，得出了时域分析法与频域分析法的关系与区别。若已知控制系统的闭环传递函数，另外系统的阶次不是很高时，采用时域分析法较合适；而如果系统的开环传递函数未知，或者系统的阶次较高，就需采用频域分析法。通过对控制系统的仿真与分析从本质上区分了时域分析法和频域分析法的利弊，从而对不同的系统可以快速的找到合

4、适的方法，达到实验的预期目的。关键词：自动控制系统；时域/频域分析；Matlab 第一章 绪论1.1 自动控制理论发展概述自动控制理论是在人类征服自然地生产实践活动中孕育、产生，并随着社会生产和科学技术的进步而不断发展、完善起来的。早在古代，劳动人民就凭借生产实践中积累的丰富经验和对反馈概念的直观认识，发明了许多闪烁控制理论智慧火花的杰作。我国北宋时代苏颂和韩公廉利用天衡装置制造的水运仪象台，就是一个按负反馈原理构成的闭环非线性自动控制理论；1681年Dennis Papin发明了用做安全调节装置的锅炉压力调节器；1765年俄国人普尔佐诺夫发明了蒸汽锅炉水位调节器。1788年，英国人瓦特在他发

5、明的蒸汽机上使用了离心调速器，解决了蒸汽机的速度控制问题，引起了人们对控制技术的重视。之后，人们曾经试图改善调速器的准确性，却常常导致系统产生振荡。1868年，英国物理学家麦克斯韦通过对调速系统线性常微分方程的建立与分析，解释了瓦特速度控制系统中出现的不稳定问题，开辟了用数学方法研究控制系统的途径。此后，英国数学家劳斯和德国数学家古尔维茨独立的建立了直接根据代数方程的系数判别系统稳定性的准则。这些方法奠定了经典控制理论中时域分析法的基础。1932年，美国物理学家乃奎斯特研究了长距离电话信号传输中出现的失真问题，运用了复变函数理论建立了以频率特性为基础的稳定性判据，奠定了频率响应法的基础。随后伯

6、德和尼克尔斯进一步将频率响应法加以发展，形成了经典控制理论的频域分析法。之后，以传递函数作为控制系统的数学模型，以时域分析法、频域分析法为主要分析设计工具，构成了经典控制理论的基本框架。到20世纪60年代初，一套以状态方程作为描述系统的数学模型，以最优控制和卡尔曼滤波为核心的控制系统分析、设计的新原理和方法基本确定，现代控制理论应运而生。控制理论目前还在向更深、更广阔的领域发展，在信息与控制学科研究中注入了蓬勃的生命力，引导人们去探讨更为深刻的运动机理。1.2 Matlab简介Matlab程序设计语言是美国MathWorks公司于20世纪80年代推出的高性能数值计算软件。其功能强大，适用范围广

7、泛，且提供了丰富的库函数（M文件），编程效率高。Matlab无论作为科学研究与工程运算的工具，还是作为计算机辅助的教学工具，都是不可多得的。由于Matlab如此强大的功能，所以它特别适合用来对控制系统进行计算与仿真。系统的设计就是基于Matlab，在正文中再做详细介绍。第二章 控制系统的时域分析与校正2.1 概述2.1.1 时域法的作用与特点时域法是一种直接在时间域中对系统进行分析与校正的方法，它可以提供系统的时间相应的全部信息，具有直观、准确的优点。但在研究系统参数改变引起系统性能指标变化的趋势这一类问题，以及对系统进行校正设计时，时域法不是非常方便的。时域法常用的典型输入信号有单位阶跃信号

8、、单位斜坡信号、等加速度信号、单位脉冲信号。系统能够稳定工作是研究系统动态性能与稳态性能的基本前提。一般情况下，阶跃输入对系统来说是最严峻的工作状态，如果系统在阶跃信号作用下的动态性能能够满足要求，那么在其他形式函数的作用下，其动态性能也是令人满意的。固有关系统的动态性能的指标均是根据系统的单位阶跃响应来定义的。2.1.2 时域性能指标对控制系统的一般要求常归纳为稳、准、快，工程上为了定量评价系统性能好坏，必须给出控制系统的性能指标的准确定义和定量计算方法。稳定是控制系统正常运行的基本条件。系统稳定，其响应工程才能收敛，研究系统的性能（包括动态性能和稳态性能）才有意义。实际物理系统都存在惯性，

9、输出量的改变是与系统所储有的能量有关的。系统所储有的能量的改变需要一个过程。在外作用激励下系统从一种稳定状态转换到另一种状态需要一定的时间。系统的动态性能指标一般有以下几个：延迟时间 阶跃响应第一次达到终值h()的50%所需的时间上升时间 阶跃响应从终值的10%上升到终值的90%所需的时间；对有振荡的系统，也可定义为从0到第一次达到终值所需的时间峰值时间 阶跃响应越过终值h()达到第一个峰值所需的时间调节时间 阶跃响应到达并保持在终值h()的5%误差带内所需的最短时间超调量% 峰值h()超出终值h()的百分比，即% =100%2.2 一阶系统的时间响应及动态性能2.2.1 一阶系统传递函数标准

10、形式及单位阶跃响应一阶系统传递函数的标准形式为(s)=式中，T=1/K称为一阶系统的时间常数，系统特征跟=1/T。2.2.2 一阶系统动态性能分析一阶系统的单位阶跃响应是单调的指数上升曲线,依据调节时间的定义,有h()=1=0.95解得=3T 时间常数是一阶系统的重要特征参数，固可用时间常数T描述一阶系统的响应特性。T越小，系统极点越远离虚轴，过渡过程越快。图2.1给出了一阶系统阶跃响应随时间常数T变化的趋势，及一阶惯性环节。图2.1 一阶系统阶跃响应随T的变化趋势图2.2为一阶积分环节的阶跃响应。图2.2 一阶积分环节的阶跃响应2.3 二阶系统的时间响应及动态性能2.3.1 二阶系统传递函数

11、标准形式及分类常见二阶系统结构图如图2.4（a）所示。其中，K，为环节参数。系统闭环传递函数为：(s)= R(s)C(s)（a）R(s)C(s)(b)图2.4 常见二阶系统结构图为分析方便起见，常将二阶系统结构图表示成如图（b）所示的标准形式，系统闭环传递函数标准形式为(s)= ，分别称为系统的阻尼比和无阻尼自然频率，这两个参数完全决定了二阶系统的响应特性，是二阶系统重要的特征参数。 若系统阻尼比取值范围不同，则特征根形式不同，响应特性也不同，由此可将二阶系统分为以下几类：当01时，系统的时域响应具有非周期特性，称为过阻尼系统当=1时，称为临界阻尼系统当=0时，系统响应为持续的等幅振荡，称为零

12、阻尼系统图2.5 =4rad/s时不同阻尼比下的单位阶跃响应图2.6 =0.7时不同的单位阶跃响应由图2.5可以验证，当01后，随阻尼比的增大，响应越来越迟钝。对于给定的，阻尼比越小，响应的速度越快，如图2.6所示，但阶跃响应的快速性指标调节时间在01时随阻尼比的减小而增大。可以看出，阶跃响应的快速性与、密切相关。对于给定的阻尼比，越大，响应越快，而超调量基本不变。2.3.2 过阻尼二阶系统动态性能指标分析 过阻尼二阶系统单位阶跃响应为=1+(t0)过阻尼二阶系统单位阶跃响应是无振荡的单调上升曲线，令取不同值，可分别求解出相应的无量纲调节时间，如图3.7所示，图中为参变量。图2.7 过阻尼二阶

13、系统%与的关系曲线 欠阻尼二阶系统动态性能指标分析欠阻尼二阶系统单位阶跃响应为=1-sin(t+arctan)如图2.8所示，响应曲线位于两条包络线1/之间，包络线收敛速度取决于（特征根实部之模），响应的阻尼振荡频率取决于（特征根虚部）。响应的初始值h(0)=0，初始斜率h（0）=0，终值h()=1。图2.8 欠阻尼二阶系统单位阶跃响应及包络线ImRe0图2.9 系统极点轨迹对典型欠阻尼二阶系统动态性能而言，当固定，增加（减小）时，系统极点在s平面按图中所示的圆弧轨迹（）移动，对应系统的超调量%减小；同时由于极点远离虚轴，增加，调节时间减小。当固定，增加时，系统极点在s平面按图所示的射线轨迹（

14、）移动，对应的系统超调量%不变；由于极点远离虚轴，增加，调节时间减小。一般实际系统中，T0是系统的固定参数，不能随意改变，而开环增益K是各环节总的传递系数，可以调节。K增大时，系统极点在s平面按图所示的垂直线（）移动，阻尼比变小，超调量%会增加。综上所述，要获得满意的系统动态性能，应当适当的选择参数，使二阶系统的闭环极点位于=45线附近，使系统具有合适的超调量，并根据情况尽量使其远离虚轴，以提高系统的快速性。2.3.4 附加闭环零、极点对系统动态性能的影响对系统附加闭环零点不会影响闭环极点，因而不会影响单位阶跃响应中的各模态，但它会改变单位阶跃响应中各模态的加权系数，由此影响系统的动态性能。附

15、加闭环零点时通过改变单位阶跃响应中各模态的加权系数影响闭环系统动态性能的。若二阶系统闭环传递函数为(s)= ，在其基础上附加闭环零、极点和同时附加闭环零、极点后，得出系统阶跃响应的变化趋势，如图2.10所示由图（a）可以看出，闭环零点的引入带来了系统超调量的增加，使系统的平稳性变差，同时上升时间缩短，响应速度加快。零点值越接近闭环极点实部，对响应的影响就越小。 由图（b）可以看出，当在闭环传递函数极点右侧增加极点时，系统的响应由周期性响应转变为非周期响应，响应平稳性变好，但过渡过程调节时间变长，快速性下降。随着增加的极点越来越靠近虚轴，对系统响应逐渐起主导作用。（a）附加闭环零点对系统阶跃响应

16、的影响（b）附加闭环极点对系统阶跃响应的影响 （c）同时附加闭环零、极点时系统的阶跃响应图2.10 附加零、极点对系统的影响2.4高阶系统的阶跃响应、动态性能及近似高阶系统传递函数一般可以表示为(s)=其单位阶跃响应为(s)=(s)=+=(0),=(s)式中，、分别是与闭环复数极点=处的留数有关的常系数.对上式进行拉式逆变换后，得到高阶系统在零初始条件下的单位阶跃响应为=A+可以看出，高阶系统的单位阶跃响应实际上是由一阶惯性环节、二阶振荡环节的响应叠加而成。当所有极点均具有负实部时，除常数项外，其他各项随时间t而衰减为零。对于系统极点而言，如果负实部远离虚轴，及或值较大，则该极点对应的响应衰减

17、快，对系统整个过渡过程的影响小，因此响应的主要特征取决于靠近虚轴的极点。经验证明，若极点与虚轴的距离大于最靠近虚轴的极点与虚轴距离的5倍以上时，该极点称为远极点，对应的瞬态分量对过渡过程的影响可忽略。对稳定的闭环系统，远离虚轴的极点对应的模态因为收敛较快，只影响阶跃响应的起始段，而距虚轴近的极点对应的模态衰减缓慢，系统动态性能主要取决于这些极点对应的响应分量。此外，各瞬态分量的具体值还与起系数的大小有关。系数大而且衰减慢的分量在瞬态响应中起主要作用。因此，距离虚轴最近而且附近没有零点的极点对系统的动态性能起主导作用，称相应极点为主导极点。图2.11 三阶系统单位阶跃响应曲线图2.12 四阶系统

18、单位阶跃响应曲线一般规定，若极点的实部大于主导极点实部的56倍以上时，则可以忽略相应分量的影响；若两相邻零、极点间的距离比他们本身的模值小一个数量级时，则称该零、极点为“偶极子”，其作用近似抵消，可以忽略相应分量的影响。闭环主导极点常取共轭复数极点，于是相应的系统近似为二阶系统。但应注意的是，应使简化后的系统与高阶系统具有相同的闭环增益，以保证阶跃响应终值相同。第三章 控制系统的频域分析与校正3.1 概述 时域响应法是一种直接法，它以传递函数为系统的数学模型，以拉氏变换为数学工具，直接求出变量的解析解。这种方法虽然直观，分析时域时十分有用，但是方法的应用需要两个前提，一是必须已知控制系统的开环

19、传递函数，另外系统的阶次不能很高。如果系统的开环传递函数未知，或者系统的阶次较高，就不能采用上述方法进行分析。频域分析法不仅是一种通过开环传递函数研究系统闭环性能的分析方法，而且当系统的数学模型未知时，还可以通过实验的方法建立。此外，大量丰富的图型方法使得频域分析法分析高阶系统时，分析的复杂性并不随阶次的增加而显著增加。 当线性系统受正弦信号作用时，其输出特性随正弦信号的频率变化而变化，这种描述系统性能与正弦信号频率之间关系的方法就称为频率特性法，或称频域响应法。之所以将正弦信号作为研究信号，是因为周期信号可以通过傅里叶级数展开成正弦信号的叠加，而非周期信号可将其看做周期T的周期信号。 为了进

20、一步弄清楚频率特性的概念，看一个实验。图3.1所示为一个线性系统，传递函数为G（s），当给系统输入一个正弦信号时，系统的输出稳定后是与输入同频率的正弦信号。如果记输入信号为r（t）=sin(wt)，则稳态输出信号可表示为c(t)=sin(wt+)，与输入信号相比，输出信号的幅值和相位发生了变化。当输入信号幅值不变而频率变化时，输出信号的幅值和相位也会随频率变化而变化，频率特新就是指输出、输入信号幅值比A=/和r(t)r(t)=Arsin(wt)G(s)C(s)c(t)=sin(wt+)0tc(t)0R(s)图3.1正弦信号作用下的系统输入输出相位差随频率变化的规律。定义 正弦信号作用下，线性定

21、常系统输出稳态分量与输入的幅值比和相位差随频率变化的规律称为频率特性，其中幅值比的变化规律A(w)称为频率特性，相位差的变化规律（w）称为相频特性。或者定义为：正弦信号作用下，线性定常系统稳态输出与输入的复数比为系统的频率特性，记为G(j)。3.2 频率特性的表示方法由于频率特性是复变函数，因此既可以表示为实部、虚部的形式：G(j)=U()+jV()也可以将幅频和相频分别表示为A（）=()=G(j)=arctan当以矢量形式表示时，有G(j)=A（）e频率特性是频率的函数，如果在相应的坐标纸上绘制成曲线，就可以直观地分析系统的输出和输入之比相位随频率变化的情况，并且可以通过分析这些曲线的某些特

22、点来判断系统的稳定性与动态品质，并对系统进行分析和综合。通常频率特性采用下面三种曲线形式表示：1、 幅相频率特性当频率由零变化到无穷大时，式表示的矢量末端在复数平面内变化的轨迹为幅相频率特性曲线，也称为极坐标图或乃奎斯特曲线。由式可知，向量G(j)的长度A()等于，由正实轴方向逆时针绕原点转动的角度()等于G(j)。2、 对数频率特性 对数频率特性是由对数幅频特性曲线与对数相频特性曲线两条曲线组成。横坐标采用对数坐标，即频率按对数分度，单位是rads。纵坐标线性分度，幅频值以L（）=20lgA()即dB为单位、相频以度（）或rad为单位，是目前应用较为广泛的一种频率响应图，又称伯德图。采用伯德

23、图表示对数频率特性时，具有以下优点：（1）化乘除运算为加减运算。当系统由多个环节构成时，利用渐进幅频的概念，系统的幅频特可以由各环节的幅频特性折线叠加而成，简洁方便；（2）对数坐标拓宽了图形所能表示的频率范围。（3）如果系统或环节的频率特性互为倒数时，其对数频率特性曲线关于零分贝线对称，相频特性曲线关于零度线对称。（4）将实验获得的频率特性数据绘制成对数频率特性曲线，可以方便地确定系统的出传递函数。3、对数幅相频率特性 在所需要的频率范围内，以频率作为参变量来表示的对数幅值和相位关系的图，称为对数幅相频率特性，也称为尼克尔斯图。3.3 频率特性的性能指标 采用频域方法进行线性控制系统设计时，时

24、域内采用的诸如超调量，调整时间等描述系统性能的指标不能直接使用，需要在频域内定义频域性能指标.1谐振峰值 谐振峰值为幅频特性曲线的A()的最大值。一般说来，的大小表明闭环控制系统相对稳定性的好坏。越大，表明系统对某个频率的正弦信号反映强烈，有共振倾向，系统的平稳性较差，相应阶跃响应的超调量越大。对应的为谐振频率。2带宽 幅频特性下降至零频幅比的70.7，或下降3dB时对应的频率称为带宽（也成为闭环截止频率）。带宽用于衡量控制系统的快速性，带宽越宽，表明系统复现快速变化信号的能力越强，阶跃响应的上升时间和调节时间就越短。带宽是控制系统及控制元件的重要性能指标。3相频宽相频宽为相频衰减90时对应的

25、频率。与一样，也用于衡量系统的快速性。相频宽高，表明输入信号的频率越高，变化较快是输出才能落后90，即系统反映快速，快速性好。4零频幅比A(0)零频幅比A(0)为频率为零时的振幅比。零频信号为直流或常值信号，A(0)=1表明系统阶跃响应的终值等于输入值，即系统的静差为0。A(0)1则表明系统有静差，其与1的差值大小反映了系统的稳态精度，因此A(0)越接近于1，系统的精度越高。谐振峰值小，带宽宽，相频宽高，系统的过渡过程性能好，A(0)越接近于1，系统的精度高，这是频域法分析系统性能的一般准则。()-90A(0)0.707A(0)图3.2 系统幅频、相频特性曲线及性能指标3.4 典型环节的频率特

26、性3.4.1 比例环节 比例环节的传递函数为（s）=K，频率特性为（j）=K于是幅相特性为（j）=K=Ke对数幅频特性与相频特性分别为L（j）=20LgK,()=0图3.3比例环节的伯德图于是伯德图如图3.3所示。 积分环节图3.4积分环节的伯德图 积分环节的传递函数为G=，其频率特性为（j）=于是幅相特性为（j）=对数幅频特性与相频特性为L（j）=-20lg,()=-90于是伯德图如图3.4所示.在伯德图中，积分环节的对数幅频特性曲线每十倍频程衰减20dB,常表示为-20dB/dec. 微分环节 微分环节的传递函数为G（s）=s，频率特性为（j）=j其幅相特性为（j）=j=对数幅频特性与相频

27、特性为L（j）=-20lg，()=90伯德图如图3.5所示。微分环节的对数幅频特性曲线每十倍频程增加20dB，故表示为+20dB/dec，与积分环节的对数幅频特性曲线关于0dB线对称。微分环节的对数相频曲线与积分环节的对数相频曲线关于0线对称。图3.5微分环节的伯德图3.4.4 惯性环节 一阶惯性环节的传递函数为G（s）=，频率特性为G（j）=其幅相特性为G（j）=A()对数幅频特性和相频特性分别为L（）=20lgA()=-20lg()=-arctan(T)由于对数幅频特性曲线L（）为曲线，实际分析中常采用简化的渐近曲线来近似。渐近的原则是以=点为界，即：当，即T，即T1时，有L（）=20lg

28、A()=-20lg-20lg(T)即L（）为lg的线性函数。可以证明，对数相频曲线关于-45线具有奇对称性。以直线代替曲线，给作图带来较大方便。对于惯性环节而言，采用渐进对数幅频曲线代替理论曲线，最大误差点出现在转折频率=处，误差值为3dB。图3.6惯性环节的伯德图3.4.5 一阶微分环节 一阶微分环节的传递函数为G（s）=s+1,频率特性为G（j）=j+1图3.7 一阶微分环节伯德图其幅相特性曲线为G（j）=j+1=对数幅频特性和相频特性分别为L（）=20lgA()=20lg()=-arctan()一阶微分环节与惯性环节的频率特性互为倒数，根据对数频率特性的特点，一阶微分环节与惯性环节的对数

29、幅频特性曲线关于0dB线对称，相频特性曲线关于0线对称。3.4.6 二阶振荡环节二阶振荡环节的传递函数为G(s)=，其频率特性为G（j）=其幅频特性和相频特性分别为A()=()=图3.8 二阶振荡环节的伯德图可以看出，当由0趋于无穷变化时，幅值A()由1衰减为0，相位()由0滞后为-180。由3.8可见，当较小时，由于曲线存在谐振，对数幅频特性渐近线与实际幅频特性曲线存在较大的误差。当渐近误差不超过3dB，可直接使用渐近线近似对数幅频特性，否则应使用准确的对数幅频曲线。3.4.7 高阶系统的伯德图 高阶系统的伯德图，随各指标的变化呈现不稳定的状态，图3.9、3.10、3.11分别表示出了三阶、

30、四阶、五阶的伯德图，以做参考。图3.9 三阶系统的伯德图图3.10 四阶系统的伯德图图3.11 五阶系统的伯德图第四章 结论 对控制系统进行分析，时域响应法是一种直接法，它以传递函数为系统的数学模型，以拉氏变换为数学工具，直接可以求出变量的解析解。这种方法虽然直观，分析时域性能十分有用，但是方法的应用需要两个前提，一是必须已知控制系统的闭环传递函数，另外系统的阶次不能很高。如果系统的开环传递函数未知，或者系统的阶次较高，就需采用频域分析法。频域分析法不仅是一种通过开环传递函数研究系统闭环传递函数性能的分析方法，而且当系统的数学模型未知时，还可以通过实验的方法建立。此外，大量丰富的图形方法使得频

31、域分析法分析高阶系统时，分析的复杂性并不随阶次的增加而显著增加。固在进行控制系统分析时，可以根据实际情况，针对不同数学模型选用最简洁、最合适的方法，从而使用相应的分析方法，达到预期的实验目的。Gm = 2.2222Pm = 38.8302wcg = 0.8164wcp = 0.4530课程设计总结 三周短暂的课程设计已然结束，此次微机测试技术综合训练是测控技术与仪器专业非常重要的一项教育环节，是对我们所学自动控制原理课程的进一步提高与总结。虽然很多知识尚还一知半解，但其中的收获对我而言还是受益匪浅。首先，非常感谢我们的指导老师张立强老师对我们的悉心指导，课设期间张老师始终寸步不离的陪伴我们进行

32、完了所有的实验课程，为我们精心的答疑解惑，使我们的课设进行的非常顺利。在张老师为我们指导的过程中，我被他的学识渊博和人格魅力所深深折服，张老师对自动控制这门课程的理解与研究甚是深入，我们对这门课程的学习只是蜻蜓点水般停留在应试的层面上，经过张老师的指点，使我对控制系统的动态性能与其影响参数的关系有了一个台阶式的飞跃。其次，此次课设的完成得益于我们小组成员的共同合作，大家明确分工，使得整个课设的过程有序而有效率。课程设计的题目是：控制系统的时域/频域分析，我们主要通过Matlab实现了分析的全过程，通过查阅网络资料，借阅相关的书籍，以及向老师请教，掌握了简单的Matlab软件的使用方法。使用Ma

33、tlab做时域分析时，使用了求取连续系统的单位阶跃响应函数step，绘制出了相应的阶跃响应曲线，从得到的曲线中研究了阻尼比和固有频率对系统动态性能的影响，同时还扩展到附加零点、极点对系统动态性能的影响，对时域分析基本有了很全面的认识与体会。在做频域分析时，使用了margin(sys)、mag,pha=bode(num,den,w)等函数，绘制出了相应的Bode图，由Bode图可以编程计算出其相对应的幅值裕度和相角裕度。通过对相同系统下时域法和频域法的对比，最终得出了对各种情况所要使用的分析方法。 邓志杰 2011年1月12日参考文献1 卢京潮自动控制理论西北工业大学出版社，2009年2 张若青

34、控制工程基础及Matlab实践高等教育出版社，2008年3 黄忠霖控制系统Matlab计算及仿真国防工业出版社，2009年4 刘振全Matlab语言与控制系统仿真实训教程化学工业出版社，2009年5 何衍庆控制系统分析、设计和应用Matlab语言的应用化学工业出版社，2004年附录：程序2.1t=0:0.1:10;T=1.0;for i=1:4num=0 1;den=i*T 1;c,x,t=step(num,den,t);plot(t,c,k-);hold on;end;xlabel(t),ylabel(h(t);title(一阶系统阶跃响应随T的变化趋势),grid on;gtext(T=1

35、)gtext(T=2)gtext(T=3)gtext(T=4)程序2.2T=0.1:0.1:1,2;figure(2)hold onfor i=Tnum=1;den=i 0;H=tf(num,den);Step(H)hold on;end;title(一阶积分环节阶跃响应),grid on;程序2.5wn=4;kosai=0.1:0.1:1,2;figure(1)hold onfor i=kosainum=wn.2;den=1,2*i*wn,wn.2;Gk=tf(num,den);sys=feedback(Gk,1,-1);step(sys)endtitle(二阶系统单位阶跃响应);gtext

36、(kosai=0.1)gtext(kosai=0.2)gtext(kosai=0.3)gtext(kosai=0.4)gtext(kosai=0.5)gtext(kosai=0.6)gtext(kosai=0.7)gtext(kosai=0.8)gtext(kosai=0.9)gtext(kosai=1.0)gtext(kosai=2.0)程序2.6wn=4;kosai=0.1:0.1:1,2;figure(1)hold onfor i=kosainum=wn.2;den=1,2*i*wn,wn.2;Gk=tf(num,den);sys=feedback(Gk,1,-1);step(sys)e

37、ndtitle(二阶系统单位阶跃响应);gtext(kosai=0.1)gtext(kosai=0.2)gtext(kosai=0.3)gtext(kosai=0.4)gtext(kosai=0.5)gtext(kosai=0.6)gtext(kosai=0.7)gtext(kosai=0.8)gtext(kosai=0.9)gtext(kosai=1.0)gtext(kosai=2.0)程序2.7Tb=;Ts=;t=0:0.01:50;T2=10; T1=T2:0.1*T2:20*T2; for j=1:length(T1)Tb=Tb T1(j)/T2; num=1/(T1(j)*T2);

38、den=1 (1/T1(j)+1/T2) 1/(T1(j)*T2);y=step(num,den,t);for k=length(y):-1:1; if(abs(y(k)-1)=0.05 Ts=Ts (k*0.01)/T1(j); break; end endendplot(Tb,Ts);grid on;xlim(1 20);xlabel(T1/T2),ylabel(Ts/T1);gtext()title(过阻尼二阶系统的调节时间特性);程序2.8wn=2.5;xi=0.4;t=0:0.05:4;t1=acos(xi)*ones(1,length(t);a1=(1/sqrt(1-xi2);h1

39、=1-a1*exp(-xi*wn*t).*sin(wn*sqrt(1-xi2)*t+t1);bu=a1*exp(-xi*wn*t)+1;b1=2-bu;plot(t,h1,k-,t,bu,-.,t,b1,:,t,ones(size(t),-);legend(阶跃输入,上包线,下包线,阶跃响应);xlabel(omega_nt),ylabel(h(t);grid on;程序2.10(a)t=0:0.1:20;r=ones(size(t);lamd=5,2,1,0.5,0.25;h,l=size(lamd);num0=1;den0=1 1 1;tf0=tf(num0,den0);c=step(tf

40、0,t);plot(t,r,b-,t,c,r-);hold on,for i=1:l tf1=tf(1/lamd(i),1,1);tfa=tf0*tf1;c=step(tfa,t);0 plot(t,c,g-);hold on; y=c; r=1;while(y(r)0.95&y(r)1.05),r=r-1;end; ts=(r+1)*0.1 plot(tr tr,0 0.1,-.);plot(tp tp,0 ymax,-.); plot(0 20,1.05 1.05,-.);plot(0 20,0.95 0.95,-.); plot(ts ts,0 0.95,-.); b=num2str(o

41、vershoot*100);ot=char(overshoot= b %); b=num2str(tp);tpchar=char(tp= b s); b=num2str(tr);trchar=char(tr= b s); b=num2str(ts);tschar=char(ts= b s); text(tp+0.2,ymax,ot); text(tp+0.5,ymax-0.1,tpchar); text(tr+0.2,0.8,trchar); text(ts+0.2,0.9,tschar);endxlabel(t/s),ylabel(h(t);title(附加闭环零点的影响);(b)t=0:0

42、.1:20;r=ones(size(t);im=1;xi=0.5;z=5,2,1,0.5,0.25;h,l=size(z);num0=1;den0=1 1 1;tf0=tf(num0,den0);c=step(tf0,t);figure; plot(t,r,b-,t,c,r-);hold on, for i=1:l tf1=tf(z(i),1,z(i);tfa=tf0*tf1;c=step(tfa,t); plot(t,c,g-);hold on,grid on; y=c; r=1;while(y(r)0.95&y(r)1.05),r=r-1;end; ts=(r+1)*0.1 plot(tr

43、 tr,0 0.1,-.);plot(tp tp,0 ymax,-.); plot(0 20,1.05 1.05,-.);plot(0 20,0.95 0.95,-.); plot(ts ts,0 0.95,-.); b=num2str(overshoot*100);ot=char(overshoot= b %); b=num2str(tp);tpchar=char(tp= b s); b=num2str(tr);trchar=char(tr= b s); b=num2str(ts);tschar=char(ts= b s); text(tp+0.2,ymax,ot); text(tp+0.5

44、,ymax-0.1,tpchar); text(tr+0.2,0.8,trchar); text(ts+0.2,0.9,tschar); end xlabel(t/s),ylabel(h(t);title(附加闭环极点的影响),grid on;程序2.11wn=2;kosai=0.1:0.1:1.0;figure(1)hold onfor i=kosai num=wn.2; den=1,2*wn*i+1,wn.2+2*wn*i,wn.2; sys=tf(num,den);step(sys)endtitle(三阶系统单位阶跃响应曲线);程序2.12wn=4;kosai=0.1:0.2:1,2;f

45、igure(1)hold onfor i=kosainum=wn.2;den=1,2*i*wn,wn.2;Gk=tf(num.2,den.2);sys=feedback(Gk,1,-1);step(sys)endtitle(四阶系统单位阶跃响应曲线);程序3.3num=100;den=1;sys=tf(num,den);margin(sys)grid程序3.4num=100;den=1 0;sys=tf(num,den);margin(sys)Grid程序3.5num=1 0;den=100;sys=tf(num,den);margin(sys)Grid程序3.6惯性环节num=1 0；den=100;sys=tf(num,den);margin(sys)gridw=logspace(-1,1,100);num=1;kosai=0.1:0.3:1.0;for T=kosai den=T,1; sys=tf(num,den);bode(sys)hold onendgridxlabel(角频率(rad/sec);title(惯性环节的伯德图);Gm,