GNU Radio 3.6.5 C++ API
|
00001 /* -*- c++ -*- */ 00002 /* 00003 * Copyright 2004,2005 Free Software Foundation, Inc. 00004 * 00005 * This file is part of GNU Radio 00006 * 00007 * GNU Radio 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 3, or (at your option) 00010 * any later version. 00011 * 00012 * GNU Radio 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 GNU Radio; see the file COPYING. If not, write to 00019 * the Free Software Foundation, Inc., 51 Franklin Street, 00020 * Boston, MA 02110-1301, USA. 00021 */ 00022 00023 #ifndef INCLUDED_GR_SINGLE_POLE_IIR_FILTER_FF_H 00024 #define INCLUDED_GR_SINGLE_POLE_IIR_FILTER_FF_H 00025 00026 #include <gr_core_api.h> 00027 #include <gr_sync_block.h> 00028 #include <gr_single_pole_iir.h> 00029 #include <stdexcept> 00030 00031 class gr_single_pole_iir_filter_ff; 00032 typedef boost::shared_ptr<gr_single_pole_iir_filter_ff> gr_single_pole_iir_filter_ff_sptr; 00033 00034 GR_CORE_API gr_single_pole_iir_filter_ff_sptr 00035 gr_make_single_pole_iir_filter_ff (double alpha, unsigned int vlen=1); 00036 00037 /*! 00038 * \brief single pole IIR filter with float input, float output 00039 * 00040 * The input and output satisfy a difference equation of the form 00041 \htmlonly 00042 \f{ 00043 y[n] - (1-alpha) y[n-1] = alpha x[n] 00044 \f} 00045 \endhtmlonly 00046 00047 \xmlonly 00048 y[n] - (1-alpha) y[n-1] = alpha x[n] 00049 \endxmlonly 00050 00051 * with the corresponding rational system function 00052 \htmlonly 00053 \f{ 00054 H(z) = \frac{alpha}{1 - (1-alpha) z^{-1}} 00055 \f} 00056 \endhtmlonly 00057 00058 \xmlonly 00059 H(z) = \ frac{alpha}{1 - (1-alpha) z^{-1}} 00060 \endxmlonly 00061 00062 * Note that some texts define the system function with a + in the denominator. 00063 * If you're using that convention, you'll need to negate the feedback tap. 00064 */ 00065 class GR_CORE_API gr_single_pole_iir_filter_ff : public gr_sync_block 00066 { 00067 private: 00068 friend GR_CORE_API gr_single_pole_iir_filter_ff_sptr 00069 gr_make_single_pole_iir_filter_ff (double alpha, unsigned int vlen); 00070 00071 unsigned int d_vlen; 00072 std::vector<gr_single_pole_iir<float,float,double> > d_iir; 00073 00074 gr_single_pole_iir_filter_ff (double alpha, unsigned int vlen); 00075 00076 public: 00077 ~gr_single_pole_iir_filter_ff (); 00078 00079 void set_taps (double alpha); 00080 00081 int work (int noutput_items, 00082 gr_vector_const_void_star &input_items, 00083 gr_vector_void_star &output_items); 00084 }; 00085 00086 #endif