%%%%   极点实部保持不变时，调节时间随阻尼比的变化规律

% 2%的调节时间；
Ts2 = [];
% 5%的调节时间；
Ts5 = [];
% 阻尼系数比
Xi = [];
% 欠阻尼极点实部 且固定为 1
re = 1;
t = 0 : 0.01 : 60;
% 欠阻尼极点虚部从10开始减小到0
for im = 10. :-0.02 :0
    Xi = [Xi cos(atan(im/re))];
    %角频率平方
    num = re * re + im * im;
    den = [ 1  2*re  re * re + im * im];
    T = 1/sqrt(num);
    y = step(num,den,t);
    for k = 6000 : -1 : 0
        % 时间是间隔 0.01 变化
        if(abs(y(k)-1)) >= 0.05, Ts2 = [ Ts2, k*0.01 ];
            break;
        end    
    end
    
    for k = 6000 : -1 : 0
        % 时间是间隔 0.01 变化
        if(abs(y(k)-1)) >= 0.02, Ts5 = [ Ts5, k*0.01 ];
            break;
        end    
    end
end
    
    plot(Xi,Ts2,'b-',Xi,Ts5,'r-');
    xlabel('Xi'),ylabel('Ts');
    title('极点实部保持不变时，调节时间随阻尼比的变化规律');grid on;