资源描述:
《matlab语言与控制系统仿真-参考答案-第5章》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、MATLAB语言与控制系统仿真实训教程-参考答案-第5章5.6控制系统的时域响应MATLAB仿真实训5.6.1实训目的1.学会利用MATLAB绘制系统的单位阶跃响应曲线,掌握读取系统动态性能指标的方法;2.学会利用MATLAB绘制系统的单位脉冲响应曲线的方法;3.掌握利用MATLAB绘制系统的零输入响应曲线的方法;4.掌握利用MATLAB绘制系统的一般输入响应曲线的方法;5.学会通过仿真曲线读取相关信息,并依据有关信息进行系统的时域分析。5.6.2实训内容1.编写程序求取下列各系统的单位阶跃响应,完成表5-5并记录相关曲线。解:>>n1=16;>>d1=[1,3.2,16];>>sys1=t
2、f(n1,d1);>>step(sys1)>>n2=16;>>d2=[1,2.4,16];>>sys2=tf(n2,d2);>>step(sys2)18MATLAB语言与控制系统仿真实训教程-参考答案-第5章>>n3=16;>>d3=[1,1.6,16];>>sys3=tf(n3,d3);>>step(sys3)>>n4=16;>>d4=[1,1,16];>>sys4=tf(n4,d4);>>step(sys4)18MATLAB语言与控制系统仿真实训教程-参考答案-第5章表5-5序号()计算值实验计算值实验计算值实验值10.441.25381.250.85690.8632.18752.120
3、.341.37231.370.82330.8282.91672.8130.241.52661.530.80160.84.37504.940.12541.67311.670.79160.8037.00007.33>>z1=0.4;w=4;cmax1=1+exp(-z1*pi/sqrt(1-z1^2));tp1=pi/(w*sqrt(1-z1^2));ts1=3.5/(z1*w);[cmax1,tp1,ts1]ans=1.25380.85692.1875>>z2=0.3;w=4;cmax2=1+exp(-z2*pi/sqrt(1-z2^2));tp2=pi/(w*sqrt(1-z2^2));ts
4、2=3.5/(z2*w);[cmax2,tp2,ts2]18MATLAB语言与控制系统仿真实训教程-参考答案-第5章ans=1.37230.82332.9167>>z3=0.2;w=4;cmax3=1+exp(-z3*pi/sqrt(1-z3^2));tp3=pi/(w*sqrt(1-z3^2));ts3=3.5/(z3*w);[cmax3,tp3,ts3]ans=1.52660.80164.3750>>z4=0.125;w=4;cmax4=1+exp(-z4*pi/sqrt(1-z4^2));tp4=pi/(w*sqrt(1-z4^2));ts4=3.5/(z4*w);[cmax4,tp4
5、,ts4]ans=1.67310.79167.0000说明:对于二阶欠阻尼系统(),若系统的闭环传递函数为则系统单位阶跃响应的输出最大值峰值时间调整时间估算值(以5%为误差带)(以2%为误差带)2.已知二阶系统的闭环传递函数如下,编程求取系统的单位阶跃响应并完成表5-6,记录相关曲线。(1)(2)(3)(4)18MATLAB语言与控制系统仿真实训教程-参考答案-第5章>>n1=8;>>d1=[1,2,8];>>sys1=tf(n1,d1);>>impulse(sys1)>>n2=10;>>d2=[1,2,10];>>sys2=tf(n2,d2);>>impulse(sys2)>>n3=12;
6、>>d3=[1,2,12];>>sys3=tf(n3,d3);>>impulse(sys3)18MATLAB语言与控制系统仿真实训教程-参考答案-第5章>>n4=16;>>d4=[1,2,16];>>sys4=tf(n4,d4);>>impulse(sys4)表5-6序号实验实验10.3542.8281.790.44220.3163.1622.080.44230.2893.4642.360.38740.2542.840.3513.已知某单位负反馈系统的开环传递函数为若系统的输入信号分别为18MATLAB语言与控制系统仿真实训教程-参考答案-第5章(1),(2),(3),(4),编写程序分别求
7、取系统的在给定的输入信号下的响应,记录相应的曲线。>>n=5;>>d=conv([1,1,0],[0.1,1]);>>[nc,dc]=cloop(n,d)>>sysc=tf(nc,dc);>>t=0:0.01:20;>>u1=t;>>u2=t.*exp(-0.5*t);>>u3=sin(2*t).*exp(-0.5*t);>>u4=cos(2*t).*exp(-3*t);>>subplot(221);lsim