控制系统计算机仿真(matlab)实验四实验报告

实验四控制系统计算机辅助分析一、实验目的1、 掌握如何使用Matlab进行系统的时域分析2、 掌握如何使用Matlab进行系统的频域分析3、 掌握如何使用Matlab进行系统的根轨迹分析二、 实验学时：2学时三、 试验原理：1、稳定性的基本概念及必要条件根据李雅普诺夫稳定性理论，线性控制系统的稳定性可定义如下：如果线性控制系统在初始扰动的影响下，其动态过程随时间的推移逐渐衰减并趋于零(原平衡工作点)，则称系统渐近稳定，简称稳定。否则，若在初始扰动影响下，系统的动态过程随时间的推移而发散，则称系统不稳定。线性系统稳定的充分必要条件是：闭环系统特征方程的所有根均具有负实部；或者说，闭环传递函数的极点均严格位于左半s平面。由上述线性系统稳定性概念及系统稳定的充分必要条件可知，判定线性系统稳定性的最直接方法就是求出闭环系统特征方程的所有根或者全部闭环极点，根据特征方程所有根是否具有负实部或闭环极点是否全部位于左半s平面来判定系统的稳定性。四、实验内容：(三题选做两题)1、时域分析(1)根据下面传递函数模型：绘制其单位阶跃响应曲线并从图上读取最大超调量，并求出单位脉冲响应曲线。=严于+)程序：s=tf('s');G=5*(s"2+5*s+6)/(s"3+6*s"2+10*s+8);step(G);grid;holdon;impulse(G);结果：超调量二（|3.75-4|）/4x100%=6.25%StepResponse(2)典型二阶系统传递函数为：%)= ——rs'-肿+cyj当Z=0.7,3n取2、4、6、8、10、12的单位阶跃响应。程序：kesi=0.7forwn=2:2:12num二wn2den二[12*kesi*wnwn"2];G=tf(num,den);t=0:0.01:10;step(G);holdon;endtitle('wn不同值下的单位阶跃响应');xlabel('t');ylabel('阶跃响应');grid;结果：1.4wn不同值下的单位阶跃响应0.40.200.511.522.533.542186100应响跃阶03）典型二阶系统传递函数为：t(sec)当3n=6，Z取0.2、0.4、0.6、0.8、1.0、1.5、2.0的单位阶跃响应。程序：wn=6;forkesi=0.2:0.2:1.0num二wn2den二[12*kesi*wnwn"2];G=tf(num，den);t=0:0.01:10;step(G);holdon;endforkesi=1.5:0.5:2.0num二wn2den二[12*kesi*wnwn"2];G=tf(num，den);t=0:0.01:10;step(G);holdon;endtitle('O不同值下的单位阶跃响应');xlabel('t');ylabel('阶跃响应');grid;结果：＜不同值下的单位阶跃响应＜不同值下的单位阶跃响应2、频域分析（1）典型二阶系统传递函数为：当z2、频域分析（1）典型二阶系统传递函数为：当z=0.7,3n取2、4、6、8、10、程序：12的伯德图kesi=0.7forwn=2:2:12num二wn2den二[12*kesi*wnwn"2];G=tf(num,den);t=0:0.01:10;bode(G);holdon;endtitle('wn不同值下的伯德图’);grid;结果wn不同值下的伯德图3Frequency(rad/sec)wn不同值下的伯德图3Frequency(rad/sec)(2)典型二阶系统传递函数为：当3n=6，Z取0.2、0.4、0.6、0.8、1.0、1.5、2.0的伯德图。程序：wn=6;forkesi=0.2:0.2:1.0num二wn2den二[12*kesi*wnwn"2];G=tf(num，den);t=0:0.01:10;bode(G);holdon;endforkesi=1.5:0.5:2.0num二wn2den二[12*kesi*wnwn"2];G=tf(num，den);t=0:0.01:10;bode(G);holdon;endtitle('6不同值下的伯德图');grid;运行结果：<5不同值下的伯德图1G(s)= (3)已知二阶系统传递函数Ks2+22s+1绘制阻尼系数E分别为1.2,1.0,0.5和0.25时系统的Nyquist曲线。程序：forkesi=0.25:0.25:0.51num=1;den=[12*kesi1];G=tf(num,den);nyquist(G);holdon;endforkesi=1.0:0.2:1.2num=1;den=[12*kesi1];G=tf(num,den);nyquist(G);holdon;endtitle('Nyquist曲线');运行结果：Nyquist曲线515.-.515.-.0-1--251510sixA-1-0.500.511.5RealAxis3、根轨迹分析根据下面负反馈系统的开环传递函数，绘制系统根轨迹，并分析系统稳定的K值范围。程序：num=[1];dun=[1320];rlocus(num,dun);[k,p]=rlocfind(num,dun)运行结果：k=6.3711p=-3.03310.0166+1.4492i0.0166-1.4492iRootLocus-1-2-3-4-5-7-6-5stem:sysGain:0.385Po