实验二
附加题：
1、积分环节
求取
Ts s G 1
)(
在不同积分时间常数T 下的单位阶跃响应，分析积分时间常数的作用。
500
10001500
Step Response
Time (sec)
A m p l i t u d e
由图可看出：积分环节强度随着T 的增加而减小
2、微分环节
在实际系统中，微分环节通常带有惯性，其传递函数为1)(21+=
s T s
T s G ，取T 2＝1，T 1
为不
同数值，分析微分时间常数T 1的作用。
不同微分常数的微分作用
Time (sec)
A m p l i t u d e
由上图可知：微分常数T 对于微分强度成正相关作用
基础题: 1、比例环节
求取K s G =)(在不同比例系数K 下的单位阶跃响应，说明比例系数对系统动态过程的影响。
程序代码： num=[0.5]; den=[1]; hold on;
step(num,den); num=[2]; den=[1]; hold on;
step(num,den); num=[10]; den=[1]; hold on;
step(num,den); 输出结果：
分析：K 越大，系统单位阶跃响应越明显。
2、一阶惯性环节
(1) 求取1
22
)(+=
s s G 的单位脉冲响应，可否用step 命令求取它的脉冲响应？ 122
)(+=
s s G 的单位脉冲响应如下图：
程序代码: num=[2]; den=[2 1]; hold on;
step(num,den);
(2) (求取1
22
)(+=
s s G 的单位脉冲响应，可否用step 命令求取它的脉冲响应？
122
)(+=
s s G 的单位脉冲响应如下图：
num=[2 0]; den=[2 1]; hold on ; step(num,den);
把传递函数乘以s 再求其单位阶跃响应，就可获得乘s 前的传递函数的脉冲响应。如下图
>> num=[0.5]; den=[0.5 1]; hold on;
step(num,den); num=[0.5]; den=[2 1]; hold on;
step(num,den); num=[0.5]; den=[10 1]; hold on;
step(num,den);
K=0.5,T 分别取0.5,2,10，G(s)单位阶跃响应曲线如下
num=[0.5];
den=[2 1];
hold on;
step(num,den);
num=[2];
den=[2 1];
hold on;
step(num,den);
num=[10];
den=[2 1];
hold on;
step(num,den);
T=0.5,K分别取0.5,2,10，G(s)单位阶跃响应曲线如下
2．说明对象模型的特征参数对系统的单位阶跃响应的影响
随着K的增大终值增大为原来的K倍，而调节时间不变。随着T的增大调节时间也随之增大，
但是终值不变。两种情况下系统的稳态误差均为0，不存在超调量，上升时间均趋近于正无穷。由此可以总结出，K 直接影响系统的终值，T 与系统的调节时间紧密相关，且均为正相关。
3、振荡环节（二阶系统）
根据传递函数2
2
22)(n n n
s s s G ωζωω++=的单位阶跃响应。
（1）
n ω=1，ζ分别取０、0.4、1.0、2；
（2）ζ=0.5，
n ω分别取0.2、0.6、1、1.4；
说明这两个特征参数对过渡过程的影响。
t=0:0.01:20 num=[1]; den=[1 0 1]; hold on ; step(num,den); num=[1]; den=[1 0.8 1]; hold on ; step(num,den); num=[1]; den=[1 2 1]; hold on ; step(num,den); num=[1]; den=[1 4 1]; hold on ; step(num,den);
2468101214161820
00.20.40.60.811.2
1.41.61.82
ωn 不变，ζ改变时的系统阶跃响应
Time (sec)
A m p l i t u d e
102030405060
00.2
0.4
0.6
0.8
1
1.2
1.4
ζ不变，ωn 改变时的系统阶跃响应
Time (sec)
A m p l i t u d e
由以上两图和两表中所列数据进行分析可得:
n ω影响二阶系统过渡过程中的峰值时间，过渡时间，上升时间（在ζ不变的情况下，峰值
时间随n ω增大而减小，过渡时间随n ω的增大而减小，上升时间随n ω的增大而减小。
）
ζ影响几乎全部过渡过程指标，其中超调量，衰减比仅与ζ有关（超调量随着ζ的增大而
减小，衰减比随着ζ的增大而增大；在
n ω不变的情况下，峰值时间随ζ增大而增大，过渡
时间随ζ的增大而减小，上升时间随ζ的增大而减小。）
自动控制原理实验 实验报告 实验三闭环电压控制系统研究 学号姓名 时间2014年10月21日 评定成绩审阅教师
实验三闭环电压控制系统研究 一、实验目的： （1）通过实例展示，认识自动控制系统的组成、功能及自动控制原理课程所要解决的问题。 （2）会正确实现闭环负反馈。 （3）通过开、闭环实验数据说明闭环控制效果。 二、预习与回答： （1）在实际控制系统调试时，如何正确实现负反馈闭环？ 答：负反馈闭环，不是单纯的加减问题，它是通过增量法实现的，具体如下： 1.系统开环； 2.输入一个增或减的变化量； 3.相应的，反馈变化量会有增减； 4.若增大，也增大，则需用减法器； 5.若增大，减小，则需用加法器，即。 （2）你认为表格中加1KΩ载后，开环的电压值与闭环的电压值，哪个更接近2V？ 答：闭环更接近。因为在开环系统下出现扰动时，系统前部分不会产生变化。故而系统不具有调节能力，对扰动的反应很大，也就会与2V相去甚远。 但在闭环系统下出现扰动时，由于有反馈的存在，扰动产生的影响会被反馈到输入端，系统就从输入部分产生了调整，经过调整后的电压值会与2V相差更小些。 因此，闭环的电压值更接近2V。 （3）学自动控制原理课程，在控制系统设计中主要设计哪一部份？ 答：应当是系统的整体框架及误差调节部分。对于一个系统，功能部分是“被控对象”部分，这部分可由对应专业设计，反馈部分大多是传感器，因此可由传感器的专业设计，而自控原理关注的是系统整体的稳定性，因此，控制系统设计中心就要集中在整个系统的协调和误差调节环节。 二、实验原理： （1）利用各种实际物理装置（如电子装置、机械装置、化工装置等）在数学上的“相似性”，将各种实际物理装置从感兴趣的角度经过简化、并抽象成相同的数学形式。我们在设计控制系统时，不必研究每一种实际装置，而用几种“等价”的数学形式来表达、研究和设计。又由于人本身的自然属性，人对数学而言，不能直接感受它的自然物理属性，这给我们分析和设计带来了困难。所以，我们又用替代、模拟、仿真的形式把数学形式再变成“模拟实物”来研究。这样，就可以“秀才不出门，遍知天下事”。实际上，在后面的课程里，不同专业的学生将面对不同的实际物理对象，而“模拟实物”的实验方式可以做到举一反三，我们就是用下列“模拟实物”——电路系统，替代各种实际物理对象。
《自动控制原理》 实验报告 姓名： 学号： 专业： 班级： 时段： 成绩： 工学院自动化系
实验一 典型环节的 MATLAB仿真 一、实验目的 1．熟悉MATLAB桌面和命令窗口，初步了解SIMULINK功能模块的使用方法。 2．通过观察典型环节在单位阶跃信号作用下的动态特性，加深对各典型环节响应曲线的理解。 3．定性了解各参数变化对典型环节动态特性的影响。 二、实验原理 1．比例环节的传递函数为 K R K R R R Z Z s G200 , 100 2 ) ( 2 1 1 2 1 2= = - = - = - = 其对应的模拟电路及SIMULINK图形如图1-3所示。 三、实验内容 按下列各典型环节的传递函数，建立相应的SIMULINK仿真模型，观察并记录其单位阶跃响应波形。 ①比例环节1 ) ( 1 = s G和2 ) ( 1 = s G； ②惯性环节 1 1 ) ( 1+ = s s G和 1 5.0 1 ) ( 2+ = s s G ③积分环节 s s G1 ) ( 1 = ④微分环节s s G= ) ( 1 ⑤比例+微分环节（PD）2 ) ( 1 + =s s G和1 ) ( 2 + =s s G ⑥比例+积分环节（PI） s s G1 1 ) ( 1 + =和s s G21 1 ) ( 2 + = 四、实验结果及分析 图1-3 比例环节的模拟电路及SIMULINK图形
① 仿真模型及波形图1)(1=s G 和2)(1=s G ② 仿真模型及波形图11)(1+= s s G 和1 5.01)(2+=s s G 11)(1+= s s G 1 5.01 )(2+=s s G ③ 积分环节s s G 1)(1= ④ 微分环节
相平面法例题解析： 要求： 1.正确求出对于非线性系统在每个线性区的相轨迹方程，也就是e e - 之间关系的方程（或c c - ）。会画相轨迹（模型中是给具体数的）。※※关键是确定开关线方 程。 2. ※※※如果发生自持振荡，计算振幅和周期。 注意相平面法一般应： 1）按照信号流向与传输关系。线性部分产生导数关系，非线性部分形成不同分区。连在一 起就形成了不同线性分区对应的运动方程，即含有c 或者e 的运动方程。 2）※※※根据不同线性分区对应的运动方程的条件方程确定开关线方程。开关线方程确定很关键。 3）※※※根据不同线性分区对应的运动方程，利用解析法（分离变量积分法或者消去t 法） 不同线性分区对应的相轨迹方程，即c c - 和e e - 之间关系。 4）※根据不同分区的初始值绘制出相轨迹，并求出稳态误差和超调、以及自持振荡的周期和振幅等。 例2 问题1. 用相平面法分析系统在输入r (t ) = 4.1(t )时的运动情况。 问题2. 如果发生自持振荡 ，求自持振荡的周期和振幅。 解：问题1：1）设系统结构图，死区特性的表达式： 0,||2 2,22,2x e x e e x e e =≤?? =->??=+<-? 2）线性部分： 2 ()1 ()C s X s s =，则微分方程为：c x = 3）绘制e e - 平面相轨迹图。因为e r c =-，c r e =-，c r e =- ，c r e =- 。代入则 e x r =-+ （1） 当0t >，0r = ，0r = 。代入，则各区的运动方程0,||2I 2,2II 2,2III e e e e e e e e =≤--?? =->---??=--<----? 由于非线性特性有3个分区，相平面e e -分为3个线性区。 注意，当相平面选好后，输入代入后，最后代入非线性特性。 4） 系统开关线：2e =±。 5） 由题意知初始条件(0)(0)(0)4e r c =-=，(0)(0)(0)0e r c =-= 在II 区，则从
自动控制原理 实验报告 姓名学号 时间地点实验楼B 院系专业 实验一系统的数学模 实验二控制系统的时域分析 实验三控制系统的频域分析
实验一系统的数学模 一、实验目的和任务 1、学会使用MATLAB的命令； 2、掌握MATLAB有关传递函数求取及其零、极点计算的函数。 3、掌握用MATLAB 求取系统的数学模型 二、实验仪器、设备及材料 1、计算机 2、MATLAB软件 三、实验原理 1、MATLAB软件的使用 2、使用MATLAB软件在计算机上求取系统的传递函数 四、实验内容 1、特征多项式的建立与特征根的求取 在命令窗口依次运行下面命令，并记录各命令运行后结果 >>p=[1,3,0,4]; p = 1 3 0 4 >>r=roots(p) r = -3.3553 + 0.0000i 0.1777 + 1.0773i 0.1777 - 1.0773i >>p=poly(r) p = 1.0000 3.0000 -0.0000 4.0000 2、求单位反馈系统的传递函数： 在命令窗口依次运行下面命令，并记录各命令运行后结果 >>numg=[1];deng=[500,0,0]; >>numc=[1,1];denc=[1,2]; >>[num1,den1]=series(numg,deng,numc,denc); >>[num,den]=cloop(num1,den1,-1) num = 0 0 1 1
den = 500 1000 1 1 >>printsys(num,den) num/den = s + 1 --------------------------- 500 s^3 + 1000 s^2 + s + 1 3、传递函数零、极点的求取 在命令窗口依次运行下面命令，并记录各命令运行后结果>>num1=[6,0,1];den1=[1,3,3,1]; >>z=roots(num1) ; >>p=roots(den1) ; >>n1=[1,1];n2=[1,2];d1=[1,2*i];d2=[1,-2*i];d3=[1,3]; >>num2=conv(n1,n2) num2 = 1 3 2 >>den2=conv(d1,conv(d2,d3)) den2 = 1 3 4 12 >>printsys(num2,den2) s^2 + 3 s + 2 ---------------------- s^3 + 3 s^2 + 4 s + 12 >>num=conv(num1,den2);den=conv(den1,num2); >>printsys(num,den) 6 s^5 + 18 s^4 + 25 s^3 + 75 s^2 + 4 s + 12 ------------------------------------------- s^5 + 6 s^4 + 14 s^3 + 16 s^2 + 9 s + 2 >>pzmap(num,den),title(‘极点－零点图’)
自动控制原理故事 小时候喜欢看杂书，没什么东西看，不正在文化大革命嘛？不过看进去了两个“化”：机械化和自动化。打小就没有弄明白，这机械化和自动化到底有什么差别，机器不是自己就会动的吗？长大了，总算稍微明白了一点，这机械化是力气活，用机器代替人的体力劳动，但还是要人管着的，不然机器是不知道该干什么不该干什么的；这自动化嘛，就是代替人的重复脑力劳动，是用来管机器的。也就是说，自动化是管着机械化的，或者说学自动化的是管着学机械的……啊，不对，不对，哪是哪啊！ 有人考证古代就有自动化的实例，但现代意义上的自动控制开始于瓦特的蒸汽机。据说纽考门比瓦特先发明蒸汽机，但是蒸汽机的转速控制问题没有解决，弄不好转速飞升，机器会坏不说，还可能说大事故。瓦特在蒸汽机的转轴上安了一个小棍，棍的一端和放汽阀连着，放气阀松开来就关闭，转速增加；按下去阀就打开，转速降低；棍的另一端是一个小重锤，棍中间某个地方通过支点和转轴连接。转轴转起来的时候，小棍由于离心力的缘故挥起来。转速太高了，小棍挥会挥得很高，放汽阀就被按下去打开，转速下降；转速太低了，小棍挥不起来，放汽阀就被松开来关闭，转速回升。这样，蒸汽机可以自动保持稳定的转速，即保证安全，又方便使用。也就是因为这个小小的转速调节器，瓦特的名字和工业革命连在一起，而纽考门的名字就要到历史书里去找了。 类似的例子在机械系统里很多，家居必备的抽水马桶是另一个例子。放水冲刷后，水箱里水位降低，浮子随水面下降，进水阀打开。随着水位的升高，进水阀逐渐关闭，直到水位达到规定高度，进水阀完全关闭，水箱的水正好准备下一次使用。这是一个非常简单但非常巧妙的水位控制系统，是一个经典的设计，但不容易用经典的控制理论来分析，不过这是题外话了. 这些机械系统设计巧妙，工作可靠，实在是巧夺天工。但是在实用中，如果每次都需要这样的创造性思维，那太累，最好有一个系统的方法，可以解决“所有”的自动控制问题，这就是控制理论的由来。 从小大人就教我们，走路要看路。为什么呢？要是不看着路，走路走歪了也不知道，结果就是东撞西撞的。要是看着路呢？走歪了，马上就看到，赶紧调整脚步，走回到正道上来。这里有自动控制里的第一个重要概念：反馈（feedback）。 反馈是一个过程： 1、设定目标，对小朋友走路的例子来说，就是前进的路线。 2、测量状态，小朋友的眼睛看着路，就是在测量自己的前进方向。 3、将测量到的状态和设定的目标比较，把眼睛看到的前进方向和心里想的前进方向作 比较，判断前进方向是否正确；如果不正确，相差有多少。 4、调整行动，在心里根据实际前进方向和设定目标的偏差，决定调整的量。 5、实际执行，也就是实际挪动脚步，重回正确的前进方向。 在整个走路的过程中，这个反馈过程周而复始，不断进行，这样，小朋友就不会走得东倒西歪了。但是，这里有一个问题：如果所有的事情都是在瞬时里同时发生的，那这个反馈过程就无法工作。要使反馈工作，一定要有一定的反应时间。还好，世上之事，都有一个过程，这就为反馈赢得了所需要的时间。 小时候，妈妈在锅里蒸东西，蒸好了，从锅里拿出来总是一个麻烦，需要抹布什么的垫着，免得烫手。但是碗和锅的间隙不大，连手带抹布伸进去颇麻烦，我常常不知天高地厚，自告奋勇地徒手把热的碗拿出来。只要动作快，手起碗落，可以不烫手。当然喽，要是
自动控制原理实验报告 一、实验名称：一、二阶系统的电子模拟及时域响应的动态测试 二、实验目的 1、了解一、二阶系统阶跃响应及其性能指标与系统参数之间的关系 2、学习在电子模拟机上建立典型环节系统模型的方法 3、学习阶跃响应的测试方法 三、实验内容 1、建立一阶系统的电子模型，观测并记录在不同时间常数T时的响应曲线，测定过渡过程时间T s 2、建立二阶系统电子模型，观测并记录不同阻尼比的响应曲线，并测定超调量及过渡过程时间T s 四、实验原理及实验数据 一阶系统 系统传递函数： 由电路图可得，取则K=1，T分别取：0.25, 0.5, 1 T 0.25 0.50 1.00 R2 0.25MΩ0.5M Ω1MΩ C 1μ1μ1μ T S 实测0.7930 1.5160 3.1050 T S 理论0.7473 1.4962 2.9927 阶跃响应曲线图1.1 图1.2 图1.3 误差计算与分析 （1）当T=0.25时，误差==6.12%； （2）当T=0.5时，误差==1.32%； （3）当T=1时，误差==3.58% 误差分析：由于T决定响应参数，而,在实验中R、C的取值上可能存在一定误差，另外，导线的连接上也存在一些误差以及干扰，使实验结果与理论值之间存在一定误差。但是本实验误差在较小范围内，响应曲线也反映了预期要求，所以本实验基本得到了预期结果。 实验结果说明 由本实验结果可看出，一阶系统阶跃响应是单调上升的指数曲线，特征有T确定，T越小，过度过程进行得越快，系统的快速性越好。 二阶系统 图1.1 图1.2 图1.3
系统传递函数： 令 二阶系统模拟线路 0.25 0.50 1.00 R4 210.5 C2 111 实测45.8% 16.9% 0.6% 理论44.5% 16.3% 0% T S实测13.9860 5.4895 4.8480 T S理论14.0065 5.3066 4.8243 阶跃响应曲线图2.1 图2.2 图2.3 注：T s理论根据matlab命令[os,ts,tr]=stepspecs(time,output,output(end),5)得出，否则误差较大。 误差计算及分析 1）当ξ=0.25时，超调量的相对误差= 调节时间的相对误差= 2）当ξ=0.5时，超调量的相对误差==3.7% 调节时间的相对误差==3.4% 4）当ξ=1时，超调量的绝对误差= 调节时间的相对误差==3.46% 误差分析：由于本试验中，用的参量比较多，有R1,R2,R3,R4;C1,C2;在它们的取值的实际调节中不免出现一些误差，误差再累加，导致最终结果出现了比较大的误差，另外，此实验用的导线要多一点，干扰和导线的传到误差也给实验结果造成了一定误差。但是在观察响应曲线方面，这些误差并不影响，这些曲线仍旧体现了它们本身应具有的特点，通过比较它们完全能够了解阶跃响应及其性能指标与系统参数之间的关系，不影响预期的效果。 实验结果说明 由本实验可以看出，当ωn一定时，超调量随着ξ的增加而减小，直到ξ达到某个值时没有了超调；而调节时间随ξ的增大，先减小，直到ξ达到某个值后又增大了。 经理论计算可知，当ξ=0.707时，调节时间最短，而此时的超调量也小于5%，此时的ξ为最佳阻尼比。此实验的ξ分布在0.707两侧，体现了超调量和调节时间随ξ的变化而变化的过程，达到了预期的效果。 图2.2 图2.1 图2.3
实验一MATLAB 仿真基础 一、实验目的： （1）熟悉MATLAB 实验环境，掌握MATLAB 命令窗口的基本操作。 （2）掌握MATLAB 建立控制系统数学模型的命令及模型相互转换的方法。 （3）掌握使用MATLAB 命令化简模型基本连接的方法。 （4）学会使用Simulink 模型结构图化简复杂控制系统模型的方法。 二、实验设备和仪器 1．计算机；2. MATLAB 软件 三、实验原理 函数tf ( ) 来建立控制系统的传递函数模型，用函数printsys ( ) 来输出控制系统的函数，用函数命令zpk ( ) 来建立系统的零极点增益模型，其函数调用格式为：sys = zpk ( z, p, k )零极点模型转换为多项式模型[num , den] = zp2tf ( z, p, k ) 多项式模型转化为零极点模型 [z , p , k] = tf2zp ( num, den ) 两个环节反馈连接后，其等效传递函数可用feedback ( ) 函数求得。 则feedback （）函数调用格式为： sys = feedback （sys1, sys2, sign ） 其中sign 是反馈极性，sign 缺省时，默认为负反馈，sign ＝-1；正反馈时，sign ＝1；单位反馈时，sys2＝1，且不能省略。 四、实验内容： 1.已知系统传递函数，建立传递函数模型 2.已知系统传递函数，建立零极点增益模型 3．将多项式模型转化为零极点模型 1 2s 2s s 3s (s)23++++=G )12()1()76()2(5)(332 2++++++= s s s s s s s s G 12s 2s s 3s (s)23++++= G )12()1()76()2(5)(3322++++++=s s s s s s s s G
自动控制原理概述及开闭环实例分析 摘要 本文简单介绍了自动控制的基本原理和发展概况，并从开环控制和闭环控制两方面对自动控制原理进行了详细介绍。列举了开环控制和闭环控制的几个实例，结合实例分析了开环控制和闭环的优缺点，并对两种控制方式进行了对比。 关键词：自动控制、基本原理、开环、闭环 1自动控制基本原理及发展概述 所谓的自动控制，就是指在没有人直接参与的情况下，利用外加的设备(称为控制器)操作被控对象（如机器、设备或生产过程）的某个状态或参数（称为被控量），使其按预先设定的规律自动运行。 一般情况下自动控制理论的发展过程可以分为以下三个阶段： 1.1经典控制理论时期 时间为20世纪40-60年代，经典控制理论主要是解决单输入单输出问题，主要采用以传递函数、频率特性、根轨迹为基础的频域分析方法。此阶段所研究的系统大多是线性定常系统，对非线性系统，分析时采用的相平面法一般不超过两个变量。 1.2现代控制理论时期 时间为20世纪60-70年代，这个时期由于计算机的飞速发展，推动了空间技术的发展。经典控制理论中的高阶常微分方程可以转化为一阶微分方程组，用以描述系统的动态过程，这种方法可以解决多输入多输出问题，系统既可以是线性的、定常的，也可以是非线性的、时变的。 1.3大系统理论、智能控制理论时期 时间为20世纪70年代末至今，控制理论向着“大系统理论”和“智能控制”方向发展。“大系统理论”是用控制和信息的观点，研究各种大系统的结构方案、总体设计中的分解方法和协调等问题的技术理论基础。而“智能控制”是研究与模拟人类智能活动及其控制与信息传递过程的规律，研究具有某些仿人智能的工程控制与信息处理系统。 2自动控制系统分类 按照控制方式和策略，系统可分为开环控制系统和闭环控制系统。 2.1开环控制系统 开环控制系统是一种简单的控制系统，在控制器和控制对象间只有正向控制作用，系统的输出量不会对控制器产生任何影响，如图1所示。在该类控制系统中，对于每一个输入量，就有一个与之对应的工作状态和输出量，系统的精度仅取决于元件的精度和执行机构的调整精度。 控制量输出量 图1 开环控制系统
上海电力学院实验报告 自动控制原理实验课程 题目：2.1.1(2.1.6课外)、2.1.4（2.1.5课内）班级：gagagagg 姓名：lalalal 学号：hahahahah 时间：zzzzzzzzzzz
实验内容一： 一、问题描述： 已知系统结构图，（1）用matlab编程计算系统的闭环传递函数；（2）用matlab转换函数表示系统状态空间模型；（3）计算其特征根。 二、理论方法分析 （1）根据系统结构图的串并联关系以及反馈关系，分别利用tf （）函数series（）函数，parallel函数以及feedback函数构建系统传递函数；（2）已求出系统传递函数G，对于线性定常系统利用函数ss（G）课得到系统的状态空间模型。（3）利用线性定常系统模型数据还原函数[num，den]=tfdata（G，‘v’）可得到系统传递函数的分子多项式num与分母多项式den，利用roots（den）函数可得到系统的特征根。 三、实验设计与实现 新建M文件，编程程序如下文所示： G1=tf([0.2],[1,1,1]); G2=tf([0.3],[1,1]); G3=tf([0.14],[2,1]); G4=series(G2,G3);%G2与G3串联 G5=0.7*feedback(G4,-1,1); G6=0.4*feedback(G1,G5,1); G7=feedback(G6,0.6)
ss(G7)%将系统传递函数转化为状态空间模型 [num den]=tfdata(G7,'v');%还原系统传递函数分子、分母系数矩阵 roots(den)%求系统传递函数特征根 点击Run运行 四、实验结果与分析 M文件如下： 运行结果如下：
实验一 MATLAB 仿真基础 、实验目的: (1) 熟悉MATLAB 实验环境，掌握MATLAB 命令窗口的基本操作。 (2) 掌握MATLAB 建立控制系统数学模型的命令及模型相互转换的方法。 (3) 掌握使用MATLAB 命令化简模型基本连接的方法。 (4) 学会使用Simulink 模型结构图化简复杂控制系统模型的方法。 二、实验设备和仪器 1 ?计算机；2. MATLAB 软件 三、实验原理 函数tf ()来建立控制系统的传递函数模型，用函数printsys ()来输出控制系 统的函数，用函数命令zpk ()来建立系统的零极点增益模型，其函数调用格式 为：sys = zpk ( z, p, k 零极点模型转换为多项式模型[num , den] = zp2tf ( z, p, k ) 多项式模型转化为零极点模型 [z , p , k] = tf2zp ( num, den ) 两个环节反馈连接后，其等效传递函数可用 feedback ()函数求得。 则 feedback ()函数调用格式为： sys = feedback (sysl, sys2, sigh 其中sign 是反馈极性，sign 缺省时，默认为负反馈，sign = -1；正反馈时, sig n = 1;单位反馈时，sys2= 1，且不能省略。 四、实验内容： 1. 已知系统传递函数，建立传递函数模型 2 2 5(s 2) (s 6s 7) 3 3 s(s 1) (s 2s 1) 2. 已知系统传递函数，建立零极点增益模型 s 3 飞 2~ s 2s 2s 1 3 ?将多项式模型转化为零极点模型 5(s 2)2(s 2 6s 7) G(s) s 3 s 3 2s 2 2s 1 G(s) G(s)
《自动控制原理》课程实验报告 实验名称系统根轨迹分析 专业班级 *********** ********* 学 号 姓名** 指导教师李离 学院名称电气信息学院 2012 年 12 月 15 日
一、实验目的 1、掌握利用MATLAB 精确绘制闭环系统根轨迹的方法； 2、了解系统参数或零极点位置变化对系统根轨迹的影响； 二、实验设备 1、硬件：个人计算机 2、软件：MATLAB 仿真软件（版本6.5或以上） 三、实验内容和步骤 1．根轨迹的绘制 利用Matlab 绘制跟轨迹的步骤如下： 1） 将系统特征方程改成为如下形式：1 + KG ( s ) = 1 + K ) () (s q s p =0， 其中，K 为我们所关心的参数。 2） 调用函数 r locus 生成根轨迹。 关于函数 rlocus 的说明见图 3.1。 不使用左边的选项也能画出根轨迹，使用左边的选项时，能 返回分别以矩阵和向量形式表征的特征根的值及与之对应的增益值。 图3.1 函数rlocus 的调用 例如，图 3.2 所示系统特征根的根轨迹及其绘制程序见图 3.3。
图3.2 闭环系统一 图3.3 闭环系统一的根轨迹及其绘制程序
图 3.4 函数 rlocfind 的使用方法 注意：在这里，构成系统 s ys 时，K 不包括在其中，且要使分子和分母中 s 最高次幂项的系数为1。 当系统开环传达函数为零、极点形式时，可调用函数 z pk 构成系统 s ys : sys = zpk([zero],[pole],1); 当系统开环传达函数无零点时，[zero]写成空集[]。 对于图 3.2 所示系统， G(s)H(s)= )2()1(++s s s K *11+s =) 3)(2() 1(+++s s s s K . 可如下式调用函数 z pk 构成系统 s ys : sys=zpk([-1],[0 -2 -3],1) 若想得到根轨迹上某个特征根及其对应的 K 的值，一种方法是在调用了函数 rlocus 并得到了根 轨迹后调用函数 rlocfind 。然后，将鼠标移至根轨迹图上会出现一个可移动的大十字。将该十字的 中心移至根轨迹上某点，再点击鼠标左键，就可在命令窗口看到该点对应的根值和 K 值了。另外一种 较为方便的做法是在调用了函数 rlocus 并得到了根轨迹后直接将鼠标移至根轨迹图中根轨迹上某
自动控制原理 实验报告 实验一典型系统的时域响应和稳定性分析 (2) 一、实验目的 (3) 二、实验原理及内容 (3) 三、实验现象分析 (5) 方法一：matlab程序 (5) 方法二：multism仿真 (12)
方法三：simulink仿真 (17) 实验二线性系统的根轨迹分析 (21) 一、确定图3系统的根轨迹的全部特征点和特征线，并绘出根轨迹 (21) 二、根据根轨迹图分析系统的闭环稳定性 (22) 三、如何通过改造根轨迹来改善系统的品质？ (25) 实验三线性系统的频率响应分析 (33) 一、绘制图1. 图3系统的奈氏图和伯德图 (33) 二、分别根据奈氏图和伯德图分析系统的稳定性 (37) 三、在图4中，任取一可使系统稳定的R值，通过实验法得到对应的伯德图，并据此导 出系统的传递函数 (38) 实验四、磁盘驱动器的读取控制 (41) 一、实验原理 (41) 二、实验内容及步骤 (41) （一）系统的阶跃响应 (41) （二) 系统动态响应、稳态误差以及扰动能力讨论 (45) 1、动态响应 (46) 2、稳态误差和扰动能力 (48) （三）引入速度传感器 (51) 1. 未加速度传感器时系统性能分析 (51) 2、加入速度传感器后的系统性能分析 (59) 五、实验总结 (64) 实验一典型系统的时域响应和稳定性分 析
一、 实验目的 1．研究二阶系统的特征参量（ξ、ωn ）对过渡过程的影响。 2．研究二阶对象的三种阻尼比下的响应曲线及系统的稳定性。 3．熟悉Routh 判据，用Routh 判据对三阶系统进行稳定性分析。 二、 实验原理及内容 1．典型的二阶系统稳定性分析 (1) 结构框图：见图1 图1 (2) 对应的模拟电路图 图2 (3) 理论分析 导出系统开环传递函数，开环增益0 1 T K K = 。 (4) 实验内容 先算出临界阻尼、欠阻尼、过阻尼时电阻R 的理论值，再将理论值应用于模拟电路中，观察二阶系统的动态性能及稳定性，应与理论分析基本吻合。在此实验中(图2)， s 1T 0=， s T 2.01=，R 200 K 1= R 200 K =?
自动控制原理实验指导书 王娜编写 电气工程与自动化学院 自动化系 2017年11月 实验一控制系统的时域分析
[实验目的] 1、熟悉并掌握Matlab 操作环境和基本方法，如数据表示、绘图等命令； 2、掌握控制信号的拉氏变换与反变换laplace 和ilaplace ，控制系统生成模型的常用函数命令sys=tf(num,den)，会绘制单位阶跃、脉冲响应曲线； 3、会构造控制系统的传递函数、会利用matlab 函数求取系统闭环特征根； 4、会分析控制系统中n ζω， 对系统阶跃、脉冲响应的影响。 [实验内容及步骤] 1、矩阵运算 a) 构建矩阵：A=[1 2;3 4]; B=[5 5;7 8]; 解： >> A=[1 2;3 4] A = 1 2 3 4 >>B=[5 5;7 8] B = 5 5 7 8 b) 已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ，求矩阵A 的特征值、特征多项式和特征向量. 解：>> A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4]; >> [V ,D]=eig(A) V = 0.4181 -0.4579 - 0.3096i -0.4579 + 0.3096i -0.6044 0.6211 -0.1757 + 0.2740i -0.1757 - 0.2740i 0.0504 0.5524 0.7474 0.7474 -0.2826 0.3665 -0.1592 - 0.0675i -0.1592 + 0.0675i 0.7432 D = 13.0527 0 0 0 0 -4.1671 + 1.9663i 0 0 0 0 -4.1671 - 1.9663i 0 0 0 0 2.1815 >> p=poly(A) p = -6.9000 -77.2600 -86.1300 604.5500 2. 基本绘图命令 a) 绘制余弦曲线y=cos(x)，x ∈[0，2π] 解：>> x=linspace(0,2*pi); >> y=cos(x); >> plot(x,y)
自动控制原理实验 实验一 典型环节的电模拟及其阶跃响应分析 一、实验目的 ⑴ 熟悉典型环节的电模拟方法。 ⑵ 掌握参数变化对动态性能的影响。 二、实验设备 ⑴ CAE2000系统（主要使用模拟机，模/数转换，微机，打印机等）。 ⑵ 数字万用表。 三、实验内容 １．比例环节的模拟及其阶跃响应 微分方程 )()(t Kr t c -= 传递函数 = )(s G ) () (s R s C K -= 负号表示比例器的反相作用。模拟机排题图如图9-1所示，分别求取K=1，K=2时的阶跃响应曲线，并打印曲线。 图9-1 比例环节排题图 图9-2 积分环节排题图 ２．积分环节的模拟及其阶跃响应 微分方程 )() (t r dt t dc T = 传递函数 s K Ts s G ==1)( 模拟机排题图如图9-2所示，分别求取K=1，K=0.5时的阶跃响应曲线，并打印曲线。 ３．一阶惯性环节的模拟及其阶跃响应 微分方程 )()() (t Kr t c dt t dc T =+ 传递函数 1 )(+=TS K S G 模拟机排题图如图３所示，分别求取K=1, T=1; K=1, T=2; K=2, T=2 时的阶跃
响应曲线，并打印曲线。 ４．二阶系统的模拟及其阶跃响应 微分方程 )()() (2)(2 22 t r t c dt t dc T dt t c d T =++ξ 传递函数 121 )(22++=Ts s T s G ξ2 2 2 2n n n s s ωξωω++= 画出二阶环节模拟机排题图，并分别求取打印： ⑴ T=1，ξ=0.1、0.5、1时的阶跃响应曲线。 ⑵ T=2，ξ=0.5 时的阶跃响应曲线。 四、实验步骤 ⑴ 接通电源，用万用表将输入阶跃信号调整为2V 。 ⑵ 调整相应系数器；按排题图接线，不用的放大器切勿断开反馈回路（接线时，阶跃开关处于关断状态）；将输出信号接至数/模转换通道。 ⑶ 检查接线无误后，开启微机、打印机电源；进入CAE2000软件，组态A/D ，运行实时仿真；开启阶跃输入信号开关，显示、打印曲线。 五．实验预习 ⑴ 一、二阶系统的瞬态响应分析；模拟机的原理及使用方法（见本章附录）。 ⑵ 写出预习报告；画出二阶系统的模拟机排题图；在理论上估计各响应曲线。 六．实验报告 ⑴ 将每个环节的实验曲线分别整理在一个坐标系上，曲线起点在坐标原点上。分析各参数变化对其阶跃响应的影响，与估计的理论曲线进行比较，不符请分析原因。 ⑵ 由二阶环节的实验曲线求得σ﹪、t s 、t p ，与理论值进行比较，并分析σ﹪、t s 、t p 等和T 、ξ的关系。 实验二 随动系统的开环控制、闭环控制及稳定性 一．实验目的 了解开环控制系统、闭环控制系统的实际结构及工作状态；控制系统稳定的概念以及系统开环比例系数与系统稳定性的关系。 二．实验要求 能按实验内容正确连接实验线路，正确使用实验所用测试仪器，在教师指导下独立
一.试建立如图所示电路的动态微分方程，并求传递函数。 解：1、建立电路的动态微分方程 根据KCL 有 2 00i 10i ) t (u )]t (u )t (d[u )t (u )t (u R dt C R =-+- 即 )t (u ) t (du )t (u )()t (du i 2i 21021021R dt C R R R R dt C R R +=++ 2、求传递函数 对微分方程进行拉氏变换得 )(U )(U )(U )()(U i 2i 21021021s R s Cs R R s R R s Cs R R +=++ 得传递函数 2 1212 21i 0)(U )(U )(R R Cs R R R Cs R R s s s G +++== 二．结构图化简。 解：1.将-G 3G 4H 4回路的综合点前移至①处 2.将反馈通道H 1的综合点也前移至①处 3.将以上两条反馈通道利用并联公式合并 4.从内到外化简回路 5.最后得到系统的传递函数为 1 432144323213324 3211)(H G G G G H G G H G G G H G G G G G G s G -+++= （注：由于电脑画图太慢，这里只给出了文字步骤，解题时请同学们用图示把各个化简步骤列写出来。）
三．系统结构图如图所示： 1、写出闭环传递函数() ()() C s s R s Φ= 表达式； 2、要使系统满足条件：707.0=ξ,2=n ω,试确定相应的参数K 和β； 3、求此时系统的动态性能指标s t ,0 0σ ； 4、t t r 2)(=时，求系统由()r t 产生的稳态误差ss e ； 解：1. 2 2222 221)()()(n n n s s K s K s K s K s K s K s R s C s ωξωωββ++=++=++==Φ 2. ? ??=====2224222 n n K K ξωβω ? ??==707 .04βK 3. 0010 32.42 ==--ξξπ σ e 12.22 3 3 == = n s t ξω 4. )1(1)(1)(2+= +=+ =s s K s s K s K s K s G βββ ? ??==11v K K β 414.12=== βK ss K A e （注：考虑再三，还是决定把时域性能指标（超调量、峰值时间、调节时间）的公式告诉大家，到时会写在黑板上。除此之外，拉氏变换表也会给大家。如有需要，也会给入射角和出射角公式。其他公式需要自行记忆。）
我们现实生活中的许许多多的事情都可以运用《自动控制原理》中的一些理论和原理来进行有效的控制和调节，进而使我们的生活更加科学合理。 就拿我们大学生的生活费用来说吧，我们就可以运用自动控制原理中的反馈调节理论对我们的生活费用进行合理有效的规划和利用，从而达到感性的利用和管理我们的生活费用。我们也看到了，我们身边的很多同学对生活费的使用是毫无计划可言的，即他们都喜欢即兴花费（今朝有酒今朝醉）。因此，往往使很多同学到学期末的时候开始出现“财政赤字”现象，导致四处借钱来维系生活。其实，与其这样还不如在每学期开始，或者每个月开始就对我们的生活费用列一项详尽的计划，日后可以照着计划开控制自己的生活费，杜绝寅吃卯粮的现象。这里就要用到自动控制原理，我们可以把钱包当作控制对象，每个月计划花费的生活费作为控制系统的输入，每天的实际花费为输出。我们对每月的花费做一个详细的安排。如用餐费、生活平花费、娱乐花费等都做一个具体的切实的计划。这就相当于我们控制系统里的控制器，即系统的核心，而控制器性能的好坏直接影响了我们控制的鲜果。对应于我们这个生活费花费系统中则应该表现在计划的具体内容、细节以及具体实现方面。计划定制的越详细具体切实，并且贯彻落实我们的计划，则控制就回去得越好的效果。在这过程中，我们自己本身就相当于系统中的执行器，我们根据计划的具体要求，具体实施我们的花费，决定是让钱包常开还是常闭，或者是快开等等的问题。这样我们就可以利用控制系统的响应原理对生活费进行了一
定的控制，而不是像有些同学那样毫无计划的挥霍。当然了，这只是相当于一个开环系统，因为过程中并没有对实际花费尽享有效的检测，控制效果无法达到理想状况，而且，有时候生活中还会出发一些意外情况，这也是避免不了的额外花费，我们姑且就叫做是扰动因素吧。我们都知道，这些扰动因素还是会影响系统的性能以及控制效果。但是，如果我们在我们上述的开环系统的基础上增加一个测量边送环节，使系统变成闭环系统，那样我们的控制效果和精度将会更加显著。而这个测量变送环节我们可以通过实际的方式（如每一天的花费进行记账，一周进行一次合计等）实现。这么我们可以调整和更改我们的计划使控制更加合理。 除此之外，运用控制原理中的一些内容方法和理论来指导我们更加美好科学的生活的例子还有很多很多。总之，我要我们学以致用，把我们所学的知识和理论云涌到我们的生活中，我相信我们的生活会更加美好！
实验二二阶系统的动态过程分析 一、 实验目的 1. 掌握二阶控制系统的电路模拟方法及其动态性能指标的测试技术。 2. 定量分析二阶系统的阻尼比ξ和无阻尼自然频率n ω对系统动态性能的影响。 3. 加深理解“线性系统的稳定性只与其结构和参数有关，而与外作用无关”的 性质。 4. 了解和学习二阶控制系统及其阶跃响应的Matlab 仿真和Simulink 实现方法。 二、 实验内容 1. 分析典型二阶系统()G s 的ξ和n ω变化时，对系统的阶跃响应的影响。 2. 用实验的方法求解以下问题： 设控制系统结构图如图2.1所示，若要求系统具有性能： %20%,1,p p t s σσ=== 试确定系统参数K 和τ，并计算单位阶跃响应的特征量d t ，r t 和s t 。 图2.1 控制系统的结构图 3. 用实验的方法求解以下问题： 设控制系统结构图如图2.2所示。图中，输入信号()r t t θ=，放大器增益A K 分别取13.5，200和1500。试分别写出系统的误差响应表达式，并估算其性能指标。 图2.2 控制系统的结构图
三、实验原理 任何一个给定的线性控制系统，都可以分解为若干个典型环节的组合。将每个典型环节的模拟电路按系统的方块图连接起来，就得到控制系统的模拟电路图。 通常，二阶控制系统 2 22 () 2 n n n G s s ω ξωω = ++ 可以分解为一个比例环节、一个惯性环节和一个积分环节，其结构原理如图 2.3所示，对应的模拟电路图如图2.4所示。 图2.3 二阶系统的结构原理图 图2.4 二阶系统的模拟电路原理图 图2.4中：()(),()() r c u t r t u t c t ==-。 比例常数（增益系数）2 1 R K R =，惯性时间常数 131 T R C =，积分时间常数242 T R C =。其闭环传递函数为： 12 2 21 112 () 1 ()(1) c r K U s TT K K U s T s T s K s s T TT == ++++ (0.1) 又：二阶控制系统的特性由两个参数来描述，即系统的阻尼比ξ和无阻尼自然频
实验一 典型环节的模拟研究及阶跃响应分析 1、比例环节 可知比例环节的传递函数为一个常数： 当Kp 分别为0.5，1，2时，输入幅值为1.84的正向阶跃信号，理论上依次输出幅值为0.92，1.84，3.68的反向阶跃信号。实验中，输出信号依次为幅值为0.94，1.88，3.70的反向阶跃信号， 相对误差分别为1.8%,2.2%,0.2%. 在误差允许范围内可认为实际输出满足理论值。 2、 积分环节 积分环节传递函数为： （1）T=0.1(0.033)时，C=1μf (0.33μf )，利用MATLAB ，模拟阶跃信号输入下的输出信号如图： T=0.1 T=0.033 与实验测得波形比较可知，实际与理论值较为吻合，理论上T=0.033时的波形斜率近似为T=0.1时的三倍，实际上为8/2.6=3.08，在误差允许范围内可认为满足理论条件。 3、 惯性环节 惯性环节传递函数为： i f i o R R U U -=TS 1 CS R 1Z Z U U i i f i 0-=-=-=1 TS K )s (R )s (C +-=
K = R f /R 1,T = R f C, （1） 保持K = R f /R 1 = 1不变，观测T = 0.1秒，0.01秒（既R 1 = 100K,C = 1μf ， 0.1μf ）时的输出波形。利用matlab 仿真得到理论波形如下： T=0.1时 t s （5%）理论值为300ms,实际测得t s =400ms 相对误差为：（400-300）/300=33.3%,读数误差较大。 K 理论值为1，实验值2.12/2.28， 相对误差为（2.28-2.12）/2.28=7%与理论值较 为接近。 T=0.01时 t s （5%）理论值为30ms,实际测得t s =40ms 相对误差为：（40-30）/30=33.3% 由于ts 较小，所以读数时误差较大。 K 理论值为1，实验值2.12/2.28， 相对误差为（2.28-2.12）/2.28=7%与理论值较为接近 （2） 保持T = R f C = 0.1s 不变，分别观测K = 1，2时的输出波形。 K=1时波形即为（1）中T0.1时波形 K=2时，利用matlab 仿真得到如下结果： t s （5%）理论值为300ms,实际测得t s =400ms 相对误差为：（400-300）/300=33.3% 读数误差较大 K 理论值为2，实验值4.30/2.28， 相对误差为（2-4.30/2.28）/2=5.7% 与理论值较为接近。