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

   1 #!/usr/bin/env python
   2 #
   3 # Copyright 2004,2007 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 3, 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
  25 class test_iir (gr_unittest.TestCase):
  26
  27     def setUp (self):
  28         self.tb = gr.top_block ()
  29
  30     def tearDown (self):
  31         self.tb = None
  32
  33     def test_iir_direct_001 (self):
  34         src_data = (1, 2, 3, 4, 5, 6, 7, 8)
  35         fftaps = ()
  36         fbtaps = ()
  37         expected_result = (0, 0, 0, 0, 0, 0, 0, 0)
  38         src = gr.vector_source_f (src_data)
  39         op = gr.iir_filter_ffd (fftaps, fbtaps)
  40         dst = gr.vector_sink_f ()
  41         self.tb.connect (src, op)
  42         self.tb.connect (op, dst)
  43         self.tb.run ()
  44         result_data = dst.data ()
  45         self.assertFloatTuplesAlmostEqual (expected_result, result_data)
  46
  47     def test_iir_direct_002 (self):
  48         src_data = (1, 2, 3, 4, 5, 6, 7, 8)
  49         fftaps = (2,)
  50         fbtaps = (0,)
  51         expected_result = (2, 4, 6, 8, 10, 12, 14, 16)
  52         src = gr.vector_source_f (src_data)
  53         op = gr.iir_filter_ffd (fftaps, fbtaps)
  54         dst = gr.vector_sink_f ()
  55         self.tb.connect (src, op)
  56         self.tb.connect (op, dst)
  57         self.tb.run ()
  58         result_data = dst.data ()
  59         self.assertFloatTuplesAlmostEqual (expected_result, result_data)
  60
  61     def test_iir_direct_003 (self):
  62         src_data = (1, 2, 3, 4, 5, 6, 7, 8)
  63         fftaps = (2, 11)
  64         fbtaps = (0,  0)
  65         expected_result = (2, 15, 28, 41, 54, 67, 80, 93)
  66         src = gr.vector_source_f (src_data)
  67         op = gr.iir_filter_ffd (fftaps, fbtaps)
  68         dst = gr.vector_sink_f ()
  69         self.tb.connect (src, op)
  70         self.tb.connect (op, dst)
  71         self.tb.run ()
  72         result_data = dst.data ()
  73         self.assertFloatTuplesAlmostEqual (expected_result, result_data)
  74
  75     def test_iir_direct_004 (self):
  76         src_data = (1, 2, 3, 4, 5, 6, 7, 8)
  77         fftaps = (2, 11)
  78         fbtaps = (0,  -1)
  79         expected_result = (2, 13, 15, 26, 28, 39, 41, 52)
  80         src = gr.vector_source_f (src_data)
  81         op = gr.iir_filter_ffd (fftaps, fbtaps)
  82         dst = gr.vector_sink_f ()
  83         self.tb.connect (src, op)
  84         self.tb.connect (op, dst)
  85         self.tb.run ()
  86         result_data = dst.data ()
  87         self.assertFloatTuplesAlmostEqual (expected_result, result_data)
  88
  89     def test_iir_direct_005 (self):
  90         src_data = (1, 2, 3, 4, 5, 6, 7, 8)
  91         fftaps = (2, 11,  0)
  92         fbtaps = (0,  -1, 3)
  93         expected_result = (2, 13, 21, 59, 58, 186, 68, 583)
  94         src = gr.vector_source_f (src_data)
  95         op = gr.iir_filter_ffd (fftaps, fbtaps)
  96         dst = gr.vector_sink_f ()
  97         self.tb.connect (src, op)
  98         self.tb.connect (op, dst)
  99         self.tb.run ()
 100         result_data = dst.data ()
 101         self.assertFloatTuplesAlmostEqual (expected_result, result_data)
 102
 103     def test_iir_direct_006 (self):
 104         src_data = (1, 2, 3, 4, 5, 6, 7, 8)
 105         expected_result = (2, 13, 21, 59, 58, 186, 68, 583)
 106         fftaps = (2, 1)
 107         fbtaps = (0, -1)
 108         src = gr.vector_source_f (src_data)
 109         op = gr.iir_filter_ffd (fftaps, fbtaps)
 110         fftaps = (2, 11,  0)
 111         fbtaps = (0,  -1, 3)
 112         op.set_taps (fftaps, fbtaps)
 113         dst = gr.vector_sink_f ()
 114         self.tb.connect (src, op)
 115         self.tb.connect (op, dst)
 116         self.tb.run ()
 117         result_data = dst.data ()
 118         self.assertFloatTuplesAlmostEqual (expected_result, result_data)
 119
 120     def test_iir_direct_007 (self):
 121         src_data = (1, 2, 3, 4, 5, 6, 7, 8)
 122         expected_result = (2,2,5,5,8,8,11,11)
 123         fftaps = (2, 1)
 124         fbtaps = (0, -1)
 125         src = gr.vector_source_f (src_data)
 126         op = gr.iir_filter_ffd (fftaps, fbtaps)
 127         fftaps = (2,0,1)
 128         fbtaps = (0, -1)
 129         op.set_taps (fftaps, fbtaps)
 130         dst = gr.vector_sink_f ()
 131         self.tb.connect (src, op)
 132         self.tb.connect (op, dst)
 133         self.tb.run ()
 134         result_data = dst.data ()
 135         self.assertFloatTuplesAlmostEqual (expected_result, result_data)
 136
 137     def test_iir_direct_008 (self):
 138         src_data = (1, 2, 3, 4, 5, 6, 7, 8)
 139         expected_result = (2,4,4,10,18,14,26,56)
 140         fftaps = (2,)
 141         fbtaps = (0, 1)
 142         src = gr.vector_source_f (src_data)
 143         op = gr.iir_filter_ffd (fftaps, fbtaps)
 144         fftaps_data = (1)
 145         fbtaps = (0,0, -1,3)
 146         op.set_taps (fftaps, fbtaps)
 147         dst = gr.vector_sink_f ()
 148         self.tb.connect (src, op)
 149         self.tb.connect (op, dst)
 150         self.tb.run ()
 151         result_data = dst.data ()
 152         self.assertFloatTuplesAlmostEqual (expected_result, result_data)
 153
 154
 155
 156
 157 if __name__ == '__main__':
 158     gr_unittest.main ()
 159