Dec 30, 2020 · Dec 29, 2020
diff --git a/control/timeresp.py b/control/timeresp.py
        The system whose time response is to be computed
    is_step : bool
        Scales the dc value by the magnitude of the nonzero mode since
        integrating the impulse response gives
        :math:`\int e^{-\lambda t} = -e^{-\lambda t}/ \lambda`
        integrating the impulse response gives
        :math:`\\int e^{-\\lambda t} = -e^{-\\lambda t}/\\lambda`
        Default is True.

    Returns
        p_u, p = p[m_u], p[~m_u]
        if p_u.size > 0:
            m_u = (p_u.real < 0) & (np.abs(p_u.imag) < sqrt_eps)
            t_emp = np.max(log_decay_percent / np.abs(np.log(p_u[~m_u])/dt))
            tfinal = max(tfinal, t_emp)
            if np.any(~m_u):
                t_emp = np.max(
                    log_decay_percent / np.abs(np.log(p_u[~m_u]) / dt))
                tfinal = max(tfinal, t_emp)

        # zero - negligible effect on tfinal
        m_z = np.abs(p) < sqrt_eps
Original file line number	Diff line number	Diff line change
Expand Up		@@ -846,8 +846,8 @@ def _ideal_tfinal_and_dt(sys, is_step=True):
		The system whose time response is to be computed
		is_step : bool
		Scales the dc value by the magnitude of the nonzero mode since
		integrating the impulse response gives
		:math:`\int e^{-\lambda t} = -e^{-\lambda t}/ \lambda`
		integrating the impulse response gives
		:math:`\\int e^{-\\lambda t} = -e^{-\\lambda t}/\\lambda`
		Default is True.

		Returns
Expand DownExpand Up		@@ -900,8 +900,10 @@ def _ideal_tfinal_and_dt(sys, is_step=True):
		p_u, p = p[m_u], p[~m_u]
		if p_u.size > 0:
		m_u = (p_u.real < 0) & (np.abs(p_u.imag) < sqrt_eps)
		t_emp = np.max(log_decay_percent / np.abs(np.log(p_u[~m_u])/dt))
		tfinal = max(tfinal, t_emp)
		if np.any(~m_u):
		t_emp = np.max(
		log_decay_percent / np.abs(np.log(p_u[~m_u]) / dt))
		tfinal = max(tfinal, t_emp)

		# zero - negligible effect on tfinal
		m_z = np.abs(p) < sqrt_eps
Expand Down