控制工程基础实验报告

1、控制工程基础实验报告北京工业大学机电学院指导教师：初红艳 学号： 姓名： 一实验目的本实验中，学生使用MATLAB语言进行控制系统的分析，可以达到以下目的：（1）通过MATLAB的分析，掌握控制系统的时域瞬态响应、频率特性，根据时域性能指标、频域性能指标评价控制系统的性能，根据系统频率特性进行稳定性分析，了解对系统进行校正的方法，从而进一步巩固、加深对课堂内容的掌握，加强对控制工程基础知识的掌握。（2）熟悉MATLAB的控制系统图形输入与仿真工具SIMULINK，能够对一些框图进行仿真或线性分析，使一个复杂系统的输入变得相当容易且直观。（3）通过本实验，使学生掌握进行控制系统计算机辅助分析的方

2、法，学会利用MATLAB语言进行复杂的实际系统的分析、校正与设计，具备解决工程实际问题的能力。二实验内容控制系统方块图如图1所式。这是一个电压转角位置随动系统，系统的功能是用电压量去控制一个设备的转角，给定值大，输出转角也就成比例地增大。+ + - - 图1 系统方块图图中，为前置放大及校正网络传递函数 为功率放大器放大倍数， 为电动机传递系数， 为电动机机电时间常数， 为电动机电磁时间常数， 为测速传递系数， 为测速反馈分压系数， 为主反馈电位计传递系数， 为输入电压 为反馈电压 为速度环输入电压 为测速机电压为电动机电压为电动机转速取 3、 实验报告1.对于二阶系统：， 其阶跃响应和单位脉

3、冲响应分别如图1-1、1-2所示： MATLAB语言为：>>num=2.83num = 2.8300>> den=0.0004,0.1,1den = 0.0004 0.1000 1.0000>> sys=tf(num,den)sys = 2.83 - 0.0004 s2 + 0.1 s + 1 Continuous-time transfer function.>> step(sys)>> impulse(sys)图1-1图1-2此时阻尼比为2.5，由其阶跃响应可知其稳态值为2.83，为过阻尼状态，瞬态响应指标如上图所示；在无阻尼自振

4、角频率不变时，通过简单计算得出如下的结论： 调整TM=0.01，Ta=0.04,使得系统处于欠阻尼状态，,其阶跃响应与单位脉冲响应1-3、1-4；MATLAB语言如下： >>num=2.83num = 2.8300>> den=0.0004,0.01,1den = 0.0004 0.0100 1.0000>> sys=tf(num,den)sys = 2.83 - 0.0004 s2 + 0.01 s + 1 Continuous-time transfer function.>> step(sys)>> impulse(sys) 其

5、瞬态响应指标在下图中标出; 图1-3图1-4调整TM=0.04，Ta=0.01,使得系统处于临界阻尼状态，其阶跃响应与单位脉冲响应如图1-5、1-6所示；MATLAB语言如下： >>num=2.83num = 2.8300>> den=0.0004,0.04,1den = 0.0004 0.0400 1.0000>> sys=tf(num,den)sys = 2.83 - 0.0004 s2 + 0.04 s + 1 Continuous-time transfer function.>> step(sys)>> impulse(sy

6、s图1-5图1-6由以上当TM=0.1，Ta=0.0004，系统处于过阻尼状态，其阶跃响应与单位脉冲响应见图1-1、1-2所示：分析：由以上响应曲线和响应指标可知，过阻尼与临界阻尼无超调、无振荡，而欠阻尼有超调和振荡，过阻尼达到平衡状态所需的时间比临界阻尼和过阻尼都要长。从快速、稳定和准确综合考虑，上诉系统性能最佳的应该是欠阻尼，即TM=0.01，Ta=0.04。对于实际的二阶振荡系统而言，考虑到系统的稳定、快速和准确性，在无阻尼自振角频率不变的情况下，应该让系统处于欠阻尼状态，即满足，且应该让尽可能大，这样系统性能最佳。2.对此系统的小闭环系统；求得：其开环传递函数为： 闭环传递函数为： 开

7、环系统的乃氏图和伯德图分别如图1-7、图1-8所示；MATLAB语言如下：>>num=32.545num = 32.5450>>den=0.0004,0.1,1den = 0.0004 0.1000 1.0000>>sys=tf(num,den)sys = 32.55 - 0.0004 s2 + 0.1 s + 1 Continuous-time transfer function.>> nyquist(sys)>> bode(sys)>> gm,pm,wcg,wcp=margin(sys)gm = Infpm = 47.

8、4634wcg = Infwcp = 239.8018图1-7图1-8根据乃奎斯特稳定判据：该系统开环传递函数无右特征根，其乃氏图不包围（-1，j0）点，故闭环系统稳定；其频域指标为（幅值裕量不是以dB为单位）：故其闭环系统稳定。3.当时，整个系统的开环传递函数如下：闭环传递函数为：开环传递函数的乃氏图和伯德图分别为图1-9、1-10所示；MATLAB语言如下;>> num=133.01num = 133.0100>> den=0.0004,0.1,33.545,0den = 0.0004 0.1000 33.5450 0>> sys=tf(num,den)

9、sys = 133 - 0.0004 s3 + 0.1 s2 + 33.55 s Continuous-time transfer function.>> nyquist(sys)>> bode(sys)>> gm,pm,wcg,wcp=margin(sys)gm = 63.0498pm = 89.3226wcg = 289.5902wcp = 3.9656图1-9 图1-10根据乃奎斯特稳定判据：该系统开环传递函数无右特征根，其乃氏图不包围（-1，j0）点，故闭环系统稳定；其频域指标为（幅值裕量不是以dB为单位）：故其闭环系统稳定。4.当时，整个系统开环传

10、递函数为=0.1，=0.3的乃氏图分别见图1-11、1-12所示；MATLAB语言如下;=0.1时:>> num1=133.01num1 = 133.0100>> den1=0.0004,0.1,33.545,0den1 = 0.0004 0.1000 33.5450 0>> num2,den2=pade(0.1,3)num2 = 1.0e+05 * -0.0000 0.0012 -0.0600 1.2000den2 = 1.0e+05 *0.0000 0.0012 0.0600 1.2000>>num=conv(num1,num2)num =

11、1.0e+07 * -0.0000 0.0016 -0.0798 1.5961>> den=conv(den1,den2)den = 1.0e+06 * 0.0000 0.0000 0.0000 0.0047 0.2133 4.0254 0>> sys=tf(num,den)sys = -133 s3 + 1.596e04 s2 - 798060 s + 1.596e07 - 0.0004 s6 + 0.148 s5 + 47.95 s4 + 4673 s3 + 213270 s2 + 4.025e06 s Continuous-time transfer functi

12、on.>> nyquist(sys)=0.3时：>> num1=133.01num1 = 133.0100>> den1=0.0004,0.1,33.545,0den1 = 0.0004 0.1000 33.5450 0>> num2 den2=pade(0.3,3)num2 = 1.0e+03 * -0.0010 0.0400 -0.6667 4.4444den2 = 1.0e+03 * 0.0010 0.0400 0.6667 4.4444>> num=conv(num1,num2)num = 1.0e+05 * -0.0013

13、0.0532 -0.8867 5.9116>> den=conv(den1,den2)den = 1.0e+05 * 0.0000 0.0000 0.0004 0.0141 0.2281 1.4909 0>> sys=tf(num,den)sys = -133 s3 + 5320 s2 - 8.867e04 s + 5.912e05 - 0.0004 s6 + 0.116 s5 + 37.81 s4 + 1410 s3 + 2.281e04 s2 + 1.491e05 s Continuous-time transfer function.>> nyquis

14、t(sys)图1-11图1-12根据那奎斯特稳定判据：该系统开环传递函数无右特征根，由上图可知=0.1，=0.3时其开环传递函数乃氏图均不包围（-1，j0）点，故=0.1，=0.3时其闭环都系统稳定；从图1-9、1-11和1-12的对比我们可以直观观察的得出增加一个延时环节会使系统的乃氏图向左移动，且越大，向左移动距离越大。增加延时环节会使系统乃氏图左移动，在相频等于-的时乃氏图上最左边的点越来越接近（-1，j0）,且越大，向左移动距离越大，故的存在与增大会降低系统的稳定性。5.校正环节的传递函数为，此校正环节的阶跃响应与单位脉冲函数分别如图1-13、1-14，伯德图与乃氏图分别如图1-15、

15、1-16所示；MATLAB语言如下：>> num=0.12,1num = 0.1200 1.0000>> den=0.01,0den = 0.0100 0>> sys=tf(num,den)sys = 0.12 s + 1 - 0.01 s Continuous-time transfer function.>> step(sys)>> impluse(sys)>> nyquist(sys)>> bode(sys)图1-13图1-14图1-15图1-16此积分环节：，由上图可知该积分环节：当输入为单位阶跃信号时

16、，输出为输入对时间的积分，输出随时间呈直线增长，当输入为脉冲信号时，输出为定值。该环节的输出在相位上滞后于输入,相角相差90°；在幅值上，经过该环节后，输出幅值明显增大，在伯德图上幅值相差40dB。此环节的稳定性比较好，可以做积分调节器来提高系统的稳定性。6.在校正环节的传递函数为时，整个系统的开环传递函数为：闭环传递函数为：此系统开环传递函数的乃氏图和伯德图分别如图1-17、1-18所示；MATLAB语言如下：>> num=15.9612,133.01num = 15.9612 133.0100>> den=0.00000,0.001,0.33545,0,0

17、den = 0 0.0010 0.3355 0 0>> den=0.000004,0.001,0.33545,0,0den = 0.0000 0.0010 0.3355 0 0>> sys=tf(num,den)sys = 15.96 s + 133 - 4e-06 s4 + 0.001 s3 + 0.3355 s2 Continuous-time transfer function.>> nyquist(sys)>> bode(sys)>> gm,pm,wcg,wcp=margin(sys)gm = 5.1236pm = 71.7979wcg = 285.9704wcp = 49.1595图1-17 图1-18其频域指标为（幅值裕量不是以dB为单位）：G1(s)=1时的频域指标为与G1(s)=1时的频域指标进行比较得到：其剪切频率、幅值裕量与相位裕量均变小，故系统快速性、稳定性均未得到改善。7.当，输入为单位阶跃函数时,SIMULINK仿真结果如下图所示：四实验总结 在完成该实验的过程中，我发现该实验运用了课本上学习的关于系统方块图的化