git.gag.com Git - debian/gnuradio/blob - gnuradio-core/src/python/gnuradio/gr/qa_goertzel.py

   1 #!/usr/bin/env python
   2 #
   3 # Copyright 2006 Free Software Foundation, Inc.
   4 #
   5 # This file is part of GNU Radio
   6 #
   7 # GNU Radio is free software; you can redistribute it and/or modify
   8 # it under the terms of the GNU General Public License as published by
   9 # the Free Software Foundation; either version 2, or (at your option)
  10 # any later version.
  11 #
  12 # GNU Radio is distributed in the hope that it will be useful,
  13 # but WITHOUT ANY WARRANTY; without even the implied warranty of
  14 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  15 # GNU General Public License for more details.
  16 #
  17 # You should have received a copy of the GNU General Public License
  18 # along with GNU Radio; see the file COPYING.  If not, write to
  19 # the Free Software Foundation, Inc., 51 Franklin Street,
  20 # Boston, MA 02110-1301, USA.
  21 #
  22
  23 from gnuradio import gr, gr_unittest
  24 from math import pi, cos
  25
  26 class test_goertzel(gr_unittest.TestCase):
  27
  28     def setUp(self):
  29         self.fg = gr.flow_graph()
  30
  31     def tearDown(self):
  32         self.fg = None
  33
  34     def make_tone_data(self, rate, freq):
  35         return [cos(2*pi*x*freq/rate) for x in range(rate)]
  36
  37     def transform(self, src_data, rate, freq):
  38         src = gr.vector_source_f(src_data, False)
  39         dft = gr.goertzel_fc(rate, rate, freq)
  40         dst = gr.vector_sink_c()
  41         self.fg.connect(src, dft, dst)
  42         self.fg.run()
  43         return dst.data()
  44
  45     def test_001(self): # Measure single tone magnitude
  46         rate = 8000
  47         freq = 100
  48         bin = freq
  49         src_data = self.make_tone_data(rate, freq)
  50         expected_result = 0.5
  51         actual_result = abs(self.transform(src_data, rate, bin)[0])
  52         self.assertAlmostEqual(expected_result, actual_result, places=5)
  53
  54     def test_002(self): # Measure off frequency magnitude
  55         rate = 8000
  56         freq = 100
  57         bin = freq/2
  58         src_data = self.make_tone_data(rate, freq)
  59         expected_result = 0.0
  60         actual_result = abs(self.transform(src_data, rate, bin)[0])
  61         self.assertAlmostEqual(expected_result, actual_result, places=5)
  62
  63 if __name__ == '__main__':
  64     gr_unittest.main()