jisuanjieyuexiangying.m

阶跃相应绘图

g1=tf(10,[1 0 -1]);
g2=tf([0.3 1],1);
g=feedback(g1,g2);
step(g,6)
[y,t]=step(g,6);
yfinal=y(length(t));
[ytop,i]=max(y);
sigma=(ytop-yfinal)/yfinal*100
tp=t(i)
y50=0.5*yfinal;
y10=0.1*yfinal;
y90=0.9*yfinal;
i=1;
while y(i)<y50
    i=i+1;
end
td=t(i)
i=1;
while y(i)<y10
    i=i+1;
end
t1=t(i)
i=1;
while y(i)<y90
    i=i+1;
end
t2=t(i);
tr=t2-t1
i=length(t);
while y(i)<yfinal*1.02&&i>=1
    i=i-1;
end
t3=t(i);
i=length(t);
while y(i)>yfinal*0.98&&i>=1
    i=i-1;
end
t4=t(i);
if t3==0
    ts=t4
else
    ts=max(t3,t4)
end