《控制工程基础》实验指导书
陈卫国
湖南文理学院
机械电子工程专业
实验一 典型环节模拟
一、实验目的
①了解、掌握计算机模拟典型环节的基本方法。
②熟悉各种典型环节的阶跃响应曲线。
③了解各种参数变化对典型环节动态特性的影响。
④了解计算机辅助分析和设计的特点与优点
二、实验要求
①通过计算机的仿真图形观测各种典型环节时域响应曲线。
②改变参数，观测参数变化时对典型环节时域响应的影响。
④对实验程序加上注释，写出实验报告。
三、实验内容
一般来讲，线性连续控制系统通常都是由一些典型环节构成的，这些典型环节有比例环节、积分环节、一阶微分环节、惯性环节、振荡环节、延迟环节等。下面分别对其性能进行仿真（建议实验程序在M文件中用单步执行的方式执行程序，以便于分析）：
1）比例环节
比例环节的传递函数为：ksG)(
编程分析当k=1～10时，比例环节在时域的情况：
①当输人信号是单位阶跃信号时，比例环节的输出曲线（单位阶跃响应曲线）是什么形状呢？实验程序如下：
for k=1:1:10
num=k;
den=1;
G=tf(num,den);
step(G);
hold on;
end
在M文件的窗口中，输入程序，录入程序完成后，保存该M文件，在弹出的“保存为”窗工中输人M文件名bl.m（也可以自己取文件名），选择存放该M文件的路径，就可以完成保存工作，然后单击“Tools”菜单中的“Run”，将在step图形窗口中显示出响应图形。试分析系统的输出信号的特点。
2）积分环节
积分环节的传递函数为： TsG1)(
① 当输人信号是单位阶跃信号时，积分环节的仿真程序如下：
num=1;
den=[1,0]; G=tf(num,den);
step(G)
执行程序，试分析系统的输出信号的特点。若G(s)=k/s，编程分析当K=1～10时，在单位阶跃信号激励下，积分环节时域响应的情况。
3）一阶微分环节
一阶微分环节的传递函数为： 1)(TssG
①当输入信号是单位阶跃信号时，一阶微分环节的输出在MATLAB的函数step（）中是无法绘制的。为了能够进行仿真，设置一个极点p，该极点｜P｜＞＞1／T，设极点P=－1000，T取值范围为：1～10。
实验程序如下：
for T=1:1:10;
num=[T,1];
den=[0.0001,1];
G=tf(num,den);
Step(G);
Hols on;
试分析一阶微分环节阶跃响应的特点以及T值的作用。
4）惯性环节
惯性环节的传递函数为:11)(TssG
①当输人信号是单位阶跃信号时，惯性环节的单位阶跃响应曲线是什么形状呢？
实验程序：
for T=1:1:10;
num=[1];
den=[T,1];
G=tf(num,den);
Step(G);
Hols on;
end
分析惯性环节的时间常数T与响应到达稳态值时间之间的关系。
5）振荡环节
振荡环节的传递函数为：121)(22TSSTsG
①当输入信号是单位阶跃信号时，振荡环节的单位阶跃响应曲线是什么形状呢？
实验参考程序如下(其中，T=6，ξ=0.1,0.4,0.7)
for T=6;
for zeta=[0.1,0.4,0.7]; num=[1];
den=[T^2,2*T*zeta,1];
G=tf(num,den);
Step(G);
Hols on;
End
end
分析阻尼系数ξ对单位阶跃响应的影响
实验二 控制系统时域仿真和稳定性研究
凡是能用二阶微分方程描述的控制系统，都称为二阶控制系统。
典型二阶控制系统的闭环传递函数为：2222)(nnnsssG
当0＜ξ＜1时，二阶控制系统被称为欠阻尼系统；当ξ＝1时，称为临界阻尼系统；当ξ＞1时，称为过阻尼系统。
典型二阶控制系统的闭环传递函数有两个可选参数：ξ和ωn，根据不同的ξ和ωn的参数值，对二阶控制系统时域仿真和稳定性进行研究。
一、实验目的
①熟悉二阶控制系统的阶跃响应曲线。
②理解ξ和ωn参数变化对系统动态特性的影响。
二、实验要求
①通过计算机的仿真图形观测二阶控制系统的时域响应曲线。
②改变ξ和ωn，观测参数变化时对典型环节时域响应的影响。
② 对实验程序加上注释，写出实验报告。
三、实验内容
1）二阶控制系统时域响应和稳定性的仿真
在程序中，阻尼比ξ用变量zeta来表示。
当ωn=1，ξ＝0.1，0.3，0.5，0.7，0.9,1.0,2.0时的单位阶跃响应实验程序：
ωn=1;
kos= [0.1:0.2:0.9,1.0,2.0]
for zeta = kos
num=ωn^2;
den=[1,2*zeta*ωn,ωn^2];
G=tf(num,den);
Step(G);
Hols on;
end title(‘step response’)
应用impulse（）函数，同样可以对二阶控制系统进行单位脉冲响应的仿真和分析。
绘制ωn＝6, ξ＝0.7时的单位脉冲响应曲线，比较这些曲线的特点
实验程序：
wn=6
zeta=0.7
figure(1)
hold on
num=wn .^2
den=[1，2*zeta*wn，wn .^2]
impulse(num,den)
title(‘impulse response’)
2）二阶控制系统稳定性的仿真
对二阶系统2222)(nnnsssG， n＝6，绘制分别为0.7, 1.0,2.0情况下系统的零极点图，判定系统的稳定性，并比较极点的分布特征。
实验程序：
wn=6
zeta=[0.7 1.0 2.0]
num=[wn*wn]
for i=1:3
figure(i)
den=[1 2*zeta(i)*wn wn*wn]
pzmap(num den)
grid
end
实验三、计算机辅助频域仿真分析
一、实验目的
①熟悉各种典型环节的频域响应的曲线
②学习采用MATLAB绘制和分析系统频率特性图
二、实验要求
①通过计算机的仿真图形观测各典型环节的频率特性曲线及系统的开环频率特性曲线
②对实验程序加上注释，写出实验报告。
三、实验内容
1） 利用频域分析函数绘制各典型环节的频率特性图
Bode（） 绘制Bode图
Nyquist（） 绘制Nyquist图 2） 已知控制系统的开环传递函数为)1(2)(sssG，绘制系统开环频率特性的极坐标图。
实验程序：
num=[2]
den=[1 1 0]
w=logspace(-1,2)
nyquist(num,den,w)
grid
3） 已知控制系统的开环传递函数为222)11.0()110()1(300)(sssssG，用SIMULINK建模，利用模型文件绘制BODE图、NYQUIST图、脉冲响应、阶跃响应。试建立模型文件和M文件
实验程序：
（1） SIMULINK窗口内构筑已知的开环系统模型，并以文件名“moxig”进行保存。
（2） 建立M文件
[a b c d]=Linmod(‘moxing’)
sys=ss (a,b,c,d)
sys=tf(sys)
figure,bode(sys)
figure,nyquist(sys)
figure,step(sys)
figure,impulse(sys) 附录
Matlab基础
关于MATLAB
MATLAB是Matrix laboratory的缩写，是美国的MATHWORK公司开发的一种进行科学和工程计算的交互式程序语言。1984年推出第一个商业版本，到现在已经到了6.1版本，功能日趋完善和强大。主要适用于矩阵运算及控制和信息处理领域的分析设计。另外还包括仿真、系统辨识、神经元网络、模糊控制等工具箱。其中Simulink仿真工具箱可以利用图形界面和框图来编程和仿真。
运行MATLAB程序的流程
• MATLAB软件通常工作在交互状态下，当键入一条命令后，MATLAB系统立即执行该命令并在屏幕上显示结果；
• MATLAB系统也可以执行储存在文件中的命令序列。这两种工作方式构成了系统的解释环境。
数字运算
MATLAB可以像计算器一样直接进行数学运算。
在MATLAB命令窗口内输入：
80*0.8+16*12.2
ans =
259.2000
帮助命令help
MATLAB提供了大量的函数和命令，如果想记住所有的函数及其调用格式几乎是不可能的。在线帮助可以由help 命令来获得。在MATLAB 的命令窗口直接键入help即可得到所有的帮助主题。如果要对某一命令或函数进行查询，直接在help后跟上该命令或函数即可。
变量的命名规则
变量名对大小写敏感。
变量名的首字符必须是字母。
每个变更名最多可包含19个字符。函数的命名规则与变量相同。
表达式
MTLAB对使用者键入的表达式进行翻译和计算，语句的形式通常为：
变量=表达式；
或简单地写作： 表达式；
其中分号是可选的，如果该语句的最后一个字符是分号，则在执行词句时不在屏幕上显示结果。
