MATLAB的根轨迹分析法及重点习题上课讲义.doc

Good is good, but better carries it.精益求精，善益求善。MATLAB的根轨迹分析法及重点习题-4.1某系统的结构如题4-1图所示，试求单位阶跃响应的调节时间ts，若要求ts=0.1秒，系统的反馈系数应调整为多少？解:(1)由系统结构图可知系统闭环传递函数为:在单位阶跃函数作用下系统输出为：为求系统单位阶跃响应，对C(s)进行拉斯反变换：根据定义调节时间等于响应曲线进入5%误差带，并保持在此误差带内所需要的最短时间，且根据响应系统单位阶跃响应的函数表达式可以看出系统单位阶跃响应的稳态值为,因此:(2)若要求ts=0.1秒,则有:即:若要求调节时间缩小为0.1秒,则需将反馈环节的反馈系数调整为0.3。4.2已知二阶系统的阶跃响应曲线如题4.2图所示，该系统为单位负反馈系统，试确定其开环传递函数。解：根据系统阶跃响应曲线可以看出：峰值时间，超调量；根据课本中对典型二阶系统暂态性能指标的推导计算可知：结合本题已知阶跃响应曲线可知：由式(2)可知：将带入式(1)中可得：回顾题意对于典型二阶系统其闭环传递函数为，且系统为单位负反馈系统，所以系统开环传递函数和闭环传递函数之间满足如下关系：4.3单位反馈控制系统开环传递函数为,若,试求(1)动态性能指标.(2)欲使,当T不变时，K应取何值。解：(1)对于单位反馈控制系统，已知开环传递函数可求出其闭环传递函数，并将其化为标准形式为：所以根据动态性能指标的计算公式将上述两参数带入后可得：(2)由于T=0.25s,所以可知:,将阻尼比带入超调量的计算公式中:4.4设控制系统如题4-4图所示，其中图(a)为无速度反馈系统，图(b)为带速度反馈系统，试确定系统阻尼比为0.5时Kt的值，并比较图(a)和图(b)系统阶跃响应的动态性能指标。解:(1)根据系统结构图可求得两系统的闭环传递函数为:(2)根据已经求得的两系统的阻尼比和无阻尼自然振荡角频率可分别计算两系统的动态性能指标:经对比可看出:采用速度反馈的b系统虽然上升时间和峰值时间稍有延长,但超调量存在明显下降,系统振荡剧烈程度下降,另外调节时间也显著降低,即使说系统能够在较快的时间内达到稳定,系统动态性能得到了提高。4.5某系统结构如题4-5图所示，试判断系统的稳定性。解：根据系统结构图可利用梅森公式求解其传递函数，结构图中前向通道有一条，回路有两个，且两回路相关，因此有：因此可得系统特征方程为：列写其劳斯表为：根据劳斯判据可知，劳斯表第一列系数符号均未发生变化，因此系统稳定。4.6已知系统特征方程如下，用劳斯判据判断系统的稳定性，如不稳定求在s右半平面的根数及虚根值。另外用MATLAB软件直接求其特征跟加以验证。(1)通过分析劳斯表中第一列系数可知并没有符号变化，所以不存在位于s右半平面的特征根，另外由于系统劳斯表中出现全零行,所以系统不稳定，存在对称于原点的根为用MATLAB软件中函数roots求特征方程根可得，系统特征方程根为如下所示，进一步验证了上述求解结果的正确性。(2)通过分析劳斯表中第一列系数可知并没有符号变化，所以不存在位于s右半平面的特征根，另外由于系统劳斯表中出现全零行,所以系统不稳定，存在对称于原点的根为用MATLAB软件中函数roots求特征方程根可得，系统特征方程根为如下所示，进一步验证了上述求解结果的正确性。(3)通过分析劳斯表中第一列系数符号变化两次，所以有两个位于s右半平面的特征根，另外由于系统劳斯表中出现全零行,所以系统不稳定，存在对称于原点的根为用MATLAB软件中函数roots求特征方程根可得，系统特征方程根为如下所示，进一步验证了上述求解结果的正确性。(4)通过分析劳斯表中第一列系数符号变化1次，所以有1个位于s右半平面的特征根，另外由于系统劳斯表中出现全零行,所以系统不稳定，存在对称于原点的根为用MATLAB软件中函数roots求特征方程根可得，系统特征方程根为如下所示，进一步验证了上述求解结果的正确性。4.7已知单位反馈系统的开环传递函数为，试确定系统稳定时的K值范围。解：根据已知单位反馈系统开环传递函数可知系统特征方程为：，对系统列写劳斯表可得：欲使系统稳定则需满足：汇总得：4.8已知系统的特征方程为，试确定系统稳定时的K值范围，若要求闭环系统极点均位于s=-1垂线之左，K值该如何调整。解：(1)对系统列写劳斯表为：欲使系统稳定则需满足：，即：。(2)将带入系统原特征方程中得：欲使系统稳定则需满足：；即4.9已知系统稳定，求的系统稳态误差。解：由系统结构图可知系统开环传递函数为:，即系统为II型系统；系统的静态位置误差系数，静态速度误差系数和静态加速度误差系数分别为：；由此可知系统在单位阶跃信号、单位斜坡信号和单位加速度信号作用下的稳态误差分别为：，根据线性系统的叠加原理可得，在信号r(t)作用下系统的稳态误差为：4.10已知单位负反馈系统开环传递函数，求作用下的稳态误差。解：由于系统开环传递函数为:，即系统为I型系统；系统的静态位置误差系数，静态速度误差系数和静态加速度误差系数分别为：；由此可知系统在单位阶跃信号、5t信号和单位加速度信号作用下的稳态误差分别为：，根据线性系统的叠加原理可得，在信号r(t)作用下系统的稳态误差为：4.11一直单位负反馈系统开环传递函数为，求作用下的稳态误差。解：由于系统开环传递函数为:，即系统为II型系统；系统的静态位置误差系数，静态速度误差系数和静态加速度误差系数分别为：；由此可知系统在单位阶跃信号、单位斜坡信号和作用下的稳态误差分别为：，根据线性系统的叠加原理可得，在信号r(t)作用下系统的稳态误差为：5.3基于MATLAB的根轨迹分析法5.3.1利用MATLAB绘制根轨迹利用伊凡斯给出的绘制根轨迹的基本规则，可以粗略地画出当系统某一参数变化时的根轨迹，但需要花费较多时间，且结果并不精确。而使用MATLAB的相关指令，绘制较为精确的根轨迹就非常方便。绘制根轨迹的常用指令为：rlocus(num,den);或rlocus(num,den,K);绘制例5-2给定单位负反馈控制系统的根轨迹图。clccleark=1;z=;p=0-1-2;num,den=zp2tf(z,p,k);%将传递函数由零极点形式转换为多项式形式rlocus(num,den);V=-32-33;axis(V);title(root-locusplotofG(s)=K/s(s+1)(s+2);xlabel(Re);ylabel(Im);返回的给定控制系统根轨迹图为：对例5-3给定的自动控制系统绘制其根轨迹图。clcclear%k(s+1)/s(s+2)(s+3)k=1;z=-1;p=0,-2,-3;n,d=zp2tf(z,p,k);rlocus(n,d)执行本程序后可得给定系统根轨迹图为：对例5-4所给定的自动控制系统绘制根轨迹图clcclear%k/s(s+3)(s2+2s+2)g=tf(1,conv(1,3,1,2,2)0);rlocus(g)执行本程序后返回给定自动控制系统根轨迹图为：对例5-5所示系统绘制根轨迹。clcclearG1=tf(1,18);G2=tf(11,150);H=tf(1,12);rlocus(G1*G2*H);V=-102-55;axis(V);gridon;xlabel(Re);ylabel(Im);执行本程序后得到系统根轨迹为5.3.2基于根轨迹的系统性能分析当做出控制系统根轨迹图之后，就可以根据根轨迹对系统进行定性的分析和定量的计算。因为系统的暂态性能和稳态性能与系统闭环极点位置密切相关，实际工程中对系统性能的要求往往可以转化为对闭环极点位置的要求。1、分析1：在对系统的分析中，一般需要确定根轨迹上某一点的根轨迹增益及其对应的闭环极点。2、分析2：在对系统的分析中，有时需要确定具有指定阻尼比的主导闭环极点及相对应的开环增益值。如下例所示：clcclearn=1;d=conv(1,1,1,2)0;kos=0.5,0.707;w=0.5,1;sgrid(kos,w);holdonrlocus(n,d);k,p=rlocfind(n,d);holdoff3、分析3-