Merge 5c37bdf into 854f5f7

aWZHY0yQH81uOYvH · web-flow · commit ae9319550556 · 2026-03-18T12:30:57.000-07:00
diff --git a/control/lti.py b/control/lti.py
@@ -211,8 +211,14 @@ def bandwidth(self, dbdrop=-3):
             # solve for the bandwidth, use scipy.optimize.root_scalar() to
             # solve using bisection
             import scipy
+
+            if self.isdtime(strict=True):
+                cvt_w = lambda w: np.exp(1j * w * self.dt)
+            else:
+                cvt_w = lambda w: 1j * w
+
             result = scipy.optimize.root_scalar(
-                lambda w: np.abs(self(w*1j)) - np.abs(dcgain)*10**(dbdrop/20),
+                lambda w: np.abs(self(cvt_w(w))) - np.abs(dcgain)*10**(dbdrop/20),
                 bracket=[omega[idx_dropped[0] - 1], omega[idx_dropped[0]]],
                 method='bisect')
 
diff --git a/control/tests/lti_test.py b/control/tests/lti_test.py
@@ -109,6 +109,11 @@ def test_bandwidth(self):
         np.testing.assert_allclose(sys1.bandwidth(), 0.099762834511098)
         np.testing.assert_allclose(bandwidth(sys1), 0.099762834511098)
 
+        # test a first-order discrete-time system, compared with matlab
+        sysd1 = tf([0.1, 0], [1, -0.9], 1)
+        np.testing.assert_allclose(sysd1.bandwidth(), 0.105207775532932)
+        np.testing.assert_allclose(bandwidth(sysd1), 0.105207775532932)
+
         # test a second-order system, compared with matlab
         wn2 = 1
         zeta2 = 0.001
