系统仿真综合实验指导书(2011.6)

系统仿真综合实验指导书电气与自动化工程学院自动化系2011年6月1前言电气与自动化工程学院为自动化专业本科生开设了控制系统仿真课程，为了使学生深入掌握MATLAB语言基本程序设计方法，运用MATLAB语言进行控制系统仿真和综合设计，同时开设了控制系统仿真综合实验，30学时。为了配合实验教学，我们编写了综合实验指导书，主要参考控制系统仿真课程的教材《自动控制系统计算机仿真》、《控制系统数字仿真与CAD》、《反馈控制系统设计与分析——MATLAB语言应用》及《基于MATLAB/Simulink的系统仿真技术与应用》。2实验一MATLAB基本操作实验目的1．熟悉MATLAB实验环境，练习MATLAB命令、m文件、Simulink的基本操作。2．利用MATLAB编写程序进行矩阵运算、图形绘制、数据处理等。3．利用Simulink建立系统的数学模型并仿真求解。实验原理MATLAB环境是一种为数值计算、数据分析和图形显示服务的交互式的环境。MATLAB有3种窗口，即：命令窗口（TheCommandWindow）、m-文件编辑窗口（TheEditWindow）和图形窗口（TheFigureWindow），而Simulink另外又有Simulink模型编辑窗口。1．命令窗口（TheCommandWindow）当MATLAB启动后，出现的最大的窗口就是命令窗口。用户可以在提示符“”后面输入交互的命令，这些命令就立即被执行。在MATLAB中，一连串命令可以放置在一个文件中，不必把它们直接在命令窗口内输入。在命令窗口中输入该文件名，这一连串命令就被执行了。因为这样的文件都是以“.m”为后缀，所以称为m-文件。2．m-文件编辑窗口（TheEditWindow）我们可以用m-文件编辑窗口来产生新的m-文件，或者编辑已经存在的m-文件。在MATLAB主界面上选择菜单“File/New/M-file”就打开了一个新的m-文件编辑窗口；选择菜单“File/Open”就可以打开一个已经存在的m-文件，并且可以在这个窗口中编辑这个m-文件。3．图形窗口（TheFigureWindow）图形窗口用来显示MATLAB程序产生的图形。图形可以是2维的、3维的数据图形，也可以是照片等。MATLAB中矩阵运算、绘图、数据处理等内容参见教材《自动控制系统计算机仿真》的相关章节。Simulink是MATLAB的一个部件，它为MATLAB用户提供了一种有效的对反馈控制系统进行建模、仿真和分析的方式。有两种方式启动Simulink:31．在Commandwindow中，键入simulink，回车。2．单击工具栏上Simulink图标。启动Simulink后，即打开了Simulink库浏览器（Simulinklibrarybrowser）。在该浏览器的窗口中单击“Createanewmodel（创建新模型）”图标，这样就打开一个尚未命名的模型窗口。把Simulink库浏览器中的单元拖拽进入这个模型窗口，构造自己需要的模型。对各个单元部件的参数进行设定，可以双击该单元部件的图标，在弹出的对话框中设置参数。实验内容1用MATLAB可以识别的格式输入下面两个矩阵12332357135732391894A144367823355422675342189543iiBi再求出它们的乘积矩阵C，并将C矩阵的右下角2×3子矩阵赋给D矩阵。赋值完成后，调用相应的命令查看MATLAB工作空间的占用情况。2分别用for和while循环结构编写程序，求出6323626302122222iiK3选择合适的步距绘制出下面的图形（1）1sin(/)t，其中0110t[.,]（2）sin(tan)tan(sin)tt，其中(,)t4对下面给出的各个矩阵求取矩阵的行列式、秩、特征多项式、范数。75350083341009103150037193......A，5765710876810957910B412345678910111213141516C，33245518118575131D5求解下面的线性代数方程，并验证得出的解真正满足原方程。(a)72124915327221151132130X，(b)1321390721264915321172211521X6假设有一组实测数据x0.10.20.30.40.50.60.70.80.91y2.32012.64702.97073.28853.60083.90904.21474.51914.82325.1275用最小二乘法拟合，求出相应的二次函数。7考虑线性微分方程(4)(3)353345sin(4/3)ttyyyyyuueet输入信号(1)(2)(3)(0)1,(0)(0)1/2,0.2,yyyy方程初值(1)试用Simulink搭建起系统的仿真模型，并绘制出仿真结果曲线。(2)将给定的微分方程转换成状态方程，并建立S函数，再利用Simulink进行仿真。8建立下图所示非线性系统的Simulink模型，并观察在单位阶跃信号输入下系统的输出曲线和误差曲线。5实验二经典控制理论实验目的以MATLAB及Simulink为工具，对控制系统进行时域、频域及根轨迹分析。实验原理1、时域分析法是根据系统的微分方程（或传递函数），利用拉普拉斯变换直接解出动态方程，并依据过程曲线及表达式分析系统的性能。时域响应指标如图所示。延迟时间td，指响应曲线第一次达到其终值一半所需要的时间。上升时间tr，指响应曲线从终值10%上升到终值90%所需要的时间；对于有振荡的系统，也可定义为响应从零第一次上升到终值所需要的时间。上升时间是系统响应速度的一种度量。峰值时间tp，指响应超过终值达到第一个峰值所需要的时间。调节时间ts，指响应达到并保持在终值±5%（或±2%）内所需要的时间。超调量σ%，指响应的最大偏离量h(tp)与终值h(∞)之差的百分比，即：%100)()()(%hhtph稳态误差，描述系统稳态性能的一种性能指标。2、频域分析法通常从频率特性出发对系统进行研究。在工程分析和设计中，通常把频率特性画成一些曲线，从频率特性曲线出发进行研究。这些曲线包括幅频特性和相频特性曲线，幅相频率特性曲线，对数频率特性曲线以及对数幅相曲线等，其中以幅相频率特性曲线，对6数频率特性曲线应用最广。对于最小相位系统，幅频特性和相频特性之间存在着唯一的对于关系，故根据对数幅频特性，可以唯一地确定相应的相频特性和传递函数。根据系统的开环频率特性去判断闭环系统的性能，并能较方便地分析系统参量对系统性能的影响，从而指出改善系统性能的途径。3、根轨迹是求解闭环系统特征根的图解方法。由于控制系统的动态性能是由系统闭环零极点共同决定，控制系统的稳定性由闭环系统极点唯一确定，利用根轨迹确定闭环系统的零极点在s平面的位置，分析控制系统的动态性能。实验内容*1．控制系统数学模型的转换《自动控制系统计算机仿真》教材第4章中的所有例题2．已知二阶系统闭环传递函数为：10210)(2sssG（1）编写程序求解系统的阶跃响应；计算系统的闭环根、阻尼比、无阻尼振荡频率；（2）编程计算并在同一幅图上绘出11和22的阶跃响应曲线；（3）编程计算并在同一幅图上绘出nn211和nn22的阶跃响应曲线。（4）绘出以下系统的阶跃响应曲线，并与原系统响应曲线比较（将5条曲线画在同一幅图上）102102)(21ssssG；102105.0)(222sssssG；1025.0)(223sssssG；42()210sGsss记录各系统的峰值和峰值时间，计算超调量。系统峰值峰值时间超调量G(s)G1(s)G2(s)G3(s)G4(s)73．已知某控制系统的开环传递函数1512(),.()()KGsKsss试绘制系统的开环频率特性曲线（奈魁斯特曲线和伯德图），判断闭环系统是否稳定？并求出系统的幅值裕量与相位裕量。4.已知开环传递函数)11.0()1()(2sssksG令k＝1作伯特图，应用频域稳定判据确定闭环系统的稳定性，并确定使系统获得最大相位裕度的增益k值。5．单位负反馈系统的开环传递函数为051051021011(.)()(.)(.)(.)KsGssss绘制根轨迹曲线，并求出使闭环系统稳定的K值范围。使用rltool（sys）函数观察不同K值的阶跃响应情况，分析K值变化对系统响应有何影响。6.非最小相位系统的开环传递函数如下32122432641060110600510111782130100ssssGsGssssssss()(),()(.)(.)绘制频率特性曲线，并解释为什么这样的系统被称为“非最小相位”系统，试从其频率特性加以解释。在Simulink环境中建立对应闭环系统的仿真模型，观察它们的阶跃响应。7．单位负反馈系统的开环传递函数为：0.50()()sGsGse，其中021()22Gsss利用Simulink工具给出Simulink仿真框图及阶跃响应曲线。8实验三PID控制器的设计实验目的研究PID控制器对系统的影响。实验原理1．模拟PID控制器控制系统结构如图所示模拟PID控制器的表达式为：PID0d()()()()ddtetxtKetKeKt式中，PK、IK和DK分别为比例系数、积分系数和微分系数。对上式进行拉普拉斯变换，整理后得到连续PID控制器的传递函数ICPDPDI1()1KGsKKsKTssTs显然，PK、IK和DK三个参数一旦确定（注意：IPI/TKK，DDP/TKK），PID控制器的性能也就确定下来。为了避免纯微分运算，通常采用近似的PID控制器，其传递函数为DCPID1()10.11TsGsKTsTs2.离散PID控制器如果采样周期为T，在第k个采样周期()et的导数可近似表示为；d()()[(1)]detekTekTtT在k个采样周期内对()et的积分可近似表示为；00()()kkTmetdtTemT9因此，离散PID控制器的表达式为；PID0()[(1)]()()()kmekTekTxkTKekTKTemTKT离散PID控制器的表达式可简化为；PID0()(1)()()()kmekekxkKekKTemKT离散PID控制器的脉冲传递函数为；1ICPD-1()(1)1-KGzKKzz式中，PK、IK和DK分别为比例系数、积分系数和微分系数。实验内容1．已知三阶对象模型3()1/(1)Gss，利用MATLAB编写程序，研究闭环系统在不同控制情况下的阶跃响应，并分析结果。(1)0,idTT时，在不同KP值下，闭环系统的阶跃响应；(2)10,pdKT时，在不同iT值下，闭环系统的阶跃响应；(3)1piKT时，在不同dT值下，闭环系统的阶跃响应；*2.已知被控对象为一电机模型，传递函数为2100067010()..Gsss，输入信号为0502().sin()rkt，采用PID控制方法设计控制器，利用MATLAB编程进行仿真，绘制具有PID控制器的控制系统正弦跟踪曲线。3.完成《自动控制系统计算机仿真》教材第7章中的例7-4题。4.完成《自动控制系统计算机仿真》教材第7章中的例7-5题。10实验四系统状态空间设计实验目的1学习系统的能控性、能观测性判别计算方法；2掌握状态反馈极点配置控制系统的设计方法。3掌握状态观察器的设计方法实验原理如果给出了对象的状态方程模型，我们希望引入状态反馈控制