The following no longer causes a Numpy ValueError:

goretkin · goretkin · commit 392f729e6ac6 · 2011-04-17T18:45:29.000Z
(cases where the order of the numerator is higher than the order of the denominator)

s = control.tf.TransferFunction([1, 0], [1])
control.bode(s)
diff --git a/src/xferfcn.py b/src/xferfcn.py
@@ -643,8 +643,10 @@ def _common_den(self):
         # are the same size as the denominator.
         for i in range(self.outputs):
             for j in range(self.inputs):
-                num[i][j] = insert(num[i][j], zeros(len(den) - len(num[i][j])),
-                    zeros(len(den) - len(num[i][j])))
+                pad = len(den) - len(num[i][j]) 
+                if(pad>0):
+                    num[i][j] = insert(num[i][j], zeros(pad),
+                        zeros(pad))
         # Finally, convert the numerator to a 3-D array.
         num = array(num)
         # Remove trivial imaginary parts.  Check for nontrivial imaginary parts.
