武科大Matlab仿真第八章系统频率响应及其仿真.ppt

《武科大Matlab仿真第八章系统频率响应及其仿真.ppt》由会员分享，可在线阅读，更多相关《武科大Matlab仿真第八章系统频率响应及其仿真.ppt（29页珍藏版）》请在装配图网上搜索。

1、第八章 系统频率响应及其仿真,本章主要内容如下： 8.1 频率特性的一般概念 8.2 频率响应的MATLAB函数 8.3 系统频域校正 8.4 系统分析图形用户界面,第八章 系统频率响应及其仿真,频率响应分析方法的基本思想是把控制系统中的各个变量看成是由许多不同频率的正弦信号叠加而成的信号；各个变量的运动就是系统对各个不同频率的信号的响应的总和。 这种源于通讯科学的分析方法，于20世纪30年代引进到控制工程后，立即得到广泛应用。这主要是由于频率响应法具有鲜明的物理意义，能够大大简化复杂机构的动力学分析和设计，更能够启发人们区分影响系统的主要因素和次要因素；其次还可以通过实验方法比较准确地求出系

2、统的数学模型并可减少手工计算量。古典控制理论实际上就是以频率响应法分析可用常系数线性微分方程描述的SISO系统。由于许多工业过程都可以近似抽象成线性定常系统，因此频率响应法在控制工程中仍然是一种重要的方法。,第八章 系统频率响应及其仿真,8.1 频率特性的一般概念 8.1.1 频率响应与频率特性 频率响应：系统对谐波输入的稳态响应。 对于线性系统，当输入为： xi(t)=Xisint 其稳态输出为同频率的正弦信号：,xo(t)=Xo() sint+(),8.1 频率特性的一般概念 8.1.1 频率响应与频率特性 频率特性：是指系统在正弦信号作用下，稳态输出与输入之比对频率的关系特性。可表示为

3、频率特性还可表示为 因此频率特性还可再分为,第八章 系统频率响应及其仿真,8.1 频率特性的一般概念 8.1.2 Nyquist图与Bode图 Nyquist图 利用封闭的Nyquist轨迹可进行系统稳定性的分析，即Nyquist稳定判据。 Nyquist图不便于分析频率特性中某个环节对频率特性的影响。,频率特性G(j)是频率的复变函数，可以在复平面上用一个矢量来表示。该矢量的幅值为 ，相角为 。当从0变化时，G(j)的矢端轨迹被称之为频率特性的极坐标图或Nyquist图。,第八章 系统频率响应及其仿真,8.1 频率特性的一般概念 8.1.2 Nyquist图与Bode图 Bode图 把频率特

4、性函数G (j)的角频率和幅频特性都取对数，则称之为对数幅频特性和对数相频特性，其中： 对数幅频特性： （单位为分贝db） 对数相频特性： （单位为度） 其频率轴采用对数分度lg。则以lg为横坐标，L(G(j)和(G(j)为纵坐标绘制的曲线分别称之为对数幅频特性图和对数相频特性图，统称为系统的Bode图。,Nyquist稳定判据引申到对数频率特性中即成为对数判据，因而也可以用Bode图分析系统的稳定性。,第八章 系统频率响应及其仿真,8.1 频率特性的一般概念 8.1.2 Nyquist图与Bode图 稳定裕度 利用系统开环频率特性的稳定裕度，可以分析闭环系统的稳定性。稳定裕度又分为幅值裕度和

5、相位裕度。在Bode图上表示为： 幅值裕度(db)： 相位裕度：,【说明】 g为相位穿越频率,即开环相频特性曲线穿越 1800线时的频率. c为幅值穿越频率,即开环幅频特性曲线穿越 0分贝线时的频率. 在工程上通常要求 kg 6db, =30o 60o,第八章 系统频率响应及其仿真,8.2 频率特性的MATLAB函数 8.2.1 频率响应计算函数 MATLAB提供了用于计算线性时不变系统的频率响应的函数，其调用格式为 h = freqs (b, a, w) 指定正实角频率向量，返回响应值。 freqs (b, a, w) 绘制对指定角频率向量的幅频和相频特性曲线. 其中 b、a均为系统传递函数

6、的分子、分母的系数向量。 在返回指令值的指令中，需调用abs（）和angle（）求取幅频和相频特性。 第2种调用可直接绘制系统的幅频和相频特性曲线，其中幅频特性曲线为全对数坐标，而相频特性曲线为半对数坐标，并且可以不指定频率向量。,【例2】频率响应演示：绘制系统 频率响应曲线。,num=11 11;den=1 15 4 0; w=0.05:0.01:0.5*pi; %产生频率向量 freqs(num,den,w) %指定频率向量 freqs(num,den,w) %不指定频率向量,指定频率向量,不指定频率向量,8.2 频率特性的MATLAB函数 8.2.2 频率特性图示法 Nyquist图 n

7、yquist (sys ) 基本调用格式绘制sys的Nyquist图 nyquist (sys, w) 指定频率范围w，绘制sys的Nyquist图 nyquist (sys1, sys2, sysn) 在同一坐标系内绘制多个模型的Nyquist图 nyquist (sys1, sys2 , sysn, w) 在同一坐标系内绘制多个模型对指定频率范围的Nyquist图,【说明】 MATLAB中频率范围w除可直接用冒号生成法生成外，还可由两个函数给定：logspace (w1, w2, N) 产生频率在w1和w2之间N个对数分布频率点；linspace (w1, w2, N) 产生频率在w1和w

8、2之间N个线性分布频率点；N可以省略。 调用nyquist()指令若指定w，则w仍然必须是正实数组,MATLAB将自动绘制与-w对应的Nyquist轨迹。 所绘Nyquist图的横坐标为系统频率响应的实部,纵坐标为虚部.,【例3】系统开环传递函数为 ，绘制当K=5、30时系统的开环频率特性Nyquist图，并判断系统的稳定性。,w=linspace(0.5,5,1000)*pi; sys1=zpk( ,0 -10 -2,100); %建立模型1,K=5 sys2=zpk( ,0 -10 -2,600); %建立模型2,K=30 figure(1), nyquist(sys1,w); %绘Nyq

9、uist图1 title(System Nyquist Charts with K=5) figure(2), nyquist(sys2,w) %绘Nyquist图2 title(System Nyquist Charts with K=30),由于系统开环稳定,因此K=5时系统是稳定的(开环Nyquist曲线没有包围(-1,j0)点,即图中的“+”号);而K=30时系统是不稳定的。,第八章 系统频率响应及其仿真,8.2 频率特性的MATLAB函数 8.2.2 频率特性图示法 Bode图 bode (sys ) 基本调用格式，绘制Bode图 bode (sys, w) 指定频率范围，绘制Bod

10、e图 bode (sys1, sys2,sysn) 在同一图内,绘制多个模型的Bode图 mag,phase,w = bode (sys ) 返回响应的幅值和相位及对应的,不绘制Bode图 bodemag (sys ) 仅绘制幅频bode图 说明 当不指定频率范围时，bode()将根据系统零极点自动确定频率范围。,【例4】系统开环传递函数为,sys1=zpk( ,0 -10 -2,100); %建立模型1，K=5 sys2=zpk( ,0 -10 -2,600); %建立模型2，K=30 figure(1),bode(sys1) %绘Bode图1 title(System Bode Chart

11、s with K=5),grid figure(2),bode(sys2) %绘Bode图2 title(System Bode Charts with K=30),grid,绘制当K=5、30时系统的开环频率特性Bode图，并判断系统的稳定性,K=5时,因为cg,所以系统闭环不稳定。,第八章 系统频率响应及其仿真,8.2 频率特性的MATLAB函数 8.2.2 频率特性图示法 计算幅值、相位裕度 margin (sys ) 为基本调用，用于绘制Bode图，并在图中标出幅值裕度和相位裕度。 Gm,Pm,Wcg,Wcp = margin (sys ) 返回幅值裕度Gm，相位裕度Pm，相位穿越频率

12、Wcg 和幅值穿越频率Wcp，不绘制Bode图。 Gm,Pm,Wcg,Wcp = margin (mag, phase, w) 根据给定幅频向量mag，相频向量phase和对应的频率向量w，计算并返回Gm，Pm，Wcg和Wcp。 说明 Gm=1/|G(jWg)|是Nyquist图对应的幅值裕度，单位不是分贝。,sys1=zpk( ,0 -10 -2,100); %建立模型1 sys2=zpk( ,0 -10 -2,600); %建立模型2 kg1,r1,wg1,wc1=margin(sys1) kg2,r2,wg2,wc2=margin(sys2),【例5】计算例4中K=5和K=30时系统的幅

13、值与相位裕度。,第八章 系统频率响应及其仿真,8.3 系统频域校正 8.3.1 频域校正的基本方法 控制系统的频域校正或频域设计是根据给定的频域性能指标,如稳定裕度、频宽、谐振频率等进行控制器的设计，是古典控制理论的一种主要设计方法。控制器的频域校正有以下几种形式 相位滞后校正 相位超前校正 相位滞后-超前校正 PID校正 本节以相位滞后校正为例，介绍应用MATLAB进行系统校正的方法。,8.3.2 相位滞后校正 相位滞后校正装置 相位滞后校正可以使系统具有希望的相位裕度和低频增益(稳态误差)，校正装置的传递函数为 相位滞后校正设计思想 先确定增益Kc使系统具有希望的稳态精度，再确定校正装置的

14、转折频率使系统具有希望的相位裕度。Kc也可合并到G(s)中。,【例】设校正前系统开环传递函数为 ，设计滞后校正控制器使系统相位裕度为600，开环增益为49。,Gp=tf(1,2,1)*tf(1,0.5,1)*tf(1,0.05,1)*4; %改变增益前系统开环传递函数模型 Gp1=Gp*49/4; %改变增益后系统开环传递函数模型 figure(1),margin(Gp) %改变增益前系统开环Bode图 figure(2),margin(Gp1) %改变增益后系统开环Bode图,计算改变增益前后，系统开环对数频率特性。,K=4,K=49,计算改变增益后，具有希望相位裕度的系统开环对数幅频穿越频

15、率wc. 考虑到校正环节在幅值穿越频率wc处相位滞后的影响，增加100的预补偿量，则,W=logspace(-1,2,100); %生成对数频率向量 mag,ph=bode(Gp1,W); %产生系统幅频、相频矩阵 mag=reshape(mag,100,1); %将幅频矩阵变为幅频向量 ph=reshape(ph,100,1); %将相频矩阵变为相频向量 Wc=interp1(ph,W,-110) %计算相位为-1100时的频率 ,Wc =1.4189 即在幅值穿越频率c处，系统开环频率特性的相位为-1100。,说明 中bode()函数返回的幅频mag、相频ph计算结果均为11100的三维矩

16、阵，需要将其转化换为向量，以便于数值计算 中interp1()为插值函数，用于确定ph中与-1100对应的角频率c.,确定校正环节参数 在c处校正环节的对数幅值应满足： 一阶微分环节转折频率1/T1可根据相位裕度的变化进行调整。,mag110=interp1(ph,mag,-110) ; %计算c处幅值（非分贝值） Beta=mag110 %求取 T1=6/Wc;BT1=Beta*T1; %计算转折频率 Gc=tf(T1,1,BT1,1) %建立校正环节模型,Beta = 13.2494 Transfer function: 4.228 s + 1 - 56.02 s + 1,系统校核,sys

17、=Gc*Gp1 %建立串联校正环节系统开环传递函数模型 figure(3) margin(sys) %计算幅值裕度、相位裕度、相位穿越频率、幅值穿越频率，并绘图,由图可知，校正后系统的幅值裕度为22.8dB，相位裕度为60.80，满足要求。,时间响应比较,figure(4) subplot(2,1,1),step(feedback(Gp1,1,-1) %绘制校正前系统阶跃响应曲线 subplot(2,1,2),step(feedback(sys,1,-1),r) %绘制校正后系统阶跃响应曲线,下图为校正后系统阶跃响应曲线，与上图未校正的系统阶跃响应相比，系统的动态特性明显改善。,第八章 系统频

18、率响应及其仿真,8.4 系统分析图形用户界面 MATLAB控制工具箱提供的线性时不变(LTI)系统仿真的图形用户分析界面LTI Viewer，可更为直观地分析系统的时域、频域响应。其使用也很简单，只需 在指令窗中建立起要分析的系统模型; 在指令窗中键入：ltiview; 即可调出LTI Viewer窗口，进行分析。,【例6】LTI-Viewer使用演示。设单位负反馈系统的开环传递函数,在指令窗中输入： ltiview 即可进入LTI-Viewer可视化仿真环境，如图。,在进入LTI-Viewer后，点击菜单FILE，选择Import项后，弹出一个装入LTI系统的窗口如图所示。该窗口将显示工作空

19、间或指定目录的文件夹内所有的系统模型对象,在LTI浏览器中选择系统“sys”后，就显出系统的阶跃响应图形窗口，在窗口内点击鼠标右键弹出现场菜单，见下图。,菜单的主要功能如下： Plot Types：选择图形类型。可选择Step(阶跃响应,缺省设置) ，Impulse(脉冲响应），Bode图，BodeMag(幅频Bode图)，Nyquist图,Pole / Zero(极点/零点图)等。 Characteristics：可对不同类型响应曲线标出相关特征值。对阶跃响应，可选择表示的特征值如。 Properties：对图形窗口进行编辑，对显示性能参数进行设置。此外，还可以选择菜单EditLinesty

20、le对曲线的线形、颜色、标志等进行选择。,此外还可进行多个图形窗口显示，其操作如下：在LTI Viewe窗口下，选择菜单Edit Plot Configurations后，弹出一个Plot Configurations（图形配置）窗口。该窗口左边显示响应图6种排列形式，通过单选按钮任选其中一种，最多有6种图形显示。该窗口右边显示响应类型，共6组，最多可选择6种（应和所选窗口数对应）。,在该界面上选择四个图形窗口，并使相应窗口分别对应阶跃、脉冲、Bode图和Nyquist图，点击 后 ，即可显示响应图形，如图7-21所示。对图中每个曲线还可分别设置相关选项，如Bode图设置显示稳定裕度，阶跃响应设置显示峰值和峰值时间，用鼠标指向图中的圆点，即可显示出相关数据。,练 习,（1）教材145页第1、2、3、4题 （2）教材166页第3题,