控制系统的根轨迹和频域特性分析样本

资料内容仅供您学习参考，如有不当之处，请联系改正或者删除。本科实验报告 课程名称:自动控制原理实验项目:控制系统的根轨迹和频域特性分析实验地点:多学科楼机房专业班级:学号:学生姓名:指导教师:年5月15日一、实验目的和要求:1．学会利用MATLAB绘制系统的根轨迹,并对系统进行分析;2．学会利用MATLAB对系统进行频域特性分析。二、实验内容和原理:1．基于MATLAB的控制系统根轨迹分析1)利用MATLAB绘制系统的根轨迹利用rlocus()函数可绘制出当根轨迹增益k由0至+∝变化时,闭环系统的特征根在s平面变化的轨迹,该函数的调用格式为[r,k]=rlocus(num,den)或[r,k]=rlocus(num,den,k)其中,返回值r为系统的闭环极点,k为相应的增益。rlocus()函数既适用于连续系统,也适用于离散系统。rlocus(num,den)绘制系统根轨迹时,增益k是自动选取的,rlocus(num,den,k)可利用指定的增益k来绘制系统的根轨迹。在不带输出变量引用函数时,rolcus()可在当前图形窗口中绘制出系统的根轨迹图。当带有输出变量引用函数时,可得到根轨迹的位置列向量r及相应的增益k列向量,再利用plot(r,‘x’)可绘制出根轨迹。2)利用MATLAB获得系统的根轨迹增益在系统分析中,常常希望确定根轨迹上某一点处的增益值k,这时可利用MATLAB中的rlocfind()函数,在使用此函数前要首先得到系统的根轨迹,然后再执行如下命令[k,poles]=rlocfind(num,den)或[k,poles]=rlocfind(num,den,p)其中,num和den分别为系统开环传递函数的分子和分母多项式的系数按降幂排列构成的系数向量;poles为所求系统的闭环极点;k为相应的根轨迹增益;p为系统给定的闭环极点。例3-1已知某反馈系统的开环传递函数为试绘制该系统根轨迹,并利用根轨迹分析系统稳定的k值范围。例3-2已知某正反馈系统的开环传递函数如例3-1所示。试绘制系统根轨迹,并计算根轨迹上点-2.3j2.02处的根轨迹增益和此时系统的稳定性。2.基于MATLAB的控制系统频域分析1)利用MATLAB绘制系统的Bode图MATLAB提供的函数bode()能够绘制系统Bode图,该函数的调用格式为[mag,phase,w]=bode(num,den)式中,num和den分别为系统开环传递函数的分子和分母多项式的系数按降幂排列构成的系数行向量;w为频率点构成的向量;mag为系统的幅值;phase为系统的相位。频率向量可由logspace()函数来构成。此函数的调用格式为ω=logspace(m,n,npts)此命令可生成一个以10为底的指数向量(10m10n),点数由npts当bode()函数带输出变量引用函数时,可得系统Bode图相应的幅值mag,相位phase及频率点ω向量,有了这些数据就可利用下面的MATLAB命令绘制系统的Bode图。>>subplot(2,1,1);semilogx(w,20*log10(mag));subplot(2,1,2);semilogx(w,phase)如果只想绘制出系统的Bode图,而对获得幅值和相位的具体数值并不感兴趣,则能够采用如下简单的调用格式bode(num,den)例3-3已知二阶系统的开环传递函数为绘制出当n=3和ζ=0.3时系统的Bode图。2)利用MATLAB绘制系统的Nyquist图MATLAB提供的函数nyquist()能够绘制系统Nyquist图,该函数的调用格式为[Re,Im,w]=nyquist(num,den)其中,num和den分别为系统开环传递函数的分子和分母多项式的系数按降幂排列构成的系数行向量;Re,Im和w分别为频率特性的实部向量、虚部向量和对应的频率向量。有了这些值就可利用命令plot(Re,Im)来直接绘出系统的奈奎斯特图。当然,Nyquist图也可采用与Bode图类似的简单命令来直接绘制。例3-4已知系统的开环传递函数为绘制Nyquist图,并判断系统的稳定性。三、主要仪器设备:安装Windows系统和MATLAB软件的计算机一台。图3-1图3-1负反馈系统的根轨迹(1)例3-1程序:>>num=1;den=conv([1,0],conv([1,1],[1,2]));>>rlocus(num,den);[k,poles]=rlocfind(num,den)执行以上命令,并移动鼠标到根轨迹与虚轴的交点处单击鼠标左键后可得如图3-1所示的根轨迹和如下结果:Selectapointinthegraphicswindowselected_point=0.0000-1.4142ik=6.0000poles=-3.00000.0000+1.4142i0.0000-1.4142i由此可见根轨迹与虚轴交点处的增益k=6,这说明当k<6时系统稳定,当k>6时,系统不稳定;利用rlocfind()函数也可找出根轨迹从实轴上的分离点处的增益k=0.38,这说明当0<k<0.38时,系统为单调衰减稳定,当0.38<k<6时系统为振荡衰减稳定的。(2)例3-2程序:>>num=1;den=conv([1,0],conv([1,1],[1,2]));>>rlocus(-num,den);[k,poles]=rlocfind(-num,den,-2.3+2.02j)执行以上命令可得如下结果和如图3-2所示的根轨迹。k=15.0166poles=-2.3011+2.0195i-2.3011-2.0195i1.6021由此可见,点-2.3j2.02确实为根轨迹上的,且该点处的增益为15.0166,而由于另一个闭环极点位于正实轴上的1.6021点处,故此时系统不稳定。实际上由于系统的一条根轨迹一直位于正实轴上,因此该系统在所有的正值增益k值下均不稳定(3)例3-3程序:>>wn=3;zeta=0.3;w=logspace(-1,2);>>num=wn.^2;den=[12*zeta*wnwn.^2];>>bode(num,den,w);grid;执行后得如图3-4所示Bode图。在曲线窗口中,经过利用鼠标单击曲线上任意一点,能够获得此点所对应的系统在该点的频率与幅值或频率与相位等有关信息。(4)例3-4程序:>>num=0.5;den=[1210.5];nyquist(num,den)执行后可得如图3-5所示的曲线,由于Nyquist曲线没有包围(-1,j0)点,且P＝0,因此由G(s