自动控制原理课程设计-高阶系统的时域分析.doc

武汉理工大学自动控制原理课程设计说明书目 录1系统稳定性分析12高阶系统的时域响应22.1系统单位阶跃响应曲线22.2系统单位斜坡响应曲线42.3系统单位加速度响应曲线52.4动态性能指标计算62.4.1动态性能指标计算62.4.2动态性能指标计算72.5 稳态性能指标计算93根轨迹图绘制103.1根轨迹数据计算103.2用matlab绘制根轨迹图114心得体会12参考文献13本科生课程设计成绩评定表高阶系统的时域分析1 系统稳定性分析给定参数系统稳定性分析：对于开环传递函数在给定条件k=15,a=2,b=4时用劳斯判据判断系统的稳定性，经过化简可得系统的特征方程为:d(s)=s4+7s3+19s2+33s+60=0其劳斯表为 s4 1 19 60 s3 7 33 0 s2 14.3 60 s1 3.6 0 s0 60从表中可以看出，第一列系数符号全部为正，故系统是稳定的。2 高阶系统的时域响应k=15，a=2，b=4时，系统的开环传递函数为:gp= 15（s+4）ss2+5s+9s+2=15s+60s4+7s3+19s2+18s系统为型系统，可以跟踪单位阶跃信号、单位斜坡信号，不能跟踪单位加速度信号。系统响应为c(s)= 15s+60s4+7s3+19s2+33s+60r(s)2.1系统单位阶跃响应曲线当输入为单位阶跃函数信号时，r(s)=1s，系统响应为c(s)= 15s+60s4+7s3+19s2+33s+601s运用matlab程序作图如图2-1，程序为：num=15 60;den=1 7 19 33 60;g=tf(num,den);step(g);grid on;xlabel(t );ylabel(c(t) );title(单位阶跃响应)图2-1 系统阶跃响应由图中数据可得:上升时间为tr=1.02s峰值时间tp=1.73s调节时间ts=33.1s超调量%=78%稳态误差为ess=02.2系统单位斜坡响应曲线 当输入为单位斜坡函数信号时，r(s)=1s2,系统响应为c(s)= 15s+60s4+7s3+19s2+33s+601s2运用matlab程序作图如图2-2，程序为：num=15 60;den=1 7 19 33 60;g=tf(num,den);t=0:0.01:10;u=t;lsim(g,u,t);grid on; xlabel(t);ylabel(c(t);title(单位斜坡响应) 图2-2 单位斜坡响应2.3系统单位加速度响应曲线当输入为单位加速度函数信号时，r(s)=1s3,系统响应为c(s)= 15s+60s4+7s3+19s2+33s+601s3运用matlab程序作图如图2-3，程序为：num=15 60;den=1 7 19 33 60;g=tf(num,den);t=0:0.01:10;u=(0.5*t.2);lsim(g,u,t)grid on;xlabel(t);ylabel(c(t) );title(单位加速度响应)图2-3 单位加速度响应2.4动态性能指标计算2.4.1主导极点法控制系统的暂态性能指标通常是零初始条件下，通过系统的阶跃响应的特征定义的，系统的暂态性能指标实际上就是刻画阶跃响应曲线特征的一些量。本系统属于高阶系统，无法用处理二阶系统的方法得到系统性能指标的精确解析表达式。即使用数值计算的方法，也会遇到求解复杂的超越方程的困难。所以对于高阶系统的暂态性能分析一般采取主导极点法。系统的闭环传递函数为：(s)= m(s)d(s)=15(s+4)s+0.085-j2.12s+0.085+j2.12s+3.4+j1.29(s+3.4-j1.29)取主导极点为s1,2=-0.0852.12j闭环零点z1=-4则高阶系统的单位阶跃响应可以近似为：c(t)=m(0)d(0)+2m(s1)s1d(s1)e-tcosdt+m(s1)s1d(s1)tp1d-i=11s1-zi+i=34s1-si计算得： tp=12.12(-28+45.8+14)=1.74sp=i=34sii=34s1-sii=11s1-zii=11zie-tp计算得： p=0.7932.4.2 matlab编程动态性能指标是指稳定的系统在单位阶跃函数作用下，动态过程随时间t的变化状况的指标，体现系统动态过程特征。用解析法求解高阶系统的动态性能指标很困难，这里用matlab编程求解。调用单位阶跃响应函数step()，获得系统的单位阶跃响应，当采用y,t=step(sys)的调用格式时，将返回值y及相应的时间t，通过对y和t进行计算，可以得到高阶系统各项动态性能指标。利用matlab编程求取系统动态性能指标程序如下：sys=tf(15 60,1 7 19 33 60)c=dcgain(sys);y,t=step(sys); y,k=max(y);tp=t(k)mp=(y-c)/cn=1;while y(n)0.98*c)&(y(i)1.02*c)i=i-1; endts=t(i)结果如图所示：2.5稳态性能指标计算e(s)=11+15s+60s4+7s3+19s2+18sr(s)=s4+7s3+19s2+18ss4+7s3+19s2+33s+60r(s)当输入信号为单位阶跃响应时，r(s)=1sse(s)= s4+7s3+19s2+18ss4+7s3+19s2+33s+60ess=lims0se(s)=0当输入为单位斜坡函数信号时，r(s)=1s2se(s)= s3+7s2+19s+18s4+7s3+19s2+33s+60ess=lims0se(s)=1860=0.3当输入为单位加速度函数信号时，r(s)=1s3se(s)= s3+7s2+19s+18s5+7s4+19s3+33s2+60sess=lims0se(s)=3 根轨迹图绘制3.1根轨迹数据计算当a=1，b=4时，gf= k（s+4）ss2+5s+9s+1根轨迹相关参数计算如下：1） 系统开环零点为z=-4；2） 开环极点为s1=0,s2=-2.5+j1.66,s3=-2.5-j1.66,s4=-1；3） 系统有max m，n =4 根分支； 4）实轴上的根轨迹为（-1，0）和（-，-4）； 5）系统的m=1，n=4，故根轨迹的渐近线为3条，渐近线与实轴的交角分 3、-3、。渐近线与横轴交点为（-23，0）；6）根轨迹的分离点：由方程 ab-ba=0 3s4+28s3+86s2+112s+36=0可计算得约为（-5，0）和（-0.462，0）；7）根轨迹与虚轴的交点：闭环特征方程：s4+6s3+14s2+（9+k）s+4k=0其劳斯表为 s4 1 14 4k s3 6 9+k 0 s2 75-k6 4k s1 9+k-144k75-k 0 s0 4k用劳斯判据可算得w=1.68 k=7.863.2用matlab绘制根轨迹图运用matlab作图如图2-4所示，程序如下num=1 4;den=1 6 14 9 0;rlocus(num,den)图3-1 根轨迹图4 心得体会对高阶系统进行时域分析，运用经典解析方法，采用拉普拉斯反变换求解瞬态响应时域表达式比较复杂，要计算出各项动态性能指标也很困难。但对于许多高阶系统，利用主导极点法可以简化系统的分析和性能指标的估算。而利用matlab软件可以方便地对高阶系统时域响应进行准确分析。本次课程设计的重要部分为系统稳定性的判断以及在各种典型信号作用下系统的响应性能分析。在求解多元不等式组方面有较大的困难。通过查阅书本资料，觉得使用函数图像表征个参变量之间的关系比较可行，但实际做的过程中，发现三变量之间的函数图像并不好作，而且涉及到高次的方程求解，十分困难。在解高阶系统的斜坡响应与加速度响应时，做时域的变换也比较困难。在课程设计中，加深了对所学自动控制原理课程知识的理解，特别是系统稳定性分析，系统各项动态性能指标，稳态误差以及系统根轨迹等相关知识的理解。设计时借助matlab软件进行控制系统分析，进一步熟悉了matlab语言及其应用。书写课程设计说明书时使用word软件，使我掌握了许多word编辑和排版技巧。通过本次课设，我不但对所学的自动控制原理的相关知识有了更深入的认识和更牢固的掌握，而且学会了运用matlab来解决高阶系统的分析和设计的相关问题。总而言之，这次课程设计又让我重新巩固了自动控制原理的相关知识，特别是关于系统稳定性分析，系统各项动态性能指标，稳态误差以及根轨迹的绘制等。通过借助强大的matlab软件，让我熟悉了matlab的语言及应用，学会运用有效的工具解决人工计算比较复杂的情况，起到事半功倍的效果。因此，在以后的学习过程中，我们不仅要在课堂上加强理论基础的学习，在实际运用上还要根据实际情况运用各种工具来进行更精准的分析，不仅是死读书，而是在各方面提高自己的能力。参考文献1 胡寿松.自动控制原理（第四版）.科学出版社,20022