典型环节的阶跃响应simulink-仿真图

1、自控原理MATLAB实验 学校：* 作者：李新华 图一 典型环节的阶跃响应simulink 仿真图图二 比例环节图三 积分环节图四 比例积分环节图五 比例微分环节图六 比例微积分环节图七 惯性环节程序一：wn=10;kosi=0,0.25,0.5,0.7,1,2; %阻尼比分别为0,0.25,0.5,0.7,1,2num=wn2;figure(1)hold onfor i=1:6; %求阻尼比分别为0,0.25,0.5,0.7,1,2时的单位阶跃响应 den=1,2*kosi(i)*wn,wn2; t=0:0.01:4; step(num,den,t)endhold offtitle(step

2、 response) %标题为step response图八 n一定变化时系统单位阶跃响应曲线程序二：wn=2,4,6,8,10,12; %以2为最小值，12为最大值，步长为2kosi=0.707; figure(1)hold onfor i=1:6 num=wn(i)2; %分别求wn为2,4,6,8,10,12时的单位阶跃响应 den=1,2*kosi*wn(i),wn(i)2; t=0:0.01:4; step(num,den,t)endhold offtitle(step response) %标题为step response图九 一定n变化时系统单位阶跃响应曲线程序一：num=10;

3、den=0.02 0.3 1 0;nc,dc=cloop(num,den,-1);sys=tf(nc,dc);step(sys) 图一 原系统的单位阶跃响应曲线程序二：G=tf(10,0.02 0.3 1 0); %建立开环系统模型figure(1) %绘制伯德图bode(G); 图二 原系统的Bode图程序三：G=tf(10,0.02 0.3 1 0); %建立开环系统模型figure(1) nyquist(G); %绘制奈奎斯特图 图三 原系统的Nyquist曲线程序四：num=10;den=0.02 0.3 1 0;sys=tf(nc,dc);rltool(sys) 图四 原系统的根轨迹

4、程序五： clearnum=10;den=conv(1,0,conv(0.1,1,0.2,1);phm=45+5;phmd=-180+phm;w=logspace(-1,2,800)mag,phase=bode(num,den,w);mag1=20*log10(mag)for i=find(phase=-132) disp(i mag1(i) phase(i) w(i)endii=input(enter index for desired abd.) t=100/(w(ii)*(sqrt(1+0.01*w(ii)2)*(sqrt(1+0.04*w(ii)2)/w(ii)beta=(w(ii)*

5、(sqrt(1+0.01*w(ii)2)*(sqrt(1+0.04*w(ii)2)/10 %运行结果： 365.0000 11.5857 -128.0499 2.3265 366.0000 11.4933 -128.3496 2.3467 367.0000 11.4006 -128.6513 2.3671 368.0000 11.3076 -128.9548 2.3877 369.0000 11.2144 -129.2603 2.4084 370.0000 11.1209 -129.5677 2.4293 371.0000 11.0272 -129.8770 2.4504 372.0000 1

6、0.9332 -130.1882 2.4717 373.0000 10.8389 -130.5013 2.4931 374.0000 10.7443 -130.8164 2.5148 375.0000 10.6494 -131.1333 2.5366 376.0000 10.5543 -131.4522 2.5586 377.0000 10.4588 -131.7730 2.5809enter index for desired abd.373ii = 373t = 13.9700beta=0.2871程序六：G1=tf(0.2871*13.97 1,13.97,1); %建立开环系统模型，绘制校正装置的bode图figure(1)bode(G1);图五 校正装置的Bode图程序七：G1=tf(0.2871*13.97 1,13.97,1);num2=10;den2=conv(1,0,conv(0.1,1,0.2,1);G2=tf(num2,den2);G=G1*G2;figure(2) %绘制校正后的bode图margin(G); %并确定校正后系统的幅值裕度和相角裕度图六 校正后系统的Bode图程序八：num3=conv(10,0.2871*13.97 1);den3=conv(conv