控制工程基础[英]实验
实验一.典型环节的模拟研究：
已知一个小车、倒单摆系统非线性系统方程为：
( 2.92)0.008x x u =-+
20.004sin 36cos n n x θωθωθθ=-+-
其中假设 (0)0;(0)0.2x x ==，
(0)0;(0); 6.781,n θθπω===
（1）要求绘出系统[0,10]t ∈的状态响应曲线
（2）并将上述系统在0θ≈的条件下线性化，并要求绘出线性化后系统
[0,10]t ∈的状态响应曲线，并与非线性系统状态响应曲线相比较。
(1)下面利用Simulink 对该系统进行仿真如下图所示。
图1.倒单摆系统仿真图
在图中已经对主要信号进行了标注下面给出每个未标注信号后加入放大器的增益：
008.092.2=
阶跃K 008
.01
-=一阶微分x K 98.45=二阶微分θK
通过示波器Scope 和Scope1观察x(t)和θ(t)的波形图如下所示。
图2.x(t)波形图3.θ(t)波形（2）将上述系统在0
θ≈的条件下线性化，则方程组改写成如下形式：
( 2.92)0.008
x x u
=-+
2
0.004sin36
n n x
θωθωθ
=-+-
在Simulink中对系统仿真如下所示。
图4.线性化后仿真系统
通过示波器模块可以观察输出信号，图形如下图所示。
图5.x(t)输出波形
图6.θ（t ）输出波形
实验二.典型系统时域响应动、静态性能和稳定性研究； 已知系统的开环传递函数为
2()11G s s s =
++
（1）利用已知的知识判断该开环系统的稳定性（系统的特征方程根、系统零极点表示法）。
（2）判别系统在单位负反馈下的稳定性，并求出闭环系统在[0,10]t ∈内的脉冲响应和单位阶跃响应，分别绘制出相应响应曲线。
（1）该系统的特征方程的根、零极点表示的求解代码如下：
输出结果如下图所示。
图7.特征方程求根结果
图8.零极点分布图
从图中可以看出两个极点在虚轴上，所以该系统处于临界稳定状
态。
（2）求解阶跃响应和脉冲响应的程序代码如下：
该系统的脉冲响应和阶跃响应图如下。
图9.系统脉冲响应和阶跃响应图
实验三.控制系统的频率特性研究; 针对如右图所示系统，已知被控系统G(s)的传递函数为：
2()1G s s s =
+
（1）要求画出该闭环系统的根轨迹曲线图（随K 变化）
（2）求出该闭环系统的传递函数T(s)，并画出当K=2该闭环系统的Bode 图，并分别在图中求出闭环系统的相角裕度与幅值裕度。
（1）由图可以看出系统的特征方程为
s
s K
s KG s ++
==∆21)()( 下面画出随K 变化的闭环系统的根轨迹曲线图，代码如下：
运行程序可以得到根轨迹如下图所示。
图10.随K 变话系统根轨迹图
系统的传递函数为
2
2
)(1)()(2++=+=
s s s KG s KG s T
下面给出计算代码如下：
程序的输出结果和波特图如下所示。
图11.系统求解传递函数输出
从波特图上的信息可以看出：
当sec /73.1rad =ω时，系统有0)(log 20=ωj G ，相角裕度为︒=60..M P 当+∞→ω时，系统有相角为180°，增益裕度为∞=..M G
图12.闭环系统波特图
实验四.线性系统校正.
针对如左图所示系统，已知系
统被控系统G(s)的传递函数为：
2()10510G s s s =
++
（1）当比例控制器增益K=1时，在Matlab/Simulink 中搭建仿真系统，当
输入为阶跃函数时，用示波器Scope 观测系统的输出，绘出响应曲线。并将对象输出y(t)与时钟输出输入Workspace 中，编程求出系统的响应时间、调节时间及超调量，并将计算结果与根据理论公式计算值相比较。
（2）调节控制器增益，使超调量%32%
e<，给出
σ<，且稳态误差0.2
ss
此时K值的范围。
（3）如果想使稳态误差0
e=，应考虑将上述的比例控制器该为何种控
ss
制器，即选用比例积分控制器，或比例微分控制器，或比例、积分与微分控制器，试从超前/滞后校正理论论证，并绘出对应的响应曲线。（1）K=1时在Simulink中搭建的仿真系统如下图所示。
图13.K=1时仿真系统
利用示波器观察输出结果可以看到如下波形
图14.K=1时仿真系统输出波形
下面将输出结果导出到Workspace中，利用To Workspace模块替换Scope将系统输出导入为矩阵y如下图所示，并编程计算相关的参数。
图14.导出系统输出数据
程序的输出结果如下:
系统超调量为pos =0.1188 系统响应时间为tr = 1.6211 系统调节时间为ts = 2.2319 系统稳态误差为ess = 0.5000
通过理论计算算出系统参数如下：
系统超调量为pos =0.1202 系统响应时间为tr = 1.801 系统调节时间为ts = 0.584
（2）调节系统放大器增益，利用上题程序计算参数。
①当K=4时，结果如下：
