三阶系统综合分析与设计.docx

目录摘要O1设计内容11*1I卜*+ii*iii*+*ii*+*i*ii+*+ii*ii*+.*+ii*+*Ei+Bi*1+*+，*+*!ii.+*+ii.*+ii.*ii.+ii.*+*ii.+ii.*:+.*+i12绘制三阶系统的根轨迹22J常规方法绘制根轨迹22.2用MATLAB绘制根轨迹33不同条件下K的取值33.1当-8为闭环系统的一个极点时，K的取值33.2主导极点阻尼比为0.7时的k值34求系统的稳态误差44,1位置误差系数44*2J度误差系数543力口j，度i吴.系数54.4输入信号为）=i十2.5f十产时的稳态误差55绘制单位阶跃响应曲线.66频域特性分析66.1绘制Bode图和Nyquist曲线66.2相角裕度和幅值裕度87加入非线性环节判断稳定性.*.*.87.1求死区特性环节的描述函数87.2根据负倒描述函数和Nyquist图判断系统的稳定性9摘要三阶系统是以三级微分方程为运动方程的控制系统。在控制工程中,三阶系统非常普遍，其动态性能指标的确定是比较复杂。在工程上常采用闭环主导极点的概念对三阶系统进行近似分析，或直接用MATLAB软件进行高阶系统分析。在课程设计中，要掌握用MATLAB绘制闭环系统根轨迹和系统响应曲线，用系统的闭环主导极点来估算三系统的动态性能，以及在比较点及开环传递函数之间加一个非线性环节判断其稳定性。1设计内容1.1设计题目三阶系统的综合分析和设计初始条件：某单位反馈系统结构图如图IT所示：图IT图1-21.2设计任务要求完成的主要任务：(包括课程设计工作量及其技术要求，以及说明书撰写等具体要求)1、试绘制随根轨迹2、当-8为闭环系统的一个极点时，K=?3、求取主导极点阻尼比为0.7时的K值(以下K取这个值)4、分别求取位置误差系数、速度误差系数、加速度误差系数及输入信号为r(f)=I(Z)+2.5?+/单位阶跃信号、斜坡信号及单位加速度信号时的稳态误差5、用MatIab绘制单位阶跃相应曲线6、绘制Bode图和Nyquist曲线，求取幅值裕度和相角裕度7、如在比较点及开环传递函数之间加1个死区非线性环节，如图1-2所示，其中%=LA=2,试求取非线性环节的描述函数，并根据负倒描述函数和Nyquist图判断系统的稳定性8、认真撰写课程设计报告。2绘制三阶系统的根轨迹图2系统结构图由图1可得，三阶系统的开环传递函数为：G(三)=2.1常规方法绘制根轨迹根据绘制根轨迹的规则，可知该系统的根轨迹绘制步骤如下：(1)根轨迹的起点和终点。根轨迹起于开环极点，终于开环零点。开环极点分别为0、-3、-6,无开环零点。(2)根轨迹的分支数。n=3,m=0,所以分支数为3。且它们是连续的并且对称于实轴。(3)根轨迹的渐进线。本系统根轨迹的渐近线有三条，据其及实轴的夹角公式：把n=3,m=0代入求得：=0-3-6=渐近线及实轴的交点为：4n-m3-0根轨迹在实轴上的分布。实轴上的某一个区域，若其右边开环零、极点的个数之和为奇数，则该区域必是根轨迹。因此实轴上-8,-6-3,0必为根轨迹。(5)确定根轨迹的分离点。该系统中没有有限零点，由法则五得：于是分离点方程为：因此可以求得分离点d=-L3,d=T.7(不合题意，舍去)(6)根轨迹及虚轴的交点。闭环特征方程式为对上式应用劳斯判据，有：?118?9kZk令劳斯表中3行的首项为零，得k=162,根据S2行的系数，得辅助方程91+k=0代k=162并令s=jw,解得交点坐标2.2用MATLAB绘制根轨迹MATLAB为绘制根轨迹编程如下：num=1;den=l9180;syms=tf(num,den);rlocus(syms)MATLAB产生的根轨迹如图2所示：图3闭环系统的根轨迹3不同条件下K的取值3.1当-8为闭环系统的一个极点时，K的取值由图1的系统的闭环传递函数为闭环特征方程式为把s=-8代入上式中解得K=803.2主导极点阻尼比为0.7时的k值在控制工程实践中，通常要求控制系统即具有较快的响应速度，又具有一定的阻尼程度，此外还要求减少死区间、间隙和库仑摩擦等非线性因素对系统性能的影响，因此高阶系统的增益常常调整到是系统具有一对闭环共学Ii主导极点。这时，可以用二阶系统的动态性能指标来估算高阶系统的性能。由于主导极点阻尼比?=0.7<1,属于欠阻尼系统。由公式,=cos/得:设系统的自然频率为%,阻尼比4=0.7,由上述用闭环主导极点分析高阶系统的方法可知，距虚轴最近的一对闭环共独主导极点为：代入数据：邑,2=-0。7%土%/Jl-(0.7r闭环特征方程式为代入S的方程化简得：0.686MT2.6%+k+j(0.686*-8.82一+12.6%)=0分别令实部和虚部为零得到两个方程：0.686-12.6w,+k=00.686*-8.82炉+12.6叫=0解得叱=L6k=17.6%=11左=-808.6(舍)所以取人=17.6。所以，开环传递函数为4求系统的稳态误差系统的误差e(t)一般定义为输出量的希望值及实际值之差，一般情况下采用从系统输入端定义的误差e(t)来进行计算分析。控制系统的稳定系统误差信号的稳态分量称为系统的稳态误差，以e.表示4.1位置误差系数对于单位阶跃输入，R(三)=I/s,求得系统的稳态误差为令KP=HmG(三)"(s)，称KP为稳态位置误差系数。sO稳态误差可表示为对于0型系统，V=O,则勺=K;当vl时，KP=8。由图1系统得：由上式得V=1,则KP=8=-04.2速度误差系数对于单位斜坡输入，此时系统的稳态误差为Ku=IimSG(三)"称为稳态速度误差系数。s0于是稳态误差可表示为由图1系统得；£=linsG(G/G)=IimI譬=0.978o.o18=l0224.3加速度误差系数对于单位抛物线输入，此时系统的稳态误差为令旦=IimfG发称K"为稳态加速度误差系数。于是稳态误差可表示为于是稳态误差可表示为则对于图1系统得：KiI=Iiml=O$>04.4输入信号为R)=Ia)+25f+"时的稳态误差系统的开环传递函数为(K=17.617)当八=l(r)时，KP=8,=0当弓=2.5/时，Kv=O.978,e“2=至=2.556当时，K"=0，=8则当输入信号为=i+2.5/+产时5绘制单位阶跃响应曲线MATLAB绘制单位阶跃响应曲线编程如下：num=17,6;den=l91817.6;syms=tf(num,den);step(tf(num,den)即系统单位阶跃相应曲线为图3所示：图4系统单位阶跃相应曲线6频域特性分析6.1绘制Bode图和Nyquist曲线6.1.1绘制Bode图Bode图又称对数频率特性曲线图，由对数幅频曲线和对数相频曲线组成，是工程中广泛使用的一组曲线。手工绘制伯德图的步骤如下：1、将传递函数写成伯德标准型，确定开环传递系数和各转折频率。系统的开环频率特性为由伯德标准型容易看书，开环传递系数为K=0.978,转折频率为0=3,<y,=6°2、确定低频段：由传递函数可知该系统为1型系统，即有微分环节，所以绘制低频段，可过=1,L=20logK=20log0.978=19作一条斜率为-20dB/dec的斜线。3、绘制开环对幅频特性的渐近线：将低频段延伸到第一个转折点频率%=3处。因为第一个转折频率是惯性环节的转折频率，所以，开环对数频率特性的渐近线下降20dB/dec,再延伸到第二个转折频率吗=6处，因为也是惯性环节，所以再下降20dB/dec。4、绘制相频特性：绘制各个环节的对数相频特性曲线，然后逐点叠加。一般在一些特征点上进行叠加，如各个转折频率处。5、修正对数幅频特性。Matlab编写程序时可以利用函数Bode和Inargin,但是margin函数可以直接读出相角裕度和幅值裕度，所以下面用margin函数绘制BODE图MATLAB绘制BODE图编程如下：num=17.6;den=l9180;syms=tf(num,den);margin(num,den)绘制出的BODE图如图2-3所示：图5开环传递函数的Bode图从图中可以看出：幅值裕度h=19.3dB,穿越频率吗=4.24radsec0相角裕度r=64deg,截止频率吗=0.924radsec°从幅值裕度和相角裕度可以看出这是一个稳定的系统。6.1.2绘制LyqUiSt曲线Nyquist曲线图又称开环幅相频率特性曲线图，是利用Nyquist判别系统稳定性的依据。MATLAB绘制NyqUiSt曲线编程如下：num=17.6;den=l9180;syms=tf(num,den);nyquist(syms)绘制出的NyqUiSt曲线如图图6开环传递函数NyqUiSt曲线6.2相角裕度和幅值裕度相角裕度7的含义是，对于闭环稳定系统的，如果系统开环相频特性再后7度，则系统将处于临界稳定状态。设也为系统的截止频率，显由上式可得：=1解得：wc=0.92rads相角裕度：=180°+GOc)(j)=64.幅值裕度人的含义是，对于闭环稳定系统，如果系统开环幅频特性再增大人倍，则系统将处于临界稳定状态。设也为系统的穿越频率，则在吗处的相角：!>(>vx)=Gjwx)H(jwx=(2k+l);左=0,土1,由上式得G(jw)H(jwc)=180叭=4.24radS幅值裕度为()=-20IgIG(JW)"(/%)=19.3d8上述为笔算结果，通过MATLAB绘制的Bode图中得到的相角裕度7为64。，幅值裕度为19.348,两种方法得到结果基本吻合。7加入非线性环节判断稳定性图7加非线性环节系统结构图在比较点及开环传递函数之间加一个死区特性非线性环节，如图6所示。7.1求死区特性环节的描述函数由正弦输入信号皿)、死区特性NX)可得死区特性环节输出的数学表达式为0f炉I由'K(Asinwr-e0)xwt!2于y(f)为奇函数，所以AO=0,A=O,而y(f)又为半周期内对称，故死区特性的描述函数为代入得：=U=2即得到描述函数为7.2根据负倒描述函数和NyqUiSt图判断系统的稳定性据图2得，对于线性环节，K=7617>7j=：5=g解得穿越频率为：非线性环节为死区特性，负倒描述函数为即-UN（八）的曲线如图7中横轴射线所示。即曲线