git.gag.com Git - debian/gnuradio/blob - gnuradio-core/src/lib/filter/gr_pfb_channelizer_ccf.h

   1 /* -*- c++ -*- */
   2 /*
   3  * Copyright 2009 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
  24 #ifndef INCLUDED_GR_PFB_CHANNELIZER_CCF_H
  25 #define INCLUDED_GR_PFB_CHANNELIZER_CCF_H
  26
  27 #include <gr_sync_block.h>
  28
  29 class gr_pfb_channelizer_ccf;
  30 typedef boost::shared_ptr<gr_pfb_channelizer_ccf> gr_pfb_channelizer_ccf_sptr;
  31 gr_pfb_channelizer_ccf_sptr gr_make_pfb_channelizer_ccf (unsigned int numchans,
  32                                                          const std::vector<float> &taps);
  33
  34 class gr_fir_ccf;
  35 class gri_fft_complex;
  36
  37
  38 /*!
  39  * \class gr_pfb_channelizer_ccf
  40  *
  41  * \brief Polyphase filterbank channelizer with
  42  *        gr_complex input, gr_complex output and float taps
  43  *
  44  * \ingroup filter_blk
  45  *
  46  * This block takes in complex inputs and channelizes it to <EM>M</EM>
  47  * channels of equal bandwidth. Each of the resulting channels is
  48  * decimated to the new rate that is the input sampling rate
  49  * <EM>fs</EM> divided by the number of channels, <EM>M</EM>.
  50  *
  51  * The PFB channelizer code takes the taps generated above and builds
  52  * a set of filters. The set contains <EM>M</EM> number of filters
  53  * and each filter contains ceil(taps.size()/decim) number of taps.
  54  * Each tap from the filter prototype is sequentially inserted into
  55  * the next filter. When all of the input taps are used, the remaining
  56  * filters in the filterbank are filled out with 0's to make sure each
  57  * filter has the same number of taps.
  58  *
  59  * Each filter operates using the gr_fir filter classs of GNU Radio,
  60  * which takes the input stream at <EM>i</EM> and performs the inner
  61  * product calculation to <EM>i+(n-1)</EM> where <EM>n</EM> is the
  62  * number of filter taps. To efficiently handle this in the GNU Radio
  63  * structure, each filter input must come from its own input
  64  * stream. So the channelizer must be provided with <EM>M</EM> streams
  65  * where the input stream has been deinterleaved. This is most easily
  66  * done using the gr_stream_to_streams block.
  67  *
  68  * The output is then produced as a vector, where index <EM>i</EM> in
  69  * the vector is the next sample from the <EM>i</EM>th channel. This
  70  * is most easily handled by sending the output to a
  71  * gr_vector_to_streams block to handle the conversion and passing
  72  * <EM>M</EM> streams out.
  73  *
  74  * The input and output formatting is done using a hier_block2 called
  75  * pfb_channelizer_ccf. This can take in a single stream and outputs
  76  * <EM>M</EM> streams based on the behavior described above.
  77  *
  78  * The filter's taps should be based on the input sampling rate.
  79  *
  80  * For example, using the GNU Radio's firdes utility to building
  81  * filters, we build a low-pass filter with a sampling rate of
  82  * <EM>fs</EM>, a 3-dB bandwidth of <EM>BW</EM> and a transition
  83  * bandwidth of <EM>TB</EM>. We can also specify the out-of-band
  84  * attenuation to use, <EM>ATT</EM>, and the filter window
  85  * function (a Blackman-harris window in this case). The first input
  86  *  is the gain of the filter, which we specify here as unity.
  87  *
  88  *      <B><EM>self._taps = gr.firdes.low_pass_2(1, fs, BW, TB,
  89  *           attenuation_dB=ATT, window=gr.firdes.WIN_BLACKMAN_hARRIS)</EM></B>
  90  *
  91  * The theory behind this block can be found in Chapter 6 of
  92  * the following book.
  93  *
  94  *    <B><EM>f. harris, "Multirate Signal Processing for Communication
  95  *       Systems," Upper Saddle River, NJ: Prentice Hall, Inc. 2004.</EM></B>
  96  *
  97  */
  98
  99 class gr_pfb_channelizer_ccf : public gr_sync_block
 100 {
 101  private:
 102   /*!
 103    * Build the polyphase filterbank decimator.
 104    * \param numchans (unsigned integer) Specifies the number of channels <EM>M</EM>
 105    * \param taps    (vector/list of floats) The prototype filter to populate the filterbank.
 106    */
 107   friend gr_pfb_channelizer_ccf_sptr gr_make_pfb_channelizer_ccf (unsigned int numchans,
 108                                                                   const std::vector<float> &taps);
 109
 110   std::vector<gr_fir_ccf*> d_filters;
 111   std::vector< std::vector<float> > d_taps;
 112   gri_fft_complex         *d_fft;
 113   unsigned int             d_numchans;
 114   unsigned int             d_taps_per_filter;
 115   bool                     d_updated;
 116
 117   /*!
 118    * Build the polyphase filterbank decimator.
 119    * \param numchans (unsigned integer) Specifies the number of channels <EM>M</EM>
 120    * \param taps    (vector/list of floats) The prototype filter to populate the filterbank.
 121    */
 122   gr_pfb_channelizer_ccf (unsigned int numchans,
 123                           const std::vector<float> &taps);
 124
 125 public:
 126   ~gr_pfb_channelizer_ccf ();
 127
 128   /*!
 129    * Resets the filterbank's filter taps with the new prototype filter
 130    * \param taps    (vector/list of floats) The prototype filter to populate the filterbank.
 131    */
 132   void set_taps (const std::vector<float> &taps);
 133
 134   /*!
 135    * Print all of the filterbank taps to screen.
 136    */
 137   void print_taps();
 138
 139   int work (int noutput_items,
 140             gr_vector_const_void_star &input_items,
 141             gr_vector_void_star &output_items);
 142 };
 143
 144 #endif