基于MATLAB的二阶系统分析

基于MATLAB的二阶系统分析凡是以二阶微分方程描述运动方程的控制系统，称为二阶系统。在控制工程中，不仅二阶系统的典型应用极为普遍，而且不少高阶系统的特性在一定条件下可用二阶系统的特性来表征。因此，着重研究二阶系统的分析和计算方法，具有较大的实际意义。1.典型二阶系统的暂态分析典型二阶系统的暂态分析是从时域方面对二阶系统进行分析。时域分析具有非常直观的分析效果，例如：给系统输入端加上阶跃信号观察系统的输出状况即二级系统的单位阶跃响应，能够很直观、很全面的对所研究的二阶系统作出全面了解。但在计算机尚未普及之前，对二阶系统单位阶跃响应曲线的绘制全依赖于人们的手工描绘，所以，对简单的、低阶系统尚能用时域法进行分析，但对于高阶系统的单位阶跃响应曲线就很难依赖手工绘制。因此，这位系统的暂态分析提出了很大挑战。然而，随着计算机技术的发展，用计算机设计的控制系统的计算机辅助设计软件层出不穷，这为控制系统的暂态分析提供了方便。因此，基于MATLAB的二阶系统分析，就是利用现在在控制系统分析、系统仿真等领域中应用非常广泛的MATLAB语言作为分析工具。1.1典型二阶系统的数学模型分析在研究典型的二阶系统时常用的数学模型有：T2血+2訂蚁+)=r(t)(l)dt2dt00①(s)二C(s)二1=馀R(s)T2s2+2gTs+1s2+2口s+32nn其中，g为系统的阻尼比，3为无阻尼自然震荡频率。n公式(1)是对二阶系统的微分方程描述，公式(2)是对二阶系2)统的传递函数描述。1.2典型二阶系统的单位阶跃响应典型二阶系统的特征方程为：D(s)=s2+2g3s+32=0nn特征根为：s=-gw±Jg2-131,2nn3)4)由公式(4)可以看出，特征根的分布主要取决于系统的阻尼比g。而系统在零初始条件下，典型二阶系统的单位阶跃响应：1321C(s)"(s)_二-ss2+2g3s+32snn单位阶跃响应的特征主要取决于特征根的分布，当3二1rads时,n5)取不同的阻尼比g时的到得阶跃响应曲线如下所示:21.81.61.41.210.80.60.40.2图1不同阻尼比下的阶跃响应因此，根据系统的阻尼比g的不同，把二阶系统分为几种不同的状态如下：1，临界阻尼状态分析当1时，特征根为重负实根S，系统的单位阶跃响应曲n线如下图所示：1987000654000302100edutilpmA1987000654000302100edutilpmA202530图2二阶系统临界阻尼状态由临界阻尼状态下系统的单位阶跃响应曲线可看出，当t二0时，响应过程的变化率为零；当t>0时，响应过程的变化率为正，响应过程单调上升；当t十时，响应过程的变化率趋于零，响应过程趋于常数1。g>1，过阻尼状态分析当g>1时，特征根为两个不相等的负实根，系统的单位阶跃响应曲线如下所示：

图3二阶系统过阻尼状态由过阻尼状态下系统的单位阶跃响应曲线可看出，系统的单位阶跃响应与临界阻尼状态相似，但响应速度比临界阻尼状态缓慢。o<g<1，欠阻尼状态分析当0<g<1时，特征根为一对共轭复数s=-»土加（耳，系nn统的单位阶跃响应曲线如下所示：99)图4二阶系统欠阻尼状态由欠阻尼状态下系统的单位阶跃响应曲线可以看出，欠阻尼状态下的响应曲线有两个特征：衰减振荡和振荡频率不变。由于二阶系统的欠阻尼状态在工程实践中具有广泛应用，因此，对二阶系统的欠阻尼状态的暂态性能指标做如下论述：对公式（5）作拉普拉斯逆变换的单位阶跃响应为：TOC\o"1-5"\h\zc（t）=1-e-sin（①t+申）（6）心2d1其中3=3J1-g2（7）dn3称为有阻尼振荡频率。d上升时间tr对于欠阻尼系统来说当系统第一次到达稳态值时所用的时间即为上升时间。所以，有c（t）=1结合公式（6）得：兀7t=r3n由公式（8）值知，在g—定时，上升时间t与有阻尼振荡频率3rd成反比。超调时间tp超调时间t为系统从初始时刻到达到最大值时的时间，所以只需p对公式（6）进行求导使导数为零得t为：冗t=-p3d由t的计算公式知t只与3有关且成反比。ppd超调量5%由计算超调时间时可知5%为:p5%=ei弋2*100%（10）p由5%的计算公式知由计算超调时间时可知5%为:p5%=ei弋2*100%（10）p由5%的计算公式知5%只与g有关。pp④调节时间ts调节时间t与调节精度有关因此有：s4T11)当A二5时t二丄二3Ts3n12)对欠阻尼状态下各性能指标在单位阶跃响应曲线中的位置如下图所示：1.41.2StepResponseSystem:G[11Pteakamplitude:1.16Overshoot(%):16.3.Attime(sec):3.62System:GSystem:G-——SettlingTime(sec):8.08FinalVabe:1edutilpmA8600System:GRiseTime:1.6551015202530Time(sec)图5欠阻尼状态下暂态性能指标1.2.4o，无阻尼状态分析当g=0时，特征根为一对纯虚数s=±jo，系统的单位阶跃响应n曲线如下所示：0051015202530Time(sec)0051015202530Time(sec)StepResponse图5二阶系统无阻尼状态由无阻尼状态下系统的单位阶跃响应曲线可看出，单位阶跃响应为等幅振荡，此时系统处于临界稳定状态。2.典型二阶系统的根轨迹分析控制系统的闭环极点在复平面上随系统参数变化的轨迹称为控制系统的根轨迹。以根轨迹放大系数k为参变量的根轨迹称为常义根轨迹。现以常义根轨迹来分析欠阻尼系统如下：通过根轨迹放大系数k的变化，通过MATLAB作图，可以很容易的看出欠阻尼系统的闭环极点在复平面上随系统参数变化，通过观察欠阻尼系统的闭环极点在复平面的根轨迹变化可以很容易判断出系统的稳定性。若控制系统的闭环传递函数为：0(s)=1(13)S2+S+1则该系统的开环传递函数为：G(s)=1(14)s2+s得该系统的根轨迹如下图所示：RealAxis图6根轨迹曲线由上图的根轨迹曲线可以看出随着开环放大系数k的变化，系统的根轨迹始终都在复平面虚轴的左侧，系统稳定，这与暂态分系统中由单位阶跃响应曲线中的到的结论一致。3.典型二阶系统的频域分析前面的时域分析法和根轨迹分析法在低阶系统中尚能灵活应用，担当系统是高阶系统时，系统的微分方程求解和根轨迹绘制都变得非常困难。因此，在经典控制理论中还用频域分析法对控制系统的稳定性和稳定裕度进行分析。所以，以下还是用以上的欠阻尼系统应用频域分析法对其分析如下：首先，利用MATLAB绘制欠阻尼系统的奈奎斯特曲线，再由奈奎斯特判据判断其稳定性如下所示：RealAxis图7奈奎斯特曲线由欠阻尼系统的奈奎斯特曲线可以看出z=0即闭环极点位于右半平面极点的个数为灵。所以，系统稳定，与以上两种方法判断结果一致。利用奈奎斯特曲线和奈奎斯特判据仅能判断出系统是否稳定，要想进一步知道稳定裕度还要用伯德图进行判定。所以，以下是用伯德图法对欠阻尼系统进行分析，从伯德图中读出系统的幅值裕度相位裕度。二阶欠阻尼系统的伯德图如下图所示：)Bd(edutingaM)ged(esah200-20-40-60-800-45-90)Bd(edutingaM)ged(esah200-20-40-60-800-45-90-135BodeDiagram-18010-2■1r10-1100101Frequency(rad/sec)102图8二