matlab7.0基础与编程入门课件第十一讲MATLAB在控制系统中的应用幻灯片.pptVIP

二 控制系统的时域分析 一个动态系统的性能常用典型输入作用下的响应来描述。响应是指零初始值条件下某种典型的输入函数作用下对象的响应，控制系统常用的输入函数为单位阶跃函数和脉冲激励函数（即冲激函数）。在MATLAB的控制系统工具箱中提供了求取这两种输入下系统响应的函数。 1、时域分析的一般方法 求取系统单位阶跃响应：step() 求取系统的冲激响应：impulse() 1、step()函数的用法 y=step(num,den,t)：其中num和den分别为系统传递函数描述中的分子和分母多项式系数，t为选定的仿真时间向量，一般可以由t=0:step:end等步长地产生出来。该函数返回值y为系统在仿真时刻各个输出所组成的矩阵。 [y,x,t]=step(A,B,C,D,iu)：其中A,B,C,D为系统的状态空间描述矩阵，iu用来指明输入变量的序号。x为系统返回的状态轨迹。 如果对具体的响应值不感兴趣，而只想绘制系统的阶跃响应曲线，可调用以下的格式： step(num,den)；step(num,den,t)；step(A,B,C,D,iu,t)；step(A,B,C,D,iu)； 线性系统的稳态值可以通过函数dcgain()来求取，其调用格式为：dc=dcgain(num,den)或dc=dcgain(a,b,c,d) [y,x,t]=step(num,den)：此时时间向量t由系统模型的特性自动生成, 状态变量x返回为空矩阵。 exp4_3.m（传递函数模型），exp4_3_.m（状态空间模型） 2、impulse()函数的用法 求取脉冲激励响应的调用方法与step()函数基本一致。 y=impulse(num,den,t)；[y,x,t]=impulse(num,den)；[y,x,t]=impulse(A,B,C,D,iu,t) impulse(num,den)；impulse(num,den,t) impulse(A,B,C,D,iu)；impulse(A,B,C,D,iu,t) 仿真时间t的选择： 对于典型二阶系统根据其响应时间的估算公式 可以确定。 对于高阶系统往往其响应时间很难估计，一般采用试探的方法，把t选大一些，看看响应曲线的结果，最后再确定其合适的仿真时间。 一般来说，先不指定仿真时间，由MATLAB自己确定，然后根据结果，最后确定合适的仿真时间。 在指定仿真时间时，步长的不同会影响到输出曲线的光滑程度，一般不易取太大。 例exp4_6_.m 2、常用时域分析函数 时间响应探究系统对输入和扰动在时域内的瞬态行为，系统特征如：上升时间、调节时间、超调量和稳态误差都能从时间响应上反映出来。MATLAB除了提供前面介绍的对系统阶跃响应、冲激响应等进行仿真的函数外，还提供了大量对控制系统进行时域分析的函数，如： covar：连续系统对白噪声的方差响应 initial：连续系统的零输入响应 lsim：连续系统对任意输入的响应 对于离散系统只需在连续系统对应函数前加d就可以，如dstep，dimpulse等。 它们的调用格式与step、impulse类似，可以通过help命令来察看自学。 三、时域分析应用实例 MATLAB的step()和impulse()函数本身可以处理多输入多输出的情况，因此编写MATLAB程序并不因为系统输入输出的增加而变得复杂。 σ超调量 三 控制系统的频域分析 频率响应是指系统对正弦输入信号的稳态响应，从频率响应中可以得出带宽、增益、转折频率、闭环稳定性等系统特征。 频率特性是指系统在正弦信号作用下，稳态输出与输入之比对频率的关系特性。频率特性函数与传递函数有直接的关系，记为： 1、频域分析的一般方法 求取系统对数频率特性图（波特图）：bode() 求取系统奈奎斯特图（幅相曲线图或极坐标图）：nyquist() 频域分析法是应用频率特性研究控制系统的一种典型方法。采用这种方法可直观地表达出系统的频率特性，分析方法比较简单，物理概念比较明确，对于诸如防止结构谐振、抑制噪声、改善系统稳定性和暂态性能等问题，都可以从系统的频率特性上明确地看出其物理实质和解决途经。通常将频率特性用曲线的形式进行表示，包括对数频率特性曲线和幅相频率特性曲线简称幅相曲线，MATLAB提供了绘制这两种曲线的函数。 1、对数频率特性图（波特图） exp4_10.m exp4_10_.m 对数频率特性图包括了对数幅频特性图和对数相频特性图。横坐标为频率w，采用对数分度，单位为弧度/秒；纵坐标均匀分度，分别为幅值函数20lgA(w)，以dB表示；相角，以度表示。 MATLAB提供了函数bode()来绘制系统的波特图，其用法如下： bode(a