第八章 系统频率响应及其仿真

,第八章 系统频率响应及其仿真,本章主要内容如下：,8.1,频率特性的一般概念,8.2,频率响应的,MATLAB,函数,8.3,系统频域校正,8.4,系统分析图形用户界面,第八章 系统频率响应及其仿真,频率响应分析方法的基本思想是把控制系统中的各个变量看成是由许多不同频率的正弦信号叠加而成的信号；各个变量的运动就是系统对各个不同频率的信号的响应的总和。,这种源于通讯科学的分析方法，于,20,世纪,30,年代引进到控制工程后，立即得到广泛应用。这主要是由于频率响应法具有鲜明的物理意义，能够大大简化复杂机构的动力学分析和设计，更能够启发人们区分影响系统的主要因素和次要因素；其次还可以通过实验方法比较准确地求出系统的数学模型并可减少手工计算量。古典控制理论实际上就是以频率响应法分析可用常系数线性微分方程描述的,SISO,系统。由于许多工业过程都可以近似抽象成线性定常系统，因此频率响应法在控制工程中仍然是一种重要的方法。,第八章 系统频率响应及其仿真,8.1,频率特性的一般概念,8.1.1,频率响应与频率特性,频率响应：,系统对谐波输入的稳态响应。,对于线性系统，当输入为：,x,i,(,t,)=,X,i,sin,t,其稳态输出为同频率的正弦信号：,x,o,(,t,)=,X,o,(,),sin,t,+,(,),8.1,频率特性的一般概念,8.1.1,频率响应与频率特性,频率特性：,是指系统在正弦信号作用下，稳态输出与输入之比对频率的关系特性。可表示为,频率特性还可表示为,因此频率特性还可再分为,实频特性：,U,（,）,虚频特性：,V,（,）,幅频特性：,相频特性：,x,o,(,t,)=,X,o,(,),sin,t,+,(,),第八章 系统频率响应及其仿真,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,图,把频率特性函数,G,(,j,),的角频率,和幅频特性都取对数，则称之为对数幅频特性和对数相频特性，其中：,对数幅频特性,：（单位为分贝,db,）,对数相频特性,：,（单位为度）,其频率轴采用对数分度,lg,。则以,lg,为横坐标，,L,(,G,(,j,),和,(,G,(,j,),为纵坐标绘制的曲线分别称之为对数幅频特性图和对数相频特性图，统称为系统的,Bode,图。,Nyquist,稳定判据引申到对数频率特性中即成为对数判据，因而也可以用,Bode,图分析系统的稳定性。,第八章 系统频率响应及其仿真,8.1,频率特性的一般概念,8.1.2,Nyquist,图与,Bode,图,稳定裕度,利用系统开环频率特性的稳定裕度，可以分析闭环系统的稳定性。稳定裕度又分为幅值裕度和相位裕度。在,Bode,图上表示为：,幅值裕度,(db),：,相位裕度,：,【,说明,】,g,为相位穿越频率,即开环相频特性曲线穿越,180,0,线时的频率,.,c,为幅值穿越频率,即开环幅频特性曲线穿越,0,分贝线时的频率,.,在工程上通常要求,kg,6db,=30,o,60,o,第八章 系统频率响应及其仿真,8.2,频率特性的,MATLAB,函数,8.2.1,频率响应计算函数,MATLAB,提供了用于计算线性时不变系统的频率响应的函数，其调用格式为,h=,freqs,(b,a,w),指定正实角频率向量，返回响应值。,freqs,(b,a,w),绘制对指定角频率向量的幅频和相频特性曲线,.,其中,b,、,a,均为系统传递函数的分子、分母的系数向量。,在返回指令值的指令中，需调用,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),%,不指定频率向量,指定频率向量,不指定频率向量,8.2,频率特性的,MATLAB,函数,8.2.2,频率,特性图示法,Nyquist,图,nyquist,(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,和,w2,之间,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);,%,绘,Nyquist,图,1,title(System,Nyquist,Charts with K=5),figure(2),nyquist,(sys2,w),%,绘,Nyquist,图,2,title(System,Nyquist,Charts with K=30),K=5,K=30,由于系统开环稳定,因此,K=5,时系统是稳定的,(,开环,Nyquist,曲线没有包围,(-1,j0),点,即图中的“,+”,号,);,而,K=30,时系统是不稳定的。,第八章 系统频率响应及其仿真,8.2,频率特性的,MATLAB,函数,8.2.2,频率,特性图示法,Bode,图,bode,(sys),基本调用格式，绘制,Bode,图,bode,(sys,w),指定频率范围，绘制,Bode,图,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 Charts with K=5),grid,figure(2),bode(sys2),%,绘,Bode,图2,title(System Bode Charts with K=30),grid,绘制当,K=5、30,时系统的开环频率特性,Bode,图，并判断系统的稳定性,K=5,K=30,K=5,时,因为,c,g,所以系统闭环不稳定,。,第八章 系统频率响应及其仿真,8.2,频率特性的,MATLAB,函数,8.2.2,频率,特性图示法,计算幅值、相位裕度,margin(sys),为基本调用，用于绘制,Bode,图，并在图中标出幅值裕度和相位裕度。,Gm,Pm,Wcg,Wcp,=margin(sys),返回幅值裕度,Gm,，相位裕度,Pm,，幅值穿越频率,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,时系统的幅值与相位裕度。,由于,K=30,时，相位裕度小于零，幅值裕度小于,1,（即小于零分贝），因此系统不稳定。,kg1=2.4000,r1=19.9079,wg1=4.4721,wc1=2.7992,kg2=0.4000,r2=-18.3711,wg2=4.4721,wc2=6.8885,（,K=5,）,（,K=30,）,第八章 系统频率响应及其仿真,8.3,系统频域校正,8.3.1,频域校正的基本方法,控制系统的频域校正或频域设计是根据给定的频域性能指标,如稳定裕度、频宽、谐振频率等进行控制器的设计，是古典控制理论的一种主要设计方法。控制器的频域校正有以下几种形式,相位滞后校正,相位超前校正,相位滞后,-,超前校正,PID,校正,本节以相位滞后校正为例，介绍应用,MATLAB,进行系统校正的方法。,8.3.2,相位滞后校正,相位滞后校正装置,相位滞后校正可以使系统具有希望的相位裕度和低频增益,(,稳态误差,),，校正装置的传递函数为,相位滞后校正设计思想,先确定增益,K,c,使系统具有希望的稳态精度，再确定校正装置的转折频率使系统具有希望的相位裕度,(,K,c,也可合并到,G,(,s,),中,),。,相位滞后校正设计步骤,根据稳态误差计算,K,c,;,根据,K,c,下原系统开环幅、相频曲线，寻找满足要求相位裕度 所对应的频率作为幅值穿越频率,c,;,根据,c,确定校正环节的转折频率：,即校正环节最大转折频率为幅值穿越频率的,1/10,【,例,】,设校正前系统开环传递函数为 ，设计滞后校正控制器使系统相位裕度为,60,0,，开环增益为,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,计算改变增益后，具有希望相位裕度的系统开环对数幅频穿越频率,wc,.,考虑到校正环节在幅值穿越频率,w,c,处相位滞后的影响，增加,10,0,的预补偿量，则,W=,logspace,(-1,2,100);,%,生成对数频率向量,mag,ph=bode(Gp1,W);,%,产生系统幅频、相频矩阵,mag,=r