控制系统的频域分析

控制系统的频域分析第1页，课件共52页，创作于2023年2月主要内容14.1控制系统的频域分析14.1.1频率特性概述14.1.2频率特性的不同表示方法14.1.3MATLAB频域分析的相关函数14.1.4MATLAB频域分析实例14.2基于频域法的控制系统稳定性能析14.2.1频域法的稳定性判定和稳定裕度概述14.2.2频域法的控制系统稳定判定相关函数14.2.3MATLAB频域法稳定性判定实例第2页，课件共52页，创作于2023年2月14.1控制系统的频域分析

14.1.1频率特性概述频域法是一种工程上广为采用的分析和综合系统的间接方法。它是一种图解分析法，所依据的是频率特性数学模型，对系统性能如稳定性、快速性和准确性进行分析。频域法因弥补了时域法的不足、使用方便、适用范围广且数学模型容易获得而得到了广泛的应用。第3页，课件共52页，创作于2023年2月14.1.2频率特性的不同表示方法1.对数坐标图对数坐标图即Bode图，由对数幅频特性和对数相频特性曲线2张图组成。第4页，课件共52页，创作于2023年2月

对数幅频特性幅度的对数值与频率的关系曲线；对数相频特性是频率特性的相角与频率的关系曲线。对数幅频特性的纵轴为第5页，课件共52页，创作于2023年2月采用线性分度；横坐标为角频率，采用对数分度。对数相频特性的纵轴为，单位为度，采用线性分度；横坐标为角频率，也采用对数分度。横坐标采用对数分度，扩展了其表示的频率范围。第6页，课件共52页，创作于2023年2月2.极坐标图

极坐标图即Nyquist曲线。系统的频率特性表示为：频率特性是输入信号频率的复变函数，当频率从连续变化时，端点的极坐标轨迹。MATLAB在绘制第7页，课件共52页，创作于2023年2月

Nyquist曲线时频率是从连续变化的。而在自动控制原理的教材中一般只绘制频率从部分曲线。可以分析得出，曲线在范围与内，是以横轴为镜像的。第8页，课件共52页，创作于2023年2月14.1.3MATLAB频域分析的相关函数bode(G)bode(G,w)bode(G1,'r--',G2,'gx',…)[mag,phase,w]=bode(G)[mag,phase]=bode(G,w)绘制系统Bode图。系统自动选取频率范围绘制系统Bode图。由用户指定选取频率范围同时绘制多系统Bode图。图形属性参数可选返回系统Bode图相应的幅值、相位和频率向量。可使用magdb=20*log10(mag)将幅值转换为分贝值返回系统Bode图与指定w相应的幅值、相位。可使用magdb=20*log10(mag)将幅值转换为分贝值表14.1频域分析的相关函数用法及功能说明第9页，课件共52页，创作于2023年2月nyquist(sys)nyquist(sys,w)nyquist(G1,'r--',G2,'gx',…)[re,im,w]=nyquist(sys)[re,im]=nyquist(sys,w)绘制系统Nyquist图。系统自动选取频率范围绘制系统Nyquist图。由用户指定选取频率范围同时绘制多系统Nyquist图。图形属性参数可选返回系统Nyquist图相应的实部、虚部和频率向量。返回系统Nyquist图与指定w相应的实部、虚部。第10页，课件共52页，创作于2023年2月14.1.4MATLAB频域分析实例例1：系统的开环传递函数为绘制系统的Bode图。第11页，课件共52页，创作于2023年2月>>s=tf('s');>>G=1000*(s+1)/(s*(s+2)*(s^2+17*s+4000))

Transferfunction:1000s+1000--------------------------------s^4+19s^3+4034s^2+8000s

>>bode(G)>>grid第12页，课件共52页，创作于2023年2月图14.1例1系统Bode图第13页，课件共52页，创作于2023年2月图14.2例1系统的Bode图(指定频率范围)第14页，课件共52页，创作于2023年2月例2：系统的开环传递函数为绘制系统的Bode图。

第15页，课件共52页，创作于2023年2月num=5;den=conv([12],[121]);w=logspace(-2,3,100);%指定频率范围[mag,phase,w]=bode(num,den,w);magdB=20*log10(mag);%进行幅值的单位转换subplot(2,1,1);semilogx(w,magdB);%绘制对数幅频特性图grid;第16页，课件共52页，创作于2023年2月title('系统Bode图');xlabel('Frequency(rad/sec)');ylabel('GaindB');subplot(2,1,2);semilogx(w,phase);%绘制对数相频特性图grid;xlabel('Frequency(rad/sec)');ylabel('Phasedeg');第17页，课件共52页，创作于2023年2月图14.3例2系统的Bode图第18页，课件共52页，创作于2023年2月例3：系统的开环传递函数为绘制K取不同值时系统的Bode图。%K分别取10，50,1000k=[105001000];forii=1:3第19页，课件共52页，创作于2023年2月G(ii)=tf(k(ii),[110500]);endbode(G(1),'r:',G(2),'b--',G(3))title('系统K/(s^2+10s+500)Bode图,K=10,500,1000','fontsize',16);grid第20页，课件共52页，创作于2023年2月图14.4例3K分别取10，50,1000的系统Bode图第21页，课件共52页，创作于2023年2月例4：单位负反馈系统的开环传递函数为

绘制系统Nyquist曲线。第22页，课件共52页，创作于2023年2月num=[202010];den=conv([110],[110]);nyquist(num,den)第23页，课件共52页，创作于2023年2月图14.5例4系统的Nyquist曲线第24页，课件共52页，创作于2023年2月

对于图14.5，如果想要看清某部分细节，也可通过设置坐标范围进行局部放大，从而得到更清晰的局部图像，如图14.6。num=[202010];den=conv([110],[110]);Nyquist(num,den)axis([-22-55])第25页，课件共52页，创作于2023年2月图14.6例4局部放大的系统Nyquist曲线第26页，课件共52页，创作于2023年2月同样，还可通过设置范围得到局部的Nyquist曲线。如只绘制系统位于的Nyquist曲线，如图14.7。num

=[202010];den=conv([110],[110]);w=0.1:0.1:100;[re,im]=nyquist(num,den,w);第27页，课件共52页，创作于2023年2月plot(re,im)axis([-22-55]);grid;title('系统(20s^2+20s+10)/[(s^2+s)(s+10)]Nyquist图(\omega>0)','fontsize',12);第28页，课件共52页，创作于2023年2月)图14.7例4系统Nyquist曲线()第29页，课件共52页，创作于2023年2月14.2基于频域法的控制系统稳定性能分析

14.2.1频域法的稳定性判定和稳定裕度概述1.Nyquist稳定判据如果开环模型含有m个正极点，则单位负反馈下单变量闭环系统稳定的充要条件是开环系统的Nyquist图逆时针围绕(-1,j0)点m周第30页，课件共52页，创作于2023年2月2.系统相对稳定性的判定（稳定裕度）

系统的稳定性固然重要，但它不是唯一刻画系统性能的准则，因为有的系统即使稳定，使其动态性能表现为很强的振荡，也是没有用途的。因为这样的系统如果出现小的变化就可能使系统不稳定。此时还应该考虑对频率响应裕度的定量分析，使系统具有一定的稳定裕度。第31页，课件共52页，创作于2023年2月

相角稳定裕度系统极坐标图上模值等于1的矢量与负轴的夹角：

相角稳定裕度表示了系统在临界稳定状态时

，系统所允许的最大相位滞后。

第32页，课件共52页，创作于2023年2月幅值稳定裕度

系统极坐标图上与负实轴交点（）的模值倒数：在对数坐标图上，采用表示的分贝值，有第33页，课件共52页，创作于2023年2月幅值稳定裕度表示了系统在临界稳定状态时，系统增益所允许的最大增大倍数。3.闭环系统频率特性性能指标

通常，描述闭环系统频率特性的性能指标主要有谐振峰值、谐振频率、带宽和带宽频率。其中：第34页，课件共52页，创作于2023年2月谐振峰值指系统闭环频率特性幅值的最大值。谐振频率指系统闭环频率特性幅值出现最大值时的频率。带宽频率指当系统的幅频特性下降到时所对应的频率。系统带宽指频率范围。第35页，课件共52页，创作于2023年2月14.2.2基于频域法的控制系统稳定判定相关MATLAB函数

除上节给出的函数可用于绘制频率响应图形并用于判定系统稳定性之外，MATLAB还提供了相关函数直接用于进一步判定系统的稳定程度，见表14.2。第36页，课件共52页，创作于2023年2月margin(G)[Gm,Pm,Wg,Wp]=margin(G)[Gm,Pm,Wg,Wp]=margin(mag,phase,w)绘制系统Bode图,带有裕量及相应频率显示给出系统相对稳定参数。分别为幅值裕度、相角裕度、幅值穿越频率、相角穿越频率给出系统相对稳定参数。由Bode函数得到的幅值、相角和频率向量计算。返回参数分别为幅值裕度、相角裕度、幅值穿越频率、相角穿越频率S=allmargin(G)返回相对稳定参数组成的结构体。包含幅值裕度、相角裕度及其相应频率，时滞幅值裕度和频率，是否稳定的标识符表14.2基于频域法的控制系统稳定判定相关MATLAB函数第37页，课件共52页，创作于2023年2月%用于求取谐振峰值,谐振频率，带宽和带宽频率

[mag,phase,w]=bode(H,w);[Mp,k]=max(mag);resonantPeak=20*log10(Mp);resonantFreq=w(k);n=1;while20*log10(mag(n))>=-3n=n+1;endbandwidth=w(n);%根据闭环Bode图求参数%求取谐振峰值%进行谐振峰值单位转换%求取谐振频率%求取带宽和带宽频率

对系统闭环频率特性的求取，MATLAB没有提供相应的函数。可根据其定义，编写程序求取：第38页，课件共52页，创作于2023年2月14.2.3MATLAB频域法稳定性判定实例例7：系统开环传递函数为绘制其极坐标图，并判定系统稳定性。第39页，课件共52页，创作于2023年2月num=10;den=conv([15],[1-1]);G=tf(num,den);nyquist(G)第40页，课件共52页，创作于2023年2月例7系统的Nyquist曲线，结合开环极点和绕行圈数，可以看出系统稳定

例7系统的闭环零极点分布图也可以得出相同的结论

第41页，课件共52页，创作于2023年2月例8：系统开环传递函数为绘制其极坐标图，并判定系统稳定性。第42页，课件共52页，创作于2023年2月num=5;den=conv([12],[125]);G=tf(num,den);nyquist(G)第43页，课件共52页，创作于2023年2月图14.15例8系统的Nyquist曲线图14.16例8系统的阶跃响应曲线第44页，课件共52页，创作于2023年2月

例9：分别判定系统

和的稳定性。如果系统稳定，进一步给出系统相对稳定参数。第45页，课件共52页，创作于2023年2月num1=5;den1=conv([12],[150]);G1=tf(num1,den1);margin(G1)figure(2)num2=200;den2=conv([12],[150]);G2=tf(num2,den2);margin(G2)第46页，课件共52页，创作于2023年2月图14.17例9系统的Bode图第47页，课件共52页，创作于2023年2月图14.18例9系统的Bode图第48页，课件共52页，创作于2023年2月例10：单位负反馈系统的开环传递函数为绘制闭环系统的Bode图，并给出闭环频率特性性能指标：谐振峰值、谐振频率和系统带宽。第49页，课件共52页，创作于2023年2月s=tf('s');G=1/s/(0.5*s+1)/(s+1);H=feedback(G,1);%闭环系统传递函数w=logspace(-1,1);