diff --git a/control/timeresp.py b/control/timeresp.py
index 946fc6270..f4f293bdf 100644
--- a/control/timeresp.py
+++ b/control/timeresp.py
@@ -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
@@ -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
pFad - Phonifier reborn
Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.
Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.
Alternative Proxies:
Alternative Proxy
pFad Proxy
pFad v3 Proxy
pFad v4 Proxy