python/test/hpf.py

   1 #!/usr/bin/env python
   2 import os, sys, math
   3
   4 sys.path.insert(0, os.getcwd())
   5 if len(sys.argv) > 1:
   6     sys.stdout = open(sys.argv[1], "w")
   7
   8 import liba  # type: ignore
   9
  10 Ts = 1.0
  11 f_0 = 1.0
  12 f_e = 1000.0
  13 c = (f_e - f_0) / Ts
  14 T = 1.0 / 10000
  15 t = 0
  16 x = []
  17 while t <= Ts:
  18     x.append(t)
  19     t += T
  20 hpf = liba.hpf(10, 0.01)
  21 for i in range(len(x)):
  22     x1 = math.cos(2 * math.pi * (f_0 * x[i] + 0.5 * c * x[i] * x[i]))
  23     print("%g,%g,%g" % (x[i], x1, hpf(x1)))
  24
  25 hpf.zero()
  26 hpf.alpha
  27 hpf.output
  28 hpf.input