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

精品范文模板 可修改删除撰写人：_日 期：_本科实验报告课程名称： 自动控制原理 实验项目：控制系统的根轨迹和频域特性分析实验地点： 多学科楼机房 专业班级： 学号： 学生姓名： 指导教师： 2012 年 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 已知二阶系统的开环传递函数为绘制出当wn=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 负反馈系统的根轨迹四、实验数据记录和处理（实验程序）（1）例3-1程序：num=1;den=conv(1,0,conv(1,1,1,2);rlocus(num,den);k,poles=rlocfind(num,den)执行以上命令，并移动鼠标到根轨迹与虚轴的交点处单击鼠标左键后可得如图3-1所示的根轨迹和如下结果： Select a point in the graphics windowselected_point = 0.0000 - 1.4142ik = 6.0000poles = -3.0000 0.0000 +1.4142i 0.0000 - 1.4142i 由此可见根轨迹与虚轴交点处的增益k=6，这说明当k6时，系统不稳定；利用rlocfind( )函数也可找出根轨迹从实轴上的分离点处的增益k =0.38, 这说明当0k0.38时，系统为单调衰减稳定，当0.38knum=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.0195i 1.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=1 2*zeta*wn wn.2;bode(num,den,w);grid;执行后得如图3-4所示Bode图。在曲线窗口中，通过利用鼠标单击曲线上任意一点，可以获得此点所对应的系统在该点的频率与幅值或频率与相位等有关信息。（4）例3-4程序：num=0.5;den=1 2 1 0.5; nyquist(num,den) 执行后可得如图3-5所示的曲线，由于Nyquist曲线没有包围(-1,j0)点，且P0，所以由G(s)H(s)构成的单位负反馈闭环系统稳定。在Nyquist曲线窗口中，也可利用鼠标通过单击曲线上任意一点，获得此点所对应的系统的开环频率特性，在该点的实部和虚部及其频率的值五、讨论、心得：通过本实验，让我回顾了matlab的使用，也增长了我对该程序编写的能力，对于控制系统的稳定性和稳态误差也有了更加深刻的理解.第 6 页 共 6 页免责声明：图文来源于网络搜集，版权归原作者所以若侵犯了您的合法权益，请作者与本上传人联系，我们将及时更正删除。