00001 // 00002 // Copyright(c) 1993-1996 Tony Kirke 00003 // author="Tony Kirke" * 00004 /* 00005 * SPUC - Signal processing using C++ - A DSP library 00006 * 00007 * This program is free software; you can redistribute it and/or modify 00008 * it under the terms of the GNU General Public License as published by 00009 * the Free Software Foundation; either version 2, or (at your option) 00010 * any later version. 00011 * 00012 * This program is distributed in the hope that it will be useful, 00013 * but WITHOUT ANY WARRANTY; without even the implied warranty of 00014 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 00015 * GNU General Public License for more details. 00016 * 00017 * You should have received a copy of the GNU General Public License 00018 * along with this program; if not, write to the Free Software 00019 * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. 00020 */ 00021 #ifndef FIRSPARSE 00022 #define FIRSPARSE 00023 #include <fir.h> 00024 namespace SPUC { 00028 00029 00030 template <class Numeric> class fir_sparse_coef : public fir<Numeric> 00031 { 00032 public: 00034 fir_sparse_coef<Numeric>(void) { } 00035 fir_sparse_coef<Numeric>(long n) : fir<Numeric>(n) { } 00036 fir_sparse_coef<Numeric>(const char* file) { read_taps(file); }; 00037 long rate; 00038 00039 void set_rate(long r) { rate = r; } 00040 void input(Numeric in) { 00041 int i; 00043 for (i=num_taps-1;i>0;i--) z[i] = z[i-1]; 00045 z[0] = in; 00046 } 00049 Numeric update(Numeric in) { 00050 input(in); 00051 // Perform FIR 00052 Numeric output; 00053 int i; 00054 for (output=0,i=0;i<num_taps/rate;i++) output += coeff[i]*z[i*rate]; 00055 return(output); 00056 } 00057 00058 }; 00059 } // namespace SPUC 00060 #endif