中科大版--现代控制系统(最新版)精品电子教案第八章频率响应法课件.ppt

《中科大版--现代控制系统(最新版)精品电子教案第八章频率响应法课件.ppt》由会员分享，可在线阅读，更多相关《中科大版--现代控制系统(最新版)精品电子教案第八章频率响应法课件.ppt（158页珍藏版）》请在三一办公上搜索。

1、自动控制原理,中国科学技术大学工业自动化研究所,第八章,频率响应法,目录8.1 引论8.2 频率响应图8.3 频率响应测量8.4 频域性能指标8.5 对数幅相图8.6 设计实例8.7 应用控制设计软件研究频率响应法8.8 循序渐进设计案例：磁盘驱动器读取系统8.9 总结,习题,8.18.3节：E8.5、E8.6、E8.14、8.1(a)(d）、P8.6、P8.17、P8.19、P8.23、AP8.18.48.8节：E8.8、P8.7、P8.8、P8.12、P8.14、P8.15、P8.24、P6.12、P6.13、AP8.5、DP8.3、DP8.5作业说明：题目中所有绘图（极坐标图，Bode图

2、，对数幅相图等）必须手工绘制，不能用Matlab绘制,频率响应frequency response：系统对正弦输入信号的稳态响应线性定常系统在正弦输入信号作用下，稳态时的输出信号、系统内部各处的信号，都是正弦信号，信号频率与输入信号相同，幅值、相位与输入信号有差别对于系统：Y(s)=T(s)R(s)，r(t)=Asint，有：,8.1 引论,8.1 引论,假定极点-pi各不相同，展开成部分分式：,系统T(s)稳定，所有极点-pi都是负实部，则：,8.1 引论,正弦输入时，稳定的线性定常系统稳态输出：,输入信号频率、振幅A确定后，稳态时的输出信号只与T(j)的幅值、相位有关,8.1 引论,频率响

3、应法的优点：各种频率、振幅的正弦测试信号容易获得，用实验方法容易确定系统的频率响应通过测量频率响应可以获得传递函数在频域中设计系统，可以有效地控制系统的带宽、抑制噪声和扰动令s=j，可以由传递函数T(s)得到T(j)用图形可以很方便地表达T(j)的幅值、相角，并深刻洞察控制系统分析、设计频率响应法的缺点：频率响应特性与时域的瞬态响应特性之间缺乏准确的定量关系,Laplace变换对Laplace transform pair：,8.1 引论,Fourier变换对Fourier transform pair：,已知f(t)的Laplace变换F(s)，令s=j，得f(t)的Fourier变换F()

4、,用Laplace变换可以得到系统传递函数T(s)，可以对T(s)在s平面的零极点分布进行分析用频率响应法可以得到传递函数T(j)，关心的是系统的幅值和相位特性研究闭环系统频率响应时，输入r(t)在频域中具有Fourier变换：,8.1 引论,对闭环系统输出Y(s)=T(s)R(s)，令s=j，可得单回路反馈系统的输出频率响应：,采用傅氏反变换，可得输出瞬态响应：,8.1 引论,但是，除了最简单的系统，傅氏反变换积分很难求解,频率响应图：系统的频率响应特性，可以在不同的坐标系中，用不同的曲线、图形表示三种常用的频率响应图：Nyquist图（极坐标图、幅相特性图）、Bode图（对数坐标图、对数频

5、率特性图）、Nichols图（对数幅相图）系统传递函数G(s)的频域表达：,8.2 频率响应图,例8.1 RC滤波器的频率响应RC滤波器的传递函数：,8.2 频率响应图,正弦稳态传递函数：,8.2 频率响应图,可以由下式绘制极坐标图：,极坐标平面：,8.2 频率响应图,第一步：确定=0和时的R()、X(),第二步：确定取其它值的图形,RC滤波器的极坐标图,8.2 频率响应图,=1时实部、虚部的幅值相等，()=-45,可以根据频率特性的幅值、相角表达来绘制极坐标图：,8.2 频率响应图,例8.2 传递函数的极坐标图,8.2 频率响应图,8.2 频率响应图,采用G(j)的实部、虚部绘制极坐标图：,

6、8.2 频率响应图,在S平面上沿虚轴s=j，在每个频率，用图解法求取向量G(j)，绘制极坐标图：,8.2 频率响应图,极坐标图的优点、缺点：采用Nyquist稳定性判据，便于分析稳定性频率响应计算繁琐单个极点、零点的影响难以明显区分增加新的极点或零点，频率响应要重新计算从试验得到的极坐标图，难以获得传递函数对数坐标图logarithmic plots，波德图、波特图、Bode图Bode plots，为纪念H.W.Bode频域的传递函数：,8.2 频率响应图,对数幅频特性图：对数角频率、对数增益值,相频特性图：对数角频率、相位(),例8.3 RC滤波器的Bode图,8.2 频率响应图,对数增益值

7、：,8.2 频率响应图,低频段（1/）对数增益值：,高频段（1/）对数增益值：,=1/处，对数增益值：,频率=1/称为转折频率break frequency，转角频率corner frequency，三分贝频率相位特性：,一阶惯性环节对数幅频特性图,一阶惯性环节相频特性图,频率轴采用线性刻度很不方便，采用对数刻度就很方便。水平轴取为log，当1/时：,8.2 频率响应图,一阶环节对数增益曲线近似为直线十倍频程decade：两个频率相差10倍，两个频率间的间隔称为十倍频程。2=101，从1到2的频率范围称十倍频程。对数增益差：,一阶环节高频渐近线的斜率为-20dB/十倍频程,Bode图采用半对数

8、坐标纸，横轴为角频率，采用以10为底的对数刻度，单位：弧度/秒rad/s；对数幅频特性图的纵轴为对数增益值，采用线性刻度，单位：分贝dB；相频特性图的纵轴为角度，采用线性刻度，单位：度2倍频程octave：两个频率相差2倍，两个频率间的间隔称为2倍频程，倍频程。2=21，从1到2的频率范围称倍频程，对数增益差：,8.2 频率响应图,一阶环节高频渐近线的斜率为-6dB/倍频程,一阶惯性环节对数幅频特性图,一阶惯性环节对数幅频特性渐近线,Bode图的主要优点：利用对数增益值，将传递函数中各个因子增益值的相乘，转化为各个因子对数增益值的相加传递函数的一般形式：,8.2 频率响应图,系统拥有Q个零点、

9、原点有N个极点、M个实极点、R对共轭复极点,频率特性的一般形式：,8.2 频率响应图,G(j)的对数增益值：,将各个因子的对数幅频特性图迭加，得到G(j)的对数幅频特性图G(j)的相角为：,8.2 频率响应图,将各个因子的相频特性图迭加，得到G(j)的相频特性图为了充分利用迭加原理，发挥Bode图的优势，必须首先掌握各种基本因子的Bode图,传递函数中共有四种因子：常数增益Kb原点的极点（零点）实极点（零点）复共轭极点（零点）,8.2 频率响应图,确定这四种因子的对数幅频特性图、相频特性图，对于一般形式的传递函数，先绘制各有关因子的Bode图，迭加后就得到系统的Bode图绘制对数幅频特性时，可

10、以先绘制渐近线，得到近似的对数幅频特性，然后获取若干关键频率的真实值，得到比较精确的对数幅频特性图Bode图的测量法,常数增益Kb,8.2 频率响应图,常数增益的对数幅频特性是一条水平线，相频特性也是一条水平线,原点的极点,8.2 频率响应图,原点的极点对数幅频特性为斜率-20dB/十倍频程的斜线，相频特性为-90的水平线原点的多重极点,原点的多重极点对数幅频特性为斜率-20NdB/十倍频程的斜线，相频特性为-90N水平线,原点的零点,8.2 频率响应图,原点的零点对数幅频特性为斜率+20dB/十倍频程的斜线，相频特性为+90的水平线,(j)N的Bode图,实极点,8.2 频率响应图,渐近线是

11、对真实对数幅频特性的近似，最大误差发生在转折频率，为-3dB,一阶惯性环节对数幅频特性,一阶惯性环节相频特性,对相频特性采用三段直线近似：,8.2 频率响应图,实际相频特性与近似相频特性在转折频率处交叉，在其他频率点上，最大误差不超过6,一阶惯性环节对数幅频特性、相频特性的真实值与近似值,实零点,8.2 频率响应图,渐近线是对真实对数幅频特性的近似，最大误差发生在转折频率，为+3dB相频特性也可以用三段直线近似,8.2 频率响应图,共轭复极点,实际的对数幅频特性和渐近线之间的误差与阻尼比有关，0.707时误差不容忽略,8.2 频率响应图,共轭复极点对数幅频特性图,共轭复极点相频特性图,谐振频率

12、resonant frequencyr处，频率响应达到最大值，称为谐振峰值Mp阻尼比0，谐振频率r无阻尼自然振荡频率n在谐振频率r处，增益达到最大值：,8.2 频率响应图,频率响应的最大值称为谐振峰值：,传递函数基本因子的渐近线,传递函数基本因子的渐近线,传递函数基本因子的渐近线,传递函数基本因子的渐近线,传递函数基本因子的渐近线,绘制Bode图S平面法：沿虚轴(s=j)取各个频率点，绘制零点、极点到该频率点的向量，确定向量的幅值、相角，可得频率响应曲线考虑具有一对共轭复极点的二阶系统：,8.2 频率响应图,令s=j，频率特性函数为：,在频率点上，G(j)的幅值、相角为：,8.2 频率响应图,

13、8.2 频率响应图,=0时的幅值、相位,8.2 频率响应图,=r时幅值、相位,8.2 频率响应图,=d时的幅值、相位,用s平面法得到的共轭复极点的Bode图,8.2 频率响应图,例8.4 双T网络Bode图双T网络twin-T network也称带阻滤波器、陷波滤波器，抑制某个频段的信号，而让其它频段的信号通过在控制中，常用于抑制机械共振,8.2 频率响应图,双T网络,双T网络可视为一个低通滤波器、一个高通滤波器并联，低通滤波器由两个电阻R、一个电容2C构成T形网络，高通滤波器由两个电容C、一个电阻R/2构成T形网络等效电路：,8.2 频率响应图,双T网络的传递函数：,8.2 频率响应图,8.

14、2 频率响应图,双T网络的频率特性：,8.2 频率响应图,8.2 频率响应图,在具有相同的幅频特性的系统中，在全部频率范围内，相角都最小的系统称为最小相位系统minimum phase system，此外则称为非最小相位系统nonminimum phase system 线性定常系统的所有极点、零点都位于闭左半S平面，并且没有纯滞后，则为最小相位系统系统有极点或者零点位于右半S平面，或者有纯滞后，则为非最小相位系统右半S平面的极点称为不稳定极点右半S平面的零点称为非最小相位零点幅频特性相同时，与非最小相位系统相比，在所有频率上最小相位系统的相角都较小,8.2 频率响应图,最小相位系统：线性定常

15、系统，如果系统本身及其逆系统都是因果的、稳定的，则为最小相位系统没有纯滞后的稳定系统，将系统的零点换成关于虚轴的对称点，就得到一个新系统。与原系统的传递函数相比，新系统的传递函数幅频特性没有改变，相频特性改变了。比较这两个系统传递函数的相移，若系统的零点全在左半S平面，则系统的相角在所有频率点上都较小具有相同幅频特性的系统，最小相位系统的相移范围小，而非最小相位系统的相移范围大最小相位系统的幅频特性、相频特性具有唯一对应关系，只需知道一个，另一个也就知道,8.2 频率响应图,以下两个系统，极点完全一致，零点不同：,8.2 频率响应图,零极点的模式表明，系统幅频特性相同，但相频特性不同,可见，二

16、个系统的幅频特性相同，而相频特性不同,G1(s)相移范围小于80，G2(s)相移范围达180,8.2 频率响应图,最小相位和非最小相位传递函数的相位特性,全通网络all-pass network：对称的格点网络,8.2 频率响应图,由零极点分布图可知，零极点的分布是对称的,8.2 频率响应图,例8.5 绘制Bode图包含所有基本因子的传递函数：,8.2 频率响应图,按特征频率递增的顺序，排列各基本因子项：常数增益K=5原点处的极点1/j=2处的极点1/(1+j0.5)=10处的零点(1+j0.1)=n=50处的共轭复极点：1/(1+j0.6(/50)+(j/50)2)，为0.3,绘制各基本因子

17、幅频特性渐近线,8.2 频率响应图,常数增益：20 log5=14dB原点处的极点：斜率-20dB/十倍频程，在=1处通过0dB线,8.2 频率响应图,=2处的极点：转折频率=2以下，幅频特性渐近线为0dB线，转折频率以上，幅频特性渐近线为斜率-20dB/十倍频程的斜线=10处的零点：转折频率=10以下，幅频特性渐近线为0dB线，转折频率以上，幅频特性渐近线为斜率+20dB/十倍频程的斜线=n=50处的共轭复极点：幅频特性渐近线斜率为-40dB/十倍频程，转折频率=n=50。由于阻尼比很小=0.3，在转折频率附近，实际的幅频特性与渐近线相比有较大偏差，必须修正将各基本因子渐近线迭加，得到近似幅

18、频特性,8.2 频率响应图,按特征频率递增的顺序迭加渐近线，就能直接绘制出G(j)的幅频特性渐近线第一段渐近线是(j)-1产生的斜率为-20dB/dec的直线，该直线在=1处与14dB线相交受=2处极点的影响，在2后，G(j)的渐近线变成斜率为-40dB/dec的直线受=10处零点的影响，在10后，G(j)的渐近线变成斜率为-20dB/dec的直线受n=50共轭复极点影响，在50后，G(j)的渐近线最终变成-60dB/dec的直线再根据各因子渐近线与实际曲线的偏差，对渐近线分段修正，得到G(j)的精确幅频特性,8.2 频率响应图,先绘制G(j)各零点、极点因子的相频特性渐近线，迭加可以得到G(

19、j)的近似相频特性常数增益：相角为0原点处的极点：相角为-90=2处的极点：相角近似为三段直线，=2处相角为-45=10处的零点：相角近似为三段直线，=10处相角为+45n=50的共轭复极点：采用准确的相频特性各基本因子相频特性渐近线迭加，得到相频特性近似曲线，可以确定实际相频特性的基本形状，估计系统在某些特殊频率点的相角,8.2 频率响应图,8.2 频率响应图,在实际使用Bode图时，通常先用G(j)的幅频和相频近似曲线确定重要的频率或频段然后再在较小的范围内，用公式准确计算系统的实际相角和幅值分析、设计控制系统时，最关心对数幅频特性为0dB、相频特性为-180的频率点从近似相频特性可知，=

20、46时，相角为-180=46时，实际相角为：,8.2 频率响应图,从图中可以精确得到某些特殊频率的幅值和相角特性,8.2 频率响应图,用正弦信号测量控制系统的开环频率响应，获得增益、相角随频率变化的曲线，利用这2条曲线，导出开环频率特性GH(j)测量控制系统的闭环频率响应T(j)，导出实际的闭环传递函数波形分析仪wave analyzer，在输入正弦波的频率改变时，测量幅值、相位的变化传递函数分析仪transfer function analyzer，测量开环、闭环传递函数典型的信号分析仪可完成从直流到100kHz范围的频率响应，内置分析和建模功能，可以由测量得到的频率响应推断极点、零点，或对

21、用户提供的模型给出频率响应,8.3 频率响应测量,8.3 频率响应测量,扫频仪 在电子测量中，经常遇到对网络的阻抗特性和传输特性进行测量的问题，其中传输特性包括增益和衰减特性、幅频特性、相频特性等。用来测量前述特性的仪器称为频率特性测试仪，简称扫频仪。它为被测网络的调整，校准及故障的排除提供了极大的方便。扫频仪一般由扫描锯齿波发生器、扫频信号发生器、宽带放大器、频标信号发生器、X轴放大、Y轴放大、显示设备、面板键盘、多路输出电源等部分组成。其基本工作过程是通过电源变压器将50Hz市电降压后送入扫描锯齿波发生器，,8.3 频率响应测量,就形成了锯齿波，这个锯齿波一方面控制扫频信号发生器，对扫频信

22、号进行调频，另一方面该锯齿波送到X轴偏转放大器放大后，去控制示波器X轴偏转板，使电子束产生水平扫描。由于这个锯齿波同时控制电子束水平扫描和扫频振荡器，因此电子束在示波管荧光屏上的每一水平位置对应于某一瞬时频率。从左向右频率逐渐增高，并且是线性变化的。扫频信号发生器产生的扫频信号送到宽带放大器放大后，送入衰减器，然后输出扫频信号到被测电路。为了消除扫频信号的寄生调幅，宽带放大器增设自动增益控制器AGC。宽带放大器输出的扫频信号送到频标混频器，在,8.3 频率响应测量,频标混频器中与1MHz和10MHz或50MHz晶振信号或外频标信号进行混频。产生的频标信号送入Y轴偏转放大器放大后输出给示波管的Y

23、轴偏转板。扫频信号通过被测电路后，经过Y轴电位器、衰减器、放大器放大后送到示波管的Y轴偏转板，得被测电路的幅频特性曲线。,8.3 频率响应测量,超低频频率特性测试仪PLT-A 一种多功能智能化仪表，能够产生较宽频带（0.000110000rad/s）的高精度超低频正弦信号，并能实现待测频率点的自动对数间距输入，能一次性输入最多达400个待测频率点，对频率特性既可单点测量、连续测量，还可以进行全自动测量、记录和运算，用LED显示或通过PP40彩描器打印测量结果，绘出各种频率特性曲线，具有全自动测量、坐标自动转换、彩色打印等功能。,8.3 频率响应测量,超低频频率特性分析仪L420526频率范围：

24、1mHz-1000Hz（0.05%）频率分辨率：1/4000扫频：对数、线性、上扫、下扫、保持、停止、任意变换停止控制：0、90、180、270谐波失真：1%分析器输入幅度范围：0V-300V、1mHz-1000Hz输入阻抗：1M测量误差：直角坐标A+jB：0.5%（满度）1个字；极坐标R：0.5%（满度）1个字，：0.11个字,传递函数未辨识的系统的Bode图,8.3 频率响应测量,8.3 频率响应测量,传递函数未辨识的系统的Bode图,从Bode图确定传递函数：电阻、电容组成的稳定的电路从100增加到1000，幅值下降约-20dB/十倍频程；=300处，相角-45，幅值-3dB。=300为

25、系统的转折频率，p1=300为极点相角突然变化+180，在n=2450穿越0；幅值的斜率由-20dB/十倍频程变成+20dB/十倍频程。n=2450处有一对共轭复零点当超过50000，幅值变成0dB；当=20000，幅值比高频0dB渐近线低-3dB，相角+45。可以推断p2=20000为系统第二个极点在转角频率n=2450处，幅值渐近线与最低幅值之差为10dB，可得=0.16,8.3 频率响应测量,传递函数：,8.3 频率响应测量,T形桥接网络bridged-T network传递函数,谐振峰值Mp：幅频特性的峰值（最大值）谐振频率r：幅频特性达到谐振峰值的频率带宽B：幅频特性从低频值下降3d

26、B的频率，称截止频率，0B的频率范围称带宽，或-3dB带宽。增益降低到0.7倍低频增益的频率考虑二阶系统，闭环传递函数：,8.4 频域性能指标,带宽B衡量系统如实复现输入信号的能力B增大，阶跃响应上升时间减少带宽与无阻尼自然振荡频率n、阻尼比有关,8.4 频域性能指标,二阶系统幅频特性,谐振峰值Mp表征系统的相对稳定性谐振峰值Mp增大，表明阻尼比减小，阶跃响应超调量增大对于二阶系统：,8.4 频域性能指标,二阶系统单位阶跃响应为：,8.4 频域性能指标,阻尼比给定，n越大，系统响应速度越快为满足时域性能指标，要求频域性能指标：谐振峰值相对较小，例如Mp1.5带宽较大，使系统时间常数=1/(n)

27、足够小以上分析的基础是，系统的主导极点是一对共轭复极点,闭环系统稳态误差由开环系统幅频特性在低频段的位置和形状决定N型开环系统传递函数为：,8.4 频域性能指标,K为相应的稳态误差系数,8.4 频域性能指标,0型系统：,幅频特性图低频段为水平线稳态位置误差系数Kp=K，在幅频特性图中表现为低频增益，影响低频段的位置幅频特性低频段或其延长线与纵轴交点为：,8.4 频域性能指标,1型系统：,幅频特性图低频段斜率为-20dB/十倍频程稳态速度误差系数Kv=K，影响低频段的位置对开环1型系统：,速度误差系数Kv=5,幅频特性低频段或其延长线与横轴交点为：,8.4 频域性能指标,1型系统幅频特性低频段或

28、其延长线与0dB轴交点的频率，近似为K=Kv,对数幅相图log-magnitude-phase diagram：在一定频率范围内，对数幅值随相角变化的曲线。横坐标：相角，单位：；纵坐标：对数幅值，单位：dB。Nichols图Nichols图的形状能够提供重要的系统信息，特别是当相角接近-180、对数幅值接近0dB时Nichols图对于研究闭环系统相对稳定性很有用考虑以下系统：,8.5 对数幅相图,例8.6 雕刻机控制系统：雕刻机采用两个驱动电机和相联的导向螺杆，在x方向定位雕刻针；另一个独立的电机在y、z两个方向定位雕刻针。本例研究x轴定位控制系统,8.6 设计实例,设计目标：用频率响应法选择

29、合适的K值，使对阶跃指令的时间响应达到要求随意选择初始值K=2，绘制系统的开环、闭环Bode图，用闭环Bode图来预估系统时域响应，并检验预估结果与实际结果的偏差闭环传递函数：,8.6 设计实例,开环系统固有部分频率响应G(j),闭环系统幅频特性图,闭环系统相频特性图,闭环幅频特性谐振频率、谐振峰值为：,8.6 设计实例,设系统有二阶主导极点，由谐振峰值、归一化谐振频率与阻尼比的关系，Mp=1.78时，阻尼比=0.29；r/n=0.91，而r=0.8，则自然振荡频率n=0.88。近似的闭环传递函数：,超调量37%，2%准则的调整时间15.7秒。实际超调量34%，调整时间17秒。近似模型合理如果

30、要减少超调量，可以降低K值,例8.7 六足机器人的单足控制Ambler是卡内基梅隆大学正在开发的六足机器人本例研究单足的定位控制,8.6 设计实例,执行机构和单足的传递函数，输入为给执行机构的电压指令，输出为单足在竖直方向位移：,8.6 设计实例,控制目标：控制机器足的位置，并在存在测量噪声时保持机器足的位置,被控变量：机器足的位置Y(s)期望机器足尽快移动到指令位置，且超调量要最小。起始目标是设计一个能移动的系统（尽管慢），初始设计时，控制系统的带宽较低设计指标：DS1：闭环带宽大于1HzDS2：阶跃输入百分比超调量小于15%DS3：阶跃输入稳态跟踪误差为零设计指标DS1、DS2保证系统跟踪

31、性能设计指标DS3很容易满足，系统固有部分（执行机构/足）传递函数是1型系统，保证对阶跃输入的稳态误差为零,8.6 设计实例,确立控制目标,确定被控变量,给出设计指标,确立系统结构,建立过程、执行机构、传感器的模型,确定控制器并选取需要调整的关键参数,优化关键参数并分析系统性能,系统性能不满足设计指标则重新选择系统结构,系统性能满足设计指标结束设计,控制机器足位置，有干扰时保持期望位置,机器足位置,如图,见方程,六足机器人控制系统设计过程中的要素,设计指标：DS1：b1HzDS2：P.O.15DS3：阶跃输入稳态误差为零,见方程,8.6 设计实例,采用控制设计软件,选择控制器：,8.6 设计实

32、例,可调参数c使控制器增加一个自由度，更加灵活选择关键调整参数：K，a，b，c也可以选择如下控制器：,闭环系统响应主要由主导极点位置决定。由性能指标确定主导极点位置，设计控制器参数，配置期望主导极点，忽略非主导极点响应,开环传递函数：,8.6 设计实例,闭环系统为四阶系统，选择参数K，a，b，c，使其中两个极点为主导极点，并配置到能够满足性能指标的位置,闭环传递函数：,期望的特征多项式：,8.6 设计实例,选非主导极点为重实极点，特征多项式为：,8.6 设计实例,对应项系数相等，K，a，b，c为：,二阶系统带宽B与自然频率n有近似关系：,8.6 设计实例,由性能指标DS1：,由百分比超调量指标

33、DS2：,尽管对调整时间没有指标，但一般希望调整时间尽可能短。所以选择=0.52,自然振荡频率、谐振峰值、谐振频率为：,8.6 设计实例,取=0.52，n=5.11，=12，可得控制器：,百分比超调量P.O.=14%，调整时间0.96秒,闭环系统阶跃响应,闭环系统幅频特性,可见，闭环系统带宽B=27.2rad/s=4.33Hz，满足性能指标DS1，由于是四阶系统，所以B大于1Hz。B 越大，调节时间越短实际的谐振峰值MP=1.21，预期值是1.125当输入为正弦信号，输入信号频率增加，输出信号幅值减小：输入为=1rad/s的正弦信号，稳态输出幅值近似为1输入为=500rad/s的正弦信号，稳态

34、输出幅值小于0.005,8.6 设计实例,输入为角频率=1rad/s正弦信号，闭环系统输出响应,输入为角频率=500rad/s正弦信号，闭环系统输出响应,本节介绍如何使用Matlab来绘制Bode图，再次讨论频率性能指标与时域性能指标的联系，并举例说明频域内的控制系统设计主要的Matlab函数:bode 和 logspacebode：用于绘制Bode图logspace:用于生成频率点数据向量，这些频率点数据是根据频率对数的相等间距生成的,8.7 应用控制设计软件研究频率响应法,Bode函数BODE(SYS)，绘制线性时不变系统sys的Bode图，频率区间及点数由Matlab自动选择BODE(S

35、YS,WMIN,WMAX)，绘制频率在wmin,wmax范围内的bode图BODE(SYS,W)（W：保存Matlab要计算幅值和相角的频率向量，可以由logspace函数指定），绘制给定W向量的Bode图MAG,PHASE=BODE(SYS,W)，将特定频率处的幅值和相角分别保存在mag和pahse向量中MAG,PHASE,W=BODE(SYS)，同时也将计算了幅值和相位的频率点保存在W向量中,8.7 应用控制设计软件研究频率响应法,对系统:,绘制Bode图的Matlab代码：g=tf(1 0.1 7.5,1 0.12 9 0 0);bode(g),8.7 应用控制设计软件研究频率响应法,需

36、要更宽的频率范围，例如0.1-100g=tf(1 0.1 7.5,1 0.12 9 0 0);bode(g,0.1,100),8.7 应用控制设计软件研究频率响应法,Logspace函数y=logspace(a,b)，产生10a-10b之间对数间距相等的50个点的行向量yy=logspace(a,b,n)，产生10a-10b之间对数间距相等的n个点的行向量y,8.7 应用控制设计软件研究频率响应法,Bode图，考虑传递函数：生成Bode图的Matlab代码：%num=5*0.1 1;f1=1 0;f2=0.5 1;f3=1/2500 6/50 1;den=conv(f1,conv(f2,f3)

37、;sys=tf(num,den);%bode(sys),8.7 应用控制设计软件研究频率响应法,Bode函数说明mag,phase,w=bode(sys,w),8.7 应用控制设计软件研究频率响应法,使用logspace函数的Bode图绘制W=logspace(-1,3,200);Bode(sys,w),8.7 应用控制设计软件研究频率响应法,通常控制系统设计的目标是使系统满足给定的时域性能指标，因此，在频域内设计控制系统时，应弄清楚频率响应和时域响应的相互关系产生2阶系统不同阻尼系数的频率响应曲线簇zeta=0.15:0.05:0.7;wr_over_wn=sqrt(1-2*zeta.2);

38、Mp=(2*zeta.*sqrt(1-zeta.2).(-1);Subplot(211),plot(zeta,Mp),gridxlabel(zeta),ylabel(M_pomega);Subplot(212),plot(zeta,wr_over_wn),gridxlabel(zeta),ylabel(omega_r/omega_n);,8.7 应用控制设计软件研究频率响应法,8.7 应用控制设计软件研究频率响应法,雕刻机系统设计目标：选择增益K的合适取值，使得闭环系统阶跃响应的各项指标保持在允许的范围内频率设计的流程图见下页图所示。基本步骤：首先确定增益K的初值，计算其闭环传递函数根据闭环传

39、递函数，绘制其Bode图，计算频率性能指标MP和r的值根据MP 和r与 和n的关系，计算和n利用和n,估计时域性能指标MP和Ts的值判断性能指标是否满足设计目标，如满足，则设计结束，否则重新给K赋初值，重复上述设计过程，直到满足设计目标,8.7 应用控制设计软件研究频率响应法,雕刻机系统的频域设计流程图,设计程序代码，根据增益K来计算时域性能指标num=K;den=1 3 2 K;sys=tf(num,den);w=logspace(-1,1,400);mag,phase,w=bode(sys,w);mp,l=max(mag);wr=w(l);zeta=sqrt(0.5*(1-sqrt(1-1

40、/mp2);wn=wr/sqrt(1-2*zeta2);ts=4/zeta/wnpo=100*exp(-zeta*pi/sqrt(1-zeta2)K=2时，计算结果：Ts=15.7，P.O.=39%如果不满足设计目标，重新选择K值再计算,8.7 应用控制设计软件研究频率响应法,满足要求系统的单位阶跃响应K=2;num=K;den=1 3 2 K;sys=tf(num,den);t=0:0.01:20;y=step(sys,t);plot(t,y);gridxlabel(Time(s),ylabel(y(t),8.7 应用控制设计软件研究频率响应法,磁盘驱动器采用弹性簧片固定读磁头基座。悬挂了读

41、磁头的弹性簧片可以用质量块M、弹簧k和摩擦b来模拟，由读磁头臂施加给弹性簧片的力为u(t),8.8 循序渐进设计案例：磁盘驱动器读取系统,弹性簧片与读磁头的弹簧、质量、摩擦力模型,弹性簧片与读磁头的弹簧-质量-阻尼传递函数:,8.8 循序渐进设计案例：磁盘驱动器读取系统,对典型的弹性簧片和读磁头：,包括弹性簧片读磁头基座的磁盘驱动器读磁头定位控制,从开环幅频特性图可见，在机械共振频率n，实际的幅频特性比渐近线高出约10dB。实践中应尽量避免激励机械共振频率,开环系统幅频特性图,闭环系统幅频特性图,闭环系统带宽B=2000rad/s如果用一对共轭复极点近似磁盘驱动器读取系统的闭环特性，则有：,8.8 循序渐进设计案例：磁盘驱动器读取系统,只要控制器增益K400，就可以将机械共振频率排除出系统带宽,闭环系统阶跃响应2%允许误差的调整时间为：,用频率响应特性表达反馈控制系统系统的频率响应是系统对正弦输入信号的稳态响应极坐标图、Bode图、对数幅相图绘制Bode图渐近线时可以采用迭加方法闭环系统频域性能指标：MP，r，B开环系统幅频特性图低频段的形状与位置和闭环系统稳态误差、稳态误差系数的关系15种典型传递函数Bode图,8.9 总结,THE END,