基于Matlab/Simulink的二阶控制系统仿真研究

　　为了研究二阶控制系统的性能，讨论了二阶控制系统参数ζ和ωn与单位阶跃响应的关系，并介绍了基于Matlab/Simulink软件仿真环境，在单位阶跃信号作用下，利用仿真实例很好地实现了对二阶控制系统进行仿真研究。该仿真具有直接观察和分析二阶控制系统输出性能变化，同时充分体现利用Matlab/Simulink进行仿真研究直观和方便的特点。

　　0 引言

　　大部分实际控制工程都是复杂的高阶系统，通常可将其分解为若干个一阶惯性环节和二阶振荡环节的叠加。在研究复杂控制系统的过程中，往往通过主导极点和偶极子相消的方法，将高阶系统简化为二阶等低阶系统的组合，再根据对二阶控制系统输出和性能的分析，对控制系统偏差进行校正，从而获得满足设计指标和使用要求的高阶系统性能。因此，掌握二阶控制系统动态特性理论，构建二阶控制系统动态特性的测试系统，可大大简化系统分析，对控制系统的分析和设计具有重要作用。

　　Simulink 是一种以Matlab 为基础的实现交互式动态系统建模、仿真与分析的软件包，可以针对控制系统等进行系统建模、仿真、分析等工作，被广泛应用于线性控制系统、非线性控制系统、数字控制及数字信号处理的建模和仿真中。本文介绍基于Matlab/Simulink 仿真环境下，通过单位阶跃信号作用，利用仿真实例对二阶控制系统进行仿真研究，直接观察和分析二阶控制系统输出性能的变化，充分体现了Matlab/Simulink仿真直观和方便的特点。
　　1 二阶控制系统模型

　　能够用二阶微分方程描述的系统称为二阶控制系统。在控制工程实践中，二阶控制系统十分常见，例如，电枢控制的直流电动机，RLC网络和弹簧-质量-阻尼器组成的机械位移系统等。此外，许多高阶系统在一定条件下，常常近似地作为二阶控制系统来研究。因此，详细讨论和分析二阶控制系统的特性，具有极为重要的实际意义。典型二阶控制系统数学模型结构如图1所示。

　　其闭环传递函数为：

　　式中：ζ为系统阻尼比;ωn为无阻尼自然振荡角频率，单位为rad/s.二阶控制系统的动态特性可由以上两个参数描述。二阶控制系统在单位阶跃信号作用下，其输出响应可分为以下几种情况：

　　（1）当ζ=0时，二阶控制系统为零阻尼状态。系统有一对共轭虚根，系统单位阶跃响应为无阻尼等幅振荡曲线。

　　（2）当0《ζ《1 时，二阶控制系统为欠阻尼状态。系统的极点为共轭复数，位于S 左半平面。系统单位阶跃响应由稳态响应和瞬态响应两部分组成，稳态响应为1，瞬态响应为振荡衰减过程，振荡角频率由阻尼比ζ和无阻尼自然振荡角频率ωn决定，并且随着ζ的减小，其振荡幅度加大。

　　（3）当ζ=1时，二阶控制系统为临界阻尼状态。系统具有两个相同的实数极点，位于S 左半平面。系统单位阶跃响应为无超调，无振荡单调上升的曲线，不存在稳态误差。

　　（4）当ζ》1时，二阶控制系统为过阻尼状态。系统具有两个不等的实数极点，位于S左半平面，系统单位阶跃响应与临界阻尼情况相似，为无超调，无振荡单调上升的曲线，但它的过渡过程时间较之临界阻尼更长。

　　2 二阶控制系统仿真设计与研究

　　2.1 二阶控制系统仿真结构设计

　　启动Matlab 7.04，进入Simulink 仿真界面，根据二阶控制系统数学模型结构，设计二阶控制系统的Simu-link仿真结构如图2所示。双击各函数模块，在出现的各参数对话框内设置相应的参数，仿真时输入单位阶跃信号，起始时间为0，分别改变ωn 和ζ的值，点击simula-tion菜单下的start命令进行仿真，双击示波器模块观察仿真结果，得到系统的阶跃响应曲线，再进行分析ωn和ζ对系统动态性能的影响。

　　2.2 二阶控制系统单位阶跃响应与参数ζ的关系

　　设定ωn=10 rad/s 不变，改变参数ζ分别为0，0.25，1和2的二阶控制系统Simulink仿真结构如图3所示，输入单位阶跃信号，其仿真响应曲线如图4所示。从图中实验数据分析可以看出，响应曲线由上至下，依次为无阻尼等幅振荡曲线，欠阻尼振荡衰减曲线，临界阻尼和过阻尼无超调单调上升曲线。当0《ζ《1 时，二阶控制系统欠阻尼状态ζ变化的阶跃响应曲线如图5所示。随着ζ 的增大，系统单位阶跃响应的超调量减少，但上升时间加长，曲线峰值较大，因此，综合考虑超调量和上升时间两个因素，应选择ζ 接近最佳阻尼比0.707。

　　2.3 二阶控制系统单位阶跃响应与参数ωn的关系

设定ζ=0.1 不变，改变参数ωn 分别为10 rad/s 和20 r