自动控制原理实验
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
实验一 典型环节的MATLAB 仿真
实验内容
按下列各典型环节的传递函数,建立相应的SIMULINK 仿真模型,观察并记录其单位阶跃响应波形。
① 比例环节1)(1=s G 和2)(1=s G ;
② 惯性环节11)(1+=
s s G 和1
5.01)(2+=s s G
③ 积分环节s
s G 1)(1
④ 微分环节s s G )(1
⑤ 比例+微分环节(PD )2)(1+=s s G 和1)(2+=s s G
⑥ 比例+积分环节(PI )s s G 11)(1+=和s
s G 211)(2+=
实验总结
通过这次接触MATLAB,真正的体会到了它强大的数值计算和符号计算功能,以及强大的数据可视化、人际智能交互能力。该工具主要处理以传递函数为主要特征的经典控制和以状态空间为主要特征的现代控制中的主要问题,它能够使图形生动形象的展现给我们,使理解更深刻。
实验二线性定常系统的瞬态响应
实验内容
1.观察函数step( )和impulse( )的调用格式,假设系统的传递函数模型为
1
4647
3)(2
342++++++=s s s s s s s G 可以用几种方法绘制出系统的阶跃响应曲线?试分别绘制。
1.阶跃响应
num=[0 0 1 3 7 ]; den=[1 4 6 4 1]; step(num,den) grid
xlabel(‘t/s ’),ylabel(‘c(t)’)
title(‘Uint-step Respinse of G(s)=(s^2+3s+7)/(s^4+4s^3+6s^2+4s+1)’)
脉冲响应1
num=[0 0 1 3 7]; den=[1 4 6 4 1]; impulse(num,den) grid
title(‘Unit-impulse Response of G(s)= (s^2+3s+7)/(s^4+4s^3+6s^2+4s+1)’)
脉冲响应2 num=[0 0 1 3 7]; den=[1 4 6 4 1]; step(num,den) grid
title(‘Unit-impulse Response of G(s)= (s^2+3s+7)/(s^4+4s^3+6s^2+4s+1)’)
斜坡响应
num=[0 0 1 3 7];
den=[1 4 6 4 1];
step(num,den)
grid
title(‘Unit-impulse Response of G(s)= (s^2+3s+7)/(s^4+4s^3+6s^2+4s+1)’)
2.对典型二阶系统
2
2
2
2)(n
n n s s s G ωζωω++= 1)分别绘出)/(2s rad n =ω,ζ分别取0,0.25,0.5,1.0和2.0时的单位阶跃响应曲线,分析参数ζ对系统的影响,并计算ζ=0.25时的时域性能指标ss s p r p e t t t ,,,,σ。
num=[ 0 0 4 ]; den1=[1 0 4 ]; den2=[1 1 4 ]; den3=[1 2 4 ]; den4=[1 4 4 ]; den5=[1 8 4]; t=0:0.1:10;
step(num,den1,t) grid
text( 4,1.7,’0’) hold step(num,den2,t) text(3.3,1.5,’0.5’) step(num,den3,t) text(3.5,1.2,’0.5 ‘ ) step(num,den4,t) text( 3.3,0.9,’1.0’ ) step(num,den5,t) text(3.3,0.6,’2.0’)
title(‘Step-Response Curves for G(s)=1/[s^2+2(zeta)s+1]’)
2)绘制出当ζ=0.25, n ω分别取1,2,4,6时单位阶跃响应曲线,分析参数n ω对系统的影响。
num1=[0 0 1]; den1=[1 0.5 1]; t=0:0.1:10;
step(num1,den1,t); grid; hold on text(3.1,1.4,'wn=1')
num2=[0 0 4]; den2=[1 1 4]; step(num2,den2,t); hold on text(1.7,1.4,'wn=2')
num3=[0 0 16]; den3=[1 2 16]; step(num3,den3,t);hold on text(0.5,1.4,'wn=3')
num4=[0 0 36];den4=[1 3 36]; step(num4,den4,t);hold on text(3.3,1.4,'wn=4')
阻尼比
Wn 超调量 上升时间 峰值时间 调节时间 0.25
1
0.44 1.9 3.2 10.2 2 0.44 0.94 1.6 5.29 4 0.44 0.95 1.6 5.4 6
0.44
0.48
0.8
2.7
实验总结
通过本次试验我熟练掌握step( )函数和impulse( )函数的使用方法,进一步加深了对线性系统在单位阶跃、单位脉冲及单位斜坡函数作用下的响应的变化的印象,并且学会了通过响应曲线观测特征参量ζ和n ω对二阶系统性能的影响。
实验三 线性系统稳态误差研究
三、实验内容
1.系统的特征方程式为010532234=++++s s s s ,试用三种判稳方式判别该系统的稳定性。 直接求根判稳
roots([2,1,3,5,10])
所得结果有实部为正数,故系统不稳定。 实验总结
通过本实验我学会了如何高阶系统的稳定性和验证稳定判据的正确性的方法,并且了解系统增益变化对系统稳定性的影响和系统结构和稳态误差之间的关系。知道了可以利用MATLAB 中的tf2zp 函数求出系统的零极点,或者利用root 函数求分母多项式的根来确定系统的闭环极点,从而判断系统的稳定性。
实验五 线性系统的频域分析
实验内容
1.典型二阶系统
2
2
22)(n
n n
s s s G ωζωω++=