典型环节的阶跃响应simulink仿真图.doc
文本预览下载声明
自控原理MATLAB实验
学校：*********
作者：李新华
图一 典型环节的阶跃响应simulink 仿真图
图二 比例环节
图三 积分环节
图四 比例积分环节
图五 比例微分环节
图六 比例微积分环节
图七 惯性环节
程序一：
wn=10;
kosi=[0,0.25,0.5,0.7,1,2]; %阻尼比分别为0,0.25,0.5,0.7,1,2
num=wn^2;
figure(1)
hold on
for i=1:6; %求阻尼比分别为0,0.25,0.5,0.7,1,2时的单位阶跃响应
den=[1,2*kosi(i)*wn,wn^2];
t=[0:0.01:4];
step(num,den,t)
end
hold off
title(step response) %标题为step response
图八 ωn一定ξ变化时系统单位阶跃响应曲线
程序二：
wn=[2,4,6,8,10,12]; %以2为最小值，12为最大值，步长为2
kosi=0.707;
figure(1)
hold on
for 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)
end
hold off
title(step response) %标题为step response
图九 ξ一定ωn变化时系统单位阶跃响应曲线
程序一：
num=[10];
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)
图四 原系统的根轨迹
程序五：
clear
num=[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=-128)(phase=-132))
disp([i mag1(i) phase(i) w(i)])
end
ii=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)*(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 -12
显示全部