资源描述：
《实验三控制系统的频率特性》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、实验三控制系统的频率特性1.实验目的1)掌握运用MATLAB软件绘制控制系统波特图的方法；2)掌握MATLAB软件绘制奈奎斯特图的方法；3)利用波特图和奈奎斯特图对控制系统性能进行分析。2.实验仪器1)PC机一台2)MATLAB软件3.实验原理1-奈華斯特稳定判据及稳定裕量(1)奈氏(Nyquist)判据：反馈控制系统稳定的充要条件是奈氏曲线逆时针包围临界点的圈数R等于开环传递函数右半s平面的极点数P,即R=P;否则闭环系统不稳定，闭环止实部特征根个数Z可按下式确定Z=P-R=P-2N(2)稳定
2、裕量利用轨迹上两个特殊点的位置来度量相角裕度和増益裕度。其中与单位圆的交点处的频率为叭(截上频率)：与负实轴的交点频率为©(穿越频率)。则相角裕度:Y=180°+ZG(皿)丹伽)=180°+0(0.)增益裕度:h1=—-—
3、Gg)"®)
4、心)(对数形式：h=-201g
5、G(j7yv)//(j7yv)
6、=-201gA(^v)2.对数频率稳定判据将系统开环频率特性曲线分为幅频特性和相频特性，分别画在两个坐标上，横轴都用频率0,纵轴一个用对数幅值和相角，这两条曲线画成的图就是Bode图，即对数频率特性
7、图。因为Bode图与奈氏图有一一对应关系，因此，奈氏稳定判据就可描述为基于Bode图的对数频率稳定判据：(1)开环系统稳定，即开环系统没有极点在止右半根平而，如杲其对数幅频曲线大于OdB的区域内，相频曲线对-180°线正负穿越次数相等，那么闭坏系统就是稳定的，否则是不稳定的。(2)开环系统不稳定，有P个极点在止右半平面，如果其对数幅频曲线人于OdB的区域内，相频曲线对-180°线正穿越次数人于负穿越次数P/2,闭环系统就是稳定的，否则是不稳定的。2.利用MATLAB绘制Nyquist图和Bode
8、图MATLAB控制系统工具箱提供了许多函数，用來绘制系统的Nyquist曲线、Bode图以及Nichols图。并町以进行增益裕度和相角裕度的分析。相关常用函数如表所示。表3-1MATLAB频率特性函数函数名函功能描述allmargin计算所有的交叉频率和稳定裕量bode计算并绘制BODE图bodemag计算并绘制BODE幅频特性图evalfr计算系统单频率点处的频率响应freqresp计算系统的频率响应interp在FRD模型频率点间插入频率响应数据linspace生成平均频率间隔的向量logs
9、pace生成平均对数频率间隔的向量margin计算增益裕度和相角裕度ngridNichols网格线nichols绘制Nichols曲线图nyquist绘制Nyquist曲线图3.实验内容绘制下例各控制系统波特图和奈奎斯特图。①系统开环传递函数为：G($)=;绘制系统Nyquist图和波特图。s(s+1)系统开坏传递函数为:G(X22(W)52(0.045+1)(0.45+1)绘制Nyquist图和波特图，并求出系统的相角稳定裕量和幅值稳定裕量。③系统开环传递函数为:500(0,01675+1)s
10、(0.05s+l)(0.0025s+1)(0.001$+1)绘制系统波特图，并求出系统的相角稳定裕呈和幅值稳定裕量。④已知控制系统开坏传递函数为：G°($)=(s+l)(0・5s+1)(02+1)定判据判定开环放人系数K为10和50时闭环系统的稳定性。2.实验步骤(1)画图程序：①k=1；z=[];p=[0,-1]；G=zpk(z,p,k);figure(1);nyquist(G);figure(2);bode(G)mp)①pnlc6e乏(6①P)①sx<100BodeDiagram0图3-2系
11、统1的伯德图①num=[22];den=conv(conv([1,00],[0.04,1]),[0.4,1]);G=tf(num,den);bode(G);[Gm,Pm,Wx]=margin(G);请改此语句，使其显示出剪切频率Wc,并记录。Gm=7.6999(单位？)Pm=20.9451Wx=5.9161②num=500*[0.0167,1];den1=conv([1,0],[0.05,1]);den2=conv([0.0025,1],[0.001,1]);den=conv(den1,den2
12、);G0=tf(num,den);w=logspace(0,4,50);bode(G0,w);margin(G0);由程序运行结果和图示知道，幅值穿越频率w=?rad/s,相角稳定裕量=?;相角穿越频率w=?,幅值稳定裕量k=?，即—dbo①⑴当KT0时G0=tf(10,conv([1,1],conv([0.5,1],[0.2,1])));nyquist(G0);(2)当K=50时G0=tf(50,conv([1,1],conv([0.5,1],[0.2,1])));nyquist(GO);上面
