python/test/tf.py

   1 #!/usr/bin/env python
   2 import os, sys
   3
   4 base = os.path.dirname(__file__)
   5 path = os.path.dirname(base)
   6 sys.path.insert(0, path)
   7 if len(sys.argv) > 1:
   8     sys.stdout = open(sys.argv[1], "w")
   9
  10 import liba  # type: ignore
  11
  12 num = [6.59492796e-05, 6.54019884e-05]
  13 den = [1.0, -1.97530991, 0.97530991]
  14
  15 tf = liba.tf(num, den[1:])
  16
  17 N = 1000
  18 T = 0.5 / N
  19 x = [0.0] * N
  20 for i in range(N):
  21     x[i] = i * T
  22
  23 y = [0.0] * N
  24 for i in range(N):
  25     y[i] = tf(x[i])
  26
  27 for i in range(N):
  28     print("%g,%g" % (x[i], y[i]))
  29
  30 tf.num = tf.num
  31 tf.den = tf.den
  32 tf.input
  33 tf.output