seems like a first working frd implementation

repagh · repagh · commit 82ac4d2a1c8f · 2012-10-24T13:49:04.000Z
diff --git a/src/frdata.py b/src/frdata.py
@@ -172,7 +172,7 @@ def __init__(self, *args):
                 self.ifunc[i,j],u = splprep(
                     u=self.omega, x=[real(self.fresp[i, j, :]), 
                                      imag(self.fresp[i, j, :])], 
-                    w=1.0/(absolute(self.fresp[i, j, :])+0.001), s=0.01)
+                    w=1.0/(absolute(self.fresp[i, j, :])+0.001), s=0.0)
 
         Lti.__init__(self, self.fresp.shape[1], self.fresp.shape[0])
         
diff --git a/tests/frd_test.py b/tests/frd_test.py
@@ -100,7 +100,7 @@ def testAuto(self):
 
     def testBode(self):
         h1 = TransferFunction([1], [1, 2, 2])
-        omega = np.logspace(-1, 2, 30)
+        omega = np.logspace(-1, 2, 40)
         f1 = FRD(h1, omega)
         control.freqplot.nyquist(f1, np.logspace(-1, 2, 100))
         plt.figure(2)
