steppoly.m

来自「这是一个用MATLAB工具设计自动控制理论中系统的仿真」· M 代码 · 共 25 行

num1=input('"num vector .eg:[0 0 1]:"');
den1=input('"den vector .eg:[1 3 2]:"');
G1=tf(num1,den1);G2=1;
G=feedback(G1,G2)
step(G);grid on;
[y,t]=step(G);
[Mp,ind]=max(y);
Yss=dcgain(G);
Tp=t(ind);
Pos=100*(Mp-Yss)/Yss;
Dimt=length(t);
    for i=1:Dimt;
    if  y(i)>1.02*Yss;
        Ts=t(i);
    elseif y(i)<0.98*Yss
        Ts=t(i);
    end
    end
disp('[Yss,Pos(%),Ts(S)]');[Yss,Pos,Ts]
disp('Get sel-define T-axit,input"steprun"');
disp('Compare in one graph,input "hold on"');