git.gag.com Git - debian/gnuradio/blob - gr-atsc/src/lib/atsci_root_raised_cosine_bandpass.cc

   1 /* -*- c++ -*- */
   2 /*
   3  * Copyright 2002 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 #include <atsci_root_raised_cosine_bandpass.h>
  24 #include <iostream>
  25 #include <cmath>
  26
  27 using std::vector;
  28 using std::cerr;
  29 using std::endl;
  30
  31
  32 vector<float>
  33 atsc_root_raised_cosine_bandpass::taps (double sampling_freq)
  34 {
  35
  36   vector<float> t = atsc_root_raised_cosine::taps (sampling_freq);
  37
  38   cerr << "atsc_root_raised_cosine_bandpass::taps -- " << t.size () << endl;
  39
  40   // heterodyne the low pass coefficients up to the specified bandpass
  41   // center frequency.  Note that when we do this, the filter bandwidth
  42   // is effectively twice the low pass (2.69 * 2 = 5.38) and hence
  43   // matches the diagram in the ATSC spec.
  44
  45   double arg = 2 * M_PI * _center_freq / sampling_freq;
  46   for (unsigned int i = 0; i < t.size (); i++)
  47     // the factor of 2 keeps the gain of the passband normalized to 1
  48     t[i] *= 2 * cos (arg * (double) i);
  49
  50   return t;
  51 }