《系统频率特性》由会员分享,可在线阅读,更多相关《系统频率特性(20页珍藏版)》请在金锄头文库上搜索。
1、第三章系统频率特性系统的时域分析时域分析是分析系统的直接方法,比较直观,但离开计算机仿真,分析高阶系统是困难的。系统频域分析频域分析是工程广为应用的系统分析和综合的间接方法。频率分析不仅可以了解系统频率特性,如截止频率、谐振频率等,而且可以间接了解系统时域特性,如快速性,稳定性等,为分析和设计系统提供更简便更可靠的方法。本章首先阐明频率响应的特点,给出计算频率响应的方法,接着介绍 Nyquist 图和 Bode 图的绘制方法、系统的稳定裕度及系统时域性能指标计算。3.1 频率响应和频率特性频率响应和频率特性3.1.1 一般概念频率响应频率响应是指系统对正弦输入的稳态响应。考虑传递函数为G(s)
2、的线性系统,若输入正弦信号(3.1-1)tXtxiisin)(根据微分方程解的理论,系统的稳态输出仍然为与输入信号系统的稳态输出仍然为与输入信号同频率的正弦信号,只是其幅值和相位发生了变化同频率的正弦信号,只是其幅值和相位发生了变化。输出幅值正比于输入的幅值,而且是输入正弦频率的函数。输出的iX相位与无关,只与输入信号产生一个相位差 ,且也是输入iX信号频率的函数。即线性系统的稳态输出为(3.1-2)(sin)()(00tXtx由此可知,输出信号与输入信号的幅值比是的函数,称为系统的幅频特性幅频特性,记为。输出信号与输入信号相位差也是)(A的函数,称为系统的相频特性相频特性,记为。)(幅频特性
3、:(3.1-3)()()(0 iXXA相频特性:(3.1-4)()()(0i频率特性频率特性是指系统在正弦信号作用下,稳态输出与输入之比对频率的关系特性,可表示为:(3.1-5)()()(0 jXjXjGi频率特性频率特性是传递函数是传递函数的一种特殊形式的一种特殊形式。任何线性)(jG)(sG连续时间系统的频率特性都可由系统传递函数中的 s 以代替j而求得。有三种表示方法:)(jG(3.1-6)()()(jeAjG(3.1-7)()()(jVUjG(3.1-8)(sin)()cos()()(jAAjG式中,实频特性:)(cos)()(AU虚频特性:)()(arctan)()()()()(si
4、n)()(22UVVUAAV一般在分析系统的结构及参数变化对系统性能的影响时,频域分析比时域分析要容易些。根据频率特性,可以较方便地判别系统的稳定性和稳定裕度,并可通过频率特性选择系统参数或对系统进行校正,使系统性能达到预期的性能指标。同时,由频率特性易于选择系统工作频率范围,或根据工作频率要求,设计具有合适的频率特性的系统。频率特性物理意义明确并且可以用实验的方法测定出来。控制系统的频率特性与其动态特性和静态性能之间存在着定性和定量的关系,因此,可以利用图表、曲线和经验公式作为辅助工具来分析和设计系统。3.1.2 频率响应的计算频率响应的计算一、连续时间系统频率响应的计算(3.1-9)01
5、101 1 )()()(asasabsbsb sUsYsGn nn nm mm m LL则系统的频率响应可以由:(3.1-10)01 101 1 )()()()()(ajajabjbjbjGn nn nm mm m LL直接求出。又设已知系统的状态方程模型为:(3.1-11) DUCXYBUAXX.则系统的频率响应可以由下式直接求出:(3.1-12)DBAIjCjG1)()(二、离散时间系统频率响应的计算若离散系统的状态空间模型为(F,G,C,D),则此系统的频率响应为:(3.1-13)DGFIeCTjGTj1)()(如离散系统以传递函数模型表示,将代入,则系统的Tjez频率响应为:(3.2-
6、14)11 2111 21 )()()()()(nTj nnTjnTjmTj mmTjmTj Tj bebebeabebebebeG LL式中,T 为采样周期。应注意应注意,离散时间系统的采样频率,而系统的频Ts2率范围应在之间。20s三、频率响应计算函数频率响应计算函数MATLAB 控制工具箱中,函数 FREQRESP 用于计算 LTI系统的频率响应,它既适用于连续时间系统,也适用于离散时间系统;既适用于 SISO 系统,也适用于 MIMO 系统。函数调用格式为:),(sysfreqrespH 其中,sys 为系统模型; 为指定的实频率向量,单位为 rad/s; 返回值 H 是系统的频率响应
7、。它是一个三维数组。例如,SISO 系统,H(1,1,5)表示频率点所对应响应值;对于 MIMO)5(系统,H(1,2,5)表示第 1 个输出和第 2 个输入之间在频率点)5(的响应值。频率响应 H 为复变量。为了说明函数 FREQRESP 所采用的计算方法,下面程序用两种方法计算一个离散的频率响应:1.采用变换;2.直接Tjez用函数 FREQRESP。例 3-1已知离散系统传递函数为:,9048. 081. 10464. 00478. 0)(2zzzzG采样周期,试计算它的频率响应并绘制其幅频图和相频sTs1 . 0图,amp301.m。图 3-1 系统的频率响应3.2 频率特性图示法频率
8、特性图示法在经典控制论中,常用图示法来描述系统的频率特性,它们是:(1)幅相频特性幅相频特性Nyquist 图图,由表示极坐标上的0的幅值和相角关系。)(jG(2)对数幅相特性对数幅相特性Bode 图图,它由两个图组成:对数幅频特性图和对数相频特性图。纵坐标分别是:幅值,)(lg20)(AL以表示;相角,以度表示。横坐标为频率,采用对数分dB)(度。(3)对数幅相特性对数幅相特性Nichols 图图,它是以为参变量来表示对数幅值和相角关系图。MATLAB 控制工具箱中,有专用的函数可方便地实现这三个图形的绘制。3.2.1 Nyfquist 图的绘制频率特性是频率的复变函数,可以在复平面上用一)
9、(jG个矢量来表示。该矢量的幅值,相角。)()(jGA)()(jG当频率从变化时,矢端的轨迹即为频率特性。0)(jG因此,把频率特性在复平面上用极坐标表示的几何图形,称为频率特性的极坐标图,或称为 Nyquist 图。Nyquist 图不便于分析频率特性中某个环节对频率特性的影响。MATLAB 控制工具箱中有绘制 Nyquist 图的函数NYQUIST,调用格式为:)()Im,Re,), 2, 1(), 2, 1(),()(sysnyquistsysNsyssysnyquistsysNsyssysnyquistsysnyquistsysnyquistLL其中,为系统模型;频率向量;为频率响应实
10、部;频sysReIm率响应虚部。MATLAB 中,频率范围可由两个函数给定:产生频率在和之间个对数分布频率点;),(log21Nspace12N产生频率在和之间个线性分布的频率点;),(21Nlinspace12NN 可以缺省。函数 NYQUIST 用于计算 LTI 系统的 Nyquist 频率响应。调用时,若不包含左边输出变量,函数 NYQUIST 绘制系统的Nyquist 图;调用时,若包含左边输出变量,则不绘图,只输出变量的向量,这常用于分析系统频率特性。同时,MATLAB 控制工具箱中还有绘制 Nichols 图的函数NICHOLS,其调用格式为:)(,), 2, 1(), 2, 1(
11、),()(sysnicholsphasemagsysNsyssysnicholssysNsyssysnicholssysnicholssysnicholsLL函数 nichols(sys)用来计算 LTI 系统的频率响应并绘制Nichols 图,分析系统的开环和闭环特性。例 3-2绘制系统的 Nyquist 图和 Nichols 图,32152)(22sssssGkamp302.m。图 3-2 Nyquist 图由图 3.2 可见,该系统的开环 Nyquist 曲线不包围点,故闭环是稳定的。)0, 1(j3.2.2 Bode 图的绘制Bode 图是由两幅图组成,分别称为对数幅频特性和对数相频特
12、性。它在频率响应法中应用最为广泛。它的横坐标是频率,对数幅频特性的纵坐标是幅值,单位 dB;)/(srad)(lg20jG对数相频特性的纵坐标为,单位 deg。)(Bode 图便于对系统中不同环节的作用以及整个系统进行分析。MATLAB 控制系统工具箱中,用于 Bode 图绘制的函数是BODE。函数函数 BODE 用于计算线性时不变系统(用于计算线性时不变系统(LTI)的频率响应、)的频率响应、幅值和相位,绘制幅值和相位,绘制 Bode 图图,调用方式为:)(,), 2, 1(), 2, 1(),()(sysbodephasemagsysNsyssysbodesysNsyssysbodesys
13、bodesysbodeLL其中,为系统模型;为幅值;为相位;频率范围。sysmagphase函数 BODE 可用于任意 LTI 系统,即单输入单输出(SISO)系统,多输入多输出(MIMO)系统 ,连续时间系统,离散时间系统。用函数绘制系统的 Bode 图时,频率范围将根据系统)(sysbode零极点自动确定。是根据给定的频率范围绘制系统的频率特性),(sysbodesys曲线。是根据给定的频率范围绘制多个系统), 2, 1(sysNsyssysbodeL的频率特性曲线。当函数调用带有左边输出变量时,函数将返回频率响应的幅值 mag,相位 phase 和频率值。例 3-3例 3-1 系统,试绘
14、制其 Bode 图,amp303.m。图 3-3 Bode 图比较图 3-3 和图 3-1 可知,在 MATLAB 中,可用不同方法求得系统的频率响应特性,函数 BODE 完成例 3-1 程序的所有计算。3.3 稳定裕度稳定裕度由 Nyquist 稳定判据可知,若系统开环的 Nyquist 轨迹不包围点,闭环系统是稳定的。当系统开环 Nyquist 轨迹离点)0, 1(j越远,闭环系统的稳定程度越高;开环 Nyquist 轨迹离点)0, 1(j越近,则其闭环系统的稳定程度越低。这称为系统的相对)0, 1(j稳定性。它通过系统开环传递函数对点的靠近程度)(jGk)0, 1(j来表征,定量表示为稳
15、定裕度:幅值裕度 kg 和相位裕度 。一、幅值裕度 kg 的求取在 Nyquist 图上,当为相位交界频率时,开环幅频特性g的倒数,称为系统的幅值 kg,即)()(jHjG(3.3-1)()(1 jHjGkg 显然,在 Nyquist 图上,Nyquist 轨迹与负实轴的交点至原点的距离即为。对于稳定系统,有;对于不稳定系统,kg/11/1kg有。1/1kg在 Bode 图上,幅值裕度改以分贝(dB)表示:(3.3-2)()(lg20)()(1lg20lg20gg ggjHjGjHjgkg对于稳定系统,kg(dB)必在 0dB 线以下,此时称0)(dBkg为正幅值裕度;对于不稳定系统,必在 0dB 线以上,)(dBkg,此时称为负幅值裕度。0)(dBkg二、相位裕度 的求取当为增益交界频率时,相频特性距线的相位差值co180为相位裕度。(3.3-3)(180co式中,一般为负值。)(c在 Nyquist 图上, 为 Nyquist 轨迹和单位圆的交点对负实轴的相位差值。对于稳定系统, 必在 Nyquist 图负实轴以下, 为正相位裕度;对于不稳定系统, 必在