控制系统计算机仿真实验指导书

《控制系统计算机仿真实验指导书》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、前言连续系统数字仿真方法是指用数字计算机对连续系统进行仿真的方法。采用这种方法时首先将连续系统的数学模型转变为适合在数字计算机上进行试验的仿真模型，实现这种转变的计算方法主要有微分方程数值解法和离散相似法。MATLAB产品家族是美国MathWorks公司开发的用于概念设计、算法开发、建模仿真、实时实现的理想的集成环境。是矩阵实验室（MatrixLaboratory）的简称，是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和SIMULINK两

2、大部分。MATLAB由于其完整的专业体系和先进的设计开发思路，使得MATLAB在多种领域都有广阔的应用空间，特别是在科学计算、建模仿真以及系统工程的设计开发上已经成为行业内的首选设计工具，它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。MATLAB软件工具在自动化专业、测控

3、技术与仪器和电气工程及其自动化等专业的本科生学习中，经常用来计算、仿真和设计，尤其是MATLAB软件的仿真功能，能使我们对所学知识有更加深入的理解和分析。40目录一、系统分析：31.绘制碾磨控制系统开环根轨迹图、BODE图和奈奎斯特图，并判断稳定性；42.当控制器为，试设计一个能满足要求的控制器（要求用根轨迹法和频率响应法进行设计）；72．1进行根轨迹校正：82.2频率校正：123.将采样周期取为，试确定与对应的数字控制器（要求用多种方法进行离散化，并进行性能比较）；174、5、6：连续，离散单位阶跃输入响应比较256.比较并讨论4和5

4、的仿真结果；277.讨论采样周期的不同选择对系统控制性能的影响；288.如控制器改为PID控制器，请确定满足性能指标的PID控制器参数。338.1用最优PID控制法设计：339.如希望尽可能的缩短系统的调节时间，请设计相对应的最小拍控制器，并画出校正后系统的阶跃响应曲线。379.1采用无纹波最少拍系统设计3710.课程设计总结：3911.参考资料：…………………………………………………………………………………………4040采样控制系统分析与仿真专业课程设计11计算机应用技术石王魁一、系统分析：某工业碾磨系统的开环传递函数为要求用数字控制

5、器D(z)来改善系统的性能，使得相角裕度大于，调节时间小于1s(2%准则)1.绘制碾磨控制系统开环根轨迹图、Bode图和奈奎斯特图，并判断稳定性；2.当控制器为，试设计一个能满足要求的控制器（要求用根轨迹法和频率响应法进行设计）；3.将采样周期取为，试确定与对应的数字控制器（要求用多种方法进行离散化，并进行性能比较）；4.仿真计算连续闭环系统对单位阶跃输入的响应；5.仿真计算数据采样系统对单位阶跃输入的响应；6.比较并讨论4和5的仿真结果；7.讨论采样周期的不同选择对系统控制性能的影响；8.如控制器改为PID控制器，请确定满足性能指标的

6、PID控制器参数。9.如希望尽可能的缩短系统的调节时间，请设计相对应的最小拍控制器，并画出校正后系统的阶跃响应曲线。40采样控制系统分析与仿真专业课程设计11计算机应用技术石王魁1.绘制碾磨控制系统开环根轨迹图、Bode图和奈奎斯特图，并判断稳定性：G=zpk([],[0-5],10);sisotool(G);margin(G);根轨迹图40采样控制系统分析与仿真专业课程设计11计算机应用技术石王魁Bode图：截止频率为1.88rad/s,相角裕度为6940采样控制系统分析与仿真专业课程设计11计算机应用技术石王魁N=0;R=0;Z=P

7、-R=0;该系统稳定。40采样控制系统分析与仿真专业课程设计11计算机应用技术石王魁2.当控制器为，试设计一个能满足要求的控制器（要求用根轨迹法和频率响应法进行设计）：调节前Gs=tf(10,[150]);Close_S=feedback(Gs,1);Step(Close_S,'b');holdon设计前调节时间为1.18s40采样控制系统分析与仿真专业课程设计11计算机应用技术石王魁设计前截止频率为1.88rad/s,相角裕度为69°（第一问中）2．1进行根轨迹校正：要使得根轨迹向左转，要加入零点。考虑到校正装置的物理可实现性，加入超

8、前校正装置。40采样控制系统分析与仿真专业课程设计11计算机应用技术石王魁检验性能：>>Ds=tf(10*[16.512],[111.499]);Gs=tf(10,[150]);Close_S=feedba