Tb=[];Ts=[];
%调节时间为5倍的最长的调节时间；所以设t 为0~20*10*5=1000
t=0:0.01:1000;
%T2时间固定，T1比T2大，为T2的（1,20)倍之间；
T2=10;

for i = 1:length(T2)
    %T1比T2大，为T2的（1,20)倍之间；
    T1 =T2(i):0.1*T2(i):20*T2(i);
    for j = 1:length(T1)
        %调节时间的横轴
        Tb = [Tb T1(j)/T2(i)];
        %分子为：((1/T1)*(1/T2))
        num = [1/T1(j)*T2(i)];
        %分母为：s^2 +(1/T1+1/T2)s + ((1/T1)*(1/T2))
        den = [1  (1/T1(j)+1/T2(i)) (1/T1(j)*1/T2(i))];
        y = step(num, den, t);
        for k = length(y) : -1 :1
            if( abs(y(k) - 1 ) >= 0.05 )
                % 刚进入稳态，记入调节时间，
                Ts = [Ts (k*0.01)/T1(j)];
                break;
            end
            
        end
        
    end
    
end
plot(Tb ,Ts);grid on; xlim([1 20]);
xlabel('T1/T2'); ylabel('ts/T1');title('过阻尼二阶系统的调节时间特性');