clc;
clear;
num = [3];
den = [1 2 10];%把分母进行了一个化简
step(num,den)%然后拖动鼠标点去查看就行了
%ltiview 打开LTI Viewer
G = tf(num,den);
ltiview
%编程获得\
t=0:0.001:10;
y = step(num,den,t);
n = length(t);
[ymax,ind] = max(y);
yss = y(n);%稳态值
mp = (ymax-yss)/yss;%超调量
for k=1:n
if y(k)<=yss & y(k+1) >=yss
tr = t(k+1);
break
end
end
for i=n:-1:1%该层循环用于计算调整时间,我们是从最终时间(无穷远的时间)往前开始遍历查找,否则从前开始查找的话会出错
if y(i)>=1.02*yss|y(i)<=0.98*yss%当得到的值与稳态值差的绝对值与稳态值之比为2%时说明到达了调整时间
ts=t(i);%ts存储调整时间
break
end
end
disp('峰值时间为:')
disp(t(ind))
disp('上升时间为:')
disp(tr)
disp('调整时间为:')
disp(ts)
disp('超调量为:')
disp(mp)
- 1
- 2
- 3
前往页