资源描述：
《控制工程基础matlab仿真分析》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、机械工程控制基础实验报告专业：班级：学号：姓名：Stepwn=1xi=[0.1,0.3,0.5,0.707,1,3,5,7]sys1=tf([wn^2],[1,2*wn*xi(1),wn^2]);step(sys1)holdonfori=2:length(xi)sys1=tf([wn^2],[1,2*wn*xi(i),wn^2]);step(sys1)endholdoffprintxklfigstep-dbitmap分析：t趋向于无穷时，二阶系统单位阶跃响应稳态值为1；随着xi的增加，最大超调量减小，响应时间变慢，振荡特性减弱；xi>=1时，最大超调量为0，过渡过程具有单调上
2、升的特性.Impulsewn=1xi=[0.1,0.3,0.5,0.707,1,3,5,7]sys1=tf([wn^2],[1,2*wn*xi(1),wn^2]);impulse(sys1)holdonfori=2:length(xi)sys1=tf([wn^2],[1,2*wn*xi(i),wn^2]);impulse(sys1)endholdoffprintxklfigimpulse-dbitmap分析：t趋向于无穷时，二阶系统单位脉冲响应稳态值为0；随着xi的增加，最大超调量减小，衰减愈快；xi>=1时，无振荡特性.Bodewn=1xi=[0.1,0.3,0.5,0.7
3、07,1,3,5,7]sys1=tf([wn^2],[1,2*wn*xi(1),wn^2]);bode(sys1)holdonfori=2:length(xi)sys1=tf([wn^2],[1,2*wn*xi(i),wn^2]);bode(sys1)endholdoffprintxklfigbode-dbitmap分析：振荡环节的对数相频特性对称于点（1，－90°）；xi越小，谐振峰值越大；xi趋于0时，最大超调量趋于无穷；w趋于无穷时，幅值变化斜率为固定值。Nyquistwn=1xi=[0.1,0.3,0.5,0.707,1,3,5,7]sys1=tf([wn^2],[1
4、,2*wn*xi(1),wn^2]);nyquist(sys1)holdonfori=2:8sys1=tf([wn^2],[1,2*wn*xi(i),wn^2]);nyquist(sys1)endholdoffprintxklfignyquist–dbitmap分析：振荡环节nyquist图始于点（1，j0），终于点（0，j0）,是关于re轴对称的封闭曲线；xi越大时nyquist曲线包围的面积越小。Magwn=1xi=[0.1,0.3,0.5,0.707,1,3,5,7]sys1=tf([wn^2],[1,2*wn*xi(1),wn^2])w1=0:0.01:10;[mag
5、1,phase1]=bode(sys1,w1);plot(w1(:),mag1(1,:))holdonfori=2:length(xi)sys1=tf([wn^2],[1,2*wn*xi(i),wn^2])w1=0:0.01:10;[mag1,phase1]=bode(sys1,w1);plot(w1(:),mag1(1,:))endholdoffxlabel('频率')ylabel('幅值')printxklfigmag-dbitmap分析：振荡环节随着xi的增加，幅值减小；xi>=0.707时，曲线单调下降，谐振峰值不存在。Phasewn=1xi=[0.1,0.3,0.5
6、,0.707,1,3,5,7]sys1=tf([wn^2],[1,2*wn*xi(1),wn^2])w1=0:0.01:10;[mag1,phase1]=bode(sys1,w1);plot(w1(:),phase1(1,:))holdonfori=2:8sys1=tf([wn^2],[1,2*wn*xi(i),wn^2])w1=0:0.01:10;[mag1,phase1]=bode(sys1,w1);plot(w1(:),phase1(1,:))endholdoffxlabel('频率')ylabel('相位')printxklfigphase-dbitmap分析：相频特
7、性对称于点（1，－90°）；随着频率增加，滞后相位角增大；xi越小，相位变化率越快。Realwn=1xi=[0.1,0.3,0.5,0.707,1,3,5,7]sys1=tf([wn^2],[1,2*wn*xi(1),wn^2])w1=0:0.01:10;[mag1,phase1]=bode(sys1,w1);plot(w1(:),mag1(1,:).*cos(phase1(1,:)*pi/180))holdonfori=2:length(xi)sys1=tf([wn^2],[1,2*wn*xi(i),
