git.gag.com Git - debian/gnuradio/blob - gr-atsc/src/lib/atsci_vsbtx_lp.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 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 #include <atsc_consts.h>
  24 #include <atsci_vsbtx_lp.h>
  25 #include <stdexcept>
  26 #include <cmath>
  27 #include <iostream>
  28
  29 using std::vector;
  30 using std::cerr;
  31 using std::endl;
  32
  33 // atsc root raised cosine filter, sampling rate = 21.52 MHz
  34 static const float atsc_vsbtx_lp2x[] = {
  35 #include "atsc_vsbtx_lp.dat"
  36 };
  37
  38 #define NELEM(x) (sizeof (x) / sizeof ((x)[0]))
  39
  40 // is a within 5% of target?
  41
  42 static bool
  43 close_enough_p (double a, double target)
  44 {
  45   double delta = fabs (target * 0.05);  //  5 percent
  46
  47   return fabs (target - a) <= delta;
  48 }
  49
  50 vector<float>
  51 atsc_vsbtx_lp::taps (double sampling_freq)
  52 {
  53   if (close_enough_p (sampling_freq, 2 * ATSC_SYMBOL_RATE)){
  54     return vector<float>(&atsc_vsbtx_lp2x[0], &atsc_vsbtx_lp2x[NELEM(atsc_vsbtx_lp2x)]);
  55   }
  56   else
  57     throw std::out_of_range (
  58      "atsc_vsbtx_lp: no pre-designed filter close enough");
  59 }